Amyloid β Protein and Alzheimer's Disease: When Computer Simulations Complement Experimental Studies

Mitigating Alzheimer’s Disease with Natural Polyphenols: A Review

:

According to Alzheimer’s Disease International (ADI), nearly 50 million people worldwide were living with dementia in 2017, and this number is expected to triple by 2050. Despite years of research in this field, the root cause and mechanisms responsible for Alzheimer’s disease (AD) have not been fully elucidated yet. Moreover, promising preclinical results have repeatedly failed to translate into patient treatments. Until now, none of the molecules targeting AD has successfully passed the Phase III trial. Although natural molecules have been extensively studied, they normally require high concentrations to be effective; alternately, they are too large to cross the blood-brain barrier (BBB).

:

In this review, we report AD treatment strategies, with a virtually exclusive focus on green chemistry (natural phenolic molecules). These include therapeutic strategies for decreasing amyloid-β (Aβ) production, preventing and/or altering Aβ aggregation, and reducing oligomers cytotoxicity such as curcumin, (-)-epigallocatechin-3-gallate (EGCG), morin, resveratrol, tannic acid, and other natural green molecules. We also examine whether consideration should be given to potential candidates used outside of medicine and nutrition, through a discussion of two intermediate-sized green molecules, with very similar molecular structures and key properties, which exhibit potential in mitigating Alzheimer’s disease.

Aggregation rate of amyloid beta peptide is controlled by beta-content in monomeric state

Understanding the key factors that govern the rate of protein aggregation is of immense interest since protein aggregation is associated with a number of neurodegenerative diseases. Previous experimental and theoretical studies have revealed that the hydrophobicity, charge, and population of the fibril-prone monomeric state control the fibril formation rate. Because the fibril structures consist of cross beta sheets, it is widely believed that those sequences that have a high beta content (β) in the monomeric state should have high aggregation rates as the monomer can serve as a template for fibril growth. However, this important fact has never been explicitly proven, motivating us to carry out this study. Using replica exchange molecular dynamics simulation with implicit water, we have computed β of 19 mutations of amyloid beta peptide of 42 residues (Aβ42) for which the aggregation rate κ has been measured experimentally. We have found that κ depends on β in such a way that the higher the propensity to aggregation, the higher the beta content in the monomeric state. Thus, we have solved a long-standing problem of the dependence of fibril formation time of the β-structure on a quantitative level.

Long Noncoding RNA NEAT1 Aggravates Aβ-Induced Neuronal Damage by Targeting miR-107 in Alzheimer's Disease

PURPOSE

Alzheimer's disease (AD) is the most common neurodegenerative disease, with a rising prevalence worldwide. Long noncoding RNAs (lncRNAs) have been found to play important roles in the development and treatment of AD. However, the exact role of lncRNA nuclear enriched abundant transcript 1 (NEAT1) in neuronal damage in AD is largely unknown.

MATERIALS AND METHODS

The AD model was established in SH-SY5Y and SK-N-SH cells via treatment with amyloid β1-42 (Aβ). The expression of NEAT1 and microRNA-107 (miR-107) was measured by quantitative real-time polymerase chain reaction. Cell viability and apoptosis were detected by MTT assay, immunocytochemistry, and flow cytometry. The expression of phosphorylated tau protein (p-Tau) was measured by Western blot. The interaction between NEAT1 and miR-107 was explored by bioinformatics analysis, luciferase activity, and RNA immunoprecipitation assays.

RESULTS

NEAT1 expression was enhanced in Aβ-treated SH-SY5Y and SK-N-SH cells, and its knockdown attenuated Aβ-induced inhibition of viability and promotion of apoptosis and p-Tau levels. NEAT1 was indicated as a decoy of miR-107. miR-107 abundance was reduced in Aβ-treated cells, and its overexpression reversed Aβ-induced injury. Moreover, interference of miR-107 abated silencing of NEAT1-mediated inhibition of neuronal damage in Aβ-treated SH-SY5Y and SK-N-SH cells.

CONCLUSION

LncRNA NEAT1 aggravated Aβ-induced neuronal damage by sponging miR-107, indicating a novel avenue for treatment of AD.

C-Terminal Plays as the Possible Nucleation of the Self-Aggregation of the S-Shape Aβ11-42 Tetramer in Solution: Intensive MD Study

Amyloid beta (Aβ) peptides are characterized as the major factors associated with neuron death in Alzheimer's disease, which is listed as the most common form of neurodegeneration. Disordered Aβ peptides are released from proteolysis of the amyloid precursor protein. The Aβ self-assembly process roughly takes place via five steps: disordered forms → oligomers → photofibrils → mature fibrils → plaques. Although Aβ fibrils are often observed in patient brains, oligomers were recently indicated to be major neurotoxic elements. In this work, the neurotoxic compound S-shape Aβ11-42 tetramer (S4Aβ11-42) was investigated over 10 μs of unbiased MD simulations. In particular, the S4Aβ11-42 oligomer adopted a high dynamics structure, resulting in unsuccessful determination of their structures in experiments. The C-terminal was suggested as the possible nucleation of the Aβ42 aggregation. The sequences 27-35 and 39-40 formed rich β-content, whereas other residues mostly adopted coil structures. The mean value of the β-content over the equilibrium interval is ∼42 ± 3%. Furthermore, the dissociation free energy of the S4Aβ11-42 peptide was predicted using a biased sampling method. The obtained free energy is ΔG US = -58.44 kcal/mol which is roughly the same level as the corresponding value of the U-shape Aβ17-42 peptide. We anticipate that the obtained S4Aβ11-42 structures could be used as targets for AD inhibitor screening over the in silico study.

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

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

Steady, Symmetric, and Reversible Growth and Dissolution of Individual Amyloid-β Fibrils

Oligomers and fibrils of the amyloid-β (Aβ) peptide are implicated in the pathology of Alzheimer's disease. Here, we monitor the growth of individual Aβ40 fibrils by time-resolved in situ atomic force microscopy and thereby directly measure fibril growth rates. The measured growth rates in a population of fibrils that includes both single protofilaments and bundles of filaments are independent of the fibril thickness, indicating that cooperation between adjacent protofilaments does not affect incorporation of monomers. The opposite ends of individual fibrils grow at similar rates. In contrast to the "stop-and-go" kinetics that has previously been observed for amyloid-forming peptides, growth and dissolution of the Aβ40 fibrils are relatively steady for peptide concentration of 0-10 μM. The fibrils readily dissolve in quiescent peptide-free solutions at a rate that is consistent with the microscopic reversibility of growth and dissolution. Importantly, the bimolecular rate coefficient for the association of a monomer to the fibril end is significantly smaller than the diffusion limit, implying that the transition state for incorporation of a monomer into a fibril is associated with a relatively high free energy.

Molecular Simulations Reveal Terminal Group Mediated Stabilization of Helical Conformers in Both Amyloid-β42 and α-Synuclein

The presence of partially structured helices in natively unfolded amyloid-β42 (Aβ42) and α-synuclein (αS) has been shown to accelerate fibrillation in the onset of Alzheimer's and Parkinson's disease, respectively. At the other extreme, folded stable helical conformers have also been reported to resist amyloid formation. Recent studies indicate that amyloidogenic aggregation can be impeded using small molecules that stabilize the α-helical monomers and switch off the neurotoxic pathway. We predict a common intrapeptide route to stabilization based on the plasticity of helical conformations of Aβ42 and αS as assessed through extensive atomistic molecular dynamics (MD) computer simulations (∼36 μs) across ten distinct protein force field and water model combinations. Computed free energies and interaction maps (not obtainable from experiments alone) show that flexible terminal groups (N-terminus of Aβ42 and C-terminus of αS) show a tendency to stabilize folded helical conformations in both peptides via primary hydrophobic interactions with central hydrophobic domains, and secondary salt bridges with other domains. These interactions confer aggregation resistance by decreasing the population of partially structured helices and are absent in control simulations of complete unfolding. Computed helical stability is also significantly reduced in terminal-deleted variants. The models suggest new strategies to tackle neurodegeneration by rationally re-engineering terminal groups to optimize their predicted ability to deactivate helical monomers.

The role of the protein-RNA recognition code in neurodegeneration

MicroRNAs are small endogenous RNAs that pair and bind to sites on mRNAs to direct post-transcriptional repression. However, there is a possibility that microRNAs directly influence protein structure and activity, and this influence can be termed post-translational riboregulation. This conceptual review explores the literature on neurodegenerative disorders. Research on the association between neurodegeneration and RNA-repeat toxicity provides data that support a protein-RNA recognition code. For example, this code explains why hnRNP H and SFPQ proteins, which are involved in amyotrophic lateral sclerosis, are sequestered by the (GGGGCC)n repeat sequence. Similarly, it explains why MNBL proteins and (CTG)n repeats in RNA, which are involved in myotonic dystrophy, are sequestered into RNA foci. Using this code, proteins involved in diseases can be identified. A simple protein BLAST search of the human genome for amino acid repeats that correspond to the nucleotide repeats reveals new proteins among already known proteins that are involved in diseases. For example, the (CAG)n repeat sequence, when transcribed into possible peptide sequences, leads to the identification of PTCD3, Rem2, MESP2, SYPL2, WDR33, COL23A1, and others. After confirming this approach on RNA repeats, in the next step, the code was used in the opposite manner. Proteins that are involved in diseases were compared with microRNAs involved in those diseases. For example, a reasonable correspondence of microRNA 9 and 107 with amyloid-β-peptide (Aβ42) was identified. In the last step, a miRBase search for micro-nucleotides, obtained by transcription of a prion amino acid sequence, revealed new microRNAs and microRNAs that have previously been identified as involved in prion diseases. This concept provides a useful key for designing RNA or peptide probes.

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.

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.

Tautomerization Effect of Histidines on Oligomer Aggregation of β-Amyloid(1-40/42) during the Early Stage: Tautomerism Hypothesis for Misfolding Protein Aggregation

As the intrinsic origin of the hypothesis for β-amyloid (Aβ) from Alzheimer's disease, histidine behaviors were found to play a crucial role in Aβ aggregation. To investigate the histidine behaviors during the early stage of aggregation, Aβ40/42 pentamers with different histidine isomer states were simulated at the atomic level. Results show that five Aβ40 (δδδ) and Aβ42 (εδδ) monomers can rapidly decrease the aggregation threshold, promote stable pentamer formation, and increase pentamer contents by 51.8% and 56.7%, respectively, as compared with the values of their wild-type (εεε) counterparts. Additionally, pentamers of Aβ40 (δδδ) and Aβ42 (εδδ) have different aggregation pathways and disassembly species, Tr+D and Te+M, during the growth of the pentamer. This work discloses the significance of histidine tautomerization in Aβ aggregation, implying a potential way to control Aβ aggregation and develop the assembly inhibitors.

Electric Field Disruption of Amyloid Aggregation: Potential Noninvasive Therapy for Alzheimer's Disease

The aggregation of β-amyloid peptides is a key event in the formative stages of Alzheimer's disease. Promoting folding and inhibiting aggregation was reported as an effective strategy in reducing Aβ-elicited toxicity. This study experimentally investigates the influence of the external electric field (EF) and magnetic field (MF) of varying strengths on the in vitro fibrillogenesis of hydrophobic core sequence, Aβ_16-22, and its parent peptide, Aβ_1-42. Biophysical methods such as ThT fluorescence, static light scattering, circular dichroism, and infrared spectroscopy suggest that EF has a stabilizing effect on the secondary structure, initiating a conformational switch of Aβ_16-22 and Aβ_1-42 from β to non-β conformation. This observation was further corroborated by dynamic light scattering and transmission electron microscopic studies. To mimic in vivo conditions, we repeated ThT fluorescence assay with Aβ_1-42 in human cerebrospinal fluid to verify EF-mediated modulation. The self-seeding of Aβ_1-42 and cross-seeding with Aβ_1-40 to verify that the autocatalytic amplification of self-assembly as a result of secondary nucleation also yields comparable results in EF-exposed and unexposed samples. Aβ-elicited toxicity of EF-treated samples in two neuroblastoma cell lines (SH-SY5Y and IMR-32) and human embryonic kidney cell line (HEK293) were found to be 15-38% less toxic than the EF untreated ones under identical conditions. Experiments with fluorescent labeled Aβ_1-42 to correlate reduced cytotoxicity and cell internalization suggest a comparatively smaller uptake of the EF-treated peptides. Our results provide a scientific roadmap for future noninvasive, therapeutic solutions for the treatment of Alzheimer's disease.

Adequate prediction for inhibitor affinity of Aβ40 protofibril using the linear interaction energy method

The search for efficient inhibitors targeting Aβ oligomers and fibrils is an important issue in Alzheimer's disease treatment. As a consequence, an accurate and computationally cheap approach to estimate the binding affinity for many ligands interacting with Aβ peptides is very important. Here, the calculated binding free energies of 30 ligands interacting with 12Aβ11-40 peptides using the linear interaction energy (LIE) approach are found to be in good correlation with experimental data (R = 0.79). The binding affinities of these complexes are also calculated by using free energy perturbation (FEP) and molecular mechanic/Poisson-Boltzmann surface area (MM/PBSA) methods. The time-consuming FEP method provides results with similar correlation (R = 0.72), whereas MM/PBSA calculations show very low correlation with experimental data (R = 0.27). In all complexes, van der Waals interactions contribute much more than electrostatic interactions. The LIE model, which is much less time-consuming than both the FEP and MM/PBSA methods, opens the door to accurate and rapid affinity prediction of ligands with Aβ peptides and the design of new ligands.

Methionine Oxidation Changes the Mechanism of Aβ Peptide Binding to the DMPC Bilayer

Using all-atom explicit solvent replica exchange molecular dynamics simulations with solute tempering, we study the effect of methionine oxidation on Aβ10–40 peptide binding to the zwitterionic DMPC bilayer. By comparing oxidized and reduced peptides, we identified changes in the binding mechanism caused by this modification. First, Met35 oxidation unravels C-terminal helix in the bound peptides. Second, oxidation destabilizes intrapeptide interactions and expands bound peptides. We explain these outcomes by the loss of amphiphilic character of the C-terminal helix due to oxidation. Third, oxidation “polarizes” Aβ binding to the DMPC bilayer by strengthening the interactions of the C-terminus with lipids while largely releasing the rest of the peptide from bilayer. Fourth, in contrast to the wild-type peptide, oxidized Aβ induces significantly smaller bilayer thinning and drop in lipid density within the binding footprint. These observations are the consequence of mixing oxidized peptide amino acids with lipids promoted by enhanced Aβ conformational fluctuations. Fifth, methionine oxidation reduces the affinity of Aβ binding to the DMPC bilayer by disrupting favorable intrapeptide interactions upon binding, which offset the gains from better hydration. Reduced binding affinity of the oxidized Aβ may represent the molecular basis for its reduced cytotoxicity.

Mechanisms of Metabolite Amyloid Formation: Computational Studies for Drug Design against Metabolic Disorders

Ordered self-organization of polypeptides into fibrillar assemblies has been associated with a number of pathological conditions linked to degenerative diseases. Recent experimental observations have demonstrated that even small-molecule metabolites can aggregate into supramolecular arrangements with structural and functional properties reminiscent of peptide-based amyloids. The molecular determinants of such mechanisms, however, are not clear yet. Herein, we examine the process of formation of ordered aggregates by adenine in aqueous solution by molecular dynamics simulations. We also investigate the effects of an inhibiting polyphenol, namely, epigallocatechin gallate (EGCG), on this mechanism. We show that, while adenine alone is able to form extended amyloid-like oligomers, EGCG interferes with the supramolecular organization process. Interestingly, acetylsalicylic acid is shown not to interfere with ordered aggregation, consistent with experiments. The results of these mechanistic studies indicate the main pharmacophoric determinants that a drug-like inhibitor should possess to effectively interfere with metabolite amyloid formation.

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.

Conformational ensemble of native α-synuclein in solution as determined by short-distance crosslinking constraint-guided discrete molecular dynamics simulations

Combining structural proteomics experimental data with computational methods is a powerful tool for protein structure prediction. Here, we apply a recently-developed approach for de novo protein structure determination based on the incorporation of short-distance crosslinking data as constraints in discrete molecular dynamics simulations (CL-DMD) for the determination of conformational ensemble of the intrinsically disordered protein α-synuclein in the solution. The predicted structures were in agreement with hydrogen-deuterium exchange, circular dichroism, surface modification, and long-distance crosslinking data. We found that α-synuclein is present in solution as an ensemble of rather compact globular conformations with distinct topology and inter-residue contacts, which is well-represented by movements of the large loops and formation of few transient secondary structure elements. Non-amyloid component and C-terminal regions were consistently found to contain β-structure elements and hairpins.

As the population ages, neurodegenerative diseases such as Parkinson’s disease will become an increasing problem in many countries. Aggregation of the protein α-synuclein is the primary cause of Parkinson’s disease, but there is still a dearth of structural information pertaining to the native, non-aggregating form of this protein. A better understanding the structural state of the native protein may prove useful for the design of new therapeutics to combat this disease. In order to obtain more structural information on this protein, we have recently modelled the native α-synuclein protein. These models were generated using a novel approach which combines protein crosslinking and discrete molecular dynamics simulations. We have found that the α-synuclein protein can adopt several shapes, all with a similar topology, resembling a three fingered closed claw. A region of the protein important for aggregation was found to be protected from the surrounding biological environment in these conformations, and the stabilization of these structures may be a fruitful avenue for future drug research into mitigating the cause and effect of Parkinson’s disease.

Free Energy Profile for Penetration of Pittsburgh Compound-B into the Amyloid β Fibril

The amyloid β (Aβ) fibril is a hallmark of Alzheimer's disease (AD) and has therefore served as an important target for early diagnosis of AD. The Pittsburgh Compound-B (PiB) is one of the most famous positron emission tomography (PET) tracers commonly used for in vivo detection of Aβ fibrils. Many theoretical studies have predicted the existence of various core binding sites with different microenvironments for probes binding to the Aβ fibril. However, little attention has been devoted to how the probes actually penetrate into the different core binding sites. In this study, an integrated molecular modeling scheme is used to study the penetration of PiB into the core binding sites of the Aβ_1-42 fibril structure recently obtained by cryogenic electron microscopy. We find that there are two core binding sites for PiB with dramatic differences in cavity size and microenvironment properties, and furthermore that the penetration of PiB into site-1 is energetically prohibitive, whereas the penetration into site-2 is much more favorable. Therefore, the binding capacity at site-2 may be larger than that at site-1 despite its lower binding affinity. Our results thus suggest that site-2 may be a major binding site for PiB binding to Aβ fibril and emphasize the importance to adopt a full dynamical picture when studying tracer-fibril binding problems in general, something that in turn can be used to guide the development of tracers with higher affinity and selectivity for the Aβ fibril.

Norepinephrine Inhibits Alzheimer's Amyloid-β Peptide Aggregation and Destabilizes Amyloid-β Protofibrils: A Molecular Dynamics Simulation Study

The abnormal self-assembly of amyloid-β (Aβ) peptides into toxic fibrillar aggregates is associated with the pathogenesis of Alzheimer's disease (AD). The inhibition of β-sheet-rich oligomer formation is considered as the primary therapeutic strategy for AD. Previous experimental studies reported that norepinephrine (NE), one of the neurotransmitters, is able to inhibit Aβ aggregation and disaggregate the preformed fibrils. Moreover, exercise can markedly increase the level of NE. However, the underlying inhibitory and disruptive mechanisms remain elusive. In this work, we performed extensive replica-exchange molecular dynamic (REMD) simulations to investigate the conformational ensemble of Aβ_1-42 dimer with and without NE molecules. Our results show that without NE molecules, Aβ_1-42 dimer transiently adopts a β-hairpin-containing structure, and the β-strand regions of this β-hairpin (residues 15QKLVFFA21 and 33GLMVGGVV40) strongly resemble those of the Aβ fibril structure (residues 15QKLVFFA21 and 30AIIGLMVG37) reported in an electron paramagnetic resonance spectroscopy study. NE molecules greatly reduce the interpeptide β-sheet content and suppress the formation of the above-mentioned β-hairpin, leading to a more disordered coil-rich Aβ dimer. Five dominant binding sites are identified, and the central hydrophobic core 16KLVFFA21 site and C-terminal 31IIGLMV36 hydrophobic site are the two most favorable ones. Our data reveal that hydrophobic, aromatic stacking, hydrogen-bonding and cation-π interactions synergistically contribute to the binding of NE molecules to Aβ peptides. MD simulations of Aβ_1-42 protofibril show that NE molecules destabilize Aβ protofibril by forming H-bonds with residues D1, A2, D23, and A42. This work reveals the molecular mechanism by which NE molecules inhibit Aβ_1-42 aggregation and disaggregate Aβ protofibrils, providing valuable information for developing new drug candidates and exercise therapy against AD.

Amyloid-β(29-42) Dimeric Conformations in Membranes Rich in Omega-3 and Omega-6 Polyunsaturated Fatty Acids

The omega-3 and omega-6 polyunsaturated fatty acids are two important components of cell membranes in human brains. When incorporated into phospholipids, omega-3 slows the progression of Alzheimer's disease (AD), whereas omega-6 is linked to increased risk of AD. Little is known on the amyloid-β (Aβ) conformations in membranes rich in omega-3 and omega-6 phospholipids. Herein, the structural properties of the Aβ_29-42 dimer embedded in both fatty acid membranes were comparatively studied to a 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine (POPC) bilayer using all-atom molecular dynamics (MD) simulations. Starting from α-helix, both omega-6 and omega-3 membranes promote new orientations and conformations of the dimer, in agreement with the observed dependence of Aβ production upon addition of these two fatty acids. This conformational result is corroborated by atomistic MD simulations of the dimer of the 99 amino acid C-terminal fragment of amyloid precursor protein spanning the residues 15-55. Starting from β-sheet, omega-6 membrane promotes helical and disordered structures of Aβ_29-42 dimer, whereas omega-3 membrane preserves the β-sheet structures differing however from those observed in POPC. Remarkably, the mixture of the two fatty acids and POPC depicts another conformational ensemble of the Aβ_29-42 dimer. This finding demonstrates that variation in the abundance of the molecular phospholipids, which changes with age, modulates membrane-embedded Aβ oligomerization.

A first-passage approach to the thermal breakage of a discrete one-dimensional chain

Using the first passage method for a Markov process, we theoretically study the fragmentation rate of a discrete one-dimensional chain (Rouse model). The fragmentation occurs due to thermal fluctuations. Assuming equilibrium initial conditions, we obtain an expression for the fragmentation rate of the one-dimensional filament as a function of the number of monomers, the position of the breaking point along the filament, the ratio of the bond energy to the thermal energy, and the Rouse relaxation time. We also obtain the fragmentation rate for a fixed initial configuration of the chain by numerically solving a Volterra equation. Our results reduce to those of previous theoretical studies at the appropriate limits, and spell out the role of the relevant time scales. The prediction of our model for the fragmentation rate of insulin fibrils under optimal growth conditions for the solution appears to be consistent with experimental data from other studies.

Protein Interaction with Charged Macromolecules: From Model Polymers to Unfolded Proteins and Post-Translational Modifications

Interaction of proteins with charged macromolecules is involved in many processes in cells. Firstly, there are many naturally occurred charged polymers such as DNA and RNA, polyphosphates, sulfated glycosaminoglycans, etc., as well as pronouncedly charged proteins such as histones or actin. Electrostatic interactions are also important for "generic" proteins, which are not generally considered as polyanions or polycations. Finally, protein behavior can be altered due to post-translational modifications such as phosphorylation, sulfation, and glycation, which change a local charge of the protein region. Herein we review molecular modeling for the investigation of such interactions, from model polyanions and polycations to unfolded proteins. We will show that electrostatic interactions are ubiquitous, and molecular dynamics simulations provide an outstanding opportunity to look inside binding and reveal the contribution of electrostatic interactions. Since a molecular dynamics simulation is only a model, we will comprehensively consider its relationship with the experimental data.

Solvent Composition Effects on the Structural Properties of the Aβ42 Monomer from the 3D-RISM-KH Molecular Theory of Solvation

Structural characterization of amyloid (A)β peptides implicated in Alzheimer's disease is a challenging problem due to their intrinsically disordered nature and their high propensity for aggregation. Only limited information is currently available from experiments on conformational properties and aggregation pathways of the peptides in cellular environments. In silico modeling complements experimental information, providing atomistic insight into structure and dynamics of different Aβ species. All-atom explicit solvent molecular dynamics (MD) simulations with a properly selected force field can deliver reliable structural and dynamic information. In the case of intrinsically disordered Aβ peptides, enhanced sampling simulations beyond the nanosecond time scale are required to obtain statistically meaningful results even for simple solvent conditions. To overcome the challenges of conformational sampling in crowded cellular environments, alternative approaches have to be used, including postprocessing of MD data. In this study, we employ the statistical-mechanical, three-dimensional reference interaction site model with the Kovalenko-Hirata closure integral equation molecular theory of solvation to describe solvent composition effects on the conformational equilibrium in a structural ensemble of the Aβ42 (covering residues 1-42) monomer based on a statistical reweighting technique. The methodology enables a computationally efficient prediction on how different factors in the cellular environment, such as solvent composition, nonpolar solvation, and macromolecular crowding, affect the structural properties of the monomer. Similarities have been identified between changes in the structural ensemble caused by nonpolar solvation and crowded environments modeled by ionic solution with large negative ions. In particular, both solvent conditions reduce the random coil content and enhance the helical structure content of the monomer. In contrast to the previous studies, which reported increased α-helical content of peptides in crowded environments, this work attributes these structural features to the difference in solvent exposure of hydrophilic residues of the monomer for different secondary structure elements, rather than to (entropic) excluded volume effects.

Frontiers in Probing Alzheimer's Disease Biomarkers with Fluorescent Small Molecules

Alzheimer's disease (AD) is the most common form of dementia. The pathogenesis of the disease is associated with aggregated amyloid-β, hyperphosphorylated tau, a high level of metal ions, abnormal enzyme activities, and reactive astrocytes. This outlook gives an overview of fluorescent small molecules targeting AD biomarkers for ex vivo and in vivo imaging. These chemical imaging probes are categorized based on the potential biomarkers, and their pros and cons are discussed. Guidelines for designing new sensing strategies as well as the desirable properties to be pursued for AD fluorescence imaging are also provided.

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.

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.

Characterizing the structural and thermodynamic properties of Aβ42 and Aβ40

The self-assembly of amyloid-beta (Aβ) proteins in aqueous extracellular environments is implicated in Alzheimer's disease. Among several alloforms of Aβ proteins differing in sequence length, the 42- and 40-residue forms (Aβ42 and Aβ40) are the most abundant ones in the human body. Although the only difference is the additional I₄₁A₄₂ residues in the C-terminus, Aβ42 exhibits more aggregation tendency and stronger neurotoxicity than Aβ40. Here, we investigate the molecular factors that confer more aggregation potential to Aβ42 than to Aβ40 based on molecular dynamics simulations combined with solvation thermodynamic analyses. It is observed that the most salient structural feature of Aβ42 relative to Aβ40 is the more enhanced β-sheet forming tendency, in particular in the C-terminal region. While such a structural characteristic of Aβ42 will certainly serve to facilitate the formation of aggregate species rich in β-sheet structure, we also detect its interesting thermodynamic consequence. Indeed, we find from the decomposition analysis that the C-terminal region substantially increases the solvation free energy (i.e., overall "hydrophobicity") of Aβ42, which is caused by the dehydration of the backbone moieties showing the enhanced tendency of forming the β-structure. Together with the two additional hydrophobic residues (I₄₁A₄₂), this leads to the higher solvation free energy of Aβ42, implying the larger water-mediated attraction toward the self-assembly. Thus, our computational results provide structural and thermodynamic grounds on why Aβ42 has more aggregation propensity than Aβ40 in aqueous environments.

Detecting lysozyme unfolding via the fluorescence of lysozyme encapsulated gold nanoclusters

Protein misfolding plays a critical role in the manifestation of amyloidosis type diseases. Therefore, understanding protein unfolding and the ability to track protein unfolding in a dynamic manner are of considerable interest. Fluorescence-based techniques are powerful tools for gaining real-time information about the local environmental conditions of a probe on the nanoscale. Fluorescent gold nanoclusters (AuNCs) are a new type of fluorescent probes which are <2 nm in diameter, incredibly robust and offer highly sensitive, wavelength tuneable emission. Their small size minimises intrusion and makes AuNCs ideal for studying protein dynamics. Lysozyme has previously been used to encapsulate AuNCs. The unfolding dynamics of lysozyme under different environmental conditions have been well-studied and being an amyloid type protein makes lysozyme an ideal candidate for encapsulating AuNCs in order to test their sensitivity to protein unfolding. In this study, we tracked the fluorescence characteristics of AuNCs encapsulated in lysozyme while inducing protein unfolding using urea, sodium dodecyl sulphate (SDS) and elevated temperature and compared them to complimentary circular dichroism spectra. It is found that AuNC fluorescence emission is quenched upon induced protein unfolding either due to a decrease in Forster Resonance Energy Transfer (FRET) efficiency between tryptophan and AuNCs or solvent exposure of the AuNC. Fluorescence lifetime measurements confirmed quenching to be collisional via oxygen dissolved in a solution which increases as the AuNC was exposed to the solvent during unfolding. Moreover, the longer decay component τ₁ was observed to decrease as the protein unfolded, due to the increased collisional quenching. It is suggested that AuNC sensitivity to solvent exposure might be utilised in the future as a new approach to studying and possibly even detecting amyloidosis type diseases.

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.

Molecular dynamics simulation on the effect of transition metal binding to the N-terminal fragment of amyloid-β

We report molecular dynamics simulations of three possible adducts of Fe(II) to the N-terminal 1-16 fragments of the amyloid-β peptide, along with analogous simulations of Cu(II) and Zn(II) adducts. We find that multiple simulations from different starting points reach pseudo-equilibration within 100-300 ns, leading to over 900 ns of equilibrated trajectory data for each system. The specifics of the coordination modes for Fe(II) have only a weak effect on peptide secondary and tertiary structures, and we therefore compare one of these with analogous models of Cu(II) and Zn(II) complexes. All share broadly similar structural features, with mixture of coil, turn and bend in the N-terminal region and helical structure for residues 11-16. Within this overall pattern, subtle effects due to changes in metal are evident: Fe(II) complexes are more compact and are more likely to occupy bridge and ribbon regions of Ramachandran maps, while Cu(II) coordination leads to greater occupancy of the poly-proline region. Analysis of representative clusters in terms of molecular mechanics energy and atoms-in-molecules properties indicates similarity of four-coordinate Cu and Zn complexes, compared to five-coordinate Fe complex that exhibits lower stability and weaker metal-ligand bonding. Communicated by Ramaswamy H. Sarma.

Selection of Membrane RNA Aptamers to Amyloid Beta Peptide: Implications for Exosome-Based Antioxidant Strategies

The distribution of amyloid beta peptide 42 (Aβ42) between model exosomal membranes and a buffer solution was measured. The model membranes contained liquid-ordered regions or phosphatidylserine. Results demonstrated that up to ca. 20% of amyloid peptide, generated in the plasma (or intracellular) membrane as a result of proteolytic cleavage of amyloid precursor proteins by β- and γ-secretases, can stay within the membrane milieu. The selection of RNA aptamers that bind to Aβ42 incorporated into phosphatidylserine-containing liposomal membranes was performed using the selection-amplification (SELEX) method. After eight selection cycles, the pool of RNA aptamers was isolated and its binding to Aβ42-containing membranes was demonstrated using the gel filtration method. Since membranes can act as a catalytic surface for Aβ42 aggregation, these RNA aptamers may inhibit the formation of toxic amyloid aggregates that can permeabilize cellular membranes or disrupt membrane receptors. Strategies are proposed for using functional exosomes, loaded with RNA aptamers specific to membrane Aβ42, to reduce the oxidative stress in Alzheimer's disease and Down's syndrome.

Alkali ion influence on structure and stability of fibrillar amyloid-β oligomers

Alzheimer's disease is characterized by the aggregation of Amyloid-β (Aβ) peptide into oligomers, fibrils and plaques. Many factors influencing this process as well as the stability of the various Aβ aggregates are known to date, and include the concentration and type of metal ions. Most experimental and theoretical studies have concentrated on heavy metal ions, like Fe²⁺, Zn²⁺, or Cu²⁺, while the smaller alkali ions Li⁺, Na⁺, and K⁺ have not gained much attention notwithstanding their role and ubiquity in physiological environments. In this work, we applied atomistic molecular dynamics simulations to investigate the potential role of these alkali ions in stabilizing fibrillar Aβ oligomers of different size and topology, i.e., single and double filament systems comprising 3-24 peptide chains per filament. We find a pronounced difference on the molecular level in the interaction behavior with free carboxylate groups of the Aβ oligomer: Li⁺ forms stable bridged interactions, whereas K⁺ interacts more transiently and lacks bridging. The behavior of Na⁺ is in between, so that this ion-protein interaction obeys the renowned Hofmeister series. These differences are also reflected in the ability of the alkali ions to stabilize the oligomer secondary structure. The stabilizing effect is most pronounced for the smaller fibrillar oligomers, suggesting that the type of alkali ion critically affects the initial stages of fibril formation. Our findings thus offer a molecular explanation for the observation that the polymorphisms of Aβ fibril structures are caused by differences in the surrounding ionic environment. Graphical abstract Influence of alkali ions on the structure and stability of fibrillar amyloid-β oligomers.