Stability of Osaka Mutant and Wild-Type Fibril Models

Single amino acid mutations in amyloid-beta (Aβ) peptides can lead to early onset and increased severity of Alzheimer's disease. An example is the Osaka mutation (Aβ1-40E22D), which is more toxic than wild-type Aβ1-40. This mutant quickly forms early stage fibrils, one of the hallmarks of the disease, and these fibrils can even seed fibrilization of wild-type monomers. Using molecular dynamic simulations, we show that because of formation of various intra- and intermolecular salt bridges the Osaka mutant fibrils are more stable than wild-type fibrils. The mutant fibril also has a wider water channel with increased water flow than the wild type. These two observations can explain the higher toxicity and aggregation rate of the Osaka mutant over the wild type.

Stability of Iowa mutant and wild type Aβ-peptide aggregates

Recent experiments indicate a connection between the structure of amyloid aggregates and their cytotoxicity as related to neurodegenerative diseases. Of particular interest is the Iowa Mutant, which causes early-onset of Alzheimer's disease. While wild-type Amyloid β-peptides form only parallel beta-sheet aggregates, the mutant also forms meta-stable antiparallel beta sheets. Since these structural variations may cause the difference in the pathological effects of the two Aβ-peptides, we have studied in silico the relative stability of the wild type and Iowa mutant in both parallel and antiparallel forms. We compare regular molecular dynamics simulations with such where the viscosity of the samples is reduced, which, we show, leads to higher sampling efficiency. By analyzing and comparing these four sets of all-atom molecular dynamics simulations, we probe the role of the various factors that could lead to the structural differences. Our analysis indicates that the parallel forms of both wild type and Iowa mutant aggregates are stable, while the antiparallel aggregates are meta-stable for the Iowa mutant and not stable for the wild type. The differences result from the direct alignment of hydrophobic interactions in the in-register parallel oligomers, making them more stable than the antiparallel aggregates. The slightly higher thermodynamic stability of the Iowa mutant fibril-like oligomers in its parallel organization over that in antiparallel form is supported by previous experimental measurements showing slow inter-conversion of antiparallel aggregates into parallel ones. Knowledge of the mechanism that selects between parallel and antiparallel conformations and determines their relative stability may open new avenues for the development of therapies targeting familial forms of early-onset Alzheimer's disease.

On the lack of polymorphism in Aβ-peptide aggregates derived from patient brains

The amyloid beta (Aβ) oligomers and fibrils that are found in neural tissues of patients suffering from Alzheimer's disease may either cause or contribute to the pathology of the disease. In vitro, these Aβ-aggregates are characterized by structural polymorphism. However, recent solid state NMR data of fibrils acquired post mortem from the brains of two Alzheimer's patients indicate presence of only a single, patient-specific structure. Using enhanced molecular dynamic simulations we investigate the factors that modulate the stability of Aβ-fibrils. We find characteristic differences in molecular flexibility, dynamics of interactions, and structural behavior between the brain-derived Aβ-fibril structure and in vitro models. These differences may help to explain the lack of polymorphism in fibrils collected from patient brains, and have to be taken into account when designing aggregation inhibitors and imaging agents for Alzheimer's disease.

Atomistic mechanism of polyphenol amyloid aggregation inhibitors: molecular dynamics study of Curcumin, Exifone, and Myricetin interaction with the segment of tau peptide oligomer

Amyloid fibrils are highly ordered protein aggregates associated with many diseases affecting millions of people worldwide. Polyphenols such as Curcumin, Exifone, and Myricetin exhibit modest inhibition toward fibril formation of tau peptide which is associated with Alzheimer's disease. However, the molecular mechanisms of this inhibition remain elusive. We investigated the binding of three polyphenol molecules to the protofibrils of an amyloidogenic fragment VQIVYK of tau peptide by molecular dynamics simulations in explicit solvent. We find that polyphenols induce conformational changes in the oligomer aggregate. These changes disrupt the amyloid H bonding, perturbing the aggregate. While the structural evolution of the control oligomer with no ligand is limited to the twisting of the β-sheets without their disassembly, the presence of polyphenol molecule pushes the β-sheets apart, and leads to a loosely packed structure where two of four β-sheets dissociate in each of the three cases considered here. The H-bonding capacity of polyphenols is responsible for the observed behavior. The calculated binding free energies and its individual components enabled better understanding of the binding. Results indicated that the contribution from Van der Waals interactions is more significant than electrostatic contribution to the binding. The findings from this study are expected to assist in the development of aggregation inhibitors. Significant binding between polyphenols and aggregate oligomer identified in our simulations confirms the previous experimental observations in which polyphenols refold the tau peptide without forming covalent bonds.

Stability of amyloid oligomers

Molecular simulations are now commonly used to complement experimental techniques in investigating amyloids and their role in human diseases. In this chapter, we will summarize techniques and approaches often used in amyloid simulations and will present recent success stories. Our examples will be focused on lessons learned from molecular dynamics simulations in aqueous environments that start from preformed aggregates. These studies explore the limitations that arise from the choice of force field, the role of mutations in the growth of amyloid aggregates, segmental polymorphism, and the importance of cross-seeding. Furthermore, they give evidence for potential toxicity mechanisms. We finally discuss the role of molecular simulations in the search for aggregation inhibitors.

Inter-species cross-seeding: stability and assembly of rat-human amylin aggregates

Diseases such as type 2 diabetes, Alzheimer's and Parkinson's share as common feature the accumulation of mis-folded disease-specific protein aggregates into fibrillar structures, or plaques. These fibrils may either be toxic by themselves, or act as reservoirs for smaller cytotoxic oligomers. This suggests to investigate molecules as potential therapeutics that either reduce fibril formation or increase fibril stability. One example is rat amylin, which can inhibit aggregation of human amylin, a hallmark of type 2 diabetes. In the present paper, we use molecular dynamics to compare the stability of various preformed aggregates, built out of either human amylin, rat amylin, or mixtures of both. We considered two types of fibril-like oligomers: a single-layer in-register conformation, and a double-layer conformation in which the first U-shaped layer consists of rat amylin and the second layer of human amylin. Our results explain the weak amyloid-inhibiting properties of rat amylin and suggest that membrane leakage due to pore formation is responsible for the toxicity of rat amylin observed in a recent experiment. Together, our results put in question the use of rat amylin or the similar FDA approved drug pramlintide as an inhibitor of human amylin aggregation. They also point to mixed human-rat amylin fibril-like oligomers as possible model-systems for studies of amyloid formation that involve cross-species transmission.

Mutations and seeding of amylin fibril-like oligomers

Seeding a protein solution with preformed fibrils can dramatically enhance the growth rate of amyloids. As the seeds do not need to be of the same protein, seeding may account for the observed correlations between amyloid diseases. In an effort to understand better the molecular mechanisms behind cross seeding we have studied in silico the effect of mutations on the seeding of amylin fibrils. Our investigations of the structural stability of decamers of wild type amylin peptides, of Y37L mutants, and of heteroassemblies of wild-type and mutant amylin molecules show that the experimentally observed efficient cross seeding can be explained based on similarity in fibril structure of components. We find that amyloids with similar side chains packing at the β-sheet interface are structurally compatible, acting as a good template for the congruent incorporation of homologues peptides. In the Y37L mutants, lack of tyrosine-specific interactions causes significant higher flexibility of the C terminal than observed in the wild-type fibril. This effects elongation of the mutant fibril leading to the longer lag times during aggregation that are observed in experiments. Our study gives guidelines for the design of ligands that could stabilize amylin fibrils.

In silico cross seeding of Aβ and amylin fibril-like oligomers

Recent epidemiological data have shown that patients suffering from Type 2 Diabetes Mellitus have an increased risk to develop Alzheimer's disease and vice versa. A possible explanation is the cross-sequence interaction between Aβ and amylin. Because the resulting amyloid oligomers are difficult to probe in experiments, we investigate stability and conformational changes of Aβ-amylin heteroassemblies through molecular dynamics simulations. We find that Aβ is a good template for the growth of amylin and vice versa. We see water molecules permeate the β-strand-turn-β-strand motif pore of the oligomers, supporting a commonly accepted mechanism for toxicity of β-rich amyloid oligomers. Aiming for a better understanding of the physical mechanisms of cross-seeding and cell toxicity of amylin and Aβ aggregates, our simulations also allow us to identify targets for the rational design of inhibitors against toxic fibril-like oligomers of Aβ and amylin oligomers.

The stability of cylindrin β-barrel amyloid oligomer models-a molecular dynamics study

Small-soluble amyloid oligomers are believed to play a significant role in the pathology of amyloid diseases. Recently, the atomic structure of a toxic oligomer formed by an 11 residue and its tandem repeat was found to have an out-off register antiparallel β-strands in the shape of a β-barrel. In the present article we investigate the effect of mutations in the hydrophobic cores on the structure and dynamic of the β-barrels using all atom multiple molecular dynamics simulations with an explicit solvent. Extending previous experiments with molecular dynamics simulations we systematically test how stability and formation of cylindrin depends on the interplay between hydrophobicity and steric effects of the core residues. We find that strong hydrophobic interactions between geometrically fitting residues keep the strands (both in register and out-off-register interface) in close proximity, which in turn stabilizes the side-chain and main-chain hydrogen bonds, and the salt bridges on the outer surface along the weak out-of-register interface. Our simulations also indicate presence of water molecules in the hydrophobic interior of the cylindrin β-barrel.Proteins 2013.

Folding and association of a homotetrameric protein complex in an all-atom Go model

The 84-residue homotetrameric BBAT1 is one of the smallest stable protein complexes and therefore is a good test system to study the self-assembly of multimeric proteins. We have researched for this protein the interplay between the folding of monomers and their assembly into tetramers. Replica exchange molecular dynamics simulations relying on a Go model are compared with earlier simulations that use the physics-based coarse-grained UNRES model.

Side-chain hydrophobicity and the stability of Aβ₁₆₋₂₂ aggregates

Recent mutagenesis studies using the hydrophobic segment of Aβ suggest that aromatic π-stacking interactions may not be critical for fibril formation. We have tested this conjecture by probing the effect of Leu, Ile, and Ala mutation of the aromatic Phe residues at positions 19 and 20, on the double-layer hexametric chains of Aβ fragment Aβ₁₆₋₂₂ using explicit solvent all-atom molecular dynamics. As these simulations rely on the accuracy of the utilized force fields, we first evaluated the dynamic and stability dependence on various force fields of small amyloid aggregates. These initial investigations led us to choose AMBER99SB-ILDN as force field in multiple long molecular dynamics simulations of 100 ns that probe the stability of the wild-type and mutants oligomers. Single-point and double-point mutants confirm that size and hydrophobicity are key for the aggregation and stability of the hydrophobic core region (Aβ₁₆₋₂₂). This suggests as a venue for designing Aβ aggregation inhibitors the substitution of residues (especially, Phe 19 and 20) in the hydrophobic region (Aβ₁₆₋₂₂) with natural and non-natural amino acids of similar size and hydrophobicity.

Similarity analysis, synthesis, and bioassay of antibacterial cyclic peptidomimetics

The chemical similarity of antibacterial cyclic peptides and peptidomimetics was studied in order to identify new promising cyclic scaffolds. A large descriptor space coupled with cluster analysis was employed to digitize known antibacterial structures and to gauge the potential of new peptidomimetic macrocycles, which were conveniently synthesized by acylbenzotriazole methodology. Some of the synthesized compounds were tested against an array of microorganisms and showed antibacterial activity against Bordetella bronchistepica, Micrococcus luteus, and Salmonella typhimurium.

Structure and dynamics of amyloid-β segmental polymorphisms

It is believed that amyloid-beta (Aβ) aggregates play a role in the pathogenesis of Alzheimer's disease. Aβ molecules form β-sheet structures with multiple interaction sites. This polymorphism gives rise to differences in morphology, physico-chemical property and level of cellular toxicity. We have investigated the conformational stability of various segmental polymorphisms using molecular dynamics simulations and find that the segmental polymorphic models of Aβ retain a U-shaped architecture. Our results demonstrate the importance of inter-sheet side chain-side chain contacts, hydrophobic contacts among the strands (β1 and β2) and of salt bridges in stabilizing the aggregates. Residues in β-sheet regions have smaller fluctuation while those at the edge and loop region are more mobile. The inter-peptide salt bridges between Asp23 and Lys28 are strong compared to intra-chain salt bridge and there is an exchange of the inter-chain salt-bridge with intra-chain salt bridge. As our results suggest that Aβ exists under physiological conditions as an ensemble of distinct segmental polymorphs, it may be necessary to account in the development of therapeutics for Alzheimer's disease the differences in structural stability and aggregation behavior of the various Aβ polymorphic forms.

Alternative packing modes leading to amyloid polymorphism in five fragments studied with molecular dynamics

Amyloid aggregates have been implicated in the pathogenesis of diseases such as type 2 diabetes, Alzheimer's, Parkinson's, and prion disease. Recently determined microcrystal structures of several short peptide segments derived from fibril-forming proteins revealed coexistence of alternative aggregation modes (amyloid polymorphism) formed by the same segment. This polymorphism may help in understanding the influence of the side chain packing on the amyloid stability. Here we use molecular dynamics (MD) simulation to analyze the stability of five pairs of polar and nonpolar polymorphic oligomers. MD simulation shows polymorphs with steric zipper interface containing large polar and/or aromatic side chains (GNNQQNY, and NNQNTF) are more stable than steric zipper interfaces made of small or hydrophobic residues (SSTNGVG, VQIVYK, and MVGGVV). Several geometric analyses revealed that larger sheet to sheet interface of the dry steric zipper through polar Q/N rich side chains holds the sheets together. Mutant simulations (Q/N→G) show substitutions with glycine disrupt the steric zipper, leading to unstable oligomers. Stability of Q/N rich oligomers was found to result from the large average number of hydrogen bonds. The molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method reports the nonpolar component of free energy to be favorable, while electrostatic solvation is unfavorable for β-sheet association. Knowledge of structural properties of these fibrils might be useful for developing therapeutic agents against amyloidoses as well as for developing biomaterials. © 2011 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 98: 131-144, 2012.

Molecular dynamic simulation of wild type and mutants of the polymorphic amyloid NNQNTF segments of elk prion: structural stability and thermodynamic of association

A hexapeptide with amino acid sequence NNQNTF from the elk prion protein forms amyloid fibrils. Here we use molecular dynamic simulations of the oligomers and their single point glycine mutants to study their stability. In an effort to probe the structural stability and association thermodynamic in a realistic environment, all wildtype of NNQNTF polymorphic forms with different size and their corresponding double layer 5 strands single point glycine mutants were subjected to a total of 500 ns of explicit-solvent molecular dynamics (MD) simulation. Our results show that the structural stability of the NNQNTF oligomers increases with increasing the number of β-strands for double layers. Our results also demonstrated that hydrophobic interaction is the principle driving force to stabilize the adjacent β-strands while the steric zipper is responsible for holding the neighboring β-sheet layers together. We used MM-PBSA approach free energy calculations to determine the role of 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 observed trend in the entropy difference between wildtype and glycine mutant. 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 may provide valuable insight for amyloid aggregation inhibitor design efforts.