Variants of 3(10)-helices in proteins

An analysis of the shortest 3(10)-helices, containing three helical residues and two flanking capping residues that participate in two consecutive i + 3 --> i hydrogen bonds, shows that not all helices belong to the classic 3(10)-helix, where the three central residues adopt the right-handed helical conformation (alpha(R)). Three variants identified are: 3L10-helix with all residues in the left-handed helical region (alpha(L)), 3EL10-helix where the first residue is in the extended region followed by two residues in the alpha(L) conformation, and its mirror-image, the 3E'R10-helix. In the context of these helices, as well as the equivalent variants of alpha-helices, the length dependence of the handedness of secondary structures in protein structure is discussed. There are considerable differences in the amino acid preferences at different positions in the various types of 3(10)-helices. Each type of 3(10)-helix can be thought to be made up of an extension of a particular type of beta-turn (made up of residues i to i + 3) such that the (i + 3)th residue assumes the same conformation as the preceding residue. Distinct residue preferences at i and i + 3 positions seem to decide whether a particular stretch of four residues will be a beta-turn or a 3(10)-helix in the folded structure.     

Solid state and solution conformations of a hybrid alphagammaalphaalphagammaalpha hexapeptide. Characterization of a backbone expanded analog of the alpha-polypeptide 3(10)-helix

The stereochemically constrained gamma amino acid residue gabapentin (1-(aminomethyl)cyclohexaneacetic acid, Gpn) has been incorporated into a host alpha-peptide sequence. The structure of a hybrid alphagammaalphaalphagammaalpha peptide, Boc-Leu-Gpn-Aib-Leu-Gpn-Aib-OMe in crystals reveals a continuous helical conformation stabilized by three intramolecular 4 --> 1 C(12) hydrogen bonds across the alphagamma/alphagamma segments and one C(10) hydrogen bond across the central alphaalpha segment. This conformation corresponds to an expanded analog of the canonical all-alpha polypeptide 3(10)-helix, with insertion of two additional backbone atoms at each gamma residue. Solvent dependence of NH chemical shifts in CDCl(3) solution are consistent with conformation in which the NH groups of Aib (3), Leu (4), Gpn (5), and Aib (6) are hydrogen bonded, a feature observed in the solid state. The nonsequential NOEs between Gpn (2) NH <--> Leu (4) NH and Gpn (2) NH <--> Gpn (5) NH support the presence of additional conformations in solution. Temperature-dependent line broadening of NH resonances confirms the occurrence of rapid exchange between multiple conformations at room temperature. Two conformational models which rationalize the observed nonsequential NOEs are presented, both of which contain three hydrogen bonds and are consistent with the known stereochemical preferences of the Gpn residue.     

Conformational features of a hexapeptide model Ac-TGAAKA-NH2 corresponding to a hydrated alpha helical segment from glyceraldehyde 3-phosphate dehydrogenase: implications for the role of turns in helix folding

Recent analysis of alpha helices in protein crystal structures, available in literature, revealed hydrated alpha helical segments in which, water molecule breaks open helix 5-->1 hydrogen bond by inserting itself, hydrogen bonds to both C=O and NH groups of helix hydrogen bond without disrupting the helix hydrogen bond, and hydrogen bonds to either C=O or NH of helix hydrogen bond. These hydrated segments display a variety of turn conformations and are thought to be 'folding intermediates' trapped during folding-unfolding of alpha helices. A role for reverse turns is implicated in the folding of alpha helices. We considered a hexapeptide model Ac-1TGAAKA6-NH2 from glyceraldehyde 3-phosphate dehydrogenase, corresponding to a hydrated helical segment to assess its role in helix folding. The sequence is a site for two 'folding intermediates'. The conformational features of the model peptide have been investigated by 1H 2D NMR techniques and quantum mechanical perturbative configuration interaction over localized orbitals (PCILO) method. Theoretical modeling largely correlates with experimental observations. Based upon the amide proton temperature coefficients, the observed d alpha n(i, i + 1), d alpha n(i, i + 2), dnn(i, i + 1), d beta n(i, i + 1) NOEs and the results from theoretical modeling, we conclude that the residues of the peptide sample alpha helical and neck regions of the Ramachandran phi, psi map with reduced conformational entropy and there is a potential for turn conformations at N and C terminal ends of the peptide. The role of reduced conformational entropy and turn potential in helix formation have been discussed. We conclude that the peptide sequence can serve as a 'folding intermediate' in the helix folding of glyceraldehyde 3-phosphate dehydrogenase.     

Structural analysis of peptide helices containing centrally positioned lactic acid residues

The effect of insertion of lactic acid (Lac) residues into peptide helices has been probed using specifically designed sequences. The crystal structures of 11-residue and 14-residue depsipeptides Boc-Val-Val-Ala-Leu-Val-Lac-Leu-Aib-Val-Ala-Leu-OMe (1) and Boc-Val-Ala-Leu-Aib-Val-Ala-Leu-Val-Lac-Leu-Aib-Val-Ala-Leu-OMe (3), containing centrally positioned Lac residues, have been determined. The structure of an 11-residue peptide Boc-Val-Ala-Leu-Aib-Val-Ala-Leu-Aib-Val-Ala-Leu-OMe (2), analog of a which is an amide previously determined Lac-containing depsipeptide, Boc-Val-Ala-Leu-Aib-Val-Lac-Leu-Aib-Val-Ala-Leu-OMe I. L. Karle, C. Das, and P. Balaram, Biopolymers, Vol. 59, (2001) pp. 276-289], is also reported. Peptide 1 adopts a helical fold, which is stabilized by mixture of 4-->1 and 5-->1 hydrogen bonds. Peptide 2 adopts a completely alpha-helical conformation stabilized by eight successive 5-->1 hydrogen bonds. Peptide 3 appears to be predominately alpha-helical, with seven 5-->1 hydrogen bonds and three 4-->1 interaction interspersed in the sequence. In the structure of peptide 3 in addition to water molecules in the head-to-tail region, hydration at an internal segment of the helix is also observed. A comparison of five related peptide helices, containing a single Lac residue, reveals that the hydroxy acid can be comfortably accommodated at interior positions in the helix, with the closest C=O...O distances lying between 2.8 and 3.3 A.     

Structural and thermodynamics characters of isolated &alpha;-syn12 peptide: long-time temperature replica-exchange molecular dynamics in aqueous solution

The structural and thermodynamics characters of α-syn12 (residues 1-12 of the human α-synuclein protein) peptide in aqueous solution were investigated through temperature replica-exchange molecular dynamics (T-REMD) simulations with the GROMOS 43A1 force field. The two independent T-REMD simulations were completed starting from an initial conformational α-helix and an irregular structure, respectively. Each replica was run for 300 ns. The structural and thermodynamics characters were studied based on parameters such as distributions of backbone dihedral angles, free energy surface, stability of folded β-hairpin structure, and favorite conformations. The results showed that the isolated α-syn12 peptide in water adopted four different conformational states: the first state was a β-hairpin ensemble with Turn(9-6) and four hydrogen bonds, the second state was a β-hairpin ensemble with two turns (Turn(9-6) and Turn(5-2)) and three hydrogen bonds, the third state was a disordered structure with both Turn(8-5) and Turn(5-2), and the last state was a π-helix ensemble. Meanwhile, we studied the free energy change of α-syn12 peptide from the unfolded state to the β-hairpin state, which was in good agreement with the experiments and molecular dynamics simulations for some other peptides. We also analyzed the driving force of the peptide transition. The results indicated that the driving forces were high solvent exposure of hydrophobic Leu8 and hydrophobic residues in secondary structure. To our knowledge, this was the first report to study the isolated α-syn12 peptide in water by T-REMD.     

Helix bending in alamethicin: molecular dynamics simulations and amide hydrogen exchange in methanol

Molecular dynamics simulations of alamethicin in methanol were carried out with either a regular alpha-helical conformation or the x-ray crystal structure as starting structures. The structures rapidly converged to a well-defined hydrogen-bonding pattern with mixed alpha-helical and 3(10)-helical hydrogen bonds, consistent with NMR structural characterization, and did not unfold throughout the 1-ns simulation, despite some sizable backbone fluctuations involving reversible breaking of helical hydrogen bonds. Bending of the helical structure around residues Aib10-Aib13 was associated with reversible flips of the peptide bonds involving G11 (Aib10-G11 or G11-L12 peptide bonds), yielding discrete structural states in which the Aib10 carbonyl or (rarely) the G11 carbonyl was oriented away from the peptide helix. These peptide bond reversals could be accommodated without greatly perturbing the adjacent helical structure, and intramolecular hydrogen bonding was generally maintained in bent states through the formation of new (non-alpha or 3[10]) hydrogen bonds with good geometries: G11 NH-V9 CO (inverse gamma turn), Aib13 NH-Aib8 CO (pi-helix) and, rarely, L12 NH- Q7 NH (pi-helix). These observations may reconcile potentially conflicting NMR structural information for alamethicin in methanol, in which evidence for conformational flexibility in the peptide sequence before P14 (G11-Aib13) contrasts with the stability of backbone amide NH groups to exchange with solvent. Similar reversible reorientation of the Thr11-Gly12 peptide bond of melittin is also observed in dynamics simulations in methanol (R. B. Sessions, N. Gibbs, and C. E. Dempsey, submitted). This phenomenon may have some role in the orientation of the peptide carbonyl in solvating the channel lumen in membrane ion channel states of these peptides.     

Addition of side-chain interactions to 3(10)-helix/coil and alpha-helix/3(10)-helix/coil theory.

An increasing number of experimental and theoretical studies have demonstrated the importance of the 3(10)-helix/ alpha-helix/coil equilibrium for the structure and folding of peptides and proteins. One way to perturb this equilibrium is to introduce side-chain interactions that stabilize or destabilize one helix. For example, an attractive i, i + 4 interaction, present only in the alpha-helix, will favor the alpha-helix over 3(10), while an i, i + 4 repulsion will favor the 3(10)-helix over alpha. To quantify the 3(10)/alpha/coil equilibrium, it is essential to use a helix/coil theory that considers the stability of every possible conformation of a peptide. We have previously developed models for the 3(10)-helix/coil and 3(10)-helix/alpha-helix/ coil equilibria. Here we extend this work by adding i, i + 3 and i, i + 4 side-chain interaction energies to the models. The theory is based on classifying residues into alpha-helical, 3(10)-helical, or nonhelical (coil) conformations. Statistical weights are assigned to residues in a helical conformation with an associated helical hydrogen bond, a helical conformation with no hydrogen bond, an N-cap position, a C-cap position, or the reference coil conformation plus i, i + 3 and i, i + 4 side-chain interactions. This work may provide a framework for quantitatively rationalizing experimental work on isolated 3(10)-helices and mixed 3(10)-/alpha-helices and for predicting the locations and stabilities of these structures in peptides and proteins. We conclude that strong i, i + 4 side-chain interactions favor alpha-helix formation, while the 3(10)-helix population is maximized when weaker i, i + 4 side-chain interactions are present.     

The structure of the ends of α-helices in globular proteins: effect of additional hydrogen bonds and implications for helix formation

We prepared a set of about 2000 α-helices from a relational database of high-resolution three-dimensional structures of globular proteins, and identified additional main chain i ← i+3 hydrogen bonds at the ends of the helices (i.e., where the hydrogen bonding potential is not fulfilled by canonical i ← i+4 hydrogen bonds). About one-third of α-helices have such additional hydrogen bonds at the N-terminus, and more than half do so at the C-terminus. Although many of these additional hydrogen bonds at the C-terminus are associated with Schellman loops, the majority are not. We compared the dihedral angles at the termini of α-helices having or lacking the additional hydrogen bonds. Significant differences were found, especially at the C-terminus, where the dihedral angles at positions C2 and C1 in the absence of additional hydrogen bonds deviate substantially from those occurring within the α-helix. Using a novel approach we show how the structure of the C-terminus of the α-helix can emerge from that of constituent overlapping α-turns and β-turns, which individually show a variation in dihedral angles at different positions. We have also considered the direction of propagation of the α-helix using this approach. If one assumes that helices start as a single α-turn and grow by successive addition of further α-turns, the paths for growth in the N → C and C → N directions differ in a way that suggests that extension in the C → N direction is favored.     

3(10)-helices in proteins are parahelices

The 3(10)-helix is characterized by having at least two consecutive hydrogen bonds between the main-chain carbonyl oxygen of residue i and the main-chain amide hydrogen of residue i + 3. The helical parameters--pitch, residues per turn, radius, and root mean square deviation (rmsd) from the best-fit helix--were determined by using the HELFIT program. All 3(10)-helices were classified as regular or irregular based on rmsd/(N - 1)1/2 where N is the helix length. For both there are systematic, position-specific shifts in the backbone dihedral angles. The average phi, psi shift systematically from approximately -58 degrees, approximately -32 degrees to approximately -90 degrees, approximately -4 degrees for helices 5, 6, and 7 residues long. The same general pattern is seen for helices, N = 8 and 9; however, in N = 9, the trend is repeated with residues 6, 7, and 8 approximately repeating the phi, psi of residues 2, 3, and 4. The residues per turn and radius of regular 3(10)-helices decrease with increasing length of helix, while the helix pitch and rise per residue increase. That is, regular 3(10)-helices become thinner and longer as N increases from 5 to 8. The fraction of regular 3(10)-helices decreases linearly with helix length. All longer helices, N > or = 9 are irregular. Energy minimizations show that regular helices become less stable with increasing helix length. These findings indicate that the definition of 3(10)-helices in terms of average, uniform dihedral angles is not appropriate and that it is inherently unstable for a polypeptide to form an extended, regular 3(10)-helix. The 3(10)-helices observed in proteins are better referred to parahelices.     

Direct computation of long time processes in peptides and proteins: reaction path study of the coil-to-helix transition in polyalanine.

The MaxFlux reaction path algorithm was used to isolate optimal transition pathways for the coil-to-helix transition in polyalanine. Eighteen transition pathways, each connecting one random coil configuration with an ideal alpha-helical configuration, were computed and analyzed. The transition pathway energetics and mechanism were analyzed in terms of the progression of the peptide nonbonded contact formation, helicity, end-to-end distance and energetics. It was found that (1) localized turns characterized by i, i + 3 hydrogen bonds form in the early stages of the coil-to-helix transition, (2) the peptide first collapses and then becomes somewhat more extended in the final stage of helix formation, and (3) 310-helix formation does not appear to be a necessary step in the transition from coil to helix. These conclusions are in agreement with the results of more computationally intensive direct molecular dynamics simulations. Proteins 1999;36:249-261.     

Conformation of di-n-propylglycine residues (Dpg) in peptides: crystal structures of a type I' beta-turn forming tetrapeptide and an alpha-helical tetradecapeptide.

The crystal structures of two oligopeptides containing di-n-propylglycine (Dpg) residues, Boc-Gly-Dpg-Gly-Leu-OMe (1) and Boc-Val-Ala-Leu-Dpg-Val-Ala-Leu-Val-Ala-Leu-Dpg-Val-Ala-Leu-OMe (2) are presented. Peptide 1 adopts a type I'beta-turn conformation with Dpg(2)-Gly(3) at the corner positions. The 14-residue peptide 2 crystallizes with two molecules in the asymmetric unit, both of which adopt alpha-helical conformations stabilized by 11 successive 5 --> 1 hydrogen bonds. In addition, a single 4 --> 1 hydrogen bond is also observed at the N-terminus. All five Dpg residues adopt backbone torsion angles (phi, psi) in the helical region of conformational space. Evaluation of the available structural data on Dpg peptides confirm the correlation between backbone bond angle N-C(alpha)-C' (tau) and the observed backbone phi,psi values. For tau > 106 degrees, helices are observed, while fully extended structures are characterized by tau < 106 degrees. The mean tau values for extended and folded conformations for the Dpg residue are 103.6 degrees +/- 1.7 degrees and 109.9 degrees +/- 2.6 degrees, respectively.     

Conformational preferences of oligopeptides rich in alpha-aminoisobutyric acid. II. A model for the 3(10)/alpha-helix transition with composition and sequence sensitivity.

The analysis of the factors that control the helical folding of Aib-rich peptides is extended to include sensitivity to sequence patterns, and in particular the presence of contiguous non-Aib alpha-mono-alkylated residues. The distinct hydrogen-bonding network of the 3(10)-helix, as contrasted with that of the competing alpha-helical structure, is explicitly incorporated into a theoretical model for the 3(10)-helix/alpha-helix equilibrium constant for a given peptide. Finite length effects and the "extra" intrahelical hydrogen bond of the 3(10) form are expressed naturally as a result of this loop analysis. This semiempirical model captures all the established features of existing empirical rules for helical conformational transitions in Aib-rich sequences, as well as the recently detected helical transition induced solely by sequence permutation.     

The structure of the Alzheimer amyloid beta 10-35 peptide probed through replica-exchange molecular dynamics simulations in explicit solvent

The conformational states sampled by the Alzheimer amyloid beta (10-35) (Abeta 10-35) peptide were probed using replica-exchange molecular dynamics (REMD) simulations in explicit solvent. The Abeta 10-35 peptide is a fragment of the full-length Abeta 40/42 peptide that possesses many of the amyloidogenic properties of its full-length counterpart. Under physiological temperature and pressure, our simulations reveal that the Abeta 10-35 peptide does not possess a single unique folded state. Rather, this peptide exists as a mixture of collapsed globular states that remain in rapid dynamic equilibrium with each other. This conformational ensemble is dominated by random coil and bend structures with insignificant presence of an alpha-helical or beta-sheet structure. The 3D structure of Abeta 10-35 is seen to be defined by a salt bridge formed between the side-chains of K28 and D23. This salt bridge is also observed in Abeta fibrils and our simulations suggest that monomeric conformations of Abeta 10-35 contain pre-folded structural motifs that promote rapid aggregation of this peptide.     

Direct observation of protein folding, aggregation, and a prion-like conformational conversion

Protein conformational transition from alpha-helices to beta-sheets precedes aggregation of proteins implicated in many diseases, including Alzheimer and prion diseases. Direct characterization of such transitions is often hindered by the complicated nature of the interaction network among amino acids. A recently engineered small protein-like peptide with a simple amino acid composition features a temperature-driven alpha-helix to beta-sheet conformational change. Here we studied the conformational transition of this peptide by molecular dynamics simulations. We observed a critical temperature, below which the peptide folds into an alpha-helical coiled-coil state and above which the peptide misfolds into beta-rich structures with a high propensity to aggregate. The structures adopted by this peptide during low temperature simulations have a backbone root mean square deviation less than 2 A from the crystal structure. At high temperatures, this peptide adopts an amyloid-like structure, which is mainly composed of coiled anti-parallel beta-sheets with the cross-beta-signature of amyloid fibrils. Most strikingly, we observed conformational conversions in which an alpha-helix is converted into a beta-strand by proximate stable beta-sheets with exposed hydrophobic surfaces and unsaturated hydrogen bonds. Our study suggested a possible generic molecular mechanism of the template-mediated aggregation process, originally proposed by Prusiner (Prusiner, S. B. (1998) Proc. Natl. Acad. Sci. U. S. A. 95, 13363-13383) to account for prion infectivity.     

Conformational investigation of alpha,beta-dehydropeptides. XVI. Beta-turn tendency in Ac-Pro-DeltaXaa-NHMe: crystallographic and theoretical studies.

The crystal structures of two diastereomeric alpha,beta-dehydrobutyrine peptides Ac-Pro-(Z)-DeltaAbu-NHMe (I) and Ac-Pro-(E)-DeltaAbu-NHMe (II) have been determined. Both dehydropeptides adopt betaI-turn conformation characterized by the pairs of (phi(i+1), psi(i+1)) and (phi(i+2), psi(i+2)) angles as -66, -19, -97, 11 degrees for I and -59, -27, -119, 29 degrees for II. In each peptide, the betaI turn is stabilized by (i + 3) --> i intramolecular hydrogen bonds with N...O distance of 3.12 A for I and 2.93 A for II. These structures have been compared to the crystal structures of homologous peptides Ac-Pro-DeltaVal-NHMe and Ac-Pro-DeltaAla-NHMe. Theoretical analyses by DFT/B3LYP/6-31 + G** method of conformers formed by these four peptides and by the saturated peptide Ac-Pro-Ala-NHMe revealed that peptides with a (Z) substituent at the C(beta) (i+2) atom of dehydroamino acid, i.e. Ac-Pro-DeltaVal-NHMe and Ac-Pro-(Z)-DeltaAbu-NHMe, predominantly form beta turns, both in vacuo and in polar environment. The tendency to adopt beta-turn conformation is much weaker for the peptides lacking the (Z) substituent, Ac-Pro-(E)-DeltaAbu-NHMe and Ac-Pro-DeltaAla-NHMe. The latter adopts a semi-extended or an extended conformation in every polar environment, including a weakly polar solvent. The saturated peptide Ac-Pro-Ala-NHMe in vacuo prefers a beta-turn conformation, but in polar environment the differences between various conformers are small. The role of pi-electron correlation and intramolecular hydrogen bonds interaction in stabilizing the hairpin structures are discussed.     

Design of peptides with branched beta-carbon dehydro-residues: syntheses, crystal structures and molecular conformations of two peptides, (I) N-Carbobenzoxy-DeltaVal-Ala-Leu-OCH3 and (II) N-Carbobenzoxy-DeltaIle-Ala-Leu-OCH3.

Highly specific structures can be designed by inserting dehydro-residues into peptide sequences. The conformational preferences of branched beta-carbon residues are known to be different from other residues. As an implication it was expected that the branched beta-carbon dehydro-residues would also induce different conformations when substituted in peptides. So far, the design of peptides with branched beta-carbon dehydro-residues at (i + 1) position has not been reported. It may be recalled that the nonbranched beta-carbon residues induced beta-turn II conformation when placed at (i + 2) position while branched beta-carbon residues induced beta-turn III conformation. However, the conformation of a peptide with a nonbranched beta-carbon residue when placed at (i + 1) position was not found to be unique as it depended on the stereochemical nature of its neighbouring residues. Therefore, in order to induce predictably unique structures with dehydro-residues at (i + 1) position, we have introduced branched beta-carbon dehydro-residues instead of nonbranched beta-carbon residues and synthesized two peptides: (I) N-Carbobenzoxy-DeltaVal-Ala-Leu-OCH3 and (II) N-Carbobenzoxy-DeltaIle-Ala-Leu-OCH3 with DeltaVal and DeltaIle, respectively. The crystal structures of peptides (I) and (II) have been determined and refined to R-factors of 0.065 and 0.063, respectively. The structures of both peptides were essentially similar. Both peptides adopted type II beta-turn conformations with torsion angles; (I): phi1 = -38.7 (4) degrees, psi1 = 126.0 (3) degrees; phi2 = 91.6 (3) degrees, psi2 = -9.5 (4) degrees and (II): phi1 = -37.0 (6) degrees, psi1 = 123.6 (4) degrees, phi2 = 93.4 (4), psi2 = -11.0(4) degrees respectively. Both peptide structures were stabilized by intramolecular 4-->1 hydrogen bonds. The molecular packing in both crystal structures were stabilized in each by two identical hydrogen bonds N1...O1' (-x, y + 1/2, -z) and N2...O2' (-x + 1, y + 1/2, -z) and van der Waals interactions.     

Helix folding simulations with various initial conformations.

Using a solvent-referenced energy calculation, a 16-residue peptide with alanine side chains folded into predominantly alpha-helical conformations during constant temperature (274 K) simulations. From different initial conformations, helical conformations were reached and the multiple energy minima did not become a serious problem. Under the same conditions, the simulation did not indiscriminately fold a sequence such as polyglycine into stable helices. Interesting observations from the simulations relate to the folding mechanism. The electrostatic interactions between the successive amides favored extended conformations (or beta strands) and caused energy barriers to helix folding. beta-bends were observed as intermediates during helix nucleation. The helix propagation toward the C-terminus seemed faster than that toward the N-terminus. In helical conformations, hydrogen bond oscillation between the i,i+ 4 and the i,i+3 patterns was observed. The i,i+3 hydrogen bonds occurred more frequently during helix propagation and deformation near both ends of the helical segment.     

Hydrogen bonding in helical polypeptides from molecular dynamics simulations and amide hydrogen exchange analysis: alamethicin and melittin in methanol.

Molecular dynamics simulations of ion channel peptides alamethicin and melittin, solvated in methanol at 27 degrees C, were run with either regular alpha-helical starting structures (alamethicin, 1 ns; melittin 500 ps either with or without chloride counterions), or with the x-ray crystal coordinates of alamethicin as a starting structure (1 ns). The hydrogen bond patterns and stabilities were characterized by analysis of the dynamics trajectories with specified hydrogen bond angle and distance criteria, and were compared with hydrogen bond patterns and stabilities previously determined from high-resolution NMR structural analysis and amide hydrogen exchange measurements in methanol. The two alamethicin simulations rapidly converged to a persistent hydrogen bond pattern with a high level of 3(10) hydrogen bonding involving the amide NH's of residues 3, 4, 9, 15, and 18. The 3(10) hydrogen bonds stabilizing amide NH's of residues C-terminal to P2 and P14 were previously proposed to explain their high amide exchange stabilities. The absence, or low levels of 3(10) hydrogen bonds at the N-terminus or for A15 NH, respectively, in the melittin simulations, is also consistent with interpretations from amide exchange analysis. Perturbation of helical hydrogen bonding in the residues before P14 (Aib10-P14, alamethicin; T11-P14, melittin) was characterized in both peptides by variable hydrogen bond patterns that included pi and gamma hydrogen bonds. The general agreement in hydrogen bond patterns determined in the simulations and from spectroscopic analysis indicates that with suitable conditions (including solvent composition and counterions where required), local hydrogen-bonded secondary structure in helical peptides may be predicted from dynamics simulations from alpha-helical starting structures. Each peptide, particularly alamethicin, underwent some large amplitude structural fluctuations in which several hydrogen bonds were cooperatively broken. The recovery of the persistent hydrogen bonding patterns after these fluctuations demonstrates the stability of intramolecular hydrogen-bonded secondary structure in methanol (consistent with spectroscopic observations), and is promising for simulations on extended timescales to characterize the nature of the backbone fluctuations that underlie amide exchange from isolated helical polypeptides.     

Molecular modeling of an antigenic complex between a viral peptide and a class I major histocompatibility glycoprotein.

Computer simulation of the conformations of short antigenic peptides (5-10 residues) either free or bound to their receptor, the major histocompatibility complex (MHC)-encoded glycoprotein H-2 Ld, was employed to explain experimentally determined differences in the antigenic activities within a set of related peptides. Starting for each sequence from the most probable conformations disclosed by a pattern-recognition technique, several energy-minimized structures were subjected to molecular dynamics simulations (MD) either in vacuo or solvated by water molecules. Notably, antigenic potencies were found to correlate to the peptides propensity to form and maintain an overall alpha-helical conformation through regular i,i + 4 hydrogen bonds. Accordingly, less active or inactive peptides showed a strong tendency to form i,i + 3 hydrogen bonds at their N-terminal end. Experimental data documented that the C-terminal residue is critical for interaction of the peptide with H-2 Ld. This finding could be satisfactorily explained by a 3-D Q.S.A.R. analysis postulating interactions between ligand and receptor by hydrophobic forces. A 3-D model is proposed for the complex between a high-affinity nonapeptide and the H-2 Ld receptor. First, the H-2 Ld molecule was built from X-ray coordinates of two homologous proteins: HLA-A2 and HLA-Aw68, energy-minimized and studied by MD simulations. With HLA-A2 as template, the only realistic simulation was achieved for a solvated model with minor deviations of the MD mean structure from the X-ray conformation. Water simulation of the H-2 Ld protein in complex with the antigenic nonapeptide was then achieved with the template-derived optimal parameters. The bound peptide retains mainly its alpha-helical conformation and binds to hydrophobic residues of H-2 Ld that correspond to highly polymorphic positions of MHC proteins. The orientation of the nonapeptide in the binding cleft is in accordance with the experimentally determined distribution of its MHC receptor-binding residues (agretope residues). Thus, computer simulation was successfully employed to explain functional data and predicts alpha-helical conformation for the bound peptide.