Stable and metastable states of human amylin in solution

Sequence-based identification of amyloidogenic β-hairpins reveals a prostatic acid phosphatase fragment promoting semen amyloid formation

β-Structure-rich amyloid fibrils are hallmarks of several diseases, including Alzheimer's (AD), Parkinson's (PD), and type 2 diabetes (T2D). While amyloid fibrils typically consist of parallel β-sheets, the anti-parallel β-hairpin is a structural motif accessible to amyloidogenic proteins in their monomeric and oligomeric states. Here, to investigate implications of β-hairpins in amyloid formation, potential β-hairpin-forming amyloidogenic segments in the human proteome were predicted based on sequence similarity with β-hairpins previously observed in Aβ, α-synuclein, and islet amyloid polypeptide, amyloidogenic proteins associated with AD, PD, and T2D, respectively. These three β-hairpins, established upon binding to the engineered binding protein β-wrapin AS10, are characterized by proximity of two sequence segments rich in hydrophobic and aromatic amino acids, with high β-aggregation scores according to the TANGO algorithm. Using these criteria, 2505 potential β-hairpin-forming amyloidogenic segments in 2098 human proteins were identified. Characterization of a test set of eight protein segments showed that seven assembled into Thioflavin T-positive aggregates and four formed β-hairpins in complex with AS10 according to NMR. One of those is a segment of prostatic acid phosphatase (PAP) comprising amino acids 185-208. PAP is naturally cleaved into fragments, including PAP(248-286) which forms functional amyloid in semen. We find that PAP(185-208) strongly decreases the protein concentrations required for fibril formation of PAP(248-286) and of another semen amyloid peptide, SEM1(86-107), indicating that it promotes nucleation of semen amyloids. In conclusion, β-hairpin-forming amyloidogenic protein segments could be identified in the human proteome with potential roles in functional or disease-related amyloid formation.

Graphitic nanoflakes modulate the structure and binding of human amylin

Human amylin is an inherently disordered protein whose ability to form amyloid fibrils is linked to the onset of type II diabetes. Graphitic nanomaterials have potential in managing amyloid diseases as they can disrupt protein aggregation processes in biological settings, but optimising these materials to prevent fibrillation is challenging. Here, we employ bias-exchange molecular dynamics simulations to systematically study the structure and adsorption preferences of amylin on graphitic nanoflakes that vary in their physical dimensions and surface functionalisation. Our findings reveal that nanoflake size and surface oxidation both influence the structure and adsorption preferences of amylin. The purely hydrophobic substrate of pristine graphene (PG) nanoflakes encourages non-specific protein adsorption, leading to unrestricted lateral mobility once amylin adheres to the surface. Particularly on larger PG nanoflakes, this induces structural changes in amylin that may promote fibril formation, such as the loss of native helical content and an increase in β-sheet character. In contrast, oxidised graphene nanoflakes form hydrogen bonds between surface oxygen sites and amylin, and as such restricting protein mobility. Reduced graphene oxide (rGO) flakes, featuring lower amounts of surface oxidation, are amphiphilic and exhibit substantial regions of bare carbon which promote protein binding and reduced conformational flexibility, leading to conservation of the native structure of amylin. In comparison, graphene oxide (GO) nanoflakes, which are predominantly hydrophilic and have a high degree of surface oxidation, facilitate considerable protein structural variability, resulting in substantial contact area between the protein and GO, and subsequent protein unfolding. Our results indicate that tailoring the size, oxygen concentration and surface patterning of graphitic nanoflakes can lead to specific and robust protein binding, ultimately influencing the likelihood of fibril formation. These atomistic insights provide key design considerations for the development of graphitic nanoflakes that can modulate protein aggregation by sequestering protein monomers in the biological environment and inhibit conformational changes linked to amyloid fibril formation.

Exploring the central region of amylin and its analogs aggregation: the influence of metal ions and residue substitutions

Human amylin (hIAPP) is found in the form of amyloid deposits within the pancreatic cells of nearly all patients diagnosed with type 2 diabetes mellitus (T2DM). However, rat amylin (rIAPP) and pramlintide - hIAPP analogs - are both non-toxic and non-amyloidogenic. Their primary sequences exhibit only slight variations in a few amino acid residues, primarily concentrated in the central region, spanning residues 20 to 29. This inspired us to study this fragment and investigate the impact on the aggregation properties of substituting residues within the central region of amylin and its analogs. Six fragments derived from amylin have undergone comprehensive testing against various metal ions by implementing a range of analytical techniques, including Nuclear Magnetic Resonance (NMR) spectroscopy, Thioflavin T (ThT) assays, Atomic Force Microscopy (AFM), and cytotoxicity assays. These methodologies serve to provide a thorough understanding of how the substitutions and interactions with metal ions impact the aggregation behavior of amylin and its analogs.

The Inhibition Effect of Epigallocatechin-3-Gallate on the Co-Aggregation of Amyloid-β and Human Islet Amyloid Polypeptide Revealed by Replica Exchange Molecular Dynamics Simulations

Alzheimer's disease and Type 2 diabetes are two epidemiologically linked diseases which are closely associated with the misfolding and aggregation of amyloid proteins amyloid-β (Aβ) and human islet amyloid polypeptide (hIAPP), respectively. The co-aggregation of the two amyloid proteins is regarded as the fundamental molecular mechanism underlying their pathological association. The green tea extract epigallocatechin-3-gallate (EGCG) has been extensively demonstrated to inhibit the amyloid aggregation of Aβ and hIAPP proteins. However, its potential role in amyloid co-aggregation has not been thoroughly investigated. In this study, we employed the enhanced-sampling replica exchange molecular dynamics simulation (REMD) method to investigate the effect of EGCG on the co-aggregation of Aβ and hIAPP. We found that EGCG molecules substantially diminish the β-sheet structures within the amyloid core regions of Aβ and hIAPP in their co-aggregates. Through hydrogen-bond, π-π and cation-π interactions targeting polar and aromatic residues of Aβ and hIAPP, EGCG effectively attenuates both inter-chain and intra-chain interactions within the co-aggregates. All these findings indicated that EGCG can effectively inhibit the co-aggregation of Aβ and hIAPP. Our study expands the potential applications of EGCG as an anti-amyloidosis agent and provides therapeutic options for the pathological association of amyloid misfolding disorders.

Nucleation of α-Synuclein Amyloid Fibrils Induced by Cross-Interaction with β-Hairpin Peptides Derived from Immunoglobulin Light Chains

Heterologous interactions between different amyloid-forming proteins, also called cross-interactions, may have a critical impact on disease-related amyloid formation. β-hairpin conformers of amyloid-forming proteins have been shown to affect homologous interactions in the amyloid self-assembly process. Here, we applied two β-hairpin-forming peptides derived from immunoglobulin light chains as models to test how heterologous β-hairpins modulate the fibril formation of Parkinson's disease-associated protein α-synuclein (αSyn). The peptides SMAhp and LENhp comprise β-strands C and C' of the κ4 antibodies SMA and LEN, which are associated with light chain amyloidosis and multiple myeloma, respectively. SMAhp and LENhp bind with high affinity to the β-hairpin-binding protein β-wrapin AS10 according to isothermal titration calorimetry and NMR spectroscopy. The addition of SMAhp and LENhp affects the kinetics of αSyn aggregation monitored by Thioflavin T (ThT) fluorescence, with the effect depending on assay conditions, salt concentration, and the applied β-hairpin peptide. In the absence of agitation, substoichiometric concentrations of the hairpin peptides strongly reduce the lag time of αSyn aggregation, suggesting that they support the nucleation of αSyn amyloid fibrils. The effect is also observed for the aggregation of αSyn fragments lacking the N-terminus or the C-terminus, indicating that the promotion of nucleation involves the interaction of hairpin peptides with the hydrophobic non-amyloid-β component (NAC) region.

Molecular Insights into the Self-Assembly of a Full-Length hIAPP Trimer: β-Protofibril Formed by β-Hairpin Lateral or Longitudinal Association

The fibrillar protein deposits of the human islet amyloid polypeptide (hIAPP) in the pancreatic islet of Langerhans are pathological hallmark of type II diabetes. Extensive experimental studies have revealed that the oligomeric formations of the hIAPP are more toxic than the mature fibrils. Exploring the oligomeric conformations in the early aggregation state is valuable for effective therapeutics. In this work, using the all-atom explicit-solvent replica exchange molecular dynamic (REMD) simulations, we investigated the structural features and the assembly mechanisms of the full-length hIAPP trimer in solution. The hIAPP trimer adopted more β-sheets than a-helix conformations, and three types of ordered conformations including open β-barrel, single-layer, and double-layer U-shaped β-sheet structures with five β-strands were captured in our simulations. A representative single-layer β-sheet conformation with a CCS value of 1400 Ų in our simulations matches exactly the experimentally ESI-IMS-MS-derived hIAPP trimer sample. These five β-strand conformations formed via the β-hairpin lateral and longitudinal association, respectively, showing two β-protofibril formation models. To the best of our knowledge, it is the first time to reveal two routes to β-sheet formation in the hIAPP trimers on the atomic level. The contact probabilities between pairs of the β-stranded residue show that the hydrophobic interactions between the residues F₁₅ ∼ V₁₇ and A₂₅ ∼ L₂₇ are responsible for the inter- and intra-peptide β-hairpin formations. All of these results indicate that the β-sheet formation is the first step in the conformational changes toward pathological aggregation and provides evidence of the β-sheet assembly mechanism into hIAPP aggregation.

Evolution of large Aβ16-22 aggregates at atomic details and potential of mean force associated to peptide unbinding and fragmentation events

Atomic characterization of large nonfibrillar aggregates of amyloid polypeptides cannot be determined by experimental means. Starting from β-rich aggregates of Y and elongated topologies predicted by coarse-grained simulations and consisting of more than 100 Aβ16-22 peptides, we performed atomistic molecular dynamics (MD), replica exchange with solute scaling (REST2), and umbrella sampling simulations using the CHARMM36m force field in explicit solvent. Here, we explored the dynamics within 3 μs, the free energy landscape, and the potential of mean force associated with either the unbinding of one single peptide in different configurations within the aggregate or fragmentation events of a large number of peptides. Within the time scale of MD and REST2, we find that the aggregates experience slow global conformational plasticity, and remain essentially random coil though we observe slow beta-strand structuring with a dominance of antiparallel beta-sheets over parallel beta-sheets. Enhanced REST2 simulation is able to capture fragmentation events, and the free energy of fragmentation of a large block of peptides is found to be similar to the free energy associated with fibril depolymerization by one chain for longer Aβ sequences.

Influence of force field choice on the conformational landscape of rat and human islet amyloid polypeptide

Human islet amyloid polypeptide (hIAPP) is a naturally occurring, intrinsically disordered protein (IDP) whose abnormal aggregation into toxic soluble oligomers and insoluble amyloid fibrils is a pathological feature in type-2 diabetes. Rat IAPP (rIAPP) differs from hIAPP by only six amino acids yet has a reduced tendency to aggregate or form fibrils. The structures of the monomeric forms of IAPP are difficult to characterize due to their intrinsically disordered nature. Molecular dynamics simulations can provide a detailed characterization of the monomeric forms of rIAPP and hIAPP in near-physiological conditions. In this work, the conformational landscapes of rIAPP and hIAPP as a function of secondary structure content were predicted using well-tempered bias exchange metadynamics simulations. Several combinations of commonly used biomolecular force fields and water models were tested. The predicted conformational preferences of both rIAPP and hIAPP are typical of IDPs, exhibiting dominant random coil structures but showing a low propensity for transient α-helical conformations. Predicted nuclear magnetic resonance Cα chemical shifts reveal different preferences with each force field towards certain conformations, with AMBERff99SBnmr2/TIP4Pd showing the best agreement with the experiment. Comparisons of secondary structure content demonstrate residue-specific differences between hIAPP and rIAPP that may reflect their different aggregation propensities.

Understanding the mechanism of amylin aggregation: From identifying crucial segments to tracing dominant sequential events to modeling potential aggregation suppressors

One of the most abundant, prevailing, and life-threatening human diseases that are currently baffling the scientific community is type 2 diabetes (T2D). The self-association of human amylin has been implicated in the pathogenesis of T2D, though with an inconclusive understanding of the mechanism. Hence, we focused on the characterization of the conformational ensembles of all the species that are believed to define the structural polymorphism of the aggregation process - the functional monomeric, the initially self-associated oligomeric, and the structured protofibril - by employing near-equilibrium, non-equilibrium, and equilibrium atomistic simulations on the sporadic, two familial variants (S20G and G33R), and their proline-substituted forms (S20P and G33P). The dynamic near-equilibrium assays hint toward - the abundance of helical conformation in the monomeric state, the retainment of the helicity in the initial self-associated oligomeric phase pointing toward the existence of the helix-helix association mechanism, the difference in preference of specific segments to have definite secondary structural features, the phase-dependent variability in the dominance of specific segments and mutation sites, and the simultaneous presence of generic and unique features among various sequences. Furthermore, the non-equilibrium pulling assays exemplify a generic sequential unzipping mechanism of the protofibrils, however, the sequence-dependent uniqueness comes from the difference in location and magnitude of the control of a specific terminus. Importantly, the equilibrium thermodynamic assays efficiently rank order the potential of aggregability among sequences and consequently suggests the probability of designing effective aggregation suppressors against sporadic and familial amylin variants incorporating proline as the mutation.

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.

Molecular Dynamics Simulations Indicate Aromaticity as a Key Factor in the Inhibition of IAPP(20-29) Aggregation

Islet amyloid polypeptide (IAPP) is a 37-residue amyloidogenic hormone implicated in the progression of Type II Diabetes (T2D). T2D affects an estimated 422 million people yearly and is a comorbidity with numerous diseases. IAPP forms toxic oligomers and amyloid fibrils that reduce pancreatic β-cell mass and exacerbate the T2D disease state. Toxic oligomer formation is attributed, in part, to the formation of interpeptide β-strands comprised of residues 20-29 (IAPP_(20-29)). Flavonoids, a class of polyphenolic natural products, have been found experimentally to inhibit IAPP aggregate formation. Many of these small flavonoids differ structurally only slightly; the influence of functional group placement on inhibiting the aggregation of the IAPP_(20-29) has yet to be explored. To probe the role of small-molecule structural features that impede IAPP aggregation, molecular dynamics simulations were performed to observe trimer formation on a model fragment of IAPP_(20-29) in the presence of morin, quercetin, dihydroquercetin, epicatechin, and myricetin. Contacts between Phe23 residues were critical to oligomer formation, and small-molecule contacts with Phe23 were a key predictor of β-strand reduction. Structural properties influencing the ability of compounds to disrupt Phe23-Phe23 contacts included aromaticity and carbonyl and hydroxyl group placement. This work provides key information on design considerations for T2D therapeutics that target IAPP aggregation.

Molecular Dynamics Simulation Studies on the Aggregation of Amyloid-β Peptides and Their Disaggregation by Ultrasonic Wave and Infrared Laser Irradiation

Alzheimer’s disease is understood to be caused by amyloid fibrils and oligomers formed by aggregated amyloid-β (Aβ) peptides. This review article presents molecular dynamics (MD) simulation studies of Aβ peptides and Aβ fragments on their aggregation, aggregation inhibition, amyloid fibril conformations in equilibrium, and disruption of the amyloid fibril by ultrasonic wave and infrared laser irradiation. In the aggregation of Aβ, a β-hairpin structure promotes the formation of intermolecular β-sheet structures. Aβ peptides tend to exist at hydrophilic/hydrophobic interfaces and form more β-hairpin structures than in bulk water. These facts are the reasons why the aggregation is accelerated at the interface. We also explain how polyphenols, which are attracting attention as aggregation inhibitors of Aβ peptides, interact with Aβ. An MD simulation study of the Aβ amyloid fibrils in equilibrium is also presented: the Aβ amyloid fibril has a different structure at one end from that at the other end. The amyloid fibrils can be destroyed by ultrasonic wave and infrared laser irradiation. The molecular mechanisms of these amyloid fibril disruptions are also explained, particularly focusing on the function of water molecules. Finally, we discuss the prospects for developing treatments for Alzheimer’s disease using MD simulations.

Sequence-dependent aggregation-prone conformations of islet amyloid polypeptide

Amyloid proteins, which aggregate to form highly ordered structures, play a crucial role in various disease pathologies. Despite many previous studies on amyloid fibrils, which are an end product of protein aggregation, the structural characteristics of amyloid proteins in the early stage of aggregation and their related aggregation mechanism still remain elusive. The role of the amino acid sequence in the aggregation-prone structures of amyloid proteins at such a stage is not understood. Here, we have studied the sequence-dependent structural characteristics of islet amyloid polypeptide based on atomistic simulations and spectroscopic experiments. We show that the amino acid sequence determines non-bonded interactions that play a leading role in the formation of aggregation-prone conformations. Specifically, a single point mutation critically changes the population of aggregation-prone conformations, resulting in a change of the aggregation mechanism. Our simulation results were supported by experimental results suggesting that mutation affects the kinetics of aggregation and the structural characteristics of amyloid aggregates. Our study provides an insight into the role of sequence-dependent aggregation-prone conformations in the underlying mechanisms of amyloid aggregation.

Distinct Modes of Action of IAPP Oligomers on Membranes

Islet amyloid polypeptide (IAPP, also known as amylin) is a peptide hormone that is co-secreted with insulin by pancreatic β-cells and forms amyloid aggregates in type II diabetes. Various lines of evidence indicate that oligomers of this peptide may induce toxicity by disrupting or forming pores in cell membranes, but the structure of these pores is unknown. Here, we create models of pores for both helical and β-structured peptides using implicit membrane modeling and test their stability using multimicrosecond all-atom simulations. We find that the helical peptides behave similarly to antimicrobial peptides; they remain stably inserted in a highly tilted or partially unfolded configuration creating a narrow water channel. Parallel helix orientation creates a somewhat larger pore. An octameric β barrel of parallel β-hairpins is highly stable in the membrane, whereas the corresponding barrel made of antiparallel hairpins is not. We propose that certain experiments probe the helical pore state while others probe the β-structured pore state; this provides a possible explanation for lack of correlation that is sometimes observed between in vivo toxicity and in vitro liposome permeabilization experiments.

Helix-coil transition and conformational deformity in Aβ42-monomer: a case study using the Zn2+ cation

The metal ions (like Fe2+, Zn²⁺, Cu2+) are known to influence the amyloid beta (Aβ) aggregation. In this study, we have examined the conformational and dynamical changes during the coordination of Aβ-monomer with the Zn2+ ion using all-atom molecular dynamics (MD) simulations using explicit solvent models. We have probed the unfolding of the full-length Aβ₄₂ monomer both inclusive and exclusive of the Zn2+ cation, with 1:1 ratio of the peptide and the Zn2+ cation. The inclusion of the Zn2+ cation shows differential intra-peptide interactions which has been probed using various analyses. The Helix - Coil transition of the wild type Aβ42 monomer is studied using the steered molecular dynamics simulations by taking the end-to-end C-α distance across the peptide. This gives an idea of the unequal intra - peptide and peptide - water interactions being found across the length of the Aβ monomer. The transition of an α-helix dominated wild-type (WT) Aβ structure to the unfolded coil structure gives significant evidence of the intra-peptide hydrogen bonding shifts in the presence of the Zn2+ cation. This accounts for the structural and the dynamical variations that take place in the Aβ monomer in the presence of the Zn2+ cation to mimic the conditions/environment at the onset of fibrillation.Communicated by Ramaswamy H. Sarma.

Role of Water Molecules and Helix Structure Stabilization in the Laser-Induced Disruption of Amyloid Fibrils Observed by Nonequilibrium Molecular Dynamics Simulations

Water plays a crucial role in the formation and destruction of biomolecular structures. The mechanism for destroying biomolecular structures was thought to be an active breaking of hydrogen bonds by water molecules. However, using nonequilibrium molecular dynamics simulations, in which an amyloid-β amyloid fibril was destroyed via infrared free-electron laser (IR-FEL) irradiation, we discovered a new mechanism, in which water molecules disrupt protein aggregates. The intermolecular hydrogen bonds formed by C═O and N-H in the fibril are broken at each pulse of laser irradiation. These bonds spontaneously re-form after the irradiation in many cases. However, when a water molecule happens to enter the gap between C═O and N-H, it inhibits the re-formation of the hydrogen bonds. Such sites become defects in the regularly aligned hydrogen bonds, from which all hydrogen bonds in the intermolecular β-sheet are broken as the fraying spreads. This role of water molecules is entirely different from other known mechanisms. This new mechanism can explain the recent experiments showing that the amyloid fibrils are not destroyed by laser irradiation under dry conditions. Additionally, we found that helix structures form more after the amyloid disruption; this is because the resonance frequency is different in a helix structure. Our findings provide a theoretical basis for the application of IR-FEL to the future treatment of amyloidosis.

Promotion and Inhibition of Amyloid-β Peptide Aggregation: Molecular Dynamics Studies

Aggregates of amyloid-β (Aβ) peptides are known to be related to Alzheimer’s disease. Their aggregation is enhanced at hydrophilic–hydrophobic interfaces, such as a cell membrane surface and air-water interface, and is inhibited by polyphenols, such as myricetin and rosmarinic acid. We review molecular dynamics (MD) simulation approaches of a full-length Aβ peptide, Aβ40, and Aβ(16–22) fragments in these environments. Since these peptides have both hydrophilic and hydrophobic amino acid residues, they tend to exist at the interfaces. The high concentration of the peptides accelerates the aggregation there. In addition, Aβ40 forms a β-hairpin structure, and this structure accelerates the aggregation. We also describe the inhibition mechanism of the Aβ(16–22) aggregation by polyphenols. The aggregation of Aβ(16–22) fragments is caused mainly by the electrostatic attraction between charged amino acid residues known as Lys16 and Glu22. Since polyphenols form hydrogen bonds between their hydroxy and carboxyl groups and these charged amino acid residues, they inhibit the aggregation.

Effect of surface chemistry on islet amyloid polypeptide conformation

The formation of dense, linear arrays (fibrils) by biomolecules is the hallmark of a number of degenerative diseases, such as Alzheimer's and type-2 diabetes. Protein fibrils have also attracted interest as building blocks for new materials. It has long been recognized that surfaces can affect the fibrillation process. Recent work on the model fibril forming protein human islet amyloid polypeptide (hIAPP) has shown that while the protein concentration is highest at hydrophobic surfaces, the rate of fibril formation is lower than on other surfaces. To understand this, replica exchange molecular dynamics simulations were used to investigate the conformations that hIAPP adopts on surfaces of different hydrophobicities. The hydrophobic surface stabilizes α-helical structures which are significantly different to those found on the hydrophilic surface and in bulk solution. There is also a greatly reduced conformational ensemble on the hydrophobic surface due to long-lived contacts between hydrophobic residues on the protein and the surface. This new microscopic information will help us determine the mechanism of the enhancement of fibril formation on surfaces and provides new insight into the effect of nanointerfaces and protein conformation.

The metastable states of proteins

The intriguing process of protein folding comprises discrete steps that stabilize the protein molecules in different conformations. The metastable state of protein is represented by specific conformational characteristics, which place the protein in a local free energy minimum state of the energy landscape. The native-to-metastable structural transitions are governed by transient or long-lived thermodynamic and kinetic fluctuations of the intrinsic interactions of the protein molecules. Depiction of the structural and functional properties of metastable proteins is not only required to understand the complexity of folding patterns but also to comprehend the mechanisms of anomalous aggregation of different proteins. In this article, we review the properties of metastable proteins in context of their stability and capability of undergoing atypical aggregation in physiological conditions.

Symmetry-breaking transitions in the early steps of protein self-assembly

Protein misfolding and subsequent self-association are complex, intertwined processes, resulting in development of a heterogeneous population of aggregates closely related to many chronic pathological conditions including Type 2 Diabetes Mellitus and Alzheimer's disease. To address this issue, here, we develop a theoretical model in the general framework of linear stability analysis. According to this model, self-assemblies of peptides with pronounced conformational flexibility may become, under particular conditions, unstable and spontaneously evolve toward an alternating array of partially ordered and disordered monomers. The predictions of the theory were verified by atomistic molecular dynamics (MD) simulations of islet amyloid polypeptide (IAPP) used as a paradigm of aggregation-prone polypeptides (proteins). Simulations of dimeric, tetrameric, and hexameric human-IAPP self-assemblies at physiological electrolyte concentration reveal an alternating distribution of the smallest domains (of the order of the peptide mean length) formed by partially ordered (mainly β-strands) and disordered (turns and coil) arrays. Periodicity disappears upon weakening of the inter-peptide binding, a result in line with the predictions of the theory. To further probe the general validity of our hypothesis, we extended the simulations to other peptides, the Aβ(1-40) amyloid peptide, and the ovine prion peptide as well as to other proteins (SOD1 dimer) that do not belong to the broad class of intrinsically disordered proteins. In all cases, the oligomeric aggregates show an alternate distribution of partially ordered and disordered monomers. We also carried out Surface Enhanced Raman Scattering (SERS) measurements of hIAPP as an experimental validation of both the theory and in silico simulations.

Theoretical and computational advances in protein misfolding

Misfolded proteins escape the cellular quality control mechanism and fail to fold properly or remain correctly folded leading to a loss in their functional specificity. Thus misfolding of proteins cause a large number of very different diseases ranging from errors in metabolism to various types of complex neurodegenerative diseases. A theoretical and computational perspective of protein misfolding is presented with a special emphasis on its salient features, mechanism and consequences. These insights quantitatively analyze different determinants of misfolding, that may be applied to design disease specific molecular targets.

Molecular Dynamics Simulations Reveal the Inhibitory Mechanism of Dopamine against Human Islet Amyloid Polypeptide (hIAPP) Aggregation and Its Destabilization Effect on hIAPP Protofibrils

The aberrant self-assembly of human islet amyloid polypeptide (hIAPP) into toxic oligomers, protofibrils, and mature fibrils is associated with the pathogenesis of type 2 diabetes (T2D). Inhibition of hIAPP aggregation and destabilization of preformed hIAPP fibrils are considered as two major therapeutic strategies for treating T2D. Previous experimental studies reported that dopamine prevented the formation of hIAPP oligomers and fibrils. However, the underlying inhibitory mechanism at the atomic level remains elusive. Herein we investigated the conformational ensembles of hIAPP dimer with and without dopamine using replica-exchange molecular dynamics simulations. The simulations demonstrated that dopamine preferentially bound to R11, L12, F15, H18, F23, I26, L27, and Y37 residues, inhibited the formation of β-sheets in the amyloidogenic regions spanning residues ¹¹RLANFLVH¹⁸, ²²NFGAIL²⁷, and ³⁰TNVGSNT³⁶, and resulted in more disordered hIAPP dimers, thus hindering the amyloid formation of hIAPP. Protonated and deprotonated dopamine molecules displayed distinct binding capabilities but bound to similar residue sites on hIAPP. Additional microsecond molecular dynamics simulations showed that dopamine mainly bound to the β1 and turn regions of hIAPP protofibril and destabilized the protofibril structure. This study not only revealed the molecular mechanism of dopamine toward the inhibition of hIAPP aggregation but also demonstrated the protofibril-destabilizing effects of dopamine, which may be helpful for the design of drug candidates to treat T2D.

Modulation of Disordered Proteins with a Focus on Neurodegenerative Diseases and Other Pathologies

Intrinsically disordered proteins (IDPs) do not have rigid 3D structures, showing changes in their folding depending on the environment or ligands. Intrinsically disordered proteins are widely spread in eukaryotic genomes, and these proteins participate in many cell regulatory metabolism processes. Some IDPs, when aberrantly folded, can be the cause of some diseases such as Alzheimer′s, Parkinson′s, and prionic, among others. In these diseases, there are modifications in parts of the protein or in its entirety. A common conformational variation of these IDPs is misfolding and aggregation, forming, for instance, neurotoxic amyloid plaques. In this review, we discuss some IDPs that are involved in neurodegenerative diseases (such as beta amyloid, alpha synuclein, tau, and the “IDP-like” PrP), cancer (p53, c-Myc), and diabetes (amylin), focusing on the structural changes of these IDPs that are linked to such pathologies. We also present the IDP modulation mechanisms that can be explored in new strategies for drug design. Lastly, we show some candidate drugs that can be used in the future for the treatment of diseases caused by misfolded IDPs, considering that cancer therapy has more advanced research in comparison to other diseases, while also discussing recent and future developments in this area of research. Therefore, we aim to provide support to the study of IDPs and their modulation mechanisms as promising approaches to combat such severe diseases.

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.

Deformation of stable and toxic hIAPP oligomers by liposomes with distinct nanomechanical features and reduced cytotoxicity

Nanoliposomes can induce hIAPP oligomers to undergo fibrillation with distinct mechanical properties and reduced cytotoxicity.

Nucleation of β-rich oligomers and β-barrels in the early aggregation of human islet amyloid polypeptide

The self-assembly of human islet amyloid polypeptide (hIAPP) into β-sheet rich amyloid aggregates is associated with pancreatic β-cell death in type 2 diabetes (T2D). Prior experimental studies of hIAPP aggregation reported the early accumulation of α-helical intermediates before the rapid conversion into β-sheet rich amyloid fibrils, as also corroborated by our experimental characterizations with transmission electron microscopy and Fourier transform infrared spectroscopy. Although increasing evidence suggests that small oligomers populating early hIAPP aggregation play crucial roles in cytotoxicity, structures of these oligomer intermediates and their conformational conversions remain unknown, hindering our understanding of T2D disease mechanism and therapeutic design targeting these early aggregation species. We further applied large-scale discrete molecule dynamics simulations to investigate the oligomerization of full-length hIAPP, employing multiple molecular systems of increasing number of peptides. We found that the oligomerization process was dynamic, involving frequent inter-oligomeric exchanges. On average, oligomers had more α-helices than β-sheets, consistent with ensemble-based experimental measurements. However, in ~4-6% independent simulations, β-rich oligomers expected as the fibrillization intermediates were observed, especially in the pentamer and hexamer simulations. These β-rich oligomers could adopt β-barrel conformations, recently postulated to be the toxic oligomer species but only observed computationally in the aggregates of short amyloid protein fragments. Free-energy analysis revealed high energies of these β-rich oligomers, supporting the nucleated conformational changes of oligomers in amyloid aggregation. β-barrel oligomers of full-length hIAPP with well-defined three-dimensional structures may play an important pathological role in T2D etiology and may be a therapeutic target for the disease.

Characterisation of the Structure and Oligomerisation of Islet Amyloid Polypeptides (IAPP): A Review of Molecular Dynamics Simulation Studies

Human islet amyloid polypeptide (hIAPP) is a naturally occurring, intrinsically disordered protein whose abnormal aggregation into amyloid fibrils is a pathological feature in type 2 diabetes, and its cross-aggregation with amyloid beta has been linked to an increased risk of Alzheimer's disease. The soluble, oligomeric forms of hIAPP are the most toxic to β-cells in the pancreas. However, the structure of these oligomeric forms is difficult to characterise because of their intrinsic disorder and their tendency to rapidly aggregate into insoluble fibrils. Experimental studies of hIAPP have generally used non-physiological conditions to prevent aggregation, and they have been unable to describe its soluble monomeric and oligomeric structure at physiological conditions. Molecular dynamics (MD) simulations offer an alternative for the detailed characterisation of the monomeric structure of hIAPP and its aggregation in aqueous solution. This paper reviews the knowledge that has been gained by the use of MD simulations, and its relationship to experimental data for both hIAPP and rat IAPP. In particular, the influence of the choice of force field and water models, the choice of initial structure, and the configurational sampling method used, are discussed in detail. Characterisation of the solution structure of hIAPP and its mechanism of oligomerisation is important to understanding its cellular toxicity and its role in disease states, and may ultimately offer new opportunities for therapeutic interventions.

Membranes as modulators of amyloid protein misfolding and target of toxicity

Abnormal protein aggregation is a hallmark of various human diseases. α-Synuclein, a protein implicated in Parkinson's disease, is found in aggregated form within Lewy bodies that are characteristically observed in the brains of PD patients. Similarly, deposits of aggregated human islet amyloid polypeptide (IAPP) are found in the pancreatic islets in individuals with type 2 diabetes mellitus. Significant number of studies have focused on how monomeric, disaggregated proteins transition into various amyloid structures leading to identification of a vast number of aggregation promoting molecules and processes over the years. Inasmuch as these factors likely enhance the formation of toxic, misfolded species, they might act as risk factors in disease. Cellular membranes, and particularly certain lipids, are considered to be among the major players for aggregation of α-synuclein and IAPP, and membranes might also be the target of toxicity. Past studies have utilized an array of biophysical tools, both in vitro and in vivo, to expound the membrane-mediated aggregation. Here, we focus on membrane interaction of α-synuclein and IAPP, and how various kinds of membranes catalyze or modulate the aggregation of these proteins and how, in turn, these proteins disrupt membrane integrity, both in vitro and in vivo. The membrane interaction and subsequent aggregation has been briefly contrasted to aggregation of α-synuclein and IAPP in solution. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy.

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.

Structural Modulation of Human Amylin Protofilaments by Naturally Occurring Mutations

Human islet amyloid polypeptide (hIAPP), also known as amylin, is a 37-amino-acid peptide, co-secreted with insulin, and widely found in fibril form in type-2 diabetes patients. By using all-atom molecular dynamics simulations, we study hIAPP fibril segments (i.e., fibrillar oligomers) formed with sequences of naturally occurring variants from cat, rat, and pig, presenting different aggregation propensities. We characterize the effect of mutations on the structural dynamics of solution-formed hIAPP fibril models built from solid-state NMR data. Results from this study are in agreement with experimental observations regarding their respective relative aggregation propensities. We analyze in detail the specific structural characteristics and infer mechanisms that modulate the conformational stability of amylin fibrils. Results provide a platform for further studies and the design of new drugs that could interfere with amylin aggregation and its cytotoxicity. One particular mutation, N31K, has fibril-destabilizing properties, and could potentially improve the solubility of therapeutic amylin analogs.

Site-specific detection of protein secondary structure using 2D IR dihedral indexing: a proposed assembly mechanism of oligomeric hIAPP

Human islet amyloid polypeptide (hIAPP) aggregates into fibrils through oligomers that have been postulated to contain α-helices as well as β-sheets. We employ a site-specific isotope labeling strategy that is capable of detecting changes in dihedral angles when used in conjunction with 2D IR spectroscopy. The method is analogous to the chemical shift index used in NMR spectroscopy for assigning protein secondary structure. We introduce isotope labels at two neighbouring residues, which results in an increased intensity and positive frequency shift if those residues are α-helical versus a negative frequency shift in β-sheets and turns. The 2D IR dihedral index approach is demonstrated for hIAPP in micelles for which the polypeptide structure is known, using pairs of 13C18O isotope labels L12A13 and L16V17, along with single labeled control experiments. Applying the approach to aggregation experiments performed in buffer, we show that about 27-38% of hIAPP peptides adopt an α-helix secondary structure in the monomeric state at L12A13, prior to aggregation, but not at L16V17 residues. At L16V17, the kinetics are described solely by the monomer and fiber conformations, but at L12A13 the kinetics exhibit a third state that is created by an oligomeric intermediate. Control experiments performed with a single isotope label at A13 exhibit two-state kinetics, indicating that a previously unknown change in dihedral angle occurs at L12A13 as hIAPP transitions from the intermediate to fiber structures. We propose a mechanism for aggregation, in which helices seed oligomer formation via structures analogous to leucine rich repeat proteins.

Effects of forcefield and sampling method in all-atom simulations of inherently disordered proteins: Application to conformational preferences of human amylin

Although several computational modelling studies have investigated the conformational behaviour of inherently disordered protein (IDP) amylin, discrepancies in identifying its preferred solution conformations still exist between various forcefields and sampling methods used. Human islet amyloid polypeptide has long been a subject of research, both experimentally and theoretically, as the aggregation of this protein is believed to be the lead cause of type-II diabetes. In this work, we present a systematic forcefield assessment using one of the most advanced non-biased sampling techniques, Replica Exchange with Solute Tempering (REST2), by comparing the secondary structure preferences of monomeric amylin in solution. This study also aims to determine the ability of common forcefields to sample a transition of the protein from a helical membrane bound conformation into the disordered solution state of amylin. Our results demonstrated that the CHARMM22* forcefield showed the best ability to sample multiple conformational states inherent for amylin. It is revealed that REST2 yielded results qualitatively consistent with experiments and in quantitative agreement with other sampling methods, however far more computationally efficiently and without any bias. Therefore, combining an unbiased sampling technique such as REST2 with a vigorous forcefield testing could be suggested as an important step in developing an efficient and robust strategy for simulating IDPs.

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.

Inhibitory Mechanism of Epigallocatechin Gallate on Fibrillation and Aggregation of Amidated Human Islet Amyloid Polypeptide

The abnormal fibrillation of human islet amyloid polypeptide (hIAPP) is associated with development of type II diabetes mellitus (T2DM). (-)-Epigallocatechin gallate (EGCG) can bind amyloid proteins to inhibit the fibrillation of these proteins. However, the mechanic detail of EGCG inhibiting amyloid formation is still unclear at the molecular level. In the present work, we sought to investigate the effect of EGCG on amidated hIAPP (hIAPP-NH₂ ) fibrillation and aggregation by using spectroscopic and microscopic techniques, and also sought to gain insights into the interaction of EGCG and hIAPP_22-27 by using spectroscopic experiments and quantum chemical calculations. ThT fluorescence, real-time NMR, and TEM studies demonstrated that EGCG inhibits the formation of hIAPP-NH₂ fibrils, while promoting the formation of hIAPP-NH₂ amorphous aggregates. Phenylalanine intrinsic fluorescence and NMR studies of the EGCG/hIAPP_22-27 complex revealed three important binding sites including the A ring of EGCG, residue Phe23, and residue Ile26. DFT calculations identified the dominant binding structures of EGCG/Phe23 and EGCG/Ile26 complexes, named structure I and structure II, respectively. Our study demonstrates the inhibitory mechanism of EGCG on fibrillation and aggregation of hIAPP-NH₂ in which EGCG interacts with hIAPP-NH₂ through hydrogen bonding and π-π interactions between the A ring and residue Phe23 as well as hydrophobic interactions between the A ring and residue Ile26, which can thus inhibit the interpeptide interaction between hIAPP-NH₂ monomers and finally inhibit fibrillation of hIAPP-NH₂ . This study agrees with and reinforces previous studies and offers an intuitive explanation at both the atomic and molecular levels. Our findings may provide an invaluable reference for the future development of new drugs in the management of diabetes.

Structural and thermodynamical properties of early human amylin oligomers using replica exchange molecular dynamics: mutation effect of three key residues F15, H18 and F23

A schematic representation of a possible oligomerization mechanism of hIAPP. β-Hairpins are proposed to self-assemble into early ordered oligomers by side-to-side association.

Effect of Post-Translational Amidation on Islet Amyloid Polypeptide Conformational Ensemble: Implications for Its Aggregation Early Steps

The human islet amyloid polypeptide (hIAPP) is an intrinsically disordered protein that can self-assemble into fibrillar aggregates that play a key role in the pathogenesis of the type II diabetes mellitus. hIAPP can transiently adopt α-helix and β-strand conformations that could be important intermediate species on the fibrillization pathway. However, experimental studies of the monomeric peptide conformations are limited due to its high aggregation propensity, and the early steps of the hIAPP association are not clearly characterized. In particular, the question of whether the aggregation-prone conformation is α-helical or β-strand-rich is still debated. In this study, combining extensive all-atom molecular dynamics (MD) and replica exchange molecular dynamics (REMD) simulations in explicit water, we shed some light on the differences between the amidated and non-amidated hIAPP conformational ensembles. Our study shows that, when compared to the amidated monomer, the non-amidation of hIAPP induces a significantly lower propensity to form β-strands, especially aggregation-prone β-hairpins. Since the fibrillization of the non-amidated hIAPP is significantly slower than that of the amidated peptide, this indicates that the early steps of the peptide oligomerization involve the association of β-hairpins or β-strands structures.