Molecular dynamics simulations of Alzheimer Abeta40 elongation and lateral association

Molecular Recognition between Aβ-Specific Single-Domain Antibody and Aβ Misfolded Aggregates

Aβ is the toxic amyloid polypeptide responsible for Alzheimer's disease (AD). Prevention and elimination of the Aβ misfolded aggregates are the promising therapeutic strategies for the AD treatments. Gammabody, the Aβ-Specific Single-domain (VH) antibody, recognizes Aβ aggregates with high affinity and specificity and reduces their toxicities. Employing the molecular dynamics simulations, we studied diverse gammabody-Aβ recognition complexes to get insights into their structural and dynamic properties and gammabody-Aβ recognitions. Among many heterogeneous binding modes, we focused on two gammabody-Aβ recognition scenarios: recognition through Aβ β-sheet backbone and on sidechain surface. We found that the gammabody primarily uses the complementarity-determining region 3 (CDR3) loop with the grafted Aβ sequence to interact with the Aβ fibril, while CDR1/CDR2 loops have very little contact. The gammabody-Aβ complexes with backbone binding mode are more stable, explaining the gammabody's specificity towards the C-terminal Aβ sequence.

Cross-Seeding Interaction between β-Amyloid and Human Islet Amyloid Polypeptide

Alzheimer's disease (AD) and type 2 diabetes (T2D) are two common protein misfolding diseases. Increasing evidence suggests that these two diseases may be correlated with each other via cross-sequence interactions between β-amyloid peptide (Aβ) associated with AD and human islet amyloid polypeptide (hIAPP) associated with T2D. However, little is known about how these two peptides work and how they interact with each other to induce amyloidogenesis. In this work, we study the effect of cross-sequence interactions between Aβ and hIAPP peptides on hybrid amyloid structures, conformational changes, and aggregation kinetics using combined experimental and simulation approaches. Experimental results confirm that Aβ and hIAPP can interact with each other to aggregate into hybrid amyloid fibrils containing β-sheet-rich structures morphologically similar to pure Aβ and hIAPP. The cross-seeding of Aβ and hIAPP leads to the coexistence of both a retarded process at the initial nucleation stage and an accelerated process at the fibrillization stage, in conjunction with a conformational transition from random structures to α-helix to β-sheet. Further molecular dynamics simulations reveal that Aβ and hIAPP oligomers can efficiently cross-seed each other via the association of two highly similar U-shaped β-sheet structures; thus, conformational compatibility between Aβ and hIAPP aggregates appears to play a key role in determining barriers to cross-seeding. The cross-seeding effects in this work may provide new insights into the molecular mechanisms of interactions between AD and T2D.

Mapping the interactions between the Alzheimer's Aβ-peptide and human serum albumin beyond domain resolution

Human serum albumin (HSA) is a potent inhibitor of Aβ self-association and this novel, to our knowledge, function of HSA is of potential therapeutic interest for the treatment of Alzheimer's disease. It is known that HSA interacts with Aβ oligomers through binding sites evenly partitioned across the three albumin domains and with comparable affinities. However, as of this writing, no information is available on the HSA-Aβ interactions beyond domain resolution. Here, we map the HSA-Aβ interactions at subdomain and peptide resolution. We show that each separate subdomain of HSA domain 3 inhibits Aβ self-association. We also show that fatty acids (FAs) compete with Aβ oligomers for binding to domain 3, but the determinant of the HSA/Aβ oligomer interactions are markedly distinct from those of FAs. Although salt bridges with the FA carboxylate determine the FA binding affinities, hydrophobic contacts are pivotal for Aβ oligomer recognition. Specifically, we identified a site of Aβ oligomer recognition that spans the HSA (494-515) region and aligns with the central hydrophobic core of Aβ. The HSA (495-515) segment includes residues affected by FA binding and this segment is prone to self-associate into β-amyloids, suggesting that sites involved in fibrilization may provide a lead to develop inhibitors of Aβ self-association.

Tanshinones inhibit amyloid aggregation by amyloid-β peptide, disaggregate amyloid fibrils, and protect cultured cells

The misfolding and aggregation of amyloid-β (Aβ) peptides into amyloid fibrils is regarded as one of the causative events in the pathogenesis of Alzheimer's disease (AD). Tanshinones extracted from Chinese herb Danshen (Salvia Miltiorrhiza Bunge) were traditionally used as anti-inflammation and cerebrovascular drugs due to their antioxidation and antiacetylcholinesterase effects. A number of studies have suggested that tanshinones could protect neuronal cells. In this work, we examine the inhibitory activity of tanshinone I (TS1) and tanshinone IIA (TS2), the two major components in the Danshen herb, on the aggregation and toxicity of Aβ1-42 using atomic force microscopy (AFM), thioflavin-T (ThT) fluorescence assay, cell viability assay, and molecular dynamics (MD) simulations. AFM and ThT results show that both TS1 and TS2 exhibit different inhibitory abilities to prevent unseeded amyloid fibril formation and to disaggregate preformed amyloid fibrils, in which TS1 shows better inhibitory potency than TS2. Live/dead assay further confirms that introduction of a very small amount of tanshinones enables protection of cultured SH-SY5Y cells against Aβ-induced cell toxicity. Comparative MD simulation results reveal a general tanshinone binding mode to prevent Aβ peptide association, showing that both TS1 and TS2 preferentially bind to a hydrophobic β-sheet groove formed by the C-terminal residues of I31-M35 and M35-V39 and several aromatic residues. Meanwhile, the differences in binding distribution, residues, sites, population, and affinity between TS1-Aβ and TS2-Aβ systems also interpret different inhibitory effects on Aβ aggregation as observed by in vitro experiments. More importantly, due to nonspecific binding mode of tanshinones, it is expected that tanshinones would have a general inhibitory efficacy of a wide range of amyloid peptides. These findings suggest that tanshinones, particularly TS1 compound, offer promising lead compounds with dual protective role in anti-inflammation and antiaggregation for further development of Aβ inhibitors to prevent and disaggregate amyloid formation.

Molecular interactions of Alzheimer amyloid-β oligomers with neutral and negatively charged lipid bilayers

Interaction of p3 (Aβ(17-42)) peptides with cell membranes is crucial for the understanding of amyloid toxicity associated with Alzheimer's disease (AD). Such p3-membrane interactions are considered to induce the disruption of membrane permeability and integrity, but the exact mechanisms of how p3 aggregates, particularly small p3 oligomers, induce receptor-independent membrane disruption are not yet completely understood. Here, we investigate the adsorption, orientation, and surface interaction of the p3 pentamer with lipid bilayers composed of both pure zwitterionic POPC (palmitoyl-oleoyl-phosphatidylcholine) and mixed anionic POPC-POPG (palmitoyl-oleoyl-phosphatidylglycerol) (3 : 1) lipids using explicit-solvent molecular dynamics (MD) simulations. MD simulation results show that the p3 pentamer has much stronger interactions with mixed POPC-POPG lipids than pure POPC lipids, consistent with experimental observation that Aβ adsorption and fibrillation are enhanced on anionic lipid bilayers. Although electrostatic interactions are main attractive forces to drive the p3 pentamer to adsorb on the bilayer surface, the adsorption of the p3 pentamer on the lipid bilayer with C-terminal β-strands facing toward the bilayer surface is a net outcome of different competitions between p3 peptides-lipid bilayer and ions-p3-bilayer interactions. More importantly, Ca(2+) ions are found to form ionic bridges to associate negatively charged residues of p3 with anionic headgroups of the lipid bilayer, resulting in Aβ-Ca(2+)-PO4(-) complexes. Intensive Ca(2+) bound to the lipid bilayer and Ca(2+) ionic bridges may lead to Ca(2+) hemostasis responsible for neuronal dysfunction and death. This work provides insights into the mutual structure, dynamics, and interactions of both Aβ peptides and lipid bilayers at the atomic level, which expand our understanding of the complex behavior of amyloid-induced membrane disruption.

Binding of Congo red to amyloid protofibrils of the Alzheimer Aβ(9-40) peptide probed by molecular dynamics simulations

Congo red (CR) is a commonly used histological amyloid dye and a weak amyloid inhibitor. There is currently no experimentally available structure of CR bound to an amyloid fibril and the binding modes, and the mechanisms governing its inhibitory and optical properties are poorly understood. In this work, we present the first, to our knowledge, atomistically detailed picture of CR binding to protofibrils of the Alzheimer Aβ(9-40) peptide. We identify three major binding modes, with the primary mode residing in the grooves formed by the β-sheets, and observe a restriction of the torsional rotation of the CR molecule upon binding. Our simulations reveal a novel, to our knowledge, electrostatic steering mechanism that plays an important role in the initial recognition and binding of CR to the positively charged surface residues of the fibril. Our simulations provide new, to our knowledge, insights into the striking spectrophotometric and inhibitory properties of CR. In particular, we show that birefringence upon CR binding is due to the anisotropic orientation of the CR dipoles resulting from the spatial ordering of these molecules in the grooves along the fibril axis. The fluorescent enhancement of the bound CR, in turn, is associated with the torsional restriction of this molecule upon binding.

Conformational basis for asymmetric seeding barrier in filaments of three- and four-repeat tau

Tau pathology in Alzheimer's disease is intimately linked to the deposition of proteinacious filaments, which akin to infectious prions, have been proposed to spread via seeded conversion. Here we use double electron-electron resonance (DEER) spectroscopy in combination with extensive computational analysis to show that filaments of three- (3R) and four-repeat (4R) tau are conformationally distinct. Distance measurements between spin labels in the third repeat, reveal tau amyloid filaments as ensembles of known β-strand-turn-β-strand U-turn motifs. Whereas filaments seeded with 3R tau are structurally homogeneous, filaments seeded with 4R tau are heterogeneous, composed of at least three distinct conformers. These findings establish a molecular basis for the seeding barrier between different tau isoforms and offer a new powerful approach for investigating the composition and dynamics of amyloid fibril ensembles.

Insight into the stability of cross-β amyloid fibril from VEALYL short peptide with molecular dynamics simulation

Amyloid fibrils are found in many fatal neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, type II diabetes, and prion disease. The VEALYL short peptide from insulin has been confirmed to aggregate amyloid-like fibrils. However, the aggregation mechanism of amyloid fibril is poorly understood. Here, we utilized molecular dynamics simulation to analyse the stability of VEALYL hexamer. The statistical results indicate that hydrophobic residues play key roles in stabilizing VEALYL hexamer. Single point and two linkage mutants confirmed that Val1, Leu4, and Tyr5 of VEALYL are key residues. The consistency of the results for the VEALYL oligomer suggests that the intermediate states might be trimer (3-0) and pentamer(3-2). These results can help us to obtain an insight into the aggregation mechanism of amyloid fibril. These methods can be used to study the stability of amyloid fibril from other short peptides.

Cross-seeding and conformational selection between three- and four-repeat human Tau proteins

In Alzheimer's disease and frontotemporal dementias, the microtubule-associated protein Tau forms intracellular paired helical filaments. The filaments can form not only by the full-length human Tau protein, but also by the three repeated (K19) or four repeated (K18) Tau segments. However, of interest, experimentally, K19 can seed K18, but not vice versa. To obtain insight into the cross-seeding between K18 and K19 aggregates, here, K18 and K19 octamers with repeat 3 (R3) in U-shaped, L-shaped, and long straight line-shaped (SL-shape) conformations are assembled into different structures. The simulation results show that K18-8/K19-8 (K18 and K19 assemblies number 8) with R3 in an L shape and K18-9/K19-9 with R3 in an SL shape are highly populated and present the highest structural similarity among all simulated K18 and K19 octamers, suggesting that similar folding of K18/K19 may serve as structural core for the K18-K19 co-assembled heterogeneous filament. We demonstrate that formation of stable R2 and R3 conformations is the critical step for K18 aggregation, and R3 is critical for K19 fibrillization. The different core units in K18 and K19 may create a cross-seeding barrier for the K18 seed to trigger K19 fibril growth because R2 is not available for K19. Our study provides insights into cross-seeding involving heterogeneous structures. The polymorphic nature of protein aggregation could be magnified in the cross-seeding process. If the seeding conformations lead to too much divergence in the energy landscape, it could impede fibril formation. Such an effect could also contribute to the asymmetric barrier between K18 and K19.

Studying interactions by molecular dynamics simulations at high concentration

Molecular dynamics simulations have been used to study molecular encounters and recognition. In recent works, simulations using high concentration of interacting molecules have been performed. In this paper, we consider the practical problems for setting up the simulation and to analyse the results of the simulation. The simulation of beta 2-microglobulin association and the simulation of the binding of hydrogen peroxide by glutathione peroxidase are provided as examples.

Polymorphic structures of Alzheimer's β-amyloid globulomers

Background

Misfolding and self-assembly of Amyloid-β (Aβ) peptides into amyloid fibrils is pathologically linked to the development of Alzheimer's disease. Polymorphic Aβ structures derived from monomers to intermediate oligomers, protofilaments, and mature fibrils have been often observed in solution. Some aggregates are on-pathway species to amyloid fibrils, while the others are off-pathway species that do not evolve into amyloid fibrils. Both on-pathway and off-pathway species could be biologically relevant species. But, the lack of atomic-level structural information for these Aβ species leads to the difficulty in the understanding of their biological roles in amyloid toxicity and amyloid formation.

Methods and Findings

Here, we model a series of molecular structures of Aβ globulomers assembled by monomer and dimer building blocks using our peptide-packing program and explicit-solvent molecular dynamics (MD) simulations. Structural and energetic analysis shows that although Aβ globulomers could adopt different energetically favorable but structurally heterogeneous conformations in a rugged energy landscape, they are still preferentially organized by dynamic dimeric subunits with a hydrophobic core formed by the C-terminal residues independence of initial peptide packing and organization. Such structural organizations offer high structural stability by maximizing peptide-peptide association and optimizing peptide-water solvation. Moreover, curved surface, compact size, and less populated β-structure in Aβ globulomers make them difficult to convert into other high-order Aβ aggregates and fibrils with dominant β-structure, suggesting that they are likely to be off-pathway species to amyloid fibrils. These Aβ globulomers are compatible with experimental data in overall size, subunit organization, and molecular weight from AFM images and H/D amide exchange NMR.

Conclusions

Our computationally modeled Aβ globulomers provide useful insights into structure, dynamics, and polymorphic nature of Aβ globulomers which are completely different from Aβ fibrils, suggesting that these globulomers are likely off-pathway species and explaining the independence of the aggregation kinetics between Aβ globulomers and fibrils.

Unique example of amyloid aggregates stabilized by main chain H-bond instead of the steric zipper: molecular dynamics study of the amyloidogenic segment of amylin wild-type and mutants

Most proteins do not aggregate while in their native functional states. However, they may be disturbed from their native conformation by certain change in the environment, and form unwanted oligomeric or polymeric aggregates. Recent experimental data demonstrate that soluble oligomers of amyloidogenic proteins are responsible for amyloidosis and its cytotoxicity. Human islet amyloid polypeptide (IAPP or amylin) is a 37-residue hormone found as fibrillar deposits in pancreatic extracts of nearly all type II diabetics. In this study we performed in silico mutation analysis to examine the stability of the double layer five strand aggregates formed by heptapeptide NNFGAIL segment from amyline peptide. This segment is one of the shortest fragments that can form amyloid fibrils similar to those formed by the full length peptide. The mutants obtained by single glycine replacement were also studied to investigate the specificity of the dry self-complementary interface between the neighboring β-sheet layers. The molecular dynamics simulations of the aggregates run for 20 ns at 330 K, the degree of the aggregate disassembly was investigated using several geometry analysis tools: the root mean square deviations of the C(α) atoms, root mean square fluctuations per residue, twist angles, interstrand distances, fraction of the secondary structure elements, and number of H-bonds. The analysis shows that most mutations make the aggregates unstable, and their stabilities were dependent to a large extent on the position of replaced residues. Our mutational simulations are in agreement with the pervious experimental observations. We also used free binding energy calculations to determine the role of different components: nonpolar effects, electrostatics and entropy in binding. Nonpolar effects remained consistently more favorable in wild type and mutants reinforcing the importance of hydrophobic effects in protein-protein binding. While entropy systematically opposed binding in all cases, there was no clear trend in the entropy difference between wildtype and glycine mutants. Free energy decomposition shows residues situated at the interface were found to make favorable contributions to the peptide-peptide association. The study of the wild type and mutants in an explicit solvent could provide valuable insight into the future computer guided design efforts for the amyloid aggregation inhibitor.

Molecular structures of quiescently grown and brain-derived polymorphic fibrils of the Alzheimer amyloid abeta9-40 peptide: a comparison to agitated fibrils

The presence of amyloid deposits consisting primarily of Amyloid-β (Aβ) fibril in the brain is a hallmark of Alzheimer's disease (AD). The morphologies of these fibrils are exquisitely sensitive to environmental conditions. Using molecular dynamics simulations combined with data from previously published solid-state NMR experiments, we propose the first atomically detailed structures of two asymmetric polymorphs of the Aβ_9-40 peptide fibril. The first corresponds to synthetic fibrils grown under quiescent conditions and the second to fibrils derived from AD patients' brain-extracts. Our core structure in both fibril structures consists of a layered structure in which three cross-β subunits are arranged in six tightly stacked β-sheet layers with an antiparallel hydrophobic-hydrophobic and an antiparallel polar-polar interface. The synthetic and brain-derived structures differ primarily in the side-chain orientation of one β-strand. The presence of a large and continually exposed hydrophobic surface (buried in the symmetric agitated Aβ fibrils) may account for the higher toxicity of the asymmetric fibrils. Our model explains the effects of external perturbations on the fibril lateral architecture as well as the fibrillogenesis inhibiting action of amphiphilic molecules.

Amyloid diseases are characterized by the presence of amyloid fibrils on organs and tissue in the body. Alzheimer's disease, Parkinson's diseases and Type II Diabetes are all examples of amyloid diseases. Determining the structure of amyloid fibrils is critical for understanding the mechanism of fibril formation as well as for the design of inhibitor molecules that can prevent aggregation. In the case of the Alzheimer Amyloid-β (Aβ) peptide, the structure of fibrils grown under conditions of mechanical agitation has been elucidated from a combination of simulation and experiments. However, the structures of the asymmetric quiescent Aβ fibrils (grown under conditions akin to physiological conditions) and of Alzheimer's brain–derived fibrils are not known. In this paper, we propose the first atomically detailed structures of these two fibrils, using molecular dynamics simulations combined with data from previously published experiments. In additions, we suggest a unifying lateral growth mechanism that explains the increased toxicity of quiescent Aβ fibrils, the effects of external perturbations on fibril lateral architecture and the inhibition mechanism of the small molecule inhibitors on fibril formation.

Atomic-scale simulations confirm that soluble beta-sheet-rich peptide self-assemblies provide amyloid mimics presenting similar conformational properties

The peptide self-assembly mimic (PSAM) from the outer surface protein A (OspA) can form highly stable but soluble beta-rich self-assembly-like structures similar to those formed by native amyloid-forming peptides. However, unlike amyloids that predominantly form insoluble aggregates, PSAMs are highly water-soluble. Here, we characterize the conformations of these soluble beta-sheet-rich assemblies. We simulate PSAMs with different-sized beta-sheets in the presence and absence of end-capping proteins using all-atom explicit-solvent molecular dynamics, comparing the structural stability, conformational dynamics, and association force. Structural and free-energy comparisons among beta-sheets with different numbers of layers and sequences indicate that in similarity to amyloids, the intersheet side chain-side chain interactions and hydrogen bonds combined with intrasheet salt bridges are the major driving forces in stabilizing the overall structural organization. A detailed structural analysis shows that in similarity to amyloid fibrils, all wild-type and mutated PSAM structures display twisted and bent beta-sheets to some extent, implying that a twisted and bent beta-sheet is a general motif of beta-rich assemblies. Thus, our studies indicate that soluble beta-sheet-rich peptide self-assemblies can provide good amyloid mimics, and as such confirm on the atomic scale that they are excellent systems for amyloid studies. These results provide further insight into the usefulness of such mimics for nanostructure design.

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

Short fragments of amyloidogenic proteins are widely used as model systems in studies of amyloid formation. Fragment 11-25 of the amyloid beta protein involved in Alzheimer's disease (Abeta11-25) was recently shown to form amyloid fibrils composed of anti-parallel beta-sheets. Interestingly, fibrils grown under neutral and acidic conditions were seen to possess different registries of their inter-beta-strand hydrogen bonds. In an effort to explain the microscopic origin of this pH dependence, we studied Abeta11-25 fibrils using methods of theoretical modeling. Several structural models were built for fibrils at low and neutral pH levels and these were examined in short molecular dynamics simulations in explicit water. The models that displayed the lowest free energy, as estimated using an implicit solvent model, were selected as representative of the true fibrillar structure. It was shown that the registry of these models agrees well with the experimental results. At neutral pH, the main contribution to the free energy difference between the two registries comes from the electrostatic interactions. The charge group of the carboxy terminus makes a large contribution to these interactions and thus appears to have a critical role in determining the registry.

Mutational analysis and allosteric effects in the HIV-1 capsid protein carboxyl-terminal dimerization domain

The carboxyl-terminal domain (CTD, residues 146-231) of the HIV-1 capsid (CA) protein plays an important role in the CA-CA dimerization and viral assembly of the human immunodeficiency virus type 1. Disrupting the native conformation of the CA is essential for blocking viral capsid formation and viral replication. Thus, it is important to identify the exact nature of the structural changes and driving forces of the CTD dimerization that take place in mutant forms. Here, we compare the structural stability, conformational dynamics, and association force of the CTD dimers for both wild-type and mutated sequences using all-atom explicit-solvent molecular dynamics (MD). The simulations show that Q155N and E159D at the major homology region (MHR) and W184A and M185A at the helix 2 region are energetically less favorable than the wild-type, imposing profound negative effects on intermolecular CA-CA dimerization. Detailed structural analysis shows that three mutants (Q155N, E159D, and W184A) display much more flexible local structures and weaker CA-CA association than the wildtype, primarily due to the loss of interactions (hydrogen bonds, side chain hydrophobic contacts, and pi-stacking) with their neighboring residues. Most interestingly, the MHR that is far from the interacting dimeric interface is more sensitive to the mutations than the helix 2 region that is located at the CA-CA dimeric interface, indicating that structural changes in the distinct motif of the CA could similarly allosterically prevent the CA capsid formation. In addition, the structural and free energy comparison of the five residues shorter CA (151-231) dimer with the CA (146-231) dimer further indicates that hydrophobic interactions, side chain packing, and hydrogen bonds are the major, dominant driving forces in stabilizing the CA interface.