Conformational Dynamics of the Human Islet Amyloid Polypeptide in a Membrane Environment: Toward the Aggregation Prone Form

Improved Highly Mobile Membrane Mimetic Model for Investigating Protein-Cholesterol Interactions

Cholesterol (CHL) plays an integral role in modulating the function and activity of various mammalian membrane proteins. Due to the slow dynamics of lipids, conventional computational studies of protein-CHL interactions rely on either long-time scale atomistic simulations or coarse-grained approximations to sample the process. A highly mobile membrane mimetic (HMMM) has been developed to enhance lipid diffusion and thus used to facilitate the investigation of lipid interactions with peripheral membrane proteins and, with customized in silico solvents to replace phospholipid tails, with integral membrane proteins. Here, we report an updated HMMM model that is able to include CHL, a nonphospholipid component of the membrane, henceforth called HMMM-CHL. To this end, we had to optimize the effect of the customized solvents on CHL behavior in the membrane. Furthermore, the new solvent is compatible with simulations using force-based switching protocols. In the HMMM-CHL, both improved CHL dynamics and accelerated lipid diffusion are integrated. To test the updated model, we have applied it to the characterization of protein-CHL interactions in two membrane protein systems, the human β2-adrenergic receptor (β2AR) and the mitochondrial voltage-dependent anion channel 1 (VDAC-1). Our HMMM-CHL simulations successfully identified CHL binding sites and captured detailed CHL interactions in excellent consistency with experimental data as well as other simulation results, indicating the utility of the improved model in applications where an enhanced sampling of protein-CHL interactions is desired.

Predicting molecular properties of α-synuclein using force fields for intrinsically disordered proteins

Independent force field validation is an essential practice to keep track of developments and for performing meaningful Molecular Dynamics simulations. In this work, atomistic force fields for intrinsically disordered proteins (IDP) are tested by simulating the archetypical IDP α-synuclein in solution for 2.5 μs. Four combinations of protein and water force fields were tested: ff19SB/OPC, ff19SB/TIP4P-D, ff03CMAP/TIP4P-D, and a99SB-disp/TIP4P-disp, with four independent repeat simulations for each combination. We compare our simulations to the results of a 73 μs simulation using the a99SB-disp/TIP4P-disp combination, provided by D. E. Shaw Research. From the trajectories, we predict a range of experimental observations of α-synuclein and compare them to literature data. This includes protein radius of gyration and hydration, intramolecular distances, NMR chemical shifts, and ³ J-couplings. Both ff19SB/TIP4P-D and a99SB-disp/TIP4P-disp produce extended conformational ensembles of α-synuclein that agree well with experimental radius of gyration and intramolecular distances while a99SB-disp/TIP4P-disp reproduces a balanced α-synuclein secondary structure content. It was found that ff19SB/OPC and ff03CMAP/TIP4P-D produce overly compact conformational ensembles and show discrepancies in the secondary structure content compared to the experimental data.

Molecular View into Preferential Binding of the Factor VII Gla Domain to Phosphatidic Acid

Factor VII (FVII) is a serine protease with a key role in initiating the coagulation cascade. It is part of a family of vitamin K-dependent clotting proteins, which require vitamin K for formation of their specialized membrane-binding domains (Gla domains). Membrane binding of the FVII Gla domain is critical to the activity of FVII, mediating the formation of its complex with other clotting factors. While Gla domains among coagulation factors are highly conserved in terms of amino acid sequence and structure, they demonstrate differential binding specificity toward anionic lipids. Although most Gla domain-containing clotting proteins display a strong preference for phosphatidylserine (PS), it has been demonstrated that FVII and protein C instead bind preferentially to phosphatidic acid (PA). We have developed the first model of the FVII Gla domain bound to PA lipids in membranes containing PA, the highly mobile membrane mimetic model, which accelerates slow diffusion of lipids in molecular dynamics simulations and therefore facilitates the membrane binding process and enhances sampling of lipid interactions. Simulations were performed using atomic level molecular dynamics, requiring a fixed charge to all atoms. The overall charge assigned to each PA lipid for this study was -1. We also developed an additional model of the FVII Gla domain bound to a 1:1 PS/PC membrane and compared the modes of binding of PS and PA lipids to FVII, allowing us to identify potential PA-specific binding sites.

Structural Dissection of the First Events Following Membrane Binding of the Islet Amyloid Polypeptide

The islet amyloid polypeptide (IAPP) is the main constituent of the amyloid fibrils found in the pancreas of type 2 diabetes patients. The aggregation of IAPP is known to cause cell death, where the cell membrane plays a dual role: being a catalyst of IAPP aggregation and being the target of IAPP toxicity. Using ATR-FTIR spectroscopy, transmission electron microscopy, and molecular dynamics simulations we investigate the very first molecular steps following IAPP binding to a lipid membrane. In particular, we assess the combined effects of the charge state of amino-acid residue 18 and the IAPP-membrane interactions on the structures of monomeric and aggregated IAPP. Distinct IAPP-membrane interaction modes for the various IAPP variants are revealed. Membrane binding causes IAPP to fold into an amphipathic α-helix, which in the case of H18K-, and H18R-IAPP readily moves beyond the headgroup region. For all IAPP variants but H18E-IAPP, the membrane-bound helix is an intermediate on the way to amyloid aggregation, while H18E-IAPP remains in a stable helical conformation. The fibrillar aggregates of wild-type IAPP and H18K-IAPP are dominated by an antiparallel β-sheet conformation, while H18R- and H18A-IAPP exhibit both antiparallel and parallel β-sheets as well as amorphous aggregates. Our results emphasize the decisive role of residue 18 for the structure and membrane interaction of IAPP. This residue is thus a good therapeutic target for destabilizing membrane-bound IAPP fibrils to inhibit their toxic actions.

Coexisting Ordered and Disordered Membrane Phases Have Distinct Modes of Interaction with Disease-Associated Oligomers

Ordered membrane domains are thought to influence the attachment and insertion of toxic amyloid oligomers, and consequently, their toxicity. However, if and how the molecular aspects of this interaction depend on the membrane order is poorly understood. Here we measure the affinity, location, and degree of insertion of the small oligomers of hIAPP (human Islet Amyloid Polypeptide, associated with Type II diabetes) at near-physiological concentrations to adjacent domains of a biphasic lipid bilayer. Using simultaneous atomic force, confocal and fluorescence lifetime microscopy (AFM-FLIM), we find that hIAPP oligomers have a nearly 8-fold higher affinity to the disordered domains over the ordered domains. To probe whether this difference indicates different modes of interaction, we measure the change of lifetime of peptide-attached fluorescent labels induced by soluble fluorescence quenchers and also measure the kinetics of localized photobleaching. We find that in the raft-like ordered domains, the oligomers primarily lie on the aqueous interface with limited membrane penetration. However, in the neighboring disordered domains, their C-termini penetrate deeper into the lipid bilayer. We conclude that local membrane order determines not only the affinity but also the mode of interaction of amyloid oligomers, which may have significant implications for disease mechanisms.

Interaction of membrane vesicles with the Pseudomonas functional amyloid protein FapC facilitates amyloid formation

Functional amyloids (FA) are proteins which are evolutionarily optimized to form highly stable fibrillar structures that strengthen the bacterial biofilm matrix. FA such as CsgA (E. coli) and FapC (Pseudomonas) are secreted to the bacterial surface where they integrate into growing fibril structures projecting from the outer membrane. FA are exposed to membrane surfaces in this process, but it remains unclear how membranes can interact with FA and potentially affect the self-assembly. Here we report the effect of different vesicles (DOPG, DMPG, DOPS, DOPC and DMPC) on the kinetics and structural endpoints of FapC fibrillation using various biophysical techniques. Particularly anionic lipids such as DMPG trigger FapC fibrillation, and the protein's second repeat sequence (R2) appears to be important for this interaction. Vesicles formed from phospholipids extracted from three different Pseudomonas strains (Δfap, ΔFapC and pfap) induce FapC fibrillation by accelerating nucleation. The general aggregation inhibitor epigallocatechin gallate (EGCG) inhibits FapC fibrillation by blocking interactions between FapC and vesicles and redirecting FapC monomers to oligomer structures. Our work indicates that biological membranes can contribute significantly to the fibrillation of functional amyloid.

The Effect of Cholesterol on Membrane-Bound Islet Amyloid Polypeptide

Islet amyloid polypeptide (IAPP) is a proposed cause of the decreased beta-cell mass in patients with type-II diabetes. The molecular composition of the cell-membrane is important for regulating IAPP cytotoxicity and aggregation. Cholesterol is present at high concentrations in the pancreatic beta-cells, and in-vitro experiments have indicated that it affects the amyloid formation of IAPP either by direct interactions or by changing the properties of the membrane. In this study we apply atomistic, unbiased molecular dynamics simulations at a microsecond timescale to investigate the effect of cholesterol on membrane bound IAPP. Simulations were performed with various combinations of cholesterol, phosphatidylcholine (PC) and phosphatidylserine (PS) lipids. In all simulations, the helical structure of monomer IAPP was stabilized by the membrane. We found that cholesterol decreased the insertion depth of IAPP compared to pure phospholipid membranes, while PS lipids counteract the effect of cholesterol. The aggregation propensity has previously been proposed to correlate with the insertion depth of IAPP, which we found to decrease with the increased ordering of the lipids induced by cholesterol. Cholesterol is depleted in the vicinity of IAPP, and thus our results suggest that the effect of cholesterol is indirect.

Binding mode of SARS-CoV-2 fusion peptide to human cellular membrane

Infection of human cells by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) relies on its binding to a specific receptor and subsequent fusion of the viral and host cell membranes. The fusion peptide (FP), a short peptide segment in the spike protein, plays a central role in the initial penetration of the virus into the host cell membrane, followed by the fusion of the two membranes. Here, we use an array of molecular dynamics simulations that take advantage of the highly mobile membrane mimetic model to investigate the interaction of the SARS-CoV2 FP with a lipid bilayer representing mammalian cellular membranes at an atomic level and to characterize the membrane-bound form of the peptide. Six independent systems were generated by changing the initial positioning and orientation of the FP with respect to the membrane, and each system was simulated in five independent replicas, each for 300 ns. In 73% of the simulations, the FP reaches a stable, membrane-bound configuration, in which the peptide deeply penetrated into the membrane. Clustering of the results reveals three major membrane-binding modes (binding modes 1-3), in which binding mode 1 populates over half of the data points. Taking into account the sequence conservation among the viral FPs and the results of mutagenesis studies establishing the role of specific residues in the helical portion of the FP in membrane association, the significant depth of penetration of the whole peptide, and the dense population of the respective cluster, we propose that the most deeply inserted membrane-bound form (binding mode 1) represents more closely the biologically relevant form. Analysis of FP-lipid interactions shows the involvement of specific residues, previously described as the "fusion-active core residues," in membrane binding. Taken together, the results shed light on a key step involved in SARS-CoV2 infection, with potential implications in designing novel inhibitors.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

Proteostasis of Islet Amyloid Polypeptide: A Molecular Perspective of Risk Factors and Protective Strategies for Type II Diabetes

The possible link between hIAPP accumulation and β-cell death in diabetic patients has inspired numerous studies focusing on amyloid structures and aggregation pathways of this hormone. Recent studies have reported on the importance of early oligomeric intermediates, the many roles of their interactions with lipid membrane, pH, insulin, and zinc on the mechanism of aggregation of hIAPP. The challenges posed by the transient nature of amyloid oligomers, their structural heterogeneity, and the complex nature of their interaction with lipid membranes have resulted in the development of a wide range of biophysical and chemical approaches to characterize the aggregation process. While the cellular processes and factors activating hIAPP-mediated cytotoxicity are still not clear, it has recently been suggested that its impaired turnover and cellular processing by proteasome and autophagy may contribute significantly toward toxic hIAPP accumulation and, eventually, β-cell death. Therefore, studies focusing on the restoration of hIAPP proteostasis may represent a promising arena for the design of effective therapies. In this review we discuss the current knowledge of the structures and pathology associated with hIAPP self-assembly and point out the opportunities for therapy that a detailed biochemical, biophysical, and cellular understanding of its aggregation may unveil.

Green tea extract EGCG plays a dual role in Aβ42 protofibril disruption and membrane protection: A molecular dynamic study

Amyloid plaques accumulated by the amyloid-β (Aβ) fibrillar aggregates are the major pathological hallmark of the Alzheimer's disease (AD). Inhibiting aggregation and disassembling preformed fibrils of Aβ by natural small molecules have developed into a promising therapeutic strategy for AD. Previous experiments reported that the green tea extract epigallocatechin-3-gallate (EGCG) can disrupt Aβ fibril and reduce Aβ cytotoxicity. The inhibitory ability of EGCG can also be affected by cellular membranes. Thus, it is essential to consider the membrane influences in the investigation of protofibril-disruptive capability of EGCG. Here, we performed multiple all-atom molecular dynamic simulations to investigate the effect of EGCG on the Aβ₄₂ protofibril in the presence of a mixed POPC/POPG (7:3) lipid bilayer and the underlying molecular mechanisms of action. Our simulations show that in the presence of membrane bilayers, EGCG has a preference to bind to the membrane, and this binding alters the binding modes between Aβ₄₂ protofibril and the lipid bilayer, leading to a reduced membrane thinning, indicative of a protective effect of EGCG on the membrane. And EGCG still displays a disruptive effect on Aβ₄₂ protofibril, albeit with a lesser extent of disruption than that in the membrane-free environment. EGCG destabilizes the two hydrophobic core regions (L17-F19-I31 and F4-L34-V36), and disrupts the intrachain K28-A42 salt bridges. Our results reveal that in the presence of lipid bilayers, EGCG plays a dual role in Aβ₄₂ protofibril disruption and membrane protection, suggesting that EGCG could be a potential effective drug candidate for the treatment of AD.

Lipid Modulation in the Formation of β-Sheet Structures. Implications for De Novo Design of Human Islet Amyloid Polypeptide and the Impact on β-Cell Homeostasis

Human islet amyloid polypeptide (hIAPP) corresponds to a 37-residue hormone present in insulin granules that maintains a high propensity to form β-sheet structures during co-secretion with insulin. Previously, employing a biomimetic approach, we proposed a panel of optimized IAPP sequences with only one residue substitution that shows the capability to reduce amyloidogenesis. Taking into account that specific membrane lipids have been considered as a key factor in the induction of cytotoxicity, in this study, following the same design strategy, we characterize the effect of a series of lipids upon several polypeptide domains that show the highest aggregation propensity. The characterization of the C-native segment of hIAPP (residues F₂₃-Y₃₇), together with novel variants F₂₃R and I₂₆A allowed us to demonstrate an effect upon the formation of β-sheet structures. Our results suggest that zwitterionic phospholipids promote adsorption of the C-native segments at the lipid-interface and β-sheet formation with the exception of the F₂₃R variant. Moreover, the presence of cholesterol did not modify this behavior, and the β-sheet structural transitions were not registered when the N-terminal domain of hIAPP (K₁-S₂₀) was characterized. Considering that insulin granules are enriched in phosphatidylserine (PS), the property of lipid vesicles containing negatively charged lipids was also evaluated. We found that these types of lipids promote β-sheet conformational transitions in both the C-native segment and the new variants. Furthermore, these PS/peptides arrangements are internalized in Langerhans islet β-cells, localized in the endoplasmic reticulum, and trigger critical pathways such as unfolded protein response (UPR), affecting insulin secretion. Since this phenomenon was associated with the presence of cytotoxicity on Langerhans islet β-cells, it can be concluded that the anionic lipid environment and degree of solvation are critical conditions for the stability of segments with the propensity to form β-sheet structures, a situation that will eventually affect the structural characteristics and stability of IAPP within insulin granules, thus modifying the insulin secretion.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Formation of α-helical and β-sheet structures in membrane-bound human IAPP monomer and the resulting membrane deformation

Upon adsorption on membrane, human IAPP monomer takes conformational changes from coils to α-helices and β-sheets. The helices inserted and β on surface cause different types of membrane deformation, implying two distinct aggregation mechanisms.

Amyloid assembly and disassembly

Amyloid fibrils are protein homopolymers that adopt diverse cross-β conformations. Some amyloid fibrils are associated with the pathogenesis of devastating neurodegenerative disorders, including Alzheimer's disease and Parkinson's disease. Conversely, functional amyloids play beneficial roles in melanosome biogenesis, long-term memory formation and release of peptide hormones. Here, we showcase advances in our understanding of amyloid assembly and structure, and how distinct amyloid strains formed by the same protein can cause distinct neurodegenerative diseases. We discuss how mutant steric zippers promote deleterious amyloidogenesis and aberrant liquid-to-gel phase transitions. We also highlight effective strategies to combat amyloidogenesis and related toxicity, including: (1) small-molecule drugs (e.g. tafamidis) to inhibit amyloid formation or (2) stimulate amyloid degradation by the proteasome and autophagy, and (3) protein disaggregases that disassemble toxic amyloid and soluble oligomers. We anticipate that these advances will inspire therapeutics for several fatal neurodegenerative diseases.

Summary: This Review showcases important advances in our understanding of amyloid structure, assembly and disassembly, which are inspiring novel therapeutic strategies for amyloid disorders.

Recent computational studies of membrane interaction and disruption of human islet amyloid polypeptide: Monomers, oligomers and protofibrils

The amyloid deposits of human islet amyloid polypeptide (hIAPP) are found in type 2 diabetes patients. hIAPP monomer is intrinsically disordered in solution, whereas it can form amyloid fibrils both in vivo and in vitro. Extensive evidence suggests that hIAPP causes the disruption of cellular membrane, and further induces cytotoxicity and the death of islet β-cells in pancreas. The presence of membrane also accelerates the hIAPP fibril formation. hIAPP oligomers and protofibrils in the early stage of aggregation were reported to be the most cytotoxic, disrupting the membrane integrity and giving rise to the pathological process. The detailed molecular mechanisms of hIAPP-membrane interactions and membrane disruption are complex and remain mostly unknown. Here in this review, we focus on recent computational studies that investigated the interactions of full length and fragmentary hIAPP monomers, oligomers and protofibrils with anionic, zwitterionic and mixed anionic-zwitterionic lipid bilayers. We mainly discuss the binding orientation of monomers at membrane surface, the conformational ensemble and the oligomerization of hIAPP inside membranes, the effect of lipid composition on hIAPP oligomers/protofibrils-membrane interactions, and the hIAPP-induced membrane perturbation. This review provides mechanistic insights into the interactions between hIAPP and lipid bilayers with different lipid composition at an atomistic level, which is helpful to understand the hIAPP cytotoxicity mediated by membrane. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy.

Amyloid growth and membrane damage: Current themes and emerging perspectives from theory and experiments on Aβ and hIAPP

Alzheimer's Disease (AD) and Type 2 diabetes mellitus (T2DM) are two incurable diseases both hallmarked by an abnormal deposition of the amyloidogenic peptides Aβ and Islet Amyloid Polypeptide (IAPP) in affected tissues. Epidemiological data demonstrate that patients suffering from diabetes are at high risk of developing AD, thus making the search for factors common to the two pathologies of special interest for the design of new therapies. Accumulating evidence suggests that the toxic properties of both Aβ or IAPP are ascribable to their ability to damage the cell membrane. However, the molecular details describing Aβ or IAPP interaction with membranes are poorly understood. This review focuses on biophysical and in silico studies addressing these topics. Effects of calcium, cholesterol and membrane lipid composition in driving aberrant Aβ or IAPP interaction with the membrane will be specifically considered. The cross correlation of all these factors appears to be a key issue not only to shed light in the countless and often controversial reports relative to this area but also to gain valuable insights into the central events leading to membrane damage caused by amyloidogenic peptides. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy.

Coupling X-Ray Reflectivity and In Silico Binding to Yield Dynamics of Membrane Recognition by Tim1

The dynamic nature of lipid membranes presents significant challenges with respect to understanding the molecular basis of protein/membrane interactions. Consequently, there is relatively little known about the structural mechanisms by which membrane-binding proteins might distinguish subtle variations in lipid membrane composition and/or structure. We have previously developed a multidisciplinary approach that combines molecular dynamics simulation with interfacial x-ray scattering experiments to produce an atomistic model for phosphatidylserine recognition by the immune receptor Tim4. However, this approach requires a previously determined protein crystal structure in a membrane-bound conformation. Tim1, a Tim4 homolog with distinct differences in both immunological function and sensitivity to membrane composition, was crystalized in a closed-loop conformation that is unlikely to support membrane binding. Here we have used a previously described highly mobile membrane mimetic membrane in combination with a conventional lipid bilayer model to generate a membrane-bound configuration of Tim1 in silico. This refined structure provided a significantly improved fit of experimental x-ray reflectivity data. Moreover, the coupling of the x-ray reflectivity analysis with both highly mobile membrane mimetic membranes and conventional lipid bilayer molecular dynamics simulations yielded a dynamic model of phosphatidylserine membrane recognition by Tim1 with atomic-level detail. In addition to providing, to our knowledge, new insights into the molecular mechanisms that distinguish the various Tim receptors, these results demonstrate that in silico membrane-binding simulations can remove the requirement that the existing crystal structure be in the membrane-bound conformation for effective x-ray reflectivity analysis. Consequently, this refined methodology has the potential for much broader applicability with respect to defining the atomistic details of membrane-binding proteins.

Identification of Key Interactions in the Initial Self-Assembly of Amylin in a Membrane Environment

Islet amyloid polypeptide, also known as amylin, forms aggregates that reduce the amount of insulin-producing cells in patients with type II diabetes mellitus. Much remains unknown about the process of aggregation and cytotoxicity, but it is known that certain cell membrane components can alter the rate of aggregation. Using atomistic molecular dynamics simulations combined with the highly mobile membrane mimetic model incorporating enhanced sampling of lipid diffusion, we investigate interaction of amylin peptides with the membrane components as well as the self-assembly of amylin. Consistent with experimental evidence, we find that an initial membrane-bound α-helical state folds into stable β-sheet structures upon self-assembly. Our results suggest the following mechanism for the initial phase of amylin self-assembly. The peptides move around on the membrane with the positively charged N-terminus interacting with the negatively charged lipid headgroups. When the peptides start to interact, they partly unfold and break some of the contacts with the membrane. The initial interactions between the peptides are dominated by aromatic and hydrophobic interactions. Oligomers are formed showing both intra- and interpeptide β-sheets, initially with interactions mainly in the C-terminal domain of the peptides. Decreasing the pH to 5.5 is known to inhibit amyloid formation. At low pH, His18 is protonated, adding a fourth positive charge at the peptide. With His18 protonated, no oligomerization is observed in the simulations. The additional charge gives a strong midpoint anchoring of the peptides to negatively charged membrane components, and the peptides experience additional interpeptide repulsion, thereby preventing interactions.

Structural Properties of Human IAPP Dimer in Membrane Environment Studied by All-Atom Molecular Dynamics Simulations

The aggregation of human islet amyloid polypeptide (hIAPP) can damage the membrane of the β-cells in the pancreatic islets and induce type 2 diabetes (T2D). Growing evidences indicated that the major toxic species are small oligomers of IAPP. Due to the fast aggregation nature, it is hard to characterize the structures of IAPP oligomers by experiments, especially in the complex membrane environment. On the other side, molecular dynamics simulation can provide atomic details of the structure and dynamics of the aggregation of IAPP. In this study, all-atom bias-exchange metadynamics (BE-Meta) and unbiased molecular dynamics simulations were employed to study the structural properties of IAPP dimer in the membranes environments. A number of intermediates, including α-helical states, β-sheet states, and fully disordered states, are identified. The formation of N-terminal β-sheet structure is prior to the C-terminal β-sheet structure towards the final fibril-like structures. The α-helical intermediates have lower propensity in the dimeric hIAPP and are off-pathway intermediates. The simulations also demonstrate that the β-sheet intermediates induce more perturbation on the membrane than the α-helical and disordered states and thus pose higher disruption ability.