Structure and Thermodynamics of Amylin Dimer Studied by Hamiltonian-Temperature Replica Exchange Molecular Dynamics Simulations

Zinc stabilization of prefibrillar oligomers of human islet amyloid polypeptide

The aggregation of human islet amyloid polypeptide (hIAPP) has been linked to beta-cell death in type II diabetes. Zinc present in secretory granules has been shown to affect this aggregation. A combination of EXAFS, NMR, and AFM experiments shows that the influence of zinc is most likely due to the stabilization of prefibrillar aggregates of hIAPP.

Structural similarities and differences between amyloidogenic and non-amyloidogenic islet amyloid polypeptide (IAPP) sequences and implications for the dual physiological and pathological activities of these peptides

IAPP, a 37 amino-acid peptide hormone belonging to the calcitonin family, is an intrinsically disordered protein that is coexpressed and cosecreted along with insulin by pancreatic islet β-cells in response to meals. IAPP plays a physiological role in glucose regulation; however, in certain species, IAPP can aggregate and this process is linked to β-cell death and Type II Diabetes. Using replica exchange molecular dynamics with extensive sampling (16 replicas per sequence and 600 ns per replica), we investigate the structure of the monomeric state of two species of aggregating peptides (human and cat IAPP) and two species of non-aggregating peptides (pig and rat IAPP). Our simulations reveal that the pig and rat conformations are very similar, and consist of helix-coil and helix-hairpin conformations. The aggregating sequences, on the other hand, populate the same helix-coil and helix-hairpin conformations as the non-aggregating sequence, but, in addition, populate a hairpin structure. Our exhaustive simulations, coupled with available peptide-activity data, leads us to a structure-activity relationship (SAR) in which we propose that the functional role of IAPP is carried out by the helix-coil conformation, a structure common to both aggregating and non-aggregating species. The pathological role of this peptide may have multiple origins, including the interaction of the helical elements with membranes. Nonetheless, our simulations suggest that the hairpin structure, only observed in the aggregating species, might be linked to the pathological role of this peptide, either as a direct precursor to amyloid fibrils, or as part of a cylindrin type of toxic oligomer. We further propose that the helix-hairpin fold is also a possible aggregation prone conformation that would lead normally non-aggregating variants of IAPP to form fibrils under conditions where an external perturbation is applied. The SAR relationship is used to suggest the rational design of therapeutics for treating diabetes.

IAPP, a 37 amino-acid peptide hormone belonging to the calcitonin family, is an intrinsically disordered peptide produced along with insulin by pancreatic islet β-cells in response to meals. In its functional form, IAPP acts as a synergic partner of insulin to reduce blood glucose. IAPP can, however, also play a pathological role, contributing to Type II diabetes (T2D). Knowledge of the structural nature of the physiological and pathological forms of IAPP will facilitate the rational design of novel drugs for therapeutic treatment of T2D. However, because IAPP does not fold to a single structure, but rather co-exists between multiple functional (and toxic) structures, it is extremely challenging for experimental methods to gain detailed structural information. Using a computational approach, we were able to obtain detailed structures of four IAPP variants and propose a novel structural hypothesis for the two opposing roles of this peptide.

Adsorption and Orientation of Human Islet Amyloid Polypeptide (hIAPP) Monomer at Anionic Lipid Bilayers: Implications for Membrane-Mediated Aggregation

Protein misfolding and aggregation cause serious degenerative diseases, such as Alzheimer's and type II diabetes. Human islet amyloid polypeptide (hIAPP) is the major component of amyloid deposits found in the pancreas of type II diabetic patients. Increasing evidence suggests that β-cell death is related to the interaction of hIAPP with the cellular membrane, which accelerates peptide aggregation. In this study, as a first step towards understanding the membrane-mediated hIAPP aggregation, we investigate the atomic details of the initial step of hIAPP-membrane interaction, including the adsorption orientation and conformation of hIAPP monomer at an anionic POPG lipid bilayer by performing all-atom molecular dynamics simulations. We found that hIAPP monomer is quickly adsorbed to bilayer surface, and the adsorption is initiated from the N-terminal residues driven by strong electrostatic interactions of the positively-charged residues K1 and R11 with negatively-charged lipid headgroups. hIAPP binds parallel to the lipid bilayer surface as a stable helix through residues 7-22, consistent with previous experimental study. Remarkably, different simulations lead to the same binding orientation stabilized by electrostatic and H-bonding interactions, with residues R11, F15 and S19 oriented towards membrane and hydrophobic residues L12, A13, L16 and V17 exposed to solvent. Implications for membrane-mediated hIAPP aggregation are discussed.

Comparative molecular dynamics study of human islet amyloid polypeptide (IAPP) and rat IAPP oligomers

Human islet amyloid polypeptide (hIAPP or amylin) is a causative agent in pancreatic amyloid deposits found in patients with type 2 diabetes. The aggregation of full-length hIAPP(1-37) into small oligomeric species is increasingly believed to be responsible for cell dysfunction and death. However, rat IAPP (rIAPP(1-37)), which differs from hIAPP in only six of 37 residues, loses its aggregation ability to form toxic amyloid species. Atomic details of the effect of sequence on the structure and toxicity between the amyloidogenic, toxic hIAPP peptide and the nonamyloidogenic, nontoxic rIAPP peptide remain unclear. Here, we probe sequence-induced differences in structural stability, conformational dynamics, and driving forces between different hIAPP and rIAPP polymorphic forms from monomer to pentamer using molecular dynamics simulations. Simulations show that hIAPP forms from trimer to pentamer exhibit high structural stability with well-preserved in-register parallel β-sheet and the U-bend conformation. The hIAPP trimer appears to be a smallest minimal seed in solution. The stabilities of parallel hIAPP oligomers increase with the number of peptides. Conversely, replacement of hIAPP sequence by rIAPP sequence causes a significant loss of favorable interpeptide interactions in all rIAPP oligomers, destabilizing the C-terminal β-sheet, turn conformation, and overall stability. A less β-sheet-rich structure and a disturbed U-shaped topology exert a large energy penalty on the self-assemble of the rIAPP peptides into highly ordered, in-register β-sheet-rich protofibrils and fibrils, which explains the nonamyloidogenic activity of rIAPP. Moreover, the absence of interior water within the U-turn region in the well-packed higher-order hIAPP oligomers, not in the poorly packed rIAPP oligomers, also stabilizes peptide association. This work provides atomic details of the sequence-structure relationship between the amyloidogenic hIAPP and its analogues such as the nonamyloidogenic rIAPP and some mutants, which could help in the development of novel therapeutic agents to block the formation of toxic hIAPP oligomeric species for type 2 diabetes.

Advanced replica-exchange sampling to study the flexibility and plasticity of peptides and proteins

Molecular dynamics (MD) simulations are ideally suited to investigate protein and peptide plasticity and flexibility simultaneously at high spatial (atomic) and high time resolution. However, the applicability is still limited by the force field accuracy and by the maximum simulation time that can be routinely achieved in current MD simulations. In order to improve the sampling the replica-exchange (REMD) methodology has become popular and is now the most widely applied advanced sampling approach. Many variants of the REMD method have been designed to reduce the computational demand or to enhance sampling along specific sets of conformational variables. An overview on recent methodological advances and discussion of specific aims and advantages of the approaches will be given. Applications in the area of free energy simulations and advanced sampling of intrinsically disordered peptides and proteins will also be discussed. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

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

The structural characterization of amyloid fibers is one of the most investigated areas in structural biology. Recently, protofibril models for amylin, i.e., the 37-residue human islet amyloid polypeptide or hIAPP were suggested by two groups based on NMR (Biochemistry 2007, 46, 13505-13522) and X-ray (Protein Sci. 2008, 17, 1467-1474) techniques. However, there are significant differences in the two models which maybe originate from the polymorphic nature of amylin fibrils. To obtain further insights into the packing and stability features of the different models, we performed a series of molecular dynamics simulations on them. Our analysis showed that even pairs of β-sheets composed of a limited number of β-strands are stable in the 100-ns simulations, which suggests that steric zipper interactions at a β-sheet-β-sheet interface strongly contribute to the stability of these amyloid aggregates. For both models, outer strands are more flexible, which might coincide with the dynamical requirement that outer strands act as growing sites facilitating conformational changes of new incoming chains. Moreover, simulation results showed that the X-ray models are structurally more compact than the NMR models and have more intimate patterns, which lead to more rigid amyloid models. As a result, the X-ray models are energetically more stable than the NMR models. Further modeling analyses verify the most likely amylin fibril model among both NMR and X-ray models. Upon further study of the force-induced dissociation of a single chain from the protofibrils, the binding energy and the mechanical stability of the fibril models are revealed. On these bases, it is possible to reconcile the crystallographic and the NMR data on the basic amylin fiber unit.

Idealized models of protofilaments of human islet amyloid polypeptide

Fibrils formed by assembly of human islet amyloid polypeptide (hIAPP) are found in most patients with type II diabetes. Structurally, these fibrils are composed of multiple protofilaments and are characterized by extended beta sheets, variable helical twists, and different morphologies. We have previously derived models for the hIAPP protofilament using simulations constrained by data from EPR spectroscopy. In the current work, these models were used as a basis for generating idealized hIAPP protofilaments with symmetrical geometrical properties using a new algorithm, MFIBRIL. We show good agreement of the idealized protofilaments with experimental data for amino acid side chain orientations and geometrical features including the inter-β sheet distance and the protofilament radius. These idealized protofilaments can be used in MFIBRIL to generate fibril models that may be experimentally testable at the molecular level. MFIBRIL can also be used for building structures of any repetitive molecular assembly starting with a single building block obtained from any source.

Distinct dimerization for various alloforms of the amyloid-beta protein: Aβ(1-40), Aβ(1-42), and Aβ(1-40)(D23N)

The Amyloid-beta protein is related to Alzheimer's disease, and various experiments have shown that oligomers as small as the dimer are cytotoxic. Two alloforms are mainly produced: Aβ(1-40) and Aβ(1-42). They have very different oligomer distributions, and it was recently suggested, from experimental studies, that this variation may originate from structural differences in their dimer structures. Little structural information is available on the Aβ dimer, however, and to complement experimental observations, we simulated the folding of the wild-type Aβ(1-40) and Aβ(1-42) dimers as well as the mutated Aβ(1-40)(D23N) dimer using an accurate coarse-grained force field coupled to Hamiltonian-temperature replica exchange molecular dynamics. The D23N variant impedes the salt-bridge formation between D23 and K28 seen in the wild-type Aβ, leading to very different fibrillation properties and final amyloid fibrils. Our results show that the Aβ(1-42) dimer has a higher propensity than the Aβ(1-40) dimer to form β-strands at the central hydrophobic core (residues 17-21) and at the C-terminal (residues 30-42), which are two segments crucial to the oligomerization of Aβ. The free energy landscape of the Aβ(1-42) dimer is also broader and more complex than that of the Aβ(1-40) dimer. Interestingly, D23N also impacts the free energy landscape by increasing the population of configurations with higher β-strand propensities when compared against Aβ(40). In addition, while Aβ(1-40)(D23N) displays a higher β-strand propensity at the C-terminal, its solvent accessibility does not change with respect to the wild-type sequence. Overall, our results show the strong impact of the two amino acids Ile41-Ala42 and the salt-bridge D23-K28 on the folding of the Aβ dimer.