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.     

β-Hairpin-Mediated Formation of Structurally Distinct Multimers of Neurotoxic Prion Peptides

Protein misfolding disorders are associated with conformational changes in specific proteins, leading to the formation of potentially neurotoxic amyloid fibrils. During pathogenesis of prion disease, the prion protein misfolds into β-sheet rich, protease-resistant isoforms. A key, hydrophobic domain within the prion protein, comprising residues 109–122, recapitulates many properties of the full protein, such as helix-to-sheet structural transition, formation of fibrils and cytotoxicity of the misfolded isoform. Using all-atom, molecular simulations, it is demonstrated that the monomeric 109–122 peptide has a preference for α-helical conformations, but that this peptide can also form β-hairpin structures resulting from turns around specific glycine residues of the peptide. Altering a single amino acid within the 109–122 peptide (A₁₁₇V, associated with familial prion disease) increases the prevalence of β-hairpin formation and these observations are replicated in a longer peptide, comprising residues 106–126. Multi-molecule simulations of aggregation yield different assemblies of peptide molecules composed of conformationally-distinct monomer units. Small molecular assemblies, consistent with oligomers, comprise peptide monomers in a β-hairpin-like conformation and in many simulations appear to exist only transiently. Conversely, larger assemblies are comprised of extended peptides in predominately antiparallel β-sheets and are stable relative to the length of the simulations. These larger assemblies are consistent with amyloid fibrils, show cross-β structure and can form through elongation of monomer units within pre-existing oligomers. In some simulations, assemblies containing both β-hairpin and linear peptides are evident. Thus, in this work oligomers are on pathway to fibril formation and a preference for β-hairpin structure should enhance oligomer formation whilst inhibiting maturation into fibrils. These simulations provide an important new atomic-level model for the formation of oligomers and fibrils of the prion protein and suggest that stabilization of β-hairpin structure may enhance cellular toxicity by altering the balance between oligomeric and fibrillar protein assemblies.     

Structural and thermodynamics characters of isolated α-syn12 peptide: long-time temperature replica-exchange molecular dynamics in aqueous solution

The structural and thermodynamics characters of α-syn12 (residues 1-12 of the human α-synuclein protein) peptide in aqueous solution were investigated through temperature replica-exchange molecular dynamics (T-REMD) simulations with the GROMOS 43A1 force field. The two independent T-REMD simulations were completed starting from an initial conformational α-helix and an irregular structure, respectively. Each replica was run for 300 ns. The structural and thermodynamics characters were studied based on parameters such as distributions of backbone dihedral angles, free energy surface, stability of folded β-hairpin structure, and favorite conformations. The results showed that the isolated α-syn12 peptide in water adopted four different conformational states: the first state was a β-hairpin ensemble with Turn(9-6) and four hydrogen bonds, the second state was a β-hairpin ensemble with two turns (Turn(9-6) and Turn(5-2)) and three hydrogen bonds, the third state was a disordered structure with both Turn(8-5) and Turn(5-2), and the last state was a π-helix ensemble. Meanwhile, we studied the free energy change of α-syn12 peptide from the unfolded state to the β-hairpin state, which was in good agreement with the experiments and molecular dynamics simulations for some other peptides. We also analyzed the driving force of the peptide transition. The results indicated that the driving forces were high solvent exposure of hydrophobic Leu8 and hydrophobic residues in secondary structure. To our knowledge, this was the first report to study the isolated α-syn12 peptide in water by T-REMD.     

Flow-Induced β-Hairpin Folding of the Glycoprotein Ibα β-Switch

Flow-induced shear has been identified as a regulatory driving force in blood clotting. Shear induces β-hairpin folding of the glycoprotein Ibα β-switch which increases affinity for binding to the von Willebrand factor, a key step in blood clot formation and wound healing. Through 2.1-μs molecular dynamics simulations, we investigate the kinetics of flow-induced β-hairpin folding. Simulations sampling different flow velocities reveal that under flow, β-hairpin folding is initiated by hydrophobic collapse, followed by interstrand hydrogen-bond formation and turn formation. Adaptive biasing force simulations are employed to determine the free energy required for extending the unfolded β-switch from a loop to an elongated state. Lattice and freely jointed chain models illustrate how the folding rate depends on the entropic and enthalpic energy, the latter controlled by flow. The results reveal that the free energy landscape of the β-switch has two stable conformations imprinted on it, namely, loop and hairpin—with flow inducing a transition between the two.     

Nascent β-Hairpin Formation of a Natively Unfolded Peptide Reveals the Role of Hydrophobic Contacts

Despite the important role of the unfolded states in protein stability, folding, and aggregation, they remain poorly understood due to the lack of residue-specific experimental data. Here, we explore features of the unfolded state of the NTL9 protein by applying all-atom replica-exchange simulations to the two fragment peptides NTL9(1–22) and NTL9(6–17). We found that while NTL9(6–17) is unstructured, NTL9(1–22) transiently folds as various β-hairpins, a fraction of which contain a native β-sheet. Interestingly, despite a large number of charged residues, the formation of backbone hydrogen bonds is concomitant with hydrophobic but not electrostatic contacts. Although the fragment peptides lack a proposed specific contact between Asp⁸ and Lys¹², the individually weak, nonspecific interactions with lysines together stabilize the charged Asp⁸, leading to a pK_a shift of nearly 0.5 units, in agreement with the NMR data. Taken together, our data suggest that the unfolded state of NTL9 likely contains a β-hairpin in segment 1–22 with sequence-distant hydrophobic contacts, thus lending support to a long-standing hypothesis that the unfolded states of proteins exhibit native-like topology with hydrophobic clusters.     

Highly Disordered Amyloid-β Monomer Probed by Single-Molecule FRET and MD Simulation

Monomers of amyloid-β (Aβ) protein are known to be disordered, but there is considerable controversy over the existence of residual or transient conformations that can potentially promote oligomerization and fibril formation. We employed single-molecule Förster resonance energy transfer (FRET) spectroscopy with site-specific dye labeling using an unnatural amino acid and molecular dynamics simulations to investigate conformations and dynamics of Aβ isoforms with 40 (Aβ40) and 42 residues (Aβ42). The FRET efficiency distributions of both proteins measured in phosphate-buffered saline at room temperature show a single peak with very similar FRET efficiencies, indicating there is apparently only one state. 2D FRET efficiency-donor lifetime analysis reveals, however, that there is a broad distribution of rapidly interconverting conformations. Using nanosecond fluorescence correlation spectroscopy, we measured the timescale of the fluctuations between these conformations to be ～35 ns, similar to that of disordered proteins. These results suggest that both Aβ40 and Aβ42 populate an ensemble of rapidly reconfiguring unfolded states, with no long-lived conformational state distinguishable from that of the disordered ensemble. To gain molecular-level insights into these observations, we performed molecular dynamics simulations with a force field optimized to describe disordered proteins. We find, as in experiments, that both peptides populate configurations consistent with random polymer chains, with the vast majority of conformations lacking significant secondary structure, giving rise to very similar ensemble-averaged FRET efficiencies.     

Effect of the aminoacid composition of model α-helical peptides on the physical properties of lipid bilayers and peptide conformation: a molecular dynamics simulation

The interaction of a model Lys flanked α-helical peptides K₂-X₂₄-K₂, (X = A,I,L,L+A,V) with lipid bilayers composed of dimyristoylphosphatidylcholine (DMPC) and dipalmitoylphosphatidylcholine (DPPC) both, in a gel and in a liquid-crystalline state, has been studied by molecular dynamics simulations. It has been shown that these peptides cause disordering of the lipid bilayer in the gel state but only small changes have been monitored in a liquid-crystalline state. The peptides affect ordering of the surrounding lipids depending on the helix stability which is determined by amino acid side chains – their volume, shape, etc. We have shown that the helix does not keep the linear shape in all simulations but often bends or breaks. During some simulations with a very small difference between hydrophobic length of peptide and membrane thickness the peptide exhibits negligible tilt. At the same time changes in peptide conformations during simulations resulted in appearance of superhelix.     

Exploring the free energy landscape of a model β‐hairpin peptide and its isoform

Secondary structural transitions from α‐helix to β‐sheet conformations are observed in several misfolding diseases including Alzheimer's and Parkinson's. Determining factors contributing favorably to the formation of each of these secondary structures is therefore essential to better understand these disease states. β‐hairpin peptides form basic components of anti‐parallel β‐sheets and are suitable model systems for characterizing the fundamental forces stabilizing β‐sheets in fibrillar structures. In this study, we explore the free energy landscape of the model β‐hairpin peptide GB1 and its E2 isoform that preferentially adopts α‐helical conformations at ambient conditions. Umbrella sampling simulations using all‐atom models and explicit solvent are performed over a large range of end‐to‐end distances. Our results show the strong preference of GB1 and the E2 isoform for β‐hairpin and α‐helical conformations, respectively, consistent with previous studies. We show that the unfolded states of GB1 are largely populated by misfolded β‐hairpin structures which differ from each other in the position of the β‐turn. We discuss the energetic factors contributing favorably to the formation of α‐helix and β‐hairpin conformations in these peptides and highlight the energetic role of hydrogen bonds and non‐bonded interactions. Proteins 2014; 82:2394–2402. © 2014 Wiley Periodicals, Inc.     

Chain length effects on helix-hairpin distribution in short peptides with Aib-DAla and Aib-Aib segments

The Aib-D Ala dipeptide segment has a tendency to form both type-I'/III' and type-I/III β-turns. The occurrence of prime turns facilitates the formation of β-hairpin conformations, while type-I/III turns can nucleate helix formation. The octapeptide Boc-Leu-Phe-Val-Aib-DAla-Leu-Phe-Val-OMe (1) has been previously shown to form a β-hairpin in the crystalline state and in solution. The effects of sequence truncation have been examined using the model peptides Boc-Phe-Val-Aib-Xxx-Leu-Phe-NHMe (2, 6), Boc-Val-Aib-Xxx-Leu-NHMe (3, 7), and Boc-Aib-Xxx-NHMe (4, 8), where Xxx=DAla, Aib. For peptides with central Aib-Aib segments, Boc-Phe-Val-Aib-Aib-Leu-Phe-NHMe (6), Boc-Val-Aib-Aib-Leu-NHMe (7), and Boc-Aib-Aib-NHMe (8) helical conformations have been established by NMR studies in both hydrogen bonding (CD3OH) and non-hydrogen bonding (CDCl3) solvents. In contrast, the corresponding hexapeptide Boc-Phe-Val-Aib-DAla-Leu-Phe-Val-NHMe (2) favors helical conformations in CDCl3 and β-hairpin conformations in CD3 OH. The β-turn conformations (type-I'/III) stabilized by intramolecular 4→1 hydrogen bonds are observed for the peptide Boc-Aib-D Ala-NHMe (4) and Boc-Aib-Aib-NHMe (8) in crystals. The tetrapeptide Boc-Val-Aib-Aib-Leu-NHMe (7) adopts an incipient 3(10)-helical conformation stabilized by three 4→1 hydrogen bonds. The peptide Boc-Val-Aib-DAla-Leu-NHMe (3) adopts a novel α-turn conformation, stabilized by three intramolecular hydrogen bonds (two 4→1 and one 5→1). The Aib-DAla segment adopts a type-I' β-turn conformation. The observation of an NOE between Val (1) NH?HNCH3 (5) in CD3OH suggests, that the solid state conformation is maintained in methanol solutions.     

Folding thermodynamics of β‐hairpins studied by replica‐exchange molecular dynamics simulations

We study the differences in folding stability of β‐hairpin peptides, including GB1 hairpin and a point mutant GB1 K10G, as well as tryptophan zippers (TrpZips): TrpZip1, TrpZip2, TrpZip3‐1, and TrpZip4. By performing replica‐exchange molecular dynamics simulations with Amber03* force field (a modified version of Amber ff03) in explicit solvent, we observe ab initio folding of all the peptides except TrpZip3‐1, which is experimentally known to be the least stable among the peptides studied here. By calculating the free energies of unfolding of the peptides at room temperature and folding midpoint temperatures for thermal unfolding of peptides, we find that TrpZip4 and GB1 K10G peptides are the most stable β‐hairpins followed by TrpZip1, GB1, and TrpZip2 in the given order. Hence, the proposed K10G mutation of GB1 peptide results in enhanced stability compared to wild‐type GB1. An important goal of our study is to test whether simulations with Amber 03* model can reproduce experimentally predicted folding stability differences between these peptides. While the stabilities of GB1 and TrpZip1 yield close agreement with experiment, TrpZip2 is found to be less stable than predicted by experiment. However, as heterogenous folding of TrpZip2 may yield divergent thermodynamic parameters by different spectroscopic methods, mismatching of results with previous experimental values are not conclusive of model shortcomings. For most of the cases, molecular simulations with Amber03* can successfully reproduce experimentally known differences between the mutated peptides, further highlighting the predictive capabilities of current state‐of‐the‐art all‐atom protein force fields. Proteins 2015; 83:1307–1315. © 2015 Wiley Periodicals, Inc.     

Direct Observations of Shifts in the β-Sheet Register of a Protein-Peptide Complex Using Explicit Solvent Simulations

Using explicit solvent molecular dynamics simulations, we were able to obtain direct observations of shifts in the hydrogen-bonding register of an intermolecular β-sheet protein-peptide complex. The β-sheet is formed between the FHA domain of cancer marker protein Ki67 (Ki67FHA) and a peptide fragment of the hNIFK signaling protein. Potential encounter complexes of the Ki67FHA receptor and hNIFK peptide are misregistered states of the β-sheet. Rearrangements of one of these misregistered states to the native state were captured in three independent simulations. All three rearrangements occurred by a common mechanism: an aromatic residue of the peptide (F263) anchors into a transient hydrophobic pocket of the receptor to facilitate the formation of native hydrogen bonds. To our knowledge, these simulations provide the first atomically detailed visualizations of a mechanism by which nature might correct for errors in the alignment of intermolecular β-sheets.     

Synthesis and characterization of β‐peptide helices as transmembrane domains in lipid model membranes

Aggregation, orientation and dynamics of transmembrane helices are of relevance for protein function and transmembrane signaling. To explore the interactions of transmembrane helices and the interdependence of peptide structure and lipid composition of the membranes, β‐peptides were explored as model transmembrane domains. Various hydrophobic β‐peptide sequences were synthesized by solid phase peptide synthesis. Conformational analyses of β‐peptide helices were performed in organic solvents (methanol and 2,2,2‐trifluoroethanol) and in large unilamellar liposomes (dimyristoylphosphatidylcholine, dipalmitoylphosphatidylcholine and dioleoylphosphatidylcholine) indicating 12‐ and 14‐helix conformations, depending on β³‐amino acid sequences. The intrinsic tryptophan fluorescence of β³‐homotryptophan units inserted in the center or near the end of the sequence was used to verify the membrane insertion of the β‐peptides. A characteristic blue shift with peripheral β³‐homotryptophan compared with β‐peptides with central tryptophan served as indication for a transmembrane orientation of the β‐peptides within the lipid bilayer. Copyright © 2016 European Peptide Society and John Wiley & Sons, Ltd.     

Mechanism of Membrane Permeation Induced by Synthetic β-Hairpin Peptides

We have investigated the membrane destabilizing properties of synthetic amphiphilic cationic peptides, MAX1 and MAX35, which have the propensity to form β-hairpin structures under certain conditions, and a control non-β-hairpin-forming peptide MAX8V16E. All three peptides bind to liposomes containing a mixture of zwitterionic POPC and negatively charged POPS lipids as determined by Zeta potential measurements. Circular dichroism measurements indicated folding of MAX1 and MAX35 in the presence of the POPC/POPS liposomes, whereas no such folding was observed with MAX8V16E. There was no binding or folding of these peptides to liposomes containing only POPC. MAX1 and MAX35 induced release of contents from negatively charged liposomes, whereas MAX8V16E failed to promote solute release under identical conditions. Thus, MAX1 and MAX35 bind to, and fold at the surface of negatively charged liposomes adopting a lytic conformation. We ruled out leaky fusion as a mechanism of release by including 2 mol % PEG-PE in the liposomes, which inhibits aggregation/fusion but not folding of MAX or MAX-induced leakage. Using a concentration-dependent quenching probe (calcein), we determined that MAX-induced leakage of liposome contents was an all-or-none process. At MAX1 concentrations, which cause release of ∼50% of the liposomes that contain small (R_h <1.5 nm) markers, only ∼15% of those liposomes release a fluorescent dextran of 40 kDa. A multimeric model of the pore is presented based on these results. Atomistic molecular dynamics simulations show that barrels consisting of 10 β-hairpin MAX1 and MAX35 peptides are relatively more stable than MAX8V16E barrels in the bilayer, suggesting that barrels of this size are responsible for the peptides lytic action.     

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.     

All-atom computer simulations of amyloid fibrils disaggregation

Amyloidlike fibrils are found in many fatal diseases, including Alzheimer's disease, type II diabetes mellitus, transmissible spongiform encephalopathies, and prion diseases. These diseases are linked to proteins that have partially unfolded, misfolded, and aggregated into amyloidlike fibrils. The kinetics of amyloidlike fibrils aggregation is still hotly debated and remains an important open question. We have utilized the GNNQQNY crystal structure and high-temperature molecular dynamics simulation in explicit solvent to study the disaggregation mechanism of the GNNQQNY fibrils and to infer its likely aggregation pathways. A hexamer model and a 12-mer model both with two parallel β-sheets separated by a dry side-chain interface were adopted in our computational analysis. A cumulative time of 1 μs was simulated for the hexamer model at five different temperatures (298 K, 348 K, 398 K, 448 K, and 498 K), and a cumulative time of 2.1 μs was simulated for the 12-mer model at four temperatures (298 K, 398 K, 448 K, and 498 K). Our disaggregation landscape and kinetics analyses indicate that tetramers probably act as the transition state in both the hexamer and the 12-mer simulations. In addition, the 12-mer simulations show that the initial aggregation nucleus is with eight peptides. Furthermore, the landscape is rather flat from 8-mers to 12-mers, indicating the absence of major barriers once the initial aggregation nucleus forms. Thus, the likely aggregation pathway is from monomers to the initial nucleus of 8-mers with tetramers as the transition state. Transition state structure analysis shows that the two dominant transition state conformations are tetramers in the 3-1 and 2-2 arrangements. The predominant nucleus conformations are in peptide arrangements maximizing dry side-chain contacts. Landscape and kinetics analyses also indicate that the parallel β-sheets form earlier than the dry side-chain contacts during aggregation. These results provide further insights in understanding the early fibrils aggregation.     

Role of β-Hairpin Formation in Aggregation: The Self-Assembly of the Amyloid-β(25-35) Peptide

The amyloid-β(25-35) peptide plays a key role in the etiology of Alzheimer's disease due to its extreme toxicity even in the absence of aging. Because of its high tendency to aggregate and its low solubility in water, the structure of this peptide is still unknown. In this work, we sought to understand the early stages of aggregation of the amyloid-β(25-35) peptide by conducting simulations of oligomers ranging from monomers to tetramers. Our simulations show that although the monomer preferentially adopts a β-hairpin conformation, larger aggregates have extended structures, and a clear transition from compact β-hairpin conformations to extended β-strand structures occurs between dimers and trimers. Even though β-hairpins are not present in the final architecture of the fibril, our simulations indicate that they play a critical role in fibril growth. Our simulations also show that β-sheet structures are stabilized when a β-hairpin is present at the edge of the sheet. The binding of the hairpin to the sheet leads to a subsequent destabilization of the hairpin, with part of the hairpin backbone dangling in solution. This free section of the peptide can then recruit an extra monomer from solution, leading to further sheet extension. Our simulations indicate that the peptide must possess sufficient conformational flexibility to switch between a hairpin and an extended conformation in order for β-sheet extension to occur, and offer a rationalization for the experimental observation that overstabilizing a hairpin conformation in the monomeric state (for example, through chemical cross-linking) significantly hampers the fibrillization process.     

Minimal folding pathways for coarse-grained biopolymer fragments

The minimal folding pathway or trajectory for a biopolymer can be defined as the transformation that minimizes the total distance traveled between a folded and an unfolded structure. This involves generalizing the usual Euclidean distance from points to one-dimensional objects such as a polymer. We apply this distance here to find minimal folding pathways for several candidate protein fragments, including the helix, the β-hairpin, and a nonplanar structure where chain noncrossing is important. Comparing the distances traveled with root mean-squared distance and mean root-squared distance, we show that chain noncrossing can have large effects on the kinetic proximity of apparently similar conformations. Structures that are aligned to the β-hairpin by minimizing mean root-squared distance, a quantity that closely approximates the true distance for long chains, show globally different orientation than structures aligned by minimizing root mean-squared distance.     

Proline-glutamate chimera’s side chain conformation directs the type of β-hairpin structure

Our aim was to study the impact of two proline chimeras, containing a glutamic acid side chain in cis- or trans-configuration, on secondary structure formation. We further investigated to what extent the configuration of the side chain contributes to the overall peptide conformation. We used a 10 residue peptide (IYSNPDGTWT) that forms a β-hairpin in water. The turn-forming proline was substituted with either a cis- or trans-proline-glutamic acid chimera, resulting in the peptides IYSNP ^{cis -E} DGTWT (P1_P ^cis-E ) and IYSNP ^{trans -E} DGTWT (P1_P ^trans-E ). We studied the conformation of the modified peptides by circular dichroism (CD) and NMR-spectroscopy, and SEC/static light scattering (SLS) analysis. NMR analysis reveals that the modified peptides maintain the β-hairpin conformation in aqueous solution. At 5 °C and pH 4.3, the peptide (P1_P ^cis-E ) was found to adopt two coexisting β-hairpin conformations (2:2 β-hairpin, and 3:5 β-hairpin). In contrast to that, the peptide (P1_P ^trans-E ) adopts a 2:2 β-hairpin that exists in equilibrium with a 4:4 β-hairpin conformation. The adoption of ordered β-hairpin structures for both modified peptides could be confirmed by CD spectroscopy, while SEC/SLS analysis showed a monomeric oligomerization state for all three investigated peptides. With the combination of several NMR methods, we were able to elucidate that even small alterations in the side chain conformation of the proline-glutamate chimera (cis or trans) can significantly influence the conformation of the adopted β-hairpin.     

Structural basis of folding cooperativity in model proteins: insights from a microcanonical perspective

Two-state cooperativity is an important characteristic in protein folding. It is defined by a depletion of states that lie energetically between folded and unfolded conformations. There are different ways to test for two-state cooperativity; however, most of these approaches probe indirect proxies of this depletion. Generalized-ensemble computer simulations allow us to unambiguously identify this transition by a microcanonical analysis on the basis of the density of states. Here, we present a detailed characterization of several helical peptides obtained by coarse-grained simulations. The level of resolution of the coarse-grained model allowed to study realistic structures ranging from small α-helices to a de novo three-helix bundle without biasing the force field toward the native state of the protein. By linking thermodynamic and structural features, we are able to show that whereas short α-helices exhibit two-state cooperativity, the type of transition changes for longer chain lengths because the chain forms multiple helix nucleation sites, stabilizing a significant population of intermediate states. The helix bundle exhibits signs of two-state cooperativity owing to favorable helix-helix interactions, as predicted from theoretical models. A detailed analysis of secondary and tertiary structure formation fits well into the framework of several folding mechanisms and confirms features that up to now have been observed only in lattice models.