Can molecular dynamics simulations assist in design of specific inhibitors and imaging agents of amyloid aggregation? Structure, stability and free energy predictions for amyloid oligomers of VQIVYK, MVGGVV and LYQLEN

Dose-dependent binding behavior of anthraquinone derivative purpurin interacting with tau-derived peptide protofibril

Alzheimer's disease is hallmarked by microtubule-associated protein tau tangles and amyloid-β plaques. The β-structure propensity of tau inclusions is closely related to the hexapeptide motif VQIVYK (termed PHF6), and disruption of this motif prevents tau aggregation. Small-molecule inhibitors are considered a promising therapeutic strategy, but the molecular mechanisms underlying the correlation between dose and inhibitory effects are still unclear. In this work, we investigated the dose-induced influence of purpurin, an anthraquinone derivative, on the structural stability of the PHF6 fibrillar nucleus by performing microsecond all-atom molecular dynamics simulations in explicit water. The stability of PHF6 protofibrils of different sizes was first examined, and it was found that the structural stability of fibrillar oligomers increases with oligomer size, and that the octamer is the minimal stable nucleus for fibril formation. When purpurin molecules were added to the protofibril octamer at a low purpurin/peptide ratio, they bound to the octamer with different coupling states, and the different states may transition to each of the other states through an uncoupling state or directly through a short-time transition. With increasing purpurin/peptide ratio, purpurins tend to self-aggregate rather than bind to the protein surface. Interestingly, the contacts between individual purpurins and the octamer as a function of the purpurin number show a power-law behavior, which may serve as a useful indicator to reflect the binding efficiency of ligands to proteins in drug screening. The interaction analysis reveals that purpurin prefers to bind to the hydrophilic and aromatic Tyr and has the lowest probability with the hydrophobic Val located in the middle of PHF6. Aromatic stacking plays a key role in the octamer-purpurin interaction, in which the three aromatic rings of purpurin have different contributions. In addition, purpurin shows a remarkable disruptive effect on the protofibril octamer when the molar ratio of purpurin to peptide is 1 : 2; above this ratio, the binding mode and disruption effect of purpurin do not change significantly. Our work provides a detailed picture of the dynamics and interactions of purpurin binding to the PHF6 protofibril and expands the understanding of the dose-induced inhibitory mechanism.

Energy Landscapes and Heat Capacity Signatures for Monomers and Dimers of Amyloid-Forming Hexapeptides

Amyloid formation is a hallmark of various neurodegenerative disorders. In this contribution, energy landscapes are explored for various hexapeptides that are known to form amyloids. Heat capacity (CV) analysis at low temperature for these hexapeptides reveals that the low energy structures contributing to the first heat capacity feature above a threshold temperature exhibit a variety of backbone conformations for amyloid-forming monomers. The corresponding control sequences do not exhibit such structural polymorphism, as diagnosed via end-to-end distance and a dihedral angle defined for the monomer. A similar heat capacity analysis for dimer conformations obtained using basin-hopping global optimisation shows clear features in end-to-end distance versus dihedral correlation plots, where amyloid-forming sequences exhibit a preference for larger end-to-end distances and larger positive dihedrals. These results hold true for sequences taken from tau, amylin, insulin A chain, a de novo designed peptide, and various control sequences. While there is a little overall correlation between the aggregation propensity and the temperature at which the low-temperature CV feature occurs, further analysis suggests that the amyloid-forming sequences exhibit the key CV feature at a lower temperature compared to control sequences derived from the same protein.

Force spectroscopy reveals the presence of structurally modified dimers in transthyretin amyloid annular oligomers

Toxicity in amyloidogenic protein misfolding disorders is thought to involve intermediate states of aggregation associated with the formation of amyloid fibrils. Despite their relevance, the heterogeneity and transience of these oligomers have placed great barriers in our understanding of their structural properties. Among amyloid intermediates, annular oligomers or annular protofibrils have raised considerable interest because they may contribute to a mechanism of cellular toxicity via membrane permeation. Here we investigated, by using AFM force spectroscopy, the structural detail of amyloid annular oligomers from transthyretin (TTR), a protein involved in systemic and neurodegenerative amyloidogenic disorders. Manipulation was performed in situ, in the absence of molecular handles and using persistence length-fit values to select relevant curves. Force curves reveal the presence of dimers in TTR annular oligomers that unfold via a series of structural intermediates. This is in contrast with the manipulation of native TTR that was more often manipulated over length scales compatible with a TTR monomer and without unfolding intermediates. Imaging and force spectroscopy data suggest that dimers are formed by the assembly of monomers in a head-to-head orientation with a nonnative interface along their β-strands. Furthermore, these dimers stack through nonnative contacts that may enhance the stability of the misfolded structure.

Molecular insight into amyloid oligomer destabilizing mechanism of flavonoid derivative 2-(4' benzyloxyphenyl)-3-hydroxy-chromen-4-one through docking and molecular dynamics simulations

Aggregation of amyloid peptide (Aβ) has been shown to be directly related to progression of Alzheimer's disease (AD). Aβ is neurotoxic and its deposition and aggregation ultimately lead to cell death. In our previous work, we reported flavonoid derivative (compound 1) showing promising result in transgenic AD model of Drosophila. Compound 1 showed prevention of Aβ-induced neurotoxicity and neuroprotective efficacy in Drosophila system. However, mechanism of action of compound 1 and its effect on the amyloid is not known. We therefore performed molecular docking and atomistic, explicit-solvent molecular dynamics simulations to investigate the process of Aβ interaction, inhibition, and destabilizing mechanism. Results showed different preferred binding sites of compound 1 and good affinity toward the target. Through the course of 35 ns molecular dynamics simulation, conformations_5 of compound 1 intercalates into the hydrophobic core near the salt bridge and showed major structural changes as compared to other conformations. Compound 1 showed interference with the salt bridge and thus reducing the inter strand hydrogen bound network. This minimizes the side chain interaction between the chains A-B leading to disorder in oligomer. Contact map analysis of amino acid residues between chains A and B also showed lesser interaction with adjacent amino acids in the presence of compound 1 (conformations_5). The study provides an insight into how compound 1 interferes and disorders the Aβ peptide. These findings will further help to design better inhibitors for aggregation of the amyloid oligomer.

Hydration water mobility is enhanced around tau amyloid fibers

The paired helical filaments (PHF) formed by the intrinsically disordered human protein tau are one of the pathological hallmarks of Alzheimer disease. PHF are fibers of amyloid nature that are composed of a rigid core and an unstructured fuzzy coat. The mechanisms of fiber formation, in particular the role that hydration water might play, remain poorly understood. We combined protein deuteration, neutron scattering, and all-atom molecular dynamics simulations to study the dynamics of hydration water at the surface of fibers formed by the full-length human protein htau40. In comparison with monomeric tau, hydration water on the surface of tau fibers is more mobile, as evidenced by an increased fraction of translationally diffusing water molecules, a higher diffusion coefficient, and increased mean-squared displacements in neutron scattering experiments. Fibers formed by the hexapeptide (306)VQIVYK(311) were taken as a model for the tau fiber core and studied by molecular dynamics simulations, revealing that hydration water dynamics around the core domain is significantly reduced after fiber formation. Thus, an increase in water dynamics around the fuzzy coat is proposed to be at the origin of the experimentally observed increase in hydration water dynamics around the entire tau fiber. The observed increase in hydration water dynamics is suggested to promote fiber formation through entropic effects. Detection of the enhanced hydration water mobility around tau fibers is conjectured to potentially contribute to the early diagnosis of Alzheimer patients by diffusion MRI.

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.

Interplay of sequence, topology and termini charge in determining the stability of the aggregates of GNNQQNY mutants: a molecular dynamics study

This study explores the stabilities of single sheet parallel systems of three sequence variants of ¹GNNQQNY⁷, N2D, N2S and N6D, with variations in aggregate size (5–8) and termini charge (charged or neutral). The aggregates were simulated at 300 and 330 K. These mutations decrease amyloid formation in the yeast prion protein Sup35. The present study finds that these mutations cause instability even in the peptide context. The protonation status of termini is found to be a key determinant of stabilities; other determinants are sequence, position of mutation and aggregate size. All systems with charged termini are unstable, whereas both stable and unstable systems are found when the termini are neutral. When termini are charged, the largest stable aggregate for the N2S and N6D systems has 3 to 4 peptides whereas N2D mutation supports oligomers of larger size (5-and 6-mers) as well. Mutation at 2^nd position (N2S and N2D) results in fewer H-bonds at the mutated as well as neighboring (Gly1/Gln4) positions. However, no such effect is found if mutation is at 6^th position (N6D). The effect of Asn→Asp mutation depends on the position and termini charge: it is more destabilizing at the 2^nd position than at the 6^th in case of neutral termini, however, the opposite is true in case of charged termini. Appearance of twist in stable systems and in smaller aggregates formed in unstable systems suggests that twist is integral to amyloid arrangement. Disorder, dissociation or rearrangement of peptides, disintegration or collapse of aggregates and formation of amorphous aggregates observed in these simulations are likely to occur during the early stages of aggregation also. The smaller aggregates formed due to such events have a variety of arrangements of peptides. This suggests polymorphic nature of oligomers and presence of a heterogeneous mixture of oligomers during early stages of aggregation.

Full length amylin oligomer aggregation: insights from molecular dynamics simulations and implications for design of aggregation inhibitors

Amyloid oligomers are considered to play essential roles in the pathogenesis of amyloid-related degenerative diseases including type 2 diabetes. Using an explicit solvent all atomic MD simulation, we explored the stability, conformational dynamics and association force of different single-layer models of the full-length wild-type and glycine mutants of amylin (pentamer) obtained from a recent high resolution fibril model. The RMSF profile shows enhanced flexibility in the disorder (Lys1-Cys7) and turn region (Ser19-Gly23), along with smallest fluctuation at the residues (Asn14-Phe15-Leu16-Val17-His18) of β1 region and (Ala25-Ile26-Leu27-Ser28-Ser29) of the β2 region. We obtained a significant difference in backbone RMSD between the wild-type and the mutants, indicating that mutations affected the stability of the peptide. The RMSD and RMSF profiles indicate the edge and loop residues are the primary contributors to the overall conformational changes. The degree of structural similarity between the oligomers in the simulation and the fibril conformation is proposed as the possible explanation for experimentally observed shortening of the nucleation lag phase of amylin with oligomer seeding. On the basis of structure-stability findings, the β1 and β2 portions are optimal target for further anti-amyloid drug design. The MM-PBSA binding energy calculation reveals the binding of amylin: amylin strands in single layer is dominated by contributions from van der Waals interactions. The non-polar solvation term is also found to be favorable. While the electrostatic interactions and polar solvation energy was found to be favorable for the interaction for the larger aggregate and unfavorable for the smaller aggregates. A per-residue decomposition of the binding free energy has been performed to identify the residues contributing most to the self-association free energy. Residues found in the β-sheet regions were found to be key residue making the largest favorable contributions to the single-layer association. The result from our simulation could be used in rational design of new amylinomimetic agent, amylin aggregation inhibitors and amylin-specific biomarkers.

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.

Interaction of glycine zipper fragments of Aβ-peptides with neuronal nitric oxide synthase: kinetic, thermodynamic and spectrofluorimetric analysis

Five peptide fragments [Aβ(17-21); Aβ(25-29); Aβ(29-33); Aβ(33-37); Aβ(25-37)] of the toxic Aβ(1-40(42)) amyloid peptide were shown to bind with neuronal nitric oxide synthase by means of hydrophobic-hydrophobic forces. The enzyme has a single site for the amyloid peptide binding, which resulted in a quenching of the intrinsic fluorescence of the enzyme. Binding constants determined from Stern-Volmer analysis were between 9×10(-3) and 1.8×10(-2) μM(-1). As temperature increased these binding constants increased reflecting that the interaction of the amyloid peptides with nNOS was endothermic and the quenching was dynamic. Kinetic analysis revealed a non-competitive interaction of the amyloid peptides to the enzyme with inhibitor constants of 5.1 μM for Aβ(17-21) to about 8-12 μM for the other peptides. According to the van't Hoff relationship the thermodynamic parameters, ΔH, ΔS and ΔG for the interaction of the amyloid peptides were all positive and between 41.28 and 77.86 kJ mol(-1)K(-1), 104.92 and 220.82 J mol(-1)K(-1) and 9.92 and 13.13 kJ mol(-1)K(-1), respectively. This suggested that the transition state, created by the amyloid peptide-nNOS complex and generated during the initial stages of Aβ aggregation had to, initially, overcome an activation barrier. Since the ΔG values decreased as temperature increased it not only implied a non-spontaneous interaction but that hydrophobic forces were operative during the binding. By FRET analysis the distance between the donor enzyme and the acceptor amyloid peptide was between 2.7 and 2.8 nm. As the temperature increased from 298 K through 313 K (and higher) the fraction of these tryptophan residues that became exposed increased, to approach a value of 1. There was strong support for the initial interaction being through the glycine zipper regions of Aβ(25-37).

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.

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.

Controlling the aggregation and rate of release in order to improve insulin formulation: molecular dynamics study of full-length insulin amyloid oligomer models

Insulin is a hormone that regulates the physiological glucose level in human blood. Insulin injections are used to treat diabetic patients. The amyloid aggregation of insulin may cause problems during the production, storage, and delivery of insulin formulations. Several modifications to the C-terminus of the B chain have been suggested in order to improve the insulin formulation. The central fragments of the A and B chains (LYQLENY and LVEALYL) have recently been identified as β-sheet-forming regions, and their microcrystalline structures have been used to build a high-resolution amyloid fibril model of insulin. Here we report on a molecular dynamics (MD) study of single-layer oligomers of the full-length insulin which aimed to identify the structural elements that are important for amyloid stability, and to suggest single glycine mutants in the β-sheet region that may improve the formulation. Structural stability, aggregation behavior and the thermodynamics of association were studied for the wild-type and mutant aggregates. A comparison of the oligomers of different sizes revealed that adding strands enhances the internal stability of the wild-type aggregates. We call this "dynamic cooperativity". The secondary structure content and clustering analysis of the MD trajectories show that the largest aggregates retain the fibril conformation, while the monomers and dimers lose their conformations. The degree of structural similarity between the oligomers in the simulation and the fibril conformation is proposed as a possible explanation for the experimentally observed shortening of the nucleation lag phase of insulin with oligomer seeding. Decomposing the free energy into electrostatic, van der Waals and solvation components demonstrated that electrostatic interactions contribute unfavorably to the binding, while the van der Waals and especially solvation effects are favorable for it. A per-atom decomposition allowed us to identify the residues that contribute most to the binding free energy. Residues in the β-sheet regions of chains A and B were found to be the key residues as they provided the largest favorable contributions to single-layer association. The positive ∆∆G (mut) values of 37.3 to 1.4 kcal mol(-1) of the mutants in the β-sheet region indicate that they have a lower tendency to aggregate than the wild type. The information obtained by identifying the parts of insulin molecules that are crucial to aggregate formation and stability can be used to design new analogs that can better control the blood glucose level. The results of our simulation may help in the rational design of new insulin analogs with a decreased propensity for self-association, thus avoiding injection amyloidosis. They may also be used to design new fast-acting and delayed-release insulin formulations.