Molecular dynamics simulations of peptides from the central domain of smooth muscle caldesmon

The central domain of smooth muscle caldesmon contains a highly charged region consisting of ten 13-residue repeats. Experimental evidence obtained from the intact protein and fragments thereof suggests that this entire region forms a single stretch of stable alpha-helix. We have carried out molecular dynamics simulations on peptides consisting of one, two and three repeats to examine the mechanism of alpha-helical stability of the central domain at the atomic level. All three peptides show high helical stability on the timescale of the MD simulations. Deviations from alpha-helical structure in all the simulations arise mainly from the formation of long stretches of pi-helix. Interconversion between alpha-helical and pi-helical conformations occurs through insertion of water molecules into alpha-helical hydrogen bonds and subsequent formation of reverse turns. The alpha-helical structure is stabilized by electrostatic interactions (salt bridges) between oppositely charged sidechains with i,i+4 spacings, while the pi-helix is stabilized by i,i+5 salt bridge interactions. Possible i,i+3 salt bridges are of minor importance. There is a strong preference for salt bridges with a Glu residue N-terminal to a basic sidechain as compared to the opposite orientation. In the double and triple repeat peptides, strong i,i+4 salt bridges exist between the last Glu residue of one repeat and the first Lys residue of the next. This demonstrates a relationship between the repetitive nature of the central domain sequence and its ability to form very long stretches of alpha-helical structure.     

Application of principal component analysis to determine the key structural features contributing to iron superoxide dismutase thermostability

Iron superoxide dismutase (Fe-SOD) is predominantly found in bacteria and mitochondria. The thermal stability of Fe-SOD from different sources can vary dramatically. We have studied the influence of structural parameters on Fe-SOD thermostability by principal component analysis (PCA). The results show that an increased α-helical and turn content, an increased α-helix and loop length, an increase in the number of main-main chains and charged-uncharged hydrogen bonds, a decrease in the 3(10) -helix content, and a decreased β-strand and loop length are all important factors for Fe-SOD thermostability. Interestingly, the use of charged residues to form salt bridges is tendentious in thermophilic Fe-SOD. Negatively charged Arg and positively charged Glu are efficiently used to form salt bridges. The cooperative action of the exposed area, the hydrogen bonds, and the secondary structure plays a crucial role in resisting high temperatures, which demonstrates that the increased stability of thermophilic Fe-SOD is provided by several structural factors acting together.     

NMR determination of pKa values in α-synuclein

The intrinsically unfolded protein α-synuclein has an N-terminal domain with seven imperfect KTKEGV sequence repeats and a C-terminal domain with a large proportion of acidic residues. We characterized pK(a) values for all 26 sites in the protein that ionize below pH 7 using 2D (1) H-(15) N HSQC and 3D C(CO)NH NMR experiments. The N-terminal domain shows systematically lowered pK(a) values, suggesting weak electrostatic interactions between acidic and basic residues in the KTKEGV repeats. By contrast, the C-terminal domain shows elevated pK(a) values due to electrostatic repulsion between like charges. The effects are smaller but persist at physiological salt concentrations. For α-synuclein in the membrane-like environment of sodium dodecylsulfate (SDS) micelles, we characterized the pK(a) of His50, a residue of particular interest since it is flanked within one turn of the α-helix structure by the Parkinson's disease-linked mutants E46K and A53T. The pK(a) of His50 is raised by 1.4 pH units in the micelle-bound state. Titrations of His50 in the micelle-bound states of the E46K and A53T mutants show that the pK(a) shift is primarily due to interactions between the histidine and the sulfate groups of SDS, with electrostatic interactions between His50 and Glu46 playing a much smaller role. Our results indicate that the pK(a) values of uncomplexed α-synuclein differ significantly from random coil model peptides even though the protein is intrinsically unfolded. Due to the long-range nature of electrostatic interactions, charged residues in the α-synuclein sequence may help nucleate the folding of the protein into an α-helical structure and confer protection from misfolding.     

Relationship of sidechain hydrophobicity and α-helical propensity on the stability of the single-stranded amphipathic α-helix

The aim of the present investigation is to determine the effect of α-helical propensity and sidechain hydrophobicity on the stability of amphipathic α-helices. Accordingly, a series of 18-residue amphipathic α-helical peptides has been synthesized as a model system where all 20 amino acid residues were substituted on the hydrophobic face of the amphipathic α-helix. In these experiments, all three parameters (sidechain hydrophobicity, α-helical propensity and helix stability) were measured on the same set of peptide analogues. For these peptide analogues that differ by only one amino acid residue, there was a 0.96 kcal/mole difference in α-helical propensity between the most (Ala) and the least (Gly) α-helical analogue, a 12.1-minute difference between the most (Phe) and the least (Asp) retentive analogue on the reversed-phase column, and a 32.3°C difference in melting temperatures between the most (Leu) and the least (Asp) stable analogue. The results show that the hydrophobicity and α-helical propensity of an amino acid sidechain are not correlated with each other, but each contributes to the stability of the amphipathic α-helix. More importantly, the combined effects of α-helical propensity and sidechain hydrophobicity at a ratio of about 2:1 had optimal correlation with α-helix stability. These results suggest that both α-helical propensity and sidechain hydrophobicity should be taken into consideration in the design of α-helical proteins with the desired stability.     

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.     

Modulation of Constitutive Activity and Signaling Bias of the Ghrelin Receptor by Conformational Constraint in the Second Extracellular Loop

Based on a rare, natural Glu for Ala-204(C+6) variant located six residues after the conserved Cys residue in extracellular loop 2b (ECL2b) associated with selective elimination of the high constitutive signaling of the ghrelin receptor, this loop was subjected to a detailed structure functional analysis. Introduction of Glu in different positions demonstrated that although the constitutive signaling was partly reduced when introduced in position 205(C+7) it was only totally eliminated in position 204(C+6). No charge-charge interaction partner could be identified for the Glu(C+6) variant despite mutational analysis of a number of potential partners in the extracellular loops and outer parts of the transmembrane segments. Systematic probing of position 204(C+6) with amino acid residues of different physicochemical properties indicated that a positively charged Lys surprisingly provided phenotypes similar to those of the negatively charged Glu residue. Computational chemistry analysis indicated that the propensity for the C-terminal segment of extracellular loop 2b to form an extended α-helix was increased from 15% in the wild type to 89 and 82% by introduction in position 204(C+6) of a Glu or a Lys residue, respectively. Moreover, the constitutive activity of the receptor was inhibited by Zn²⁺ binding in an engineered metal ion site, stabilizing an α-helical conformation of this loop segment. It is concluded that the high constitutive activity of the ghrelin receptor is dependent upon flexibility in the C-terminal segment of extracellular loop 2 and that mutations or ligand binding that constrains this segment and thereby conceivably the movements of transmembrane domain V relative to transmembrane domain III inhibits the high constitutive signaling.     

Positional dependence of the effects of negatively charged Glu side chains on the stability of two-stranded α-helical coiled-coils

The effects on protein stability of negatively charged Glu side chains at different positions along the length of the α-helix were investigated in the two-stranded α-helical coiled-coil. A native coiled-coil has been designed which consists of two identical 35 residue polypeptide chains with a heptad repeat Q_gV_aG_bA_cL_dQ_eK_f and a Cys residue at position 2 to allow the formation of an interchain 2-2′ disulphide bridge. This coiled-coil contains no intra- or interchain electrostatic interactions and served as a control for peptides in which Glu was substituted for Gln in the e or g heptad positions. The effect of the substitutions on stability was determined by urea denaturation at 20°C with the degree of unfolding monitored by circular dichroism spectroscopy. A Glu substituted for Gln near the N-terminus in each chain of the coiled-coil stabilizes the coiled-coil at pH 7, consistent with the charge–helix dipole interaction model. This stability increase is modulated by pH change and the addition of salt (KCl or guanidine hydrochloride), confirming the electrostatic nature of the effect. In contrast, Glu substitution in the middle of the helix destabilizes the coiled-coil because of the lower helical propensity and hydrophobicity of Glu compared with Gln at pH 7. Taking the intrinsic differences into account, the apparent charge–helix dipole interaction at the N-terminus is approximately 0.35 kcal/mol per Glu substitution. A Glu substitution at the C-terminus destabilizes the coiled-coil more than in the middle owing to the combined effects of intrinsic destabilization and unfavourable charge–helix dipole interaction with the negative pole of the helix dipole. The estimated destabilizing charge–helix dipole interaction of 0.08 kcal/mol is smaller than the stabilizing interaction at the N-terminus. The presence of a 2-2′disulphide bridge appears to have little influence on the magnitude of the charge–helix dipole interactions at either end of the coiled-coil. © 1997 European Peptide Society and John Wiley & Sons, Ltd.     

Peptide Binding to a PDZ Domain by Electrostatic Steering via Nonnative Salt Bridges

We have captured the binding of a peptide to a PDZ domain by unbiased molecular dynamics simulations. Analysis of the trajectories reveals on-pathway encounter complex formation, which is driven by electrostatic interactions between negatively charged carboxylate groups in the peptide and positively charged side chains surrounding the binding site. In contrast, the final stereospecific complex, which matches the crystal structure, features completely different interactions, namely the burial of the hydrophobic side chain of the peptide C-terminal residue and backbone hydrogen bonds. The simulations show that nonnative salt bridges stabilize kinetically the encounter complex during binding. Unbinding follows the inverse sequence of events with the same nonnative salt bridges in the encounter complex. Thus, in contrast to protein folding, which is driven by native interactions, the binding of charged peptides can be steered by nonnative interactions, which might be a general mechanism, e.g., in the recognition of histone tails by bromodomains.     

Proline-induced constraints in alpha-helices

The disrupting effect of a prolyl residue on an α-helix has been analyzed by means of conformational energy computations. In the preferred, nearly α-helical conformations of Ac-Ala₄-Pro-NHMe and of Ac-Ala₇-Pro-Ala₇-NHMe, only the residue preceding Pro is not α-helical, while all other residues can occur in the α-helical A conformation; i.e., it is sufficient to introduce a conformational change of only one residue in order to accommodate proline in a distorted α-helix. Other low-energy conformations exist in which the conformational state of three residues preceding proline is altered considerably; on the other hand, another conformation in which these three residues retain the near-α-helical A-conformational state (with up to 26° changes of their dihedral angles ? and ψ, and a 48° change in one ω from those of the ideal α-helix) has a considerably higher energy. These conclusions are not altered by the substitution of other residues in the place of the Ala preceding Pro. The conformations of the peptide chain next to prolyl residues in or near an α-helix have been analyzed in 58 proteins of known structure, based on published atomic coordinates. Of 331 α-helices, 61 have a Pro at or next to their N-terminus, 21 have a Pro next to their C-terminus, and 30 contain a Pro inside the helix. Of the latter, 16 correspond to a break in the helix, 9 are located inside distorted first turns of the helix, and 5 are parts of irregular helices. Thus, the reported occurrence of prolyl residues next to or inside observed α-helices in proteins is consistent with the computed steric and energetic requirements of prolyl peptides.     

Probing the structure of membrane proteins with electron spin echo envelope modulation spectroscopy

A new approach has been developed to probe the structural properties of membrane peptides and proteins using the pulsed electron paramagnetic resonance technique of electron spin echo envelope modulation (ESEEM) spectroscopy and the α-helical M2δ subunit of the acetylcholine receptor incorporated into phospholipid bicelles. To demonstrate the practicality of this method, a cysteine-mutated nitroxide spin label (SL) is positioned 1, 2, 3, and 4 residues away from a fully deuterated Val side chain (denoted i + 1 to i + 4). The characteristic periodicity of the α-helical structure gives rise to a unique pattern in the ESEEM spectra. In the i + 1 and i + 2 samples, the 2H nuclei are too far away to be detected. However, with the 3.6 residue per turn pattern of an α-helix, the i + 3 and i + 4 samples reveal a strong signal from the 2H nuclei of the Val side chain. Modeling studies verify these data suggesting that the closest 2H-labeled Val to SL distance would in fact be expected in the i + 3 and i + 4 samples. This technique is very advantageous, because it provides pertinent qualitative structural information on an inherently difficult system like membrane proteins in a short period of time (minutes) with small amounts of protein (μg).     

Solution structure of P01, a natural scorpion peptide structurally analogous to scorpion toxins specific for apamin-sensitive potassium channel

The venom of the North African scorpion Androctonus mauretanicus mauretanicus possesses numerous highly active neurotoxins that specifically bind to various ion channels. One of these, P05, has been found to bind specifically to calcium-activated potassium channels and also to compete with apamin, a toxin extracted from bee venom. Besides the highly potent ones, several of these peptides (including that of P01) have been purified and been found to possess only a very weak, although significant, activity in competition with apamin. The amino acid sequence of P01 shows that it is shorter than P05 by two residues. This deletion occurs within an α-helix stretch (residues 5–12). This α-helix has been shown to be involved in the interaction of P05 with its receptor via two arginine residues. These two arginines are absent in the P01 sequence. Furthermore, a proline residue in position 7 of the P01 sequence may act as an α-helix breaker. We have determined the solution structure of P01 by conventional two-dimensional ¹H nuclear magnetic resonance and show that 1) the proline residue does not disturb the α-helix running from residues 5 to 12; 2) the two arginines are topologically replaced by two acidic residues, which explains the drop in activity; 3) the residual binding activity may be due to the histidine residue in position 9; and 4) the overall secondary structure is conserved, i.e., an α-helix running from residues 5 to 12, two antiparallel stretches of β-sheet (residues 15–20 and 23–27) connected by a type I′ β-turn, and three disulfide bridges connecting the α-helix to the β-sheet.     

Evolutionary origin of a secondary structure: π-helices as cryptic but widespread insertional variations of α-helices that enhance protein functionality

Formally annotated π-helices are rare in protein structures but have been correlated with functional sites. Here, we analyze protein structures to show that π-helices are the same as structures known as α-bulges, α-aneurisms, π-bulges, and looping outs, and are evolutionarily derived by the insertion of a single residue into an α-helix. This newly discovered evolutionary origin explains both why π-helices are cryptic, being rarely annotated despite occurring in 15% of known proteins, and why they tend to be associated with function. An analysis of π-helices in the diverse ferritin-like superfamily illustrates their tendency to be conserved in protein families and identifies a putative π-helix-containing primordial precursor, a “missing link” intermediary form of the ribonucleotide reductase family, vestigial π-helices, and a novel function for π-helices that we term a “peristaltic-like shift.” This new understanding of π-helices paves the way for this generally overlooked motif to become a noteworthy feature that will aid in tracing the evolution of many protein families, guide investigations of protein and π-helix functionality, and contribute additional tools to the protein engineering toolkit.     

Protein-protein interactions in actin-myosin binding and structural effects of R405Q mutation: a molecular dynamics study

Detailed residue-wise interactions involved in the binding of myosin to actin in the rigor conformation without nucleotides have been examined using molecular dynamics simulations of the chicken skeletal myosin head complexed with two actin monomers, based on the cryo-microscopic model of Holmes et al. (Nature 2003;425:423-427). The overall interaction is largely electrostatic in nature, because of the charged residues in the four loops surrounding the central primary binding site. The 50k/20k loop, disordered in crystal structures and in simulations of free myosin in solution, was found to be in a conformation stabilized with 1 - 2 internal salt bridges. The cardiomyopathy loop forms 2 - 3 interprotein salt bridges with actin monomers upon binding, whereas its Arg405 residue, the mutation site associated with the hypertrophic cardiomyopathy, forms a strong salt bridge with Glu605 in the neighboring helix away from actin in the actin-bound myosin. The myopathy loop of the R405Q mutant maintains a high degree of two-strand beta-sheet character when bound to actin with the corresponding salt bridges broken.     

Conformational choreography of a molecular switch region in myelin basic protein--molecular dynamics shows induced folding and secondary structure type conversion upon threonyl phosphorylation in both aqueous and membrane-associated environments

The 18.5 kDa isoform of myelin basic protein is essential to maintaining the close apposition of myelin membranes in central nervous system myelin, but its intrinsic disorder (conformational dependence on environment), a variety of post-translational modifications, and a diversity of protein ligands (e.g., actin and tubulin) all indicate it to be multifunctional. We have performed molecular dynamics simulations of a conserved central segment of 18.5 kDa myelin basic protein (residues Glu80-Gly103, murine sequence numbering) in aqueous and membrane-associated environments to ascertain the stability of constituent secondary structure elements (α-helix from Glu80-Val91 and extended poly-proline type II from Thr92-Gly103) and the effects of phosphorylation of residues Thr92 and Thr95, individually and together. In aqueous solution, all four forms of the peptide bent in the middle to form a hydrophobic cluster. The phosphorylated variants were stabilized further by electrostatic interactions and formation of β-structures, in agreement with previous spectroscopic data. In simulations performed with the peptide in association with a dimyristoylphosphatidylcholine bilayer, the amphipathic α-helical segment remained stable and membrane-associated, although the degree of penetration was less in the phosphorylated variants, and the tilt of the α-helix with respect to the plane of the membrane also changed significantly with the modifications. The extended segment adjacent to this α-helix represents a putative SH3-ligand and remained exposed to the cytoplasm (and thus accessible to binding partners). The results of these simulations demonstrate how this segment of the protein can act as a molecular switch: an amphipathic α-helical segment of the protein is membrane-associated and presents a subsequent proline-rich segment to the cytoplasm for interaction with other proteins. Phosphorylation of threonyl residues alters the degree of membrane penetration of the α-helix and the accessibility of the proline-rich ligand and can stabilize a β-bend. A bend in this region of 18.5 kDa myelin basic protein suggests that the N- and C-termini of the proteins can interact with different leaflets of the myelin membrane and explain how a single protein can bring them close together.     

An amino acid packing code for α-helical structure and protein design

This work demonstrates that all packing in α-helices can be simplified to repetitive patterns of a single motif: the knob-socket. Using the precision of Voronoi Polyhedra/Delauney Tessellations to identify contacts, the knob-socket is a four-residue tetrahedral motif: a knob residue on one α-helix packs into the three-residue socket on another α-helix. The principle of the knob-socket model relates the packing between levels of protein structure: the intra-helical packing arrangements within secondary structure that permit inter-helix tertiary packing interactions. Within an α-helix, the three-residue sockets arrange residues into a uniform packing lattice. Inter-helix packing results from a definable pattern of interdigitated knob-socket motifs between two α-helices. Furthermore, the knob-socket model classifies three types of sockets: (1) free, favoring only intra-helical packing; (2) filled, favoring inter-helical interactions; and (3) non, disfavoring α-helical structure. The amino acid propensities in these three socket classes essentially represent an amino acid code for structure in α-helical packing. Using this code, we used a novel yet straightforward approach for the design of α-helical structure to validate the knob-socket model. Unique sequences for three peptides were created to produce a predicted amount of α-helical structure: mostly helical, some helical, and no helix. These three peptides were synthesized, and helical content was assessed using CD spectroscopy. The measured α-helicity of each peptide was consistent with the expected predictions. These results and analysis demonstrate that the knob-socket motif functions as the basic unit of packing and presents an intuitive tool to decipher the rules governing packing in protein structure.     

Selection of a buried salt bridge by phage display

The α-helical coiled coil is a valuable folding motif for protein design and engineering. By means of phage display technology, we selected a capable binding partner for one strand of a coiled coil bearing a charged amino acid in a central hydrophobic core position. This procedure resulted in a novel coiled coil pair featuring an opposed Glu-Lys pair arranged staggered within the hydrophobic core of a coiled coil structure. Structural investigation of the selected coiled coil dimer by CD spectroscopy and MD simulations suggest that a buried salt bridge within the hydrophobic core enables the specific dimerization of two peptides.     

Sequence preference of α-helix N-terminal tetrapeptide

The α-helix is the most abundant secondary structure in proteins. Due to the specific i, i+4 hydrogen bond pattern, the two termini have unsatisfied hydrogen bonds, and are less constrained; in order to compensate for this, specific residues are preferred for the terminal positions. However, a naive combination of the statistically-preferred residues for each position may not result in a stable N-terminal helical sequence. In order to provide a set of preferable N-terminal peptides for α-helix design, we have studied the N-terminal tetrapeptide sequence motifs that are favorable for helix formation using statistical analysis and atomistic simulations. A set of tetrapeptide sequences including TEEE and TPEE were found to be favorable motifs. In addition to forming more hydrogen bonds in the helical conformation, the favorable motifs also tended to form more capping boxes. To empirically test our predictions, we obtained 10 peptides with different N-terminal motifs and measured their α-helical content by circular dichroism spectroscopy. The experimental results agreed qualitatively with the statistical and simulation results. Furthermore, some of the suggested preferable tetrapeptide sequences have been successfully applied in de novo protein design.     

Conformational analysis of peptide analogues of silkmoth chorion protein segments using CD, NMR and molecular modelling.

Silkmoth proteins secreted from the follicular cells that surround the oocyte form a large extracellular assembly which is important for protecting and sustaining the structure of the oocyte and the developing embryo. These proteins have been classified into two major families (A and B). Sequence analysis showed conservation of a central domain containing long stretches of six amino acid residue repeats in both families, which have been suggested to be organized in beta-sheet structures. In this work NMR and CD spectra, as well as molecular calculations, have been used to investigate the conformational properties of two synthetic peptides (A and B), analogues of parts of the central domain of silkmoth chorion proteins of the A and B families, respectively. These peptides consist of three tandem repeats of the six-residue basic motif. Analysis of CD spectra of the two peptides in aqueous solutions and mixtures with organic solvents revealed beta-sheet and turn structural elements with a percentage higher than 40%. NOESY spectra at low temperatures (263-273 K) show sequential nOe connectivities (i, i + 1), indicative of a relative flexibility. The presence of HNi-HNi+1 cross-peaks and medium Halphai-HNi+1 connectivities, chemical shift deviations and temperature coefficient data provide, for the first time, experimental evidence that local folded structures around Gly residues occur in peptide segments of chorion proteins in solution. Simulated annealing calculations were used to examine the conformational space of the peptides and to probe the initial steps of amyloid fibril formation in the case of chorion proteins.     

Molecular structures and interactions of repetitive peptides based on wheat glutenin subunits depend on chain length

Synthetic and recombinant peptide models of the central repetitive domain of the high molecular weight subunits of wheat glutenin with different numbers of the consensus repeat motifs PGQGQQ + GYYPTSLQQ (21, 45, 110, and 203 residues long) and a recombinant 58,000-Da relative molecular mass (M(r) 58,000) repetitive peptide from a single subunit (1Dx5) are studied using Fourier transform IR spectroscopy. The spectra of the dry peptides are very similar; at low water contents (<76% relative humidity) there is an increase in beta-sheet structures in all peptides. However, on further hydration the content of beta sheets decrease and more beta turns are observed. The changes during the second step of hydration are very marked in the 21 and 45 residue peptides, but they are less apparent in the longer perfect repeat peptides. In the 110 and 203 residue peptides hydration results in increased contents of intermolecular beta-sheets and less beta-turn formation. In contrast, the beta-turn content of the M(r) 58,000 peptide increases during the second hydration step. The decreased extent of structure changes with increasing chain length indicates that cumulative intermolecular interactions, in particular hydrogen bonds, are an important factor in determining the structures in the solid state. The regularity of the perfect repeat sequences in the 21, 45, 110, and 203 residue peptides may favor the formation of larger stretches of intermolecular beta sheets. In contrast, the M(r) 58,000 peptide contains imperfect repeats (in common with native glutenin subunits), which may limit its ability to form intermolecular beta sheets.     

Kinetic steps for α-helix formation

The kinetics of α-helix formation in polyalanine and polyglycine eicosamers (20-mers) were examined using torsional-coordinate molecular dynamics (MD). Of one hundred fifty-five MD experiments on extended (Ala)₂₀ carried out for 0.5 ns each, 129 (83%) formed a persistent α-helix. In contrast, the extended state of (Gly)₂₀ only formed a right-handed α-helix in two of the 20 MD experiments (10%), and these helices were not as long or as persistent as those of polyalanine. These simulations show helix formation to be a competition between the rates of (a) forming local hydrogen bonds (i.e. hydrogen bonds between any residue i and its i + 2, i + 3, i + 4, or i + 5th neighbor) and (b) forming nonlocal hydrogen bonds (HBs) between residues widely separated in sequence. Local HBs grow rapidly into an α-helix; but nonlocal HBs usually retard helix formation by “trapping” the polymer in irregular, “balled-up” structures. Most trajectories formed some nonlocal HBs, sometimes as many as eight. But, for (Ala)₂₀, most of these eventually rearranged to form local HBs that lead to α-helices. A simple kinetic model describes the rate of converting nonlocal HBs into α-helices. Torsional-coordinate MD speeds folding by eliminating bond and angle degrees of freedom and reducing dynamical friction. Thus, the observed 210 ps half-life for helix formation is likely to be a lower bound on the real rate. However, we believe the sequential steps observed here mirror those of real systems. Proteins 33:343–357, 1998. © 1998 Wiley-Liss, Inc.