Comparison of the structural characteristics of Cu2+-bound and unbound α-syn12 peptide obtained in simulations using different force fields

The effects of Cu²⁺ binding and the utilization of different force fields when modeling the structural characteristics of α-syn12 peptide were investigated. To this end, we performed extensive temperature replica exchange molecular dynamics (T-REMD) simulations on Cu²⁺-bound and unbound α-syn12 peptide using the GROMOS 43A1, OPLS-AA, and AMBER03 force fields. Each replica was run for 300 ns. The structural characteristics of α-syn12 peptide were studied based on backbone dihedral angle distributions, free-energy surfaces obtained with different reaction coordinates, favored conformations, the formation of different Turn structures, and the solvent exposure of the hydrophobic residues. The findings show that AMBER03 prefers to sample helical structures for the unbound α-syn12 peptide and does not sample any β-hairpin structure for the Cu²⁺-bound α-syn12 peptide. In contrast, the central structure of the major conformational clusters for the Cu²⁺-bound and unbound α-syn12 peptide according to simulations performed using the GROMOS 43A1 and OPLS-AA force fields is a β-hairpin with Turn_9-6. Cu²⁺ can also promote the formation of the β-hairpin and increase the solvent exposure of hydrophobic residues, which promotes the aggregation of α-syn12 peptide. This study can help us to understand the mechanisms through which Cu²⁺ participates in the fibrillation of α-syn12 peptide at the atomic level, which in turn represents a step towards elucidating the nosogenesis of Parkinson’s disease.

Peptides containing the sulfonamide junction: Synthesis,structure, and conformation of Z-Tau-Pro-Phe-NHiPr

The taurine (Tau) containing tripeptide derivative Z-Tau-Pro-Phe-NHiPr (1) has been synthesized as suitable sulfonamido-pseudopeptide model to investigate formation and conformational properties of folded secondary structures stabilized by intramolecular H bonds directly involving the sulfonamide junction. In the crystal the pseudopeptide 1 adopts a type I β-turn with the Pro and Phe residues located at the (i + 1) and (i + 2) corner positions, respectively. The turn is stabilized by a 4 → 1 H bond engaging one of the SO₂ oxygen atoms and the isopropylamide NH. In CDCl₃ solution the β-turn folding is accompanied by a γ-turn centered at the Pro and involving a 3 → 1 H bond between the SO₂ and the Phe NH. A comparison of the structural and conformational properties found in 1 with those of the already known sulfonamido-pseudopeptides, with particular reference to the models containing the Tau-Pro junction, is also reported. © 1997 John Wiley & Sons, Inc. Biopoly 41: 555–567, 1997.     

β-Turn conformations in crystal structures of model peptides containing α,α-Di-n-propylglycine and α,α-Di-n-butylglycine

The crystal state conformations of three peptides containing the α,α-dialkylated residues. α,α-di-n-propylglycine (Dpg) and α,α-di-n-butylglycine (Dbg), have been established by x-ray diffraction. Boc-Ala-Dpg-Alu-OMe (I) and Boc-Ala-Dbg-Ala-OMe (III) adopt distorted type II β-turn conformations with Ala (1) and Dpg/Dbg (2) as the corner residues. In both peptides the conformational angles at the Dxg residue (I: ? = 66.2°, ψ = 19.3°; III: ? = 66.5°. ψ = 21.1°) deviate appreciably from ideal values for the i + 2 residue in a type II β-turn. In both peptides the observed (N…O) distances between the Boc CO and Ala (3) NH groups are far too long (1: 3.44 Å: III: 3.63 Å) for an intramolecular 4 → 1 hydrogen bond. Boc-Ala-Dpg-Ata-NHMe (II) crystallizes with two independent molecules in the asymmetric unit. Both molecules HA and HB adopt consecutive β-turn (type III-III in HA and type III-I in IIB) or incipient 3₁₀-helical structures, stabilized by two intramolecular 4 → 1 hydrogen bonds. In all four molecules the bond angle N-C^α-C′ (τ) at the Dxg residues are ≥ 110°. The observation of conformational angles in the helical region of ?,ψ space at these residues is consistent with theoretical predictions. © 1995 John Wiley & Sons, Inc.     

Understanding the folding and stability of a designed WW domain protein with replica exchange molecular dynamics simulations

WW domain proteins are usually regarded as simple models for understanding the folding mechanism of β-sheet. CC45 is an artificial protein that is capable of folding into the same structure as WW domain. In this article, the replica exchange molecular dynamics simulations are performed to investigate the folding mechanism of CC45. The analysis of thermal stability shows that β-hairpin 1 is more stable than β-hairpin 2 during the unfolding process. Free energy analysis shows that the unfolding of this protein substantially proceeds through solvating the smaller β-hairpin 2, followed by the unfolding of β-hairpin 1. We further propose the unfolding process of CC45 and the folding mechanism of two β-hairpins. These results are similar to the previous folding studies of formin binding protein 28 (FBP28). Compared with FBP28, it is found that CC45 has more aromatic residues in N-terminal loop, and these residues contact with C-terminal loop to form the outer hydrophobic core, which increases the stability of CC45. Knowledge about the stability and folding behaviour of CC45 may help in understanding the folding mechanisms of the β-sheet and in designing new WW domains.     

Analysis and prediction of the packing of α-helices against a β-sheet in the tertiary structure of globular proteins

The packing of α-helices and β-sheets in six $\text{α}\text{β}$ proteins (e.g. flavodoxin) has been analysed. The results provide the basis for a computer algorithm to predict the tertiary structure of an $\text{α}\text{β}$ protein from its amino acid sequence and actual assignment of secondary structure.The packing of an individual α-helix against a β-sheet generally involves two adjacent ± 4 rows of non-polar residues on the α-helix at the positions i, i + 4, i + 8, i + 1, i + 5, i + 9. The pattern of interacting β-sheet residues results from the twisted nature of the sheet surface and the attendant rotation of the side-chains. At a more detailed level, four of the α-helical residues (i + 1, i + 4, i + 5 and i + 8) form a diamond that surrounds one particular β-sheet residue, generally isoleucine, leucine or valine. In general, the α-helix sits 10 Å above the sheet and lies parallel to the strand direction.The prediction follows a combinational approach. First, a list of possible β-sheet structures (10⁶ to 10¹⁴) is constructed by the generation of all β-sheet topologies and β-strand alignments. This list is reduced by constraints on topology and the location of non-polar residues to mediate the sheet/helix packing, and then rank-ordered on the extent of hydrogen bonding. This algorithm was uniformly applied to 16 $\text{α}\text{β}$ domains in 13 proteins. For every structure, one member of the reduced list was close to the crystal structure; the root-mean-square deviation between equivalenced C^α atoms averaged 5.6 Å for 100 residues. For the $\text{α}\text{β}$ proteins with pure parallel β-sheets, the total number of structures comparable to or better than the native in terms of hydrogen bonds was between 1 and 148. For proteins with mixed β-sheets, the worst case is glyceraldehyde-3-phosphate dehydrogenase, where as many as 3800 structures would have to be sampled. The evolutionary significance of these results as well as the potential use of a combinatorial approach to the protein folding problem are discussed.     

Conformational interconversions in peptide β-turns: Discrimination between enantiomeric conformations by chiral perturbation

The crystal structure of the model tripeptide Boc-Aib-Gly-Leu-OMe ( 1 ) reveals two independent molecules in the asymmetric unit that adopt “enantiomeric” type I and type I′ β-turn conformations with the Aib and Gly residues occupying the corner (i + 1 and i + 2) positions. ¹³C cross polarization and magic angle sample spinning spectra in the solid state also support the coexistence of two conformational species. ¹³C-nmr in CDCl₃ establishes the presence of a single species or rapid exchange between conformations. 400 MHz ¹H-nmr provides evidence for conformational exchange involving a major and minor species, with β-turn conformations supported by the low solvent exposure of Leu(3) NH and the observation of N_iH ↔ N_i+1H nuclear Overhauser effects. CD bands in the region 190–230 nm are positive, supporting a major population of type I′ β-turns. The isomeric peptide, Boc-Gly-Leu-Aib-OMe ( 2 ), adopts an “open” type II′ β-turn conformation in crystals. Solid state and solution nmr support population of a single conformational species. Chiral perturbation introduced outside the folded region of peptides may provide a means of modulating screw sense in achiral sequences. © 1998 John Wiley & Sons, Inc. Biopoly 45: 191–202, 1998     

Role of different β-turns in β-hairpin conformation and stability studied by optical spectroscopy

Model β-hairpin peptides based on variations in the turn sequence of Cochran's tryptophan zipper peptide, SWTWENGKWTWK, were studied using electronic circular dichroism (ECD), fluorescence, and infrared (IR) spectroscopies. The trpzip2 Asn-Gly turn sequence was substituted with Thr-Gly, Aib-Gly, (D)Pro-Gly, and Gly-Asn (trpzip1) to study the impact of turn stability on β-hairpin formation. Stability and conformational changes of these hairpins were monitored by thermodynamic analyses of the temperature variation of both FTIR (amide I') and ECD spectral intensities. These changes were fit to a two-state model which yielded different T(m) values, representing the folding/unfolding process, for hairpins with different β-turns. Different β-turns show systematic contributions to hairpin structure formation, and their inclusion in hairpin design can modify the folding pathways. Aib-Gly or (D)Pro-Gly sequences stabilize the turn resulting in residual Trp-Trp interaction at high temperatures, but at the same time the β-structure (cross strand H-bonds) can become less stable due to constraints of the turn, as seen for (D)Pro-Gly. The structure of the Aib-Gly turn containing hairpin was determined by NMR and was shown to be like trpzip2 (Asn-Gly turn) as regards turn and strand geometries, but to differ from trpzip1 (Gly-Asn turn). The Munoz and Eaton statistical mechanically derived multistate model, tested as an alternate point of view, represented contributions from H-bonds and hydrophobic interactions as well as conformational change as interdependent. Use of different spectral methods that vary in dependence on these physical interactions along with the structural variations provided insight to the complex folding pathways of these small, well-folded peptides.     

Synthesis and β-sheet propensity of constrained N-amino peptides

The stabilization of β-sheet secondary structure through peptide backbone modification represents an attractive approach to protein mimicry. Here, we present strategies toward stable β-hairpin folds based on peptide strand N-amination. Novel pyrazolidinone and tetrahydropyridazinone dipeptide constraints were introduced via on-resin Mitsunobu cyclization between α-hydrazino acid residues and a serine or homoserine side chain. Acyclic and cyclic N-amino peptide building blocks were then evaluated for their effect on β-hairpin stability in water using a GB1-derived model system. Our results demonstrate the strong β-sheet stabilizing effect of the peptide N-amino substituent, and provide useful insights into the impact of covalent dipeptide constraint on β-sheet folding.     

Residue-specific α-helix propensities from molecular simulation

Formation of α-helices is a fundamental process in protein folding and assembly. By studying helix formation in molecular simulations of a series of alanine-based peptides, we obtain the temperature-dependent α-helix propensities of all 20 naturally occurring residues with two recent additive force fields, Amber ff03w and Amber ff99SB¹. Encouragingly, we find that the overall helix propensity of many residues is captured well by both energy functions, with Amber ff99SB¹ being more accurate. Nonetheless, there are some residues that deviate considerably from experiment, which can be attributed to two aspects of the energy function: i), variations of the charge model used to determine the atomic partial charges, with residues whose backbone charges differ most from alanine tending to have the largest error; ii), side-chain torsion potentials, as illustrated by the effect of modifications to the torsion angles of I, L, D, N. We find that constrained refitting of residue charges for charged residues in Amber ff99SB¹ significantly improves their helix propensity. The resulting parameters should more faithfully reproduce helix propensities in simulations of protein folding and disordered proteins.     

Conformational features of a peptide model Ac-DTVKLMYKGQPMTFR-NH2, corresponding to an early folding beta hairpin region of staphylococcal nuclease.

Recent H-D exchange 1H NMR studies of the refolding of Staphylococcal nuclease (P117G) variant suggest that, a region of the protein corresponding to a beta hairpin in the native structure folded early in the refolding process. In order to investigate whether the formation of beta hairpin is an early folding event, we investigated the conformational features of the beta hairpin peptide model Ac-DTVKLMYKGQPMTFR-NH2 from Staphylococcal nuclease with 1H NMR techniques. It appears that the peptide aggregates even at a low concentration. However, based on the observation of weak dnn(i, i + 1) NOEs between K8-G9, G9-Q10, an upfield shift of Gly9 NH and a low temperature coefficient (-d delta/dT) for Gly9 NH, we suggest that the sequence YKGQP as part of the beta hairpin peptide model samples conformational forms with reduced conformational entropy and turn potential. The presence of aggregation could be restricting the population of folded conformational forms and formation of beta hairpin at detectable concentrations. We suggest that, formation of beta hairpin could be an early event in the folding of Staphylococcal nuclease and this observation correlates with H-D exchange 1H NMR results and also with the prediction of a protein folding model proposed in literature.     

The roles of the conserved tyrosine in the β2-α2 loop of the prion protein

ABSTRACT

Prions cause neurodegenerative diseases for which no cure exists. Despite decades of research activities the function of the prion protein (PrP) in mammalians is not known. Moreover, little is known on the molecular mechanisms of the self-assembly of the PrP from its monomeric state (cellular PrP, PrP^C) to the multimeric state. The latter state includes the toxic species (scrapie PrP, PrP^Sc) knowledge of which would facilitate the development of drugs against prion diseases. Here we analyze the role of a tyrosine residue (Y169) which is strictly conserved in mammalian PrPs. Nuclear magnetic resonance (NMR) spectroscopy studies of many mammalian PrP^C proteins have provided evidence of a conformational equilibrium between a 3₁₀-helical turn and a type I β turn conformation in the β2-α2 loop (residues 165–175). In vitro cell-free experiments of the seeded conversion of PrP^C indicate that non-aromatic residues at position 169 reduce the formation of proteinase K-resistant PrP. Recent molecular dynamics (MD) simulations of monomeric PrP and several single-point mutants show that Y169 stabilizes the 3₁₀-helical turn conformation more than single-point mutants at position 169 or residues in contact with it. In the 3₁₀-helical turn conformation the hydrophobic and aggregation-prone segment 169-YSNQNNF-175 is buried and thus not-available for self-assembly. From the combined analysis of simulation and experimental results it emerges that Y169 is an aggregation gatekeeper with a twofold role. Mutations related to 3 human prion diseases are interpreted on the basis of the gatekeeper role in the monomeric state. Another potential role of the Y169 side chain is the stabilization of the ordered aggregates, i.e., reduction of frangibility of filamentous protofibrils and fibrils, which is likely to reduce the generation of toxic species.     

NMR investigations of structural and dynamics features of natively unstructured drug peptide – salmon calcitonin: implication to rational design of potent sCT analogs

Backbone dynamics and conformational properties of drug peptide salmon calcitonin have been studied in aqueous solution using nuclear magnetic resonance (NMR). Although salmon calcitonin (sCT) is largely unfolded in solution (as has been reported in several circular dichroism studies), the secondary H^α chemical shifts and three bond H^N–H^α coupling constants indicated that most of the residues of the peptide are populating the α‐helical region of the Ramachandran (?, ψ) map. Further, the peptide in solution has been found to exhibit multiple conformational states exchanging slowly on the NMR timescale (10²–10³ s^?1), inferred by the multiple chemical shift assignments in the region Leu4–Leu12 and around Pro23 (for residues Gln20–Tyr22 and Arg24). Possibly, these slowly exchanging multiple conformational states might inhibit symmetric self‐association of the peptide and, in part, may account for its reduced aggregation propensity compared with human calcitonin (which lacks this property). The ¹⁵N NMR‐relaxation data revealed (i) the presence of slow (microsecond‐to‐millisecond) timescale dynamics in the N‐terminal region (Cys1–Ser5) and core residues His17 and Asn26 and (ii) the presence of high frequency (nanosecond‐to‐picosecond) motions in the C‐terminal arm. Put together, the various results suggested that (i) the flexible C‐terminal of sCT (from Thr25–Thr31) is involved in identification of specific target receptors, (ii) whereas the N‐terminal of sCT (from Cys1–Gln20) in solution – exhibiting significant amount of conformational plasticity and strong bias towards biologically active α‐helical structure – facilitates favorable conformational adaptations while interacting with the intermembrane domains of these target receptors. Thus, we believe that the structural and dynamics features of sCT presented here will be useful guiding attributes for the rational design of biologically active sCT analogs. Copyright © 2012 European Peptide Society and John Wiley & Sons, Ltd.     

Molecular dynamics simulations of certain mutant peptide models from staphylococcal nuclease reveal that initial hydrophobic collapse associated with turn propensity drives β‐hairpin folding

An important nucleation event during the folding of staphylococcal nuclease involves the formation of a β‐hairpin by the sequence ²¹DTVKLMYKGQPMTFR³⁵. Earlier studies show that the turn sequence ‘YKGQP’ has an important role in the folding of this β‐hairpin. To understand the active or passive nature of the turn sequence ‘YKGQP’ in the folding of the aforementioned β‐hairpin sequence, we studied glycine mutant peptides Ac‐²DTVKLMYGGQPMTFR¹⁶‐NMe (K9G:15), Ac‐²DTVKLMYKGGPMTFR¹⁶‐NMe (Q11G:15), Ac‐²DTVKLMYGGGPMTFR¹⁶‐NMe (K9G/Q11G:15), and Ac‐²DTVKLMGGGGGMTFR¹⁶‐NMe (penta‐G:15) by using molecular dynamics simulations, starting with two different unfolded states, polyproline II and extended conformational forms. Further, 5mer mutant turn peptides Ac‐²YGGQP⁶‐NMe (K3G:5), Ac‐²YKGGP⁶‐NMe (Q5G:5), Ac‐²YGGGP⁶‐NMe (K3G/Q5G:5), and Ac‐²GGGGG⁶‐NMe (penta‐G:5) were also studied individually. Our results show that an initial hydrophobic collapse and loop closure occurs in all 15mer mutants, but only K9G:15 mutant forms a stable native‐like β‐hairpin. In the other 15mer mutants, the hydrophobic collapsed state would not proceed to β‐hairpin formation. Of the different simulations performed for the penta‐G:15 mutant, in only one simulation a nonnative β‐hairpin conformation is sampled with highly flexible loop region (⁸GGGGG¹²), which has no specific conformational preference as a 5mer. While the sequence ‘YGGQP’ in the K3G:5 simulation shows relatively higher β‐turn propensity, the presence of this sequence in K9G:15 peptide seems to be driving the β‐hairpin formation. Thus, these results seem to suggest that for the formation of a stable β‐hairpin, the initial hydrophobic collapse is to be assisted by a turn propensity. Initial hydrophobic collapse alone is not sufficient to guide β‐hairpin formation. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Molecular Origin of Gerstmann-Sträussler-Scheinker Syndrome: Insight from Computer Simulation of an Amyloidogenic Prion Peptide

Prion proteins become pathogenic through misfolding. Here, we characterize the folding of a peptide consisting of residues 109–122 of the Syrian hamster prion protein (the H1 peptide) and of a more amyloidogenic A117V point mutant that leads in humans to an inheritable form of the Gerstmann-Sträussler-Scheinker syndrome. Atomistic molecular dynamics simulations are performed for 2.5 μs. Both peptides lose their α-helical starting conformations and assume a β-hairpin that is structurally similar in both systems. In each simulation several unfolding/refolding events occur, leading to convergence of the thermodynamics of the conformational states to within 1 kJ/mol. The similar stability of the β-hairpin relative to the unfolded state is observed in the two peptides. However, substantial differences are found between the two unfolded states. A local minimum is found within the free energy unfolded basin of the A117V mutant populated by misfolded collapsed conformations of comparable stability to the β-hairpin state, consistent with increased amyloidogenicity. This population, in which V117 stabilizes a hydrophobic core, is absent in the wild-type peptide. These results are supported by simulations of oligomers showing a slightly higher stability of the associated structures and a lower barrier to association for the mutated peptide. Hence, a single point mutation carrying only two additional methyl groups is here shown to be responsible for rather dramatic differences of structuring within the unfolded (misfolded) state.     

Equilibrium folding and unfolding pathways for a model protein

The folding–unfolding process of reduced bovine pancreatic trypsin inhibitor was investigated with an idealized model employing approximate free energies. The protein is regarded to consist of only C^α and C^β atoms. The backbone dihedral angles are the only conformational variables and are permitted to take discrete values at every 10°. Intraresidue energies consist of two terms: an empirical part taken from the observed frequency distributions of (?,ψ) and an additional favorable energy assigned to the native conformation of each residue. Interresidue interactions are simplified by assuming that there is an attractive energy operative only between residue pairs in close contact in the native structure. A total of 230,000 molecular conformations, with no atomic overlaps, ranging from the native state to the denatured state, are randomly generated by changing the sampling bias. Each conformation is classified according to its conformational energy, F; a conformational entropy, S(F) is estimated for each value of F from the number of samples. The dependence of S(F) on energy reveals that the folding–unfolding transition for this idealized model is an “all-or-none” type; this is attributable to the specific long-range interactions. Interresidue contact probabilities, averaged over samples representing various stages of folding, serve to characterize folding intermediates. Most probable equilibrium pathways for the folding–unfolding transition are constructed by connecting conformationally similar intermediates. The specific details obtained for bovine pancreatic trypsin inhibitor are as follows: (1) Folding begins with the appearance of nativelike medium-range contacts at a β-turn and at the α-helix. (2) These grow to include the native pair of interacting β-strands. This state includes intact regular secondary conformations, as well as the interstrand sheet contacts, and corresponds to an activated state with the highest free energy on the pathway. (3) Additional native long-range contacts are completely formed either toward the amino terminus or toward the carboxyl terminus. (4) In a final step, the missing contacts appear. Although these folding pathways for this model are not consistent with experimental reports, it does indicate multiple folding pathways. The method is general and can be applied to any set of calculated conformational energies and furthermore permits investigation of gross folding features.     

Stable and metastable states of human amylin in solution

Patients with type II diabetes exhibit fibrillar deposits of human amylin protein in the pancreas. It has been proposed that amylin oligomers arising along the aggregation or fibril-formation pathways are important in the genesis of the disease. In a step toward understanding these aggregation pathways, in this work we report the conformational preferences of human amylin monomer in solution using molecular simulations and infrared experiments. In particular, we identify a stable conformer that could play a key role in aggregation. We find that amylin adopts three stable conformations: one with an α-helical segment comprising residues 9-17 and a short antiparallel β-sheet comprising residues 24-28 and 31-35; one with an extended antiparallel β-hairpin with the turn region comprising residues 20-23; and one with no particular structure. Using detailed calculations, we determine the relative stability of these various conformations, finding that the β-hairpin conformation is the most stable, followed by the α-helical conformation, and then the unstructured coil. To test our predicted structure, we calculate its infrared spectrum in the amide I stretch regime, which is sensitive to secondary structure through vibrational couplings and linewidths, and compare it to experiment. We find that theoretically predicted spectra are in good agreement with the experimental line shapes presented herein. The implications of the monomer secondary structures on its aggregation pathway and on its interaction with cell membranes are discussed.     

Short-range conformational energies,secondary structure propensities,and recognition of correct sequence-structure matches

A statistical analysis of known structures is made for an assessment of the utility of short-range energy considerations. For each type of amino acid, the potentials governing (1) the torsions and bond angle changes of virtual C^α-C^α bonds and (2) the coupling between torsion and bond angle changes are derived. These contribute approximately −2 RT per residue to the stability of native proteins, approximately half of which is due to coupling effects. The torsional potentials for the α-helical states of different residues are verified to be strongly correlated with the free-energy change measurements made upon single-site mutations at solvent-exposed regions. Likewise, a satisfactory correlation is shown between the β-sheet potentials of different amino acids and the scales from free-energy measurements, despite the role of tertiary context in stabilizing β-sheets. Furthermore, there is excellent agreement between our residue-specific potentials for α-helical state and other thermodynamic based scales. Threading experiments performed by using an inverse folding protocol show that 50 of 62 test structures correctly recognize their native sequence on the basis of short-range potentials. The performance is improved to 55, upon simultaneous consideration of short-range potentials and the nonbonded interaction potentials between sequentially distant residues. Interactions between near residues along the primary structure, i.e., the local or short-range interactions, are known to be insufficient, alone, for understanding the tertiary structural preferences of proteins alone. Yet, knowledge of short-range conformational potentials permits rationalizing the secondary structure propensities and aids in the discrimination between correct and incorrect tertiary folds. Proteins 29:292–308, 1997. © 1997 Wiley-Liss, Inc.     

Importance of context in protein folding: secondary structural propensities versus tertiary contact-assisted secondary structure formation

Molecular dynamics simulations can be used to reveal the detailed conformational behaviors of peptides and proteins. By comparing fragment and full-length protein simulations, we can investigate the role of each peptide segment in the folding process. Here, we take advantage of information regarding the helix formation process from our previous simulations of barnase and protein A as well as new simulations of four helical fragments from these proteins at three different temperatures, starting with both helical and extended structures. Segments with high helical propensity began the folding process by tethering the chain through side chain interactions involving either polar interactions, such as salt bridges, or hydrophobic staples. These tethers were frequently nonnative (i.e., not i --> i + 4 spacing) and provided a scaffold for other residues, thereby limiting the conformational search. The helical structure then propagated on both sides of the tether. Segments with low stability and propensity formed later in the folding process and utilized contacts with other portions of the protein when folding. These helices formed via a tertiary contact-assisted mechanism, primarily via hydrophobic contacts between residues distant in sequence. Thus, segments with different helical propensities appear to play different roles during protein folding. Furthermore, the active role of nonlocal side chains in helix formation highlights why we must move beyond simple hierarchical models of protein folding.     

Structural versatility of peptides from Cα, α-disubstituted glycines: Crystal-state conformational analysis of homopeptides from Cα-methyl,Cα-benzylglycine [(αMe)Phe]n

The molecular and crystal structures of one derivative and three homopeptides (from the di-to the tetrapeptide level) of the chiral, C^{α, α}-disubstituted glycine C^α-methyl, C^α-benzylglycine [(αMe)Phe], have been determined by x-ray diffraction. The derivative is mClAc-D -(αMe)Phe-OH, and the peptides are pBrBz-[D -(αMe)Phe]₂-NHMe, pBrBz-[D -(αMe)Phe]₃-OH hemihydrate, and pBrBz-[D -(αMe)Phe]₄-OtBu sesquihydrate. All (αMe)Phe residues prefer ?,ψ torsion angles in the helical region of the conformational map. The dipeptide methylamide and the tripeptide carboxylic acid adopt a β-turn conformation with a 1 ← 4 C?O…?H? N intramolecular H bond. The structure of the tripeptide carboxylic acid is further stabilized by a 1 ← 4 C?O…?H? O intramolecular H bond, forming an “oxy-analogue” of a β-turn. The tetrapeptide ester is folded in a regular (incipient) 3₁₀-helix. In general, the relationship between (αMe)Phe chirality and helix screw sense is opposite to that exhibited by protein amino acids. A comparison is made with the conclusions extracted from published work on homopeptides from other C^α-methylated α-amino acids. © 1993 John Wiley & Sons, Inc.