Oligomerization of amyloid Abeta16-22 peptides using hydrogen bonds and hydrophobicity forces

Structures of soluble amyloid oligomers from computer simulations

Alzheimer's, Parkinson's, and Creutzfeldt-Jakob's neurodegenerative diseases are all linked with the assembly of normally soluble proteins into amyloid fibrils. Because of experimental limitations, structural characterization of the soluble oligomers, which form early in the process of fibrillogenesis and are cytotoxic, remains to be determined. In this article, we study the aggregation paths of seven chains of the shortest amyloid-forming peptide, using an activitated method and a reduced atomic representation. Our simulations show that disordered KFFE monomers ultimately form three distinct topologies of similar energy: amorphous oligomers, incomplete rings with beta-barrel character, and cross-beta-sheet structures with the meridional but not the equatorial X-ray fiber reflections. The simulations also shed light on the pathways from misfolded aggregates to fibrillar-like structures. They also underline the multiplicity of building blocks that can lead to the formation of the critical nucleus from which rapid growth of the fibril occurs.

Thermal versus mechanical unfolding of ubiquitin

The authors studied the temperature-induced unfolding of ubiquitin by all-atom Monte Carlo simulations. The unfolding behavior is compared with that seen in previous simulations of the mechanical unfolding of this protein, based on the same model. In mechanical unfolding, secondary-structure elements were found to break in a quite well-defined order. In thermal unfolding, the authors saw somewhat larger event-to-event fluctuations, but the unfolding pathway was still far from random. Two long-lived secondary-structure elements could be identified in the simulations. These two elements have been found experimentally to be the thermally most stable ones. Interestingly, one of these long-lived elements, the first beta-hairpin, was found to break early in the mechanical unfolding simulations. Their combined simulation results thus enable the authors to predict in detail important differences between the thermal and mechanical unfolding behaviors of ubiquitin.

Folding Landscapes of the Alzheimer Amyloid-β(12-28) Peptide

The energy landscape for folding of the 12-28 fragment of the Alzheimer amyloid beta (Abeta) peptide is characterized using replica-exchange molecular dynamics simulations with an all-atom peptide model and explicit solvent. At physiological temperatures, the peptide exists mostly as a collapsed random coil, populating a small fraction (less than 10%) of hairpins with a beta-turn at position V18F19, with another 10% of hairpin-like conformations possessing a bend rather than a turn in the central VFFA positions. A small fraction of the populated states, approximately 14%, adopt polyproline II (PPII) conformations. Folding of the structured hairpin states proceeds through the assembly of two locally stable segments, VFFAE and EDVGS. The interactions stabilizing these locally folded structural motifs are in conflict with those stabilizing the global fold of A12-28, a signature of underlying residual frustration in this peptide. At increased temperature, the population of both beta-strand and PPII conformations diminishes in favor of beta-turn and random-coil states. On the basis of the conformational preferences of Abeta 12-28 monomers, two models for the molecular structure of amyloid fibrils formed by this peptide are proposed.

Stability and structure of oligomers of the Alzheimer peptide Abeta16-22: from the dimer to the 32-mer

Several neurodegenerative diseases such as Alzheimer's, Parkinson's, and Huntington's diseases are associated with amyloid fibrils formed by different polypeptides. We probe the structure and stability of oligomers of different sizes of the fragment Abeta(16-22) of the Alzheimer beta-amyloid peptide using atomic-detail molecular dynamics simulations with explicit solvent. We find that only large oligomers form a stable beta-sheet aggregate, the minimum nucleus size being of the order of 8-16 peptides. This effect is attributed to better hydrophobic contacts and a better shielding of backbone-backbone hydrogen bonds from the solvent in bigger assemblies. Moreover, the observed stability of beta-sheet aggregates with a different number of layers can be explained on the basis of their solvent-accessible surface area. Depending on the stacking interface between the sheets, we observe straight or twisted structures, which could be linked to the experimentally observed polymorphism of amyloid fibrils. To compare our 32-mer structure to experimental data, we calculate its x-ray diffraction pattern. Good agreement is found between experimentally and theoretically determined reflections, suggesting that our model indeed closely resembles the structures found in vitro.

Effects of solvent on the structure of the Alzheimer amyloid-beta(25-35) peptide

The free energy landscape for folding of the Alzheimer's amyloid-beta(25-35) peptide is explored using replica exchange molecular dynamics in both pure water and in HFIP/water cosolvent. This amphiphilic peptide is a natural by-product of the Alzheimer's amyloid-beta(1-40) peptide and retains the toxicity of its full-length counterpart as well as the ability to aggregate into beta-sheet-rich fibrils. Our simulations reveal that the peptide preferentially populates a helical structure in apolar organic solvent, while in pure water, the peptide adopts collapsed coil conformations and to a lesser extent beta-hairpin conformations. The beta-hairpin is characterized by a type II' beta-turn involving residues G29 and A30 and two short beta-strands involving residues N27, K28, I31, and I32. The hairpin is stabilized by backbone hydrogen-bonding interactions between residues K28 and I31; S26 and G33; and by side-chain-to-side-chain interactions between N27 and I32. Implications regarding the mechanism of aggregation of this peptide into fibrils and the role of the environment in modulating secondary structure are discussed.

Atomic-Level Description of Amyloid β-Dimer Formation

All-atom simulations have been carried out on a monomer and dimer of the aggregation-prone fragment (16-22) of amyloid beta peptide, which is implicated in Alzheimer's disease. The replica exchange molecular dynamics method, which has been successfully applied to peptide folding, is utilized as a means to sample the configurational space with proper Boltzmann weighting so that the structural, motional, and thermodynamic description of self-assembly can be obtained. The free energy landscape showing the delicate balance between different monomer and dimer conformations is mapped along carefully chosen reaction coordinates. The canonical ensembles at 38 different temperatures are used to describe the thermodynamics and the relative stabilities of at least six different dimer conformations including that of parallel and antiparallel orientations. We also delineate the nature of the molecular forces that activate and stabilize these different dimer conformations as a function of temperature, especially as related to secondary structural propensity of monomer. We identify parallel loop dimer conformations that are stabilized due to specific interactions with water molecules.

PROFASI: A Monte Carlo simulation package for protein folding and aggregation

We present a flexible and efficient program package written in C++, PROFASI, for simulating protein folding and aggregation. The systems are modeled using an all-atom description of the protein chains with only torsional degrees of freedom, and implicit water. The program package has a modular structure that makes the interaction potential easy to modify. The currently implemented potential is able to fold several peptides with about 20 residues, and has also been used to study aggregation and force-induced unfolding. The simulation methods implemented in PROFASI are Monte Carlo-based and include a semilocal move and simulated tempering. Adding new updates is easy. The code runs fast in both single- and multi-chain applications, as is illustrated by several examples.

Exploring the early steps of amyloid peptide aggregation by computers

The assembly of normally soluble proteins into amyloid fibrils is a hallmark of neurodegenerative diseases. Because protein aggregation is very complex, involving a variety of oligomeric metastable intermediates, the detailed aggregation paths and structural characterization of the intermediates remain to be determined. Yet, there is strong evidence that these oligomers, which form early in the process of fibrillogenesis, are cytotoxic. In this paper, we review our current understanding of the underlying factors that promote the aggregation of peptides into amyloid fibrils. We focus here on the structural and dynamic aspects of the aggregation as observed in state-of-the-art computer simulations of amyloid-forming peptides with an emphasis on the activation-relaxation technique.

Elongation of Ordered Peptide Aggregate of an Amyloidogenic Hexapeptide NFGAIL Observed in Molecular Dynamics Simulations with Explicit Solvent

The mechanisms by which amyloidogenic peptides and proteins form soluble toxic oligomers remain elusive. We have studied the formation of partially ordered tetramers and well-ordered octamers of an amyloidogenic hexapeptide NFGAIL (residues 22-27 of the human islet amyloid polypeptide) in our previous work. Continuing the effort, we here probe the beta-sheet elongation process by a combined total of 2.0 micros molecular dynamics simulations with explicit solvent. In a set of 10 simulations with the peptides restrained to the extended conformation, we observed that the main growth mode was elongation along the beta-sheet hydrogen bonds through primarily a two-stage process. Driven by hydrophobic forces, the peptides initially attached to the surface of the ordered oligomer, moved quickly to the beta-sheet edges, and formed stable beta-sheet hydrogen bonds. Addition of peptides to the existing oligomer notably improved the order of the peptide aggregate in which labile outer layer beta-sheets were stabilized, which provides good templates for further elongation. These simulations suggested that elongation along the beta-sheet hydrogen bonds occurs at the intermediate stage when low-weight oligomers start to form. We did not observe significant preference toward either parallel or antiparallel beta-sheets at the elongation stage for this peptide. In another set of 10 unrestrained simulations, the dominant growth mode was disordered aggregation. Taken together, these results offered a glimpse at the molecular events leading to the formation of ordered and disordered low-weight oligomers.

Dissecting the mechanical unfolding of ubiquitin

The unfolding behavior of ubiquitin under the influence of a stretching force recently was investigated experimentally by single-molecule constant-force methods. Many observed unfolding traces had a simple two-state character, whereas others showed clear evidence of intermediate states. Here, we use Monte Carlo simulations to investigate the force-induced unfolding of ubiquitin at the atomic level. In agreement with experimental data, we find that the unfolding process can occur either in a single step or through intermediate states. In addition to this randomness, we find that many quantities, such as the frequency of occurrence of intermediates, show a clear systematic dependence on the strength of the applied force. Despite this diversity, one common feature can be identified in the simulated unfolding events, which is the order in which the secondary-structure elements break. This order is the same in two- and three-state events and at the different forces studied. The observed order remains to be verified experimentally but appears physically reasonable.

Lattice models of peptide aggregation: Evaluation of conformational search algorithms

We present a series of conformational search calculations on the aggregation of short peptide fragments that form fibrils similar to those seen in many protein mis-folding diseases. The proteins were represented by a face-centered cubic lattice model with the conformational energies calculated using the Miyazawa-Jernigan potential. The searches were performed using algorithms based on the Metropolis Monte Carlo method, including simulated annealing and replica exchange. We also present the results of searches using the tabu search method, an algorithm that has been used for many optimization problems, but has rarely been used in protein conformational searches. The replica exchange algorithm consistently found more stable structures then the other algorithms, and was particularly effective for the octamers and larger systems.

Folding thermodynamics of peptides

A simplified interaction potential for protein folding studies at the atomic level is discussed and tested on a set of peptides with approximately 20 residues each. The test set contains both alpha-helical (Trp cage, F(s)) and beta-sheet (GB1p, GB1m2, GB1m3, Betanova, LLM) peptides. The model, which is entirely sequence-based, is able to fold these different peptides for one and the same choice of model parameters. Furthermore, the melting behavior of the peptides is in good quantitative agreement with experimental data. Apparent folded populations obtained using different observables are compared, and are found to be very different for some of the peptides (e.g., Betanova). In other cases (in particular, GB1m2 and GB1m3), the different estimates agree reasonably well, indicating a more two-state-like melting behavior.