Microscopic factors that control beta-sheet registry in amyloid fibrils formed by fragment 11-25 of amyloid beta peptide: insights from computer simulations

Molecular dynamics simulations of amyloid-β(16-22) peptide aggregation at air-water interfaces

Oligomers of amyloid-β (Aβ) peptides are known to be related to Alzheimer's disease, and their formation is accelerated at hydrophilic-hydrophobic interfaces, such as the cell membrane surface and air-water interface. Here, we report molecular dynamics simulations of aggregation of Aβ(16-22) peptides at air-water interfaces. First, 100 randomly distributed Aβ(16-22) peptides moved to the interface. The high concentration of peptides then accelerated their aggregation and formation of antiparallel β-sheets. Two layers of oligomers were observed near the interface. In the first layer from the interface, the oligomer with less β-bridges exposed the hydrophobic residues to the air. The second layer consisted of oligomers with more β-bridges that protruded into water. They are more soluble in water because the hydrophobic residues are covered by N- and C-terminal hydrophilic residues that are aligned well along the oligomer edge. These results indicate that amyloid protofibril formation mainly occurs in the second layer.

A theoretical study of polymorphism in VQIVYK fibrils

The VQIVYK fragment from the Tau protein, also known as PHF6, is essential for aggregation of Tau into neurofibrillary lesions associated with neurodegenerative diseases. VQIVYK itself forms amyloid fibrils composed of paired β-sheets. Therefore, the full Tau protein and VQIVYK fibrils have been intensively investigated. A central issue in these studies is polymorphism, the ability of a protein to fold into more than one structure. Using all-atom molecular simulations, we generate five stable polymorphs of VQIVYK fibrils, establish their relative free energy with umbrella sampling methods, and identify the side chain interactions that provide stability. The two most stable polymorphs, which have nearly equal free energy, are formed by interdigitation of the mostly hydrophobic VIY "face" sides of the β-sheets. Another stable polymorph is formed by interdigitation of the QVK "back" sides. When we turn to examine structures from cryo-electron microscopy experiments on Tau filaments taken from diseased patients or generated in vitro, we find that the pattern of side chain interactions found in the two most stable face-to-face as well as the back-to-back polymorphs are recapitulated in amyloid structures of the full protein. Thus, our studies suggest that the interactions stabilizing PHF6 fibrils explain the amyloidogenicity of the VQIVYK motif within the full Tau protein and provide justification for the use of VQIVYK fibrils as a test bed for the design of molecules that identify or inhibit amyloid structures.

pH-dependent conformational ensemble and polymorphism of amyloid-β core fragment

Characterization of amyloid oligomeric species is important due to its possible responsibility for the toxicity of amyloid proteins, whereas it is difficult to detect by current spectroscopic techniques. The pH-dependent tetramerization and fibrillation of the central hydrophobic segment of Alzheimer amyloid β-peptide (Aβ(12-24)) were respectively explored by all-atom replica exchange molecular dynamics simulations and by fluorescence and atomic force microscopy measurements. Our combined study shows that more β-sheet structures in the early event of tetramerization is linked directly to the high propensity to form amyloid fibrils in the consequent fibrillation. Both tetramerization and fibrillation are strongly regulated by pH. At pH 5.0, Aβ(12-24) has two opposite terminal charges. The electrostatic attraction between the side-chains of His13/His14 and Glu22/Asp23 thus acts as a "pattern keeper", resulting in high propensity of amyloid formation. These results suggest that pH effects most likely by affecting the ionization properties of the Aβ(12-24) peptide. Specifically, the pH-dependent equilibrium conformational distribution of different aggregate species are well-investigated in detail. Our findings also give hints to other experimental findings that the kinetics and morphologies of Aβ fibril formation are strongly pH-dependent.

Dynamical stability and assembly cooperativity of β-sheet amyloid oligomers--effect of polarization

The soluble intermediate oligomers of amyloidogenic proteins are suspected to be more cytotoxic than the mature fibrils in neurodegenerative disorders. Here, the dynamic stability and assembly cooperativity of a model oligomer of human islet amyloid polypeptide (hIAPP) segments were explored by means of all-atom molecular dynamics (MD) simulations under different force fields including AMBER99SB, OPLS, and polarized protein-specific charge (PPC) model. Simulation results show that the dynamic stability of β-sheet oligomers is seriously impacted by electrostatic polarization. Without inclusion of polarization (simulation under standard AMBER and OPLS force field), the β-sheet oligomers are dynamically unstable during MD simulation. For comparison, simulation results under PPC give significantly more stable dynamical structures of the oligomers. Furthermore, calculation of electrostatic interaction energy between the neighboring β strands with an approximate polarizable method produces energetic evidence for cooperative assembly of β-strand oligomers. This result supports a picture of downhill-like cooperative assembly of β strands during fibrillation process. The present study demonstrates the critical role of polarization in dynamic stability and assembly cooperativity of β-sheet-rich amyloid oligomers.

Nanomaterials design and tests for neural tissue engineering

Nanostructured scaffolds recently showed great promise in tissue engineering: nanomaterials can be tailored at the molecular level and scaffold morphology may more closely resemble features of extracellular matrix components in terms of porosity, framing and biofunctionalities. As a consequence, both biomechanical properties of scaffold microenvironments and biomaterial-protein interactions can be tuned, allowing for improved transplanted cell engraftment and better controlled diffusion of drugs. Easier said than done, a nanotech-based regenerative approach encompasses different fields of know-how, ranging from in silico simulations, nanomaterial synthesis and characterization at the nano-, micro- and mesoscales to random library screening methods (e.g. phage display), in vitro cellular-based experiments and validation in animal models of the target injury. All of these steps of the "assembly line" of nanostructured scaffolds are tightly interconnected both in their standard analysis techniques and in their most recent breakthroughs: indeed their efforts have to jointly provide the deepest possible analyses of the diverse facets of the challenging field of neural tissue engineering. The purpose of this review is therefore to provide a critical overview of the recent advances in and drawbacks and potential of each mentioned field, contributing to the realization of effective nanotech-based therapies for the regeneration of peripheral nerve transections, spinal cord injuries and brain traumatic injuries. Far from being the ultimate overview of such a number of topics, the reader will acknowledge the intrinsic complexity of the goal of nanotech tissue engineering for a conscious approach to the development of a regenerative therapy and, by deciphering the thread connecting all steps of the research, will gain the necessary view of its tremendous potential if each piece of stone is correctly placed to work synergically in this impressive mosaic.

Biochemistry of amyloid β-protein and amyloid deposits in Alzheimer disease

Progressive cerebral deposition of the amyloid β-protein (Aβ) in brain regions serving memory and cognition is an invariant and defining feature of Alzheimer disease. A highly similar but less robust process accompanies brain aging in many nondemented humans, lower primates, and some other mammals. The discovery of Aβ as the subunit of the amyloid fibrils in meningocerebral blood vessels and parenchymal plaques has led to innumerable studies of its biochemistry and potential cytotoxic properties. Here we will review the discovery of Aβ, numerous aspects of its complex biochemistry, and current attempts to understand how a range of Aβ assemblies, including soluble oligomers and insoluble fibrils, may precipitate and promote neuronal and glial alterations that underlie the development of dementia. Although the role of Aβ as a key molecular factor in the etiology of Alzheimer disease remains controversial, clinical trials of amyloid-lowering agents, reviewed elsewhere in this book, are poised to resolve the question of its pathogenic primacy.

Intrinsic determinants of Aβ(12-24) pH-dependent self-assembly revealed by combined computational and experimental studies

The propensity of amyloid- (A) peptide to self-assemble into highly ordered amyloid structures lies at the core of their accumulation in the brain during Alzheimer's disease. By using all-atom explicit solvent replica exchange molecular dynamics simulations, we elucidated at the atomic level the intrinsic determinants of the pH-dependent dimerization of the central hydrophobic segment A and related these with the propensity to form amyloid fibrils measured by experimental tools such as atomic force microscopy and fluorescence. The process of A dimerization was evaluated in terms of free energy landscape, side-chain two-dimensional contact probability maps, -sheet registries, potential mean force as a function of inter-chain distances, secondary structure development and radial solvation distributions. We showed that dimerization is a key event in A amyloid formation; it is highly prompted in the order of pH 5.02.98.4 and determines further amyloid growth. The dimerization is governed by a dynamic interplay of hydrophobic, electrostatic and solvation interactions permitting some variability of -sheets at each pH. These results provide atomistic insight into the complex process of molecular recognition detrimental for amyloid growth and pave the way for better understanding of the molecular basis of amyloid diseases.

Effect of atom- and group-based truncations on biomolecules simulated with reaction-field electrostatics

The performance of the reaction-field method of electrostatics is tested in molecular dynamics simulations of protein human interleukin-4 and a short DNA fragment in explicit solvent. Two truncation schemes are considered: one based on the position of atomic charges in water molecules and the other on the position of groups of charges. The group-based truncation leads to the melting of the DNA double helix. In contrast, the atom-based truncation maintains the helical structure intact. Similarly for the protein, the group-based truncation leads to an unfolding at pH 2 while the atom-based truncation produces stable trajectories at low and normal pH, in agreement with experiment. Artificial repulsion between charged residues associated with the group-based truncation is identified as the microscopic reason behind unfolding of the protein. Implications of different truncation schemes in reaction-field simulations of biomolecules are discussed.

Probing the strand orientation and registry alignment in the propagation of amyloid fibrils

Detailed knowledge of the structure and growth mechanism of amyloid fibrils is important for understanding the disease process. Recently, solid-state NMR and other spectroscopic data have revealed the equilibrium organization of the tertiary structure of fibrils formed by various segments of beta-amyloid peptides. A three-step "dock-and-lock" mechanism for fibril growth has been proposed on the basis of the kinetic data. Here we use all-atom replica-exchange molecular dynamics simulations in generalized-Born implicit solvent to probe the mechanism of tertiary structure propagation in fibrils of Abeta(16-22) modeled as an oligomer consisting of two beta-sheets each having four strands. The data show that following association with the oligomer, but before being fully locked onto the existing beta-sheet, the added monomer predominantly samples states with the antiparallel strand orientation, but both in- and one-residue shifted backbone hydrogen bond alignments. The in-register state, which is the experimentally observed equilibrium alignment, is marginally more stable than the registry-shifted one. These results suggest that, following the fast docking step, the added monomer dynamically slides in the backbone registry, and stabilization of the preferential alignment must occur in the second locking step as the monomer becomes fully integrated with the fibril. We also delineate the electrostatic and hydrophobic effects in directing the registry alignment during monomer addition. Surprisingly, the in-register alignment provides both increased cross-strand electrostatic attraction and hydrophobic surface burial. Finally, our data support the notion that side chain hydrophobic burial is a major driving force for beta-sheet assembly.