Self-Assembly of the β2-Microglobulin NHVTLSQ Peptide Using a Coarse-Grained Protein Model Reveals a β-Barrel Species

A multiscale approach to characterize the early aggregation steps of the amyloid-forming peptide GNNQQNY from the yeast prion sup-35

The self-organization of peptides into amyloidogenic oligomers is one of the key events for a wide range of molecular and degenerative diseases. Atomic-resolution characterization of the mechanisms responsible for the aggregation process and the resulting structures is thus a necessary step to improve our understanding of the determinants of these pathologies. To address this issue, we combine the accelerated sampling properties of replica exchange molecular dynamics simulations based on the OPEP coarse-grained potential with the atomic resolution description of interactions provided by all-atom MD simulations, and investigate the oligomerization process of the GNNQQNY for three system sizes: 3-mers, 12-mers and 20-mers. Results for our integrated simulations show a rich variety of structural arrangements for aggregates of all sizes. Elongated fibril-like structures can form transiently in the 20-mer case, but they are not stable and easily interconvert in more globular and disordered forms. Our extensive characterization of the intermediate structures and their physico-chemical determinants points to a high degree of polymorphism for the GNNQQNY sequence that can be reflected at the macroscopic scale. Detailed mechanisms and structures that underlie amyloid aggregation are also provided.

The formation of amyloid fibrils is associated with many neurodegenerative diseases such as Alzheimer's, Creutzfeld-Jakob, Parkinson's, the Prion disease and diabetes mellitus. In all cases, proteins misfold to form highly ordered insoluble aggregates called amyloid fibrils that deposit intra- and extracellularly and are resistant to proteases. All proteins are believed to have the instrinsic capability of forming amyloid fibrils that share common specific structural properties that have been observed by X-ray crystallography and by NMR. However, little is known about the aggregation dynamics of amyloid assemblies, and their toxicity mechanism is therefore poorly understood. It is believed that small amyloid oligomers, formed on the aggregation pathway of full amyloid fibrils, are the toxic species. A detailed atomic characterization of the oligomerization process is thus necessary to further our understanding of the amyloid oligomer's toxicity. Our approach here is to study the aggregation dynamics of a 7-residue amyloid peptide GNNQQNY through a combination of numerical techniques. Our results suggest that this amyloid sequence can form fibril-like structures and is polymorphic, which agrees with recent experimental observations. The ability to fully characterize and describe the aggregation pathway of amyloid sequences numerically is key to the development of future drugs to target amyloid oligomers.

HiRE-RNA: a high resolution coarse-grained energy model for RNA

Although RNAs play many cellular functions, the gap between their sequences and 3D structures is increasing and our knowledge of RNA thermodynamics and long time scale dynamics is still limited at an atomic level of detail. In principle, all-atom molecular dynamics (MD) and replica exchange molecular dynamics (REMD) simulations can investigate these issues, but with current computer facilities, these simulations in explicit solvent have been limited to small RNAs and to short times. To move to larger systems, we can resort to coarse-graining. In this study, we present HiRE-RNA, a generic high resolution coarse-grained model for RNA, and report MD and REMD simulations on two RNAs of 22 and 36 nucleotides. Starting from unfolded structures, the 22-mer folds within 1.8 A rmsd from the NMR structure, while the 36-mer folds within 4.6 A rmsd. Current results suggest that further optimization of the HiRE-RNA force field should open the door to a relevant model for studying large RNAs, such as riboswitches, and for predicting 3D structures from secondary structure information.

Effects of G33A and G33I mutations on the structures of monomer and dimer of the amyloid-β fragment 29-42 by replica exchange molecular dynamics simulations

The early formed oligomers of amyloid-β proteins with 40 and 42 amino acids are believed to be the culprits of Alzheimer's disease. Aβ1-42 peptides with alanine and isoleucine mutations of glycine 33 are known to be much less toxic than the wild-type Aβ1-42 and promote the aggregation process in vitro. The fragment Aβ29-42 has also been shown to form fibrils, disrupt Aβ1-42 oligomerization, and inhibit Aβ1-42-induced neurotoxicity. As a first step toward understanding the impact of G33A and G33I mutations on the earliest steps along the Aβ1-42 aggregation pathway, we have studied the structures of the monomer and dimer of Aβ29-42 and its two G33 variants using coarse-grained replica exchange molecular dynamics simulations. These simulations, totaling 15 μs, indicate that both substitutions impact the conformational ensemble of Aβ29-42. For the monomer, the population of the β-hairpin is high for wild-type Aβ29-42, but marginal for Aβ29-42 G33I mutant. The three dimers are also stabilized by different patterns of interaction. The data are discussed in terms of the differences in the aggregation characteristics between wild-type Aβ1-42 and its two G33A and G33I variants.

Molecular dynamics simulation of β₂-microglobulin in denaturing and stabilizing conditions

β₂-Microglobulin has been a model system for the study of fibril formation for 20 years. The experimental study of β₂-microglobulin structure, dynamics, and thermodynamics in solution, at atomic detail, along the pathway leading to fibril formation is difficult because the onset of disorder and aggregation prevents signal resolution in Nuclear Magnetic Resonance experiments. Moreover, it is difficult to characterize conformers in exchange equilibrium. To gain insight (at atomic level) on processes for which experimental information is available at molecular or supramolecular level, molecular dynamics simulations have been widely used in the last decade. Here, we use molecular dynamics to address three key aspects of β₂-microglobulin, which are known to be relevant to amyloid formation: (1) 60 ns molecular dynamics simulations of β₂-microglobulin in trifluoroethanol and in conditions mimicking low pH are used to study the behavior of the protein in environmental conditions that are able to trigger amyloid formation; (2) adaptive biasing force molecular dynamics simulation is used to force cis-trans isomerization at Proline 32 and to calculate the relative free energy in the folded and unfolded state. The native-like trans-conformer (known as intermediate 2 and determining the slow phase of refolding), is simulated for 10 ns, detailing the possible link between cis-trans isomerization and conformational disorder; (3) molecular dynamics simulation of highly concentrated doxycycline (a molecule able to suppress fibril formation) in the presence of β₂-microglobulin provides details of the binding modes of the drug and a rationale for its effect.

A fast method for large-scale de novo peptide and miniprotein structure prediction

Although peptides have many biological and biomedical implications, an accurate method predicting their equilibrium structural ensembles from amino acid sequences and suitable for large-scale experiments is still missing. We introduce a new approach-PEP-FOLD-to the de novo prediction of peptides and miniproteins. It first predicts, in the terms of a Hidden Markov Model-derived structural alphabet, a limited number of local conformations at each position of the structure. It then performs their assembly using a greedy procedure driven by a coarse-grained energy score. On a benchmark of 52 peptides with 9-23 amino acids, PEP-FOLD generates lowest-energy conformations within 2.8 and 2.3 A Calpha root-mean-square deviation from the full nuclear magnetic resonance structures (NMR) and the NMR rigid cores, respectively, outperforming previous approaches. For 13 miniproteins with 27-49 amino acids, PEP-FOLD reaches an accuracy of 3.6 and 4.6 A Calpha root-mean-square deviation for the most-native and lowest-energy conformations, using the nonflexible regions identified by NMR. PEP-FOLD simulations are fast-a few minutes only-opening therefore, the door to in silico large-scale rational design of new bioactive peptides and miniproteins.

Structures and thermodynamics of Alzheimer's amyloid-beta Abeta(16-35) monomer and dimer by replica exchange molecular dynamics simulations: implication for full-length Abeta fibrillation

Many proteins display a strand-loop-strand motif in their amyloid fibrillar states. For instance, the amyloid beta-protein, Abeta1-40, associated with Alzheimer's disease, displays a loop at positions 22-28 in its amyloid fibril state. It has been suggested that this loop could appear early in the aggregation process, but quantitative information regarding its presence in small oligomers remains scant. Because residues 1-15 are disordered in Abeta1-42 fibrils and Abeta10-35 forms fibrils in vitro, we select the peptide Abeta16-35, centered on residues 22-28 and determine the structures and thermodynamics of the monomer and dimer using coarse-grained implicit solvent replica exchange molecular dynamics simulations. Our simulations totalling 5 mus for the monomer and 12 micros for the dimer show no sign of strong secondary structure signals in both instances and the significant impact of dimerization on the global structure of Abeta16-35. They reveal however that the loop 22-28 acts as a quasi-independent unit in both species. The loop structure ensemble we report in Abeta16-35 monomer and dimer has high similarity to the loop formed by the Abeta21-30 peptide in solution and, to a lesser extent, to the loop found in Abeta1-40 fibrils. We discuss the implications of our findings on the assembly of full-length Abeta.

Self-assembly of phenylalanine oligopeptides: insights from experiments and simulations

Studies of peptide-based nanostructures provide general insights into biomolecular self-assembly and can lead material engineering toward technological applications. The diphenylalanine peptide (FF) self-assembles into discrete, hollow, well ordered nanotubes, and its derivatives form nanoassemblies of various morphologies. Here we demonstrate for the first time, to our knowledge, the formation of planar nanostructures with beta-sheet content by the triphenylalanine peptide (FFF). We characterize these structures using various microscopy and spectroscopy techniques. We also obtain insights into the interactions and structural properties of the FF and FFF nanostructures by 0.4-micros, implicit-solvent, replica-exchange, molecular-dynamics simulations of aqueous FF and FFF solutions. In the simulations the peptides form aggregates, which often contain open or ring-like peptide networks, as well as elementary and network-containing structures with beta-sheet characteristics. The networks are stabilized by polar and nonpolar interactions, and by the surrounding aggregate. In particular, the charged termini of neighbor peptides are involved in hydrogen-bonding interactions and their aromatic side chains form "T-shaped" contacts, as in three-dimensional FF crystals. These interactions may assist the FF and FFF self-assembly at the early stage, and may also stabilize the mature nanostructures. The FFF peptides have higher network propensities and increased aggregate stabilities with respect to FF, which can be interpreted energetically.

A model of amyloid's role in disease based on fibril fracture

Although the correlative evidence relating the presence of amyloid fibrils and certain disease states (e.g. Alzheimer's disease and Type 2 Diabetes) is overwhelming, a direct causative role for amyloid has proved harder to establish. Current thinking links a narrow region of the aggregate protein size distribution, the so called 'early aggregate' domain to cellular toxicity. A troubling feature of this theory however is that the nucleated reaction mechanism by which amyloid formation is believed to occur results in a very low number concentration of early aggregates which are rapidly extended to form amyloid fibrils. This situation means that the concentration of early aggregates is very low at the time when they are supposedly at their most toxic. Here we adopt a novel explicit simulation strategy to examine a kinetic regime involving nucleated growth combined with fibril fragmentation under which this situation can be reversed so as to produce a high number concentration of small on-pathway toxic aggregates. Dependent upon the rate of fragmentation, the time scale for generation of toxic early aggregates may be coupled, uncoupled or disassociated from the time scale for the appearance of amyloid fibrils. Furthermore the model presents itself as a biochemical 'switch' transitioning between modes of amyloid induced cell death dependent upon either specific amyloid toxicity or non-specific solid mass induced tissue damage.

Structural diversity of the soluble trimers of the human amylin(20-29) peptide revealed by molecular dynamics simulations

The human islet amyloid polypeptide (hIAPP) or amylin is a 37-residue hormone found as amyloid deposits in pancreatic extracts of nearly all type 2 diabetes patients. The fragment 20-29 of sequence SNNFGAILSS (hIAPP20-29) has been shown to be responsible for the amyloidogenic propensities of the full length protein. Various polymorphic forms of hIAPP20-29 fibrils were described by using Fourier transform infrared (FTIR) and solid-state NMR experiments: unseeded hIAPP20-29 fibril with out-of-register antiparallel beta-strands, and two forms of seeded hIAPP20-29 fibril, with in-register antiparallel or in-register parallel beta-strands. As a first step toward understanding this polymorphism, we explore the equilibrium structures of the soluble hIAPP20-29 trimer, using multiple molecular dynamics (MD) simulations with the Optimized Potential for Efficient structure Prediction (OPEP) coarse-grained implicit solvent force field for a total length of 3.2 micros. Although, the trimer is found mainly random coil, consistent with the signal measured experimentally during the lag phase of hIAPP20-29 fibril formation, the central FGAIL residues have a relative high propensity to form interpeptide beta-sheets and antiparallel beta-strands are more probable than parallel beta-strands. One MD-predicted out-of-register antiparallel three-stranded beta-sheet matches exactly the FTIR-derived unseeded hIAPP20-29 fibril model. Our simulations, however, do not reveal any evidence of in-register parallel or in-register antiparallel beta-sheets as reported for seeded hIAPP20-29 fibrils. All these results indicate that fibril polymorphism is partially encoded in a trimer.

Targeting the early steps of Abeta16-22 protofibril disassembly by N-methylated inhibitors: a numerical study

Aggregation of the Abeta1-40/Abeta1-42 peptides is a key factor in Alzheimer's disease. Though the inhibitory effect of N-methylated Abeta16-22 (mAbeta16-22) peptides is well characterized in vitro, there is little information on how they disassemble full-length Abeta fibrils or block fibril formation. Here, we report coarse-grained implicit solvent molecular dynamics (MD) and replica exchange molecular dynamics (REMD) simulations on Abeta16-22 and mAbeta16-22 monomers, and then a preformed six-chain Abeta16-22 bilayer with either four copies of Abeta16-22 or four copies of mAbeta16-22. Our simulations show that the effect of N-methylation on mAbeta16-22 monomer is to reduce the density of compact forms. While 100 ns MD trajectories do not reveal any significant differences between the two ten-chain systems, the REMD simulations totaling 1 micros help understand the first steps of Abeta16-22 protofibril disassembly by N-methylated inhibitors. Notably, we find that mAbeta16-22 preferentially interacts with Abeta16-22 by blocking both beta-sheet extension and lateral association of layers, but also by intercalation of the inhibitors allowing sequestration of Abeta16-22 peptides. This third binding mode is particularly appealing for blocking Abeta fibrillogenesis.

Replica exchange molecular dynamics simulations of coarse-grained proteins in implicit solvent

Current approaches aimed at determining the free energy surface of all-atom medium-size proteins in explicit solvent are slow and are not sufficient to converge to equilibrium properties. To ensure a proper sampling of the configurational space, it is preferable to use reduced representations such as implicit solvent and/or coarse-grained protein models, which are much lighter computationally. Each model must be verified, however, to ensure that it can recover experimental structures and thermodynamics. Here we test the coarse-grained implicit solvent OPEP model with replica exchange molecular dynamics (REMD) on six peptides ranging in length from 10 to 28 residues: two alanine-based peptides, the second beta-hairpin from protein G, the Trp-cage and zinc-finger motif, and a dimer of a coiled coil peptide. We show that REMD-OPEP recovers the proper thermodynamics of the systems studied, with accurate structural description of the beta-hairpin and Trp-cage peptides (within 1-2 A from experiments). The light computational burden of REMD-OPEP, which enables us to generate many hundred nanoseconds at each temperature and fully assess convergence to equilibrium ensemble, opens the door to the determination of the free energy surface of larger proteins and assemblies.

Formation and growth of oligomers: a Monte Carlo study of an amyloid tau fragment

Small oligomers formed early in the process of amyloid fibril formation may be the major toxic species in Alzheimer's disease. We investigate the early stages of amyloid aggregation for the tau fragment AcPHF6 (Ac-VQIVYK-NH2) using an implicit solvent all-atom model and extensive Monte Carlo simulations of 12, 24, and 36 chains. A variety of small metastable aggregates form and dissolve until an aggregate of a critical size and conformation arises. However, the stable oligomers, which are β-sheet-rich and feature many hydrophobic contacts, are not always growth-ready. The simulations indicate instead that these supercritical oligomers spend a lengthy period in equilibrium in which considerable reorganization takes place accompanied by exchange of chains with the solution. Growth competence of the stable oligomers correlates with the alignment of the strands in the β-sheets. The larger aggregates seen in our simulations are all composed of two twisted β-sheets, packed against each other with hydrophobic side chains at the sheet–sheet interface. These β-sandwiches show similarities with the proposed steric zipper structure for PHF6 fibrils but have a mixed parallel/antiparallel β-strand organization as opposed to the parallel organization found in experiments on fibrils. Interestingly, we find that the fraction of parallel β-sheet structure increases with aggregate size. We speculate that the reorganization of the β-sheets into parallel ones is an important rate-limiting step in the formation of PHF6 fibrils.

It is believed that the self association of certain protein molecules into aggregated structures, known as amyloid fibrils, plays an important role in a variety of human diseases, such as Alzheimer's disease and Parkinson's disease. Although the ability to form such amyloid fibrils is a common property for proteins, the process leading to these fibrils is incompletely understood. The early stages of the process involve small transient heterogeneous structures made of a few protein chains and are especially difficult to characterize. Here we use atomic-level simulations to explore the early part of the aggregation process for a fibril-forming fragment of the protein tau associated with Alzheimer's disease. We find that a multitude of small aggregates, rich in sheetlike structures, form through a nucleation process. Interestingly, a statistically preferred type of aggregate, consisting of two tightly packed sheets, emerges with increasing aggregate size. Growth of these larger aggregates seems to be a slow process that correlates with the emergence of more uniformly ordered sheets. We speculate that reorganization of the protein chains leading to that ordered arrangement is an important bottleneck to amyloid fibril formation for this peptide.

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 beta-sheets separated by a dry side-chain interface were adopted in our computational analysis. A cumulative time of 1 micros 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 micros 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 beta-sheets form earlier than the dry side-chain contacts during aggregation. These results provide further insights in understanding the early fibrils aggregation.