Mechanism of amyloid-β fibril elongation

S100A9 inhibits and redirects prion protein 89-230 fragment amyloid aggregation

Protein aggregation in the form of amyloid fibrils has long been associated with the onset and development of various amyloidoses, including Alzheimer's, Parkinson's or prion diseases. Recent studies of their fibril formation process have revealed that amyloidogenic protein cross-interactions may impact aggregation pathways and kinetic parameters, as well as the structure of the resulting aggregates. Despite a growing number of reports exploring this type of interaction, they only cover just a small number of possible amyloidogenic protein pairings. One such pair is between two neurodegeneration-associated proteins: the pro-inflammatory S100A9 and prion protein, which are known to co-localize in vivo. In this study, we examined their cross-interaction in vitro and discovered that the fibrillar form of S100A9 modulated the aggregation pathway of mouse prion protein 89-230 fragment, while non-aggregated S100A9 also significantly inhibited its primary nucleation process. These results complement previous observations of the pro-inflammatory protein's role in amyloid aggregation and highlight its potential role against neurodegenerative disorders.

Molecular Insights into the Dynamics of Amyloid Fibril Growth: Elongation and Lateral Assembly of GNNQQNY Protofibrils

The self-assembly of peptides and proteins into β-sheet rich amyloid fibrils is linked to both functional and pathological states. In this study, the growth of fibrillar structures of the short peptide GNNQQNY, a fragment from the yeast prion Sup35 protein, was examined. Molecular dynamics simulations were used to study alternative mechanisms of fibril growth, including elongation through binding of monomers as well as fibril self-assembly into larger, more mature structures. It was found that after binding, monomers diffused along preformed fibrils toward the ends, supporting the mechanism of fibril growth via elongation. Lateral assembly of protofibrils was found to occur readily, suggesting that this could be the key to transitioning from isolated fibrils to mature multilayer structures. Overall, the work provides mechanistic insights into the competitive pathways that govern amyloid fibril growth.

Investigating lysozyme amyloid fibril formation and structural variability dependence on its initial folding state under different pH conditions

Protein fibril formation and accumulation are associated with dozens of amyloidoses, including the widespread and yet-incurable Alzheimer's and Parkinson's diseases. Currently, there are still several aspects of amyloid aggregation that are not fully understood, which negatively contributes to the development of disease-altering drugs and treatments. One factor which requires a more in-depth analysis is the effect of the environment on both the initial state of amyloidogenic proteins and their aggregation process and resulting fibril characteristics. In this work, we examine how lysozyme's folding state influences its amyloid formation kinetics and resulting aggregate structural characteristics under several different pH conditions, ranging from acidic to neutral. We demonstrate that both the initial state of the protein and the solution's pH value have a significant combined effect on the variability of the resulting aggregate secondary structures, as well as their stabilities, interactions with amyloid-specific dye molecules, and self-replication properties.

Unlocking novel therapies: cyclic peptide design for amyloidogenic targets through synergies of experiments, simulations, and machine learning

Existing therapies for neurodegenerative diseases like Parkinson's and Alzheimer's address only their symptoms and do not prevent disease onset. Common therapeutic agents, such as small molecules and antibodies struggle with insufficient selectivity, stability and bioavailability, leading to poor performance in clinical trials. Peptide-based therapeutics are emerging as promising candidates, with successful applications for cardiovascular diseases and cancers due to their high bioavailability, good efficacy and specificity. In particular, cyclic peptides have a long in vivo stability, while maintaining a robust antibody-like binding affinity. However, the de novo design of cyclic peptides is challenging due to the lack of long-lived druggable pockets of the target polypeptide, absence of exhaustive conformational distributions of the target and/or the binder, unknown binding site, methodological limitations, associated constraints (failed trials, time, money) and the vast combinatorial sequence space. Hence, efficient alignment and cooperation between disciplines, and synergies between experiments and simulations complemented by popular techniques like machine-learning can significantly speed up the therapeutic cyclic-peptide development for neurodegenerative diseases. We review the latest advancements in cyclic peptide design against amyloidogenic targets from a computational perspective in light of recent advancements and potential of machine learning to optimize the design process. We discuss the difficulties encountered when designing novel peptide-based inhibitors and we propose new strategies incorporating experiments, simulations and machine learning to design cyclic peptides to inhibit the toxic propagation of amyloidogenic polypeptides. Importantly, these strategies extend beyond the mere design of cyclic peptides and serve as template for the de novo generation of (bio)materials with programmable properties.

A Light Scattering Investigation of Enzymatic Gelation in Self-Assembling Peptides

Self-assembling peptides (SAPs) have been increasingly studied as hydrogel-former gelators because they can create biocompatible environments. A common strategy to trigger gelation, is to use a pH variation, but most methods result in a change in pH that is too rapid, leading to gels with hardly reproducible properties. Here, we use the urea-urease reaction to tune gel properties, by a slow and uniform pH increase. We were able to produce very homogeneous and transparent gels at several SAP concentrations, ranging from c=1g/L to c=10g/L. In addition, by exploiting such a pH control strategy, and combining photon correlation imaging with dynamic light scattering measurements, we managed to unravel the mechanism by which gelation occurs in solutions of (LDLK)3-based SAPs. We found that, in diluted and concentrated solutions, gelation follows different pathways. This leads to gels with different microscopic dynamics and capability of trapping nanoparticles. At high concentrations, a strong gel is formed, made of relatively thick and rigid branches that firmly entrap nanoparticles. By contrast, the gel formed in dilute conditions is weaker, characterized by entanglements and crosslinks of very thin and flexible filaments. The gel is still able to entrap nanoparticles, but their motion is not completely arrested. These different gel morphologies can potentially be exploited for controlled multiple drug release.

Combination of structure-based virtual screening, molecular docking and molecular dynamics approaches for the discovery of anti-prion fibril flavonoids

Prion diseases are a group of rare neurodegenerative diseases caused by the structural conversion of cellular prion into Scrapie prion resulting aggregated fibrils. Therapy of prion diseases has been developed for several decades, especially drug designs based on the structure of prion monomers. Unfortunately, none of the designed anti-prion drugs function well clinically. To fight against prion fibrils, a drug design based on the precise structure of mammalian prion fibrils is highly required. Fortunately, based on the advantage of newly advanced cryo-electron microscopy (cryo-EM) in the deconvolution of large complexes, three prion fibril structures were resolved in the last 2 years. Based on the cryo-EM solved prion fibril structures, we are able to find some molecules fighting against prion fibrils. Quercetin, one flavonoid molecule in the polyphenol group, has been found to disaggregate the prion fibrils in vitro. In this study, we performed the molecular docking and molecular dynamics simulation on quercetin-like molecules possessing pharmacological properties to evaluate the anti-prion ability of tested molecules. As a result, four quercetin-like molecules interact with prion fibril and decrease the β-strand content by converting some β-strands into loop and helical structures to disintegrate the existing fibril structure. The results of this study are significant in the treatment of prion diseases, and the approaches used in this study are applicable to other amyloid diseases.

Higher-order interactions of Bcr-Abl can broaden chronic myeloid leukemia (CML) drug repertoire

Bcr-Abl, a nonreceptor tyrosine kinase, is associated with leukemias, especially chronic myeloid leukemia (CML). Deletion of Abl's N-terminal region, to which myristoyl is linked, renders the Bcr-Abl fusion oncoprotein constitutively active. The substitution of Abl's N-terminal region by Bcr enables Bcr-Abl oligomerization. Oligomerization is critical: it promotes clustering on the membrane, which is essential for potent MAPK signaling and cell proliferation. Here we decipher the Bcr-Abl specific, step-by-step oligomerization process, identify a specific packing surface, determine exactly how the process is structured and identify its key elements. Bcr's coiled coil (CC) domain at the N-terminal controls Bcr-Abl oligomerization. Crystallography validated oligomerization via Bcr-Abl dimerization between two Bcr CC domains, with tetramerization via tight packing between two binary assemblies. However, the structural principles guiding Bcr CC domain oligomerization are unknown, hindering mechanistic understanding and drugs exploiting it. Using molecular dynamics (MD) simulations, we determine that the binary complex of the Bcr CC domain serves as a basic unit in the quaternary complex providing a specific surface for dimer-dimer packing and higher-order oligomerization. We discover that the small α1-helix is the key. In the binary assembly, the helix forms interchain aromatic dimeric packing, and in the quaternary assembly, it contributes to the specific dimer-dimer packing. Our mechanism is supported by the experimental literature. It offers the key elements controlling this process which can expand the drug discovery strategy, including by Bcr CC-derived peptides, and candidate residues for small covalent drugs, toward quenching oligomerization, supplementing competitive and allosteric tyrosine kinase inhibitors.

Suppression of amyloid-β fibril growth by drug-engineered polymorph transformation

Fibrillization of the protein amyloid β is assumed to trigger Alzheimer's pathology. Approaches that target amyloid plaques, however, have garnered limited clinical success, and their failures may relate to the scarce understanding of the impact of potential drugs on the intertwined stages of fibrillization. Here, we demonstrate that bexarotene, a T-cell lymphoma medication with known antiamyloid activity both in vitro and in vivo, suppresses amyloid fibrillization by promoting an alternative fibril structure. We employ time-resolved in situ atomic force microscopy to quantify the kinetics of growth of individual fibrils and supplement it with structure characterization by cryo-EM. We show that fibrils with structure engineered by the drug nucleate and grow substantially slower than "normal" fibrils; remarkably, growth remains stunted even in drug-free solutions. We find that the suppression of fibril growth by bexarotene is not because of the drug binding to the fibril tips or to the peptides in the solution. Kinetic analyses attribute the slow growth of drug-enforced fibril polymorph to the distinctive dynamics of peptide chain association to their tips. As an additional benefit, the bexarotene fibrils kill primary rat hippocampal neurons less efficiently than normal fibrils. In conclusion, the suggested drug-driven polymorph transformation presents a mode of action to irreversibly suppress toxic aggregates not only in Alzheimer's but also potentially in myriad diverse pathologies that originate with protein condensation.

Proteinopathies: Deciphering Physiology and Mechanisms to Develop Effective Therapies for Neurodegenerative Diseases

Neurodegenerative diseases (NDs) are a cluster of diseases marked by progressive neuronal loss, axonal transport blockage, mitochondrial dysfunction, oxidative stress, neuroinflammation, and aggregation of misfolded proteins. NDs are more prevalent beyond the age of 50, and their symptoms often include motor and cognitive impairment. Even though various proteins are involved in different NDs, the mechanisms of protein misfolding and aggregation are very similar. Recently, several studies have discovered that, like prions, these misfolded proteins have the inherent capability of translocation from one neuron to another, thus having far-reaching implications for understanding the processes involved in the onset and progression of NDs, as well as the development of innovative therapy and diagnostic options. These misfolded proteins can also influence the transcription of other proteins and form aggregates, tangles, plaques, and inclusion bodies, which then accumulate in the CNS, leading to neuronal dysfunction and neurodegeneration. This review demonstrates protein misfolding and aggregation in NDs, and similarities and differences between different protein aggregates have been discussed. Furthermore, we have also reviewed the disposal of protein aggregates, the various molecular machinery involved in the process, their regulation, and how these molecular mechanisms are targeted to build innovative therapeutic and diagnostic procedures. In addition, the landscape of various therapeutic interventions for targeting protein aggregation for the effective prevention or treatment of NDs has also been discussed.

Rapid restructurization of conformationally-distinct alpha-synuclein amyloid fibrils at an elevated temperature

Protein aggregation in the form of amyloid fibrils is linked with the onset and progression of more than 30 amyloidoses, including multiple neurodegenerative disorders, such as Alzheimer's or Parkinson's disease. Despite countless studies and years of research, the process of such aggregate formation is still not fully understood. One peculiar aspect of amyloids is that they appear to be capable of undergoing structural rearrangements even after the fibrils have already formed. Such a phenomenon was reported to occur in the case of alpha-synuclein and amyloid beta aggregates after a long period of incubation. In this work, we examine whether incubation at an elevated temperature can induce the restructurization of four different conformation alpha-synuclein amyloid fibrils. We show that this structural alteration occurs in a relatively brief time period, when the aggregates are incubated at 60 °C. Additionally, it appears that during this process multiple conformationally-distinct alpha-synuclein fibrils all shift towards an identical secondary structure.

Lysozyme Amyloid Fibril Structural Variability Dependence on Initial Protein Folding State

Amyloid fibril formation is associated with several amyloidoses, including neurodegenerative Alzheimer's or Parkinson's diseases. The process of such fibrillar structure formation is still not fully understood, with new mechanistic insights appearing on a regular basis. This, in turn, has limited the development of potential anti-amyloid compounds, with only a handful of effective cures or treatment modalities available. One of the multiple amyloid aggregation factors that requires further examination is the ability of proteins to form multiple, structurally distinct aggregates, based on the environmental conditions. In this work, we examine how the initial folding state affects the fibrilization of lysozyme-an amyloidogenic protein, often used in protein aggregation studies. We show that there is a correlation between the initial state of the protein and the aggregate formation lag time, rate of elongation, resulting aggregate structural variability and dye-binding properties, as well as formation lag time and rate of elongation.

Energy Landscapes for Base-Flipping in a Model DNA Duplex

We explore the process of base-flipping for four central bases, adenine, guanine, cytosine, and thymine, in a deoxyribonucleic acid (DNA) duplex using the energy landscape perspective. NMR imino-proton exchange and fluorescence correlation spectroscopy studies have been used in previous experiments to obtain lifetimes for bases in paired and extrahelical states. However, the difference of almost 4 orders of magnitude in the base-flipping rates obtained by the two methods implies that they are exploring different pathways and possibly different open states. Our results support the previous suggestion that minor groove opening may be favored by distortions in the DNA backbone and reveal links between sequence effects and the direction of opening, i.e., whether the base flips toward the major or the minor groove side. In particular, base flipping along the minor groove pathway was found to align toward the 5' side of the backbone. We find that bases align toward the 3' side of the backbone when flipping along the major groove pathway. However, in some cases for cytosine and thymine, the base flipping along the major groove pathway also aligns toward the 5' side. The sequence effect may be caused by the polar interactions between the flipping-base and its neighboring bases on either of the strands. For guanine flipping toward the minor groove side, we find that the equilibrium constant for opening is large compared to flipping via the major groove. We find that the estimated rates of base opening, and hence the lifetimes of the closed state, obtained for thymine flipping through small and large angles along the major groove differ by 6 orders of magnitude, whereas for thymine flipping through small angles along the minor groove and large angles along the major groove, the rates differ by 3 orders of magnitude.

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.

Mechanistic Kinetic Model Reveals How Amyloidogenic Hydrophobic Patches Facilitate the Amyloid-β Fibril Elongation

Abnormal aggregation of amyloid β (Aβ) peptides into fibrils plays a critical role in the development of Alzheimer's disease. A two-stage "dock-lock" model has been proposed for the Aβ fibril elongation process. However, the mechanisms of the Aβ monomer-fibril binding process have not been elucidated with the necessary molecular-level precision, so it remains unclear how the lock phase dynamics leads to the overall in-register binding of the Aβ monomer onto the fibril. To gain mechanistic insights into this critical step during the fibril elongation process, we used molecular dynamics (MD) simulations with a physics-based coarse-grained UNited-RESidue (UNRES) force field and sampled extensively the dynamics of the lock phase process, in which a fibril-bound Aβ_(9-40) peptide rearranged to establish the native docking conformation. Analysis of the MD trajectories with Markov state models was used to quantify the kinetics of the rearrangement process and the most probable pathways leading to the overall native docking conformation of the incoming peptide. These revealed a key intermediate state in which an intra-monomer hairpin is formed between the central core amyloidogenic patch ₁₈VFFA₂₁ and the C-terminal hydrophobic patch ₃₄LMVG₃₇. This hairpin structure is highly favored as a transition state during the lock phase of the fibril elongation. We propose a molecular mechanism for facilitation of the Aβ fibril elongation by amyloidogenic hydrophobic patches.

Superoxide dismutase-1 alters the rate of prion protein aggregation and resulting fibril conformation

The assembly of amyloidogenic proteins into highly-structured fibrillar aggregates is related to the onset and progression of several amyloidoses, including neurodegenerative Alzheimer's or Parkinson's diseases. Despite years of research and a general understanding of the process of such aggregate formation, there are currently still very few drugs and treatment modalities available. One of the factors that is relatively insufficiently understood is the cross-interaction between different amyloid-forming proteins. In recent years, it has been shown that several of these proteins or their aggregates can alter each other's fibrillization properties, however, there are still many unknowns in the amyloid interactome. In this work, we examine the interaction between amyloid disease-related prion protein and superoxide dismutase-1. We show that not only does superoxide dismutase-1 increase the lag time of prion protein fibril formation, but it also changes the conformation of the resulting aggregates.

Probing Protein Aggregation Using the Coarse-Grained UNRES Force Field

Protein aggregation is the cause of many, often lethal, diseases, including the Alzheimer's, Parkinson's, and Huntington's diseases, and familial amyloidosis. Theoretical investigation of the mechanism of this process, including the structures of the oligomeric intermediates which are the most toxic, is difficult because of long time scale of aggregation. Coarse-grained models, which enable us to extend the simulation time scale by three or more orders of magnitude, are, therefore, of great advantage in such studies. In this chapter, we describe the application of the physics-based UNited RESidue (UNRES) force field developed in our laboratory to study protein aggregation, in both free simulations and simulations of aggregation propagation from an existing template (seed), and illustrate it with the examples of Aβ-peptide aggregation and Aβ-peptide-assisted aggregation of the peptides derived from the repeat domains of tau (TauRD).

Frustrated peptide chains at the fibril tip control the kinetics of growth of amyloid-β fibrils

Amyloid fibrillization is an exceedingly complex process in which incoming peptide chains bind to the fibril while concertedly folding. The coupling between folding and binding is not fully understood. We explore the molecular pathways of association of Aβ40 monomers to fibril tips by combining time-resolved in situ scanning probe microscopy with molecular modeling. The comparison between experimental and simulation results shows that a complex supported by nonnative contacts is present in the equilibrium structure of the fibril tip and impedes fibril growth in a supersaturated solution. The unraveling of this frustrated state determines the rate of fibril growth. The kinetics of growth of freshly cut fibrils, in which the bulk fibril structure persists at the tip, complemented by molecular simulations, indicate that this frustrated complex comprises three or four monomers in nonnative conformations and likely is contained on the top of a single stack of peptide chains in the fibril structure. This pathway of fibril growth strongly deviates from the common view that the conformational transformation of each captured peptide chain is templated by the previously arrived peptide. The insights into the ensemble structure of the frustrated complex may guide the search for suppressors of Aβ fibrillization. The uncovered dynamics of coupled structuring and assembly during fibril growth are more complex than during the folding of most globular proteins, as they involve the collective motions of several peptide chains that are not guided by a funneled energy landscape.

Multiscale Models for Fibril Formation: Rare Events Methods, Microkinetic Models, and Population Balances

Amyloid fibrils are thought to grow by a two-step dock-lock mechanism. However, previous simulations of fibril formation (i) overlook the bi-molecular nature of the docking step and obtain rates with first-order units, or (ii) superimpose the docked and locked states when computing the potential of mean force for association and thereby muddle the docking and locking steps. Here, we developed a simple microkinetic model with separate locking and docking steps and with the appropriate concentration dependences for each step. We constructed a simple model comprised of chiral dumbbells that retains qualitative aspects of fibril formation. We used rare events methods to predict separate docking and locking rate constants for the model. The rate constants were embedded in the microkinetic model, with the microkinetic model embedded in a population balance model for “bottom-up” multiscale fibril growth rate predictions. These were compared to “top-down” results using simulation data with the same model and multiscale framework to obtain maximum likelihood estimates of the separate lock and dock rate constants. We used the same procedures to extract separate docking and locking rate constants from experimental fibril growth data. Our multiscale strategy, embedding rate theories, and kinetic models in conservation laws should help to extract docking and locking rate constants from experimental data or long molecular simulations with correct units and without compromising the molecular description.

Molecular dynamics simulations of amyloid-β(16-22) peptide aggregation at air-water interfaces

Oligomers of amyloid-β (Aβ) peptides are known to be related to Alzheimer's disease, and their formation is accelerated at hydrophilic-hydrophobic interfaces, such as the cell membrane surface and air-water interface. Here, we report molecular dynamics simulations of aggregation of Aβ(16-22) peptides at air-water interfaces. First, 100 randomly distributed Aβ(16-22) peptides moved to the interface. The high concentration of peptides then accelerated their aggregation and formation of antiparallel β-sheets. Two layers of oligomers were observed near the interface. In the first layer from the interface, the oligomer with less β-bridges exposed the hydrophobic residues to the air. The second layer consisted of oligomers with more β-bridges that protruded into water. They are more soluble in water because the hydrophobic residues are covered by N- and C-terminal hydrophilic residues that are aligned well along the oligomer edge. These results indicate that amyloid protofibril formation mainly occurs in the second layer.

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.

The role of amyloid oligomers in neurodegenerative pathologies

Many neurodegenerative diseases are rooted in the activities of amyloid-like proteins which possess conformations that spread to healthy proteins. These include Alzheimer's disease (AD), Parkinson's disease (PD), Huntington's disease (HD) and amyotrophic lateral sclerosis (ALS). While their clinical manifestations vary, their protein-level mechanisms are remarkably similar. Aberrant monomeric proteins undergo conformational shifts, facilitating aggregation and formation of solid fibrils. However, there is growing evidence that intermediate oligomeric stages are key drivers of neuronal toxicity. Analysis of protein dynamics is complicated by the fact that nucleation and growth of amyloid-like proteins is not a linear pathway. Feedback within this pathway results in exponential acceleration of aggregation, but activities exerted by oligomers and fibrils can alter cellular interactions and the cellular environment as a whole. The resulting cascade of effects likely contributes to the late onset and accelerating progression of amyloid-like protein disorders and the widespread effects they have on the body. In this review we explore the amyloid-like proteins associated with AD, PD, HD and ALS, as well as the common mechanisms of amyloid-like protein nucleation and aggregation. From this, we identify core elements of pathological progression which have been targeted for therapies, and which may become future therapeutic targets.

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

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

Structural analysis of the β-sheet edge of peptide self-assembly using a model protein

Peptides and proteins self-assemble into β-sheet-rich fibrils, amyloid, which extends its structure by incorporating peptide/protein molecules from solution. At the elongation edge, the peptide/protein molecule binds to the edge of the amyloid β-sheet. Such processes are transient and elusive when observing molecular details by experimental methods. We used a model protein system, peptide self-assembly mimic (PSAM), which mimics an amyloid-like structure within a globular protein by capping both edges of single-layer β sheet (SLB) with certain domains. We constructed a PSAM variant that lacks the capping domain on the C-terminal side to observe the structure of the β-sheet edge of the peptide self-assembly. This variant, which we termed PSAM-edge, proved to be soluble with a monomeric form. Urea-induced unfolding experiments revealed that PSAM-edge displayed two-state cooperative unfolding, indicating the N-terminal capping domain and extended SLB folded as one unit. The crystal structure showed that SLB was almost completely structured except for a few terminal residues. A molecular dynamics simulation results revealed that the SLB structure was retained while the C-terminal four residues fluctuated, which was consistent with the crystal structure. Our findings indicate that SLB is stable even when one side of the β-sheet edge is exposed to a solvent. This stability may prevent the dissociation of the attached peptide from the peptide self-assembly. Because of the scarcity of SLB proteins with exposed β-sheet edges in nature, successful construction of the PSAM-edge expands our understanding of protein folding and design.

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.

Lysozyme Fibrils Alter the Mechanism of Insulin Amyloid Aggregation

Protein aggregation into amyloid fibrils is linked to multiple disorders. The understanding of how natively non-harmful proteins convert to these highly cytotoxic amyloid aggregates is still not sufficient, with new ideas and hypotheses being presented each year. Recently it has been shown that more than one type of protein aggregates may co-exist in the affected tissue of patients suffering from amyloid-related disorders, sparking the idea that amyloid aggregates formed by one protein may induce another protein’s fibrillization. In this work, we examine the effect that lysozyme fibrils have on insulin amyloid aggregation. We show that not only do lysozyme fibrils affect insulin nucleation, but they also alter the mechanism of its aggregation.

Exploring the occurrence of thioflavin-T-positive insulin amyloid aggregation intermediates

The aggregation of proteins is considered to be the main cause of several neurodegenerative diseases. Despite much progress in amyloid research, the process of fibrillization is still not fully understood, which is one of the main reasons why there are still very few effective treatments available. When the aggregation of insulin, a model amyloidogenic protein, is tracked using thioflavin-T (ThT), an amyloid specific dye, there is an anomalous occurrence of double-sigmoidal aggregation kinetics. Such an event is likely related to the formation of ThT-positive intermediates, which may affect the outcome of both aggregation kinetic data, as well as final fibril structure. In this work we explore insulin fibrillization under conditions, where both normal and double-sigmoidal kinetics are observed and show that, despite their dye-binding properties and random occurrence, the ThT-positive intermediates do not significantly alter the overall aggregation process.

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.

Amyloid-Like Peptide Aggregates

Self-assembly of proteins and peptides into the amyloid fold is a widespread phenomenon in the natural world. The structural hallmark of self-assembly into amyloid fibrillar assemblies is the cross-beta motif, which conveys distinct morphological and mechanical properties. The amyloid fibril formation has contrasting results depending on the organism, in the sense that it can bestow an organism with the advantages of mechanical strength and improved functionality or, on the contrary, could give rise to pathological states. In this chapter we review the existing information on amyloid-like peptide aggregates, which could either be derived from protein sequences, but also could be rationally or de novo designed in order to self-assemble into amyloid fibrils under physiological conditions. Moreover, the development of self-assembled fibrillar biomaterials that are tailored for the desired properties towards applications in biomedical or environmental areas is extensively analyzed. We also review computational studies predicting the amyloid propensity of the natural amino acid sequences and the structure of amyloids, as well as designing novel functional amyloid materials.

On the Mechanism of Self-Assembly by a Hydrogel-Forming Peptide

Self-assembling peptide-based hydrogels are a class of tunable soft materials that have been shown to be highly useful for a number of biomedical applications. The dynamic formation of the supramolecular fibrils that compose these materials has heretofore remained poorly characterized. A better understanding of this process would provide important insights into the behavior of these systems and could aid in the rational design of new peptide hydrogels. Here, we report the determination of the microscopic steps that underpin the self-assembly of a hydrogel-forming peptide, SgI_37-49. Using theoretical models of linear polymerization to analyze the kinetic self-assembly data, we show that SgI_37-49 fibril formation is driven by fibril-catalyzed secondary nucleation and that all the microscopic processes involved in SgI_37-49 self-assembly display an enzyme-like saturation behavior. Moreover, this analysis allows us to quantify the rates of the underlying processes at different peptide concentrations and to calculate the time evolution of these reaction rates over the time course of self-assembly. We demonstrate here a new mechanistic approach for the study of self-assembling hydrogel-forming peptides, which is complementary to commonly used materials science characterization techniques.

High-resolution structure of a partially folded insulin aggregation intermediate

Insulin has long been served as a model for protein aggregation, both due to the importance of aggregation in the manufacture of insulin and because the structural biology of insulin has been extensively characterized. Despite intensive study, details about the initial triggers for aggregation have remained elusive at the molecular level. We show here that at acidic pH, the aggregation of insulin is likely initiated by a partially folded monomeric intermediate. High-resolution structures of the partially folded intermediate show that it is coarsely similar to the initial monomeric structure but differs in subtle details-the A chain helices on the receptor interface are more disordered and the B chain helix is displaced from the C-terminal A chain helix when compared to the stable monomer. The result of these movements is the creation of a hydrophobic cavity in the center of the protein that may serve as nucleation site for oligomer formation. Knowledge of this transition may aid in the engineering of insulin variants that retain the favorable pharamacokinetic properties of monomeric insulin but are more resistant to aggregation.

SAMPL7 TrimerTrip host-guest binding poses and binding affinities from spherical-coordinates-biased simulations

Host-guest binding remains a major challenge in modern computational modelling. The newest 7th statistical assessment of the modeling of proteins and ligands (SAMPL) challenge contains a new series of host-guest systems. The TrimerTrip host binds to 16 structurally diverse guests. Previously, we have successfully employed the spherical coordinates as the collective variables coupled with the enhanced sampling technique metadynamics to enhance the sampling of the binding/unbinding event, search for possible binding poses and calculate the binding affinities in all three host-guest binding cases of the 6th SAMPL challenge. In this work, we report a retrospective study on the TrimerTrip host-guest systems by employing the same protocol to investigate the TrimerTrip host in the SAMPL7 challenge. As no binding pose is provided by the SAMPL7 host, our simulations initiate from randomly selected configurations and are proceeded long enough to obtain converged free energy estimates and search for possible binding poses. The calculated binding affinities are in good agreement with the experimental reference, and the obtained binding poses serve as a nice starting point for end-point or alchemical free energy calculations. Note that as the work is performed after the close of the SAMPL7 challenge, we do not participate in the challenge and the results are not formally submitted to the SAMPL7 challenge.

N-Terminal Modified Aβ Variants Enable Modulations to the Structures and Cytotoxicity Levels of Wild-Type Aβ Fibrils through Cross-Seeding

Post-translational modifications (PTMs) of β-amyloid (Aβ) peptides are considered as triggering factors in sporadic Alzheimer's disease. However, studies to show the influence of pre-existing PTM-Aβ fibrils on wild-type Aβ peptides, which directly mimic the triggering scenarios, are rare. Here we show that three types of pathologically relevant PTM-Aβ variants with modifications in a particular segment (from D7 to V12) of the primary sequence lead to distinct impacts on the fibrillization of wild-type Aβ peptides. In general, the triggering effects are observed through cross-seeding between the PTM-Aβ seeds and wild-type peptides, which consequently induce modulations in the resultant wild-type fibril structures and elevations in the fibrillar cytotoxicity levels. Modifications with a similar chemical nature, such as the S8-phosphorylation and Y10-nitration, both of which introduce additional side-chain negative charges, show comparable structural-modulation and cytotoxicity-elevation effects. The results imply the biological influences of PTM-Aβ variants on the formation of amyloid deposits through cross-seeded fibrillization.

Redox-Dependent Copper Ion Modulation of Amyloid-β (1-42) Aggregation In Vitro

Plaque deposits composed of amyloid-β (Aβ) fibrils are pathological hallmarks of Alzheimer's disease (AD). Although copper ion dyshomeostasis is apparent in AD brains and copper ions are found co-deposited with Aβ peptides in patients' plaques, the molecular effects of copper ion interactions and redox-state dependence on Aβ aggregation remain elusive. By combining biophysical and theoretical approaches, we here show that Cu2+ (oxidized) and Cu+ (reduced) ions have opposite effects on the assembly kinetics of recombinant Aβ(1-42) into amyloid fibrils in vitro. Cu2+ inhibits both the unseeded and seeded aggregation of Aβ(1-42) at pH 8.0. Using mathematical models to fit the kinetic data, we find that Cu2+ prevents fibril elongation. The Cu2+-mediated inhibition of Aβ aggregation shows the largest effect around pH 6.0 but is lost at pH 5.0, which corresponds to the pH in lysosomes. In contrast to Cu2+, Cu+ ion binding mildly catalyzes the Aβ(1-42) aggregation via a mechanism that accelerates primary nucleation, possibly via the formation of Cu+-bridged Aβ(1-42) dimers. Taken together, our study emphasizes redox-dependent copper ion effects on Aβ(1-42) aggregation and thereby provides further knowledge of putative copper-dependent mechanisms resulting in AD.

Interactions between Curcumin Derivatives and Amyloid-β Fibrils: Insights from Molecular Dynamics Simulations

The aggregation of amyloid-β (Aβ) peptides into senile plaques is a hallmark of Alzheimer's disease (AD) and is hypothesized to be the primary cause of AD related neurodegeneration. Previous studies have shown the ability of curcumin to both inhibit the aggregation of Aβ peptides into oligomers or fibrils and reduce amyloids in vivo. Despite the promise of curcumin and its derivatives to serve as diagnostic, preventative, and potentially therapeutic AD molecules, the mechanism by which curcumin and its derivatives bind to and inhibit Aβ fibrils' formation remains elusive. Here, we investigated curcumin and a set of curcumin derivatives in complex with a hexamer peptide model of the Aβ1-42 fibril using nearly exhaustive docking, followed by multi-ns molecular dynamics simulations, to provide atomistic-detail insights into the molecules' binding and inhibitory properties. In the vast majority of the simulations, curcumin and its derivatives remain firmly bound in complex with the fibril through primarily three different principle binding modes, in which the molecules interact with residue domain 17LVFFA21, in line with previous experiments. In a small subset of these simulations, the molecules partly dissociate the outermost peptide of the Aβ1-42 fibril by disrupting β-sheets within the residue domain 12VHHQKLVFF20. A comparison between binding modes leading or not leading to partial dissociation of the outermost peptide suggests that the latter is attributed to a few subtle key structural and energetic interaction-based differences. Interestingly, partial dissociation appears to be either an outcome of high affinity interactions or a cause leading to high affinity interactions between the molecules and the fibril, which could partly serve as a compensation for the energy loss in the fibril due to partial dissociation. In conjunction with this, we suggest a potential inhibition mechanism of Αβ1-42 aggregation by the molecules, where the partially dissociated 16KLVFF20 domain of the outermost peptide could either remain unstructured or wrap around to form intramolecular interactions with the same peptide's 29GAIIG33 domain, while the molecules could additionally act as a patch against the external edge of the second outermost peptide's 16KLVFF20 domain. Thereby, individually or concurrently, these could prohibit fibril elongation.

Folding dynamics of Aβ42 monomer at pH 4.0–7.5 with and without physiological salt conditions – does the β1 or β2 region fold first?

The seeding region of Aβ42 monomer is jointly affected by the solution acidity, ionic distribution of the salt, and charged residues.

Structural and Binding Properties on Aβ Mature Fibrils Due to the Histidine Tautomeric Effect

Many studies have focused on histidine behaviors in misfolding diseases. However, histidine behaviors on mature fibrils are still unknown. In the current study, we investigated mature fibrils with various histidine states to understand the structural properties of the histidine tautomeric effect on mature fibrils. Our results show that substituting chain 1 with different histidine states affects Aβ structural properties in A2, D7-G9, H14-Q15, S26-N27, and G33-G37 regions. The binding free energies with substituted fibrils were influenced not only along the axial direction, but also between duplex fibrils. Our results suggest that substituted (εδδ) preferentially disturbed the stability among the current mature fibrils. Further, H-bonded network differences indicate that twisted morphologies in mature fibrils are derived from the position and orientation of the imidazole ring in histidines. Our current study helps to elucidate histidine behaviors on mature fibrils, which will present opportunities to understand the misfolding mechanisms.

The Formation Mechanism of Segmented Ring-Shaped Aβ Oligomers and Protofibrils

A clear understanding of amyloid formation with diverse morphologies is critical to overcoming the fatal disease amyloidosis. Studies have revealed that monomer concentration is a crucial factor for determining amyloid morphologies, such as protofibrils, annular, or spherical oligomers. However, gaining a complete understanding of the mechanism of formation of the various amyloid morphologies has been limited by the lack of experimental devices and insufficient knowledge. In this study, we demonstrate that the monomer concentration is an essential factor in determining the morphology of beta-amyloid (Aβ) oligomers or protofibrils. By computational and experimental approaches, we investigated the strategies for structural stabilization of amyloid protein, the morphological changes, and amyloid aggregation. In particular, we found unprecedented conformations, e.g., single bent oligomers and segmented ring-shaped protofibrils, the formation of which was explained by the computational analysis. Our findings provide insight into the structural features of amyloid molecules formed at low concentrations of monomer, which will help determine the clinical targets (in therapy) to effectively inhibit amyloid formation in the early stages of the amyloid growth phase.

Early mechanisms of amyloid fibril nucleation in model and disease-related proteins

Protein amyloid aggregation is a hallmark in neuropathologies and other diseases of tremendous impact such as Alzheimer's or Parkinson's diseases. During the last decade, it has become increasingly evident that neuronal death is mainly induced by proteinaceous oligomers rather than the mature amyloid fibrils. Therefore, the earliest molecular events occurring during the amyloid aggregation cascade represent a growing interest of study. Important breakthroughs have been achieved using experimental data from different proteins, used as models, as well as systems related to diseases. Here, we summarize the structural properties of amyloid oligomeric and fibrillar aggregates and review the recent advances on how biophysical techniques can be combined with quantitative kinetic analysis and theoretical models to study the detailed mechanism of oligomer formation and nucleation of fibrils. These insights into the mechanism of early oligomerization and amyloid nucleation are of relevant interest in drug discovery and in the design of preventive strategies against neurodegenerative diseases.

Interference with Amyloid-β Nucleation by Transient Ligand Interaction

Amyloid-β peptide (Aβ) is an intrinsically disordered protein (IDP) associated with Alzheimer's disease. The structural flexibility and aggregation propensity of Aβ pose major challenges for elucidating the interaction between Aβ monomers and ligands. All-D-peptides consisting solely of D-enantiomeric amino acid residues are interesting drug candidates that combine high binding specificity with high metabolic stability. Here we characterized the interaction between the 12-residue all-D-peptide D3 and Aβ42 monomers, and how the interaction influences Aβ42 aggregation. We demonstrate for the first time that D3 binds to Aβ42 monomers with submicromolar affinities. These two highly unstructured molecules are able to form complexes with 1:1 and other stoichiometries. Further, D3 at substoichiometric concentrations effectively slows down the β-sheet formation and Aβ42 fibrillation by modulating the nucleation process. The study provides new insights into the molecular mechanism of how D3 affects Aβ assemblies and contributes to our knowledge on the interaction between two IDPs.

Molecular insights into the surface-catalyzed secondary nucleation of amyloid-β40 (Aβ40) by the peptide fragment Aβ16-22

Understanding the structural mechanism by which proteins and peptides aggregate is crucial, given the role of fibrillar aggregates in debilitating amyloid diseases and bioinspired materials. Yet, this is a major challenge as the assembly involves multiple heterogeneous and transient intermediates. Here, we analyze the co-aggregation of Aβ₄₀ and Aβ_16-22, two widely studied peptide fragments of Aβ₄₂ implicated in Alzheimer's disease. We demonstrate that Aβ_16-22 increases the aggregation rate of Aβ₄₀ through a surface-catalyzed secondary nucleation mechanism. Discontinuous molecular dynamics simulations allowed aggregation to be tracked from the initial random coil monomer to the catalysis of nucleation on the fibril surface. Together, the results provide insight into how dynamic interactions between Aβ₄₀ monomers/oligomers on the surface of preformed Aβ_16-22 fibrils nucleate Aβ₄₀ amyloid assembly. This new understanding may facilitate development of surfaces designed to enhance or suppress secondary nucleation and hence to control the rates and products of fibril assembly.

Understanding Aldose Reductase-Inhibitors interactions with free energy simulation

Aldose Reductase (AR) reduces a variety of substrates, such as aldehydes, aldoses and corticosteroids. It is the first and rate-limiting enzyme of the polyol pathway and is an important target enzyme for diabetes-associated complications, including retinopathy, neuropathy, and nephropathy. Inhibitors targeting this enzyme are structurally different and some of them have side effects. In existing publications, computational techniques are applied to investigate the binding affinities of existing inhibitors and predicting the affinities of newly designed ligands. However, these calculations only employ coarse and approximated methods such as docking and MM/PBSA. Brute force simulations are employed to study the dynamics of the system but no converged statistics is obtained. As a result, these computations provide results not consistent with experimental values and large discrepancies exist. In the current work, we employ the enhanced sampling technique of alchemical free energy simulation to calculate the binding affinities of several ligands targeting AR. The statistical error is defined with care and the correlation in the time-series data is fully considered. The statistically optimal estimators are used to extract the free energy estimates and the predicted results are in agreement with the experimental values. Less computationally demanding end-point free energy methods are also performed to compare their efficiency with the alchemical methods. As is expected, the end-point methods are of less accuracy and reliability compared with the alchemical free energy methods. The decomposition of the free energy difference in each alchemical transformation into the enthalpic and entropic components gives further insights on the thermodynamics. The enthalpy-entropy compensation is observed in this case. As the structural data obtained from experiments are only snapshots and more details are needed to understand the dynamics of the protein-ligand system, the conformational ensemble is analyzed. We identify important residues involved in the protein-ligand binding case and short-lived interactions formed due to fluctuations in the conformational ensemble. The current work shed light on the atomic detailed understanding of the dynamics of AR-inhibitors interactions.

Probing transient non-native states in amyloid beta fiber elongation by NMR

Using NMR to probe transient binding of Aβ1-40 monomers to fibers, we find partially bound conformations with the highest degree of interaction near F19-K28 and a lesser degree of interaction near the C-terminus (L34-G37). This represents a shift away from the KLVFFA recognition sequence (residues 16-21) currently used for inhibitor design.

Amyloid Fibril Design: Limiting Structural Polymorphism in Alzheimer's Aβ Protofilaments

Nanoscale fibrils formed by amyloid peptides have a polymorphic character, adopting several types of molecular structures in similar growth conditions. As shown by experimental (e.g., solid-state NMR) and computational studies, amyloid fibril polymorphism hinders both the structural characterization of Alzheimer's Aβ amyloid protofilaments and fibrils at a molecular level, as well as the possible applications (e.g., development of drugs or biomarkers) that rely on similar, controlled molecular arrangements of the Aβ peptides in amyloid fibril structures. We have explored the use of several contact potentials for the efficient identification of minimal sequence mutations that could enhance the stability of specific fibril structures while simultaneously destabilizing competing topologies, controlling thus the amount of structural polymorphism in a rational way. We found that different types of contact potentials, while having only partial accuracy on their own, lead to similar results regarding ranking the compatibility of wild-type (WT) and mutated amyloid sequences with different fibril morphologies. This approach allows exhaustive screening and assessment of possible mutations and the identification of minimal consensus mutations that could stabilize fibrils with the desired topology at the expense of other topology types, a prediction that is further validated using atomistic molecular dynamics with explicit water molecules. We apply this two-step multiscale (i.e., residue and atomistic-level) approach to predict and validate mutations that could bias either parallel or antiparallel packing in the core Alzheimer's Aβ_9-40 amyloid fibril models based on solid-state NMR experiments. Besides shedding new light on the molecular origins of structural polymorphism in WT Aβ fibrils, our study could also lead to efficient tools for assisting future experimental approaches for amyloid fibril determination, and for the development of biomarkers or drugs aimed at interfering with the stability of amyloid fibrils, as well as for the future design of amyloid fibrils with a controlled (e.g., reduced) level of structural polymorphism.

Mechanistic and Structural Origins of the Asymmetric Barrier to Prion-like Cross-Seeding between Tau-3R and Tau-4R

The spread and deposition of infectious fibrillar protein aggregates in the brain via a prion-like mechanism is a critical component in the patho-physiology of various neurodegenerative diseases, including the tauopathies. In tauopathies, two isoforms of tau, containing three and four microtubule binding repeats, are found to aggregate, and the type of isoform present in aggregates determines the type of tauopathy. Cross-seeding between the two tau isoforms is limited by an asymmetric barrier similar to the species barrier that restricts prion transmission across species, whose origin has remained unclear. In this study, the growth of the tau fibrils is shown to be describable by a two-step Michaelis-Menten-like model. Delineation of the mechanism as a Michaelis-Menten-like mechanism has enabled a quantitative understanding of the asymmetric seeding barrier that exists between two isoforms of tau, tau-K18 and tau-K19 (which differ in containing four and three microtubule binding repeats, respectively), wherein tau-K18 fibrils cannot seed tau-K19 monomer. Furthermore, high-resolution structural analysis of the two isoforms shows that the structural core is more ordered in tau-K19 than in tau-K18. Hence, the current work provides kinetic and structural rationales for asymmetric seeding barriers in general and for the two tau isoforms in particular.

Nonproductive Binding Modes as a Prominent Feature of Aβ40 Fiber Elongation: Insights from Molecular Dynamics Simulation

The formation of amyloid fibers has been implicated in a number of neurodegenerative diseases. The growth of amyloid fibers is strongly thermodynamically favorable, but kinetic traps exist where the incoming monomer binds in an incompatible conformation that blocks further elongation. Unfortunately, this process is difficult to follow experimentally at the atomic level. It is also too complex to simulate in full detail and to date has been explored either through coarse-grained simulations, which may miss many important interactions, or full atomic simulations, in which the incoming peptide is constrained to be near the ideal fiber geometry. Here we use an alternate approach starting from a docked complex in which the monomer is from an experimental NMR structure of one of the major conformations in the unbound ensemble, a largely unstructured peptide with the central hydrophobic region in a 3₁₀ helix. A 1000 ns full atomic simulation in explicit solvent shows the formation of a metastable intermediate by sequential, concerted movements of both the fiber and the monomer. A Markov state model shows that the unfolded monomer is trapped at the end of the fiber in a set of interconverting antiparallel β-hairpin conformations. The simulation here may serve as a model for the binding of other non-β-sheet conformations to amyloid fibers.