Mechanism for the alpha-helix to beta-hairpin transition

The aggregation of alpha-helix-rich proteins into beta-sheet-rich amyloid fibrils is associated with fatal diseases, such as Alzheimer's disease and prion disease. During an aggregation process, protein secondary structure elements-alpha-helices-undergo conformational changes to beta-sheets. The fact that proteins with different sequences and structures undergo a similar transition on aggregation suggests that the sequence nonspecific hydrogen bond interaction among protein backbones is an important factor. We perform molecular dynamics simulations of a polyalanine model, which is an alpha-helix in its native state and observe a metastable beta-hairpin intermediate. Although a beta-hairpin has larger potential energy than an alpha-helix, the entropy of a beta-hairpin is larger because of fewer constraints imposed by the hydrogen bonds. In the vicinity of the transition temperature, we observe the interconversion of the alpha-helix and beta-sheet states via a random coil state. We also study the effect of the environment by varying the relative strength of side-chain interactions for a designed peptide-an alpha-helix in its native state. For a certain range of side-chain interaction strengths, we find that the intermediate beta-hairpin state is destabilized and even disappears, suggesting an important role of the environment in the aggregation propensity of a peptide.     

Assembly of a tetrameric alpha-helical bundle: computer simulations on an intermediate-resolution protein model.

Discontinuous molecular dynamics (DMD) simulation on an intermediate-resolution protein model is used to study the folding of an isolated, small model peptide to an amphipathic alpha-helix and the assembly of four of these model peptides into a four-helix bundle. A total of 129 simulations were performed on the isolated peptide, and 50 simulations were performed on the four-peptide system. Simulations efficiently sample conformational space allowing complete folding trajectories from random initial configurations to be observed within 15 min for the one-peptide system and within 15 h for the four-peptide system on a 500-MHz workstation. The native structures of both the alpha-helix and the four-helix bundle are consistent with experimental characterization studies and with results from previous simulations on these model peptides. In both the one- and four-peptide systems, the native state is achieved during simulations within an optimal temperature range, a phenomenon also observed experimentally. The ease with which our simulations yield reasonable estimates of folded structures demonstrates the power of the intermediate-resolution model developed for this work and the DMD algorithm and suggests that simulations of very long times and of multiprotein systems may be possible with this model.     

One gene, two diseases and three conformations: molecular dynamics simulations of mutants of human prion protein at room temperature and elevated temperatures

Fatal familial insomnia (FFI) and Creutzfeldt-Jakob disease (CJD) are associated to the same mutation at codon 178 but differentiate into clinicopathologically distinct diseases determined by this mutation and a naturally occurring methionine-valine polymorphism at codon 129 of the prion protein gene. It has been suggested that the clinical and pathological difference between FFI and CJD is caused by different conformations of the prion protein. Using molecular dynamics (MD), we investigated the effect of the mutation at codon 178 and the polymorphism at codon 129 on prion protein dynamics and conformation at normal and elevated temperatures. Four model structures were examined with a focus on their dynamics and conformational changes. The results showed differences in stability and dynamics between polymorphic variants. Methionine variants demonstrated a higher stability than valine variants. Elongation of existing beta-sheets and formation of new beta-sheets was found to occur more readily in valine polymorphic variants. We also discovered the inhibitory effect of proline residue on existing beta-sheet elongation.     

Molecular dynamics simulations of alanine rich beta-sheet oligomers: Insight into amyloid formation

The aggregation observed in protein conformational diseases is the outcome of significant new beta-sheet structure not present in the native state. Peptide model systems have been useful in studies of fibril aggregate formation. Experimentally, it was found that a short peptide AGAAAAGA is one of the most highly amyloidogenic peptides. This peptide corresponds to the Syrian hamster prion protein (ShPrP) residues 113-120. The peptide was observed to be conserved in all species for which the PrP sequence has been determined. We have simulated the stabilities of oligomeric AGAAAAGA and AAAAAAAA (A8) by molecular dynamic simulations. Oligomers of both AGAAAAGA and AAAAAAAA were found to be stable when the size is 6 to 8 (hexamer to octamer). Subsequent simulation of an additional alpha-helical AAAAAAAA placed on the A8-octamer surface has revealed molecular events related to conformational change and oligomer growth. Our study addresses both the minimal oligomeric size of an aggregate seed and the mechanism of seed growth. Our simulations of the prion-derived 8-residue amyloidogenic peptide and its variant have indicated that an octamer is stable enough to be a seed and that the driving force for stabilization is the hydrophobic effect.     

β-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.     

Correlations among morphology, beta-sheet stability, and molecular structure in prion peptide aggregates

The misfolding of proteins into beta-sheets and the subsequent aggregation of these sheets into fibrous networks underlies many diseases. In this paper, the role of peptide structure in determining the ordering of beta-sheet aggregates and the morphology of fibrils and protofibrils is dissected. Using a series of peptides based on residues 109-122 of the Syrian hamster prion protein (H1) with a range of substitutions at position 117, the link between side chain interactions and beta-sheet thermal stability has been investigated. The thermal stability of beta-sheets is associated with the peptides' ability to adopt the same alignment as wild-type H1, with residue 117 in register across all beta-strands [Silva, R. A. G. D., Barber-Armstrong, W., and Decatur, S. M. (2003) J. Am. Chem. Soc. 125, 13674-13675]. These aligned strands are capable of forming long, rigid, and twisted fibrils (as visualized by atomic force microscopy) which are thermostable. Peptides which do not adopt this strand alignment aggregate to form thin, flexible, and smooth protofibrils. The ability to form ordered aggregates, and thus to form twisted fibrils, is modulated by the structure of the side chain of residue 117.     

Exploring structural and thermodynamic stabilities of human prion protein pathogenic mutants D202N, E211Q and Q217R

The central event in the pathogenesis of prion protein (PrP) is a profound conformational change from its α-helical (PrP(C)) to its β-sheet-rich isoform (PrP(Sc)). Many single amino acid mutations of PrP are associated with familial prion diseases, such as D202N, E211Q, and Q217R mutations located at the third native α-helix of human PrP. In order to explore the underlying structural and dynamic effects of these mutations, we performed all-atom molecular dynamics (MD) simulations for the wild-type (WT) PrP and its mutants. The obtained results indicate that these amino acid substitutions have subtle effects on the protein structures, but show large changes of the overall electrostatic potential distributions. We can infer that the changes of PrP electrostatic surface due to the studied mutations may influence the intermolecular interactions during the aggregation process. In addition, the mutations also affect the thermodynamic stabilities of PrP.     

Structure and orientation study of fusion peptide FP23 of gp41 from HIV-1 alone or inserted into various lipid membrane models (mono-, bi- and multibi-layers) by FT-IR spectroscopies and Brewster angle microscopy

In the present work, we study the structure and the orientation of the 23 N-terminal peptide of the HIV-1 gp 41 protein (AVGIGALFLGFLGAAGSTMGARS) called FP23. The behaviour of FP23 was investigated alone at the air/water interface and inserted into various lipid model systems: in monolayer or multibilayers of a DOPC/cholesterol/DOPE/DOPG (6/5/3/2) and in a DMPC bilayer. PMIRRAS and polarized ATR spectroscopy coupled with Brewster angle microscopy and spectral simulations were used to precisely determine the structure and the orientation of the peptide in its environment as well as the lipid perturbations induced by the FP23 insertion. The infra-red results show the structural polymorphism of the FP23 and its ability to transit quasi irreversibly from an alpha-helix to antiparallel beta-sheets. At the air/water interface, the transition is induced by compression of the peptide alone and is modulated by compression and lipid to peptide ratio (Ri) when FP23 is inserted into a lipid monolayer. In multibilayers and in a single bilayer, there is coexistence in quasi equal proportions of alpha-helix and antiparallel beta-sheets of FP23 at low peptide content (Ri=100, 200) while antiparallel beta-sheets are predominant at high FP23 concentration (Ri=50). In (multi)bilayer systems, evaluation of dichroic ratios and sprectral simulations show that both the alpha-helix and the antiparallel beta-sheets are tilted at diluted FP23 concentrations (tilt angle of alpha-helix with respect to the normal of the interface=36.5+/-3.0 degrees for FP23 in multibilayers of DOPC/Chol/DOPE/DOPG at Ri=200 and 39.0+/-5.0 degrees in a single bilayer of DMPC at Ri=100 and tilt angle of the beta-sheets=36.0+/-2.0 degrees for the beta-sheets in multibilayers and 30.0+/-2.0 degrees in the lipid bilayer). In parallel, the FP23 induces an increase of the lipid chain disorder which shows both by an increase of the methylene stretching frequencies and an increase of the average C-C-C angle of the acyl chains. At high FP23 content (Ri=50), the antiparallel beta-sheets induce a complete disorganization of the lipid chains in (multi)bilayers.     

Aggregation of PrP106–126 on surfaces of neutral and negatively charged membranes studied by molecular dynamics simulations

Prion diseases are neurodegenerative disorders characterized by the aggregation of an abnormal form of prion protein. The interaction of prion protein and cellular membrane is crucial to elucidate the occurrence and development of prion diseases. Its fragment, residues 106–126, has been proven to maintain the pathological properties of misfolded prion and was used as a model peptide. In this study, explicit solvent molecular dynamics (MD) simulations were carried out to investigate the adsorption, folding and aggregation of PrP106–126 with different sizes (2-peptides, 4-peptides and 6-peptides) on the surface of both pure neutral POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and negatively charged POPC/POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol) (3:1) lipids. MD simulation results show that PrP106–126 display strong affinity with POPC/POPG but does not interact with pure POPC. The positively charged and polar residues participating hydrogen bonding with membrane promote the adsorption of PrP106–126. The presence of POPC and POPC/POPG exert limited influence on the secondary structures of PrP106–126 and random coil structures are predominant in all simulation systems. Upon the adsorption on the POPC/POPG surface, the aggregation states of PrP106–126 have been changed and more small oligomers were observed. This work provides insights into the interactions of PrP106–126 and membranes with different compositions in atomic level, which expand our understanding the role membrane plays in the development of prion diseases. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy.     

Computational studies on prion proteins: effect of Ala(117)-->Val mutation

Molecular dynamics calculations demonstrated the conformational change in the prion protein due to Ala(117)-->Val mutation, which is related to Gerstmann-Str?ussler-Sheinker disease, one of the familial prion diseases. Three kinds of model structures of human and mouse prion proteins were examined: (model 1) nuclear magnetic resonance structures of human prion protein HuPrP (125-228) and mouse prion protein MoPrP (124-224), each having a globular domain consisting of three alpha-helices and an antiparallel beta-sheet; (model 2) extra peptides including Ala(117) (109-124 in HuPrP and 109-123 in MoPrP) plus the nuclear magnetic resonance structures of model 1; and (model 3) extra peptides including Val(117) (109-124 in HuPrP and 109-123 in MoPrP) plus the nuclear magnetic resonance structures of model 1. The results of molecular dynamics calculations indicated that the globular domains of models 1 and 2 were stable and that the extra peptide in model 2 tended to form a new alpha-helix. On the other hand, the globular domain of model 3 was unstable, and the beta-sheet region increased especially in HuPrP.     

Conformation of the cecropin-melittin hybrid, cecropin A(1-8)-melittin (1-18) antibacterial Peptide

Cecropin A (1-8)-Melittin (1-18) is a synthetic cecropin A-melittin hybrid peptide with leishmanicidal activity. The primary sequence of the peptide is as follows: KWKLPKKIGIGAVLKVLTTGLPALIS-NH2. 1H and 13C 2D NMR techniques were used to deduce the conformational parameters of chemical shift, 3JNHalpha coupling constants, temperature coefficients of NH chemical shifts and the pattern of intra and inter-residue nOe's. NMR studies were carried out in water (pH 6.0) and hexafluoroacetone (HFA). The peptide was found in a beta-pleated structure in water, and in HFA it adopts a right-handed alpha-helix conformation. Solution structures generated using restrained molecular dynamics simulations were refined by Mardigras to R factors ranging from 0.5 to 0.6.     

Concealment of epitope by reduction and alkylation in prion protein

Conversion of the cellular prion protein (PrP(C)) into its pathological isoform (PrP(Sc)), the key molecular event in the pathogenesis of prion diseases, is accompanied by a conformational transition of alpha-helix into beta-sheet structures involving alpha-helix 1 (alpha1) domain from residues 144 to 154 of the protein. Reduction and alkylation of PrP(C) have been found to inhibit the conversion of PrP(C) into PrP(Sc) in vitro. Here we report that while antibody affinity of epitopes in the N- and C-terminal domains remained unchanged, reduction and alkylation of the PrP molecule induced complete concealment of an epitope in alpha1 for anti-PrP antibody 6H4 that is able to cure prion infection in the cell model. Mass spectrometric analysis of recombinant PrP showed that the alkylation reaction takes place at reduced cysteines but no modification was observed in this cryptic epitope. Our study suggests that reduction and alkylation result in local or global rearrangement of PrP tertiary structure that is maintained in both liquid and solid phases. The implications in the conversion of PrP(C) into PrP(Sc) and the therapeutics of prion diseases are discussed.     

Optimization of hydrophobic domains in peptides that undergo transformation from alpha-helix to beta-fibril

Recent studies on peptide fibrillogenesis by the de novo method as well as amyloidogenic proteins including prion proteins and Alzheimer's beta-peptides have provided insights into the conformational changes, such as alpha-helix to beta-structure, involved in folding and misfolding processes. We have found that an exposed hydrophobic nucleation domain at N-terminal causes a structural transition of a peptide from alpha-helix to beta-fibril. It became clear that N-terminal acyl groups of particular lengths in a 2alpha-helix peptide caused the peptide to undergo an alpha-to-beta transition. The peptide with the octanoyl group (C8-2alpha) showed the highest rate of transformation. The study of the designed peptides revealed that these alpha-to-beta transitions were closely related to the initial alpha-helix conformation and its stability. Engineering peptides that undergo alpha-to-beta transitions are attractive not only to the study of pathogenic proteins such as prion proteins, but also to the control of self-assembly of peptides, which will lead to the development of peptidyl self-assembling materials.     

Possible role of region 152-156 in the structural duality of a peptide fragment from sheep prion protein

The conformational conversion of the nonpathogenic "cellular" prion isoform into a pathogenic "scrapie" protease-resistant isoform is a fundamental event in the onset of transmissible spongiform encephalopathies (TSE). During this pathogenic conversion, helix H1 and its two flanking loops of the normal prion protein are thought to undergo a conformational transition into a beta-like structure. A peptide spanning helix H1 and beta-strand S2 (residues 142-166 in human numbering) was studied by circular dichroism and nuclear magnetic resonance spectroscopies. This peptide in aqueous solution, in contrast to many prion fragments studied earlier (1) is highly soluble and (2) does not aggregate until the millimolar concentration range, and (3) exhibits an intrinsic propensity to a beta-hairpin-like conformation at neutral pH. We found that this peptide can also fold into a helix H1 conformation when dissolved in a TFE/PB mixture. The structures of the peptide calculated by MD showed solvent-dependent internal stabilizing forces of the structures and evidenced a higher mobility of the residues following the end of helix H1. These data suggest that the molecular rearrangement of this peptide in region 152-156, particularly in position 155, could be associated with the pathogenic conversion of the prion protein.     

Beta-hairpin conformation of fibrillogenic peptides: structure and alpha-beta transition mechanism revealed by molecular dynamics simulations

Understanding the conformational transitions that trigger the aggregation and amyloidogenesis of otherwise soluble peptides at atomic resolution is of fundamental relevance for the design of effective therapeutic agents against amyloid-related disorders. In the present study the transition from ideal alpha-helical to beta-hairpin conformations is revealed by long timescale molecular dynamics simulations in explicit water solvent, for two well-known amyloidogenic peptides: the H1 peptide from prion protein and the Abeta(12-28) fragment from the Abeta(1-42) peptide responsible for Alzheimer's disease. The simulations highlight the unfolding of alpha-helices, followed by the formation of bent conformations and a final convergence to ordered in register beta-hairpin conformations. The beta-hairpins observed, despite different sequences, exhibit a common dynamic behavior and the presence of a peculiar pattern of the hydrophobic side-chains, in particular in the region of the turns. These observations hint at a possible common aggregation mechanism for the onset of different amyloid diseases and a common mechanism in the transition to the beta-hairpin structures. Furthermore the simulations presented herein evidence the stabilization of the alpha-helical conformations induced by the presence of an organic fluorinated cosolvent. The results of MD simulation in 2,2,2-trifluoroethanol (TFE)/water mixture provide further evidence that the peptide coating effect of TFE molecules is responsible for the stabilization of the soluble helical conformation.     

Energy landscape of a peptide consisting of alpha-helix, 3(10)-helix, beta-turn, beta-hairpin, and other disordered conformations

The energy landscape of a peptide [Ace-Lys-Gln-Cys-Arg-Glu-Arg-Ala-Nme] in explicit water was studied with a multicanonical molecular dynamics simulation, and the AMBER parm96 force field was used for the energy calculation. The peptide was taken from the recognition helix of the DNA-binding protein, c-MYB: A rugged energy landscape was obtained, in which the random-coil conformations were dominant at room temperature. The CD spectra of the synthesized peptide revealed that it is in the random state at room temperature. However, the 300 K canonical ensemble, Q(300K), contained alpha-helix, 3(10)-helix, beta-turn, and beta-hairpin structures with small but notable probabilities of existence. The complete alpha-helix, imperfect alpha-helix, and random-coil conformations were separated from one another in the conformational space. This means that the peptide must overcome energy barriers to form the alpha-helix. The overcoming process may correspond to the hydrogen-bond rearrangements from peptide-water to peptide-peptide interactions. The beta-turn, imperfect 3(10)-helix, and beta-hairpin structures, among which there are no energy barriers at 300 K, were embedded in the ensemble of the random-coil conformations. Two types of beta-hairpin with different beta-turn regions were observed in Q(300K). The two beta-hairpin structures may have different mechanisms for the beta-hairpin formation. The current study proposes a scheme that the random state of this peptide consists of both ordered and disordered conformations. In contrast, the energy landscape obtained from the parm94 force field was funnel like, in which the peptide formed the helical conformation at room temperature and random coil at high temperature.     

Checking the pH-induced conformational transition of prion protein by molecular dynamics simulations: effect of protonation of histidine residues

The role of acidic pH in the conversion of human prion protein to the pathogenic isoform is investigated by means of molecular dynamics simulations, focusing the attention on the effect of protonation of histidine residues on the conformational behavior of human PrPC globular domain. Our simulations reveal a significant loss of alpha-helix content under mildly acidic conditions, due to destructuration of the C-terminal part of HB (thus suggesting a possible involvement of HB into the conformational transition leading to the pathogenic isoform) and a transient lengthening of the native beta-sheet. Protonation of His-187 and His-155 seems to be crucial for the onset of the conformational rearrangement. This finding can be related to the existence of a pathogenic mutation, H187R, which is associated with GSS syndrome. Finally, the relevance of our results for the location of a Cu2+-binding pocket in the C-terminal part of the prion is discussed.     

Characterization by Raman microspectroscopy of the strain-induced conformational transition in fibroin fibers from the silkworm Samia cynthia ricini

Raman microspectroscopy has been used to quantitatively study the effect of a mechanical deformation on the conformation and orientation of Samia cynthia ricini (S. c. ricini) silk fibroin. Samples were obtained from the aqueous solution stored in the silk gland and stretched at draw ratios (lambda) ranging from 0 to 11. Using an appropriate band decomposition procedure, polarized and orientation-insensitive spectra have been analyzed to determine order parameters and the content of secondary structures, respectively. The data unambiguously show that, in response to mechanical deformation, S. c. ricini fibroin undergoes a cooperative alpha-helix to beta-sheet conformational transition above a critical draw ratio of 4. The alpha-helix content decreases from 33 to 13% when lambda increases from 0 to 11, while the amount of beta-sheets increases from 15 to 37%. In comparison, cocoon silk is devoid of alpha-helical structure and always contains a larger amount of beta-sheets. Although the presence of isosbestic points in different spectral regions reveals that the conformational change induced by mechanical deformation is a two-state process, our results suggest that part of the glycine residues might be incorporated into beta-poly(alanine) structures. The beta-sheets are initially isotropically distributed and orient along the fiber axis as lambda increases, but do not reach the high level of orientation found in the cocoon fiber. The increase in the orientation level of the beta-sheets is found to be concomitant with the alpha --> beta conformational conversion, whereas alpha-helices do not orient under the applied strain but are rather readily converted into beta-sheets. The components assigned to turns exhibit a small orientation perpendicular to the fiber axis in stretched samples, showing that, overall, the polypeptide chains are aligned along the stretching direction. Our results suggest that, in nature, factors other than stretching contribute to the optimization of the amount of beta-sheets and the high degree of orientation found in natural cocoon silk.     

Dimensionality of Carbon Nanomaterials Determines the Binding and Dynamics of Amyloidogenic Peptides: Multiscale Theoretical Simulations

Experimental studies have demonstrated that nanoparticles can affect the rate of protein self-assembly, possibly interfering with the development of protein misfolding diseases such as Alzheimer''s, Parkinson''s and prion disease caused by aggregation and fibril formation of amyloid-prone proteins. We employ classical molecular dynamics simulations and large-scale density functional theory calculations to investigate the effects of nanomaterials on the structure, dynamics and binding of an amyloidogenic peptide apoC-II(60-70). We show that the binding affinity of this peptide to carbonaceous nanomaterials such as C60, nanotubes and graphene decreases with increasing nanoparticle curvature. Strong binding is facilitated by the large contact area available for π-stacking between the aromatic residues of the peptide and the extended surfaces of graphene and the nanotube. The highly curved fullerene surface exhibits reduced efficiency for π-stacking but promotes increased peptide dynamics. We postulate that the increase in conformational dynamics of the amyloid peptide can be unfavorable for the formation of fibril competent structures. In contrast, extended fibril forming peptide conformations are promoted by the nanotube and graphene surfaces which can provide a template for fibril-growth.