Development of an extended simulated annealing method: application to the modeling of complementary determining regions of immunoglobulins.

An extended simulated annealing process (ESAP) has been developed in order to obtain an ensemble of conformations of a peptide segment from a protein fluctuating at a given temperature. The annealing process was performed with a fast Monte Carlo method using the scaled collective variables developed by Noguti and Go. The system was divided into two parts: one consists of one or more peptide segments and is flexible around the main-chain and side-chain torsional angles; the other represents the rest of the molecule and was maintained fixed at the atomic positions determined by x-ray experiments. The target function included the nonbonding atomic interactions and a distance function to anchor the N and C terminal ends of each segment to the fixed part. Three systems of complementary determining regions (CDR) of antibodies were tested and compared to x-ray data: L2 loop (7 residues) of the light chain of lambda-type Bence-Jones protein, H1 and the H2 loops (14 residues) of McPC603, and H1 and H2 loops (12 residues) of HyHEL-5. Each state of CDR conformations was characterized at room temperature by the average of their coordinates (average conformation) and the internal energy. With a limited number of annealing processes (10), starting from the extended conformation, we have obtained states with conformations close to the observed x-ray structures, from 1.1 to 1.7 A root mean square deviation (rmsd) of main-chain atoms depending on the system. These states were identical or within 0.25 A rmsd of those of lowest internal energy. For unknown CDR structures the criteria of lowest internal energies from ESAP can be used to predict hypervariable loop structures in antibodies with an accuracy comparable to other methods.     

Canonical antigen-binding loop structures in immunoglobulins: more structures, more canonical classes?

Grafting the antigen-binding loops onto a human antibody scaffold is a widely used technique to humanise murine antibodies. The success of this approach depends largely on the observation that the antigen-binding loops adopt only a limited number of canonical structures. Identification of the correct canonical structure is therefore essential. Algorithms that predict the main-chain conformation of the hypervariable loops using only the amino acid sequence often provide this information. Here, we describe new canonical loop conformations for the hypervariable regions H1 and H2 as found in single-domain antibody fragments of dromedaries or llama. Although the occurrence of these new loop conformations was not predicted by the algorithms used, it seems that they could occur in human or mouse antigen-binding loops. Their discovery indicates that the currently used set of canonical structures is incomplete and that the prediction algorithms should be extended to include these new structures.     

Ab initio construction of all-atom loop conformations

In this study, a new ab initio method named CLOOP has been developed to build all-atom loop conformations. In this method, a loop main-chain conformation is generated by sampling main-chain dihedral angles from a restrained varphi/psi set, and the side-chain conformations are built randomly. The CHARMM all-atom force field was used to evaluate the loop conformations. Soft core potentials were used to treat the non-bond interactions, and a designed energy-minimization technique was used to close and optimize the loop conformations. It is shown that the two strategies improve the computational efficiency and the loop-closure rate substantially compared to normal minimization methods. CLOOP was used to construct the conformations of 4-, 8-, and 12-residue loops in Fiser's test set. The average main-chain root-mean-square deviations obtained in 1,000 trials for the 10 different loops of each size are 0.33, 1.27, and 2.77 A, respectively. CLOOP can build all-atom loop conformations with a sampling accuracy comparable with previous loop main-chain construction algorithms. [Figure: see text].     

Molecular recognition. Conformational analysis of limited proteolytic sites and serine proteinase protein inhibitors

The conformations of known tryptic limited proteolytic sites have been analysed and compared to the structures of the binding regions of serine proteinase inhibitors, as they are found when complexed to a serine proteinase. Conformational parameters studied include main-chain torsion angles, root-mean-square fits, accessibility, mobility and protrusion indices. As observed before, the inhibitors share a common main-chain conformation at the binding loop from P3-P'3 (Schechter & Berger notation), which is maintained throughout all the serine proteinase inhibitor families for which X-ray data is available, despite lack of similarity in the rest of the protein. This canonical structure is not found amongst the limited proteolytic sites (or nicksites), which differ markedly from the inhibitor binding loop conformation, and also amongst themselves. The experimentally determined nicksites are in general both accessible and protruding; as are the inhibitor binding loops, as well as being typically flexible regions of structure, as denoted by elevated temperature factors from crystallographic determinations. For cleavage by serine proteinases these loops must radically alter their local conformations and a large motion of the loop relative to the structure, in some cases, would be required to orientate these sites for cleavage.     

Structural determinants of the conformations of medium-sized loops in proteins

Loops are integral components of protein structures, providing links between elements of secondary structure, and in many cases contributing to catalytic and binding sites. The conformations of short loops are now understood to depend primarily on their amino acid sequences. In contrast, the structural determinants of longer loops involve hydrogen-bonding and packing interactions within the loop and with other parts of the protein. By searching solved protein structures for regions similar in main chain conformation to the antigen-binding loops in immunoglobulins, we identified medium-sized loops of similar structure in unrelated proteins, and compared the determinants of their conformations. For loops that form compact substructures the major determinant of the conformation is the formation of hydrogen bonds to inward-pointing main chain atoms. For loops that have more extended conformations, the major determinant of their structure is the packing of a particular residue or residues against the rest of the protein. The following picture emerges: Medium-sized loops of similar conformation are stabilized by similar interactions. The groups that interact with the loop have very similar spatial dispositions with respect to the loop. However, the residues that provide these interactions may arise from dissimilar parts of the protein: The conformation of the loop requires certain interactions that the protein may provide in a variety of ways.     

Conformational analysis and clustering of short and medium size loops connecting regular secondary structures: a database for modeling and prediction.

Loops are regions of nonrepetitive conformation connecting regular secondary structures. We identified 2,024 loops of one to eight residues in length, with acceptable main-chain bond lengths and peptide bond angles, from a database of 223 protein and protein-domain structures. Each loop is characterized by its sequence, main-chain conformation, and relative disposition of its bounding secondary structures as described by the separation between the tips of their axes and the angle between them. Loops, grouped according to their length and type of their bounding secondary structures, were superposed and clustered into 161 conformational classes, corresponding to 63% of all loops. Of these, 109 (51% of the loops) were populated by at least four nonhomologous loops or four loops sharing a low sequence identity. Another 52 classes, including 12% of the loops, were populated by at least three loops of low sequence similarity from three or fewer nonhomologous groups. Loop class suprafamilies resulting from variations in the termini of secondary structures are discussed in this article. Most previously described loop conformations were found among the classes. New classes included a 2:4 type IV hairpin, a helix-capping loop, and a loop that mediates dinucleotide-binding. The relative disposition of bounding secondary structures varies among loop classes, with some classes such as beta-hairpins being very restrictive. For each class, sequence preferences as key residues were identified; those most frequently at these conserved positions than in proteins were Gly, Asp, Pro, Phe, and Cys. Most of these residues are involved in stabilizing loop conformation, often through a positive phi conformation or secondary structure capping. Identification of helix-capping residues and beta-breakers among the highly conserved positions supported our decision to group loops according to their bounding secondary structures. Several of the identified loop classes were associated with specific functions, and all of the member loops had the same function; key residues were conserved for this purpose, as is the case for the parvalbumin-like calcium-binding loops. A significant number, but not all, of the member loops of other loop classes had the same function, as is the case for the helix-turn-helix DNA-binding loops. This article provides a systematic and coherent conformational classification of loops, covering a broad range of lengths and all four combinations of bounding secondary structure types, and supplies a useful basis for modelling of loop conformations where the bounding secondary structures are known or reliably predicted.     

Anti-insulin antibody structure and conformation. II. Molecular dynamics with explicit solvent.

Molecular dynamics at 300 K was used as a conformation searching tool to analyze a knowledge-based structure prediction of an anti-insulin antibody. Solvation effects were modeled by packing water molecules around the antigen binding loops. Some loops underwent backbone and side-chain conformational changes during the 95-ps equilibration, and most of these new, lower potential energy conformations were stable during the subsequent 200-ps simulation. Alterations to the model include changes in the intraloop, main-chain hydrogen bonding network of loop H3, and adjustments of Tyr and Lys side chains of H3 induced by hydrogen bonding to water molecules. The structures observed during molecular dynamics support the conclusion of the previous paper that hydrogen bonding will play the dominant role in antibody-insulin recognition. Determination of the structure of the antibody by x-ray crystallography is currently being pursued to provide an experimental test of these results. The simulation appears to improve the model, but longer simulations at higher temperatures should be performed.     

Electronic and vibrational circular dichroism of aromatic amino acids by density functional theory

Structures of model compounds mimicking aromatic amino acid residues in proteins are optimized by density functional theory (DFT), assuming that the main-chain conformation was a random coil. Excitation energies and dipole and rotational strengths for the optimized structures were calculated based on time-dependent DFT (TD-DFT). The electronic circular dichroism (ECD) bands of the models were significantly affected by side-chain conformations. Hydration models of the aromatic residues were also subjected to TD-DFT calculations, and the ECD bands of these models were found to be highly perturbed by the hydration of the main-chain amide groups. In addition to calculating the random-coil conformation, we also performed TD-DFT calculations of the aromatic residue models, assuming that the main-chain conformation was an alpha-helix or beta-strand. As expected, the overall feature of the ECD bands was also perturbed by the main-chain conformations. Moreover, vibrational circular dichroism (VCD) spectra of the hydration models in a random-coil structure were simulated by DFT, which showed that the VCD spectra are more sensitive to the side-chain conformations than the ECD spectra. The present results show that analyses combining ECD and VCD spectroscopy and using DFT calculations can elucidate the main- and side-chain conformations of aromatic residues in proteins.     

Ab initio construction of polypeptide fragments: efficient generation of accurate,representative ensembles

We describe a novel method to generate ensembles of conformations of the main-chain atoms [N, C(alpha), C, O, Cbeta] for a sequence of amino acids within the context of a fixed protein framework. Each conformation satisfies fundamental stereo-chemical restraints such as idealized geometry, favorable phi/psi angles, and excluded volume. The ensembles include conformations both near and far from the native structure. Algorithms for effective conformational sampling and constant time overlap detection permit the generation of thousands of distinct conformations in minutes. Unlike previous approaches, our method samples dihedral angles from fine-grained phi/psi state sets, which we demonstrate is superior to exhaustive enumeration from coarse phi/psi sets. Applied to a large set of loop structures, our method samples consistently near-native conformations, averaging 0.4, 1.1, and 2.2 A main-chain root-mean-square deviations for four, eight, and twelve residue long loops, respectively. The ensembles make ideal decoy sets to assess the discriminatory power of a selection method. Using these decoy sets, we conclude that quality of anchor geometry cannot reliably identify near-native conformations, though the selection results are comparable to previous loop prediction methods. In a subsequent study (de Bakker et al.: Proteins 2003;51:21-40), we demonstrate that the AMBER forcefield with the Generalized Born solvation model identifies near-native conformations significantly better than previous methods.     

Neutron diffraction study of carbonmonoxymyoglobin.

Neutron diffraction data from a crystal of carbonmonoxymyoglobin were refined by PROLSQ, a modern restrained least-squares procedure in reciprocal space, in conjunction with a solvent analysis technique, to a final R-factor of 11.3%. The ligand CO occupies two sites and its binding conformations are distorted from the linear conformation. The N epsilon atom of the distal histidine residue is deprotonated (not deuterated), and a water molecule is bound to the N delta atom of the distal histidine. The side-chain of Lys56 (D6) exists in two alternative charge-binding sites. His24 (B5) and His119 (GH1) share a hydrogen atom. His12 (A10) and His36 (C1) are deprotonated. The deprotonated imidazole ring of His12 (A10) may act as a hydrogen-bond acceptor. The heme group is planar within 0.09 A root-mean-square (r.m.s.) deviation from planarity. The solvent environments for the two propionic acid groups are different. The side-chain of Arg45 (CD3) forms hydrogen bonds with the side-chain of Asp60 (E3) and one of the two propionic acid groups. An average N-2H . . . O angle in helical regions is 147 (+/- 11) degrees. Eleven main-chain amide hydrogen atoms from hydrophobic residues do not exchange with deuterium. The overall atomic occupancy factors for the main-chain and side-chain atoms are quite uniform, at 0.97 (+/- 0.07) and 0.93 (+/- 0.10), respectively, as shown by an occupancy analysis made at the end of the refinement procedure.     

Context-dependent conformation of diethylglycine residues in peptides.

Diethylglycine (Deg) residues incorporated into peptides can stabilize fully extended (C5) or helical conformations. The conformations of three tetrapeptides Boc-Xxx-Deg-Xxx-Deg-OMe (Xxx=Gly, GD4; Leu, LD4 and Pro, PD4) have been investigated by NMR. In the Gly and Leu peptides, NOE data suggest that the local conformations at the Deg residues are fully extended. Low temperature coefficients for the Deg(2) and Deg(4) NH groups are consistent with their inaccessibility to solvent, in a C5 conformation. NMR evidence supports a folded beta-turn conformation involving Deg(2)-Gly(3), stabilized by a 4-->1 intramolecular hydrogen bond between Pro(1) CO and Deg(4) NH in the proline containing peptide (PD4). The crystal structure of GD4 reveals a hydrated multiple turn conformation with Gly(1)-Deg(2) adopting a distorted type II/II' conformation, while the Deg(2)-Pro(3) segment adopts a type III/III' structure. A lone water molecule is inserted into the potential 4-->1 hydrogen bond of the Gly(1)-Deg(2) beta-turn.     

Dihydrofolate reductase from Lactobacillus casei. Stereochemistry of NADPH binding.

The NADPH molecule binds to dihydrofolate reductase in an extended conformation. Several of the individual dihedral angles, especially in the adenine mononucleotide portion of the coenzyme, differ from their minimum energy conformations. The ribose phosphate portions of the coenzyme are involved in numerous specific hydrogen-bonded and charge-charge interactions. The adenine ring resides in an apparently nonspecific hydrophobic cleft and the nicotinamide ring is bound within an intricately constructed cavity, one wall of which includes the pyrazine ring of bound methotrexate. Two rather extended loops (residues 10 to 24 and 117 to 135) connecting beta A to alpha B and beta F to beta G, respectively, move 2 to 3 A when NADPH binds to dihydrofolate reductase. No overall structural homology is evident between the dinucleotide binding domains of dihydrofolate reductase on the one hand and the four NAD+-dependent dehydrogenases of known structure on the other. However, binding does occur in both cases at the carboxyl edge of a region of parallel beta sheet flanked by a pair of alpha helices.     

Properties of polyproline II, a secondary structure element implicated in protein-protein interactions

The polyproline II (PPII) conformation of protein backbone is an important secondary structure type. It is unusual in that, due to steric constraints, its main-chain hydrogen-bond donors and acceptors cannot easily be satisfied. It is unable to make local hydrogen bonds, in a manner similar to that of alpha-helices, and it cannot easily satisfy the hydrogen-bonding potential of neighboring residues in polyproline conformation in a manner analogous to beta-strands. Here we describe an analysis of polyproline conformations using the HOMSTRAD database of structurally aligned proteins. This allows us not only to determine amino acid propensities from a much larger database than previously but also to investigate conservation of amino acids in polyproline conformations, and the conservation of the conformation itself. Although proline is common in polyproline helices, helices without proline represent 46% of the total. No other amino acid appears to be greatly preferred; glycine and aromatic amino acids have low propensities for PPII. Accordingly, the hydrogen-bonding potential of PPII main-chain is mainly satisfied by water molecules and by other parts of the main-chain. Side-chain to main-chain interactions are mostly nonlocal. Interestingly, the increased number of nonsatisfied H-bond donors and acceptors (as compared with alpha-helices and beta-strands) makes PPII conformers well suited to take part in protein-protein interactions.     

Modeling studies of the change in conformation required for cleavage of limited proteolytic sites.

Previous analyses of limited proteolytic sites within native, folded protein structures have shown that a significant conformational change is required in order to facilitate binding into the active site of the attacking proteinase. For the serine proteinases, the optimum conformation to match the proteinase binding-site geometry has been well characterized crystallographically by the conserved main-chain geometry of the reactive site loops of their protein inhibitors. A good substrate must adopt a conformation very similar to this "target" main-chain conformation prior to cleavage. Using a "loop-closure" modeling approach, we have tested the ability of a set of tryptic-limited proteolytic sites to achieve this target conformation and further tested their suitability for cleavage. The results show that in most cases, significant changes in the conformation of at least 12 residues are required. All the putative tryptic cleavage sites in 1 protein, elastase, were also modeled and tested to compare the results to the actual nicksite in that protein. These results strongly suggest that large local motions proximate to the scissile bond are required for proteolysis, and it is this ability to unfold locally without perturbing the overall protein conformation that is the prime determinant for limited proteolysis.     

Conformational studies of irreversible HIV-1 protease inhibitors containing cis-epoxide as an amide isostere.

We have carried out NMR and molecular modeling studies of peptidomimetic HIV-1 protease inhibitors, LB71116: Qc-Asn-Phepsi[(1R,2S)-cis-epoxide]Gly-NH-CH(isopropyl)2 where Qc stands for quinaldic acid and LB71148: Qc-(SMe)Pen(O)2-Phepsi[(1R,2S)-cis-epoxide]Gly-NH-CH(isoprop yl)2 where (SMe)Pen(O)2 stands for S-methyl-S-dioxo-penicillamine. Through conformational calculations and NMR data analysis, we have obtained preferred conformations of the two inhibitors in solution. To our knowledge, this work is one of the first extensive conformational studies of peptidomimetics containing cis-epoxide amide isostere. The resulting preferred conformations contain extended structures. In these conformations, the psi of Phe(cep) is maintained about 130 degrees and the phi angle of (cep)Gly prefers +/- 150 degrees [where Phe(cep) and (cep)Gly are the residues generated by the replacement of the Phe-Gly peptide bond with cis-epoxide]. Two conformations were commonly observed in the preferred conformations of each inhibitor. Through restrained molecular dynamics simulating the hydrogen bond formation between our inhibitor and a water molecule ('flap water'), one of the conformations is assumed as the conformation which can bind to the enzyme without large conformational changes. Recently, we had the opportunity to compare the selected preferred conformation with the binding conformation of LB71116 observed from the X-ray studies of the complex between LB71116 and HIV-1 protease. These two conformations are surprisingly similar to each other. Thus, we can explain high activity and selectivity of our inhibitors to the HIV-1 protease by the similarity between the preferred conformations in solution and the binding conformation.     

Molecular dynamics studies of the conformation of sorbitol

Molecular dynamics simulations of a 3 molal aqueous solution of d-sorbitol (also called d-glucitol) have been performed at 300 K, as well as at two elevated temperatures to promote conformational transitions. In principle, sorbitol is more flexible than glucose since it does not contain a constraining ring. However, a conformational analysis revealed that the sorbitol chain remains extended in solution, in contrast to the bent conformation found experimentally in the crystalline form. While there are 243 staggered conformations of the backbone possible for this open-chain polyol, only a very limited number were found to be stable in the simulations. Although many conformers were briefly sampled, only eight were significantly populated in the simulation. The carbon backbones of all but two of these eight conformers were completely extended, unlike the bent crystal conformation. These extended conformers were stabilized by a quite persistent intramolecular hydrogen bond between the hydroxyl groups of carbon C-2 and C-4. The conformational populations were found to be in good agreement with the limited available NMR data except for the C-2–C-3 torsion (spanned by the O-2–O-4 hydrogen bond), where the NMR data support a more bent structure.     

Conformational studies on somatostatin. II. The C-terminal hexapeptide fragment

The results of a conformational study on the C terminal hexapeptide of Somatostatin are presented. Semi-empirical energy calculations and high resolution NMR methods have been used to obtain information on the conformational properties of SRIF9-14 in [2H6]dimethylsulfoxide and 2H2O. It is concluded from the energy calculations that the peptide has an averaged conformation in which semi extended and folded structures are important. Only some of the folded conformations can explain the chemical shift differences between the amino acid residues Thr10 and Thr12 as a ring current shift by the Phe11 aromatic ring on Thr10. The nonequivalence is more pronounced in dimethyl-sulfoxide (0.23--0.15 ppm) where it decreases with increasing temperature towards the temperature independent value in 2H2O (0.03 ppm). This suggests that the folded conformations are somewhat predominant in dimethylsulfoxide solutions. In 2H2O the semi extended and folded structures are statistically equally important and the peptide is more flexible. A comparison with a study on the smaller fragments SRIF10-12 and SRIF10-13 which have similar conformational properties, demonstrates the usefulness of the fragment approach in conformational studies of peptides.     

Refined 1.83 A structure of trypanosomal triosephosphate isomerase crystallized in the presence of 2.4 M-ammonium sulphate. A comparison with the structure of the trypanosomal triosephosphate isomerase-glycerol-3-phosphate complex

Triosephosphate isomerase (TIM) is a dimeric glycolytic enzyme. TIM from Trypanosoma brucei brucei has been crystallized at pH 7.0 in 2.4 M-ammonium sulphate. The well-diffracting crystals have one dimer per asymmetric unit. The structure has been refined at 1.83 A resolution with an R-factor of 18.3% for all data between 6 A and 1.83 A (37,568 reflections). The model consists of 3778 protein atoms and 297 solvent atoms. Subunit 1 is involved in considerably more crystal contacts than subunit 2. Correlated with these differences in crystal packing is the observation that only in the active site of subunit 2 is a sulphate ion bound. Furthermore, significant differences with respect to structure and flexibility are observed in three loops near the active site. In particular, there is a 7 A positional difference of the tip of the flexible loop (loop 6) when comparing subunit 1 and subunit 2. Also, the neighbouring loops (loop 5 and loop 7) have significantly different conformations and flexibility. In subunit 1, loop 6 is in an "open" conformation, in subunit 2, loop 6 is in an "almost closed" conformation. Only in the presence of a phosphate-containing ligand, such as glycerol-3-phosphate, does loop 6 take up the "closed" conformation. Loop 6 and loop 7 (and also to some extent loop 5) are rather flexible in the almost closed conformation, but well defined in the open and closed conformations. The closing of loop 6 (167 to 180), as observed in the almost closed conformation, slightly changes the main-chain conformation of the catalytic glutamate, Glu167, leading to a change of the chi 1 angle of this residue from approximately -60 degrees to approximately 60 degrees and the weakening of the hydrogen bonds between its polar side-chain atoms and Ser96. In the closed conformation, in the presence of glycerol-3-phosphate, the main-chain atoms of Glu167 remain in the same position as in the almost closed conformation, but the side-chain has rotated around the CA-CB bond changing chi 1 from approximately 60 degrees to approximately -60 degrees. In this new position the hydrogen bonding to Ser96 is completely lost and also a water-mediated salt bridge between OE2(Glu167) and NE(Arg99) is lost. Comparison of the two independently refined subunits, showed that the root-mean-square deviation for all 249 CA atoms is 0.9 A; for the CA atoms of the beta-strands this is only 0.2 A. The average B-factor for all subunit 1 and subunit 2 atoms is 20 A2 and 25 A2, respectively.(ABSTRACT TRUNCATED AT 400 WORDS)