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.

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.

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.

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.

The aggregation modes of hexapeptide fragments of Tau, Insulin and Aß peptide (VQIVYK, MVGGVV and LYQLEN) were found from their microcrystalline structures that had been recently resolved by X-ray analysis. The atomic structures reveal a dry self-complementary interface between the neighboring ß-sheet layers, termed "steric zipper". In this study we perform several all-atom molecular dynamics simulations with explicit water to analyze stability of the crystalline fragments of 2-10 hexapeptides each and their analogs with single glycine replacement mutations to investigate the structural stability, aggregation behavior and thermodynamic of the amyloid oligomers. Upon comparing single and double layer models, our results reveal that additional strands contribute significantly to the structural stability of the peptide oligomers for double layer model, while in the case of single layer model the stability decreases (or remains the same in the case of LYQLEN). This is in agreement with the previous studies performed on different types of amyloid models. We also replaced the side-chains participating in the steric zipper interfaces with glycine. None of the mutants were structurally stable compared to the respective wild type model, except for mutants V2G and V6G in MVGGVV2 case. The exception can be explained by structural features of this particular polymorph. The double layer decamer and dodecamer aggregates of the wild type hexapeptides appear to be stable at 300K, which is confirmed by the conservation of high anti-parallel ß-sheet content throughout the whole simulation time. Deletions of the side chains resulted in decline of secondary structure content compared to corresponding wild type indicating that the role of the replaced amino acid in stabilizing the structure. Detailed analysis of the binding energy reveals that stability of these peptide aggregates is determined mainly by the van der Waals and hydrophobic forces that can serve as quantitative measure of shape complementarities between the side chains. This observation implies that interactions among side chains forming the dehydrated steric zipper, rather than among those exposed to water, are the major structural determinant. The electrostatic repulsion destabilizes the studied double layer aggregates in two cases, while stabilizes the other two. Negative total binding free energy indicates that both wild type and mutants complex formation is favorable. However, the mutants complexation is less favorable than the wild type's. The present study provides the atomic level understanding of the aggregation behavior and the driving force for the amyloid aggregates, and could be useful for rational design of amyloid inhibitors and amyloid-specific biomarkers for diagnostic purposes.

Are density functional theory predictions of the Raman spectra accurate enough to distinguish conformational transitions during amyloid formation?

We report density functional theory (DFT) calculations of the Raman spectra for hexapeptides of glutamic acid and lysine in three different conformations (alpha, beta and PPII). The wave numbers of amide I, amide II and amide III bands of all three conformations predicted at B3LYP/6-31G and B3LYP/6-31G* are in good agreement with previously reported experimental values of polyglutamic acid and polylysine. Agreement with experiment improves when polarization functions are included in the basis set. Explicit water molecules, H-bonded to the backbone amide groups were found to be absolutely necessary to obtain this agreement. Our results indicate that DFT is a promising tool for assignment of the spectral data on kinetics of conformational changes for peptides during amyloid formation.