Exploring the inter-molecular interactions in amyloid-β protofibril with molecular dynamics simulations and molecular mechanics Poisson-Boltzmann surface area free energy calculations

Effects of Nitro-Oxidative Stress on Biomolecules: Part 1-Non-Reactive Molecular Dynamics Simulations

Plasma medicine, or the biomedical application of cold atmospheric plasma (CAP), is an expanding field within plasma research. CAP has demonstrated remarkable versatility in diverse biological applications, including cancer treatment, wound healing, microorganism inactivation, and skin disease therapy. However, the precise mechanisms underlying the effects of CAP remain incompletely understood. The therapeutic effects of CAP are largely attributed to the generation of reactive oxygen and nitrogen species (RONS), which play a crucial role in the biological responses induced by CAP. Specifically, RONS produced during CAP treatment have the ability to chemically modify cell membranes and membrane proteins, causing nitro-oxidative stress, thereby leading to changes in membrane permeability and disruption of cellular processes. To gain atomic-level insights into these interactions, non-reactive molecular dynamics (MD) simulations have emerged as a valuable tool. These simulations facilitate the examination of larger-scale system dynamics, including protein-protein and protein-membrane interactions. In this comprehensive review, we focus on the applications of non-reactive MD simulations in studying the effects of CAP on cellular components and interactions at the atomic level, providing a detailed overview of the potential of CAP in medicine. We also review the results of other MD studies that are not related to plasma medicine but explore the effects of nitro-oxidative stress on cellular components and are therefore important for a broader understanding of the underlying processes.

Destabilization of Aβ fibrils by omega-3 polyunsaturated fatty acids: a molecular dynamics study

The senile plaques of neurotoxic aggregates of Aβ protein, deposited extraneuronally, mark the pathological hallmark of Alzheimer's disease (AD). The natural compounds such as omega-3 (ω-3) polyunsaturated fatty acids (PUFAs), which can access blood-brain barrier, are believed to be potential disruptors of preformed Aβ fibrils to cure AD with unknown mechanism. Herein, we present the destabilization potential of three ω-3 PUFAs, viz. Eicosapentaenoic acid (EPA), Docosahexaenoic acid (HXA), and α-linolenic acid (LNL) by molecular dynamics simulation. After an initial testing of 300 ns, EPA and HXA have been considered further for extended production run time, 500 ns. The increased value of root mean square deviation (RMSD), radius of gyration, and solvent-accessible surface area (SASA), the reduced number of H-bonds and β-sheet content, and disruption of salt bridges and hydrophobic contacts establish the binding of these ligands to Aβ fibril leading to destabilization. The polar head was found to interact with positively charged lysine (K28) residue in the fibril. However, the hydrophobicity of the long aliphatic tail competes with the intrinsic hydrophobic interactions of Aβ fibril. This amphiphilic nature of EPA and HXA led to the breaking of inherent hydrophobic contacts and formation of new bonds between the tail of PUFA and hydrophobic residues of Aβ fibril, leading to the destabilization of fibril. The Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) results explain the binding of EPA and HXA to Aβ fibril by interacting with different residues. The destabilization potential of EPA and HXA establishes them as promising drug leads to cure AD, and encourages prospecting of other fatty acids for therapeutic intervention in AD.Communicated by Ramaswamy H. Sarma.

Bionics design of affinity peptide inhibitors for SARS-CoV-2 RBD to block SARS-CoV-2 RBD-ACE2 interactions

Coronavirus Disease 2019 (COVID-19), has already posed serious threats and impacts on the health of the population and the country's economy. Therefore, it is of great theoretical significance and practical application value to better understand the process of COVID-19 infection and develop effective therapeutic drugs. It is known that the receptor-binding structural domain (SARS-CoV-2 RBD) on the spike protein of the novel coronavirus directly mediates its interaction with the host receptor angiotensin-converting enzyme 2 (ACE2), and thus blocking SARS-CoV-2 RBD-ACE2 interaction is capable of inhibiting SARS-CoV-2 infection. Firstly, the interaction mechanism between SARS-CoV-2RBD-ACE2 was explored using molecular dynamics simulation (MD) coupled with molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) free energy calculation method. The results of energy analysis showed that the key residues R403, R408, K417, and Y505 of SARS-CoV-2 RBD and the key residues D30, E37, D38, and Y41 of ACE2 were identified. Therefore, according to the hotspot residues of ACE2 and their distribution, a short peptide library of high-affinity SARS-CoV-2 RBD was constructed. And by using molecular docking virtual screening, six short peptides including DDFEDY, DEFEDY, DEYEDY, DFVEDY, DFHEDY, and DSFEDY with high affinity for SARS-CoV-2 RBD were identified. The results of MD simulation further confirmed that DDFEDY, DEYEDY, and DFVEDY are expected to be effective inhibitors. Finally, the allergenicity, toxicity and solubility properties of the three peptide inhibitors were validated.

Potential Anti-Alzheimer Agents from Guanidinyl Tryptophan Derivatives with Activities of Membrane Adhesion and Conformational Transition Inhibitions

Guanidinyl tryptophan derivatives TGN1, TGN2, TGN3, and TGN4 were synthesized, and these compounds were shown to possess in vitro inhibitory activity for amyloid aggregation in a previous study. Nevertheless, the influence of the TGN series of compounds on the binding and permeation behaviors of an Aβ monomer to the cell membranes was not elucidated. In this study, we investigated the effect of compounds in the TGN series on the behavior of an Aβ monomer regarding its toxicity toward the bilayer lipid membrane using molecular dynamics (MD) simulation. MD simulations suggest that TGN4 is a potential agent that can interfere with the movement of the Aβ monomer into the membrane. The MM-GBSA result demonstrated that TGN4 exhibits the highest affinity to the Aβ_1–42 monomer but has the lowest affinity to the bilayer. Moreover, TGN4 also contributes to a decrease in the binding affinity between the Aβ_1–42 monomer and the POPC membrane. Regarding the results of the binding mode and conformational analyses, a high number of amino-acid residues were shown to provide the binding interactions between TGN4 and the Aβ_1–42 monomer. TGN4 also reduces the conformational transition of the Aβ_1–42 monomer by means of interacting with the monomer. The present study presents molecular-level insights into how the TGN series of compounds affect the membrane adsorption and the conformational transition of the Aβ_1–42 monomer, which could be valuable for the further development of new anti-Alzheimer agents.

A theoretical study of polymorphism in VQIVYK fibrils

The VQIVYK fragment from the Tau protein, also known as PHF6, is essential for aggregation of Tau into neurofibrillary lesions associated with neurodegenerative diseases. VQIVYK itself forms amyloid fibrils composed of paired β-sheets. Therefore, the full Tau protein and VQIVYK fibrils have been intensively investigated. A central issue in these studies is polymorphism, the ability of a protein to fold into more than one structure. Using all-atom molecular simulations, we generate five stable polymorphs of VQIVYK fibrils, establish their relative free energy with umbrella sampling methods, and identify the side chain interactions that provide stability. The two most stable polymorphs, which have nearly equal free energy, are formed by interdigitation of the mostly hydrophobic VIY "face" sides of the β-sheets. Another stable polymorph is formed by interdigitation of the QVK "back" sides. When we turn to examine structures from cryo-electron microscopy experiments on Tau filaments taken from diseased patients or generated in vitro, we find that the pattern of side chain interactions found in the two most stable face-to-face as well as the back-to-back polymorphs are recapitulated in amyloid structures of the full protein. Thus, our studies suggest that the interactions stabilizing PHF6 fibrils explain the amyloidogenicity of the VQIVYK motif within the full Tau protein and provide justification for the use of VQIVYK fibrils as a test bed for the design of molecules that identify or inhibit amyloid structures.

Benefits under the Sea: The Role of Marine Compounds in Neurodegenerative Disorders

Marine habitats offer a rich reservoir of new bioactive compounds with great pharmaceutical potential; the variety of these molecules is unique, and its production is favored by the chemical and physical conditions of the sea. It is known that marine organisms can synthesize bioactive molecules to survive from atypical environmental conditions, such as oxidative stress, photodynamic damage, and extreme temperature. Recent evidence proposed a beneficial role of these compounds for human health. In particular, xanthines, bryostatin, and 11-dehydrosinulariolide displayed encouraging neuroprotective effects in neurodegenerative disorders. This review will focus on the most promising marine drugs' neuroprotective potential for neurodegenerative disorders, such as Parkinson's and Alzheimer's diseases. We will describe these marine compounds' potential as adjuvant therapies for neurodegenerative diseases, based on their antioxidant, anti-inflammatory, and anti-apoptotic properties.

Destabilization of Alzheimer's Aβ42 protofibrils with acyclovir, carmustine, curcumin, and tetracycline: insights from molecular dynamics simulations

Among the neurodegenerative diseases, one of the most common dementia is Alzheimer's disease (AD).

Expression and purification of amyloid β-protein, tau, and α-synuclein in Escherichia coli: a review

Misfolding and accumulation of amyloidogenic proteins into various forms of aggregated intermediates and insoluble amyloid fibrils is associated with more than 50 human diseases. Large amounts of high-quality amyloid proteins are required for better probing of their aggregation and neurotoxicity. Due to their intrinsic hydrophobicity, it is a challenge to obtain amyloid proteins with high yield and purity, and they have attracted the attention of researchers from all over the world. The rapid development of bioengineering technology provides technical support for obtaining large amounts of recombinant amyloidogenic proteins. This review discusses the available expression and purification methods for three amyloid proteins including amyloid β-protein, tau, and α-synuclein in microbial expression systems, especially Escherichia coli, and discusses the advantages and disadvantages of these methods. Importantly, these protocols can also be referred to for the expression and purification of other hydrophobic proteins.

Destabilization potential of phenolics on Aβ fibrils: mechanistic insights from molecular dynamics simulation

Ellagic acid from pomegranate and walnuts is found to destabilize Aβ fibrils. It can be a potential drug to treat AD.

Inhibitory Effect of a Flavonoid Dihydromyricetin against Aβ40 Amyloidogenesis and Its Associated Cytotoxicity

Misfolding and fibrillogenesis of amyloid-β protein (Aβ) play a key role in the onset and progression of Alzheimer's disease (AD). Screening for inhibitors against Aβ amyloidogenesis is helpful for rational designing and developing new anti-AD drugs and therapeutic strategies. Dihydromyricetin, a natural flavonoid extracted from a Chinese herb, Ampelopsis grossedentata, has been proven with antioxidative, anti-inflammatory, and neuroprotective effects against neurodegenerative disease. Herein, we found that dihydromyricetin could inhibit Aβ40 aggregation, impede the protofibril formation, disassemble preformed Aβ40 fibrils, and protect PC12 cells from the Aβ40-induced cytotoxicity using a series of biochemical and biophysical assays, including thioflavin T fluorescence, atomic force microscopy, and cell toxicity assays. Circular dichroism spectroscopy data proved that dihydromyricetin delayed the Aβ40 conformational conversion. In addition, the results of molecular dynamics simulations indicated that the interaction between dihydromyricetin and Aβ40 trimer is mainly nonpolar interactions. Key residues (i.e., V18, A21, and D23) of the Aβ40 interacting with dihydromyricetin were also identified. This study suggested that dihydromyricetin shows great potential to be developed as a novel Aβ40 inhibitor.

Edaravone inhibits the conformational transition of amyloid-β42: insights from molecular dynamics simulations

Previous work has shown that edaravone inhibits fibrillogenesis of amyloid-β protein (Aβ). However, the detailed mechanism by which edaravone inhibits the conformational transition of the Aβ42 monomer is not known at the molecular level. Here, explicit-solvent molecular dynamics (MD) simulations were coupled with molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method to address the issue. MD simulations confirmed that edaravone inhibits the conformational transition of the Aβ42 monomer in a dose-dependent manner. It was found that the direct interactions between edaravone and Aβ42 are responsible for its inhibiting effects. The analysis of binding free energy using the MM-PBSA method demonstrated that the nonpolar interactions provide favourable contributions (about -71.7 kcal/mol). Conversely, the polar interactions are unfavourable for the binding process. A total of 14 residues were identified as greatly contributing to the binding free energy between edaravone and the Aβ42 monomer. In addition, the intra-peptide hydrophobic interactions were weakened and the salt bridge D23-K28 was interrupted by edaravone. Therefore, the conformational transition was inhibited. Our studies provide molecular-level insights into how edaravone molecules inhibit the conformational transition of the Aβ42 monomer, which may be useful for designing amyloid inhibitors.Communicated by Ramaswamy H. Sarma.

Oxidation destabilizes toxic amyloid beta peptide aggregation

The aggregation of insoluble amyloid beta (Aβ) peptides in the brain is known to trigger the onset of neurodegenerative diseases, such as Alzheimer’s disease. In spite of the massive number of investigations, the underlying mechanisms to destabilize the Aβ aggregates are still poorly understood. Some studies indicate the importance of oxidation to destabilize the Aβ aggregates. In particular, oxidation induced by cold atmospheric plasma (CAP) has demonstrated promising results in eliminating these toxic aggregates. In this paper, we investigate the effect of oxidation on the stability of an Aβ pentamer. By means of molecular dynamics simulations and umbrella sampling, we elucidate the conformational changes of Aβ pentamer in the presence of oxidized residues, and we estimate the dissociation free energy of the terminal peptide out of the pentamer form. The calculated dissociation free energy of the terminal peptide is also found to decrease with increasing oxidation. This indicates that Aβ pentamer aggregation becomes less favorable upon oxidation. Our study contributes to a better insight in one of the potential mechanisms for inhibition of toxic Aβ peptide aggregation, which is considered to be the main culprit to Alzheimer’s disease.

Amyloidogenicity and Cytotoxicity of a Recombinant C-Terminal His6-Tagged Aβ1-42

Aggregation of amyloid β peptide (Aβ) is closely associated with the occurrence and development of Alzheimer's disease (AD). Reproducible and detailed studies on the aggregation kinetics and structure of various aggregates have been conducted using recombinant Aβ peptides. While the His₆-tag is commonly used in the purification of recombinant proteins due to its great simplicity and affinity, there is little information on the aggregation of His₆-tagged Aβ and its corresponding cytotoxicity. Moreover, it is also unclear whether there is an effect of the His₆-tag on the amyloidogenicity and cytotoxicity of recombinant Aβ_1-42. Herein, a method to express and purify a mutant C-terminally His₆-tagged Aβ_1-42 (named as Aβ_1-42-His₆) from Escherichia coli was described. Aβ_1-42-His₆ aggregated into β-sheet-rich fibrils as shown by thioflavin T fluorescence, atomic force microscopy and circular dichroism spectroscopy. Moreover, the fibrillar recombinant Aβ_1-42-His₆ showed strong toxicity toward PC12 cells in vitro. Molecular dynamics simulations revealed that the His₆-tag contributed little to the secondary structure and intermolecular interactions, including hydrophobic interactions, salt bridges, and hydrogen bonding of the fibrillar pentamer of Aβ_1-42. This highlights the biological importance of modification on the molecular structure of Aβ. Thus, the easily purified high-quality Aβ_1-42-His₆ offers great advantages for screening aggregation inhibitors or in vitro confirmation of rationally designed drugs for the treatment of AD.

9-Methylfascaplysin Is a More Potent Aβ Aggregation Inhibitor than the Marine-Derived Alkaloid, Fascaplysin, and Produces Nanomolar Neuroprotective Effects in SH-SY5Y Cells

β-Amyloid (Aβ) is regarded as an important pathogenic target for Alzheimer’s disease (AD), the most prevalent neurodegenerative disease. Aβ can assemble into oligomers and fibrils, and produce neurotoxicity. Therefore, Aβ aggregation inhibitors may have anti-AD therapeutic efficacies. It was found, here, that the marine-derived alkaloid, fascaplysin, inhibits Aβ fibrillization in vitro. Moreover, the new analogue, 9-methylfascaplysin, was designed and synthesized from 5-methyltryptamine. Interestingly, 9-methylfascaplysin is a more potent inhibitor of Aβ fibril formation than fascaplysin. Incubation of 9-methylfascaplysin with Aβ directly reduced Aβ oligomer formation. Molecular dynamics simulations revealed that 9-methylfascaplysin might interact with negatively charged residues of Aβ₄₂ with polar binding energy. Hydrogen bonds and π–π interactions between the key amino acid residues of Aβ₄₂ and 9-methylfascaplysin were also suggested. Most importantly, compared with the typical Aβ oligomer, Aβ modified by nanomolar 9-methylfascaplysin produced less neuronal toxicity in SH-SY5Y cells. 9-Methylfascaplysin appears to be one of the most potent marine-derived compounds that produces anti-Aβ neuroprotective effects. Given previous reports that fascaplysin inhibits acetylcholinesterase and induces P-glycoprotein, the current study results suggest that fascaplysin derivatives can be developed as novel anti-AD drugs that possibly act via inhibition of Aβ aggregation along with other target mechanisms.

Hydroxylated Single-Walled Carbon Nanotubes Inhibit Aβ42 Fibrillogenesis, Disaggregate Mature Fibrils, and Protect against Aβ42-Induced Cytotoxicity

The fibrillogenesis of amyloid-β protein (Aβ) is considered a crucial factor in the pathogenesis of Alzheimer's disease (AD). Hence, inhibiting Aβ fibrillogenesis is regarded as the primary therapeutic strategy for the prevention and treatment of AD. However, the development of effective inhibitors against Aβ fibrillogenesis has faced significant challenges. Previous studies have shown that pristine single-walled carbon nanotubes (SWNTs) can inhibit fibrillogenesis of some amyloid proteins. However, the poor dispersibility of SWNTs in an aqueous environment greatly hinders their inhibitory efficacy. Here, we examined the inhibitory activity of hydroxylated SWNTs (SWNT-OH) on the aggregation and cytotoxicity of Aβ₄₂ using thioflavin T (ThT) fluorescence, atomic force microscopy (AFM), cellular viability assays, and molecular dynamics (MD) simulations. ThT and AFM results showed that SWNT-OH inhibits Aβ₄₂ fibrillogenesis and disaggregates preformed amyloid fibrils in a dose-dependent manner. Furthermore, the ratio of hydroxyl groups in SWNT-OH is crucial for their effect against Aβ₄₂ aggregation. SWNT-OH exerted cytoprotective effects against Aβ₄₂ fibrillation-induced cytotoxicity. The results of free-energy decomposition studies based on MD simulations revealed that nonpolar interactions, and especially van der Waals forces, contributed most of the free energy of binding in the SWNT-OH-Aβ complex. Two regions of the Aβ pentamer were identified to interact with SWNT-OH, spanning H13-Q15 and V36-G38. The findings presented here will contribute to a comprehensive understanding of the inhibitory effect of hydroxylated nanoparticles against Aβ fibrillogenesis, which is critical for the search for more effective agents that can counteract amyloid-mediated pathologies.

Understanding the microscopic origin behind heterogeneous properties of water confined in and around Aβ17-42 protofilaments

Aggregation of amyloid beta (Aβ) peptides in the brain is responsible for one of the most devastating neurodegenerative diseases, namely, Alzheimer's disease. In this study, we have carried out atomistic molecular dynamics simulations to explore the effects of non-uniform structural distortions of Aβ_17-42 pre-fibrillar aggregates of different sizes on the microscopic structure and ordering of water molecules confined within their amphiphilic nanocores. The calculations revealed non-uniform peptide-water interactions resulting in simultaneous existence of both highly ordered and disordered water molecules within the spatially heterogeneous confined environment of the protofilament cores. It is found that the high degree of ordering originates from a sizable fraction of doubly coordinated core water molecules, while the randomly oriented ones are those that are coordinated with three neighbors in their first coordination shells. Furthermore, it is quantitatively demonstrated that relative fractions of these two types of water molecules are correlated with the protofilament core topology and the degree of confinement within that. It is proposed that the ordered core waters are likely to stabilize the Aβ protofilaments by screening the residue charges and favoring water-mediated salt bridge formations, while the randomly oriented ones can drive further growth of the protofilaments by being displaced easily during the docking of additional peptides. In that way, both types of core water molecules can play equally important roles in controlling the growth and stability of the Aβ-aggregates.

Nature-Inspired Construction of Two-Dimensionally Self-Assembled Peptide on Pristine Graphene

Peptide assemblies have received significant attention because of their important role in biology and applications in bionanotechnology. Despite recent efforts to elucidate the principles of peptide self-assembly for developing novel functional devices, peptide self-assembly on two-dimensional nanomaterials has remained challenging. Here, we report nature-inspired two-dimensional peptide self-assembly on pristine graphene via optimization of peptide-peptide and peptide-graphene interactions. Two-dimensional peptide self-assembly was designed based on statistical analyses of >10⁴ protein structures existing in nature and atomistic simulation-based structure predictions. We characterized the structures and surface properties of the self-assembled peptide formed on pristine graphene. Our study provides insights into the formation of peptide assemblies coupled with two-dimensional nanomaterials for further development of nanobiocomposite devices.

Scrutiny of the mechanism of small molecule inhibitor preventing conformational transition of amyloid-β42 monomer: insights from molecular dynamics simulations

Alzheimer's disease (AD) is a progressive neurodegenerative disorder that is characterized by loss of intellectual functioning of brain and memory loss. According to amyloid cascade hypothesis, aggregation of amyloid-β₄₂ (Aβ₄₂) peptide can generate toxic oligomers and their accumulation in the brain is responsible for the onset of AD. In spite of carrying out a large number of experimental studies on inhibition of Aβ₄₂ aggregation by small molecules, the detailed inhibitory mechanism remains elusive. In the present study, comparable molecular dynamics (MD) simulations were performed to elucidate the inhibitory mechanism of a sulfonamide inhibitor C1 (2,5-dichloro-N-(4-piperidinophenyl)-3-thiophenesulfonamide), reported for its in vitro and in vivo anti-aggregation activity against Aβ₄₂. MD simulations reveal that C1 stabilizes native α-helix conformation of Aβ₄₂ by interacting with key residues in the central helix region (13-26) with hydrogen bonds and π-π interactions. C1 lowers the solvent-accessible surface area of the central hydrophobic core (CHC), KLVFF (16-20), that confirms burial of hydrophobic residues leading to the dominance of helical conformation in the CHC region. The binding free energy analysis with MM-PBSA demonstrates that Ala2, Phe4, Tyr10, Gln15, Lys16, Leu17, Val18, Phe19, Phe20, Glu22, and Met35 contribute maximum to binding free energy (-43.1 kcal/mol) between C1 and Aβ₄₂ monomer. Overall, MD simulations reveal that C1 inhibits Aβ₄₂ aggregation by stabilizing native helical conformation and inhibiting the formation of aggregation-prone β-sheet conformation. The present results will shed light on the underlying inhibitory mechanism of small molecules that show potential in vitro anti-aggregation activity against Aβ₄₂.

Dimerization Mechanism of Alzheimer Aβ40 Peptides: The High Content of Intrapeptide-Stabilized Conformations in A2V and A2T Heterozygous Dimers Retards Amyloid Fibril Formation

Amyloid beta (Aβ) oligomerization is associated with the origin and progression of Alzheimer's disease (AD). While the A2V mutation enhances aggregation kinetics and toxicity, mixtures of wild-type (WT) and A2V, and also WT and A2T, peptides retard fibril formation and protect against AD. In this study, we simulate the equilibrium ensemble of WT:A2T Aβ₄₀ dimer by means of extensive atomistic replica exchange molecular dynamics and compare our results with previous equivalent simulations of A2V:A2V, WT:WT, and WT:A2V Aβ₄₀ dimers for a total time scale of nearly 0.1 ms. Qualitative comparison of the resulting thermodynamic properties, such as the relative binding free energies, with the reported experimental kinetic and thermodynamic data affords us important insight into the conversion from slow-pathway to fast-pathway dimer conformations. The crucial reaction coordinate or driving force of such transformation turns out to be related to hydrophobic interpeptide interactions. Analysis of the equilibrium ensembles shows that the fast-pathway conformations contain interpeptide out-of-register antiparallel β-sheet structures at short interpeptide distances. In contrast, the slow-pathway conformations are formed by the association of peptides at large interpeptide distances and high intrapeptide compactness, such as conformations containing intramolecular three-stranded β-sheets which sharply distinguish fast (A2V:A2V and WT:WT) and slow (WT:A2T and WT:A2V) amyloid-forming sequences. Also, this analysis leads us to predict that a molecule stabilizing the intramolecular three-stranded β-sheet or inhibiting the formation of an interpeptide β-sheet spanning residues 17-20 and 31-37 would further reduce fibril formation and probably the cytotoxicity of Aβ species.

Role of sequence and structural polymorphism on the mechanical properties of amyloid fibrils

Amyloid fibrils playing a critical role in disease expression, have recently been found to exhibit the excellent mechanical properties such as elastic modulus in the order of 10 GPa, which is comparable to that of other mechanical proteins such as microtubule, actin filament, and spider silk. These remarkable mechanical properties of amyloid fibrils are correlated with their functional role in disease expression. This suggests the importance in understanding how these excellent mechanical properties are originated through self-assembly process that may depend on the amino acid sequence. However, the sequence-structure-property relationship of amyloid fibrils has not been fully understood yet. In this work, we characterize the mechanical properties of human islet amyloid polypeptide (hIAPP) fibrils with respect to their molecular structures as well as their amino acid sequence by using all-atom explicit water molecular dynamics (MD) simulation. The simulation result suggests that the remarkable bending rigidity of amyloid fibrils can be achieved through a specific self-aggregation pattern such as antiparallel stacking of β strands (peptide chain). Moreover, we have shown that a single point mutation of hIAPP chain constituting a hIAPP fibril significantly affects the thermodynamic stability of hIAPP fibril formed by parallel stacking of peptide chain, and that a single point mutation results in a significant change in the bending rigidity of hIAPP fibrils formed by antiparallel stacking of β strands. This clearly elucidates the role of amino acid sequence on not only the equilibrium conformations of amyloid fibrils but also their mechanical properties. Our study sheds light on sequence-structure-property relationships of amyloid fibrils, which suggests that the mechanical properties of amyloid fibrils are encoded in their sequence-dependent molecular architecture.

In silico cross seeding of Aβ and amylin fibril-like oligomers

Recent epidemiological data have shown that patients suffering from Type 2 Diabetes Mellitus have an increased risk to develop Alzheimer's disease and vice versa. A possible explanation is the cross-sequence interaction between Aβ and amylin. Because the resulting amyloid oligomers are difficult to probe in experiments, we investigate stability and conformational changes of Aβ-amylin heteroassemblies through molecular dynamics simulations. We find that Aβ is a good template for the growth of amylin and vice versa. We see water molecules permeate the β-strand-turn-β-strand motif pore of the oligomers, supporting a commonly accepted mechanism for toxicity of β-rich amyloid oligomers. Aiming for a better understanding of the physical mechanisms of cross-seeding and cell toxicity of amylin and Aβ aggregates, our simulations also allow us to identify targets for the rational design of inhibitors against toxic fibril-like oligomers of Aβ and amylin oligomers.