Alanine-based Peptides Research Articles

Fast folding dynamics of α-helical peptides – Effect of solvent additives and pH

The fast folding/unfolding dynamics of an alanine-based α-helical model peptide was investigated using nanosecond temperature jump experiments in D 2O. Temperature jumps were induced either by indirect heating, using visible laser pulses and a heat-transducing triphenylmethane dye, or by direct heating, using IR laser pulses. Although direct heating improves the quality, reliability and speed of the experiments, the heat-transducing dye was shown to not affect the helix-coil dynamics of alanine-based model peptides, thus validating the indirect heating method for this class of peptides. On the other hand, the presence of high concentrations of salt ions was found to alter the helix folding dynamics by screening of charged residues. We conclude that electrostatic interactions between residues have significant effects on the dynamics of helix folding. Similarly, the solvent pH was observed to affect the dynamics of helix folding, even in a range where no protonation/deprotonation is expected to occur.

Molecular Dynamics Simulation of Folding of a Short Helical Peptide with Many Charged Residues

A molecular dynamics simulation of the folding of conantokin-T (con-T), a short helical peptide with 5 helical turns of 21 amino acids with 10 charged residues, was carried out to examine folding pathways for this peptide and to predict the folding rate. In the 18 trajectories run at 300 K, 16 trajectories folded, with an averaged folding time of approximately 50 ns. Two trajectories did not fold in up to 200 ns simulation. The folded structure in folded trajectories is in good agreement with experimental structure. An analysis of the trajectories showed that, at the beginning of a few nanoseconds, helix formation started from residues 5-9 with assistance of a hydrophobic clustering involving Tyr5, Met8, and Leu9. The peptide formed a U-shape mainly due to charge-charge interactions between charged residues at the N- and C-terminus segments. In the next approximately 10 ns, several nonnative charge-charge interactions were broken and nonnative Gla10-Lys18 (this denotes a salt bridge between Gal10 and Lys18) and/or Gla10-Lys19 interactions appeared more frequently in this folding step and the peptide became a fishhook J-shape. From this structure, the peptide folded to the folded state in 7 of all 16 folded trajectories in approximately 15 ns. Alternatively, in approximately 30 ns, the con-T went to a conformation in an L-shape with 4 helical turns and a kink at the Arg13 and Gla14 segment in the other 9 trajectories. Con-T in the L-shape then required another approximately 15 ns to fold into the folded state. In addition, in overall folding times, the former 7 trajectories folded faster with the total folding times all shorter than 45 ns, while the latter 9 trajectories folded at a time longer than 45 ns, resulting in an average folding time of approximately 50 ns. Two major folding intermediates found in 2 nonfolded trajectories are stabilized by charge clusters of 5 and 6 charged residues, respectively. With inclusion of friction and solvent-solvent interactions, which were ignored in the present GB/SA solvation model, the folding time obtained above should be multiplied by a factor of 1.25-1.7 according to a previous, similar simulation study. This results in a folding time of 65-105 ns, slightly shorter than the folding time of 127 ns for an alanine-based peptide of the same length. This suggests that the energy barrier of folding for this type of peptides with many charged residues is slightly lower than alanine-based helical peptides by less than 1 kcal/mol.

Unusual compactness of a polyproline type II structure

Polyproline type II (PPII) helix has emerged recently as the dominant paradigm for describing the conformation of unfolded polypeptides. However, most experimental observables used to characterize unfolded proteins typically provide only short-range, sequence-local structural information that is both time- and ensemble-averaged, giving limited detail about the long-range structure of the chain. Here, we report a study of a long-range property: the radius of gyration of an alanine-based peptide, Ace-(diaminobutyric acid)2-(Ala)7-(ornithine)2-NH2. This molecule has previously been studied as a model for the unfolded state of proteins under folding conditions and is believed to adopt a PPII fold based on short-range techniques such as NMR and CD. By using synchrotron radiation and small-angle x-ray scattering, we have determined the radius of gyration of this peptide to be 7.4 +/- 0.5 angstroms, which is significantly less than the value expected from an ideal PPII helix in solution (13.1 angstroms). To further study this contradiction, we have used molecular dynamics simulations using six variants of the AMBER force field and the GROMOS 53A6 force field. However, in all cases, the simulated ensembles underestimate the PPII content while overestimating the experimental radius of gyration. The conformational model that we propose, based on our small angle x-ray scattering results and what is known about this molecule from before, is that of a very flexible, fluctuating structure that on the level of individual residues explores a wide basin around the ideal PPII geometry but is never, or only rarely, in the ideal extended PPII helical conformation.

Interactions of Peptides with a Protein Pore

The partitioning of polypeptides into nanoscale transmembrane pores is of fundamental importance in biology. Examples include protein translocation in the endoplasmic reticulum and the passage of proteins through the nuclear pore complex. Here we examine the exchange of cationic α-helical peptides between the bulk aqueous phase and the transmembrane β-barrel of the α-hemolysin ( αHL) protein pore at the single-molecule level. The experimental kinetic data suggest a two-barrier, single-well free energy profile for peptide transit through the αHL pore. This free energy profile is strongly voltage- and peptide-length-dependent. We used the Woodhull-Eyring formalism to rationalize the values measured for the association and dissociation rate constants k on and k off, and to separate k off into individual rate constants for exit through each of the openings of the protein pore. The rate constants k on, k off cis , and k off trans decreased with the length of the peptide. At high transmembrane potentials, the alanine-based peptides, which include bulky lysine side chains, bind more strongly (formation constants K f ∼ tens of M −1) than highly flexible polyethylene glycols ( K f ∼ M −1) to the lumen of the αHL protein pore. In contrast, at zero transmembrane potential, the peptides bind weakly to the lumen of the pore, and the affinity decreases with the peptide length, similar to the case of the polyethylene glycols. The binding is enhanced at increased transmembrane potentials, because the free energy contribution Δ G = − ζδFV/R T predominates with the peptides. We suggest that the αHL protein may serve as a robust and versatile model for examining the interactions between positively charged signal peptides and a β-barrel pore.

Left-Handed and Ambidextrous Helices in the Gas Phase

We have used ion mobility mass spectrometry to study the effect of d-residues on helix formation in unsolvated alanine-based peptides. The right-handed helix of AC-A15K + H+ is significantly disrupted when five or more of the natural L-residues are randomly replaced with D-residues. On the other hand, when a block of L-residues is replaced with D-residues, an unusual ambidextrous structure with helical segments of opposite chirality is formed. A peptide with all D-residues forms a left-handed helix.

Monte Carlo Studies of Folding, Dynamics, and Stability in α-Helices

Folding simulations of polyalanine peptides were carried out using an off-lattice Monte Carlo simulation technique. The peptide was represented as a chain of residues, each of which contains two interaction sites: one corresponding to the C α atom and the other to the side chain. A statistical potential was used to describe the interaction between these sites. The preferred conformations of the peptide chain on the energy surface, starting from several initial conditions, were searched by perturbations on its generalized coordinates with the Metropolis criterion. We observed that, at low temperatures, the effective energy was low and the helix content high. The calculated helix propagation ( s) and nucleation ( σ) parameters of the Zimm-Bragg model were in reasonable agreement with the empirical data. Exploration of the energy surface of the alanine-based peptides (AAQAA) 3 and AAAAA(AAARA) 3A demonstrated that their behavior is similar to that of polyalanine, in regard to their effective energy, helix content, and the temperature-dependence of their helicity. In contrast, stable secondary structures were not observed for (Gly) 20 at similar temperatures, which is consistent with the nonfolder nature of this peptide. The fluctuations in the slowest dynamics mode, which describe the elastic behavior of the chain, showed that as the temperature decreases, the polyalanine peptides become stiffer and retain conformations with higher helix content. Clustering of conformations during the folding phase implied that polyalanine folds into a helix through fewer numbers of intermediate conformations as the temperature decreases.

Pairwise Coupling in an Arg-Phe-Met Triplet Stabilizes α-Helical Peptide via Shared Rotamer Preferences

The hydrophobic Arg-Phe and Phe-Met side chain interactions stabilize the alpha-helix by -0.29 and -0.59 kcal/mol, respectively, when placed i, i + 4 in an alanine-based peptide. When both interactions are present simultaneously, however, they stabilize the helix by an additional -0.75 kcal/mol, nearly as much as the sum of its parts. We attribute this coupling to a shared rotamer preference, as the central Phe is t in both bonds. The energetic cost of restricting the Phe residue into a t conformation is only paid once in the triplet, rather than twice when the interactions are separate. Coupling is thus demonstrated to have large effects on protein stability.

Entropic Stabilization of Isolated β-Sheets

Temperature-dependent electric deflection measurements have been performed for a series of unsolvated alanine-based peptides (Ac-WA(n)-NH(2), where Ac = acetyl, W = tryptophan, A = alanine, and n = 3, 5, 10, 13, and 15). The measurements are interpreted using Monte Carlo simulations performed with a parallel tempering algorithm. Despite alanine's high helix propensity in solution, the results suggest that unsolvated Ac-WA(n)-NH(2) peptides with n > 10 adopt beta-sheet conformations at room temperature. Previous studies have shown that protonated alanine-based peptides adopt helical or globular conformations in the gas phase, depending on the location of the charge. Thus, the charge more than anything else controls the structure.

Non-Covalent Interactions between Unsolvated Peptides: Helical Complexes Based on Acid−Base Interactions

We have used ion-mobility mass spectrometry to examine the conformations of the protonated complex formed between AcA(7)KA(6)KK and AcEA(7)EA(7), helical alanine-based peptides that incorporate glutamic acid (E) and lysine (K). Designed interactions between the acidic E and basic K residues help to stabilize the complex, which is generated by electrospray and studied in the gas phase. There are two main conformations: (1) a coaxial linear arrangement where the helices are tethered together by an EKK interaction between the pair of lysines at the C-terminus of the AcA(7)KA(6)KK peptide and a glutamic acid at the N-terminus of the AcEA(7)EA(7) peptide and (2) a coiled-coil arrangement with side-by-side antiparallel helices where there is an additional EK interaction between the E and K residues in the middle of the helices. The coiled-coil opens up to the coaxial linear structure as the temperature is raised. Entropy and enthalpy changes for the opening of the coiled-coil were derived from the measurements. The enthalpy change indicates that the interaction between the E and K residues in the middle of the helices is a weak neutral hydrogen bond. The EKK interaction is significantly stronger.

Comparative Study of Generalized Born Models: Born Radii and Peptide Folding

In this study, we have implemented four analytical generalized Born (GB) models and investigated their performance in conjunction with the GROMOS96 force field. The four models include that of Still and co-workers, the HCT model of Cramer, Truhlar, and co-workers, a modified form of the AGB model of Levy and co-workers, and the GBMV2 model of Brooks and co-workers. The models were coded independently and implemented in the GROMOS software package and in TINKER. They were compared in terms of their ability to reproduce the results of Poisson-Boltzmann (PB) calculations and in their performance in the ab initio peptide folding of two peptides, one that forms a beta-hairpin in solution and one that forms an alpha-helix. In agreement with previous work, the GBMV2 model is most successful in reproducing PB results while the other models tend to underestimate the effective Born radii of buried atoms. In contrast, stochastic dynamics simulations on the folding of the two peptides, the C-terminus beta-hairpin of the B1 domain of protein G and the alanine-based alpha-helical peptide 3K(I), suggest that the simpler GB models are more effective in sampling conformational space. Indeed, the Still model used in conjunction with the GROMOS96 force field is able to fold the hairpin peptide to a native-like structure without the benefit of enhanced sampling techniques. This is due in part to the properties of the united-atom GROMOS96 force field which appears to be more flexible, and hence to sample more efficiently, than force fields such as OPLSAA. Our results suggest a general strategy which involves using different combinations of force fields and solvent models in different applications, for example, using GROMOS96 and a simple GB model in sampling and OPLSAA and a more accurate GB model in refinement. The fact that various methods have been implemented in a unified way should facilitate the testing and subsequent use of different methods to evaluate conformational free energies in different applications. Our results also bear on some general issues involved in peptide folding and structure prediction which are addressed in the Discussion.

The effects of individual amino acids on the fast folding dynamics of α-helical peptides

Nanosecond temperature jump experiments coupled to time-resolved infrared spectroscopy were carried out on a series of alanine-based peptides containing different guest amino acids to study the effects of residues with different helix propensities on the helix-coil dynamics.

Dissecting Contributions to the Denaturant Sensitivities of Proteins

Understanding the molecular basis for protein denaturation by urea and guanidinium chloride (GdmCl) should accommodate the observation that, on a molar basis, GdmCl is generally 2-2.5-fold more effective as a protein denaturant than urea. Previous studies [Smith, J. S., and Scholtz, J. M. (1996) Biochemistry 35, 7292-7297] have suggested that the effects of GdmCl on the stability of alanine-based helical peptides can be separated into denaturant and salt effects, since adding equimolar NaCl to urea enhanced urea-induced unfolding to an extent that was close to that of Gdm. We reinvestigated this observation using an alanine-based helical peptide (alahel) that lacks side chain electrostatic contributions to stability, and compared the relative denaturant sensitivities of this peptide with that of tryptophan zipper peptides (trpzip) whose native conformations are stabilized largely by cross-strand indole ring interactions. In contrast to the observations of Smith and Scholtz, GdmCl was only slightly more powerful as a denaturant of alahel than urea in salt-free buffer (the denaturant m value m(GdmCl)/m(urea) ratio = 1.4), and the denaturation of alahel by urea exhibited only a small dependence on NaCl or KCl. The trpzip peptides were much more sensitive to GdmCl than to urea (m(GdmCl)/m(urea) = 3.5-4). These observations indicate that the m(GdmCl)/m(urea) ratio of 2-2.5 for proteins results from a combination of effects on the multiple contributions to protein stability, for which GdmCl may be only slightly more effective than urea (e.g., hydrogen bonds) or considerably more effective than urea (e.g., indole-indole interactions).

Cutoff Size Need Not Strongly Influence Molecular Dynamics Results for Solvated Polypeptides

The correct treatment of van der Waals and electrostatic nonbonded interactions in molecular force fields is essential for performing realistic molecular dynamics (MD) simulations of solvated polypeptides. The most computationally tractable treatment of nonbonded interactions in MD utilizes a spherical distance cutoff (typically, 8-12 A) to reduce the number of pairwise interactions. In this work, we assess three spherical atom-based cutoff approaches for use with all-atom explicit solvent MD: abrupt truncation, a CHARMM-style electrostatic shift truncation, and our own force-shifted truncation. The chosen system for this study is an end-capped 17-residue alanine-based alpha-helical peptide, selected because of its use in previous computational and experimental studies. We compare the time-averaged helical content calculated from these MD trajectories with experiment. We also examine the effect of varying the cutoff treatment and distance on energy conservation. We find that the abrupt truncation approach is pathological in its inability to conserve energy. The CHARMM-style shift truncation performs quite well but suffers from energetic instability. On the other hand, the force-shifted spherical cutoff method conserves energy, correctly predicts the experimental helical content, and shows convergence in simulation statistics as the cutoff is increased. This work demonstrates that by using proper and rigorous techniques, it is possible to correctly model polypeptide dynamics in solution with a spherical cutoff. The inherent computational advantage of spherical cutoffs over Ewald summation (and related) techniques is essential in accessing longer MD time scales.

Improved prediction for N-termini of alpha-helices using empirical information.

The prediction of the secondary structure of proteins from their amino acid sequences remains a key component of many approaches to the protein folding problem. The most abundant form of regular secondary structure in proteins is the alpha-helix, in which specific residue preferences exist at the N-terminal locations. Propensities derived from these observed amino acid frequencies in the Protein Data Bank (PDB) database correlate well with experimental free energies measured for residues at different N-terminal positions in alanine-based peptides. We report a novel method to exploit this data to improve protein secondary structure prediction through identification of the correct N-terminal sequences in alpha-helices, based on existing popular methods for secondary structure prediction. With this algorithm, the number of correctly predicted alpha-helix start positions was improved from 30% to 38%, while the overall prediction accuracy (Q3) remained the same, using cross-validated testing. Although the algorithm was developed and tested on multiple sequence alignment-based secondary structure predictions, it was also able to improve the predictions of start locations by methods that use single sequences to make their predictions. Furthermore, the residue frequencies at N-terminal positions of the improved predictions better reflect those seen at the N-terminal positions of alpha-helices in proteins. This has implications for areas such as comparative modeling, where a more accurate prediction of the N-terminal regions of alpha-helices should benefit attempts to model adjacent loop regions. The algorithm is available as a Web tool, located at http://rocky.bms.umist.ac.uk/elephant.

Fs-21 Peptides Can Form Both Single Helix and Helix−Turn−Helix

Detailed folding processes and mechanisms of two alanine-based peptides (Fs-21 and MABA-Fs) were investigated by all atom molecular dynamic simulation with a new AMBER force field and Generalized Born continuum solvent model. Both peptides showed multiphase folding processes. Much like what has been envisaged by the folding funnel theory, the number of accessible conformations descended quickly as folding progresses. Interestingly, MABA-Fs and Fs-21 peptides exhibited notably different folding kinetics; the Fs-21-folding was a two-phase process while MABA-Fs went through three phases and folded more slowly than the Fs-21 peptide by four times. These difference highlights the contribution of the bulky N-terminal MABA group. Furthermore, it is found that helix−turn−helix conformation was the most stable state at 300 K, instead of the expected full helix conformation. At 273 K, however, the full helix became the most stable state. The turn structure was found to be stabilized mainly by the hydrophobic interactions. Statistical analysis of high-resolution PDB structures indicated that most helices are shorter than 16 amino acids. Taken together, we suggest that the intrinsic property of polypeptide chain dictates the formation of short helices in proteins.

Role of solvent in determining conformational preferences of alanine dipeptide in water.

Evidence from a variety of spectroscopic probes indicates that (phi, psi) values corresponding to the left-handed polyproline II helix (P(II)) are preferred for short alanine-based peptides in water. On the basis of results from theoretical studies, it is believed that the observed preference is dictated by favorable peptide-solvent interactions, which are realized through formation of optimal hydrogen-bonding water bridges between peptide donor and acceptor groups. In the present study, we address this issue explicitly by analyzing the hydration structure and thermodynamics of 16 low-energy conformers of the alanine dipeptide (N-acetylalanine-N'-methylamide) in liquid water. Monte Carlo simulations in the canonical ensemble were performed under ambient conditions with all-atom OPLS parameters for the alanine dipeptide and the TIP5P model for water. We find that the number of hydrogen-bonded water molecules connecting the peptide group donor and acceptor atoms has no effect on the solvation thermodynamics. Instead, the latter are determined by the work done to fully hydrate the peptide. This work is minimal for conformations that are characterized by a minimal overlap of the primary hydration shells around the peptide donor and acceptor atoms. As a result, peptide-solvent interactions favor "compact" conformations that do not include P(II)-like geometries. Our main conclusion is that the experimentally observed preference for P(II) does not arise due to favorable direct interactions between the peptide and water molecules. Instead, the latter act to unmask underlying conformational preferences that are a consequence of minimizing intrapeptide steric conflicts.

Polyproline II Helix Conformation in a Proline-Rich Environment: A Theoretical Study

Interest centers here on whether a polyproline II helix can propagate through adjacent non-proline residues, and on shedding light on recent experimental observations suggesting the presence of significant PP II structure in a short alanine-based peptide with no proline in the sequence. For this purpose, we explored the formation of polyproline II helices in proline-rich peptides with the sequences Ac-(Pro) 3-X-(Pro) 3-Gly-Tyr-NH 2, with X = Pro (PPP), Ala (PAP), Gln (PQP), Gly (PGP), and Val (PVP), and Ac-(Pro) 3-Ala-Ala-(Pro) 3-Gly-Tyr-NH 2 (PAAP), by using a theoretical approach that includes a solvent effect as well as cis ↔ trans isomerization of the peptide groups and puckering conformations of the pyrrolidine ring of the proline residues. Since 13C chemical shifts have proven to be useful for identifying secondary-structure preferences in proteins and peptides, and because values of the dihedral angles ( ϕ, ψ) are the main determinants of their magnitudes, we have, therefore, computed the Boltzmann-averaged 13C chemical shifts for the guest residues in the PXP peptide ( X = Pro, Ala, Gln, Gly, and Val) with a combination of approaches, involving molecular mechanics, statistical mechanics, and quantum mechanics. In addition, an improved procedure was used to carry out the conformational searches and to compute the solvent polarization effects faster and more accurately than in previous work. The current theoretical work and additional experimental evidence show that, in short proline-rich peptides, alanine decreases the polyproline II helix content. In particular, the theoretical evidence accumulated in this work calls into question the proposal that alanine has a strong preference to adopt conformations in the polyproline II region of the Ramachandran map.

All-atom generalized-ensemble simulations of small proteins

We give an overview of some generalized-ensemble techniques that have proven successful in all-atom simulations of proteins. We show that these techniques enable efficient investigations of secondary structure formation and folding in peptides and small proteins. Results are presented for various alanine-based artificial peptides and a small protein, the 36-residued villin headpiece subdomain (HP-36). Our results indicate that all-atom simulations of proteins may be more restricted by the accuracy of the present energy functions than by the efficiency of the search algorithms.

Fast folding dynamics of an α-helical peptide with bulky side chains

The fast folding/unfolding dynamics of the α-helical polypeptide poly-N5-(3-hydroxypropyl)-L-glutamine, PHPG, was observed in D2O after a nanosecond laser-induced temperature jump of 1 °C. The sudden rise in temperature disturbs the helix–coil equilibrium of the peptide and the ensuing structural relaxation was monitored by time-resolved IR-spectroscopy of the amide I′ band with a time resolution of 12 ns. After an initial relaxation on the time scale of a few nanoseconds, probably arising from a local rearrangement of random coil or side chain structures, the helix–coil relaxation proceeds on the time scale of a few hundred nanoseconds. The observed time constant of 227 ns is very similar to the helix–coil relaxation time constants observed for short alanine-based peptides, in spite of the greater length and more complex and bulky side chain of PHPG.

Role of Backbone Hydration and Salt-Bridge Formation in Stability of α-Helix in Solution

We test molecular level hypotheses for the high thermal stability of α-helical conformations of alanine-based peptides by performing detailed atomistic simulations of a 20-amino-acid peptide with explicit treatment of water. To assess the contribution of large side chains to α-helix stability through backbone desolvation and salt-bridge formation, we simulate the alanine-rich peptide, Ac-YAEAAKAAEAAKAAEAAKAF-Nme, referred to as the EK peptide, that has three pairs of “ i, i + 3” glutamic acid(−) and lysine(+) substitutions. Efficient configurational sampling of the EK peptide over a wide temperature range enabled by the replica exchange molecular dynamics technique allows characterization of the stability of α-helix with respect to heat-induced unfolding. We find that near ambient temperatures, the EK peptide predominately samples α-helical configurations with 80% fractional helicity at 300 K. The helix melts over a broad range of temperatures with melting temperature, T m, equal to 350 K, that is significantly higher than the T m of a 21-residue polyalanine peptide, A 21. Salt-bridges between oppositely charged Glu − and Lys + side chains can, in principle, provide thermal stability to α-helical conformers. For the specific EK peptide sequence, we observe infrequent formation of Glu-Lys salt-bridges (with ∼10–20% probability) and therefore we conclude that salt-bridge formation does not contribute significantly to the EK peptide's helical stability. However, lysine side chains are found to shield specifi c “ i, i + 4” backbone hydrogen bonds from water, indicating that large side-chain substituents can play an important role in stabilizing α-helical configurations of short peptides in aqueous solution through mediation of water access to backbone hydrogen bonds. These observations have implications on molecular engineering of peptides and biomolecules in the design of their thermostable variants where the shielding mechanism can act in concert with other factors such as salt-bridge formation, thereby increasing thermal stability considerably.