Helix forming tendency of valine substituted poly-alanine: a molecular dynamics investigation

In this study, classical molecular dynamics simulations have been carried out on the valine (guest) substituted poly alanine (host) using the host-guest peptide approach to understand the role of valine in the formation and stabilization of helix. Valine has been substituted in the host peptide starting from N terminal to C terminal. Various structural parameters have been obtained from the molecular dynamics simulation to understand the tolerance of helical motif to valine. Depending on the position of valine in the host peptide, it stabilizes (or destabilizes) the formation of the helical structure. The substitution of valine in the poly alanine at some positions has no effect on the helix formation (deformation). It is interesting to observe the coexistence of 3 10 and alpha-helix in the peptides due to the dynamical nature of the hydrogen bonding interaction and sterical interactions.     

Helices and Sheets in vacuo

The structures and properties of unsolvated peptides large enough to possess secondary structure have been examined by experiments and simulations. Some of the factors that stabilize unsolvated helices and sheets have been identified. The charge, in particular, plays a critical role in stabilizing alpha-helices and destabilizing beta-sheets. Some helices are much more stable in vacuum than in aqueous solution. Factors like helix propensity, context, and the incorporation of specific stabilizing interactions have been examined. The helix propensities in vacuum differ from those found in solution. Studies of the hydration of unsolvated peptides can be performed one water molecule at a time. The first few water molecules only bind weakly to unsolvated peptides, and they bind much more strongly to some conformations than to others. The most favorable binding locations are not the protonation sites, but clefts or pockets where a water molecule can establish a network of hydrogen bonds. Non-covalent interactions between secondary structure elements leads to the formation of tertiary structure. Helical peptides assemble into complexes with a variety of intriguing structures. The intramolecular coupling of helices to make antiparallel coiled-coil geometries has also been investigated with model peptides.     

Effect of thioxopeptide bonds on alpha-helix structure and stability

Thioxoamide (thioamide) bonds are nearly isosteric substitutions for amides but have altered hydrogen-bonding and photophysical properties. They are thus well-suited backbone modifications for physicochemical studies on peptides and proteins. The effect of thioxoamides on protein structure and stability has not been subject to detailed experimental investigations up to date. We used alanine-based model peptides to test the influence of single thioxoamide bonds on alpha-helix structure and stability. The results from circular dichroism measurements show that thioxoamides are strongly helix-destabilizing. The effect of an oxo-to-thioxoamide backbone substitution is of similar magnitude as an alanine-to-glycine substitution resulting in a helix destabilization of about 7 kJ/mol. NMR characterization of a helical peptide with a thioxopeptide bond near the N-terminus indicates that the thioxopeptide moiety is tolerated in helical structures. The thioxoamide group is engaged in an i, i+4 hydrogen bond, arguing against the formation of a 3(10)-helical structure as suggested for the N-termini of alpha-helices in general and for thioxopeptides in particular.     

Helix formation in alpha,gamma- and beta,gamma-hybrid peptides: theoretical insights into mimicry of alpha- and beta-peptides

Alpha,gamma- and beta,gamma-hybrid peptides, which are composed of two different homologous amino acid constituents in alternate order, are suggested as novel classes of peptide foldamers. On the basis of a systematic conformational search employing the methods of ab initio MO theory, the possibilities for the formation of periodic secondary structures in these systems are described. The conformational analysis provides a great number of helix conformers widely differing in energy, which can be arranged into three groups: (i) helices with all hydrogen bonds formed in forward direction along the sequence, (ii) helices with all hydrogen bonds in backward direction, and (iii) helices with alternate hydrogen-bond directions (mixed or beta-helices). Most stable are representatives of beta-helices, but their stability decreases considerably in more polar environments in comparison to helix conformers from the other two classes. There is a great similarity between the overall topology of the most stable hybrid peptide helices and typical helices of peptides which are exclusively composed of a single type of homologous amino acids. Thus, the helices of the beta,gamma-hybrid peptides mimic perfectly those of the native alpha-peptides as, for instance, the well-known alpha-helix, whereas the most stable helix conformers of alpha,gamma-hybrid peptides correspond well to the overall structure of beta-peptide helices. The two suggested novel hybrid peptide classes expand considerably the pool of peptide foldamers and may be promising tools in peptide design and in material sciences.     

The Effects of Various Alcohols on the Stability of Melittin: A Molecular Dynamics Study

In this study, 200 ps molecular dynamics simulations were conducted to investigate the effects of various alcohols on the structural stability of melittin. The averaged helicity of melittin remained 80% in pure butanol, whereas it was below 60% both in pure water and in pure methanol. The α‐helix propensity of melittin increased with the aliphatic chain length of the alcohol. Charge‐charge interaction between Lys21 and Arg24 and polar‐nonpolar interaction between Trp19 and Arg22 are probably responsible for the higher structural integrity of the C‐terminal α‐helix over the N‐terminal one. The weaker dielectric constant of longer aliphatic chain length of alcohol possibly reduces the hydrogen bonding between amide protons and surrounding solvent molecules and simultaneously promotes the intramolecular hydrogen bonding in melittin and therefore stabilizes the secondary structure of melittin. The effect of various alcohols on stabilizing melittin is most likely due to their ability to form clusters on the surface of melittin effectively, favoring the formation of intramolecular hydrogen bonds instead of intermolecular hydrogen bonds and promoting the formation of stable α‐helices.     

Modeling of folding and unfolding mechanisms in alanine-based alpha-helical polypeptides

alpha-Helix formation is known to be opposed by the entropy loss due to the folding and favored by the energy of molecular interactions. However, the underlying mechanism of these factors is still being discussed. Here we have used the experimental and calculation data for short alanine-based peptides embedded in water to model the mechanism of helix folding and unfolding and to calculate microscopically the free energy factors of alanine in the frame of helix coil conformational integrals. Classical helix-coil transition theories take into account the interactions in a peptide chain only if the i, i + 3 peptide bond participates in hydrogen bonding. But quantum mechanical calculations showed that interactions of the i, i + 2 peptide bond play an important role in helix folding too. We also included the short-range repulsive interactions due to molecular steric clashes and the end effects due to polar/hydrogen-bonding interactions at the N and C termini. The helix and coil regions of peptide conformational space were defined using an experimental steric criterion for hydrogen bonding. Arginine helix propensity was discussed and estimated. Monte Carlo numerical simulations of thermodynamics and kinetics for the 21 amino acid alpha-helical polypeptide Ac-A5(AAARA)3A-NMe were carried out and found to be in an agreement with the experimental results.     

Specific and nonspecific effects of glycosylation

Glycosylation regulates vital cellular processes and dramatically influences protein folding and stability. In particular, experiments have demonstrated that asparagine (N)-linked disaccharides drive a "conformational switch" in a model peptide. The present work investigates this conformational switch via extensive atomically detailed replica exchange molecular dynamics simulations in explicit solvent. To distinguish the effects of specific and nonspecific interactions upon the peptide conformational ensemble, these simulations considered model peptides that were N-linked to a disaccharide and to a steric crowder of the same shape. The simulations are remarkably consistent with experiment and provide detailed insight into the peptide structure ensemble. They suggest that steric crowding by N-linked disaccharides excludes extended conformations, but does not significantly impact the tetrahedral structure of the surrounding solvent or otherwise alter the peptide free energy surface. However, the combination of steric crowding with specific hydrogen bonds and hydrophobic stacking interactions more dramatically impacts the peptide ensemble and stabilizes new structures.     

Identifying residual structure in intrinsically disordered systems: a 2D IR spectroscopic study of the GVGXPGVG peptide

The peptide amide-I vibration of a proline turn encodes information on the turn structure. In this study, FTIR, two-dimensional IR spectroscopy and molecular dynamics simulations were employed to characterize the varying turn conformations that exist in the GVGX(L)PGVG family of disordered peptides. This analysis revealed that changing the size of the side chain at the X amino acid site from Gly to Ala to Val substantially alters the conformation of the peptide. To quantify this effect, proline peak shifts and intensity changes were compared to a structure-based spectroscopic model. These simulated spectra were used to assign the population of type-II β turns, bulged turns, and irregular β turns for each peptide. Of particular interest was the Val variant commonly found in the protein elastin, which contained a 25% population of irregular β turns containing two peptide hydrogen bonds to the proline C═O.     

Solvent effect on the folding dynamics and structure of E6-associated protein characterized from ab initio protein folding simulations

Solvent effect on protein conformation and folding mechanism of E6-associated protein (E6ap) peptide are investigated using a recently developed charge update scheme termed as adaptive hydrogen bond-specific charge (AHBC). On the basis of the close agreement between the calculated helix contents from AHBC simulations and experimental results, we observed based on the presented simulations that the two ends of the peptide may simultaneously take part in the formation of the helical structure at the early stage of folding and finally merge to form a helix with lowest backbone RMSD of about 0.9 A? in 40% 2,2,2-trifluoroethanol solution. However, in pure water, the folding may start at the center of the peptide sequence instead of at the two opposite ends. The analysis of the free energy landscape indicates that the solvent may determine the folding clusters of E6ap, which subsequently leads to the different final folded structure. The current study demonstrates new insight to the role of solvent in the determination of protein structure and folding dynamics.     

Secondary structure bias in generalized Born solvent models: comparison of conformational ensembles and free energy of solvent polarization from explicit and implicit solvation

The effects of the use of three generalized Born (GB) implicit solvent models on the thermodynamics of a simple polyalanine peptide are studied via comparing several hundred nanoseconds of well-converged replica exchange molecular dynamics (REMD) simulations using explicit TIP3P solvent to REMD simulations with the GB solvent models. It is found that when compared to REMD simulations using TIP3P the GB REMD simulations contain significant differences in secondary structure populations, most notably an overabundance of alpha-helical secondary structure. This discrepancy is explored via comparison of the differences in the electrostatic component of the free energy of solvation (DeltaDeltaG(pol)) between TIP3P (via thermodynamic Integration calculations), the GB models, and an implicit solvent model based on the Poisson equation (PE). The electrostatic components of the solvation free energies are calculated using each solvent model for four representative conformations of Ala10. Since the PE model is found to have the best performance with respect to reproducing TIP3P DeltaDeltaG(pol) values, effective Born radii from the GB models are compared to effective Born radii calculated with PE (so-called perfect radii), and significant and numerous deviations in GB radii from perfect radii are found in all GB models. The effect of these deviations on the solvation free energy is discussed, and it is shown that even when perfect radii are used the agreement of GB with TIP3P DeltaDeltaG(pol) values does not improve. This suggests a limit to the optimization of the effective Born radius calculation and that future efforts to improve the accuracy of GB models must extend beyond such optimizations.     

Energetic switch of the proline effect in collision‐induced dissociation of singly and doubly protonated peptide Ala‐Ala‐Arg‐Pro‐Ala‐Ala

Suppression of the selective cleavage at N‐terminal of proline is observed in the peptide cleavage by proteolytic enzyme trypsin and in the fragment ion mass spectra of peptides containing Arg‐Pro sequence. An insight into the fragmentation mechanism of the influence of arginine residue on the proline effect can help in prediction of mass spectra and in protein structure analysis. In this work, collision‐induced dissociation spectra of singly and doubly charged peptide AARPAA were studied by ESI MS/MS and theoretical calculation methods. The proline effect was evaluated by comparing the experimental ratio of fragments originated from cleavage of different amide bonds. The results revealed that the backbone amide bond cleavage was selected by the energy barrier height of the fragmentation pathway although the strong proton affinity of the Arg side chain affected the stereostructure of the peptide and the dissociation mechanism. The thermodynamic stability of the fragment ions played a secondary role in the abundance ratio of fragments generated via different pathways. Fragmentation studies of protonated peptide AACitPAA supported the energy‐dependent hypothesis. The results provide an explanation to the long‐term arguments between the steric conflict and the proton mobility mechanisms of proline effect.     

Elongation of ordered peptide aggregate of an amyloidogenic hexapeptide NFGAIL observed in molecular dynamics simulations with explicit solvent

The mechanisms by which amyloidogenic peptides and proteins form soluble toxic oligomers remain elusive. We have studied the formation of partially ordered tetramers and well-ordered octamers of an amyloidogenic hexapeptide NFGAIL (residues 22-27 of the human islet amyloid polypeptide) in our previous work. Continuing the effort, we here probe the beta-sheet elongation process by a combined total of 2.0 micros molecular dynamics simulations with explicit solvent. In a set of 10 simulations with the peptides restrained to the extended conformation, we observed that the main growth mode was elongation along the beta-sheet hydrogen bonds through primarily a two-stage process. Driven by hydrophobic forces, the peptides initially attached to the surface of the ordered oligomer, moved quickly to the beta-sheet edges, and formed stable beta-sheet hydrogen bonds. Addition of peptides to the existing oligomer notably improved the order of the peptide aggregate in which labile outer layer beta-sheets were stabilized, which provides good templates for further elongation. These simulations suggested that elongation along the beta-sheet hydrogen bonds occurs at the intermediate stage when low-weight oligomers start to form. We did not observe significant preference toward either parallel or antiparallel beta-sheets at the elongation stage for this peptide. In another set of 10 unrestrained simulations, the dominant growth mode was disordered aggregation. Taken together, these results offered a glimpse at the molecular events leading to the formation of ordered and disordered low-weight oligomers.     

Molecular dynamics and thermodynamics of protein-RNA interactions: mutation of a conserved aromatic residue modifies stacking interactions and structural adaptation in the U1A-stem loop 2 RNA complex.

Molecular dynamics (MD) simulations and free energy component analysis have been performed to evaluate the molecular origins of the 5.5 kcal/mol destabilization of the complex formed between the N-terminal RNP domain of U1A and stem loop 2 of U1 snRNA upon mutation of a conserved aromatic residue, Phe56, to Ala. MD simulations, including counterions and water, have been carried out on the wild type and Phe56Ala peptide-stem loop 2 RNA complexes, the free wild type and Phe56Ala peptides, and the free stem loop 2 RNA. The MD structure of the Phe56Ala-stem loop 2 complex is similar to that of the wild type complex except the stacking interaction between Phe56 and A6 of stem loop 2 is absent and loop 3 of the peptide is more dynamic. However, the MD simulations predict large changes in the structure and dynamics of helix C and increased dynamic range of loop 3 for the free Phe56Ala peptide compared to the wild type peptide. Since helix C and loop 3 are highly variable regions of RNP domains, this indicates that a significant contribution to the reduced affinity of the Phe56Ala peptide for RNA results from cooperation between highly conserved and highly variable regions of the RNP domain of U1A. Surprisingly, these structural effects, which are manifested as cooperative free energy changes, occur in the free peptide, rather than in the complex, and are revealed only by study of both the initial and final states of the complexation process. Free energy component analysis correctly accounts for the destabilization of the Phe56Ala-stem loop 2 complex, and indicates that approximately 80% of the destabilization is due to the loss of the stacking interaction and approximately 20% is due to differences in U1A adaptation.     

A 10-A spectroscopic ruler applied to short polyprolines

Fluorescence resonance energy transfer (FRET) from the amino acid tryptophan (Trp) as donor and a 2,3-diazabicyclo[2.2.2]oct-2-ene-labeled asparagine (Dbo) as acceptor in peptides of the general structure Trp-(Pro)n-Dbo-NH2 (n = 1-6) was investigated by steady-state and time-resolved fluorescence, CD, and NMR spectroscopy as well as by molecular dynamics (MD) simulations (GROMOS96 force field). The Trp/Dbo FRET pair is characterized by a very short F?rster radius (R0 ca. 9 A), which allowed distance determinations in such short peptides. Water and propylene glycol were investigated as solvents. The peptides were designed to show an early nucleation of the poly(Pro)II (PPII) secondary helix structure for n > or = 2, which was confirmed by their CD spectra. The shortest peptide (n = 1) adopts preferentially the trans conformation about the Trp-Pro bond, as confirmed by NMR spectra. The FRET efficiencies ranged 2-72% and were found to depend sensitively on the peptide length, i.e., the number of intervening proline residues. The analysis of the FRET data at different levels of theory (assuming either a fixed distance or distance distributions according to a wormlike chain or Gaussian model) afforded donor-acceptor distances between ca. 8 A (n = 1) and ca. 16 A (n = 6) in water, which were found to be similar or slightly higher in propylene glycol. The distances afforded by the Trp/Dbo FRET pair were found to be reasonable in comparison to literature data, expectations from the PPII helix structure, and the results from MD simulations. The persistence lengths for the longer peptides were found to lie at 30-70 A in water and 220 +/- 40 A in propylene glycol, suggesting a more rigid PPII helical structure in propylene glycol. A detailed comparison with literature data on FRET in polyprolines demonstrates that the donor-acceptor distances extracted by FRET are correlated with the F?rster radii of the employed FRET pairs. This demonstrates the limitations of using FRET as a spectroscopic ruler for short polyprolines, which is presumably due to the breakdown of the point dipole approximation in F?rster theory, when the size of the chromophores becomes comparable or larger than the distances under investigation.     

Molecular Dynamics Simulation of Ester‐Linked Hen Egg White Lysozyme Reveals the Effect of Missing Backbone Hydrogen Bond Donors on the Protein Structure

The three‐dimensional structure of a protein is stabilized by a number of different atomic interactions. One of these is hydrogen bonding. Its influence on the spatial structure of the hen egg white lysozyme is investigated by replacing peptide bonds (except those of the two proline residues) by ester bonds. Molecular dynamics simulations of native and ester‐linked lysozyme are compared with the native crystal structure and with NOE distance bounds derived from solution NMR experiments. The ester‐linked protein shows a slight compaction while losing its native structure. However, it does not unfold completely. The structure remains compact due to its hydrophobic core and a changed network of hydrogen bonds involving side chains.     

Molecular dynamics study of the solvation of an alpha-helical transmembrane peptide by DMSO

A 10-ns molecular dynamics study of the solvation of a hydrophobic transmembrane helical peptide in dimethyl sulfoxide (DMSO) is presented. The objective is to analyze how this aprotic polar solvent is able to solvate three groups of amino acid residues (i.e., polar, apolar, and charged) that are located in a stable helical region of a transmembrane peptide. The 25-residue peptide (sMTM7) used mimics the cytoplasmic proton hemichannel domain of the seventh transmembrane segment (TM7) from subunit a of H(+)-V-ATPase from Saccharomyces cerevisiae. The three-dimensional structure of peptide sMTM7 in DMSO has been previously solved by NMR spectroscopy. The radial and spatial distributions of the DMSO molecules surrounding the peptide as well as the number of hydrogen bonds between DMSO and the side chains of the amino acid residues involved are extracted from the molecular dynamics simulations. Analysis of the molecular dynamics trajectories shows that the amino acid side chains are fully embedded in DMSO. Polar and positively charged amino acid side chains have dipole-dipole interactions with the oxygen atom of DMSO and form hydrogen bonds. Apolar residues become solvated by DMSO through the formation of a hydrophobic pocket in which the methyl groups of DMSO are pointing toward the hydrophobic side chains of the residues involved. The dual solvation properties of DMSO cause it to be a good membrane-mimicking solvent for transmembrane peptides that do not unfold due to the presence of DMSO.     

Helix Nucleation by the Smallest Known α‐Helix in Water

Cyclic pentapeptides (e.g. Ac‐(cyclo‐1,5)‐[KAXAD]‐NH₂; X=Ala, 1 ; Arg, 2 ) in water adopt one α‐helical turn defined by three hydrogen bonds. NMR structure analysis reveals a slight distortion from α‐helicity at the C‐terminal aspartate caused by torsional restraints imposed by the K(i)–D(i+4) lactam bridge. To investigate this effect on helix nucleation, the more water‐soluble 2 was appended to N‐, C‐, or both termini of a palindromic peptide ARAARAARA (≤5 % helicity), resulting in 67, 92, or 100 % relative α‐helicity, as calculated from CD spectra. From the C‐terminus of peptides, 2 can nucleate at least six α‐helical turns. From the N‐terminus, imperfect alignment of the Asp5 backbone amide in 2 reduces helix nucleation, but is corrected by a second unit of 2 separated by 0–9 residues from the first. These cyclic peptides are extremely versatile helix nucleators that can be placed anywhere in 5–25 residue peptides, which correspond to most helix lengths in protein–protein interactions.     

Mechanism of helix nucleation and propagation: microscopic view from microsecond time scale MD simulations

Microsecond time scale molecular dynamics simulations of the 13-residue peptide RN24 were carried out to investigate the mechanism of helix nucleation and propagation. An extended and an ideal alpha-helical conformation were used as starting structures. NOE-derived interatomic distances were compared with distances calculated from the simulations, showing good agreement between experimental and simulation results. Based on almost 200 helix nucleation events observed, beta-turn and 3(10)-helix play an important role in helix nucleation; in most cases, helix nucleation is preceded by the formation of a short-lived beta-turn (60% probability) or 3(10)-helix (20% probability), and the conversion from beta-turn to alpha-turn involves bifurcated hydrogen bonds. Helix propagation in RN24 appears to occur preferentially from the N-terminus to the C-terminus, and helix unfolding preferentially in the opposite direction.     

Conformational equilibrium in alanine-rich peptides probed by reversible stretching simulations

Reversible stretching of the alanine-rich peptide 3K (Proc. Natl. Acad. Sci. USA 1989, 86, 5286-5290) and its analogue MW (Nature 1992, 359, 653-655) is examined using molecular dynamics simulations in explicit water. In both cases, sampling of the extension pathway is obtained on the 10 ns time scale by applying an adaptive biasing force. The free energy profile reveals a single minimum associated with a contiguous alpha-helix. Short 3(10)-helical motifs are observed in folded as well as extended conformations, in accordance with their proposed role as folding intermediates. The native 3(10)-helical content of both peptides is found, however, to be no higher than a few percent. Difficulties in both the definition and the detection of secondary structure motifs, most notably in relation to bifurcated hydrogen bonds, are proposed to account for the discrepancy between 3(10)-helical propensities reported by several authors, based on experimental and computational results.