An atomistic characterization of high-density lipoproteins and the conserved "LN" region of apoA-I.

The physicochemical characteristics of the various subpopulations of high-density lipoproteins (HDLs) and, in particular, their surface properties determine their ability to scavenge lipids and interact with specific receptors and peptides. Five representative spheroidal HDL subpopulation models were mapped from a previously reported equilibrated coarse-grained (CG) description to an atomistic representation for subsequent molecular dynamics simulation. For each HDL model a range of finer-level analyses was undertaken, including the component-wise characterization of HDL surfaces, the average size and composition of hydrophobic surface patches, dynamic protein secondary structure monitoring, and the proclivity for solvent exposure of the proposed ß-amyloid (Aß) binding region of apolipoprotein A-I (apoA-I), "LN." This study reveals that previously characterized ellipsoidal HDL3a and HDL2a models revert to a more spherical geometry in an atomistic representation due to the enhanced conformational flexibility afforded to the apoA-I protein secondary structure, allowing for enhanced surface lipid packing and lower overall surface hydrophobicity. Indeed, the proportional surface hydrophobicity and apoA-I exposure reduced with increasing HDL size, consistent with previous characterizations. Furthermore, solvent exposure of the "LN" region of apoA-I was exclusively limited to the smallest HDL3c model within the timescale of the simulations, and typically corresponded to a distinct loss in secondary structure across the "LN" region to form part of a significant contiguous hydrophobic patch on the HDL surface. Taken together, these findings provide preliminary evidence for a subpopulation-specific interaction between HDL3c particles and circulating hydrophobic species such as Aß via the exposed "LN" region of apoA-I.

What determines sub-diffusive behavior in crowded protein solutions?

The aqueous environment inside cells is densely packed. A typical cell has a macromolecular concentration in the range 90-450 g/L, with 5%-40% of its volume being occupied by macromolecules, resulting in what is known as macromolecular crowding. The space available for the free diffusion of metabolites and other macromolecules is thus greatly reduced, leading to so-called excluded volume effects. The slow diffusion of macromolecules under crowded conditions has been explained using transient complex formation. However, sub-diffusion noted in earlier works is not well characterized, particularly the role played by transient complex formation and excluded volume effects. We have used Brownian dynamics simulations to characterize the diffusion of chymotrypsin inhibitor 2 in protein solutions of bovine serum albumin and lysozyme at concentrations ranging from 50 to 300 g/L. The predicted changes in diffusion coefficient as a function of crowder concentration are consistent with NMR experiments. The sub-diffusive behavior observed in the sub-microsecond timescale can be explained in terms of a so-called cage effect, arising from rattling motion in a local molecular cage as a consequence of excluded volume effects. By selectively manipulating the nature of interactions between protein molecules, we determined that excluded volume effects induce sub-diffusive dynamics at sub-microsecond timescales. These findings may help to explain the diffusion-mediated effects of protein crowding on cellular processes.

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.

The Molecular Docking of MAX Fungal Effectors with Plant HMA Domain-Binding Proteins.

Fungal effector proteins are important in mediating disease infections in agriculturally important crops. These secreted small proteins are known to interact with their respective host receptor binding partners in the host, either inside the cells or in the apoplastic space, depending on the localisation of the effector proteins. Consequently, it is important to understand the interactions between fungal effector proteins and their target host receptor binding partners, particularly since this can be used for the selection of potential plant resistance or susceptibility-related proteins that can be applied to the breeding of new cultivars with disease resistance. In this study, molecular docking simulations were used to characterise protein-protein interactions between effector and plant receptors. Benchmarking was undertaken using available experimental structures of effector-host receptor complexes to optimise simulation parameters, which were then used to predict the structures and mediating interactions of effector proteins with host receptor binding partners that have not yet been characterised experimentally. Rigid docking was applied for both the so-called bound and unbound docking of MAX effectors with plant HMA domain protein partners. All bound complexes used for benchmarking were correctly predicted, with 84% being ranked as the top docking pose using the ZDOCK scoring function. In the case of unbound complexes, a minimum of 95% of known residues were predicted to be part of the interacting interface on the host receptor binding partner, and at least 87% of known residues were predicted to be part of the interacting interface on the effector protein. Hydrophobic interactions were found to dominate the formation of effector-plant protein complexes. An optimised set of docking parameters based on the use of ZDOCK and ZRANK scoring functions were established to enable the prediction of near-native docking poses involving different binding interfaces on plant HMA domain proteins. Whilst this study was limited by the availability of the experimentally determined complexed structures of effectors and host receptor binding partners, we demonstrated the potential of molecular docking simulations to predict the likely interactions between effectors and their respective host receptor binding partners. This computational approach may accelerate the process of the discovery of putative interacting plant partners of effector proteins and contribute to effector-assisted marker discovery, thereby supporting the breeding of disease-resistant crops.

Ab Initio Modelling of the Structure of ToxA-like and MAX Fungal Effector Proteins.

Pathogenic fungal diseases in crops are mediated by the release of effector proteins that facilitate infection. Characterising the structure of these fungal effectors is vital to understanding their virulence mechanisms and interactions with their hosts, which is crucial in the breeding of plant cultivars for disease resistance. Several effectors have been identified and validated experimentally; however, their lack of sequence conservation often impedes the identification and prediction of their structure using sequence similarity approaches. Structural similarity has, nonetheless, been observed within fungal effector protein families, creating interest in validating the use of computational methods to predict their tertiary structure from their sequence. We used Rosetta ab initio modelling to predict the structures of members of the ToxA-like and MAX effector families for which experimental structures are known to validate this method. An optimised approach was then used to predict the structures of phenotypically validated effectors lacking known structures. Rosetta was found to successfully predict the structure of fungal effectors in the ToxA-like and MAX families, as well as phenotypically validated but structurally unconfirmed effector sequences. Interestingly, potential new effector structural families were identified on the basis of comparisons with structural homologues and the identification of associated protein domains.

Functional and structural consequences of TBK1 missense variants in frontotemporal lobar degeneration and amyotrophic lateral sclerosis.

Mutations in the Tank-binding kinase 1 (TBK1) gene were identified in 2015 in individuals with frontotemporal lobar degeneration (FTLD) and amyotrophic lateral sclerosis (ALS). They account for ∼3-4% of cases. To date, over 100 distinct mutations, including missense, nonsense, deletion, insertion, duplication, and splice-site mutations have been reported. While nonsense mutations are predicted to cause disease via haploinsufficiency, the mechanisms underlying disease pathogenesis with missense mutations is not fully elucidated. TBK1 is a kinase involved in neuroinflammation, which is commonly observed in these diseases. TBK1 also phosphorylates key autophagy mediators, thereby regulating proteostasis, a pathway that is dysregulated in ALS-FTLD. Recently, several groups have characterised various missense mutations with respect to their effects on the phosphorylation of known TBK1 substrates, TBK1 homodimerization, interaction with optineurin, and the regulation of autophagy and neuroinflammatory pathways. Further, the effects of either global or conditional heterozygous knock-out of Tbk1, or the heterozygous or homozygous knock-in of ALS-FTLD associated mutations, alone or when crossed with the SOD1G93A classical ALS mouse model or a TDP-43 mouse model, have been reported. In this review we summarise the known functional effects of TBK1 missense mutations. We also present novel modelling data that predicts the structural effects of missense mutations and discuss how they correlate with the known functional effects of these mutations.

Review: Investigating the aggregation of amyloid beta with surface plasmon resonance: Do different approaches yield different results?

Aggregation of amyloid beta into amyloid plaques in the brain is a hallmark characteristic of Alzheimer's disease. Therapeutics aimed at preventing or retarding amyloid formation often rely on detailed characterization of the underlying mechanism and kinetics of protein aggregation. Surface plasmon resonance (SPR) spectroscopy is a robust technique used to determine binding affinity and kinetics of biomolecular interactions. This approach has been used to characterize the mechanism of aggregation of amyloid beta but there are multiple pitfalls that need to be addressed when working with this and other amyloidogenic proteins. The choice of method for analyte preparation and ligand immobilization to a sensor chip can lead to different theoretical and practical implications in terms of the mathematical modelling of binding data, different mechanisms of binding and the presence of different interacting species. This review examines preparation methods for SPR characterisation of the aggregation of amyloid beta and their influence on the findings derived from such studies.

How does metabolic rate in plant shoot tips change after cryopreservation?

Cryopreservation allows the long-term storage of plant germplasm, but can cause damage to plant tissues, which must be repaired for survival to occur. This repair process is fuelled by the metabolic function of mitochondria; however, little is known about how metabolic function is affected by the cryopreservation process in plants. We compared metabolic rates of shoot tips of two Australian native species, Androcalva perlaria and Anigozanthos viridis. Overall, cryopreservation resulted in a significant reduction in the metabolic rates of shoot tips from both species, even in tissues that regenerated after cryopreservation. Metabolic rate did not increase within 48 h after of thawing, even in shoot tips which later regenerated. When examined in isolation, both pre-treatment on desiccation medium and exposure to cryoprotective agents significantly decreased metabolic rates in regenerating shoot tips of A. viridis, however both caused a significant increase in shoot tips of A. perlaria, suggesting diversity of response to cryopreservation stresses across species. Measurements of shoot tip metabolic rate during cryopreservation will inform investigations into cellular energy production and provide critical information on the state of shoot health after exposure to different cryoprotective treatments, which could play a useful role in guiding protocol optimisation for threatened species to maximise post-cryopreservation regeneration.

Characterisation of the Molecular Mechanism of Permeation of the Prodrug Me-5ALA across the Human Stratum Corneum Using Molecular Dynamics Simulations.

The barrier imposed by the outer layer of the skin, the stratum corneum, creates an almost impermeable environment for exogenous substances. Few lipophilic drugs with low molecular mass can passively diffuse through this layer, highlighting the need to develop methods to enable the delivery of more drugs via the transdermal route. The prodrug approach involves modifying the structure of a drug molecule to enhance its permeability across the skin, but it is often difficult to predict how exactly changes in chemical structure affect permeation. This study uses molecular dynamics simulations to predict permeability values and adequately characterise the molecular mechanism of permeation of the prodrugs Me-5ALA and its parent compound 5ALA across a molecular model of the lipid bilayers of the human stratum corneum. The influence of increased hydrophobicity in Me-5ALA on its permeation revealed a reduction in hydrogen bonding capability that enables it to interact more favourably with the hydrophobic region of the bilayer and diffuse at a faster rate with less resistance, thus making it a better permeant compared to its more hydrophilic parent compound. This molecular simulation approach offers a promising route for the rational design of drug molecules that can permeate effectively across the stratum corneum.

Pharmacological and structure-activity relationship studies of oleoyl-lysophosphatidylinositol synthetic mimetics.

Metabolic diseases, such as obesity and type 2 diabetes, are relentlessly spreading worldwide. The beginning of the 21st century has seen the introduction of mechanistically novel types of drugs, aimed primarily at keeping these pathologies under control. In particular, an important family of therapeutics exploits the beneficial physiology of the gut-derived glucagon-like peptide-1 (GLP-1), with important clinical benefits, from glycaemic control to cardioprotection. Nonetheless, these protein-based drugs act systemically as exogenous GLP-1 mimetics and are not exempt from side effects. The food-derived lipid oleoyl-lysophosphatidylinositol (LPI) is a potent GPR119-dependent GLP-1 secreting agent. Here we present a structure-activity relationship (SAR) study of a synthetic library of oleoyl-LPI mimetics capable to induce the physiological release of GLP-1 from gastrointestinal enteroendocrine cells (EECs). The best lead compounds have shown potent and efficient release of GLP-1 in vitro from human and murine cells, and in vivo in diabetic db/db mice. We have also generated a molecular model of oleoyl-LPI, as well as its best performing analogues, interacting with the orthosteric site of GPR119, laying foundational evidence for their pharmacological activity.

Does the SARS-CoV-2 Spike Protein Receptor Binding Domain Interact Effectively with the DPP4 (CD26) Receptor? A Molecular Docking Study.

ACE2 has been established as the main receptor for SARS-CoV-2. Since other human coronaviruses are known to use co-receptors for viral cell entry, it has been suggested that DPP4 (CD26) could be a potential additional binding target or co-receptor, supported by early molecular docking simulation studies. However, recent biophysical studies have shown this interaction to be very weak. We have conducted detailed molecular docking simulations to predict the potential binding interactions between the receptor binding domain (RBD) of the spike protein of SARS-CoV-2 and DPP4 and compare them with the interactions observed in the experimentally determined structure of the complex of MERS-CoV with DPP4. Whilst the overall binding mode of the RBD of SARS-CoV-2 to DPP4 is predicted to be similar to that observed in the MERS-CoV-DPP4 complex, including a number of equivalent interactions, important differences in the amino acid sequences of SARS-CoV-2 and MERS-CoV result in substantially weakened interactions with DPP4. This is shown to arise from differences in the predicted proximity, nature and secondary structure at the binding interface on the RBD of SARS-CoV-2. These findings do not support DPP4 being a significant receptor for SARS-CoV-2.

Design and Characterization of a Cell-Penetrating Peptide Derived from the SOX2 Transcription Factor.

SOX2 is an oncogenic transcription factor overexpressed in nearly half of the basal-like triple-negative breast cancers associated with very poor outcomes. Targeting and inhibiting SOX2 is clinically relevant as high SOX2 mRNA levels are positively correlated with decreased overall survival and progression-free survival in patients affected with breast cancer. Given its key role as a master regulator of cell proliferation, SOX2 represents an important scaffold for the engineering of dominant-negative synthetic DNA-binding domains (DBDs) that act by blocking or interfering with the oncogenic activity of the endogenous transcription factor in cancer cells. We have synthesized an interference peptide (iPep) encompassing a truncated 24 amino acid long C-terminus of SOX2 containing a potential SOX-specific nuclear localization sequence, and the determinants of the binding of SOX2 to the DNA and to its transcription factor binding partners. We found that the resulting peptide (SOX2-iPep) possessed intrinsic cell penetration and promising nuclear localization into breast cancer cells, and decreased cellular proliferation of SOX2 overexpressing cell lines. The novel SOX2-iPep was found to exhibit a random coil conformation predominantly in solution. Molecular dynamics simulations were used to characterize the interactions of both the SOX2 transcription factor and the SOX2-iPep with FGF4-enhancer DNA in the presence of the POU domain of the partner transcription factor OCT4. Predictions of the free energy of binding revealed that the iPep largely retained the binding affinity for DNA of parental SOX2. This work will enable the future engineering of novel dominant interference peptides to transport different therapeutic cargo molecules such as anti-cancer drugs into cells.

The unusual conformation of cross-strand disulfide bonds is critical to the stability of ß-hairpin peptides.

The cross-strand disulfides (CSDs) found in ß-hairpin antimicrobial peptides (ß-AMPs) show a unique disulfide geometry that is characterized by unusual torsion angles and a short Cα-Cα distance. While the sequence and disulfide bond connectivity of disulfide-rich peptides is well studied, much less is known about the disulfide geometry found in CSDs and their role in the stability of ß-AMPs. To address this, we solved the nuclear magnetic resonance (NMR) structure of the ß-AMP gomesin (Gm) at 278, 298, and 310 K, examined the disulfide bond geometry of over 800 disulfide-rich peptides, and carried out extensive molecular dynamics (MD) simulation of the peptides Gm and protegrin. The NMR data suggests Cα-Cα distances characteristic for CSDs are independent of temperature. Analysis of disulfide-rich peptides from the Protein Data Bank revealed that right-handed and left-handed rotamers are equally likely in CSDs. The previously reported preference for right-handed rotamers was likely biased by restricting the analysis to peptides and proteins solved using X-ray crystallography. Furthermore, data from MD simulations showed that the short Cα-Cα distance is critical for the stability of these peptides. The unique disulfide geometry of CSDs poses a challenge to biomolecular force fields and to retain the stability of ß-hairpin fold over long simulation times, restraints on the torsion angles might be required.

Molecular Dynamics Simulation of Small Molecules Interacting with Biological Membranes.

Cell membranes protect and compartmentalise cells and their organelles. The semi-permeable nature of these membranes controls the exchange of solutes across their structure. Characterising the interaction of small molecules with biological membranes is critical to understanding of physiological processes, drug action and permeation, and many biotechnological applications. This review provides an overview of how molecular simulations are used to study the interaction of small molecules with biological membranes, with a particular focus on the interactions of water, organic compounds, drugs and short peptides with models of plasma cell membrane and stratum corneum lipid bilayers. This review will not delve on other types of membranes which might have different composition and arrangement, such as thylakoid or mitochondrial membranes. The application of unbiased molecular dynamics simulations and enhanced sampling methods such as umbrella sampling, metadynamics and replica exchange are described using key examples. This review demonstrates how state-of-the-art molecular simulations have been used successfully to describe the mechanism of binding and permeation of small molecules with biological membranes, as well as associated changes to the structure and dynamics of these membranes. The review concludes with an outlook on future directions in this field.

Revisiting the Interaction of Melittin with Phospholipid Bilayers: The Effects of Concentration and Ionic Strength.

Melittin is an anti-microbial peptide (AMP) and one of the most studied membrane-disrupting peptides. There is, however, a lack of accurate measurements of the concentration-dependent kinetics and affinity of binding of melittin to phospholipid membranes. In this study, we used surface plasmon resonance spectroscopy to determine the concentration-dependent effect on the binding of melittin to 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) bilayers in vesicles. Three concentration ranges were considered, and when combined, covered two orders of magnitudes (0.04 µM to 8 µM), corresponding to concentrations relevant to the membrane-disrupting and anti-microbial activities of melittin. Binding kinetics data were analysed using a 1:1 Langmuir-binding model and a two-state reaction model. Using in-depth quantitative analysis, we characterised the effect of peptide concentration, the addition of NaCl at physiological ionic strength and the choice of kinetic binding model on the reliability of the calculated kinetics and affinity of binding parameters. The apparent binding affinity of melittin for POPC bilayers was observed to decrease with increasing peptide/lipid (P/L) ratio, primarily due to the marked decrease in the association rate. At all concentration ranges, the two-state reaction model provided a better fit to the data and, thus, a more reliable estimate of binding affinity. Addition of NaCl significantly reduced the signal response during the association phase; however, no substantial effect on the binding affinity of melittin to the POPC bilayers was observed. These findings based on POPC bilayers could have important implications for our understanding of the mechanism of action of melittin on more complex model cell membranes of higher physiological relevance.

Does Sucrose Change Its Mechanism of Stabilization of Lipid Bilayers during Desiccation? Influences of Hydration and Concentration.

Sugar-membrane interactions are believed to be responsible for cell preservation during desiccation and freezing, but the molecular mechanism by which they achieve this is still not well understood. The associated decrease of the main phase transition temperature of phospholipid bilayers is explained by two opposing views on the matter: the direct sugar-phospholipid interaction at the bilayer interface (water replacement hypothesis) and an entropy-driven phase transition with sugar molecules concentrating away from the lipid interface (hydration forces explanation). Both mechanisms are supported by experiments but molecular dynamics (MD) simulations have overwhelmingly shown the occurrence of direct sugar-phospholipid interactions. We have performed MD simulations of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) bilayers at different water and sucrose contents. The behavior of sucrose was found to depend on both the sucrose and water contents: at high sucrose concentration and at low hydration, it is best described by the hydration forces explanation model, whereas at low sucrose concentration, it is consistent with the water replacement hypothesis model. These simulations reveal that at low concentration, sucrose molecules preferentially interact directly with the membrane interface while at high concentration, they preferentially accumulate in the intermembrane solution. The transition between the two modes of interaction is revealed for the first time as being governed by the saturation of the lipid bilayer interface with sucrose molecules, and this occurs more rapidly as the level of hydration decreases.

Influence of Bilayer Size and Number in Multi-Bilayer DOPC Simulations at Full and Low Hydration.

Biophysical studies of model cell membranes at full and low hydration are usually carried out using scattering measurements on multi-bilayer systems. Molecular simulations of lipid bilayers aimed at reproducing those experimental conditions are usually conducted using single bilayers with different amounts of water. These simulation conditions may lead to artifacts arising from size effects and self-interactions because of periodic boundary conditions. We have tested the influence of the size and number of bilayers on membrane properties using the Lipid14 force field for lipids in molecular dynamics simulations of 1,2-dioleoyl- sn-glycero-3-phosphocholine bilayers at full hydration (44 water molecules per lipid), low hydration (18 water molecules per lipid), and dehydration (9 water molecules per lipid). A number of additional simulations were conducted with the Slipids force field for comparison. We have found that the average area per lipid (APL), thickness, mass density profiles, and acyl tail order parameters are insensitive to the size and the number of bilayers for all hydration states. The Lipid14 force field can also successfully reproduce the experimentally observed decrease in APL and corresponding increase in bilayer thickness upon dehydration, reflecting the increase in ordering as the system becomes more gel-like. Additionally, decreasing hydration levels were associated with a trend away from normal lateral diffusion and toward more subdiffusive regimes across both force fields. In summary, at least for the Lipid14 force field, the use of a single bilayer with 128 phospholipid molecules provides an adequate representation of multi-bilayer systems at varying levels of hydration.

Secondary Structural Changes in Proteins as a Result of Electroadsorption at Aqueous-Organogel Interfaces.

The electroadsorption of proteins at aqueous-organic interfaces offers the possibility to examine protein structural rearrangements upon interaction with lipophilic phases, without modifying the bulk protein or relying on a solid support. The aqueous-organic interface has already provided a simple means of electrochemical protein detection, often involving adsorption and ion complexation; however, little is yet known about the protein structure at these electrified interfaces. This work focuses on the interaction between proteins and an electrified aqueous-organic interface via controlled protein electroadsorption. Four proteins known to be electroactive at such interfaces were studied: lysozyme, myoglobin, cytochrome c, and hemoglobin. Following controlled protein electroadsorption onto the interface, ex situ structural characterization of the proteins by FTIR spectroscopy was undertaken, focusing on secondary structural traits within the amide I band. The structural variations observed included unfolding to form aggregated antiparallel ß-sheets, where the rearrangement was specifically dependent on the interaction with the organic phase. This was supported by MALDI ToF MS measurements, which showed the formation of protein-anion complexes for three of these proteins, and molecular dynamic simulations, which modeled the structure of lysozyme at an aqueous-organic interface. On the basis of these findings, the modulation of protein secondary structure by interfacial electrochemistry opens up unique prospects to selectively modify proteins.

The Biological and Biophysical Properties of the Spider Peptide Gomesin.

This review summarises the current knowledge of Gomesin (Gm), an 18-residue long, cationic anti-microbial peptide originally isolated from the haemocytes of the Brazilian tarantula Acanthoscurria gomesiana. The peptide shows potent cytotoxic activity against clinically relevant microbes including Gram-positive and Gram-negative bacteria, fungi, and parasites. In addition, Gm shows in-vitro and in-vivo anti-cancer activities against several human and murine cancers. The peptide exerts its cytotoxic activity by permeabilising cell membranes, but the underlying molecular mechanism of action is still unclear. Due to its potential as a therapeutic agent, the structure and membrane-binding properties, as well as the leakage and cytotoxic activities of Gm have been studied using a range of techniques. This review provides a summary of these studies, with a particular focus on biophysical characterisation studies of peptide variants that have attempted to establish a structure-activity relationship. Future studies are still needed to rationalise the binding affinity and cell-type-specific selectivity of Gm and its variants, while more pre-clinical studies are required to develop Gm into a therapeutically useful peptide.

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.