Molecular dynamics simulation study on the molecular structures of the amylin fibril models

Molecular simulations of proteins: From simplified physical interactions to complex biological phenomena

Molecular dynamics simulation is the most popular computational technique for investigating the structural and dynamical behaviour of proteins, in search of the molecular basis of their function. Far from being a completely settled field of research, simulations are still evolving to best capture the essential features of the atomic interactions that govern a protein's inner motions. Modern force fields are becoming increasingly accurate in providing a physical description adequate to this purpose, and allow us to model complex biological systems under fairly realistic conditions. Furthermore, the use of accelerated sampling techniques is improving our access to the observation of progressively larger molecular structures, longer time scales, and more hidden functional events. In this review, the basic principles of molecular dynamics simulations and a number of key applications in the area of protein science are summarized, and some of the most important results are discussed. Examples include the study of the structure, dynamics and binding properties of 'difficult' targets, such as intrinsically disordered proteins and membrane receptors, and the investigation of challenging phenomena like hydration-driven processes and protein aggregation. The findings described provide an overall picture of the current state of this research field, and indicate new perspectives on the road ahead to the upcoming future of molecular simulations.

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.

Structural characterization and cryo-electron tomography analysis of human islet amyloid polypeptide suggest a synchronous process of the hIAPP1-37 amyloid fibrillation

Revealing the aggregation and fibrillation process of variant amyloid proteins is critical for understanding the molecular mechanism of related amyloidosis diseases. Here we characterized the fibrillation morphology and kinetics of type 2 diabetes (T2D) related human islet amyloid polypeptide (hIAPP_1-37) fibril formation process using negative staining transmission electron microscopy (NS-TEM), cryo-electron microscopy (cryo-EM) analysis, and 3D cryo-electron tomography (cryo-ET) reconstruction, together with circular dichroism (CD) and Thioflavin-T (ThT) assays. Our results showed that various amyloid fibrils can be observed at different time points of hIAPP_1-37 fibrillization process, while the winding of protofibrils presents in different growth stages, which suggests a synchronous process of hIAPP_1-37 amyloid fibrillization. This work provides insights into the understanding of hIAPP_1-37 amyloid aggregation process and the pathogenesis of Type 2 diabetes disease.

Mechanistic perspective and functional activity of insulin in amylin aggregation

This work provides the first-ever complete atomic model of insulin–amylin aggregates, identifying the specific interactions that stabilize the insulin–amylin complex.

Insulin is a key regulatory polypeptide that is secreted from pancreatic β-cells and has several important effects on the synthesis of lipids, regulation of enzymatic activities, blood glucose levels and the prevention of hyperglycemia. Insulin was demonstrated to self-assemble into ordered amyloid fibrils upon repeated injections, although the possible biological significance of the supramolecular structures is enigmatic. Amylin is also an amyloidogenic polypeptide that is secreted from pancreatic β-cells and plays an important role in glycemic regulation preventing post-prandial spikes in blood glucose levels. These two amyloidogenic proteins are secreted together from the pancreas and have the ability to interact and produce insulin–amylin aggregates. So far, the molecular architecture of insulin–amylin complexes at the atomic resolution has been unknown. The current work identifies for the first time the specific π–π interactions between Y16 in insulin and F19 in amylin that contribute to the stability of the insulin–amylin complex, by using experimental and molecular modeling techniques. We performed additional experiments that verify the functional activity of insulin in amylin aggregation. Our findings illustrate for the first time the specific interactions between insulin and amylin aggregates at the atomic resolution and provide a new mechanistic perspective on the effect of insulin on amylin aggregation and may pave the way towards pharmacological intervention in this process.

Insight into a New Binding Site of Zinc Ions in Fibrillar Amylin

Amylin peptides are secreted together with insulin and zinc ions from pancreatic β-cells. Under unknown conditions, the amylin peptides aggregate to produce oligomers and fibrils, and in some cases Zn²⁺ ions can bind to amylin peptides to form Zn²⁺-aggregate complexes. Consequently, these aggregates lead to the death of the β-cells and a decrease in insulin, which is one of the symptoms of type-2 diabetes (T2D). Therefore, it is crucial to investigate the binding sites of the Zn²⁺ ions in fibrillary amylin. It was previously found by in vitro and simulation studies that Zn²⁺ ion binds to two or four His residues in the turn domain of fibrillary amylin. In the current study, we present a new Zn²⁺ binding site in the N-terminus of fibrillary amylin with three different coordination modes. Our simulations showed that Zn²⁺ ions bind to polymorphic amylin fibrils with a preference to bind to four Cys residues rather than two Cys residues of two neighboring amylin monomers. The new binding site leads to conformational changes, increases the number of polymorphic states, and demonstrates the existence of competition between various binding sites. Our study provides insight into the molecular mechanisms through which Zn²⁺ ions that play a critical role in amylin aggregation can bind to amylin and promote amylin aggregation in T2D.

Aspects of structural landscape of human islet amyloid polypeptide

The human islet amyloid polypeptide (hIAPP) co-operates with insulin to maintain glycemic balance. It also constitutes the amyloid plaques that aggregate in the pancreas of type-II diabetic patients. We have performed extensive in silico investigations to analyse the structural landscape of monomeric hIAPP, which is presumed to be intrinsically disordered. For this, we construct from first principles a highly predictive energy function that describes a monomeric hIAPP observed in a nuclear magnetic resonance experiment, as a local energy minimum. We subject our theoretical model of hIAPP to repeated heating and cooling simulations, back and forth between a high temperature regime where the conformation resembles a random walker and a low temperature limit where no thermal motions prevail. We find that the final low temperature conformations display a high level of degeneracy, in a manner which is fully in line with the presumed intrinsically disordered character of hIAPP. In particular, we identify an isolated family of α-helical conformations that might cause the transition to amyloidosis, by nucleation.

Orientations of residues along the β-arch of self-assembled amylin fibril-like structures lead to polymorphism

Amylin is an endocrine hormone peptide that consists of 37 residues and is the main component of extracellular amyloid deposits found in the pancreas of most type 2 diabetes patients. Amylin peptides are self-assembled to form oligomers and fibrils. So far, four different molecular structures of the self-assembled amylin fibrils have been observed experimentally: two ssNMR models and two crystal models. This study reveals, for the first time, that there are four self-assembled amylin forms that differ in the orientations of the side chains along the β-arch and are all derived from the two ssNMR models. The two ssNMR models are composed of these four different self-assembled forms of amylin, and the two crystal models are composed of two different self-assembled forms of amylin. This study illustrates at the atomic level the differences among the four experimental models and proposes eight new models of self-assembled amylin that are also composed of the four different self-assembled forms of amylin. Our results show polymorphism of the self-assembled fibril-like amylin, with a slight preference of some of the newly constructed models over the experimental models. Finally, we propose that two different self-assembled fibril-like forms of amylin can interact to form a new fibril-like amylin. We investigated this argument and found that some fibril-like amylin prefers to interact to form stable fibril-like structures, whereas others disfavor it. Our work provides new insights that may suggest strategies for future pharmacological studies that aim to find ways to ameliorate the interactions between polymorphic oligomers and fibrils of amylin.

Simulations of Protein Aggregation: Insights from Atomistic and Coarse-Grained Models

This Perspective highlights recent computational approaches to protein aggregation, from coarse-grained models to atomistic simulations, using the islet amyloid polypeptide (IAPP) as a case study. We review salient open questions where simulations can make an impact, discuss the successes and challenges met by simulations, and explore new directions.

Cross-sequence interactions between human and rat islet amyloid polypeptides

Human islet amyloid polypeptide (hIAPP) can assemble into toxic oligomers and fibrils, which are associated with cell degeneration and the pathogenesis of type 2 diabetes. Cross-interaction of hIAPP with rat IAPP (rIAPP)--a non-amyloidogenic peptide with high sequence similarity to hIAPP--might influence the aggregation and toxicity of hIAPP. However, the exact role of rIAPP in hIAPP aggregation and toxicity still remains unclear. In this work, we investigated the effect of cross-sequence interactions between full-length hIAPP(1-37) and rIAPP(1-37) on hybrid amyloid structures, aggregation kinetics, and cell toxicity using combined computational and experimental approaches. Experimental results indicate a contrasting role of rIAPP in hIAPP aggregation, in which rIAPP initially inhibits the early aggregation and nuclei formation of hIAPP, but hIAPP seeds can also recruit both hIAPP and rIAPP to form more hybrid fibrils, thus promoting amyloid fibrillation ultimately. The coincubation of hIAPP and rIAPP also decreases cell viability, presumably due to the formation of more toxic hybrid oligomers at the prolonged lag phase. Comparative MD simulations confirm that the cross-sequence interactions between hIAPP and rIAPP stabilize β-sheet structure and thus likely promote their fibrillization. This work provides valuable insights into a critical role of cross-amyloid interactions in protein aggregation.

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.

Effect of the amyloid β hairpin's structure on the handedness of helices formed by its aggregates

Various structural models for amyloid β fibrils have been derived from a variety of experimental techniques. However, these models cannot differentiate between the relative position of the two arms of the β hairpin called the stagger. Amyloid fibrils of various hierarchical levels form left-handed helices composed of β sheets. However it is unclear if positive, negative and zero staggers all form the macroscopic left-handed helices. To address this issue we have conducted extensive molecular dynamics simulations of amyloid β sheets of various staggers and shown that only negative staggers lead to the experimentally observed left-handed helices while positive staggers generate the incorrect right-handed helices. This result suggests that the negative staggers are physiologically relevant structure of the amyloid β fibrils.