Temperature coefficients of peptides dissolved in hexafluoroisopropanol monitor distortions of helices

Summary Temperature coefficients are widely used as an indication of solvent accessibility to amide protons. Low temperature coefficients are related to low accessibility and are often interpreted as evidence for intramolecular hydrogen bonding. Conformational shifts, i.e. the difference between chemical shifts of a particular residue in a structured and in a random-coil conformation, provide information on secondary structure. In particular, negative CH^α conformational shifts are often used to delineate the extent of helical stretches. NH conformational shifts show large oscillations within a helix that have been interpreted as the result of helix distortions affecting hydrogen bond lengths. In the ocurse of the study of differnet peptides that adopt a helical structure in the presence of the structure-inducing solvent hexafluoroisopropanol (HFIP), we have found a strong correlation between temperature coefficients and amide conformational shifts. However, contrary to the initial expectations, lower temperature coefficients were associated to amide protons involved in longer, and presumably weaker, hydrogen bonds. The correlation can be explained, however, assuming that, in helical peptides dissolved in HFIP, temperature affects the chemical shift of amide protons mainly by changing the average length of intramolecular hydrogen bonds and changes in solvent accessibility play only a secondary role under these experimental conditions. The pattern of temperature coefficients in helical peptides can therefore be used to identify short or long hydragen bonds causing bending of the helix axis.     

310-Helices,Helix Screw Sense and Screw Sense Reversal in the Dehydro-peptide Boc-Val-ΔPhe-Gly-ΔPhe-Val-OMe

The pentapeptide Boc-Val-ΔPhe-Gly-ΔPhe-Val-OMe, containing two dehydro-phenylalanine (ΔPhe) residues, has been synthesized and its structure investigated. In the crystalline state, the molecule adopts a right-handed 3₁₀-helical conformation stabilized by two intramolecular hydrogen bonds between CO of Val¹ and NH of ΔPhe⁴, and between CO of ΔPhe² and NH of Val⁵, respectively. NMR measurements are consistent with the presence of 3₁₀-helical structures also in acetonitrile and dimethylsulphoxide solution: the distances between backbone protons estimated from NOE connectivities are in overall agreement with those observed in the solid state; the chemical shifts of the amide protons show the smaller temperature coefficients for the NHs that in solid state are involved in intramolecular hydrogen bonds. The CD spectra in acetonitrile, chloroform, methanol and dimethylsulphoxide display exciton couplets of bands corresponding to the ΔPhe electronic transition at 280nm; the sign of the bands is consistent with the presence of helical structures having a prevalent left-handed screw sense. Addition of 1,1,1,3,3,3-hexafluoro- propan-2-ol gives rise to the gradual appearance of a couplet of opposite sign, suggesting the helix reversal from left-handed sense to right-handed sense. The conformational behaviour is discussed on the basis of the specific sequence of the peptide.     

Temperature coefficients of amide proton NMR resonance frequencies in trifluoroethanol: A monitor of intramolecular hydrogen bonds in helical peptides?

Summary 2D ¹H NMR spectroscopy of two -helical peptides which differ in their amphipathicity has been used to investigate the relationships between amide-proton chemical shifts, amide-proton exchange rates, temperature, and trifluoroethanol (TFE) concentration. In 50% TFE, in which the peptides are maximally helical, the amide-proton chemical shift and temperature coefficient patterns are very similar to each other in each peptide. Temperature coefficients from –10 to –6 ppb/K, usually indicative of the lack of intramolecular hydrogen bonds, were observed even for hydrophobic amino acids in the center of the -helices. However, slow hydrogen isotope exchange for residues from 4 to 16 in both 18-mer helices indicates intact intramolecular hydrogen bonds over most of the length of these peptides. Based on these anomalous observations, we suggest that the pattern of amide-proton shifts in -helices in H₂O/TFE solvents is dominated by bifurcated intermolecular hydrogen-bond formation between the backbone carbonyl groups and TFE. The amide-proton chemical shift changes with increasing temperature may be interpreted by a disruption of intermolecular hydrogen bonds between carbonyl groups and the TFE in TFE/water rather than by the length of intramolecular hydrogen bonds in -helices. Supplementary Material is available upon request, comprising seven pages with listings of experimental details and the NMR shift data for the two peptides.     

Comparative (1)H NMR and molecular modeling study of hydroxy protons of beta-D-Galp-(1-->4)-beta-D-GlcpNAc-(1-->2)-alpha-D-Manp-(1-->O)(CH(2))(7)CH(3) analogues in aqueous solution

The (1)H chemical shifts, coupling constants, temperature coefficients, exchange rates, and inter-residual ROEs have been measured, in aqueous solution, for the hydroxy and amine/amide proton resonances of a set of beta-D-Galp-(1-->4)-beta-D-GlcpNAc-(1-->2)-alpha-D-Manp-(1-->O)(CH(2))(7)CH(3) analogues. From the structural data, a few significant structural features could be ascertained, such as a preferential anti-conformation for the amide protons of the N-acetyl and N-propionyl groups. The introduction of systematic modifications at Gal 2-C and Gal 6-C resulted in alterations of the Gal 4-OH, Gal 3-OH, and GlcNAc 3-OH areas, since variations in chemical shifts and temperature coefficient were observed. In order to verify the possibility of hydrogen bonds, molecular dynamics simulations in the gas phase and explicit solvent were performed and correlated with the experimental data. A network of hydrogen bonds to solvent molecules was observed, but no strong intramolecular hydrogen bonding was observed.     

Immunogenic peptides corresponding to the dominant antigenic region alanine-597 to cysteine-619 in the transmembrane protein of simian immunodeficiency virus have a propensity to fold in aqueous solution.

Two synthetic peptides corresponding to the N- and C-terminal halves of a 23 amino acid sequence representing an immunodominant domain of the simian immunodeficiency virus of macaque origin (SIVmac) were examined for conformational preferences in aqueous solution by proton nuclear magnetic resonance methods. The two constituent peptides, termed A12-7 (Ala597-Ile-Glu-Lys-Tyr-Leu-Glu-Asp-Gln-Ala-Gln607) and A12-9 (Leu608-Asn-Ala-Trp-Gly-Cys-Ala-Phe-Arg-Gln-Val-Ser619), were found to contain a considerable conformational preference for states in which the backbone phi and psi angles populate the alpha region of the Ramachandran plot. Further, for peptide A12-9, the types and intensities of the nuclear Overhauser effect (NOE) connectivities between protons in the polypeptide backbone suggest that these states appear to include helical turns. The temperature dependence of the amide proton chemical shifts indicates that some degree of intramolecular hydrogen bonding occurs in these peptides. These results are consistent with a model in which immunogenic peptides which induce antibodies reactive with the intact protein from which the peptide sequence was derived contain conformational preferences in water solution for states other than the extended-chain forms typically found in "random coil" peptides.     

Amide proton temperature coefficients as hydrogen bond indicators in proteins

Correlations between amide proton temperature coefficients (_HN/T) and hydrogen bonds were investigated for a data set of 793 amides derived from 14 proteins. For amide protons showing temperature gradients more positive than –4.6 ppb/K there is a hydrogen bond predictivity value exceeding 85%. It increases to over 93% for amides within the range between –4 and –1 ppb/K. Detailed analysis shows an inverse proportionality between amide proton temperature coefficients and hydrogen bond lengths. Furthermore, for hydrogen bonds of similar bond lengths, values of temperature gradients in -helices are on average 1 ppb/K more negative than in -sheets. In consequence, a number of amide protons in -helices involved in hydrogen bonds shorter than 2 Å show _HN/T < –4.6 ppb/K. Due to longer hydrogen bonds, 90% of amides in 3₁₀ helices and 98% in -turns have temperature coefficients more positive than –4.6ppb/K. Ring current effect also significantly influences temperature coefficients of amide protons. In seven out of eight cases non-hydrogen bonded amides strongly deshielded by neighboring aromatic rings show temperature coefficients more positive than –2 ppb/K. In general, amide proton temperature gradients do not change with pH unless they correspond to conformational changes. Three examples of pH dependent equilibrium showing hydrogen bond formation at higher pH were found. In conclusion, amide proton temperature coefficients offer an attractive and simple way to confirm existence of hydrogen bonds in NMR determined structures.     

400 MHz NMR study on the C-terminal fragment 21–28 of vasoactive intestinal peptide

A conformational study of the C-terminal fragment 21–28 of vasoactive intestinal peptide (VIP) was carried out by high resolution NMR spectroscopy. All spectral data were recorded with a Brüker 400 MHz spectrometer. The correct assignment of peaks was determined by specific homonuclear decoupling and by titration. The chemical shifts of amide protons were measured as a function of temperature in order to detect the presence of intramolecular hydrogen bonds. A tridimensional structure is proposed for the peptide.     

The hydration of amides in helices; a comprehensive picture from molecular dynamics,IR, and NMR

We examined the hydration of amides of alpha(3)D, a simple, designed three-helix bundle protein. Molecular dynamics calculations show that the amide carbonyls on the surface of the protein tilt away from the helical axis to interact with solvent water, resulting in a lengthening of the hydrogen bonds on this face of the helix. Water molecules are bonded to these carbonyl groups with partial occupancy ( approximately 50%-70%), and their interaction geometries show a large variation in their hydrogen bond lengths and angles on the nsec time scale. This heterogeneity is reflected in the carbonyl stretching vibration (amide I' band) of a group of surface Ala residues. The surface-exposed amides are broad, and shift to lower frequency (reflecting strengthening of the hydrogen bonds) as the temperature is decreased. By contrast, the amide I' bands of the buried (13)C-labeled Leu residues are significantly sharper and their frequencies are consistent with the formation of strong hydrogen bonds, independent of temperature. The rates of hydrogen-deuterium exchange and the proton NMR chemical shifts of the helical amide groups also depend on environment. The partial occupancy of the hydration sites on the surface of helices suggests that the interaction is relatively weak, on the order of thermal energy at room temperature. One unexpected feature that emerged from the dynamics calculations was that a Thr side chain subtly disrupted the helical geometry 4-7 residues N-terminal in sequence, which was reflected in the proton chemical shifts and the rates of amide proton exchange for several amides that engage in a mixed 3(10)/alpha/pi-helical conformation.     

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.     

Structure and dynamics of the neutrophil defensins NP-2, NP-5, and HNP-1: NMR studies of amide hydrogen exchange kinetics

The exchange kinetics for the slowly exchanging amide hydrogens in three defensins, rabbit NP-2, rabbit NP-5, and human HNP-1, have been measured over a range of pH at 25°C using 1D and 2D NMR methods. These NHs have exchange rates 10² to 10⁵ times slower than rates from unstructured model peptides. The observed distribution of exchange rates under these conditions can be rationalized by intramolecular hydrogen bonding of the individual NHs, solvent accessibility of the NHs, and local fluctuations in structure. The temperature dependencies of NH chemical shifts (NH temperature coefficients) were measured for the defensins and these values are consistent with the defensin structure. A comparison is made between NH exchange kinetics, NH solvent accessibility, and NH temperature coefficients of the defensins and other globular proteins. Titration of the histidine side chain in NP-2 was examined and the results are mapped to the three-dimensional structure. © 1994 Wiley-Liss, Inc.     

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.     

The conformational characteristics in solution of the synthetic cyclic hexapeptide 5L-ala-D-ala

An attempt to elucidate the solution conformation(s) of the synthetic cyclic hexapeptide 5L -ala·D-ala is described. Nuclear magnetic resonance (nmr) spectra are recorded for the purpose of measuring the vicinal coupling constant between the amide and α-protons in each residue and to observe the deuterium exchange rate and temperature dependence of the chemical shift of each amide proton. Low-energy cyclic conformations, whose individual residues are in conformations consistent with the observed amide to α-proton coupling constant, are searched for in an approximate theoretical treatment. The two lowest energy, all trans peptide bond conformations generated are distinguishable by the presence or absence of a single intramolecular hydrogen bond. The observed temperature independence of the chemical shift of one of the amide protons is consistent with the presence of a single intramolecular hydrogen bond, while the observation of similar deuterium exchange rates for each of the amide protons indicates their comparable availability to solvent. Consequently, it is concluded that 5L -ala·D-ala is in rapid equilibrium between conformations with and without a single internal hydrogen bond and possesses considerable conformational flexibility in solution.     

Synthetic amatoxin analogues. II. A proton n.m.r. study of S-deoxo-Ile3-(L)Ala5 and S-deoxo-Ile3-(D)Ala5-amaninamide

Amatoxin analogues with D and L-Ala substitutions in position 5 have been studied by means of 1- and 2-dimensional n.m.r. spectroscopy at 500 MHz. The assignment of all resonances for both analogues has been carried out mostly with the use of COSY and NOESY type experiments. Temperature coefficients for the amide NH protons have been measured and the data compared to known amatoxin structures. The results obtained demonstrate that the rigidity of the bicyclic amatoxin framework is preserved in the D and L-Ala5 analogues, although the temperature coefficients point to intramolecular hydrogen bonds stronger in the case of the L-Ala analogue. The 10-fold decrease of biological activity is discussed in terms of structural features involving also the Trp4 indole accessibility.     

Temperature-dependence of protein hydrogen bond properties as studied by high-resolution NMR

The temperature-dependence of a large number of NMR parameters describing hydrogen bond properties in the protein ubiquitin was followed over a range from 5 to 65 degrees C. The parameters comprise hydrogen bond (H-bond) scalar couplings, h3JNC', chemical shifts, amide proton exchange rates, 15N relaxation parameters as well as covalent 1JNC' and 1JNH couplings. A global weakening of the h3JNC' coupling with increasing temperature is accompanied by a global upfield shift of the amide protons and a decrease of the sequential 1JNC' couplings. If interpreted as a linear increase of the N...O distance, the change in h3JNC' corresponds to an average linear thermal expansion coefficient for the NH-->O hydrogen bonds of 1.7 x 10(-4)/K, which is in good agreement with overall volume expansion coefficients observed for proteins. A residue-specific analysis reveals that not all hydrogen bonds are affected to the same extent by the thermal expansion. The end of beta-sheet beta1/beta5 at hydrogen bond E64-->Q2 appears as the most thermolabile, whereas the adjacent hydrogen bond I3-->L15 connecting beta-strands beta1 and beta2 is even stabilized slightly at higher temperatures. Additional evidence for the stabilization of the beta1/beta2 beta-hairpin at higher temperatures is found in reduced hydrogen exchange rates for strand end residue V17. This reduction corresponds to a stabilizing change in free energy of 9.7 kJ/mol for the beta1/beta2 hairpin. The result can be linked to the finding that the beta1/beta2 hairpin behaves as an autonomously folding unit in the A-state of ubiquitin under changed solvent conditions. For several amide groups the temperature-dependencies of the amide exchange rates and H-bond scalar couplings are uncorrelated. Therefore, amide exchange rates are not a sole function of the hydrogen bond "strength" as given by the electronic overlap of donors and acceptors, but are clearly dependent on other blocking mechanisms.     

pH dependence of hydrogen exchange from backbone peptide amides in apamin

The kinetics of hydrogen exchange of the 11 most protected backbone amides of bee venom apamin have been measured between pH 1 and pH 8.5 by using time-resolved and saturation-transfer NMR spectroscopy. The five amides most protected from base-catalyzed exchange, those of residues 5 and 12-15, show highly correlated exchange behavior in the base-catalyzed regime. It is proposed that the intramolecular hydrogen bonds stabilizing these amides define a stable cooperative unit of secondary structure in apamin (a C-terminal helix and an N-terminal beta-turn). This conformational unit is further stabilized (by 5-6 kJ mol-1) on titration of the Glu-7 side-chain carboxyl group. The relative contributions of specific intramolecular interactions to this conformational stabilization are estimated. The pHminima in the pH-dependent single amide exchange curves are compared with values predicted by correcting for sequence-dependent contributions to amide exchange rates [Molday, R. S., Englander, S. W., & Kallen, R. G. (1972) Biochemistry 11, 150-158]. The lack of correlation suggests that the "open" conformers from which amide exchange occurs are nonrandom. This conclusion is dependent on the assumption that acid-catalyzed exchange occurs via N-protonation so that residual conformational effects on exchange rates in the open conformers will affect acid- and base-catalyzed rates in approximately equal and opposite ways. A strong correlation between the measured pHminima and the amide proton chemical shifts is observed, however, and this may be most easily accommodated if acid-catalyzed exchange occurs by the imidic acid mechanism (via amide O-protonation).     

Solution conformation of a cyclic neurokinin antagonist: A nmr and molecular dynamics study

The solution structure of a hexapeptide, cyclo(Gln-Trp-Phe-Gly-Leu-Met), which is a selective NK-2 antagonist, has been studied by a combination of two-dimensional nmr and molecular dynamics (MD) techniques. The simulation based on nmr and MD data resulted in the convergence to a family of structures. Free molecular dynamics for 50 ps in the presence of DMSO solvent molecules shows that the structure is energetically stable. One intramolecular hydrogen bond between the amide proton of Gin and the carbonyl oxygen of Gly was revealed. This result is consistent with the results from the measurement of the temperature coefficient of the amide protons. The extent of intermolecular hydrogen bonding between the amide protons of the peptide and DMSO was also revealed by the free MD simulation. The resulting structure of the cyclic peptide contains a variation type I′ β-turn in the Gly-Leu-Met-Gln segment. Comparison of the structure of this peptide with that of other NK-2 antagonist cyclic hexapeptides was made, and the activity of cyclic antagonists appears to be inversely related to the conformational rigidity of the cyclic peptides. © 1994 John Wiley & Sons, Inc.     

Solvent accessibility in folded proteins. Studies of hydrogen exchange in trypsin.

In a native protein, the exchange of a peptide amide proton with solvent occurs by one of two pathways, either directly from the folded protein, or via unfolding, exchange taking place from the unfolded protein. From the thermal unfolding rate constants, the contribution of unfolding to the over-all kinetics as a function of solvent and temperature has been determined. Exchange involving unfolding of the protein is characterized by a high activation energy, in the range of 50 to 60 Cal per mol. The activiation energy (Eapp) of the rates of exchange directly from the folded protein is approximately 20 to 25 Cal per mol. Because for the proton transfer step, Eapp approximately equal to 20 Cal per mol, the activation energy for any contributing protein conformational process(es) is approximately equal to 0 to 5 Cal per mol. Most, if not all, of the peptide amide protons in a folded protein can exchange directly with solvent without the protein unfolding. The number of "slowly" exchanging protons at a given condition of pH and temperature is not related to a discrete structural unit, but rather to the distribution of observed rates within the broader distribution of actual rates. The large attenuation of hydrogen exchange rates in folded proteins, resulting in a distribution of first order rates over 6 orders of magnitude, is primarily due to the effects of restricted solvent accessibility of labile protons in the three-dimensional structure. Any protein conformational process, such as protein fluctuations, invoked to explain the solvent accessibility must be of low activation energy and attenuated by ethanol and other co-solvents (Woodward, C. K., Ellis, L. M., and Rosenberg, A. (1974) J. Biol. Chem. 250, 440-444).     

Nuclear magnetic resonance studies of the internal dynamics in Apo, (Cd2+)1 and (Ca2+)2 calbindin D9k. The rates of amide proton exchange with solvent.

The backbone dynamics of the EF-hand Ca(2+)-binding protein, calbindin D9k, has been investigated in the apo, (Cd2+)1 and (Ca2+)2 states by measuring the rate constants for amide proton exchange with solvent. 15N-1H correlation spectroscopy was utilized to follow direct 1H-->2H exchange of the slowly exchanging amide protons and to follow indirect proton exchange via saturation transfer from water to the rapidly exchanging amide protons. Plots of experimental rate constants versus intrinsic rate constants have been analyzed to give qualitative insight into the opening modes of the protein that lead to exchange. These results have been interpreted within the context of a progressive unfolding model, wherein hydrophobic interactions and metal chelation serve to anchor portions of the protein, thereby damping fluctuations and retarding amide proton exchange. The addition of Ca2+ or Cd2+ was found to retard the exchange of many amide protons observed to be in hydrogen-bonding environments in the crystal structure of the (Ca2+)2 state, but not of those amide protons that were not involved in hydrogen bonds. The largest changes in rate constant occur for residues in the ion-binding loops, with substantial effects also found for the adjacent residues in helices I, II and III, but not helix IV. The results are consistent with a reorganization of the hydrogen-bonding networks in the metal ion-binding loops, accompanied by a change in the conformation of helix IV, as metal ions are chelated. Further analysis of the results obtained for the three states of metal occupancy provides insight into the nature of the changes in conformational fluctuations induced by ion binding.     

Periodic properties of proton conformational shifts in isolated protein helices. An experimental study.

In this work, the helix-forming residues in fragments of several proteins (ribonuclease, thermolysin, tendamistat and angiogenin) were identified by NOE and the helix proton shifts were measured as delta changes associated with helix-population increments driven by trifluoroethanol addition. When estimated in this way, a regular pattern of helix conformational shifts was clearly seen in the delta delta versus sequence profiles of all the peptides studied. The helix periodicity of the H alpha and H beta resonances was especially clear, an observation that earlier statistical studies of protein delta values failed to predict. Amide protons showed the largest helix shifts, but with a less-sharply defined periodic character. Aromatic residues considerably distorted the periodicity of the helix amide shifts in some peptides, as evidenced by the delta shifts of a RNase A fragment 1-15 analog in which the two aromatic residues were replaced by Ala. The relationship between helix periodicity and peptide amphiphatic character is discussed.     

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.