An Account of Amyloid Oligomers: Facts and Figures Obtained from Experiments and Simulations

A turn for the worse: Aβ β-hairpins in Alzheimer's disease

Amyloid-β (Aβ) oligomers are a cause of neurodegeneration in Alzheimer's disease (AD). These soluble aggregates of the Aβ peptide have proven difficult to study due to their inherent metastability and heterogeneity. Strategies to isolate and stabilize homogenous Aβ oligomer populations have emerged such as mutations, covalent cross-linking, and protein fusions. These strategies along with molecular dynamics simulations have provided a variety of proposed structures of Aβ oligomers, many of which consist of molecules of Aβ in β-hairpin conformations. β-Hairpins are intramolecular antiparallel β-sheets composed of two β-strands connected by a loop or turn. Three decades of research suggests that Aβ peptides form several different β-hairpin conformations, some of which are building blocks of toxic Aβ oligomers. The insights from these studies are currently being used to design anti-Aβ antibodies and vaccines to treat AD. Research suggests that antibody therapies designed to target oligomeric Aβ may be more successful at treating AD than antibodies designed to target linear epitopes of Aβ or fibrillar Aβ. Aβ β-hairpins are good epitopes to use in antibody development to selectively target oligomeric Aβ. This review summarizes the research on β-hairpins in Aβ peptides and discusses the relevance of this conformation in AD pathogenesis and drug development.

Perspective for Molecular Dynamics Simulation Studies of Amyloid-β Aggregates

The cause of Alzheimer's disease is related to aggregates such as oligomers and amyloid fibrils consisting of amyloid-β (Aβ) peptides. Molecular dynamics (MD) simulation studies have been conducted to understand the molecular mechanism of the formation and disruption of Aβ aggregates. In this Perspective, the MD simulation studies are classified into four categories, focusing on the target systems: aggregation of Aβ peptides in bulk solution, Aβ aggregation at the interface, aggregation inhibitor against Aβ peptides, and nonequilibrium MD simulation of Aβ aggregates. MD simulation studies in these categories are first reviewed. Future perspectives in each category are then presented. Finally, the overall perspective is presented on how MD simulations of Aβ aggregates can be utilized for developing Alzheimer's disease treatment.

Aggregation of an Amyloidogenic Peptide on Gold Surfaces

Solid surfaces have been shown to affect the aggregation and assembly of many biomolecular systems. One important example is the formation of protein fibrils, which can occur on a range of biological and synthetic surfaces. The rate of fibrillation depends on both the protein structure and the surface chemistry, with the different molecular and oligomer structures adopted by proteins on surfaces likely to be crucial. In this paper, the aggregation of the model amyloidogenic peptide, Aβ(16-22), corresponding to a hydrophobic segment of the amyloid beta protein on a gold surface is studied using molecular dynamics simulation. Previous simulations of this peptide on gold surfaces have shown that it adopts conformations on surfaces that are quite different from those in bulk solution. These simulations show that this then leads to significant differences in the oligomer structures formed in solution and on gold surfaces. In particular, oligomers formed on the surface are low in beta-strands so are unlike the structures formed in bulk solution. When oligomers formed in solution adsorb onto gold surfaces they can then restructure themselves. This can then help explain the inhibition of Aβ(16-22) fibrillation by gold surfaces and nanoparticles seen experimentally.

Differentiation of subnucleus-sized oligomers and nucleation-competent assemblies of the Aβ peptide

A significant feature of Alzheimer's disease is the formation of amyloid deposits in the brain consisting mainly of misfolded derivatives of proteolytic cleavage products of the amyloid precursor protein amyloid-β (Aβ) peptide. While high-resolution structures already exist for both the monomer and the amyloid fibril of the Aβ peptide, the mechanism of amyloid formation itself still defies precise characterization. In this study, low and high molecular weight oligomers (LMWOs and HMWOs) were identified by sedimentation velocity analysis, and for the first time, the temporal evolution of oligomer size distributions was correlated with the kinetics of amyloid formation as determined by thioflavin T-binding studies. LMWOs of subnucleus size contain fewer than seven monomer units and exist alongside a heterogeneous group of HMWOs with 20-160 monomer units that represent potential centers of nucleus formation due to high local monomer concentrations. These HMWOs already have slightly increased β-strand content and appear structurally similar regardless of size, as shown by examination with a range of fluorescent dyes. Once fibril nuclei are formed, the monomer concentration begins to decrease, followed by a decrease in oligomer concentration, starting with LMWOs, which are the least stable species. The observed behavior classifies the two LMWOs as off pathway. In contrast, we consider HMWOs to be on-pathway, prefibrillar intermediates, representing structures in which nucleated conformational conversion is facilitated by high local concentrations. Aβ40 and Aβ42 M35ox take much longer to form nuclei and enter the growth phase than Aβ42 under identical reaction conditions, presumably because both the size and the concentration of HMWOs formed are much smaller.

Understanding the mechanism of amylin aggregation: From identifying crucial segments to tracing dominant sequential events to modeling potential aggregation suppressors

One of the most abundant, prevailing, and life-threatening human diseases that are currently baffling the scientific community is type 2 diabetes (T2D). The self-association of human amylin has been implicated in the pathogenesis of T2D, though with an inconclusive understanding of the mechanism. Hence, we focused on the characterization of the conformational ensembles of all the species that are believed to define the structural polymorphism of the aggregation process - the functional monomeric, the initially self-associated oligomeric, and the structured protofibril - by employing near-equilibrium, non-equilibrium, and equilibrium atomistic simulations on the sporadic, two familial variants (S20G and G33R), and their proline-substituted forms (S20P and G33P). The dynamic near-equilibrium assays hint toward - the abundance of helical conformation in the monomeric state, the retainment of the helicity in the initial self-associated oligomeric phase pointing toward the existence of the helix-helix association mechanism, the difference in preference of specific segments to have definite secondary structural features, the phase-dependent variability in the dominance of specific segments and mutation sites, and the simultaneous presence of generic and unique features among various sequences. Furthermore, the non-equilibrium pulling assays exemplify a generic sequential unzipping mechanism of the protofibrils, however, the sequence-dependent uniqueness comes from the difference in location and magnitude of the control of a specific terminus. Importantly, the equilibrium thermodynamic assays efficiently rank order the potential of aggregability among sequences and consequently suggests the probability of designing effective aggregation suppressors against sporadic and familial amylin variants incorporating proline as the mutation.

Dendrimers as Antiamyloid Agents

Dendrimer–protein conjugates have significant prospects for biological applications. The complexation changes the biophysical behavior of both proteins and dendrimers. The dendrimers could influence the secondary structure of proteins, zeta-potential, distribution of charged regions on the surface, the protein–protein interactions, etc. These changes offer significant possibilities for the application of these features in nanotheranostics and biomedicine. Based on the dendrimer–protein interactions, several therapeutic applications of dendrimers have emerged. Thus, the formation of stable complexes retains the disordered proteins on the aggregation, which is especially important in neurodegenerative diseases. To clarify the origin of these properties and assess the efficiency of action, the mechanism of protein–dendrimer interaction and the nature and driving force of binding are considered in this review. The review outlines the antiamyloid activity of dendrimers and discusses the effect of dendrimer structures and external factors on their antiamyloid properties.

Disaggregation of Islet Amyloid Polypeptide Fibrils as a Potential Anti-Fibrillation Mechanism of Tetrapeptide TNGQ

Islet amyloid polypeptide (IAPP) fibrillation has been commonly associated with the exacerbation of type 2 diabetes prognosis. Consequently, inhibition of IAPP fibrillation to minimize β-cell cytotoxicity is an important approach towards β-cell preservation and type 2 diabetes management. In this study, we identified three tetrapeptides, TNGQ, MANT, and YMSV, that inhibited IAPP fibrillation. Using thioflavin T (ThT) fluorescence assay, circular dichroism (CD) spectroscopy, dynamic light scattering (DLS), and molecular docking, we evaluated the potential anti-fibrillation mechanism of the tetrapeptides. ThT fluorescence kinetics and microscopy as well as transmission electron microscopy showed that TNGQ was the most effective inhibitor based on the absence of normal IAPP fibrillar morphology. CD spectroscopy showed that TNGQ maintained the α-helical conformation of monomeric IAPP, while DLS confirmed the presence of varying fibrillation species. Molecular docking showed that TNGQ and MANT interact with monomeric IAPP mainly by hydrogen bonding and electrostatic interaction, with TNGQ binding at IAPP surface compared to YMSV, which had the highest docking score, but interact mainly through hydrophobic interaction in IAPP core. The highly polar TNGQ was the most active and appeared to inhibit IAPP fibrillation by disaggregation of preformed IAPP fibrils. These findings indicate the potential of TNGQ in the development of peptide-based anti-fibrillation and antidiabetic nutraceuticals.

Molecular insights into the primary nucleation of polymorphic amyloid β dimers in DOPC lipid bilayer membrane

Alzheimer's disease (AD) pathology is characterized by loss of memory cognitive and behavioral deterioration. One of the hallmarks of AD is amyloid β (Aβ) plaques in the brain that consists of Aβ oligomers and fibrils. It is accepted that oligomers, particularly dimers, are toxic species that are produced extracellularly and intracellularly in membranes. It is believed that the disruption of membranes by polymorphic Aβ oligomers is the key for the pathology of AD. This is a first study that investigate the effect of polymorphic “α‐helix/random coil” and “fibril‐like” Aβ dimers on 1,2‐dioleoyl‐sn‐glycero‐3‐phosphocholine (DOPC) membrane. It has been found that the DOPC membrane promotes Aβ_1–42 “fibril‐like” dimers and impedes Aβ_1–42 “α‐helix/random coil” dimers. The N‐termini domains within Aβ_1–42 dimers play a role in Aβ aggregation in membrane milieus. In addition, the aromatic π–π interactions (involving residues F19 and F20 in Aβ_1–42) are the driving forces for the hydrophobic interactions that initiate the primary nucleation of polymorphic Aβ_1–42 dimers within DOPC membrane. Finally, the DOPC bilayer membrane thickness is locally decreased, and it is disrupted by an embedded distinct Aβ_1–42 dimer, due to relatively large contacts between Aβ_1–42 monomers and the DOPC membrane. This study reveals insights into the molecular mechanisms by which polymorphic early‐stage Aβ_1–42 dimers have distinct impacts on DOPC membrane.

Early aggregation mechanism of Aβ16-22 revealed by Markov state models

Aβ_16-22 is believed to have critical role in early aggregation of full length amyloids that are associated with the Alzheimer's disease and can aggregate to form amyloid fibrils. However, the early aggregation mechanism is still unsolved. Here, multiple long-term molecular dynamics simulations combining with Markov state model were used to probe the early oligomerization mechanism of Aβ_16-22 peptides. The identified dimeric form adopted either globular random-coil or extended β-strand like conformations. The observed dimers of these variants shared many overall conformational characteristics but differed in several aspects at detailed level. In all cases, the most common type of secondary structure was intermolecular antiparallel β-sheets. The inter-state transitions were very frequent ranges from few to hundred nanoseconds. More strikingly, those states which contain fraction of β secondary structure and significant amount of extended coiled structures, therefore exposed to the solvent, were majorly participated in aggregation. The assembly of low-energy dimers, in which the peptides form antiparallel β sheets, occurred by multiple pathways with the formation of an obligatory intermediates. We proposed that these states might facilitate the Aβ_16-22 aggregation through a significant component of the conformational selection mechanism, because they might increase the aggregates population by promoting the inter-chain hydrophobic and the hydrogen bond contacts. The formation of early stage antiparallel β sheet structures is critical for oligomerization, and at the same time provided a flat geometry to seed the ordered β-strand packing of the fibrils. Our findings hint at reorganization of this part of the molecule as a potentially critical step in Aβ aggregation and will insight into early oligomerization for large β amyloids.

Site specific NMR characterization of abeta-40 oligomers cross seeded by abeta-42 oligomers

Extracellular accumulation of β amyloid peptides of 40 (Aβ₄₀) and 42 residues (Aβ₄₂) has been considered as one of the hallmarks in the pathology of Alzheimer's disease. In this work, we are able to prepare oligomeric aggregates of Aβ with uniform size and monomorphic structure. Our experimental design is to incubate Aβ peptides in reverse micelles (RMs) so that the peptides could aggregate only through a single nucleation process and the size of the oligomers is confined by the physical dimension of the reverse micelles. The hence obtained Aβ oligomers (AβOs) are 23 nm in diameter and they belong to the category of high molecular-weight (MW) oligomers. The solid-state NMR data revealed that Aβ₄₀Os adopt the structural motif of β-loop-β but the chemical shifts manifested that they may be structurally different from low-MW AβOs and mature fibrils. From the thioflavin-T results, we found that high-MW Aβ₄₂Os can accelerate the fibrillization of Aβ₄₀ monomers. Our protocol allows performing cross-seeding experiments among oligomeric species. By comparing the chemical shifts of Aβ₄₀Os cross seeded by Aβ₄₂Os and those of Aβ₄₀Os prepared in the absence of Aβ₄₂Os, we observed that the chemical states of E11, K16, and E22 were altered, whereas the backbone conformation of the β-sheet region near the C-terminus was structurally invariant. The use of reverse micelles allows hitherto the most detailed characterization of the structural variability of Aβ₄₀Os.

Extracellular accumulation of β amyloid peptides of 40 (Aβ₄₀) and 42 residues (Aβ₄₂) has been considered as one of the hallmarks in the pathology of Alzheimer's disease.

Anatomy and Formation Mechanisms of Early Amyloid-β Oligomers with Lateral Branching: Graph Network Analysis on Large-Scale Simulations

Oligomeric amyloid-β aggregates (AβOs) effectively trigger Alzheimer's disease-related toxicity, generating great interest in understanding their structures and formation mechanisms. However, AβOs are heterogeneous and transient, making their structure and formation difficult to study. Here, we performed graph network analysis of tens of microsecond massive simulations of early amyloid-β (Aβ) aggregations at near-atomic resolution to characterize AβO structures with sizes up to 20-mers. We found that AβOs exhibit highly curvilinear, irregular shapes with occasional lateral branches, consistent with recent cryo-electron tomography experiments. We also found that Aβ40 oligomers were more likely to develop branches than Aβ42 oligomers, explaining an experimental observation that only Aβ40 was trapped in network-like aggregates and exhibited slower fibrillization kinetics. Moreover, AβO architecture dissection revealed that their curvilinear appearance is related to the local packing geometries of neighboring peptides and that Aβ40's greater branching ability originates from specific C-terminal interactions at branching interfaces. Finally, we demonstrate that whether Aβ oligomerization causes oligomers to elongate or to branch depends on the sizes and shapes of colliding aggregates. Collectively, this study provides bottom-up structural information for understanding early Aβ aggregation and AβO toxicity.

Graph network analysis on large-scale simulations uncovers the differential branching behaviours of large Aβ40 and Aβ42 oligomers.

Computer Simulations Aimed at Exploring Protein Aggregation and Dissociation

Protein aggregation can lead to well-defined structures that are functional, but is also the cause of the death of neuron cells in many neurodegenerative diseases. The complexity of the molecular events involved in the aggregation kinetics of amyloid proteins and the transient and heterogeneous characters of all oligomers prevent high-resolution structural experiments. As a result, computer simulations have been used to determine the atomic structures of amyloid proteins at different association stages as well as to understand fibril dissociation. In this chapter, we first review the current computer simulation methods used for aggregation with some atomistic and coarse-grained results aimed at better characterizing the early formed oligomers and amyloid fibril formation. Then we present the applications of non-equilibrium molecular dynamics simulations to comprehend the dissociation of protein assemblies.

Amyloid-β peptide dimers undergo a random coil to β-sheet transition in the aqueous phase but not at the neuronal membrane

Mounting evidence suggests that the neuronal cell membrane is the main site of oligomer-mediated neuronal toxicity of amyloid-β peptides in Alzheimer's disease. To gain a detailed understanding of the mutual interference of amyloid-β oligomers and the neuronal membrane, we carried out microseconds of all-atom molecular dynamics (MD) simulations on the dimerization of amyloid-β (Aβ)42 in the aqueous phase and in the presence of a lipid bilayer mimicking the in vivo composition of neuronal membranes. The dimerization in solution is characterized by a random coil to β-sheet transition that seems on pathway to amyloid aggregation, while the interactions with the neuronal membrane decrease the order of the Aβ42 dimer by attenuating its propensity to form a β-sheet structure. The main lipid interaction partners of Aβ42 are the surface-exposed sugar groups of the gangliosides GM1. As the neurotoxic activity of amyloid oligomers increases with oligomer order, these results suggest that GM1 is neuroprotective against Aβ-mediated toxicity.

Energy landscapes of protein aggregation and conformation switching in intrinsically disordered proteins

The protein folding problem was apparently solved recently by the advent of a deep learning method for protein structure prediction called AlphaFold. However, this program is not able to make predictions about the protein folding pathways. Moreover, it only treats about half of the human proteome, as the remaining proteins are intrinsically disordered or contain disordered regions. By definition these proteins differ from natively folded proteins and do not adopt a properly folded structure in solution. However these intrinsically disordered proteins (IDPs) also systematically differ in amino acid composition and uniquely often become folded upon binding to an interaction partner. These factors preclude solving IDP structures by current machine-learning methods like AlphaFold, which also cannot solve the protein aggregation problem, since this meta-folding process can give rise to different aggregate sizes and structures. An alternative computational method is provided by molecular dynamics simulations that already successfully explored the energy landscapes of IDP conformational switching and protein aggregation in multiple cases. These energy landscapes are very different from those of 'simple' protein folding, where one energy funnel leads to a unique protein structure. Instead, the energy landscapes of IDP conformational switching and protein aggregation feature a number of minima for different competing low-energy structures. In this review, I discuss the characteristics of these multifunneled energy landscapes in detail, illustrated by molecular dynamics simulations that elucidated the underlying conformational transitions and aggregation processes.

Thermodynamics and kinetics of the amyloid-β peptide revealed by Markov state models based on MD data in agreement with experiment

The amlyoid-β peptide (Aβ) is closely linked to the development of Alzheimer's disease. Molecular dynamics (MD) simulations have become an indispensable tool for studying the behavior of this peptide at the atomistic level. General key aspects of MD simulations are the force field used for modeling the peptide and its environment, which is important for accurate modeling of the system of interest, and the length of the simulations, which determines whether or not equilibrium is reached. In this study we address these points by analyzing 30-μs MD simulations acquired for Aβ40 using seven different force fields. We assess the convergence of these simulations based on the convergence of various structural properties and of NMR and fluorescence spectroscopic observables. Moreover, we calculate Markov state models for the different MD simulations, which provide an unprecedented view of the thermodynamics and kinetics of the amyloid-β peptide. This further allows us to provide answers for pertinent questions, like: which force fields are suitable for modeling Aβ? (a99SB-UCB and a99SB-ILDN/TIP4P-D); what does Aβ peptide really look like? (mostly extended and disordered) and; how long does it take MD simulations of Aβ to attain equilibrium? (at least 20-30 μs). We believe the analyses presented in this study will provide a useful reference guide for important questions relating to the structure and dynamics of Aβ in particular, and by extension other similar disordered proteins.

Aggregation of Aβ40/42 chains in the presence of cyclic neuropeptides investigated by molecular dynamics simulations

Alzheimer’s disease is associated with the formation of toxic aggregates of amyloid beta (Aβ) peptides. Despite tremendous efforts, our understanding of the molecular mechanisms of aggregation, as well as cofactors that might influence it, remains incomplete. The small cyclic neuropeptide somatostatin-14 (SST₁₄) was recently found to be the most selectively enriched protein in human frontal lobe extracts that binds Aβ₄₂ aggregates. Furthermore, SST₁₄’s presence was also found to promote the formation of toxic Aβ₄₂ oligomers in vitro. In order to elucidate how SST₁₄ influences the onset of Aβ oligomerization, we performed all-atom molecular dynamics simulations of model mixtures of Aβ₄₂ or Aβ₄₀ peptides with SST₁₄ molecules and analyzed the structure and dynamics of early-stage aggregates. For comparison we also analyzed the aggregation of Aβ₄₂ in the presence of arginine vasopressin (AVP), a different cyclic neuropeptide. We observed the formation of self-assembled aggregates containing the Aβ chains and small cyclic peptides in all mixtures of Aβ₄₂–SST₁₄, Aβ₄₂–AVP, and Aβ₄₀–SST₁₄. The Aβ₄₂–SST₁₄ mixtures were found to develop compact, dynamically stable, but small aggregates with the highest exposure of hydrophobic residues to the solvent. Differences in the morphology and dynamics of aggregates that comprise SST₁₄ or AVP appear to reflect distinct (1) regions of the Aβ chains they interact with; (2) propensities to engage in hydrogen bonds with Aβ peptides; and (3) solvent exposures of hydrophilic and hydrophobic groups. The presence of SST₁₄ was found to impede aggregation in the Aβ₄₂–SST₁₄ system despite a high hydrophobicity, producing a stronger “sticky surface” effect in the aggregates at the onset of Aβ₄₂–SST₁₄ oligomerization.

Improper folding of proteins causes disorders known as protein misfolding diseases. Under normal conditions most proteins adopt particular folds, which allow them functioning properly. However, for reasons that are not yet fully understood, proteins may misfold and aggregate, forming deposits known as amyloid fibrils, which accumulate in the brain or other tissues. This process affects functioning of the nervous system, gradually causing loss of cognitive abilities. Alzheimer’s disease is one of the most common diseases from this group. A better understanding of the aggregation of peptides implicated in Alzheimer’s disease, known as amyloid beta (Aβ) peptides, may facilitate the development of treatments that ameliorate or prevent the disease. We use detailed molecular dynamics simulations to investigate the influence of somatostatin-14 (SST₁₄), a cyclic neuropeptide that might be involved in the amyloidogenic aggregation of Aβ, on molecular processes occurring at the onset of Aβ aggregation. Results of these simulations explain how the presence of SST₁₄ might alter pathways of aggregation of Aβ, shedding light upon the possible role of extrinsic factors in the aggregation at a molecular level.

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.

Real-Time Analysis of Specific Binding between Apolipoprotein E Isoforms and Amyloid β-Peptide by Dual Polarization Interferometry

One of the pathogenesis hypotheses of Alzheimer's disease (AD) is amyloid depositions and neurofibrillary tangles. Apolipoprotein E (Apo E) acts a vital part in the development of AD by affecting the aggregation and clearance of amyloid-β (Aβ). In this paper, a dual polarization interferometry (DPI) technique was employed for a real-time investigation toward the binding events of Apo E isoforms, for instance, Apo E2, Apo E3, and Apo E4, with Aβ_1-40. By evaluation of detailed binding information provided by DPI, the affinities between Apo E isoforms and Aβ_1-40 follow the order of E4 > E3 > E2, and the dissociation constants (K_D) of Aβ_1-40 with Apo E2, Apo E3, and Apo E4 were determined to be 251 ± 37, 40 ± 0.65, and 24.6 ± 2.42 nM, respectively. Our findings reveal the isoform-specific binding behaviors from a kinetics perspective, which can help us understand that Apo E4 has a higher risk of causing AD because of its promoting effect on Aβ aggregation and fibrillation and inefficient clearance of Aβ. Remarkably, this work provides a promising method for exploring the dynamics of interactions between biomolecules and expectantly contributes to the development of AD drugs and therapies targeting Apo E and Aβ.

The F19W mutation reduces the binding affinity of the transmembrane Aβ11-40 trimer to the membrane bilayer

Alzheimer's disease is linked to the aggregation of the amyloid-β protein (Aβ) of 40 or 42 amino acids. Lipid membranes are known to modulate the rate and mechanisms of the Aβ aggregation. Point mutations in Aβ can alter these rates and mechanisms. In particular, experiments show that F19 mutations influence the aggregation rate, but maintain the fibril structures. Here, we used molecular dynamics simulations to examine the effect of the F19W mutation in the 3Aβ11-40 trimer immersed in DPPC lipid bilayers submerged in aqueous solution. Substituting Phe by its closest (non-polar) aromatic amino acid Trp has a dramatic reduction in binding affinity to the phospholipid membrane (measured with respect to the solvated protein) compared to the wild type: the binding free energy of the protein-DPPC lipid bilayer increases by 40-50 kcal mol-1 over the wild-type. This is accompanied by conformational changes and loss of salt bridges, as well as a more complex free energy surface, all indicative of a more flexible and less stable mutated trimer. These results suggest that the impact of mutations can be assessed, at least partially, by evaluating the interaction of the mutated peptides with the lipid membranes.

Protein aggregation detection with fluorescent macromolecular and nanostructured probes: challenges and opportunities

Nanoprobes based on various nanomaterials, polymers or AIEgens are overcoming previous limitations for diagnosis and therapy of early-stage protein aggregation.

Visualizing and trapping transient oligomers in amyloid assembly pathways

Oligomers which form during amyloid fibril assembly are considered to be key contributors towards amyloid disease. However, understanding how such intermediates form, their structure, and mechanisms of toxicity presents significant challenges due to their transient and heterogeneous nature. Here, we discuss two different strategies for addressing these challenges: use of (1) methods capable of detecting lowly-populated species within complex mixtures, such as NMR, single particle methods (including fluorescence and force spectroscopy), and mass spectrometry; and (2) chemical and biological tools to bias the amyloid energy landscape towards specific oligomeric states. While the former methods are well suited to following the kinetics of amyloid assembly and obtaining low-resolution structural information, the latter are capable of producing oligomer samples for high-resolution structural studies and inferring structure-toxicity relationships. Together, these different approaches should enable a clearer picture to be gained of the nature and role of oligomeric intermediates in amyloid formation and disease.

Shear-induced aggregation of amyloid β (1-40) in a parallel plate geometry

Protein aggregation is induced by various environmental or external factors and associated with various neurodegenerative diseases. Among various external factors, shear stress is inevitable for both in vivo and in vitro applications of proteins. In this study, Aβ (1-40) peptide, a derivative of the amyloid precursor protein, was subjected to constant (300, 500, 700 s^-1) and varying (ramp) shear in a parallel plate geometry to explore the implications of shear in terms of macro (viscosity) and micro (secondary structure, morphology) characteristics. Aβ (1-40) solution followed a shear thickening flow behaviour with performance index value 'n' of 2.12. The fibrillation process resulting from the shear force was evaluated in terms of dissipation energy, which was found to exceed the free energy of unfolding. This resulted in the formation of β-sheet rich structures, which were confirmed by CD and FTIR analyses and enhanced Th-T fluorescence. The apparent rate of aggregation (k) was found to increase with the shear rate, and inversely related to the solution viscosity. The maximum k value was 0.21 ± 0.3 min^-1 at 700 s^-1. The molecular weights of aggregates were determined using gel filtration, which were proportionally related to the solution viscosity. The average molecular weights were estimated to be 70, 62 and 52 KDa for samples sheared at 300, 500 and 700 s^-1, respectively. The present study has deciphered the interplay of viscosity, a fluid property, with the aggregation process and its corresponding change in the secondary structures of the peptide. These findings provide useful insights for understanding various proteopathies under shear force.Communicated by Ramaswamy H. Sarma.

Long-range Regulation of Partially Folded Amyloidogenic Peptides

Neurodegeneration involves abnormal aggregation of intrinsically disordered amyloidogenic peptides (IDPs), usually mediated by hydrophobic protein-protein interactions. There is mounting evidence that formation of α-helical intermediates is an early event during self-assembly of amyloid-β42 (Aβ42) and α-synuclein (αS) IDPs in Alzheimer’s and Parkinson’s disease pathogenesis, respectively. However, the driving force behind on-pathway molecular assembly of partially folded helical monomers into helical oligomers assembly remains unknown. Here, we employ extensive molecular dynamics simulations to sample the helical conformational sub-spaces of monomeric peptides of both Aβ42 and αS. Our computed free energies, population shifts, and dynamic cross-correlation network analyses reveal a common feature of long-range intra-peptide modulation of partial helical folds of the amyloidogenic central hydrophobic domains via concerted coupling with their charged terminal tails (N-terminus of Aβ42 and C-terminus of αS). The absence of such inter-domain fluctuations in both fully helical and completely unfolded (disordered) states suggests that long-range coupling regulates the dynamicity of partially folded helices, in both Aβ42 and αS peptides. The inter-domain coupling suggests a form of intra-molecular allosteric regulation of the aggregation trigger in partially folded helical monomers. This approach could be applied to study the broad range of amyloidogenic peptides, which could provide a new path to curbing pathogenic aggregation of partially folded conformers into oligomers, by inhibition of sites far from the hydrophobic core.

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.

A mechanistic hypothesis for the impairment of synaptic plasticity by soluble Aβ oligomers from Alzheimer's brain

It is increasingly accepted that early cognitive impairment in Alzheimer's disease results in considerable part from synaptic dysfunction caused by the accumulation of a range of oligomeric assemblies of amyloid β-protein (Aβ). Most studies have used synthetic Aβ peptides to explore the mechanisms of memory deficits in rodent models, but recent work suggests that Aβ assemblies isolated from human (AD) brain tissue are far more potent and disease-relevant. Although reductionist experiments show Aβ oligomers to impair synaptic plasticity and neuronal viability, the responsible mechanisms are only partly understood. Glutamatergic receptors, GABAergic receptors, nicotinic receptors, insulin receptors, the cellular prion protein, inflammatory mediators, and diverse signaling pathways have all been suggested. Studies using AD brain-derived soluble Aβ oligomers suggest that only certain bioactive forms (principally small, diffusible oligomers) can disrupt synaptic plasticity, including by binding to plasma membranes and changing excitatory-inhibitory balance, perturbing mGluR, PrP, and other neuronal surface proteins, down-regulating glutamate transporters, causing glutamate spillover, and activating extrasynaptic GluN2B-containing NMDA receptors. We synthesize these emerging data into a mechanistic hypothesis for synaptic failure in Alzheimer's disease that can be modified as new knowledge is added and specific therapeutics are developed.

Heat-induced degradation of fibrils: Exponential vs logistic kinetics

The degradation of fibrils under the influence of thermal fluctuations was studied experimentally by various groups around the world. In the first set of experiments, it was shown that the decay of fibril content, which can be measured by the ThT fluorescence assay, obeys a bi-exponential function. In the second series of experiments, it was demonstrated that when the monomers separated from the aggregate are not recyclable, the time dependence of the number of monomers belonging to the dominant cluster is described by a single-exponential function if the fraction of bound chains becomes less than a certain threshold. Note that the time dependence of the fraction of bound chains can be measured by tryptophan fluorescence. To understand these interesting experimental results, we developed a phenomenological theory and performed molecular simulation. According to our theory and simulations using the lattice and all-atom models, the time dependence of bound chains is described by a logistic function, which slowly decreases at short time scales but becomes a single exponential function at large time scales. The results, obtained by using lattice and all-atom simulations, ascertained that the time dependence of the fibril content can be described by a bi-exponential function that decays faster than the logistic function on short time scales. We have uncovered the molecular mechanism for the distinction between the logistic and bi-exponential behavior. Since the dissociation of the chain from the fibrils requires the breaking of a greater number of inter-chain contacts as compared to the breaking of the beta sheet structure, the decrease in the number of connected chains is slower than the fibril content. Therefore, the time dependence of the aggregate size is logistic, while the two-exponential behavior is preserved for the content of fibrils. Our results are in agreement with the results obtained in both sets of experiments.

Combining molecular dynamics simulations and experimental analyses in protein misfolding

The fold of a protein determines its function and its misfolding can result in loss-of-function defects. In addition, for certain proteins their misfolding can lead to gain-of-function toxicities resulting in protein misfolding diseases such as Alzheimer's, Parkinson's, or the prion diseases. In all of these diseases one or more proteins misfold and aggregate into disease-specific assemblies, often in the form of fibrillar amyloid deposits. Most, if not all, protein misfolding diseases share a fundamental molecular mechanism that governs the misfolding and subsequent aggregation. A wide variety of experimental methods have contributed to our knowledge about misfolded protein aggregates, some of which are briefly described in this review. The misfolding mechanism itself is difficult to investigate, as the necessary timescale and resolution of the misfolding events often lie outside of the observable parameter space. Molecular dynamics simulations fill this gap by virtue of their intrinsic, molecular perspective and the step-by-step iterative process that forms the basis of the simulations. This review focuses on molecular dynamics simulations and how they combine with experimental analyses to provide detailed insights into protein misfolding and the ensuing diseases.

Conjugation of RTHLVFFARK to human lysozyme creates a potent multifunctional modulator for Cu2+-mediated amyloid β-protein aggregation and cytotoxicity

Conjugation of alkaline decapeptide (RTHLVFFARK) to lysozyme creates a potent multifunctional modulator (R-hLys) for Cu²⁺-mediated amyloid β-protein aggregation and cytotoxicity.

Conformational dynamics of amyloid-β (16–22) peptide in aqueous ionic liquids

Molecular dynamics simulations of amyloid-β (16–22) peptide dimer in water as well as at two different experimentally studied concentrations of hydrated ionic liquids (ILs), ethylammonium mesylate (EAM), ethylammonium nitrate (EAN), and triethylammonium mesylate (TEAM), were carried out employing an umbrella sampling method. We used the average Ψ angle of the peptide backbone as the reaction coordinate to observe the conformational changes of a peptide dimer. Secondary structural element values were calculated for the peptide dimer along the reaction coordinate to see the transition of the peptide dimer between β-sheet and α-helix conformations. We observe the β-sheet conformation as the global minimum on the free energy surfaces in both EAM and EAN ILs at both the concentrations and at a low concentration of TEAM. However, we observe α-helix conformation as the global minimum at a high concentration of TEAM. Our results are in good correlation with the experimental findings. We calculated the average number of intramolecular and intermolecular hydrogen bonds of α-helix and β-sheet conformations in all solutions, and they are in correlation with the secondary structure element values. To understand the peptide–IL interactions, atom–atom radial distribution functions of cation, anion, and water around amide oxygen and hydrogen atoms were calculated. The solvent-accessible surface area of the peptide dimer was calculated to understand the exposure of the peptide towards the solvent during conformational changes. Finally, van der Waals (vdW) and Coulomb interaction energies were calculated between peptide–cation, peptide–anion, and peptide–water to understand the stability of conformations in different concentrations. We find that the TEA cation has more vdW interaction energy compared to Coulomb interaction energy with peptide in 70% (w/w) TEAM, which mimics a membrane-like environment to induce α-helix conformation rather than β-sheet conformation.

Molecular dynamics simulations of amyloid-β (16–22) peptide dimer at two different experimentally studied concentrations of hydrated ethylammonium mesylate, ethylammonium nitrate, and triethylammonium mesylate were carried out employing an umbrella sampling method.

Effects of in vivo conditions on amyloid aggregation

One of the grand challenges of biophysical chemistry is to understand the principles that govern protein misfolding and aggregation, which is a highly complex process that is sensitive to initial conditions, operates on a huge range of length- and timescales, and has products that range from protein dimers to macroscopic amyloid fibrils. Aberrant aggregation is associated with more than 25 diseases, which include Alzheimer's, Parkinson's, Huntington's, and type II diabetes. Amyloid aggregation has been extensively studied in the test tube, therefore under conditions that are far from physiological relevance. Hence, there is dire need to extend these investigations to in vivo conditions where amyloid formation is affected by a myriad of biochemical interactions. As a hallmark of neurodegenerative diseases, these interactions need to be understood in detail to develop novel therapeutic interventions, as millions of people globally suffer from neurodegenerative disorders and type II diabetes. The aim of this review is to document the progress in the research on amyloid formation from a physicochemical perspective with a special focus on the physiological factors influencing the aggregation of the amyloid-β peptide, the islet amyloid polypeptide, α-synuclein, and the hungingtin protein.

The interactions of an Aβ protofibril with a cholesterol-enriched membrane and involvement of neuroprotective carbazolium-based substances

Recent studies have shown that the aggregation of the amyloid-beta peptide (Aβ) in the brain cell membrane is responsible for the emergence of Alzheimer's disease (AD); the exploration of effective factors involved in the extension of the aggregation process and alternatively the examination of an effective inhibitor via theoretical and experimental tools are among the main research topics in the field of AD treatment. Therefore, in this study, we used all-atom molecular dynamics (MD) simulations to clarify the impact of cell membrane cholesterol on the interaction of Aβ with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) as a membrane model. Moreover, the effect of the P7C3-S243 molecule on the abovementioned process was investigated. The simulation results disclosed the neuroprotective property of the P7C3-S243 molecule. The MD simulation results indicate that the interaction of cholesterol molecules with the Aβ oligomer is negligible and cannot enhance membrane rupture. However, strong hydrogen bonding between the POPC molecules and the oligomers led to membrane perturbation. According to our modellings, the P7C3-S243 molecular layer can protect the cell membrane by inhibiting the direct interaction between the bilayer and Aβ. In addition, free-energy calculations were conducted to determine the possible penetration of Aβ fibrils into the cholesterol-enriched membrane.

Network-Based Classification and Modeling of Amyloid Fibrils

Amyloid fibrils are locally ordered protein aggregates that self-assemble under a variety of physiological and in vitro conditions. Their formation is of fundamental interest as a physical chemistry problem and plays a central role in Alzheimer's disease, Type II diabetes, and other human diseases. As the number of known amyloid fibril structures has grown, the need has arisen for a nomenclature for describing and classifying fibril types, as well as a theoretical description of the physics that gives rise to the self-assembly of these structures. Here, we introduce a systematic nomenclature and coarse-graining methodology for describing the topology of fibrils and other protein aggregates, along with a computational methodology for simulating protein aggregation. Both have mathematical underpinnings in graph theory and statistical mechanics and are consistent with available experimental data on the fibril structure and aggregation kinetics. Our graph representation of the fibril topology enables us to define a network Hamiltonian based on connectivity patterns among monomers rather than detailed intermolecular interactions, greatly speeding up the simulation of large ensembles. Our simulation strategy is capable of recapitulating the formation of all currently known amyloid fibril topologies found in the Protein Data Bank, as well as the formation kinetics of fibrils and oligomers.

Exploring Amyloid-β Dimer Structure Using Molecular Dynamics Simulations

A major hallmark of Alzheimer's disease (AD) is the aggregation of amyloid-β peptides in the brains of people afflicted by the disease. The exact pathway to this catastrophic event is unknown. In this work, a total of 9.5 μs molecular dynamics simulations have been performed to investigate the structure and dynamics of the smallest form of toxic Aβ oligomers, i.e., the Aβ dimers. This study suggests that specific hydrophobic regions are vital in the aggregation process. Different possible structures for Aβ dimers are reported along with their relative binding affinity. These data may be used to design better Aβ-aggregation inhibitors. The diversity of the dimer structures suggests several aggregation pathways.

d-Enantiomeric RTHLVFFARK-NH2: A Potent Multifunctional Decapeptide Inhibiting Cu2+-Mediated Amyloid β-Protein Aggregation and Remodeling Cu2+-Mediated Amyloid β Aggregates

Aggregation of amyloid β-protein (Aβ) into β-sheet-rich plaques is a general feature of Alzheimer's disease (AD). Homeostasis dysregulation of Cu²⁺ mediates Aβ to form high cytotoxic aggregates, which causes cell damage by generation of reactive oxygen species (ROS). To improve the inhibitory potency and explore the multifaceted functions of our previously designed decapeptide, RTHLVFFARK-NH₂ (RK10), we have herein reformulated the decapeptide into its d-enantiomer, rk10, and the effects of chirality on Aβ aggregation, Cu²⁺-mediated Aβ aggregations, and aggregate-remodeling effects were investigated. The results revealed the following: (1) The d-enantiomer presented enhanced inhibitory potency on Aβ fibrillogenesis in comparison to RK10; rk10 and RK10 increased the cell viability from 60% to 91% and 71%, respectively, at an equimolar concentration to Aβ. (2) The enantiomers were chemically equivalent to Cu²⁺ chelation, ROS suppression and oxidative damage rescue. (3) The d-enantiomer exhibited higher performance to inhibit Cu²⁺-mediated Aβ aggregation, and more significantly attenuated the cytotoxicity caused by Aβ₄₂-Cu²⁺ complex than RK10. Cell viability was rescued from 51% to 89% and 74% by coincubating with rk10 and RK10 at 50 μM, respectively. Intracellular ROS levels generated by Aβ₄₂ and Aβ₄₂-Cu²⁺ species were also remarkably decreased by treating with rk10. (4) The enantiomers could remodel mature Aβ₄₂-Cu²⁺ aggregates by Cu²⁺ chelation, and rk10 showed higher performance than RK10, as evidenced by the enhanced cell viability from 57% to 86% by RK10 and to 96% by rk10. The d-enantiomer also showed higher ability than RK10 on protecting the disrupted species from reaggregation. Taken together, D-chiral derivatization of the decapeptide resulted in a potent multifunctional agent in inhibiting Cu²⁺-mediated Aβ aggregation and remodeling mature Aβ-Cu²⁺ species. To the best of our knowledge, this is the first investigation on the chirality effect of a multifunctional peptide inhibitor on Cu²⁺-mediated Aβ aggregation and on the remodeling effect of mature Aβ-Cu²⁺ aggregates. The work provides new insights into the critical role of chirality in the multifaceted functions of peptide inhibitors against amyloid formation and its toxicity.

Obstacles to translating the promise of nanoparticles into viable amyloid disease therapeutics

Nanoparticles (NPs) constitute a powerful therapeutic platform with exciting prospects as potential inhibitors of amyloid-[Formula: see text] (Aβ) aggregation, a process associated with Alzheimer's disease (AD). Researchers have synthesized and tested a large collection of NPs with disparate sizes, shapes, electrostatic properties and surface ligands that evoke a variety of responses on Aβ aggregation. In spite of a decade of research on the NP-Aβ system and many promising experimental results, NPs have failed to progress to any level of clinical trials for AD. A theoretical framework with which to approach this physical system is presented featuring two simple metrics, (1) the extent to which NPs adsorb Aβ, and (2) the degree to which interaction with a NP alters Aβ conformation relative to aggregation propensity. Most of our current understanding of these two interactions has been gained through experimentation, and many of these studies are reviewed herein. We also provide a potential roadmap for studies that we believe could produce viable NPs as an effective AD therapeutic platform.

Basified Human Lysozyme: A Potent Inhibitor against Amyloid β-Protein Fibrillogenesis

The aggregation of amyloid β-proteins (Aβ) has been recognized as a key process in the pathogenesis of Alzheimer's disease (AD), so inhibiting Aβ aggregation is an important strategy to prevent the onset and treatment of AD. Our recent work indicated that decreasing the positive charges (or introducing negative charges) on human lysozyme (hLys) was unfavorable in keeping the inhibiting capability of hLys on Aβ aggregation. Therefore, we have herein proposed to basify hLys by conversion of the carboxyl groups into amino groups by modification with ethylene diamine. Basified hLys (Lys-B) preparations of three modification degrees (MDs), denoted as hLys-B1 (MD, 1.5), hLys-B2 (MD, 3.3), and hLys-B3 (MD, 4.4), were synthesized for modulating Aβ fibrillogenesis. The hLys-B preparations kept the stability and biocompatibility as native hLys did, whereas the inhibitory potency of hLys-B on Aβ fibrillogenesis increased with increasing MD. Cytotoxicity analysis showed that cell viability with 2.5 μM hLys-B3 increased from 62.5% (with 25 μM Aβ only) to 76.1%, similar to the case with 12.5 μM hLys (75.5%); cell viability with 6.25 μM hLys-B3 increased to 82.0%, similar to the case with 25 μM hLys (80.9%). The results indicate about four- to fivefold increase in the inhibition efficiency of hLys by the amino modification. Mechanistic analysis suggests that such a superior inhibitory capability of hLys-B was attributed to its more widely distributed positive charges, which promoted broad electrostatic interactions between Aβ and hLys-B. Thus, hLys-B suppressed the conformational transition of Aβ to β-sheet structures at low concentrations (e.g., 2.5 μM hLys-B3), leading to changes in the aggregation pathway and the formation of Aβ species with less cytotoxicity. The findings provided new insights into the development of more potent protein-based inhibitors against Aβ fibrillogenesis.

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.

High pressure NMR reveals conformational perturbations by disease-causing mutations in amyloid β-peptide

Here we present the high pressure NMR characterization of Aβ42 and two Aβ40 variants with Alzheimer-causing mutations E22G and D23N. While chemical shifts only identified localized changes at ambient pressure compared with Aβ40, high pressure NMR revealed a common site with heightened pressure sensitivity at Q15, K16 and L17 in all three variants, which correlates to higher β-propensity at central hydrophobic cluster (CHC) and faster aggregation.