Predicting antibody hypervariable loop conformations. II: Minimization and molecular dynamics studies of MCPC603 from many randomly generated loop conformations.

Abstract

We describe a method for predicting the conformations of loops in proteins and its application to four of the complementarity determining regions [CDRs] in the crystallographically determined structure of MCPC603. The method is based on the generation of a large number of randomly generated conformations for the backbone of the loop being studied, followed by either minimization or molecular dynamics followed by minimization starting from these random structures. The details of the algorithm for the generation of the loops are presented in the first paper in this series (Shenkin et al. [submitted]). The results of minimization and molecular dynamics applied to these loops is presented here. For the two shortest CDRs studied (H1 and L2, which are five and seven amino acids long), minimizations and dynamics simulations which ignore interactions of the loop amino acids beyond the carbon beta replicate the conformation of the crystal structure closely. This suggests that these loops fold independently of sequence variation. For the third CDR (L3, which is nine amino acids), those portions of the CDR near its base which are hydrogen bonded to framework are well replicated by our procedures, but the top of the loop shows significant conformational variability. This variability persists when side chain interactions for the MCPC603 sequence are included. For a fourth CDR (H3, which is 11 amino acids long), new low-energy backbone conformations are found; however, only those which are close to the crystal are compatible with the sequence when side chain interactions are taken into account. Results from minimization and dynamics on single CDRs with all other CDRs removed are presented. These allow us to explore the extent to which individual CDR conformations are determined by interactions with framework only.