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

Converting the Amyloidogenic Islet Amyloid Polypeptide into a Potent Nonaggregating Peptide Ligand by Side Chain-to-Side Chain Macrocyclization

The islet amyloid polypeptide (IAPP), also known as amylin, is a hormone playing key physiological roles. However, its aggregation and deposition in the pancreatic islets are associated with type 2 diabetes. While this peptide adopts mainly a random coil structure in solution, its secondary conformational conversion into α-helix represents a critical step for receptor activation and contributes to amyloid formation and associated cytotoxicity. Considering the large conformational landscape and high amyloidogenicity of the peptide, as well as the complexity of the self-assembly process, it is challenging to delineate the delicate interplay between helical folding, peptide aggregation, and receptor activation. In the present study, we probed the roles of helical folding on the function-toxicity duality of IAPP by restricting its conformational ensemble through side chain-to-side chain stapling via azide-alkyne cycloaddition. Intramolecular macrocyclization (i; i + 4) constrained IAPP into α-helix and inhibited its aggregation into amyloid fibrils. These helical derivatives slowed down the self-assembly of unmodified IAPP. Site-specific macrocyclization modulated the capacity of IAPP to perturb lipid bilayers and cell plasma membrane and reduced, or even fully inhibited, the cytotoxicity associated with aggregation. Furthermore, the α-helical IAPP analogs showed moderate to high potency toward cognate G protein-coupled receptors. Overall, these results indicate that macrocyclization represents a promising strategy to protect an amyloidogenic peptide hormone from aggregation and associated toxicity, while maintaining high receptor activity.

In silico analysis of SOD1 aggregation inhibition modes of tertiary amine pyrazolone and pyrano coumarin ferulate as ALS drug candidates

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease, the familial form (fALS) of which is often cognate to mutations in the antioxidant enzyme Cu/Zn superoxide dismutase 1 (SOD1) leading to misfolding and aggregation. Two small molecules, a tertiary amine pyrazolone (TAP) and a pyrano coumarin ferulate (PCF) were suggested to be ALS drug candidates following experimental observation of their ability to inhibit SOD1 protein misfolding and aggregation. The present work aims at computational investigation of these experimentally proposed drug candidates to gain insight into their mechanism of SOD1 misfolding and aggregation inhibition. On the basis of molecular docking, molecular dynamics simulation, MM-PBSA and per-residue energy decomposition analysis, we examined the specific interactions of TAP and PCF with three probable binding sites of SOD1, namely, dimeric interface cavity, W32 and, UMP binding sites. Results suggest that the binding of TAP at W32 and at UMP sites are least probable due to absence of any favorable interaction. The binding of TAP to dimeric cavity is also unstable due to strong unfavorable interactions. In case of PCF, binding at the UMP site is least probable while binding at dimeric cavity is accompanied by unfavorable interactions. PCF, however, exhibits stable binding with the W32 binding site of SOD1 by stabilizing the solvent accessible hydrophobic residues, which otherwise would have acted as contact points for aggregation. Thus the results imply that compound PCF functions as an inhibitior of SOD1 misfolding/aggregation through direct interaction with the protein SOD1 at the W32 binding site. However, TAP is likely to act as an inhibitor through a different mechanism rather than direct interaction with the protein SOD1. These results apart from reinforcing previous experimental findings, shed light on the probable mechanism of action of the proposed drug candidates.

Molecular basis of the anchoring and stabilization of human islet amyloid polypeptide in lipid hydroperoxidized bilayers

The molecular structure of membrane lipids is formed by mono- or polyunsaturations on their aliphatic tails that make them susceptible to oxidation, facilitating the incorporation of hydroperoxide (R-OOH) functional groups. Such groups promote changes in both composition and complexity of the membrane significantly modifying its physicochemical properties. Human Langerhans islets amyloid polypeptide (hIAPP) is the main component of amyloid deposits found in the pancreas of patients with type-2 diabetes (T2D). hIAPP in the presence of membranes with oxidized lipid species accelerates the formation of amyloid fibrils or the formation of intermediate oligomeric structures. However, the molecular bases at the initial stage of the anchoring and stabilization of the hIAPP in a hydroperoxidized membrane are not yet well understood. To shed some light on this matter, in this contribution, three bilayer models were modeled: neutral (POPC), anionic (POPS), and oxidized (POPC_OOH), and full atom Molecular Dynamics (MD) simulations were performed. Our results show that the POPC_OOH bilayer increases the helicity in hIAPP when compared to POPC or POPS bilayer. The modification in the secondary structure covers the residues of the so-called amyloidogenic core of the hIAPP. Overall, the hydroperoxidation of the neutral lipids modifies both the anchoring and the stabilization of the peptide hIAPP by reducing the random conformations of the peptide and increasing of hydrogen bond population with the hydroperoxidized lipids.

Aggregation of Lysozyme in the Presence of a Mixed Bilayer of POPC and POPG

Understanding the molecular mechanisms by which amyloidogenic proteins interact with membranes is a challenging task. Amyloid accumulates from many human diseases have been observed to contain membrane lipids. In this work, coarse-grained molecular dynamics simulations have been used to inspect hen egg white lysozyme (HEWL) aggregation and membrane association in the presence of a pure POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) bilayer and a POPC and POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol) mixed bilayer. It was observed that, in both cases, two HEWLs formed aggregates. In the presence of a mixed bilayer, after aggregation, the aggregated system started to interact with the membrane. It has been found that one of the lysozymes which came closer to the mixed bilayer unfolded more. The process of the initial insertion of an aggregated system in the mixed bilayer has been analyzed. The structural rearrangements of the protein and lipids were analyzed as well along the course of the simulation. Although with a pure POPC bilayer, aggregation was observed, the aggregated system moved away from the membrane. We believe that our study will provide considerable insights into lysozyme aggregation in the presence of a membrane environment.

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.

Targeting Human Islet Amyloid Polypeptide Aggregation and Toxicity in Type 2 Diabetes: An Overview of Peptide-Based Inhibitors

Type 2 diabetes (T2D) is a chronic metabolic disease characterized by insulin resistance and a progressive loss of pancreatic islet β-cell mass, which leads to insufficient secretion of insulin and hyperglycemia. Emerging evidence suggests that toxic oligomers and fibrils of human islet amyloid polypeptide (hIAPP) contribute to the death of β-cells and lead to T2D pathogenesis. These observations have opened new avenues for the development of islet amyloid therapies for the treatment of T2D. The peptide-based inhibitors are of great value as therapeutic agents against hIAPP aggregation in T2D owing to their biocompatibility, feasibility of synthesis and modification, high specificity, low toxicity, proteolytic stability (modified peptides), and weak immunogenicity as well as the large size of involved interfaces during self-aggregation of hIAPP. An understanding of what has been done and achieved will provide key insights into T2D pathology and assist in the discovery of more potent drug candidates for the treatment of T2D. In this article, we review various peptide-based inhibitors of hIAPP aggregation, including those derived from the hIAPP sequence and those not based on the sequence, consisting of both natural as well as unnatural amino acids and their derivatives. The present review will be beneficial in advancing the field of peptide medicine for the treatment of T2D.

Human islet amyloid polypeptide (hIAPP) - a curse in type II diabetes mellitus: insights from structure and toxicity studies

The human islet amyloid polypeptide (hIAPP) or amylin, a neuroendocrine peptide hormone, is known to misfold and form amyloidogenic aggregates that have been observed in the pancreas of 90% subjects with Type 2 Diabetes Mellitus (T2DM). Under normal physiological conditions, hIAPP is co-stored and co-secreted with insulin; however, under chronic hyperglycemic conditions associated with T2DM, the overexpression of hIAPP occurs that has been associated with the formation of amyloid deposits; as well as the death and dysfunction of pancreatic β-islets in T2DM. Hitherto, various biophysical and structural studies have shown that during this process of aggregation, the peptide conformation changes from random structure to helix, then to β-sheet, subsequently to cross β-sheets, which finally form left-handed helical aggregates. The intermediates, formed during this process, have been shown to induce higher cytotoxicity in the β-cells by inducing cell membrane disruption, endoplasmic reticulum stress, mitochondrial dysfunction, oxidative stress, islet inflammation, and DNA damage. As a result, several research groups have attempted to target both hIAPP aggregation phenomenon and the destabilization of preformed fibrils as a therapeutic intervention for T2DM management. In this review, we have summarized structural aspects of various forms of hIAPP viz. monomer, oligomers, proto-filaments, and fibrils of hIAPP. Subsequently, cellular toxicity caused by toxic conformations of hIAPP has been elaborated upon. Finally, the need for performing structural and toxicity studies in vivo to fill in the gap between the structural and cellular aspects has been discussed.

Computer simulations of protein-membrane systems

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

Effects of in vivo conditions on amyloid aggregation

One of the grand challenges of biophysical chemistry is to understand the principles that govern protein misfolding and aggregation, which is a highly complex process that is sensitive to initial conditions, operates on a huge range of length- and timescales, and has products that range from protein dimers to macroscopic amyloid fibrils. Aberrant aggregation is associated with more than 25 diseases, which include Alzheimer's, Parkinson's, Huntington's, and type II diabetes. Amyloid aggregation has been extensively studied in the test tube, therefore under conditions that are far from physiological relevance. Hence, there is dire need to extend these investigations to in vivo conditions where amyloid formation is affected by a myriad of biochemical interactions. As a hallmark of neurodegenerative diseases, these interactions need to be understood in detail to develop novel therapeutic interventions, as millions of people globally suffer from neurodegenerative disorders and type II diabetes. The aim of this review is to document the progress in the research on amyloid formation from a physicochemical perspective with a special focus on the physiological factors influencing the aggregation of the amyloid-β peptide, the islet amyloid polypeptide, α-synuclein, and the hungingtin protein.

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.

Atomistic-level study of the interactions between hIAPP protofibrils and membranes: Influence of pH and lipid composition

The pathology of type 2 diabetes mellitus is associated with the aggregation of human islet amyloid polypeptide (hIAPP) and aggregation-mediated membrane disruption. The interactions of hIAPP aggregates with lipid membrane, as well as the effects of pH and lipid composition at the atomic level, remain elusive. Herein, using molecular dynamics simulations, we investigate the interactions of hIAPP protofibrillar oligomers with lipids, and the membrane perturbation that they induce, when they are partially inserted in an anionic dipalmitoyl-phosphatidylglycerol (DPPG) membrane or a mixed dipalmitoyl-phosphatidylcholine (DPPC)/DPPG (7:3) lipid bilayer under acidic/neutral pH conditions. We observed that the tilt angles and insertion depths of the hIAPP protofibril are strongly correlated with the pH and lipid composition. At neutral pH, the tilt angle and insertion depth of hIAPP protofibrils at a DPPG bilayer reach ~52° and ~1.62 nm with respect to the membrane surface, while they become ~77° and ~1.75 nm at a mixed DPPC/DPPG membrane. The calculated tilt angle of hIAPP at DPPG membrane is consistent with a recent chiral sum frequency generation spectroscopic study. The acidic pH induces a smaller tilt angle of ~40° and a shallower insertion depth (~1.24 nm) of hIAPP at the DPPG membrane surface, mainly due to protonation of His18 near the turn region. These differences mainly result from a combination of distinct electrostatic, van der Waals, hydrogen bonding and salt-bridge interactions between hIAPP and lipid bilayers. The hIAPP-membrane interaction energy analysis reveals that besides charged residues K1, R11 and H18, aromatic residues Phe15 and Phe23 also exhibit strong interactions with lipid bilayers, revealing the crucial role of aromatic residues in stabilizing the membrane-bound hIAPP protofibrils. hIAPP-membrane interactions disturb the lipid ordering and the local bilayer thickness around the peptides. Our results provide atomic-level information of membrane interaction of hIAPP protofibrils, revealing pH-dependent and membrane-modulated hIAPP aggregation at the early stage.

Molecular dynamics simulations of β2-microglobulin interaction with hydrophobic surfaces

Hydrophobic surfaces are known to adsorb and unfold proteins, a process that has been studied only for a few proteins. Here we address the interaction of β2-microglobulin, a paradigmatic protein for the study of amyloidogenesis, with hydrophobic surfaces. A system with 27 copies of the protein surrounded by a model cubic hydrophobic box is studied by implicit solvent molecular dynamics simulations. Most proteins adsorb on the walls of the box without major distortions in local geometry, whereas free molecules maintain proper structures and fluctuations as observed in explicit solvent molecular dynamics simulations. The major conclusions from the simulations are as follows: (i) the adopted implicit solvent model is adequate to describe protein dynamics and thermodynamics; (ii) adsorption occurs readily and is irreversible on the simulated timescale; (iii) the regions most involved in molecular encounters and stable interactions with the walls are the same as those that are important in protein-protein and protein-nanoparticle interactions; (iv) unfolding following adsorption occurs at regions found to be flexible by both experiments and simulations; (v) thermodynamic analysis suggests a very large contribution from van der Waals interactions, whereas unfavorable electrostatic interactions are not found to contribute much to adsorption energy. Surfaces with different degrees of hydrophobicity may occur in vivo. Our simulations show that adsorption is a fast and irreversible process which is accompanied by partial unfolding. The results and the thermodynamic analysis presented here are consistent with and rationalize previous experimental work.

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.

Membrane Interactions of hIAPP Monomer and Oligomer with Lipid Membranes by Molecular Dynamics Simulations

Interaction of human islet amyloid polypeptide (hIAPP) peptides with cell membrane is crucial for the understanding of amyloid toxicity associated with Type II diabetes (T2D). While it is known that the hIAPP-membrane interactions are considered to promote hIAPP aggregation into fibrils and induce membrane disruption, the membrane-induced conformation, orientation, aggregation, and adsorption behaviors of hIAPP peptides have not been well understood at the atomic level. Herein, we perform all-atom explicit-water molecular dynamics (MD) simulations to study the adsorption, orientation, and surface interaction of hIAPP aggregates with different sizes (monomer to tetramer) and conformations (monomer with α-helix and tetramer with β-sheet-rich U-turn) upon adsorption on the lipid bilayers composed of both pure zwitterionic POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and mixed anionic POPC/POPE (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine) (3:1) lipids. MD simulation results show that hIAPP monomer with α-helical conformation and hIAPP pentamer with β-sheet conformation can adsorb on both POPC and POPC/POPE bilayers via a preferential orientation of N-terminal residues facing toward the bilayer surface. The hIAPP aggregates show stronger interactions with mixed POPC/POPE lipids than pure POPC lipids, consistent with experimental observation that hIAPP adsorption and fibrililation are enhanced on mixed lipid bilayers. While electrostatic interactions are main attractive forces to drive the hIAPP aggregates to adsorb on both bilayers, the introduction of the more hydrophilic head groups of POPE lipids further promote the formation of the interfacial hydrogen bonds. Complement to our previous studies of hIAPP aggregates in bulk solution, this computational work increases our knowledge about the mechanism of amyloid peptide-membrane interactions that is central to the understanding the progression of all amyloid diseases.

Molecular dynamics simulation of the membrane binding and disruption mechanisms by antimicrobial scorpion venom-derived peptides

Pandinin 2 (Pin2) is an alpha-helical polycationic peptide, identified and characterized from venom of the African scorpion Pandinus imperator with high antimicrobial activity against Gram-positive bacteria and less active against Gram-negative bacteria, however it has demonstrated strong hemolytic activity against sheep red blood cells. In the chemically synthesized Pin2GVG analog, the GVG motif grants it low hemolytic activity while keeping its antimicrobial activity. In this work, we performed 12 μs all-atom molecular dynamics simulation of the antimicrobial peptides (AMPs) Pin2 and Pin2GVG to explore their adsorption mechanism and the role of their constituent amino acid residues when interacting with pure POPC and pure POPG membrane bilayers. Starting from an α-helical conformation, both AMPs are attracted at different rates to the POPC and POPG bilayer surfaces due to the electrostatic interaction between the positively charged amino acid residues and the charged moieties of the membranes. Since POPG is an anionic membrane, the PAMs adhesion is stronger to the POPG membrane than to the POPC membrane and they are stabilized more rapidly. This study reveals that, before the insertion begins, Pin2 and Pin2GVG remained partially folded in the POPC surface during the first 300 and 600 ns, respectively, while they are mostly unfolded in the POPG surface during most of the simulation time. The unfolded structures provide for a large number of intermolecular hydrogen bonds and stronger electrostatic interactions with the POPG surface. The results show that the aromatic residues at the N-terminus of Pin2 initiate the insertion process in both POPC and POPG bilayers. As for Pin2GVG in POPC the C-terminus residues seem to initiate the insertion process while in POPG this process seems to be slowed down due to a strong electrostatic attraction. The membrane conformational effects upon PAMs binding are measured in terms of the area per lipid and the contact surface area. Several replicas of the systems lead to the same observations.

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

Human islet amyloid polypeptide (hIAPP) is a 37-residue peptide hormone, which upon misfolding changes from the physiologically active monomer into pathological amyloid fibril aggregates in the pancreas of type 2 diabetes mellitus patients. During this process, the insulin-producing pancreatic β-cells are damaged; however, the underlying mechanism of this mode of cytotoxicity remains elusive. It is known that anionic lipids accelerate amyloid fibril formation, implicating the importance of the cellular membrane in the process, and that a pH close to the level in the β-cell secretory granules (pH 5.5) inhibits amyloid fibril formation. Using all-atom molecular dynamics simulations, we have investigated the membrane-associated monomer state of α-helical hIAPP, analyzed specific interactions of hIAPP with a mixed anionic-zwitterionic lipid membrane and examined the influence of pH on the structure and dynamics of hIAPP and its interaction with the membrane. We find that hIAPP primarily interacts with the membrane by forming favorable interactions between anionic lipids and the positively charged residues in the N-terminal part of the peptide. Rationalizing experimental findings, the simulations show that the N-terminal part of the peptide interacts with the membrane in the lipid headgroup region. At neutral pH, the C-terminal part of the peptide, which contains the residues that initiate fibril formation, displays a highly dynamic, unfolded state, which interacts with the membrane significantly less than the N-terminal part. Such an unfolded form can be proposed to contribute to the acceleration of fibril formation. At low pH, protonation of His18 mediates a stronger interaction of the C-terminal part with the membrane, resulting in the immobilization of the C-terminal part on the membrane surface that might constitute a mechanism by which low pH inhibits fibril formation.

Role of charged lipids in membrane structures - Insight given by simulations

Lipids and proteins are the main components of cell membranes. It is becoming increasingly clear that lipids, in addition to providing an environment for proteins to work in, are in many cases also able to modulate the structure and function of those proteins. Particularly charged lipids such as phosphatidylinositols and phosphatidylserines are involved in several examples of such effects. Molecular dynamics simulations have proved an invaluable tool in exploring these aspects. This so-called computational microscope can provide both complementing explanations for the experimental results and guide experiments to fruitful directions. In this paper, we review studies that have utilized molecular dynamics simulations to unravel the roles of charged lipids in membrane structures. We focus on lipids as active constituents of the membranes, affecting both general membrane properties as well as non-lipid membrane components, mainly proteins. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Binding Orientations and Lipid Interactions of Human Amylin at Zwitterionic and Anionic Lipid Bilayers

Increasing evidence suggests that the interaction of human islet amyloid polypeptide (hIAPP) with lipids may facilitate hIAPP aggregation and cause the death of pancreatic islet β-cells. However, the detailed hIAPP-membrane interactions and the influences of lipid compositions are unclear. In this study, as a first step to understand the mechanism of membrane-mediated hIAPP aggregation, we investigate the binding behaviors of hIAPP monomer at zwitterionic palmitoyloleoyl-phosphatidylcholine (POPC) bilayer by performing atomistic molecular dynamics simulations. The results are compared with those of hIAPP at anionic palmitoyloleoyl-phosphatidylglycerol (POPG) bilayers. We find that the adsorption of hIAPP to POPC bilayer is mainly initiated from the C-terminal region and the peptide adopts a helical structure with multiple binding orientations, while the adsorption to POPG bilayer is mostly initiated from the N-terminal region and hIAPP displays one preferential binding orientation, with its hydrophobic residues exposed to water. hIAPP monomer inserts into POPC lipid bilayers more readily than into POPG bilayers. Peptide-lipid interaction analyses show that the different binding features of hIAPP at POPC and POPG bilayers are attributed to different magnitudes of electrostatic and hydrogen-bonding interactions with lipids. This study provides mechanistic insights into the different interaction behaviors of hIAPP with zwitterionic and anionic lipid bilayers.

Spontaneous adsorption of coiled-coil model peptides K and E to a mixed lipid bilayer

A molecular description of the lipid-protein interactions underlying the adsorption of proteins to membranes is crucial for understanding, for example, the specificity of adsorption or the binding strength of a protein to a bilayer, or for characterizing protein-induced changes of membrane properties. In this paper, we extend an automated in silico assay (DAFT) for binding studies and apply it to characterize the adsorption of the model fusion peptides E and K to a mixed phospholipid/cholesterol membrane using coarse-grained molecular dynamics simulations. In addition, we couple the coarse-grained protocol to reverse transformation to atomistic resolution, thereby allowing to study molecular interactions with high detail. The experimentally observed differential binding of the peptides E and K to membranes, as well as the increased binding affinity of helical over unstructered peptides, could be well reproduced using the polarizable Martini coarse-grained (CG) force field. Binding to neutral membranes is shown to be dominated by initial binding of the positively charged N-terminus to the phospholipid headgroup region, followed by membrane surface-aligned insertion of the peptide at the interface between the hydrophobic core of the membrane and its polar headgroup region. Both coarse-grained and atomistic simulations confirm a before hypothesized snorkeling of lysine side chains for the membrane-bound state of the peptide K. Cholesterol was found to be enriched in peptide vicinity, which is probably of importance for the mechanism of membrane fusion. The applied sequential multiscale method, using coarse-grained simulations for the slow adsorption process of peptides to membranes followed by backward transformation to atomistic detail and subsequent atomistic simulations of the preformed peptide-lipid complexes, is shown to be a versatile approach to study the interactions of peptides or proteins with biomembranes.

Distinct helix propensities and membrane interactions of human and rat IAPP(1-19) monomers in anionic lipid bilayers

Islet amyloid polypeptide, IAPP or amylin, is a 37-residue peptide hormone coexpressed with insulin by pancreatic β-cells. The aggregation of human IAPP (hIAPP) into amyloid deposits is associated with type II diabetes. Substantial evidence suggests that the interaction of anionic membranes with hIAPP may facilitate peptide aggregation and the N-terminal 1-19 fragment (IAPP(1-19)) plays an important role in peptide-membrane interaction. As a first step to understand how structural differences between human and rat IAPP peptides in membranes may influence the later oligomerization process, we have investigated the structures and orientations of hIAPP(1-19) and the less toxic rIAPP(1-19) (i.e., the H18R mutant of hIAPP(1-19)) monomers in anionic POPG bilayers by performing replica exchange molecular dynamics (REMD) simulations. On the basis of ∼20 μs REMD simulations started from a random coil conformation of the peptide placed in water, we find that unfolded h(r)IAPP(1-19) can insert into the bilayers and the membrane-bound peptide stays mainly at the lipid head-tail interface. hIAPP(1-19) displays higher helix propensity than rIAPP(1-19), especially in the L12-L16 region. The helix is oriented parallel to the bilayer surface and buried in the membrane 0.3-0.8 nm below the phosphorus atoms, consistent with previous electron paramagnetic resonance data. The helical conformation is an amphiphilic helix with its hydrophilic and hydrophobic faces pointing, respectively, to the lipid head and tail regions. The H18R substitution enhances the electrostatic interactions of IAPP(1-19) with the membrane, while it weakens the intrapeptide interactions crucial for helix formation, thus leading to lower helix propensity of rIAPP(1-19). Implications of our simulation results on the membrane-mediated IAPP(1-19) oligomerization are discussed.

The Effects of Lipid Membranes, Crowding and Osmolytes on the Aggregation, and Fibrillation Propensity of Human IAPP

Type 2 diabetes mellitus (T2DM) is an age-related and metabolic disease. Its development is hallmarked, among others, by the dysfunction and degeneration of β-cells of the pancreatic islets of Langerhans. The major pathological characteristic thereby is the formation of extracellular amyloid deposits consisting of the islet amyloid polypeptide (IAPP). The process of human IAPP (hIAPP) self-association, and the intermediate structures formed as well as the interaction of hIAPP with membrane systems seem to be, at least to a major extent, responsible for the cytotoxicity. Here we present a summary and comparison of the amyloidogenic propensities of hIAPP in bulk solution and in the presence of various neutral and charged lipid bilayer systems as well as biological membranes. We also discuss the cellular effects of macromolecular crowding and osmolytes on the aggregation pathway of hIAPP. Understanding the influence of different cellular factors on hIAPP aggregation will provide more insight into the onset of T2DM and help to develop novel therapeutic strategies.

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.

Atomistic mechanisms of huntingtin N-terminal fragment insertion on a phospholipid bilayer revealed by molecular dynamics simulations

The huntingtin protein is characterized by a segment of consecutive glutamines (Q(N)) that is responsible for its fibrillation. As with other amyloid proteins, misfolding of huntingtin is related to Huntington's disease through pathways that can involve interactions with phospholipid membranes. Experimental results suggest that the N-terminal 17-amino-acid sequence (htt(NT)) positioned just before the Q(N) region is important for the binding of huntingtin to membranes. Through all-atom explicit solvent molecular dynamics simulations, we unveil the structure and dynamics of the htt(NT)Q(N) fragment on a phospholipid membrane at the atomic level. We observe that the insertion dynamics of this peptide can be described by four main steps-approach, reorganization, anchoring, and insertion-that are very diverse at the atomic level. On the membrane, the htt(NT) peptide forms a stable α-helix essentially parallel to the membrane with its nonpolar side-chains-mainly Leu-4, Leu-7, Phe-11 and Leu-14-positioned in the hydrophobic core of the membrane. Salt-bridges involving Glu-5, Glu-12, Lys-6, and Lys-15, as well as hydrogen bonds involving Thr-3 and Ser-13 with the phospholipids also stabilize the structure and orientation of the htt(NT) peptide. These observations do not significantly change upon adding the Q(N) region whose role is rather to provide, through its hydrogen bonds with the phospholipids' head group, a stable scaffold facilitating the partitioning of the htt(NT) region in the membrane. Moreover, by staying accessible to the solvent, the amyloidogenic Q(N) region could also play a key role for the oligomerization of htt(NT)Q(N) on phospholipid membranes.

Misfolding of amyloidogenic proteins and their interactions with membranes

In this paper, we discuss amyloidogenic proteins, their misfolding, resulting structures, and interactions with membranes, which lead to membrane damage and subsequent cell death. Many of these proteins are implicated in serious illnesses such as Alzheimer’s disease and Parkinson’s disease. Misfolding of amyloidogenic proteins leads to the formation of polymorphic oligomers and fibrils. Oligomeric aggregates are widely thought to be the toxic species, however, fibrils also play a role in membrane damage. We focus on the structure of these aggregates and their interactions with model membranes. Study of interactions of amlyoidogenic proteins with model and natural membranes has shown the importance of the lipid bilayer in protein misfolding and aggregation and has led to the development of several models for membrane permeabilization by the resulting amyloid aggregates. We discuss several of these models: formation of structured pores by misfolded amyloidogenic proteins, extraction of lipids, interactions with receptors in biological membranes, and membrane destabilization by amyloid aggregates perhaps analogous to that caused by antimicrobial peptides.

Membrane binding and insertion of a pHLIP peptide studied by all-atom molecular dynamics simulations

Recent experiments in function mechanism study reported that a pH low-insertion peptide (pHLIP) can insert into a zwitterionic palmitoyloleoylphosphatidylcholine (POPC) lipid bilayer at acidic pH while binding to the bilayer surface at basic pH. However, the atomic details of the pH-dependent interaction of pHLIP with a POPC bilayer are not well understood. In this study, we investigate the detailed interactions of pHLIP with a POPC bilayer at acidic and basic pH conditions as those used in function mechanism study, using all-atom molecular dynamics (MD) simulations. Simulations have been performed by employing the initial configurations, where pHLIP is placed in aqueous solution, parallel to bilayer surface (system S), partially-inserted (system P), or fully-inserted (system F) in POPC bilayers. On the basis of multiple 200-ns MD simulations, we found (1) pHLIP in system S can spontaneously insert into a POPC bilayer at acidic pH, while binding to the membrane surface at basic pH; (2) pHLIP in system P can insert deep into a POPC bilayer at acidic pH, while it has a tendency to exit, and stays at bilayer surface at basic pH; (3) pHLIP in system F keeps in an α-helical structure at acidic pH while partially unfolding at basic pH. This study provides at atomic-level the pH-induced insertion of pHLIP into POPC bilayer.