The conformations of the amyloid-beta (21-30) fragment can be described by three families in solution

Modeling amyloid fibril formation: a free-energy approach

Amyloid fibrils are structures consisting of many proteins with a well-defined conformation. The formation of these fibrils has been the subject of intense research, largely due to their connection to several diseases. We focus here on the computational studies and discuss these from a free-energy point of view. The fibrillogenic properties of many proteins can be predicted and understood by taking the relevant free energies into account in an appropriate way. This is because both the equilibrium and the kinetic properties of the protein system depend on its free-energy landscape. Advanced simulation techniques can be used to understand the relationship between the free-energy landscape of a protein and its three-dimensional structure and propensity to form amyloid fibrils. We give an overview of existing simulation techniques that operate at a molecular level of detail and that are capable of generating relevant free-energy values. The free energies obtained with these methods can be inserted into a statistical-mechanical or kinetic framework to predict mean fibril properties on length scales and time scales that are inaccessible by molecular-scale simulation methods.

Familial Alzheimer's disease mutations alter the stability of the amyloid beta-protein monomer folding nucleus

Amyloid beta-protein (Abeta) oligomers may be the proximate neurotoxins in Alzheimer's disease (AD). Recently, to elucidate the oligomerization pathway, we studied Abeta monomer folding and identified a decapeptide segment of Abeta, (21)Ala-(22)Glu-(23)Asp-(24)Val-(25)Gly-(26)Ser-(27)Asn-(28)Lys-(29)Gly-(30)Ala, within which turn formation appears to nucleate monomer folding. The turn is stabilized by hydrophobic interactions between Val-24 and Lys-28 and by long-range electrostatic interactions between Lys-28 and either Glu-22 or Asp-23. We hypothesized that turn destabilization might explain the effects of amino acid substitutions at Glu-22 and Asp-23 that cause familial forms of AD and cerebral amyloid angiopathy. To test this hypothesis, limited proteolysis, mass spectrometry, and solution-state NMR spectroscopy were used here to determine and compare the structure and stability of the Abeta(21-30) turn within wild-type Abeta and seven clinically relevant homologues. In addition, we determined the relative differences in folding free energies (DeltaDeltaG(f)) among the mutant peptides. We observed that all of the disease-associated amino acid substitutions at Glu-22 or Asp-23 destabilized the turn and that the magnitude of the destabilization correlated with oligomerization propensity. The Ala21Gly (Flemish) substitution, outside the turn proper (Glu-22-Lys-28), displayed a stability similar to that of the wild-type peptide. The implications of these findings for understanding Abeta monomer folding and disease causation are discussed.

Probing amyloid fibril formation of the NFGAIL peptide by computer simulations

Amyloid fibril formation, as observed in Alzheimer's disease and type II diabetes, is currently described by a nucleation-condensation mechanism, but the details of the process preceding the formation of the nucleus are still lacking. In this study, using an activation-relaxation technique coupled to a generic energy model, we explore the aggregation pathways of 12 chains of the hexapeptide NFGAIL. The simulations show, starting from a preformed parallel dimer and ten disordered chains, that the peptides form essentially amorphous oligomers or more rarely ordered beta-sheet structures where the peptides adopt a parallel orientation within the sheets. Comparison between the simulations indicates that a dimer is not a sufficient seed for avoiding amorphous aggregates and that there is a critical threshold in the number of connections between the chains above which exploration of amorphous aggregates is preferred.

Coarse-grained protein molecular dynamics simulations

A limiting factor in biological science is the time-scale gap between experimental and computational trajectories. At this point, all-atom explicit solvent molecular dynamics (MD) are clearly too expensive to explore long-range protein motions and extract accurate thermodynamics of proteins in isolated or multimeric forms. To reach the appropriate time scale, we must then resort to coarse graining. Here we couple the coarse-grained OPEP model, which has already been used with activated methods, to MD simulations. Two test cases are studied: the stability of three proteins around their experimental structures and the aggregation mechanisms of the Alzheimer's Abeta16-22 peptides. We find that coarse-grained isolated proteins are stable at room temperature within 50 ns time scale. Based on two 220 ns trajectories starting from disordered chains, we find that four Abeta16-22 peptides can form a three-stranded beta sheet. We also demonstrate that the reptation move of one chain over the others, first observed using the activation-relaxation technique, is a kinetically important mechanism during aggregation. These results show that MD-OPEP is a particularly appropriate tool to study qualitatively the dynamics of long biological processes and the thermodynamics of molecular assemblies.