Amyloid fibril formation by A beta 16-22, a seven-residue fragment of the Alzheimer's beta-amyloid peptide, and structural characterization by solid state NMR

Solvent-Assisted Tyrosine-Based Dipeptide Forms Low-Molecular Weight Gel: Preparation and Its Potential Use in Dye Removal and Oil Spillage Separation from Water

Low-molecular weight gelators (supramolecular, or simply molecular gels) are highly important molecular frameworks because of their potential application in drug delivery, catalysis, pollutant removal, sensing materials, and so forth. Herein, a small dipeptide composed of N-(tert-butoxycarbonyl)pentafluoro-l-phenylalanine and O-benzyl-l-tyrosine methyl ester was synthesized, and its gelation ability was investigated in different solvent systems. It was found that the dipeptide was unable to form gel with a single solvent, but a mixture of solvent systems was found to be suitable for the gelation of this dipeptide. Interestingly, water was found to be essential for gelation with the polar protic solvent, and long-chain hydrocarbon units such as, petroleum ether, kerosene, and diesel, were important for gelation with aromatic solvents. The structural insights of these gels were characterized by field-emission scanning electronic microscopy, atomic force microscopy, Fourier transform infrared analysis, and X-ray diffraction studies, and their mechanical strengths were characterized by rheological experiments. Both of the gels obtained from these two solvent systems were thermoreversible in nature, and these translucent gels had potential application for the treatment of waste water. The gel obtained from dipeptides with methanol-water was used to remove toxic dyes (crystal violet, Eriochrome Black T, and rhodamine B) from water. Furthermore, the gel obtained from dipeptide with assistance from toluene-petroleum ether was used as a phase-selective gelator for oil-spill recovery.

Guiding the Morphology of Amyloid Assemblies by Electrostatic Capping: from Polymorphic Twisted Fibrils to Uniform Nanorods

Peptides that self-assemble into cross-β-sheet amyloid structures constitute promising building blocks to construct highly ordered proteinaceous materials and nanoparticles. Nevertheless, the intrinsic polymorphism of amyloids and the difficulty of controlling self-assembly currently limit their usage. In this study, the effect of electrostatic interactions on the supramolecular organization of peptide assemblies is investigated to gain insights into the structural basis of the morphological diversities of amyloids. Different charged capping units are introduced at the N-terminus of a potent β-sheet-forming sequence derived from the 20-29 segment of islet amyloid polypeptide, known to self-assemble into polymorphic fibrils. By tuning the charge and the electrostatic strength, different mesoscopic morphologies are obtained, including nanorods, rope-like fibrils, and twisted ribbons. Particularly, the addition of positive capping units leads to the formation of uniform rod-like assemblies, with lengths that can be modulated by the charge number. It is proposed that electrostatic repulsions between N-terminal positive charges hinder β-sheet tape twisting, leading to a unique control over the size of these cytocompatible nanorods by protofilament growth frustration. This study reveals the high susceptibility of amyloid formation to subtle chemical modifications and opens to promising strategies to control the final architecture of proteinaceous assemblies from the peptide sequence.

Limpid hydrogels from β-turn motif-connected tandem repeats of Aβ16-22

Self-assembling peptides constitute an important class of functional biomaterials. A number of short amyloidogenic stretches have been identified from amyloid proteins. Such peptides, as such or through subtle modifications, can turn out to be promising candidates for functional biomaterials. End-capped Aβ16-22, the well-studied amyloidogenic stretch from β-amyloid, is reported to be non-hydrogelating up to 20 mM concentration. Here we investigated the hydrogelation propensity of Aβ16-22 repeats connected through β-turn-supporting motifs. The peptide repeats connected through Asn-Gly, Aib-DPro, and DPro-Gly formed transparent hydrogels at concentrations ≥2 mM. The repeats of the aromatic analog Aβ16-22(F20Y) also resulted in similar hydrogels. Like other peptide-based gels reported earlier, these gels could trap the anticancer drug doxorubicin and displayed steady release in water. In addition, the gels supported the growth of mammalian cell lines, HEK-293 and RIN-5F. These data show that turn-inducing motifs can have marked effects on the hydrogelating propensity of self-assembling peptides.

Designed Two- and Three-Dimensional Protein Nanocage Networks Driven by Hydrophobic Interactions Contributed by Amyloidogenic Motifs

Precise manipulation of protein self-assembly by noncovalent interactions into programmed networks to mimic naturally occurring nanoarchitectures in living organisms is a challenge due to its structural heterogeneity, flexibility, and complexity. Herein, by taking advantage of both the hydrophobic forces contributed by the "GLMVG" motif, a kind of amyloidogenic motif (AM), and the high symmetry of protein nanocages, we have built an effective protein self-assembly strategy for the construction of two-dimensional (2D) or three-dimensional (3D) protein nanocage arrays. According to this strategy, "GLMVG" AMs from β-amyloid 42 were grafted onto the outer surface of a 24-mer ferritin nanocage close to its C₄ symmetry channels, initially resulting in the production of subgrade 2D nanocage arrays and ultimately generating 3D highly ordered arrays with a simple cubic packing pattern as the reaction time increases. More importantly, the reversibility and the formation rate of these protein arrays can be modulated by pH. This work provides a de novo design strategy for accurate control over 2D or 3D protein self-assemblies.

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.

Effects of Metal Ions on Aβ42 Peptide Conformations from Molecular Simulation Studies

In this study, we investigate the conformational characteristics of full-length Aβ₄₂ peptide monomers in the presence of Na⁺ and Zn²⁺ metal ions using atomistic molecular dynamics (MD) simulations with an aim to explore the possible driving forces behind enhanced aggregation rates of the peptides in the presence of salts. The calculations reveal that the presence of metal ions shifts the conformational equilibrium more toward the compact ordered Aβ structures. Such compact ordered structures stabilized by distant nonlocal contacts between two crucial hydrophobic segments, hp1 and hp2, primarily through two important hydrophobic aromatic residues, Phe-19 and Phe-20, are expected to trigger the aggregation process at a faster rate by populating and stabilizing the aggregation prone structures. Formation of a significant number of such distant contacts in the presence of Na⁺ ions has also been found to result in breaking of the N-terminal helix. On the contrary, binding of Zn²⁺ ion to Aβ peptide is highly specific, which stabilizes the N-terminal helix instead of breaking it. This explains why the aggregation rate of Aβ peptides is higher in the presence of divalent Zn²⁺ ions than monovalent Na⁺ ions. Relatively higher overall stability of the most populated Aβ peptide monomers in the presence of Zn²⁺ ions has been found to be associated with specific Zn²⁺-Aβ binding and significant free energy gain.

Rippled β-Sheet Formation by an Amyloid-β Fragment Indicates Expanded Scope of Sequence Space for Enantiomeric β-Sheet Peptide Coassembly

In 1953, Pauling and Corey predicted that enantiomeric β-sheet peptides would coassemble into so-called "rippled" β-sheets, in which the β-sheets would consist of alternating l- and d-peptides. To date, this phenomenon has been investigated primarily with amphipathic peptide sequences composed of alternating hydrophilic and hydrophobic amino acid residues. Here, we show that enantiomers of a fragment of the amyloid-β (Aβ) peptide that does not follow this sequence pattern, amyloid-β (16-22), readily coassembles into rippled β-sheets. Equimolar mixtures of enantiomeric amyloid-β (16-22) peptides assemble into supramolecular structures that exhibit distinct morphologies from those observed by self-assembly of the single enantiomer pleated β-sheet fibrils. Formation of rippled β-sheets composed of alternating l- and d-amyloid-β (16-22) is confirmed by isotope-edited infrared spectroscopy and solid-state NMR spectroscopy. Sedimentation analysis reveals that rippled β-sheet formation by l- and d-amyloid-β (16-22) is energetically favorable relative to self-assembly into corresponding pleated β-sheets. This work illustrates that coassembly of enantiomeric β-sheet peptides into rippled β-sheets is not limited to peptides with alternating hydrophobic/hydrophilic sequence patterns, but that a broader range of sequence space is available for the design and preparation of rippled β-sheet materials.

Multiscale Aggregation of the Amyloid Aβ16-22 Peptide: From Disordered Coagulation and Lateral Branching to Amorphous Prefibrils

In this work we investigate the multiscale dynamics of the aggregation process of an amyloid peptide, Aβ_16-22. By performing massive coarse-grained simulations at the quasi-atomistic resolution and including hydrodynamic effects, we followed the formation and growth of a large elongated aggregate and its slow structuring. The elongation proceeds via a two-step nucleation mechanism with disordered aggregates formed initially and subsequently fusing to elongate the amorphous prefibril. A variety of coagulation events coexist, including lateral growth. The latter mechanism, sustained by long-range hydrodynamics correlations, actually can create a large branched structure spanning a few tens of nanometers. Our findings confirm the experimental hypothesis of a critical contribution of lateral growth to the amyloid aggregation kinetics and the capability of our model to sample critical structures like prefibril hosting annular pores.

Effects of All-Atom Molecular Mechanics Force Fields on Amyloid Peptide Assembly: The Case of Aβ16-22 Dimer

We investigated the effects of 17 widely used atomistic molecular mechanics force fields (MMFFs) on the structures and kinetics of amyloid peptide assembly. To this end, we performed large-scale all-atom molecular dynamics simulations in explicit water on the dimer of the seven-residue fragment of the Alzheimer's amyloid-β peptide, Aβ16-22, for a total time of 0.34 ms. We compared the effects of these MMFFs by analyzing various global reaction coordinates, secondary structure contents, the fibril population, the in-register and out-of-register architectures, and the fibril formation time at 310 K. While the AMBER94, AMBER99, and AMBER12SB force fields do not predict any β-sheets, the seven force fields, AMBER96, GROMOS45a3, GROMOS53a5, GROMOS53a6, GROMOS43a1, GROMOS43a2, and GROMOS54a7, form β-sheets rapidly. In contrast, the following five force fields, AMBER99-ILDN, AMBER14SB, CHARMM22*, CHARMM36, and CHARMM36m, are the best candidates for studying amyloid peptide assembly, as they provide good balances in terms of structures and kinetics. We also investigated the assembly mechanisms of dimeric Aβ16-22 and found that the fibril formation rate is predominantly controlled by the total β-strand content.

Thermodynamic phase diagram of amyloid-β (16-22) peptide

The aggregation of monomeric amyloid β protein (Aβ) peptide into oligomers and amyloid fibrils in the mammalian brain is associated with Alzheimer's disease. Insight into the thermodynamic stability of the Aβ peptide in different polymeric states is fundamental to defining and predicting the aggregation process. Experimental determination of Aβ thermodynamic behavior is challenging due to the transient nature of Aβ oligomers and the low peptide solubility. Furthermore, quantitative calculation of a thermodynamic phase diagram for a specific peptide requires extremely long computational times. Here, using a coarse-grained protein model, molecular dynamics (MD) simulations are performed to determine an equilibrium concentration and temperature phase diagram for the amyloidogenic peptide fragment Aβ_16-22 Our results reveal that the only thermodynamically stable phases are the solution phase and the macroscopic fibrillar phase, and that there also exists a hierarchy of metastable phases. The boundary line between the solution phase and fibril phase is found by calculating the temperature-dependent solubility of a macroscopic Aβ_16-22 fibril consisting of an infinite number of β-sheet layers. This in silico determination of an equilibrium (solubility) phase diagram for a real amyloid-forming peptide, Aβ_16-22, over the temperature range of 277-330 K agrees well with fibrillation experiments and transmission electron microscopy (TEM) measurements of the fibril morphologies formed. This in silico approach of predicting peptide solubility is also potentially useful for optimizing biopharmaceutical production and manufacturing nanofiber scaffolds for tissue engineering.

Structural Analysis of a Trimer of β2-Microgloblin Fragment by Molecular Dynamics Simulations

A peptide β₂-m_21-31, which is a fragment from residue 21 to residue 31 of β₂-microgloblin, is experimentally known to self-assemble and form amyloid fibrils. In order to understand the mechanism of amyloid fibril formations, we applied the replica-exchange molecular dynamics method to the system consisting of three fragments of β₂-m_21-31. From the analyses on the temperature dependence, we found that there is a clear phase transition temperature in which the peptides aggregate with each other. Moreover, we found by the free energy analyses that there are two major stable states: One of them is like amyloid fibrils and the other is amorphous aggregates.

Mechanistic insight into E22Q-mutation-induced antiparallel-to-parallel β-sheet transition of Aβ16−22fibrils: an all-atom simulation study

E22Q mutation of Aβ₁₆₋₂₂fibrils facilitates parallel β-sheet formation by enhancing Q22–Q22 hydrogen-bonding interaction and A21–A21, F20–F20, F19–F19 and V18–V18 hydrophobic interaction.

Pathways of amyloid-beta absorption and aggregation in a membranous environment

Aggregation of misfolded oligomeric amyloid-beta (Aβ) peptides on lipid membranes has been identified as a primary event in Alzheimer's pathogenesis. However, the structural and dynamical features of this membrane assisted Aβ aggregation have not been well characterized. The microscopic characterization of dynamic molecular-level interactions in peptide aggregation pathways has been challenging both computationally and experimentally. In this work, we explore differential patterns of membrane-induced Aβ 16-22 (K-L-V-F-F-A-E) aggregation from the microscopic perspective of molecular interactions. Physics-based coarse-grained molecular dynamics (CG-MD) simulations were employed to investigate the effect of lipid headgroup charge - zwitterionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine: POPC) and anionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine: POPS) - on Aβ 16-22 peptide aggregation. Our analyses present an extensive overview of multiple pathways for peptide absorption and biomechanical forces governing peptide folding and aggregation. In agreement with experimental observations, anionic POPS molecules promote extended configurations in Aβ peptides that contribute towards faster emergence of ordered β-sheet-rich peptide assemblies compared to POPC, suggesting faster fibrillation. In addition, lower cumulative rates of peptide aggregation in POPS due to higher peptide-lipid interactions and slower lipid diffusion result in multiple distinct ordered peptide aggregates that can serve as nucleation seeds for subsequent Aβ aggregation. This study provides an in-silico assessment of experimentally observed aggregation patterns, presents new morphological insights and highlights the importance of lipid headgroup chemistry in modulating the peptide absorption and aggregation process.

Morphological Transformation of Peptide Nanoassemblies through Conformational Transition of Core-forming Peptides

Morphological control of nanostructures that are composed of amphiphilic di- or tri-block molecules by external stimuli broadens their applications for molecular containers, nanoreactors, and controlled release materials. In this study, triblock amphiphiles comprising oligo(ethylene glycol), oligo(l-lysine), and tetra(l-phenylalanine) were prepared for the construction of nanostructures that can transform accompanying α-to-β transition of core-forming peptides. Circular dichroic (CD) measurements showed that the triblock amphiphiles adopted different secondary structures depending on the solvent environment: they adopt β-sheet structures in aqueous solution, while α-helix structures in 25% 2,2,2-trifluoroethanol (TFE) solution under basic pH conditions. Transmission electron microscopic (TEM) observation revealed that the triblock amphiphiles formed vesicle structures in 25% TFE aq. Solvent exchange from 25% TFE to water induced morphological transformation from vesicles to arc-shaped nanostructures accompanying α-β conformational transition. The transformable nanostructures may be useful as novel smart nanomaterials for molecular containers and micro reactors.

RETRACTED: Peptide-induced formation of protein aggregates and amyloid fibrils in human and guinea pig αA-crystallins under physiological conditions of temperature and pH

This article has been retracted: please see Elsevier Policy on Article Withdrawal (https://www.elsevier.com/about/our-business/policies/article-withdrawal). This article has been retracted at the request of the authors. The senior author contacted the journal in a forthright manner, in an effort to preserve the scientific integrity of the literature, after discovering a significant error in the results reported in the article. The authors were recently made aware of a paper by Kim et al. (Nature Commun. 2019) which shows a spirosome structure (the enzyme aldehyde-alcohol dehydrogenase) present in E. coli (Fig. 5a) that is very similar to the structure the authors thought formed when synthetic alpha A crystallin (66-80) peptide was incubated for 24 h with recombinant guinea pig alpha A insert crystallin (see Kumarasamy et al., Figs. 7C and F, and Fig. 9). Subsequent to publication of their report, the authors later found a number of images that showed what appeared to be the same structure present in samples of their presumably purified recombinant guinea pig alpha A insert crystallin which had been incubated without peptide for 24 h. Hence, the authors now conclude that the structures shown in Figs. 7C and F, and Fig. 9 of their article published in this journal are actually due to E. coli contaminant aldehyde-alcohol dehydrogenase. The authors deeply regret this error and any inconvenience it may have caused.

On the Applicability of Force Fields To Study the Aggregation of Amyloidogenic Peptides Using Molecular Dynamics Simulations

Molecular dynamics simulations play an essential role in understanding biomolecular processes such as protein aggregation at temporal and spatial resolutions which are not attainable by experimental methods. For a correct modeling of protein aggregation, force fields must accurately represent molecular interactions. Here, we study the effect of five different force fields on the oligomer formation of Alzheimer's Aβ_16-22 peptide and two of its mutants: Aβ_16-22(F19V,F20V), which does not form fibrils, and Aβ_16-22(F19L) which forms fibrils faster than the wild type. We observe that while oligomer formation kinetics depends strongly on the force field, structural properties, such as the most relevant protein-protein contacts, are similar between them. The oligomer formation kinetics obtained with different force fields differ more from each other than the kinetics between aggregating and nonaggregating peptides simulated with a single force field. We discuss the difficulties in comparing atomistic simulations of amyloid oligomer formation with experimental observables.

A Novel Peroxidase Mimics and Ameliorates Alzheimer's Disease-Related Pathology and Cognitive Decline in Mice

Alzheimer's disease (AD) is the most common neurodegenerative disorder in the elderly, which is characterized by the accumulation of amyloid β (Aβ) plaques, oxidative stress, and neuronal loss. Therefore, clearing Aβ aggregates and reducing oxidative stress could be an effective therapeutic strategy for AD. Deuterohemin-AlaHisThrValGluLys (DhHP-6), a novel deuterohemin-containing peptide mimetic of the natural microperoxidase-11 (MP-11), shows higher antioxidant activity and stability compared to the natural microperoxidases. DhHP-6 possesses the ability of extending lifespan and alleviating paralysis in the Aβ1-42 transgenic Caenorhabditis elegans CL4176 model of AD, as shown in our previous study. Therefore, this study was aimed at exploring the neuroprotective effect of DhHP-6 in the APPswe/PSEN1dE9 transgenic mouse model of AD. DhHP-6 reduced the diameter and fiber structure of Aβ1-42 aggregation in vitro, as shown by dynamic light scattering and transmission electron microscope. DhHP-6 exerted its neuroprotective effect by inhibiting Aβ aggregation and plaque formation, and by reducing Aβ1-42 oligomers-induced neurotoxicity on HT22 (mouse hippocampal neuronal) and SH-SY5Y (human neuroblastoma) cells. In the AD mouse model, DhHP-6 significantly ameliorated cognitive decline and improved spatial learning ability in behavioral tests including the Morris water maze, Y-maze, novel object recognition, open field, and nest-building test. Moreover, DhHP-6 reduced the deposition of Aβ plaques in the cerebral cortex and hippocampus. More importantly, DhHP-6 restored the morphology of astrocytes and microglia, and significantly reduced the levels of pro-inflammatory cytokines. Our findings provide a basis for considering the non-toxic, peroxidase mimetic DhHP-6 as a new candidate drug against AD.

AggScore: Prediction of aggregation-prone regions in proteins based on the distribution of surface patches

Protein aggregation is a phenomenon that has attracted considerable attention within the pharmaceutical industry from both a developability standpoint (to ensure stability of protein formulations) and from a research perspective for neurodegenerative diseases. Experimental identification of aggregation behavior in proteins can be expensive; and hence, the development of accurate computational approaches is crucial. The existing methods for predicting protein aggregation rely mostly on the primary sequence and are typically trained on amyloid-like proteins. However, the training bias toward beta amyloid peptides may worsen prediction accuracy of such models when applied to larger protein systems. Here, we present a novel algorithm to identify aggregation-prone regions in proteins termed "AggScore" that is based entirely on three-dimensional structure input. The method uses the distribution of hydrophobic and electrostatic patches on the surface of the protein, factoring in the intensity and relative orientation of the respective surface patches into an aggregation propensity function that has been trained on a benchmark set of 31 adnectin proteins. AggScore can accurately identify aggregation-prone regions in several well-studied proteins and also reliably predict changes in aggregation behavior upon residue mutation. The method is agnostic to an amyloid-specific aggregation context and thus may be applied to globular proteins, small peptides and antibodies.

Etersalate prevents the formations of 6Aβ16-22 oligomer: An in silico study

Oligomerization of amyloid beta (Aβ) peptides has been considered as the crucially causative agent in the development of Alzheimer's disease. Etersalate, a nonsteroidal anti-inflammatory oral drug (United State Food and Drug Administration—Unique Ingredient Identifier: 653GN04T2G) was previously suggested to bind well to proto-fibrils of Aβ peptides in silico. Here, the effect of etersalate on the oligomerization of soluble Aβ_16–22 hexamer (6Aβ_16–22) were extensively investigated using temperature replica exchange molecular dynamics (REMD) simulations over ~16.8 μs in total for 48 replicas (350 ns per replica). The results reveal that etersalate can enter the inner space or bind on the surface of 6Aβ_16–22 conformations, which destabilizes the hexamer. Etersalate was predicted to able to cross the blood brain barrier using prediction of absorption, distribution, metabolism, and excretion—toxicity (preADMET) tools. Overall, although the investigation was performed with the low concentration of trial inhibitor, the obtained results indicate that etersalate is a potential drug candidate for AD through inhibiting formation of Aβ oligomers with the average binding free energy of -11.7 kcal/mol.

Understanding amyloid fibril formation using protein fragments: structural investigations via vibrational spectroscopy and solid-state NMR

It is well established that amyloid proteins play a primary role in neurodegenerative diseases. Alzheimer's, Parkinson's, type II diabetes, and Creutzfeldt-Jakob's diseases are part of a wider family encompassing more than 50 human pathologies related to aggregation of proteins. Although this field of research is thoroughly investigated, several aspects of fibrillization remain misunderstood, which in turn slows down, or even impedes, advances in treating and curing amyloidoses. To solve this problem, several research groups have chosen to focus on short fragments of amyloid proteins, sequences that have been found to be of great importance for the amyloid formation process. Studying short peptides allows bypassing the complexity of working with full-length proteins and may provide important information relative to critical segments of amyloid proteins. To this end, efficient biophysical tools are required. In this review, we focus on two essential types of spectroscopic techniques, i.e., vibrational spectroscopy and its derivatives (conventional Raman scattering, deep-UV resonance Raman (DUVRR), Raman optical activity (ROA), surface-enhanced Raman spectroscopy (SERS), tip-enhanced Raman spectroscopy (TERS), infrared (IR) absorption spectroscopy, vibrational circular dichroism (VCD)) and solid-state nuclear magnetic resonance (ssNMR). These techniques revealed powerful to provide a better atomic and molecular comprehension of the amyloidogenic process and fibril structure. This review aims at underlining the information that these techniques can provide and at highlighting their strengths and weaknesses when studying amyloid fragments. Meaningful examples from the literature are provided for each technique, and their complementarity is stressed for the kinetic and structural characterization of amyloid fibril formation.

Characterization of Interfacial Structure in Polymer-Fullerene Bulk Heterojunctions via ^{13}C {^{2}H} Rotational Echo Double Resonance NMR

We introduce a new application of solid state NMR measurements towards characterizing the donor-acceptor interfaces within bulk heterojunction (BHJ) films. Rotational echo double resonance (REDOR) is used to measure dipolar couplings between ^{13}C nuclei on the acceptor phenyl-C_{61}-butyric acid methyl ester (PCBM) fullerene cage, which is ≈18% isotopically enriched with ^{13}C, and beta hydrogens on the donor poly(3-hexyl thiophene) (P3HT) main chain, which are >95% isotopically enriched with ^{2}H. The ^{13}C-^{2}H dipolar couplings are used for constraining possible models of molecular packing in the amorphous mixed phase of a P3HT/PCBM BHJ. The films studied are highly mixed (>80%) and have a maximum length scale of composition nonuniformity of ≈6  nm in the mixed phase, as demonstrated by ^{1}H spin diffusion NMR and supported by TEM. The REDOR results show that despite the lack of phase separation at length scales greater than ≈6  nm, neat P3HT and PCBM clusters exist on ≈3  nm size scales, and, for the average PCBM molecule, the number of nearest neighbors P3HTs is two.

Kinetics and mechanical stability of the fibril state control fibril formation time of polypeptide chains: A computational study

Fibril formation resulting from protein misfolding and aggregation is a hallmark of several neurodegenerative diseases such as Alzheimer's and Parkinson's diseases. Despite much progress in the understanding of the protein aggregation process, the factors governing fibril formation rates and fibril stability have not been fully understood. Using lattice models, we have shown that the fibril formation time is controlled by the kinetic stability of the fibril state but not by its energy. Having performed all-atom explicit solvent molecular dynamics simulations with the GROMOS43a1 force field for full-length amyloid beta peptides Aβ₄₀ and Aβ₄₂ and truncated peptides, we demonstrated that kinetic stability can be accessed via mechanical stability in such a way that the higher the mechanical stability or the kinetic stability, the faster the fibril formation. This result opens up a new way for predicting fibril formation rates based on mechanical stability that may be easily estimated by steered molecular dynamics.

Ultrafast Hydrogen-Bonding Dynamics in Amyloid Fibrils

While there are many studies on the subject of hydrogen-bonding dynamics in biological systems, few, if any, have investigated this fundamental process in amyloid fibrils. Herein, we seek to add insight into this topic by assessing the dynamics of a hydrogen bond buried in the dry interface of amyloid fibrils. To prepare a suitable model peptide system for this purpose, we introduce two mutations into the amyloid-forming Aβ_16-22 peptide. The first one is a lysine analogue at position 19, which is used to help form structurally homogeneous fibrils, and the second one is an aspartic acid derivative (D_M) at position 17, which is intended (1) to be used as a site-specific infrared probe and (2) to serve as a hydrogen-bond acceptor to lysine so that an inter-β-sheet hydrogen bond can be formed in the fibrils. Using both infrared spectroscopy and atomic force microscopy, we show that (1) this mutant peptide indeed forms well-defined fibrils, (2) when bulk solvent is removed, there is no detectable water present in the fibrils, (3) infrared results obtained with the D_M probe are consistent with a protofibril structure that is composed of two antiparallel β-sheets stacked in a parallel fashion, leading to formation of the expected hydrogen bond. Using two-dimensional infrared spectroscopy, we further show that the dynamics of this hydrogen bond occur on a time scale of ∼2.3 ps, which is attributed to the rapid rotation of the -NH₃⁺ group of lysine around its Cε-Nζ bond. Taken together, these results suggest that (1) D_M is a useful infrared marker in facilitating structure determination of amyloid fibrils and (2) even in the tightly packed core of amyloid fibrils certain amino acid side chains can undergo ultrafast motions, hence contributing to the thermodynamic stability of the system.

Exploring the Mechanism of Inhibition of Au Nanoparticles on the Aggregation of Amyloid-β(16-22) Peptides at the Atom Level by All-Atom Molecular Dynamics

The abnormal self-assembly of the amyloid-β peptide into toxic β-rich aggregates can cause Alzheimer’s disease. Recently, it has been shown that small gold nanoparticles (AuNPs) inhibit Aβ aggregation and fibrillation by slowing down the nucleation process in experimental studies. However, the effects of AuNPs on Aβ oligomeric structures are still unclear. In this study, we investigate the conformation of Aβ(16-22) tetramers/octamers in the absence and presence of AuNPs using extensive all-atom molecular-dynamics simulations in explicit solvent. Our studies demonstrate that the addition of AuNPs into Aβ(16-22) solution prevents β-sheet formation, and the inhibition depends on the concentration of Aβ(16-22) peptides. A detailed analysis of the Aβ(16-22)/Aβ(16-22)/water/AuNPs interactions reveals that AuNPs inhibit the β-sheet formation resulting from the same physical forces: hydrophobic interactions. Overall, our computational study provides evidence that AuNPs are likely to inhibit Aβ(16-22) and full-length Aβ fibrillation. Thus, this work provides theoretical insights into the development of inorganic nanoparticles as drug candidates for treatment of AD.

Guiding principles for peptide nanotechnology through directed discovery

Life's diverse molecular functions are largely based on only a small number of highly conserved building blocks - the twenty canonical amino acids. These building blocks are chemically simple, but when they are organized in three-dimensional structures of tremendous complexity, new properties emerge. This review explores recent efforts in the directed discovery of functional nanoscale systems and materials based on these same amino acids, but that are not guided by copying or editing biological systems. The review summarises insights obtained using three complementary approaches of searching the sequence space to explore sequence-structure relationships for assembly, reactivity and complexation, namely: (i) strategic editing of short peptide sequences; (ii) computational approaches to predicting and comparing assembly behaviours; (iii) dynamic peptide libraries that explore the free energy landscape. These approaches give rise to guiding principles on controlling order/disorder, complexation and reactivity by peptide sequence design.

Identification and application of self-binding zipper-like sequences in SARS-CoV spike protein

Self-binding peptides containing zipper-like sequences, such as the Leu/Ile zipper sequence within the coiled coil regions of proteins and the cross-β spine steric zippers within the amyloid-like fibrils, could bind to the protein-of-origin through homophilic sequence-specific zipper motifs. These self-binding sequences represent opportunities for the development of biochemical tools and/or therapeutics. Here, we report on the identification of a putative self-binding β-zipper-forming peptide within the severe acute respiratory syndrome-associated coronavirus spike (S) protein and its application in viral detection. Peptide array scanning of overlapping peptides covering the entire length of S protein identified 34 putative self-binding peptides of six clusters, five of which contained octapeptide core consensus sequences. The Cluster I consensus octapeptide sequence GINITNFR was predicted by the Eisenberg’s 3D profile method to have high amyloid-like fibrillation potential through steric β-zipper formation. Peptide C6 containing the Cluster I consensus sequence was shown to oligomerize and form amyloid-like fibrils. Taking advantage of this, C6 was further applied to detect the S protein expression in vitro by fluorescence staining. Meanwhile, the coiled-coil-forming Leu/Ile heptad repeat sequences within the S protein were under-represented during peptide array scanning, in agreement with that long peptide lengths were required to attain high helix-mediated interaction avidity. The data suggest that short β-zipper-like self-binding peptides within the S protein could be identified through combining the peptide scanning and predictive methods, and could be exploited as biochemical detection reagents for viral infection.

Multi-Target Anti-Alzheimer Activities of Four Prenylated Compounds from Psoralea Fructus

Alzheimer’s disease (AD) is an age-related neurodegenerative disease that is mediated by multiple signaling pathways. In recent years, the components of Psoralea Fructus (PF) have demonstrated some anti-Alzheimer effects both in vitro and in vivo. To further reveal the active compounds of PF and their mechanisms regulating key targets of AD, in this study, we identified four prenylated compounds from the 70% ethanolic aqueous extract of PF, namely bavachin, bavachinin, bavachalcone, and isobavachalcone. Multi-target bioactivity analysis showed that these compounds could differentially inhibit neuroinflammation, oxidative damage, and key AD-related protein targets, such as amyloid β-peptide 42, β-secretase, glycogen synthase kinase 3β, and acetylcholinesterase. These compounds may generate beneficial effects in AD prevention and treatment.

Solid-State NMR Structural Characterization of Self-Assembled Peptides with Selective 13C and 15N Isotopic Labels

For the structural characterization methods discussed here, information on molecular conformation and intermolecular organization within nanostructured peptide assemblies is discerned through analysis of solid-state NMR spectral features. This chapter reviews general NMR methodologies, requirements for sample preparation, and specific descriptions of key experiments. An attempt is made to explain choices of solid-state NMR experiments and interpretation of results in a way that is approachable to a nonspecialist. Measurements are designed to determine precise NMR peak positions and line widths, which are correlated with secondary structures, and probe nuclear spin-spin interactions that report on three-dimensional organization of atoms. The formulation of molecular structural models requires rationalization of data sets obtained from multiple NMR experiments on samples with carefully chosen 13C and 15N isotopic labels. The information content of solid-state NMR data has been illustrated mostly through the use of simulated data sets and references to recent structural work on amyloid fibril-forming peptides and designer self-assembling peptides.

Methods to Characterize the Nanostructure and Molecular Organization of Amphiphilic Peptide Assemblies

Methods to characterize the nanostructure and molecular organization of aggregates of peptides such as amyloid or amphiphilic peptide assemblies are reviewed. We discuss techniques to characterize conformation and secondary structure including circular and linear dichroism and FTIR and Raman spectroscopies, as well as fluorescence methods to detect aggregation. NMR spectroscopy methods, especially solid-state NMR measurements to probe beta-sheet packing motifs, are also briefly outlined. Also discussed are scattering methods including X-ray diffraction and small-angle scattering techniques including SAXS (small-angle X-ray scattering) and SANS (small-angle neutron scattering) and dynamic light scattering. Imaging methods are direct methods to uncover features of peptide nanostructures, and we provide a summary of electron microscopy and atomic force microscopy techniques. Selected examples are highlighted showing data obtained using these techniques, which provide a powerful suite of methods to probe ordering from the molecular scale to the aggregate superstructure.

Boc-Val-Val-OMe (Aβ39–40) and Boc-Ile-Ala-OMe (Aβ41–42) crystallize in a parallel β-sheet arrangement but generate a different morphology

Crystal structures and morphologies of Boc-Val-Val-OMe (Aβ_39–40) and Boc-Ile-Ala-OMe (Aβ_41–42), the two consecutive dipeptides of the C-terminus of Aβ, are reported.

Amyloid-like aggregation of designer bolaamphiphilic peptides: Effect of hydrophobic section and hydrophilic heads

Amyloid-like aggregation of natural proteins or polypeptides is an important process involved in many human diseases as well as some normal biological functions. Plenty of works have been done on this ubiquitous phenomenon, but the molecular mechanism of amyloid-like aggregation has not been fully understood yet. In this study, we showed that a series of designer bolaamphiphilic peptides could undergo amyloid-like aggregation even though they didn't possess typical β-sheet secondary structure. Through systematic amino acid substitution, we found that for the self-assembling ability, the number and species of amino acid in hydrophobic section could be variable as long as enough hydrophobic interaction is provided, while different polar amino acids as the hydrophilic heads could change the self-assembling nanostructures with their aggregating behaviors affected by pH value change. Based on these results, novel self-assembling models and aggregating mechanisms were proposed, which might provide new insight into the molecular basis of amyloid-like aggregation.

β-Sheet Structure within the Extracellular Domain of C99 Regulates Amyloidogenic Processing

Familial mutations in C99 can increase the total level of the soluble Aβ peptides produced by proteolysis, as well as the Aβ42/Aβ40 ratio, both of which are linked to the progression of Alzheimer’s disease. We show that the extracellular sequence of C99 forms β-sheet structure upon interaction with membrane bilayers. Mutations that disrupt this structure result in a significant increase in Aβ production and, in specific cases, result in an increase in the amount of Aβ42 relative to Aβ40. Fourier transform infrared and solid-state NMR spectroscopic studies reveal a central β-hairpin within the extracellular sequence comprising Y10-E11-V12 and L17-V18-F19 connected by a loop involving H13-H14-Q15. These results suggest how familial mutations in the extracellular sequence influence C99 processing and provide a structural basis for the development of small molecule modulators that would reduce Aβ production.

Assemblies of amyloid-β30-36 hexamer and its G33V/L34T mutants by replica-exchange molecular dynamics simulation

The aggregation of amyloid-β peptides is associated with the pathogenesis of Alzheimer’s disease, in which the 30–36 fragments play an important part as a fiber-forming hydrophobic region. The fibrillar structure of Aβ_30–36 has been detected by means of X-ray diffraction, but its oligomeric structural determination, biophysical characterization, and pathological mechanism remain elusive. In this study, we have investigated the structures of Aβ_30–36 hexamer as well as its G33V and L34T mutants in explicit water environment using replica-exchange molecular dynamics (REMD) simulations. Our results show that the wild-type (WT) Aβ_30–36 hexamer has a preference to form β-barrel and bilayer β-sheet conformations, while the G33V or L34T mutation disrupts the β-barrel structures: the G33V mutant is homogenized to adopt β-sheet-rich bilayers, and the structures of L34T mutant on the contrary get more diverse. The hydrophobic interaction plays a critical role in the formation and stability of oligomeric assemblies among all the three systems. In addition, the substitution of G33 by V reduces the β-sheet content in the most populated conformations of Aβ_30–36 oligomers through a steric effect. The L34T mutation disturbs the interpeptide hydrogen bonding network, and results in the increased coil content and morphological diversity. Our REMD runs provide structural details of WT and G33V/L34T mutant Aβ_30–36 oligomers, and molecular insight into the aggregation mechanism, which will be helpful for designing novel inhibitors or amyloid-based materials.

Multistep Conformation Selection in Amyloid Assembly

Defining pathways for amyloid assembly could impact therapeutic strategies for as many as 50 disease states. Here we show that amyloid assembly is subject to different forces regulating nucleation and propagation steps and provide evidence that the more global β-sheet/β-sheet facial complementarity is a critical determinant for amyloid nucleation and structural selection.

Distinct oligomerization and fibrillization dynamics of amyloid core sequences of amyloid-beta and islet amyloid polypeptide

A direct observation of amyloid aggregation from isolated peptides to cross-β fibrils is crucial for understanding the nucleation-dependence process, but the corresponding macroscopic timescales impose a major computational challenge. Using rapid all-atom discrete molecular dynamics simulations, we capture the oligomerization and fibrillization dynamics of the amyloid core sequences of amyloid-β (Aβ) in Alzheimer's disease and islet amyloid polypeptide (IAPP) in type-2 diabetes, namely Aβ16-22 and IAPP22-28. Both peptides and their mixture spontaneously assemble into cross-β aggregates in silico, but follow distinct pathways. Aβ16-22 is highly aggregation-prone with a funneled free energy basin toward multi-layer β-sheet aggregates. IAPP22-28, on the other hand, features the accumulation of unstructured oligomers before the nucleation of β-sheets and growth into double-layer β-sheet aggregates. In the presence of Aβ16-22, the aggregation of IAPP22-28 is promoted by forming co-aggregated multi-layer β-sheets. Our study offers a detailed molecular insight to the long-postulated oligomerization-nucleation process in the amyloid aggregations.

1,2,3,4,6-penta-O-galloyl-β-D-glucopyranose Binds to the N-terminal Metal Binding Region to Inhibit Amyloid β-protein Oligomer and Fibril Formation

The early oligomerization of amyloid β-protein (Aβ) is a crucial step in the etiology of Alzheimer's disease (AD), in which soluble and highly neurotoxic oligomers are produced and accumulated inside neurons. In search of therapeutic solutions for AD treatment and prevention, potent inhibitors that remodel Aβ assembly and prevent neurotoxic oligomer formation offer a promising approach. In particular, several polyphenolic compounds have shown anti-aggregation properties and good efficacy on inhibiting oligomeric amyloid formation. 1,2,3,4,6-penta-O-galloyl-β-D-glucopyranose is a large polyphenol that has been shown to be effective at inhibiting aggregation of full-length Aβ_1-40 and Aβ_1-42, but has the opposite effect on the C-terminal fragment Aβ_25-35. Here, we use a combination of ion mobility coupled to mass spectrometry (IMS-MS), transmission electron microscopy (TEM) and molecular dynamics (MD) simulations to elucidate the inhibitory effect of PGG on aggregation of full-length Aβ_1-40 and Aβ_1-42. We show that PGG interacts strongly with these two peptides, especially in their N-terminal metal binding regions, and suppresses the formation of Aβ_1-40 tetramer and Aβ_1-42 dodecamer. By exploring multiple facets of polyphenol-amyloid interactions, we provide a molecular basis for the opposing effects of PGG on full-length Aβ and its C-terminal fragments.

Amyloid scaffolds as alternative chlorosomes

Living systems contain remarkable functional capability built within sophisticated self-organizing frameworks. Defining the assembly codes that coordinate these systems could greatly extend nanobiotechnology. To that end, we have highlighted the self-assembling architecture of the chlorosome antenna arrays and report the emulation and extension of their features for the development of cell-compatible photoredox materials. We specifically review work on amyloid peptide scaffolds able to (1) organize light-harvesting chromophores, (2) break peptide bilayer symmetry for directional energy and electron transfer, and (3) incorporate redox active metal ions at high density for energy storage.

Structural and Spectroscopic Properties of Assemblies of Self-Replicating Peptide Macrocycles

Self-replication at the molecular level is often seen as essential to the early origins of life. Recently a mechanism of self-replication has been discovered in which replicator self-assembly drives the process. We have studied one of the examples of such self-assembling self-replicating molecules to a high level of structural detail using a combination of computational and spectroscopic techniques. Molecular Dynamics simulations of self-assembled stacks of peptide-derived replicators provide insights into the structural characteristics of the system and serve as the basis for semiempirical calculations of the UV-vis, circular dichroism (CD) and infrared (IR) absorption spectra that reflect the chiral organization and peptide secondary structure of the stacks. Two proposed structural models are tested by comparing calculated spectra to experimental data from electron microscopy, CD and IR spectroscopy, resulting in a better insight into the specific supramolecular interactions that lead to self-replication. Specifically, we find a cooperative self-assembly process in which β-sheet formation leads to well-organized structures, while also the aromatic core of the macrocycles plays an important role in the stability of the resulting fibers.

Theory of amyloid fibril nucleation from folded proteins

We present a theoretical model for the nucleation of amyloid fibrils. In our model we use helix-coil theory to describe the equilibrium between a soluble native state and an aggregation-prone unfolded state. We then extend the theory to include oligomers with β-sheet cores and calculate the free energy of these states using estimates for the energies of H-bonds, steric zipper interactions, and the conformational entropy cost of forming secondary structure. We find that states with fewer than ~10 β-strands are unstable relative to the dissociated state and three β-strands is the highest free energy state. We then use a modified version of Classical Nucleation Theory to compute the nucleation rate of fibrils from a supersaturated solution of monomers, dimers, and trimers. The nucleation rate has a non-monotonic dependence on denaturant concentration reflecting the competing effects of destabilizing the fibril and increasing the concentration of unfolded monomers. We estimate heterogeneous nucleation rates and discuss the application of our model to secondary nucleation.

Scaleable two-component gelator from phthalic acid derivatives and primary alkyl amines: acid-base interaction in the cooperative assembly

A series of phthalic acid derivatives (P) with a carbon-chain tail was designed and synthesized as single-component gelators. A combination of the single-component gelator P and a non-gelling additive n-alkylamine A through acid-base interaction brought about a series of novel phase-selective two-component gelators PA. The gelation capabilities of P and PA, and the structural, morphological, thermo-dynamic and rheological properties of the corresponding gels were investigated. A molecular dynamics simulation showed that the H-bonding network in PA formed between the NH of A and the carbonyl oxygen of P altered the assembly process of gelator P. Crude PA could be synthesized through a one-step process without any purification and could selectively gel the oil phase without a typical heating-cooling process. Moreover, such a crude PA and its gelation process could be amplified to the kilogram scale with high efficiency, which offers a practical economically viable solution to marine oil-spill recovery.

Lattice model for amyloid peptides: OPEP force field parametrization and applications to the nucleus size of Alzheimer's peptides

Coarse-grained protein lattice models approximate atomistic details and keep the essential interactions. They are, therefore, suitable for capturing generic features of protein folding and amyloid formation at low computational cost. As our aim is to study the critical nucleus sizes of two experimentally well-characterized peptide fragments Aβ16-22 and Aβ37-42 of the full length Aβ1-42 Alzheimer's peptide, it is important that simulations with the lattice model reproduce all-atom simulations. In this study, we present a comprehensive force field parameterization based on the OPEP (Optimized Potential for Efficient protein structure Prediction) force field for an on-lattice protein model, which incorporates explicitly the formation of hydrogen bonds and directions of side-chains. Our bottom-up approach starts with the determination of the best lattice force parameters for the Aβ16-22 dimer by fitting its equilibrium parallel and anti-parallel β-sheet populations to all-atom simulation results. Surprisingly, the calibrated force field is transferable to the trimer of Aβ16-22 and the dimer and trimer of Aβ37-42. Encouraged by this finding, we characterized the free energy landscapes of the two decamers. The dominant structure of the Aβ16-22 decamer matches the microcrystal structure. Pushing the simulations for aggregates between 4-mer and 12-mer suggests a nucleus size for fibril formation of 10 chains. In contrast, the Aβ37-42 decamer is largely disordered with mixed by parallel and antiparallel chains, suggesting that the nucleus size is >10 peptides. Our refined force field coupled to this on-lattice model should provide useful insights into the critical nucleation number associated with neurodegenerative diseases.

Signature of an aggregation-prone conformation of tau

The self-assembly of the microtubule associated tau protein into fibrillar cell inclusions is linked to a number of devastating neurodegenerative disorders collectively known as tauopathies. The mechanism by which tau self-assembles into pathological entities is a matter of much debate, largely due to the lack of direct experimental insights into the earliest stages of aggregation. We present pulsed double electron-electron resonance measurements of two key fibril-forming regions of tau, PHF6 and PHF6*, in transient as aggregation happens. By monitoring the end-to-end distance distribution of these segments as a function of aggregation time, we show that the PHF6^(*) regions dramatically extend to distances commensurate with extended β-strand structures within the earliest stages of aggregation, well before fibril formation. Combined with simulations, our experiments show that the extended β-strand conformational state of PHF6^(*) is readily populated under aggregating conditions, constituting a defining signature of aggregation-prone tau, and as such, a possible target for therapeutic interventions.

Self-Assembly and Anti-Amyloid Cytotoxicity Activity of Amyloid beta Peptide Derivatives

The self-assembly of two derivatives of KLVFF, a fragment Aβ(16-20) of the amyloid beta (Aβ) peptide, is investigated and recovery of viability of neuroblastoma cells exposed to Aβ (1-42) is observed at sub-stoichiometric peptide concentrations. Fluorescence assays show that NH2-KLVFF-CONH2 undergoes hydrophobic collapse and amyloid formation at the same critical aggregation concentration (cac). In contrast, NH2-K(Boc)LVFF-CONH2 undergoes hydrophobic collapse at a low concentration, followed by amyloid formation at a higher cac. These findings are supported by the β-sheet features observed by FTIR. Electrospray ionization mass spectrometry indicates that NH2-K(Boc)LVFF-CONH2 forms a significant population of oligomeric species above the cac. Cryo-TEM, used together with SAXS to determine fibril dimensions, shows that the length and degree of twisting of peptide fibrils seem to be influenced by the net peptide charge. Grazing incidence X-ray scattering from thin peptide films shows features of β-sheet ordering for both peptides, along with evidence for lamellar ordering of NH2-KLVFF-CONH2. This work provides a comprehensive picture of the aggregation properties of these two KLVFF derivatives and shows their utility, in unaggregated form, in restoring the viability of neuroblastoma cells against Aβ-induced toxicity.

The Levinthal Problem in Amyloid Aggregation: Sampling of a Flat Reaction Space

The formation of amyloid fibrils has been associated with many neurodegenerative disorders, yet the mechanism of aggregation remains elusive, partly because aggregation time scales are too long to probe with atomistic simulations. A microscopic theory of fibril elongation was recently developed that could recapitulate experimental results with respect to the effects of temperature, denaturants, and protein concentration on fibril growth kinetics (Schmit, J. D. J. Chem. Phys. 2013, 138 (18), 185102). The theory identifies the conformational search over H-bonding states as the slowest step in the aggregation process and suggests that this search can be efficiently modeled as a random walk on a rugged one-dimensional energy landscape. This insight motivated the multiscale computational algorithm for simulating fibril growth presented in this paper. Briefly, a large number of short atomistic simulations are performed to compute the system diffusion tensor in the reaction coordinate space predicted by the analytic theory. Ensemble aggregation pathways and growth kinetics are then computed from Markov state model (MSM) trajectories. The algorithm is deployed here to understand the fibril growth mechanism and kinetics of Aβ_16-22 and three of its mutants. The order of growth rates of the wild-type and two single mutation peptides (CHA19 and CHA20) predicted by the MSM trajectories is consistent with experimental results. The simulation also correctly predicts that the double mutation (CHA19/CHA20) would reduce the fibril growth rate, even though the degree of rate reduction with respect to either single mutation is overestimated. This artifact may be attributed to the simplistic implicit solvent model. These trends in the growth rate are not apparent from inspection of the rate constants of individual bonds or the lifetimes of the mis-registered states that are the primary kinetic traps but only emerge in the ensemble of trajectories generated by the MSM.