Dynamics of locking of peptides onto growing amyloid fibrils

Multiscale modeling of protofilament structures: A case study on insulin amyloid aggregates

Under certain conditions, proteins may undergo misfolding and form long insoluble aggregates called amyloid fibrils. The presence of these aggregates is often associated with various diseases. The molecular mechanisms governing the aggregation process are yet to be fully understood. The self-assembly of amyloid protofilaments occurs over extended time frames, making the simulation of such events problematic. In this work, we describe a pipeline for multiscale modeling protofilament structures. In the first stage, the self-assembly of short fibrillar oligomers occurs during coarse-grained docking simulations of multiple copies of aggregating peptides. Subsequently, symmetry criteria are used to select the highest-ranked oligomer structures. Selected models are then reconstructed to an all-atom representation and used for the assembly of longer protofilaments. Models are optimized using molecular dynamics. Final structures are selected using various scoring protocols. We evaluated this modeling procedure through the test prediction of insulin amyloid protofilaments whose experimental structures have been published recently. The resulting insulin protofilament models closely resemble the experimental structures. This work provides a proof of concept for the proposed modeling procedure aiming to predict amyloid protofilament structures that exhibit in-register and parallel arrangement of β-sheets based solely on the amino acid sequence of aggregating peptides.

The energy landscape of Aβ42: a funnel to disorder for the monomer becomes a folding funnel for self-assembly

The aggregation of amyloid-β (Aβ) peptides, particularly Aβ1-42, plays a key role in Alzheimer's disease pathogenesis. In this study, we investigate how dimerisation transforms the free energy surface (FES) of the Aβ1-42 monomer when it interacts with another Aβ1-42 peptide. We find that the monomer FES is a structurally inverted funnel with a disordered state at the global minimum. However, in the presence of a second Aβ1-42 peptide, the landscape becomes a folding funnel, leading to a β-hairpin state. Using first passage time analysis, we analyse the pathway for the transition from disordered to the β-hairpin state, which highlights the initial formation of a D23-K28 salt bridge as the driving force, together with hydrophobic contacts.

Molecular Integrative Study on Inhibitory Effects of Pentapeptides on Polymerization and Cell Toxicity of Amyloid-β Peptide (1-42)

Alzheimer's Disease (AD) is a multifaceted neurodegenerative disease predominantly defined by the extracellular accumulation of amyloid-β (Aβ) peptide. In light of this, in the past decade, several clinical approaches have been used aiming at developing peptides for therapeutic use in AD. The use of cationic arginine-rich peptides (CARPs) in targeting protein aggregations has been on the rise. Also, the process of peptide development employing computational approaches has attracted a lot of attention recently. Using a structure database containing pentapeptides made from 20 L-α amino acids, we employed molecular docking to sort pentapeptides that can bind to Aβ42, then performed molecular dynamics (MD) analyses, including analysis of the binding stability, interaction energy, and binding free energy to screen ligands. Transmission electron microscopy (TEM), circular dichroism (CD), thioflavin T (ThT) fluorescence detection of Aβ42 polymerization, MTT (3-[4,5-dimethylthiazol-2-yl]-2,5 diphenyl tetrazolium bromide) assay, and the flow cytometry of reactive oxygen species (ROS) were carried out to evaluate the influence of pentapeptides on the aggregation and cell toxicity of Aβ42. Two pentapeptides (TRRRR and ARRGR) were found to have strong effects on inhibiting the aggregation of Aβ42 and reducing the toxicity of Aβ42 secreted by SH-SY5Y cells, including cell death, reactive oxygen species (ROS) production, and apoptosis.

Insights into the mechanism of peptide fibril growth on gold surface

Understanding the formation of β-fibrils over the gold surface is of paramount interest in nano-bio-medicinal Chemistry. The intricate mechanism of self-assembly of neurofibrillogenic peptides and their growth over the gold surface remains elusive, as experiments are limited in unveiling the microscopic dynamic details, in particular, at the early stage of the peptide aggregation. In this work, we carried out equilibrium molecular dynamics and enhanced sampling simulations to elucidate the underlying mechanism of the growth of an amyloid-forming sequence of tau fragments over the gold surface. Our results disclose that the collective intermolecular interactions between the peptide chains and peptides with the gold surface facilitate the peptide adsorption, followed by integration, finally leading to the fibril formation.

Mechanistic Modeling of Amyloid Oligomer and Protofibril Formation in Bovine Insulin

Early phase of amyloid formation, where prefibrillar aggregates such as oligomers and protofibrils are often observed, is crucial for understanding pathogenesis. However, the detailed mechanisms of their formation have been difficult to elucidate because they tend to form transiently and heterogeneously. Here, we found that bovine insulin protofibril formation proceeds in a monodisperse manner, which allowed us to characterize the detailed early aggregation process by light scattering in combination with thioflavin T fluorescence and Fourier transform infrared spectroscopy. The protofibril formation was specific to bovine insulin, whereas no significant aggregation was observed in human insulin. The kinetic analysis combining static and dynamic light scattering data revealed that the protofibril formation process in bovine insulin can be divided into two steps based on fractal dimension. When modeling the experimental data based on Smoluchowski aggregation kinetics, an aggregation scheme consisting of initial fractal aggregation forming spherical oligomers and their subsequent end-to-end association forming protofibrils was clarified. Furthermore, the analysis of temperature and salt concentration dependencies showed that the end-to-end association is the rate-limiting step, involving dehydration. The established model for protofibril formation, wherein oligomers are incorporated as a precursor, provides insight into the molecular mechanism by which protein molecules assemble during the early stage of amyloid formation.

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.

A Molecular Integrative Study on the Inhibitory Effects of WRR and ERW on Amyloid β Peptide (1-42) Polymerization and Cell Toxicity

Alzheimer's disease (AD) is a neurodegenerative disease and the main pathological characteristic of AD is the deposition of Aβ42 in the brain. Inhibition of Aβ42 polymerization is one of the important research directions. Due to the pathological complexity of Alzheimer's disease, studies on Aβ42 polymerization inhibitors have not made significant progress worldwide. Using an independently constructed structure database of oligopeptides, in this study, molecular docking, umbrella sampling analysis of free energy, ThT fluorescence detection of Aβ42 polymerization, transmission electron microscopy, and flow cytometry detection of reactive oxygen species (ROS) and apoptosis were performed to screen tripeptides and pentapeptides that inhibit polymerization. It was found that two tripeptides, i.e., WRR and ERW, bind stably to the core of Aβ42 polymerization in the molecular dynamics analysis, and they significantly inhibited the aggregation of Aβ42 and reduced their cell toxicity in vitro.

Flanking regions, amyloid cores, and polymorphism: the potential interplay underlying structural diversity

The β-sheet-rich amyloid core is the defining feature of protein aggregates associated with neurodegenerative disorders. Recent investigations have revealed that there exist multiple examples of the same protein, with the same sequence, forming a variety of amyloid cores with distinct structural characteristics. These structural variants, termed as polymorphs, are hypothesized to influence the pathological profile and the progression of different neurodegenerative diseases, giving rise to unique phenotypic differences. Thus, identifying the origin and properties of these structural variants remain a focus of studies, as a preliminary step in the development of therapeutic strategies. Here, we review the potential role of the flanking regions of amyloid cores in inducing polymorphism. These regions, adjacent to the amyloid cores, show a preponderance for being structurally disordered, imbuing them with functional promiscuity. The dynamic nature of the flanking regions can then manifest in the form of conformational polymorphism of the aggregates. We take a closer look at the sequences flanking the amyloid cores, followed by a review of the polymorphic aggregates of the well-characterized proteins amyloid-β, α-synuclein, Tau, and TDP-43. We also consider different factors that can potentially influence aggregate structure and how these regions can be viewed as novel targets for therapeutic strategies by utilizing their unique structural properties.

Antibody binding modulates the dynamics of the membrane-bound prion protein

Misfolding of the cellular prion protein (PrPC) is associated with lethal neurodegeneration. PrPC consists of a flexible tail (residues 23-123) and a globular domain (residues 124-231) whose C-terminal end is anchored to the cell membrane. The neurotoxic antibody POM1 and the innocuous antibody POM6 recognize the globular domain. Experimental evidence indicates that POM1 binding to PrPC emulates the influence on PrPC of the misfolded prion protein (PrPSc) while the binding of POM6 has the opposite biological response. Little is known about the potential interactions between flexible tail, globular domain, and the membrane. Here, we used atomistic simulations to investigate how these interactions are modulated by the binding of the Fab fragments of POM1 and POM6 to PrPC and by interstitial sequence truncations to the flexible tail. The simulations show that the binding of the antibodies restricts the range of orientations of the globular domain with respect to the membrane and decreases the distance between tail and membrane. Five of the six sequence truncations influence only marginally this distance and the contact patterns between tail and globular domain. The only exception is a truncation coupled to a charge inversion mutation of four N-terminal residues, which increases the distance of the flexible tail from the membrane. The interactions of the flexible tail and globular domain are modulated differently by the two antibodies.

Exploring Heterogeneous Dynamical Environment around an Ensemble of Aβ42 Peptide Monomer Conformations

Exploring the conformational properties of amyloid β (Aβ) peptides and the role of solvent (water) in guiding the dynamical environment at their interfaces is crucial for microscopic understanding of Aβ misfolding, which is involved in causing the most common neurodegenerative disorder, i.e., Alzheimer's disease. While numerous studies in the past have emphasized examining the conformational states of Aβ peptides, the role of water has not received much attention. Here, we have performed all-atom molecular dynamics simulations of several full-length Aβ₄₂ peptide monomers with different initial configurations. Our efforts are directed toward probing the origin of the heterogeneous dynamics of water around various segments of the Aβ peptide, identified as the two terminal segments (N-term and C-term) and the two hydrophobic segments (hp1 and hp2), along with the central turn region interconnecting hp1 and hp2. Our results revealed that water hydrating hp1, hp2, and turn (nonterminal segments) and C-term segments exhibit nonuniformly restricted translational as well as rotational motions. The degree of such restriction has been found to be correlated with the hydrogen bond relaxation time scales at the interface. Importantly, it is revealed that the water molecules around hp1 and, to some extent, around hp2, form relatively rigid hydration layers, compared to that around the other segments. Such rigid hydration layers arise due to relatively more solid-like caging motions resulting in relatively lesser hydration entropy. As hp1 and hp2 have been demonstrated to play a central role in Aβ aggregation, we believe that distinct water dynamics in the vicinity of these two segments, as outlined in this study, can provide vital information in understanding the early stages of the onset of the aggregation process of such peptides at higher concentration that can further aid toward advances in AD therapeutics.

Scan-Find-Scan-Model: Discrete Site-Targeted Suppressor Design Strategy for Amyloid-β

Alzheimer's disease is undoubtedly the most well-studied neurodegenerative disease. Consequently, the amyloid-β (Aβ) protein ranks at the top in terms of getting attention from the scientific community for structural property-based characterization. Even after decades of extensive research, there is existing volatility in terms of understanding and hence the effective tackling procedures against the disease that arises due to the lack of knowledge of both specific target- and site-specific drugs. Here, we develop a multidimensional approach based on the characterization of the common static-dynamic-thermodynamic trait of the monomeric protein, which efficiently identifies a small target sequence that contains an inherent tendency to misfold and consequently aggregate. The robustness of the identification of the target sequence comes with an abundance of a priori knowledge about the length and sequence of the target and hence guides toward effective designing of the target-specific drug with a very low probability of bottleneck and failure. Based on the target sequence information, we further identified a specific mutant that showed the maximum potential to act as a destabilizer of the monomeric protein as well as enormous success as an aggregation suppressor. We eventually tested the drug efficacy by estimating the extent of modulation of binding affinity existing within the fibrillar form of the Aβ protein due to a single-point mutation and hence provided a proof of concept of the entire protocol.

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.

Single-molecule Observation of Self-Propagating Amyloid Fibrils

The assembly of misfolded proteins into amyloid fibrils is associated with amyloidosis, including neurodegenerative diseases, such as Alzheimer’s, Parkinson’s, and prion diseases. The self-propagation of amyloid fibrils is widely observed in the aggregation pathways of numerous amyloidogenic proteins. This propensity with plasticity in primary nucleation allows amyloid fibril polymorphism, which is correlated with the pathology/phenotypes of patients. Because the interference with the nucleation and replication processes of amyloid fibrils can alter the amyloid structure and the outcome of the disease, these processes can be a target for developing clinical drugs. Single-molecule observation of amyloid fibril replication can be an experimental system to provide the kinetic parameters for simulation studies and confirm the effect of clinical drugs. Here, we review single-molecule observation of the amyloid fibril replication process using fluorescence microscopy and time-lapse atomic force microscopy, including high-speed atomic force microscopy. We discussed the amyloid fibril replication process and combined single-molecule observation results with molecular dynamics simulations.

Mini Abstract Structural dynamics in amyloid aggregation is related with various Alzheimer’s and Parkinson’s disease symptoms. Single-molecule observation using high-speed atomic force microscopy can directly visualize the structural dynamics of individual amyloid aggregate assemblies. Here, we review historical and recent studies of single-molecule observation of amyloid aggregation with supportive molecular dynamics simulation.

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).

Molecular dynamics study of water channels in natural and synthetic amyloid-β fibrils

We compared all-atom explicit solvent molecular dynamics simulations of three types of Aβ(1-40) fibrils: brain-seeded fibrils (2M4J, with a threefold axial symmetry) and the other two, all-synthetic fibril polymorphs (2LMN and 2LMP, made under different fibrillization conditions). Fibril models were constructed using either a finite or an infinite number of layers made using periodic images. These studies yielded four conclusions. First, finite fibrils tend to unravel in a manner reminiscent of fibril dissolution, while infinite fibrils were more stable during simulations. Second, salt bridges in these fibrils remained stable in those fibrils that contained them initially, and those without salt bridges did not develop them over the time course of the simulations. Third, all fibrils tended to develop a "stagger" or register shift of β-strands along the fibril axis. Fourth and most importantly, the brain-seeded, 2M4J, infinite fibrils allowed bidirectional transport of water in and out of the central longitudinal core of the fibril by rapidly developing gaps at the fibril vertices. 2LMP fibrils also showed this behavior, although to a lesser extent. The diffusion of water molecules in the fibril core region involved two dynamical states: a localized state and directed diffusion in the presence of obstacles. These observations provided support for the hypothesis that Aβ fibrils could act as nanotubes. At least some Aβ oligomers resembled fibrils structurally in having parallel, in-register β-sheets and a sheet-turn-sheet motif. Thus, our findings could have implications for Aβ cytotoxicity, which may occur through the ability of oligomers to form abnormal water and ion channels in cell membranes.

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.

The role of water in the primary nucleation of protein amyloid aggregation

The understanding of the complex conformational landscape of amyloid aggregation and its modulation by relevant physicochemical and cellular factors is a prerequisite for elucidating some of the molecular basis of pathology in amyloid related diseases, and for developing and evaluating effective disease-specific therapeutics to reduce or eliminate the underlying sources of toxicity in these diseases. Interactions of proteins with solvating water have been long considered to be fundamental in mediating their function and folding; however, the relevance of water in the process of protein amyloid aggregation has been largely overlooked. Here, we provide a perspective on the role water plays in triggering primary amyloid nucleation of intrinsically disordered proteins (IDPs) based on recent experimental evidences. The initiation of amyloid aggregation likely results from the synergistic effect between both protein intermolecular interactions and the properties of the water hydration layer of the protein surface. While the self-assembly of both hydrophobic and hydrophilic IDPs would be thermodynamically favoured due to large water entropy contributions, large desolvation energy barriers are expected, particularly for the nucleation of hydrophilic IDPs. Under highly hydrating conditions, primary nucleation is slow, being facilitated by the presence of nucleation-active surfaces (heterogeneous nucleation). Under conditions of poor water activity, such as those found in the interior of protein droplets generated by liquid-liquid phase separation, however, the desolvation energy barrier is significantly reduced, and nucleation can occur very rapidly in the bulk of the solution (homogeneous nucleation), giving rise to structurally distinct amyloid polymorphs. Water, therefore, plays a key role in modulating the transition free energy of amyloid nucleation, thus governing the initiation of the process, and dictating the type of preferred primary nucleation and the type of amyloid polymorph generated, which could vary depending on the particular microenvironment that the protein molecules encounter in the cell.

An explicitly designed paratope of amyloid-β prevents neuronal apoptosis in vitro and hippocampal damage in rat brain

Synthetic antibodies hold great promise in combating diseases, diagnosis, and a wide range of biomedical applications. However, designing a therapeutically amenable, synthetic antibody that can arrest the aggregation of amyloid-β (Aβ) remains challenging. Here, we report a flexible, hairpin-like synthetic paratope (SP1, ∼2 kDa), which prevents the aggregation of Aβ monomers and reverses the preformed amyloid fibril to a non-toxic species. Structural and biophysical studies further allowed dissecting the mode and affinity of molecular recognition events between SP1 and Aβ. Subsequently, SP1 reduces Aβ-induced neurotoxicity, neuronal apoptosis, and ROS-mediated oxidative damage in human neuroblastoma cells (SH-SY5Y). The non-toxic nature of SP1 and its ability to ameliorate hippocampal neurodegeneration in a rat model of AD demonstrate its therapeutic potential. This paratope engineering module could readily implement discoveries of cost-effective molecular probes to nurture the basic principles of protein misfolding, thus combating related diseases.

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.

The Dynamism of Intrinsically Disordered Proteins: Binding-Induced Folding, Amyloid Formation, and Phase Separation

Intrinsically disordered proteins (IDPs) or natively unfolded proteins do not undergo autonomous folding into a well-defined 3-D structure and challenge the conventional structure-function paradigm. They are involved in a multitude of critical physiological functions by adopting various structural states via order-to-disorder transitions or by maintaining their disordered characteristics in functional complexes. In recent times, there has been a burgeoning interest in the investigation of intriguing behavior of IDPs using highly multidisciplinary and complementary approaches due to the pivotal role of this unique class of protein chameleons in physiology and disease. Over the past decade or so, our laboratory has been actively investigating the unique physicochemical properties of this class of highly dynamic, flexible, rapidly interconverting proteins. We have utilized a diverse array of existing and emerging tools involving steady-state and time-resolved fluorescence, Raman spectroscopy, circular dichroism, light scattering, fluorescence microscopy, and atomic force microscopy coupled with site-directed mutagenesis and other biochemical and biophysical tools to study a variety of interesting and important aspects of IDPs. In this Feature Article, I describe our work on the conformational characteristics, solvation dynamics, binding-induced folding, amyloid formation, and liquid-liquid phase separation of a number of amyloidogenic IDPs. A series of these studies described here captures the role of conformational plasticity and dynamics in directing binding, folding, assembly, aggregation, and phase transitions implicated in physiology and pathology.

Structural Characteristics of Monomeric Aβ42 on Fibril in the Early Stage of Secondary Nucleation Process

Amyloid-β (Aβ) aggregates are believed to be one of the main causes of Alzheimer's disease. Aβ peptides form fibrils having cross β-sheet structures mainly through primary nucleation, secondary nucleation, and elongation. In particular, self-catalyzed secondary nucleation is of great interest. Here, we investigate the adsorption of Aβ42 peptides to the Aβ42 fibril to reveal a role of adsorption as a part of secondary nucleation. We performed extensive molecular dynamics simulations based on replica exchange with solute tempering 2 (REST2) to two systems: a monomeric Aβ42 in solution and a complex of an Aβ42 peptide and Aβ42 fibril. Results of our simulations show that the Aβ42 monomer is extended on the fibril. Furthermore, we find that the hairpin structure of the Aβ42 monomer decreases but the helix structure increases by adsorption to the fibril surface. These structural changes are preferable for forming fibril-like aggregates, suggesting that the fibril surface serves as a catalyst in the secondary nucleation process. In addition, the stabilization of the helix structure of the Aβ42 monomer on the fibril indicates that the strategy of a secondary nucleation inhibitor design for Aβ40 can also be used for Aβ42.

Clustering and Fibril Formation during GNNQQNY Aggregation: A Molecular Dynamics Study

The precise kinetic pathways of peptide clustering and fibril formation are not fully understood. Here we study the initial clustering kinetics and transient cluster morphologies during aggregation of the heptapeptide fragment GNNQQNY from the yeast prion protein Sup35. We use a mid-resolution coarse-grained molecular dynamics model of Bereau and Deserno to explore the aggregation pathways from the initial random distribution of free monomers to the formation of large clusters. By increasing the system size to 72 peptides we could follow directly the molecular events leading to the formation of stable fibril-like structures. To quantify those structures we developed a new cluster helicity parameter. We found that the formation of fibril-like structures is a cooperative processes that requires a critical number of monomers, M⋆≈25, in a cluster. The terminal tyrosine residue is the structural determinant in the formation of helical fibril-like structures. This work supports and quantifies the two-step aggregation model where the initially formed amorphous clusters grow and, when they are large enough, rearrange into mature twisted structures. However, in addition to the nucleated fibrillation, growing aggregates undergo further internal reorganization, which leads to more compact structures of large aggregates.

Latent Models of Molecular Dynamics Data: Automatic Order Parameter Generation for Peptide Fibrillization

Variational autoencoders are artificial neural networks with the capability to reduce highly dimensional sets of data to smaller dimensional, latent representations. In this work, these models are applied to molecular dynamics simulations of the self-assembly of coarse-grained peptides to obtain a singled-valued order parameter for amyloid aggregation. This automatically learned order parameter is constructed by time-averaging the latent parametrizations of internal coordinate representations and compared to the nematic order parameter which is commonly used to study ordering of similar systems in literature. It is found that the latent space value provides more tailored insight into the aggregation mechanism's details, correctly identifying fibril formation in instances where the nematic order parameter fails to do so. A means is provided by which the latent space value can be analyzed so that the major contributing internal coordinates are identified, allowing for a direct interpretation of the latent space order parameter in terms of the behavior of the system. The latent model is found to be an effective and convenient way of representing the data from the dynamic ensemble and provides a means of reducing the dimensionality of a system whose scale exceeds molecular systems so-far considered with similar tools. This bypasses a need for researcher speculation on what elements of a system best contribute to summarizing major transitions and suggests latent models are effective and insightful when applied to large systems with a diversity of complex behaviors.

High-resolution probing of early events in amyloid-β aggregation related to Alzheimer's disease

In Alzheimer's disease (AD), soluble oligomers of amyloid-β (Aβ) are emerging as a crucial entity in driving disease progression as compared to insoluble amyloid deposits. The lacuna in establishing the structure to function relationship for Aβ oligomers prevents the development of an effective treatment for AD. While the transient and heterogeneous properties of Aβ oligomers impose many challenges for structural investigation, an effective use of a combination of NMR techniques has successfully identified and characterized them at atomic-resolution. Here, we review the successful utilization of solution and solid-state NMR techniques to probe the aggregation and structures of small and large oligomers of Aβ. Biophysical studies utilizing the commonly used solution and 19F based NMR experiments to identify the formation of small size early intermediates and to obtain their structures, and dock-lock mechanism of fiber growth at atomic-resolution are discussed. In addition, the use of proton-detected magic angle spinning (MAS) solid-state NMR experiments to obtain high-resolution insights into the aggregation pathways and structures of large oligomers and other aggregates is also presented. We expect these NMR based studies to be valuable for real-time monitoring of the depletion of monomers and the formation of toxic oligomers and high-order aggregates under a variety of conditions, and to solve the high-resolution structures of small and large size oligomers for most amyloid proteins, and therefore to develop inhibitors and drugs.

Thermo- and pH-responsive fibrillization of squid suckerin A1H1 peptide

Stimuli-responsive smart materials have attracted considerable attention with numerous applications in nanotechnology, sensing, and biomedicine. Suckerin family proteins found in squid ring teeth represent such a class of peptide-based smart materials with their self-assemblies featuring excellent thermo-plasticity and pH-dependence. Similar to block copolymers, suckerin proteins are comprised of two repeating sequence motifs, where M1 motifs are abundant in alanine and histidine residues and M2 are rich in glycine. Experimental studies of suckerin assemblies suggested that M1 regions mainly formed nano-confined β-sheets within an amorphous matrix made of M2 modules stabilizing these β-rich nano-assemblies. The histidine-containing M1 modules are believed to govern the pH- and temperature-sensitive properties of suckerin assemblies. To better understand the stimuli-responsive properties of suckerin assemblies at the molecular level, we systematically studied the self-assembly dynamics of A1H1 peptides - a representative M1 sequence - at different temperatures and pH conditions with atomistic discrete molecular dynamic simulations. Our simulations with twenty A1H1 peptides demonstrated that below the transition temperature Tagg, they could readily self-assemble from isolated monomers into well-defined β-sheet nanostructures by both primary and secondary nucleation of β-sheets and subsequent aggregation growth via elongation and coagulation. Interestingly, the dissociation of pre-formed A1H1 β-sheet nanostructures featured a melting temperature Tm higher than Tagg, exhibiting the thermal hysteresis that is characteristic of first-order phase transitions with high energy barriers. In acidic environments where all histidine residues were protonated, the stability of the A1H1 β-sheet nano-assemblies was reduced and the β-rich assemblies easily dissociated into unstructured monomers at significantly lower temperatures than in the neutral solution. The computationally derived molecular mechanisms for pH- and temperature-dependent A1H1 self-assembly will help to understand the supramolecular assembly structures and functions of the large suckerin family and aid in the future design of peptide-based stimuli-responsive smart materials.

Spontaneous self-assembly of amyloid β (1-40) into dimers

The self-assembly and fibrillation of amyloid β (Aβ) proteins is the neuropathological hallmark of Alzheimer's disease. However, the molecular mechanism of how disordered monomers assemble into aggregates remains largely unknown. In this work, we characterize the assembly of Aβ (1-40) monomers into dimers using long-time molecular dynamics simulations. Upon interaction, the monomers undergo conformational transitions, accompanied by change of the structure, leading to the formation of a stable dimer. The dimers are stabilized by interactions in the N-terminal region (residues 5-12), in the central hydrophobic region (residues 16-23), and in the C-terminal region (residues 30-40); with inter-peptide interactions focused around the N- and C-termini. The dimers do not contain long β-strands that are usually found in fibrils.

Aggregation kinetics of short peptides: All-atom and coarse-grained molecular dynamics study

Peptides can aggregate into ordered structures with different morphologies. The aggregation mechanism and evolving structures are the subject of intense research. In this paper we have used molecular dynamics to examine the sequence-dependence of aggregation kinetics for three short peptides: octaalanine (Ala8), octaasparagine (Asn8), and the heptapeptide GNNQQNY (abbreviated as GNN). First, we compared the aggregation of 20 randomly distributed peptides using the coarse-grained MARTINI force field and the atomistic OPLS-AA force field. We found that the MARTINI and OPLS-AA aggregation kinetics are similar for Ala8, Asn8, and GNN. Second, we used the MARTINI force field to study the early stages of aggregation kinetics for a larger system with 72 peptides. In the initial stage of aggregation small clusters grow by monomer addition. In the second stage, when the free monomers are depleted, the dominant cluster growth path is cluster-cluster coalescence. We quantified the aggregation kinetics in terms of rate equations. Our study shows that the initial aggregation kinetics are similar for Ala8, Asn8, and GNN but the molecular details can be different, especially for MARTINI Ala8. We hypothesize that peptide aggregation proceed in two steps. In the first step amorphous aggregates are formed, and then, in the second step, they reorganize into ordered structures. We conclude that sequence-specific differences show up in the second step of aggregation.

Top-down Multiscale Approach To Simulate Peptide Self-Assembly from Monomers

Modeling peptide assembly from monomers on large time and length scales is often intractable at the atomistic resolution. To address this challenge, we present a new approach which integrates coarse-grained (CG), mixed-resolution, and all-atom (AA) modeling in a single simulation. We simulate the initial encounter stage with the CG model, while the further assembly and reorganization stages are simulated with the mixed-resolution and AA models. We have implemented this top-down approach with new tools to automate model transformations and to monitor oligomer formations. Further, a theory was developed to estimate the optimal simulation length for each stage using a model peptide, melittin. The assembly level, the oligomer distribution, and the secondary structures of melittin simulated by the optimal protocol show good agreement with prior experiments and AA simulations. Finally, our approach and theory have been successfully validated with three amyloid peptides (β-amyloid 16-22, GNNQQNY fragment from the yeast prion protein SUP35, and α-synuclein fibril 35-55), which highlight the synergy from modeling at multiple resolutions. This work not only serves as proof of concept for multiresolution simulation studies but also presents practical guidelines for further self-assembly simulations at more physically and chemically relevant scales.

Hyperosmotic stress induces cell-dependent aggregation of α-synuclein

The aggregation of alpha-synuclein (α-syn) is a pathological feature of a number of neurodegenerative conditions, including Parkinson’s disease. Genetic mutations, abnormal protein synthesis, environmental stress, and aging have all been implicated as causative factors in this process. The importance of water in the polymerisation of monomers, however, has largely been overlooked. In the present study, we highlight the role of hyperosmotic stress in inducing human α-syn to aggregate in cells in vitro, through rapid treatment of the cells with three different osmolytes: sugar, salt and alcohol. This effect is cell-dependent and not due to direct protein-osmolyte interaction, and is specific for α-syn when compared to other neurodegeneration-related proteins, such as Tau or Huntingtin. This new property of α-syn not only highlights a unique aspect of its behaviour which may have some relevance for disease states, but may also be useful as a screening test for compounds to inhibit the aggregation of α-syn in vitro.

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.

Paradoxical Effect of Trehalose on the Aggregation of α-Synuclein: Expedites Onset of Aggregation yet Reduces Fibril Load

Aggregation of α-synuclein is closely connected to the pathology of Parkinson's disease. The phenomenon involves multiple steps, commenced by partial misfolding and eventually leading to mature amyloid fibril formation. Trehalose, a widely accepted osmolyte, has been shown previously to inhibit aggregation of various globular proteins owing to its ability to prevent the initial unfolding of protein. In this study, we have examined if it behaves in a similar fashion with intrinsically disordered protein α-synuclein and possesses the potential to act as therapeutic agent against Parkinson's disease. It was observed experimentally that samples coincubated with trehalose fibrillate faster compared to the case in its absence. Molecular dynamics simulations suggested that this initial acceleration is manifestation of trehalose's tendency to perturb the conformational transitions between different conformers of monomeric protein. It stabilizes the aggregation prone "extended" conformer of α-synuclein, by binding to its exposed acidic residues of the C terminus. It also favors the β-rich oligomers once formed. Interestingly, the total fibrils formed are still promisingly less since it accelerates the competing pathway toward formation of amorphous aggregates.

The attachment of α-synuclein to a fiber: A coarse-grain approach

We present simulations of the amyloidogenic core of α-synuclein, the protein causing Parkinson's disease, as a short chain of coarse-grain patchy particles. Each particle represents a sequence of about a dozen amino acids. The fluctuating secondary structure of this intrinsically disordered protein is modelled by dynamic variations of the shape and interaction characteristics of the patchy particles, ranging from spherical with weak isotropic attractions for the disordered state to spherocylindrical with strong directional interactions for a β-sheet. Flexible linkers between the particles enable sampling of the tertiary structure. This novel model is applied here to study the growth of an amyloid fibril, by calculating the free energy profile of a protein attaching to the end of a fibril. The simulation results suggest that the attaching protein readily becomes trapped in a mis-folded state, thereby inhibiting further growth of the fibril until the protein has readjusted to conform to the fibril structure, in line with experimental findings and previous simulations on small fragments of other proteins.

Intrinsic Conformational Preferences and Interactions in α-Synuclein Fibrils: Insights from Molecular Dynamics Simulations

Amyloid formation by the intrinsically disordered α-synuclein protein is the hallmark of Parkinson's disease. We present atomistic Molecular Dynamics simulations of the core of α-synuclein using enhanced sampling techniques to describe the conformational and binding free energy landscapes of fragments implicated in fibril stabilization. The theoretical framework is derived to combine the free energy profiles of the fragments into the reaction free energy of a protein binding to a fibril. Our study shows that individual fragments in solution have a propensity toward attaining non-β conformations, indicating that in a fibril β-strands are stabilized by interactions with other strands. We show that most dimers of hydrogen-bonded fragments are unstable in solution, while hydrogen bonding stabilizes the collective binding of five fragments to the end of a fibril. Hydrophobic effects make further contributions to the stability of fibrils. This study is the first of its kind where structural and binding preferences of the five major fragments of the hydrophobic core of α-synuclein have been investigated. This approach improves sampling of intrinsically disordered proteins, provides information on the binding mechanism between the core sequences of α-synuclein, and enables the parametrization of coarse grained models.

Energy Landscapes for the Aggregation of Aβ17-42

The aggregation of the Aβ peptide (Aβ_1-42) to form fibrils is a key feature of Alzheimer's disease. The mechanism is thought to be a nucleation stage followed by an elongation process. The elongation stage involves the consecutive addition of monomers to one end of the growing fibril. The aggregation process proceeds in a stop-and-go fashion and may involve off-pathway aggregates, complicating experimental and computational studies. Here we present exploration of a well-defined region in the free and potential energy landscapes for the Aβ_17-42 pentamer. We find that the ideal aggregation process agrees with the previously reported dock-lock mechanism. We also analyze a large number of additional stable structures located on the multifunnel energy landscape, which constitute kinetic traps. The key contributors to the formation of such traps are misaligned strong interactions, for example the stacking of F19 and F20, as well as entropic contributions. Our results suggest that folding templates for aggregation are a necessity and that aggregation studies could employ such species to obtain a more detailed description of the process.

An adaptive bias - hybrid MD/kMC algorithm for protein folding and aggregation

In this paper, we present a novel hybrid Molecular Dynamics/kinetic Monte Carlo (MD/kMC) algorithm and apply it to protein folding and aggregation in explicit solvent. The new algorithm uses a dynamical definition of biases throughout the MD component of the simulation, normalized in relation to the unbiased forces. The algorithm guarantees sampling of the underlying ensemble in dependency of one average linear coupling factor 〈α〉_τ. We test the validity of the kinetics in simulations of dialanine and compare dihedral transition kinetics with long-time MD-simulations. We find that for low 〈α〉_τ values, kinetics are in good quantitative agreement. In folding simulations of TrpCage and TrpZip4 in explicit solvent, we also find good quantitative agreement with experimental results and prior MD/kMC simulations. Finally, we apply our algorithm to study growth of the Alzheimer Amyloid Aβ 16-22 fibril by monomer addition. We observe two possible binding modes, one at the extremity of the fibril (elongation) and one on the surface of the fibril (lateral growth), on timescales ranging from ns to 8 μs.

Inhibition of GNNQQNY prion peptide aggregation by trehalose: a mechanistic view

Deposition of amyloid fibrils is the seminal event in the pathogenesis of numerous neurodegenerative diseases. The formation of this amyloid assembly is the manifestation of a cascade of structural transitions including toxic oligomer formation in the early stages of aggregation. Thus a viable therapeutic strategy involves the use of small molecular ligands to interfere with this assembly. In this perspective, we have explored the kinetics of aggregate formation of the fibril forming GNNQQNY peptide fragment from the yeast prion protein SUP35 using multiple all atom MD simulations with explicit solvent and provided mechanistic insights into the way trehalose, an experimentally known aggregation inhibitor, modulates the aggregation pathway. The results suggest that the assimilation process is impeded by different barriers at smaller and larger oligomeric sizes: the initial one being easily surpassed at higher temperatures and peptide concentrations. The kinetic profile demonstrates that trehalose delays the aggregation process by increasing both these activation barriers, specifically the latter one. It increases the sampling of small-sized aggregates that lack the beta sheet conformation. Analysis reveals that the barrier in the growth of larger stable oligomers causes the formation of multiple stable small oligomers which then fuse together bimolecularly. The PCA of 26 properties was carried out to deconvolute the events within the temporary lag phases, which suggested dynamism in lags involving an increase in interchain contacts and burial of SASA. The predominant growth route is monomer addition, which changes to condensation on account of a large number of depolymerisation events in the presence of trehalose. The favourable interaction of trehalose specifically with the sidechain of the peptide promotes crowding of trehalose molecules in its vicinity - the combination of both these factors imparts the observed behaviour. Furthermore, increasing trehalose concentration leads to faster expulsion of water molecules than interpeptide interactions. These expelled water molecules have larger translational movement, suggesting an entropy factor to favor the assembly process. Different conformations observed under this condition suggest the role of water molecules in guiding the morphology of the aggregates as well. A similar scenario exists on increasing peptide concentration.

Asparagine and glutamine ladders promote cross-species prion conversion

Prion transmission between species is governed in part by primary sequence similarity between the infectious prion aggregate, PrPSc, and the cellular prion protein of the host, PrPC A puzzling feature of prion formation is that certain PrPC sequences, such as that of bank vole, can be converted by a remarkably broad array of different mammalian prions, whereas others, such as rabbit, show robust resistance to cross-species prion conversion. To examine the structural determinants that confer susceptibility or resistance to prion conversion, we systematically tested over 40 PrPC variants of susceptible and resistant PrPC sequences in a prion conversion assay. Five key residue positions markedly impacted prion conversion, four of which were in steric zipper segments where side chains from amino acids tightly interdigitate in a dry interface. Strikingly, all five residue substitutions modulating prion conversion involved the gain or loss of an asparagine or glutamine residue. For two of the four positions, Asn and Gln residues were not interchangeable, revealing a strict requirement for either an Asn or Gln residue. Bank voles have a high number of Asn and Gln residues and a high Asn:Gln ratio. These findings suggest that a high number of Asn and Gln residues at specific positions may stabilize β-sheets and lower the energy barrier for cross-species prion transmission, potentially because of hydrogen bond networks from side chain amides forming extended Asn/Gln ladders. These data also suggest that multiple PrPC segments containing Asn/Gln residues may act in concert along a replicative interface to promote prion conversion.

Probing oligomerization of amyloid beta peptide in silico

Aggregation of amyloid β (Aβ) peptide is implicated in fatal Alzheimer's disease, for which no cure is available. Understanding the mechanisms responsible for this aggregation is required in order for therapies to be developed. In an effort to better understand the molecular mechanisms involved in spontaneous aggregation of Aβ peptide, extensive molecular dynamics simulations are reported, and the results are analyzed through a combination of structural biology tools and a novel essential collective dynamics method. Several model systems composed of ten or twelve Aβ_17-42 chains in water are investigated, and the influence of metal ions is probed. The results suggest that Aβ monomers tend to aggregate into stable globular-like oligomers with 13-23% of β-sheet content. Two stages of oligomer formation have been identified: quick collapse within the first 40 ns of the simulation, characterized by a decrease in inter-chain separation and build-up of β-sheets, and the subsequent slow relaxation of the oligomer structure. The resulting oligomers comprise a stable, coherently moving sub-aggregate of 6-9 strongly inter-correlated chains. Cu²⁺ and Fe²⁺ ions have been found to develop coordination bonds with carboxylate groups of E22, D23 and A42, which remain stable during 200 ns simulations. The presence of Fe²⁺, and particularly Cu²⁺ ions, in negatively charged cavities has been found to cause significant changes in the structure and dynamics of the oligomers. The results indicate, in particular, that formation of non-fibrillar oligomers might be involved in early template-free aggregation of Aβ_17-42 monomers, with charged species such as Cu²⁺ or Fe²⁺ ions playing an important role.

Role of Hydrophobic/Aromatic Residues on the Stability of Double-Wall β-Sheet Structures Formed by a Triblock Peptide

Bioinspired self-assembling peptides serve as powerful building blocks in the manufacturing of nanomaterials with tailored features. Because of their ease of synthesis, biocompatibility, and tunable activity, this emerging branch of biomolecules has become very popular. The triblock peptide architecture designed by the Hartgerink group is a versatile system that allows control over its assembly and has been shown to demonstrate tunable bioactivity. Three main forces, Coulomb repulsion, hydrogen bonding and hydrophobicity act together to guide the triblock peptides' assembly into one-dimensional objects and hydrogels. It was shown previously that both the nanofiber morphology (e.g., intersheet spacing, formation of antiparallel/parallel β-sheets) and hydrogel rheology strictly depend on the choice of the core residue where the triblock peptide fibers with aromatic cores in general form shorter fibers and yield poor hydrogels with respect to the ones with aliphatic cores. However, an elaborate understanding of the molecular reasons behind these changes remained unclear. In this study, by using carefully designed computer based free energy calculations, we analyzed the influence of the core residue on the formation of double-wall fibers and single-wall β-sheets. Our results demonstrate that the aromatic substitution impairs the fiber cores and this impairment is mainly associated with a reduced hydrophobic character of the aromatic side chains. Such weakening is most obvious in tryptophan containing peptides where the fiber core absorbs a significant amount of water. We also show that the ability of tyrosine to form side chain hydrogen bonds plays an indispensable role in the fiber stability. As opposed to the impairment of the fiber cores, single-wall β-sheets with aromatic faces become more stable compared to the ones with aliphatic faces suggesting that the choice of the core residue can also affect the underlying assembly mechanism. We also provide an in-depth comparison of competing structures (zero-dimensional aggregates, short and long fibers) in the triblock peptides' assembly and show that by adjusting the length of the terminal blocks, the fiber growth can be turned on or off while keeping the nanofiber morphology intact.

Elucidating Important Sites and the Mechanism for Amyloid Fibril Formation by Coarse-Grained Molecular Dynamics

Fibrils formed by the β-amyloid (Aβ) peptide play a central role in the development of Alzheimer's disease. In this study, the principles governing their growth and stability are investigated by analyzing canonical and replica-exchange molecular dynamics trajectories of Aβ_(9-40) fibrils. In particular, an unstructured monomer was allowed to interact freely with an Aβ fibril template. Trajectories were generated with the coarse-grained united-residue force field, and one- and two-dimensional free-energy landscapes (FELs) along the backbone virtual-bond angle θ and backbone virtual-bond-dihedral angle γ of each residue and principal components, respectively, were analyzed. Also, thermal unbinding (unfolding) of an Aβ peptide from the fibril template was investigated. These analyses enable us to illustrate the entire process of Aβ fibril elongation and to elucidate the key residues involved in it. Several different pathways were identified during the search for the fibril conformation by the monomer, which finally follows a dock-lock mechanism with two distinct locking stages. However, it was found that the correct binding, with native hydrogen bonds, of the free monomer to the fibril template at both stages is crucial for fibril elongation. In other words, if the monomer is incorrectly bound (with nonnative hydrogen bonds) to the fibril template during the first "docking" stage, it can remain attached to it for a long time before it dissociates and either attempts a different binding or allows another monomer to bind. This finding is consistent with an experimentally observed "stop-and-go" mechanism of fibril growth.

Inhibitory effect of hydrophobic fullerenes on the β-sheet-rich oligomers of a hydrophilic GNNQQNY peptide revealed by atomistic simulations

Fullerenes suppress fibril-like β-sheet oligomers by interacting strongly with the nonpolar aliphatic groups of polar residues of GNNQQNY peptide, thus inhibit peptide aggregation.