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

The molecular behavior of a single β-amyloid inside a dipalmitoylphosphatidylcholine bilayer at three different temperatures: An atomistic simulation study: Aβ interaction with DPPC: Atomistic simulation

The behavior of a single Aβ40 molecule within a dipalmitoylphosphatidylcholine (DPPC) bilayer was studied by all-atom molecular dynamics simulations. The effect of membrane structure was investigated on Aβ40 behavior, secondary structure, and insertion depth. Simulations were performed at three temperatures (323, 310, and 300 K) to probe three different bilayer fluidities. Results show that at all above temperatures, the peptide contains two short helices, coil, bend, and turn structures. At 300 K, the peptide contains a region with β structure in C-terminal region. Our results also show that Aβ decreases the bilayer thickness and the order of lipids in its vicinity which leads to water insertion into the bilayer and concomitant increase in the local fluidity. The peptide remains embedded in the bilayer at all temperatures, and become inserted into the bilayer up to several residues at 323 and 310 K. At 310 and 300 K, the dominant interaction energy between Aβ and bilayer changes from electrostatic to van der Waals. It can be proposed that at higher temperatures (e.g., 323 K), Lys28 and the C-terminal region of the peptide play the role of two anchors that keep Aβ inside the top leaflet. This study demonstrates that Aβ molecule can perturb the integrity of cellular membranes. Proteins 2017; 85:1298-1310. © 2017 Wiley Periodicals, Inc.

Rational design, synthesis and biological screening of triazine-triazolopyrimidine hybrids as multitarget anti-Alzheimer agents

In our endeavor towards the development of potent multitarget ligands for the treatment of Alzheimer's disease, a series of triazine-triazolopyrimidine hybrids were designed, synthesized and characterized by various spectral techniques. Docking and scoring techniques were used to design the inhibitors and to display their interaction with key residues of active site. Organic synthesis relied upon convergent synthetic routes were mono and di-substituted triazines were connected with triazolopyrimidine using piperazine as a linker. In total, seventeen compounds were synthesized in which the di-substituted triazine-triazolopyrimidine derivatives 9a-d showed better acetylcholinesterase (AChE) inhibitory activity than the corresponding tri-substituted triazine-triazolopyrimidine derivatives 10a-f. Out of the disubstituted triazine-triazolopyrimidine based compounds, 9a and 9b showed encouraging inhibitory activity on AChE with IC₅₀ values 0.065 and 0.092 μM, respectively. Interestingly, 9a and 9b also demonstrated good inhibition selectivity towards AChE over BuChE by ∼28 folds. Furthermore, kinetic analysis and molecular modeling studies showed that 9a and 9b target both catalytic active site as well as peripheral anionic site of AChE. In addition, these derivatives effectively modulated Aβ self-aggregation as investigated through CD spectroscopy, ThT fluorescence assay and electron microscopy. Besides, these compounds exhibited potential antioxidants (2.15 and 2.91 trolox equivalent by ORAC assay) and metal chelating properties. In silico ADMET profiling highlighted that, these novel triazine derivatives have appropriate drug like properties and possess very low toxic effects in the primarily pharmacokinetic study. Overall, the multitarget profile exerted by these novel triazine molecules qualified them as potential anti-Alzheimer drug candidates in AD therapy.

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

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

Conformational Ensembles of the Wild-Type and S8C Aβ1-42 Dimers

We characterized the dimer of the amyloid-β wild-type (WT) peptide, Aβ, of 42 residues and its disulfide-bond-locked double mutant (S8C) by replica exchange molecular dynamics simulations. Aβ dimers are known to be the smallest toxic species in Alzheimer's disease, and the S8C mutant has been shown experimentally to form an exclusive homogeneous and neurotoxic dimer. Our 50 μs all-atom simulations reveal similar secondary structures and collision cross-sections but very different intramolecular and intermolecular conformations upon double S8C mutation. Both dimers are very dynamic with hundreds of free-energy minima that differ from the U-shape and S-shape conformations of the peptides in the fibrils. The only common structural feature, shared by both species with a probability of 4% in WT and 12% in S8C-S8C, is a three-stranded β-sheet spanning the 17-23, 29-36, and 39-41 residues, which does not exist in the Aβ40 WT dimers.

Structural properties of amyloid β(1-40) dimer explored by replica exchange molecular dynamics simulations

Replica exchange molecular dynamics simulations (300 ns) were used to study the dimerization of amyloid β(1-40) (Aβ(1-40)) polypeptide. Configurational entropy calculations revealed that at physiological temperature (310 K, 37°C) dynamic dimers are formed by randomly docked monomers. Free energy of binding of the two chains to each other was -93.56 ± 6.341 kJ mol^-1 . Prevalence of random coil conformations was found for both chains with the exceptions of increased β-sheet content from residues 16-21 and 29-32 of chain A and residues 15-21 and 30-33 of chain B with β-turn/β-bend conformations in both chains from residues 1-16, 21-29 of chain A, 1-16, and 21-29 of chain B. There is a mixed β-turn/β-sheet region from residues 33-38 of both chains. Analysis of intra- and interchain residue distances shows that, although the individual chains are highly flexible, the dimer system stays in a loosely packed antiparallel β-sheet configuration with contacts between residues 17-21 of chain A with residues 17-21 and 31-36 of chain B as well as residues 31-36 of chain A with residues 17-21 and 31-36 of chain B. Based on dihedral principal component analysis, the antiparallel β-sheet-loop-β-sheet conformational motif is favored for many low energy sampled conformations. Our results show that Aβ(1-40) can form dynamic dimers in aqueous solution that have significant conformational flexibility and are stabilized by collapse of the central and C-terminal hydrophobic cores with the expected β-sheet-loop-β-sheet conformational motif. Proteins 2017; 85:1024-1045. © 2017 Wiley Periodicals, Inc.

Alzheimer's Protective Cross-Interaction between Wild-Type and A2T Variants Alters Aβ42 Dimer Structure

Whole genome sequencing has recently revealed the protective effect of a single A2T mutation in heterozygous carriers against Alzheimer's disease (AD) and age-related cognitive decline. The impact of the protective cross-interaction between the wild-type (WT) and A2T variants on the dimer structure is therefore of high interest, as the Aβ dimers are the smallest known neurotoxic species. Toward this goal, extensive atomistic replica exchange molecular dynamics simulations of the solvated WT homo- and A2T hetero- Aβ_1-42 dimers have been performed, resulting into a total of 51 μs of sampling for each system. Weakening of a set of transient, intrachain contacts formed between the central and C-terminal hydrophobic residues is observed in the heterodimeric system. The majority of the heterodimers with reduced interaction between central and C-terminal regions lack any significant secondary structure and display a weak interchain interface. Interestingly, the A2T N-terminus, particularly residue F4, is frequently engaged in tertiary and quaternary interactions with central and C-terminal hydrophobic residues in those distinct structures, leading to hydrophobic burial. This atypical involvement of the N-terminus within A2T heterodimer revealed in our simulations implies possible interference on Aβ₄₂ aggregation and toxic oligomer formation, which is consistent with experiments. In conclusion, the present study provides detailed structural insights onto A2T Aβ₄₂ heterodimer, which might provide molecular insights onto the AD protective effect of the A2T mutation in the heterozygous state.

Stereoisomers Probe Steric Zippers in Amyloid-β

Shape complementarity between close-packed residues plays a critical role in the amyloid aggregation process. Here, we probe such "steric zipper" interactions in amyloid-β (Aβ₄₀), whose aggregation is linked to Alzheimer's disease, by replacing natural residues by their stereoisomers. Such mutations are expected to specifically destabilize the shape sensitive "packing" interactions, which may potentially increase their solubility and change other properties. We study the stereomutants DF19 and DL34 and also the DA2/DF4/DH6/DS8 mutant of Aβ₄₀. F19-L34 is a critical contact in a tightly packed region of Aβ, while residues 1-9 are known to be disordered. While both DF19 and DL34 slow down the kinetics of aggregation and form amyloid fibrils efficiently, only DL34 increases the final solubility. DF19 gives rise to additional off-pathway aggregation which results in large, kinetically stable aggregates, and has lower net solubility. DA2/DF4/DH6/DS8 does not have an effect on the kinetics or the solubility. Notably, both DF19 and DL34 oligomers have a significantly lower level of interactions with lipid vesicles and live cells. We conclude that stereoisomers can cause complex site dependent changes in amyloid properties, and provide an effective tool to determine the role of individual residues in shaping the packed interiors of amyloid aggregates.

Fibrillation-prone conformations of the amyloid-β-42 peptide at the gold/water interface

Proteins in the proximity of inorganic surfaces and nanoparticles may undergo profound adjustments that trigger biomedically relevant processes, such as protein fibrillation. The mechanisms that govern protein-surface interactions at the molecular level are still poorly understood. In this work, we investigate the adsorption onto a gold surface, in water, of an amyloid-β (Aβ) peptide, which is the amyloidogenic peptide involved in Alzheimer's disease. The entire adsorption process, from the peptide in bulk water to its conformational relaxation on the surface, is explored by large-scale atomistic molecular dynamics (MD) simulations. We start by providing a description of the conformational ensemble of Aβ in solution by a 22 μs temperature replica exchange MD simulation, which is consistent with previous results. Then, we obtain a statistical description of how the peptide approaches the gold surface by multiple MD simulations, identifying the preferential gold-binding sites and giving a kinetic picture of the association process. Finally, relaxation of the Aβ conformations at the gold/water interface is performed by a 19 μs Hamiltonian-temperature replica exchange MD simulation. We find that the conformational ensemble of Aβ is strongly perturbed by the presence of the surface. In particular, at the gold/water interface the population of the conformers akin to amyloid fibrils is significantly enriched, suggesting that this extended contact geometry may promote fibrillation.

Computer Simulations of Intrinsically Disordered Proteins

The investigation of intrinsically disordered proteins (IDPs) is a new frontier in structural and molecular biology that requires a new paradigm to connect structural disorder to function. Molecular dynamics simulations and statistical thermodynamics potentially offer ideal tools for atomic-level characterizations and thermodynamic descriptions of this fascinating class of proteins that will complement experimental studies. However, IDPs display sensitivity to inaccuracies in the underlying molecular mechanics force fields. Thus, achieving an accurate structural characterization of IDPs via simulations is a challenge. It is also daunting to perform a configuration-space integration over heterogeneous structural ensembles sampled by IDPs to extract, in particular, protein configurational entropy. In this review, we summarize recent efforts devoted to the development of force fields and the critical evaluations of their performance when applied to IDPs. We also survey recent advances in computational methods for protein configurational entropy that aim to provide a thermodynamic link between structural disorder and protein activity.

Hydrodynamic effects on β-amyloid (16-22) peptide aggregation

Computer simulations based on simplified representations are routinely used to explore the early steps of amyloid aggregation. However, when protein models with implicit solvent are employed, these simulations miss the effect of solvent induced correlations on the aggregation kinetics and lifetimes of metastable states. In this work, we apply the multi-scale Lattice Boltzmann Molecular Dynamics technique (LBMD) to investigate the initial aggregation phases of the amyloid Aβ16-22 peptide. LBMD includes naturally hydrodynamic interactions (HIs) via a kinetic on-lattice representation of the fluid kinetics. The peptides are represented by the flexible OPEP coarse-grained force field. First, we have tuned the essential parameters that control the coupling between the molecular and fluid evolutions in order to reproduce the experimental diffusivity of elementary species. The method is then deployed to investigate the effect of HIs on the aggregation of 100 and 1000 Aβ16-22 peptides. We show that HIs clearly impact the aggregation process and the fluctuations of the oligomer sizes by favouring the fusion and exchange dynamics of oligomers between aggregates. HIs also guide the growth of the leading largest cluster. For the 100 Aβ16-22 peptide system, the simulation of ∼300 ns allowed us to observe the transition from ellipsoidal assemblies to an elongated and slightly twisted aggregate involving almost the totality of the peptides. For the 1000 Aβ16-22 peptides, a system of unprecedented size at quasi-atomistic resolution, we were able to explore a branched disordered fibril-like structure that has never been described by other computer simulations, but has been observed experimentally.

A computational study on how structure influences the optical properties in model crystal structures of amyloid fibrils

Amyloid fibrils have been shown to have peculiar optical properties since they can exhibit fluorescence in the absence of aromatic residues. In a recent study, we have shown that proton transfer (PT) events along hydrogen bonds (HBs) are coupled to absorption in the near UV range. Here, we gain more insights into the different types of hydrogen bonding interactions that occur in our model systems and the molecular factors that control the susceptibility of the protons to undergo PT and how this couples to the optical properties. In the case of the strong N-C termini interactions, a nearby methionine residue stabilizes the non-zwitterionic NH₂-COOH pair, while zwitterionic NH₃⁺-COO- is stabilized by the proximity of nearby crystallographic water molecules. Proton motion along the hydrogen bonds in the fibril is intimately coupled to the compression of the heavier atoms, similar to what is observed in bulk water. Small changes in the compression of the hydrogen bonds in the protein can lead to significant changes in both the ground and excited state potential energy surfaces associated with PT. Finally, we also reinforce the importance of nuclear quantum fluctuations of protons in the HBs of the amyloid proteins.

Comparison of force fields for Alzheimer's A β42: A case study for intrinsically disordered proteins

Intrinsically disordered proteins are essential for biological processes such as cell signalling, but are also associated to devastating diseases including Alzheimer's disease, Parkinson's disease or type II diabetes. Because of their lack of a stable three-dimensional structure, molecular dynamics simulations are often used to obtain atomistic details that cannot be observed experimentally. The applicability of molecular dynamics simulations depends on the accuracy of the force field chosen to represent the underlying free energy surface of the system. Here, we use replica exchange molecular dynamics simulations to test five modern force fields, OPLS, AMBER99SB, AMBER99SB*ILDN, AMBER99SBILDN-NMR and CHARMM22*, in their ability to model Aβ42 , an intrinsically disordered peptide associated with Alzheimer's disease, and compare our results to nuclear magnetic resonance (NMR) experimental data. We observe that all force fields except AMBER99SBILDN-NMR successfully reproduce local NMR observables, with CHARMM22* being slightly better than the other force fields.

Numerical Simulations Reveal Randomness of Cu(II) Induced Aβ Peptide Dimerization under Conditions Present in Glutamatergic Synapses

The interactions between the Aβ1-40 molecules species and the copper ions (Cu(II)) were intensively investigated due to their potential role in the development of the Alzheimer Disease (AD). The rate and the mechanism of the Cu(II)-Aβ complexes formation determines the aggregation pathway of the Aβ species, starting from smaller but more cytotoxic oligomers and ending up in large Aβ plaques, being the main hallmark of the AD. In our study we exploit the existing knowledge on the Cu(II)-Aβ interactions and create the theoretical model of the initial phase of the copper- driven Aβ aggregation mechanism. The model is based on the direct solution of the Chemical Master Equations, which capture the inherent stochastics of the considered system. In our work we argue that due to a strong Cu(II) affinity to Aβ and temporal accessibility of the Cu(II) ions during normal synaptic activity the aggregation driven by Cu(II) dominates the pure Aβ aggregation. We also demonstrate the dependence of the formation of different Cu(II)-Aβ complexes on the concentrations of reagents and the synaptic activity. Our findings correspond to recent experimental results and give a sound hypothesis on the AD development mechanisms.

Oligomerization of FVFLM peptides and their ability to inhibit beta amyloid peptides aggregation: consideration as a possible model

Preeclampsia, a pregnancy-specific disorder, shares typical pathophysiological features with protein misfolding disorders including Alzheimer's disease. Characteristic for preeclampsia is the involvement of multiple proteins of which fragments of SERPINA1 and β-amyloid co-aggregate in urine and placenta of preeclamptic women. To explore the biophysical basis of this interaction, we investigated the multidimensional efficacy of the FVFLM sequence in SERPINA1, as a model inhibitory agent of β-amyloid aggregation. After studying the oligomerization of FVFLM peptides using all-atom molecular dynamics simulations with the GROMOS43a1 force field and explicit water, we report that FVFLM can aggregate and its aggregation is spontaneous with a remarkably faster rate than that recorded for KLVFF (aggregation "hot-spot" from β-amyloid). The fast kinetics of FVFLM aggregation was found to be driven primarily by core-like aromatic interactions originating from the anti-parallel orientation of complementarily uncharged strands. The conspicuously stable aggregation mechanism observed for FVFLM peptides is found not to conform to the popular 'dock-lock' scheme. We also found high propensity of FVFLM for KLVFF binding. When present, FVFLM disrupts the β-amyloid aggregation pathway and we propose that FVFLM-like peptides might be used to prevent the assembly of full-length Aβ or other pro-amyloidogenic peptides into amyloid fibrils.

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

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

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

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

Replica exchange molecular dynamics study of the amyloid beta (11–40) trimer penetrating a membrane

The transmembrane Aβ_11–40 trimer is investigated for the first time using REMD and FEP.

Curcumin-like compounds designed to modify amyloid beta peptide aggregation patterns

This study suggests new concepts and potential difficulties in the design of novel drugs against diverse amyloidoses, including Alzheimer’s disease.

Molecular modeling studies of quinazolinone derivatives as novel PI3Kδ selective inhibitors

The main molecular modeling method, the docking results of newly designed compoundD04and the best pharmacophore model are reported herein.

The GOR Method of Protein Secondary Structure Prediction and Its Application as a Protein Aggregation Prediction Tool

The GOR method of protein secondary structure prediction is described. The original method was published by Garnier, Osguthorpe, and Robson in 1978 and was one of the first successful methods to predict protein secondary structure from amino acid sequence. The method is based on information theory, and an assumption that information function of a protein chain can be approximated by a sum of information from single residues and pairs of residues. The analysis of frequencies of occurrence of secondary structure for singlets and doublets of residues in a protein database enables prediction of secondary structure for new amino acid sequences. Because of these simple physical assumptions the GOR method has a conceptual advantage over other later developed methods such as PHD, PSIPRED, and others that are based on Machine Learning methods (like Neural Networks), give slightly better predictions, but have a "black box" nature. The GOR method has been continuously improved and modified for 30 years with the last GOR V version published in 2002, and the GOR V server developed in 2005. We discuss here the original GOR method and the GOR V program and the web server. Additionally we discuss new highly interesting and important applications of the GOR method to chameleon sequences in protein folding simulations, and for prediction of protein aggregation propensities. Our preliminary studies show that the GOR method is a promising and efficient alternative to other protein aggregation predicting tools. This shows that the GOR method despite being almost 40 years old is still important and has significant potential in application to new scientific problems.

Calycosin improves cognitive function in a transgenic mouse model of Alzheimer's disease by activating the protein kinase C pathway

The major pathological changes in Alzheimer's disease are beta amyloid deposits and cognitive impairment. Calycosin is a typical phytoestrogen derived from radix astragali that binds to estrogen receptors to produce estrogen-like effects. Radix astragali Calycosin has been shown to relieve cognitive impairment induced by diabetes mellitus, suggesting calycosin may improve the cognitive function of Alzheimer's disease patients. The protein kinase C pathway is upstream of the mitogen-activated protein kinase pathway and exerts a neuroprotective effect by regulating Alzheimer's disease-related beta amyloid degradation. We hypothesized that calycosin improves the cognitive function of a transgenic mouse model of Alzheimer's disease by activating the protein kinase C pathway. Various doses of calycosin (10, 20 and 40 mg/kg) were intraperitoneally injected into APP/PS1 transgenic mice that model Alzheimer's disease. Calycosin diminished hippocampal beta amyloid, Tau protein, interleukin-1beta, tumor necrosis factor-alpha, acetylcholinesterase and malondialdehyde levels in a dose-dependent manner, and increased acetylcholine and glutathione activities. The administration of a protein kinase C inhibitor, calphostin C, abolished the neuroprotective effects of calycosin including improving cognitive ability, and anti-oxidative and anti-inflammatory effects. Our data demonstrated that calycosin mitigated oxidative stress and inflammatory responses in the hippocampus of Alzheimer's disease model mice by activating the protein kinase C pathway, and thereby improving cognitive function.

Replica exchange molecular dynamics study of the truncated amyloid beta (11–40) trimer in solution

The structure of the 3Aβ_11–40 oligomer is determined for the first time using T-REMD simulations.

Structural and thermodynamical properties of early human amylin oligomers using replica exchange molecular dynamics: mutation effect of three key residues F15, H18 and F23

A schematic representation of a possible oligomerization mechanism of hIAPP. β-Hairpins are proposed to self-assemble into early ordered oligomers by side-to-side association.

AChE Inhibition-based Multi-target-directed Ligands, a Novel Pharmacological Approach for the Symptomatic and Disease-modifying Therapy of Alzheimer's Disease

Alzheimer's disease (AD) is the most common form of dementia in elder people, characterised by a progressive decline in memory as a result of an impairment of cholinergic neurotransmission. To date acetylcholinesterase inhibitors (AChEIs) have become the most prescribed drugs for the symptomatic treatment of mild to moderate AD. However, the traditional “one molecule-one target” paradigm is not sufficient and appropriate to yield the desired therapeutic efficacy since multiple factors, such as amyloid-β (Aβ) deposits, neuroinflammation, oxidative stress, and decreased levels of acetylcholine (ACh) have been thought to play significant roles in the AD pathogenesis. New generation of multi-target drugs is earnestly demanded not only for ameliorating symptoms but also for modifying the disease. Herein, we delineated the catalytic and non-catalytic functions of AChE, and summarized the works of our group and others in research and development of novel AChEI-based multi-target-directed ligands (MTDLs), such as dual binding site AChEIs and multi-target AChEIs inhibiting Aβ aggregation, regulating Aβ procession, antagonizing platelet-activating factor (PAF) receptor, scavenging oxygen radical, chelating metal ions, inhibiting monoamine oxidase B (MAO-B), blocking N-methyl-D-aspartic acid (NMDA) receptor and others.

Characterization of Mn(II) ion binding to the amyloid-β peptide in Alzheimer's disease

Growing evidence links neurodegenerative diseases to metal exposure. Aberrant metal ion concentrations have been noted in Alzheimer's disease (AD) brains, yet the role of metals in AD pathogenesis remains unresolved. A major factor in AD pathogenesis is considered to be aggregation of and amyloid formation by amyloid-β (Aβ) peptides. Previous studies have shown that Aβ displays specific binding to Cu(II) and Zn(II) ions, and such binding has been shown to modulate Aβ aggregation. Here, we use nuclear magnetic resonance (NMR) spectroscopy to show that Mn(II) ions also bind to the N-terminal part of the Aβ(1-40) peptide, with a weak binding affinity in the milli- to micromolar range. Circular dichroism (CD) spectroscopy, solid state atomic force microscopy (AFM), fluorescence spectroscopy, and molecular modeling suggest that the weak binding of Mn(II) to Aβ may not have a large effect on the peptide's aggregation into amyloid fibrils. However, identification of an additional metal ion displaying Aβ binding reveals more complex AD metal chemistry than has been previously considered in the literature.

Polymorphism of fibrillar structures depending on the size of assembled Aβ17-42 peptides

The size of assembled Aβ_17-42 peptides can determine polymorphism during oligomerization and fibrillization, but the mechanism of this effect is unknown. Starting from separate random monomers, various fibrillar oligomers with distinct structural characteristics were identified using discontinuous molecular dynamics simulations based on a coarse-grained protein model. From the structures observed in the simulations, two characteristic oligomer sizes emerged, trimer and paranuclei, which generated distinct structural patterns during fibrillization. A majority of the simulations for trimers and tetramers formed non-fibrillar oligomers, which primarily progress to off-pathway oligomers. Pentamers and hexamers were significantly converted into U-shape fibrillar structures, meaning that these oligomers, called paranuclei, might be potent on-pathway intermediates in fibril formation. Fibrillar oligomers larger than hexamers generated substantial polymorphism in which hybrid structures were readily formed and homogeneous fibrillar structures appeared infrequently.

Dimerization Mechanism of Alzheimer Aβ40 Peptides: The High Content of Intrapeptide-Stabilized Conformations in A2V and A2T Heterozygous Dimers Retards Amyloid Fibril Formation

Amyloid beta (Aβ) oligomerization is associated with the origin and progression of Alzheimer's disease (AD). While the A2V mutation enhances aggregation kinetics and toxicity, mixtures of wild-type (WT) and A2V, and also WT and A2T, peptides retard fibril formation and protect against AD. In this study, we simulate the equilibrium ensemble of WT:A2T Aβ₄₀ dimer by means of extensive atomistic replica exchange molecular dynamics and compare our results with previous equivalent simulations of A2V:A2V, WT:WT, and WT:A2V Aβ₄₀ dimers for a total time scale of nearly 0.1 ms. Qualitative comparison of the resulting thermodynamic properties, such as the relative binding free energies, with the reported experimental kinetic and thermodynamic data affords us important insight into the conversion from slow-pathway to fast-pathway dimer conformations. The crucial reaction coordinate or driving force of such transformation turns out to be related to hydrophobic interpeptide interactions. Analysis of the equilibrium ensembles shows that the fast-pathway conformations contain interpeptide out-of-register antiparallel β-sheet structures at short interpeptide distances. In contrast, the slow-pathway conformations are formed by the association of peptides at large interpeptide distances and high intrapeptide compactness, such as conformations containing intramolecular three-stranded β-sheets which sharply distinguish fast (A2V:A2V and WT:WT) and slow (WT:A2T and WT:A2V) amyloid-forming sequences. Also, this analysis leads us to predict that a molecule stabilizing the intramolecular three-stranded β-sheet or inhibiting the formation of an interpeptide β-sheet spanning residues 17-20 and 31-37 would further reduce fibril formation and probably the cytotoxicity of Aβ species.

Multiscale simulation of molecular processes in cellular environments

We describe the recent advances in studying biological systems via multiscale simulations. Our scheme is based on a coarse-grained representation of the macromolecules and a mesoscopic description of the solvent. The dual technique handles particles, the aqueous solvent and their mutual exchange of forces resulting in a stable and accurate methodology allowing biosystems of unprecedented size to be simulated.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.

Nonequilibrium all-atom molecular dynamics simulation of the bubble cavitation and application to dissociate amyloid fibrils

The cavitation of gas bubbles in liquids has been applied to different disciplines in life and natural sciences, and in technologies. To obtain an appropriate theoretical description of effects induced by the bubble cavitation, we develop an all-atom nonequilibrium molecular-dynamics simulation method to simulate bubbles undergoing harmonic oscillation in size. This allows us to understand the mechanism of the bubble cavitation-induced liquid shear stress on surrounding objects. The method is then employed to simulate an Aβ fibril model in the presence of bubbles, and the results show that the bubble expansion and contraction exert water pressure on the fibril. This yields to the deceleration and acceleration of the fibril kinetic energy, facilitating the conformational transition between local free energy minima, and leading to the dissociation of the fibril. Our work, which is a proof-of-concept, may open a new, efficient way to dissociate amyloid fibrils using the bubble cavitation technique, and new venues to investigate the complex phenomena associated with amyloidogenesis.

Recent progress in repositioning Alzheimer's disease drugs based on a multitarget strategy

Alzheimer's disease (AD) is a serious progressive neurological disorder, characterized by impaired cognition and profound irreversible memory loss. The multifactorial nature of AD and the absence of a cure so far have stimulated medicinal chemists worldwide to follow multitarget drug-design strategies based on repositioning approved drugs. This review describes a summary of recently published works focused on tailoring new derivatives of US FDA-approved acetylcholinesterase inhibitors, in addition to huperzine (a drug approved in China), either by hybridization with other pharmacophore elements (to hit more AD targets), or by combination of two FDA-approved drugs. Besides the capacity for improving the cholinergic activity, these polyfunctional derivatives are also able to tackle other important neuroprotective properties, such as anti-β-amyloid aggregation, scavenging of radical oxygen species, modulation of redox-active metals or inhibition of monoamine oxidase, thereby resulting in potentially novel and more effective therapeutics for the treatment of AD.

Protein Aggregation and Molecular Crowding: Perspectives From Multiscale Simulations

Cells are extremely crowded environments, thus the use of diluted salted aqueous solutions containing a single protein is too simplistic to mimic the real situation. Macromolecular crowding might affect protein structure, folding, shape, conformational stability, binding of small molecules, enzymatic activity, interactions with cognate biomolecules, and pathological aggregation. The latter phenomenon typically leads to the formation of amyloid fibrils that are linked to several lethal neurodegenerative diseases, but that can also play a functional role in certain organisms. The majority of molecular simulations performed before the last few years were conducted in diluted solutions and were restricted both in the timescales and in the system dimensions by the available computational resources. In recent years, several computational solutions were developed to get close to physiological conditions. In this review we summarize the main computational techniques used to tackle the issue of protein aggregation both in a diluted and in a crowded environment.

Impact of Cu(II) Binding on Structures and Dynamics of Aβ42 Monomer and Dimer: Molecular Dynamics Study

The classical force field, which is compatible with the Amber force field 99SB, has been obtained for the interaction of Cu(II) with monomer and dimers of amyloid-β peptides using the coordination where Cu(II) is bound to His6, His13 (or His14), and Asp1 with distorted planar geometry. The newly developed force field and molecular dynamics simulation were employed to study the impact of Cu(II) binding on structures and dynamics of Aβ₄₂ monomer and dimers. It was shown that in the presence of Cu(II) the β content of monomer is reduced substantially compared with the wild-type Aβ₄₂ suggesting that, in accord with experiments, metal ions facilitate formation of amorphous aggregates rather than amyloid fibrils with cross-β structures. In addition, one possible mechanism for amorphous assembly is that the Asp23-Lys28 salt bridge, which plays a crucial role in β sheet formation, becomes more flexible upon copper ion binding to the Aβ N-terminus. The simulation of dimers was conducted with the Cu(II)/Aβ stoichiometric ratios of 1:1 and 1:2. For the 1:1 ratio Cu(II) delays the Aβ dimerization process as observed in a number of experiments. The mechanism underlying this phenomenon is associated with slow formation of interchain salt bridges in dimer as well as with decreased hydrophobicity of monomer upon Cu-binding.

Amyloid-β 1-24 C-terminal truncated fragment promotes amyloid-β 1-42 aggregate formation in the healthy brain

Substantial data indicate that amyloid-β (Aβ), the major component of senile plaques, plays a central role in Alzheimer’s Disease and indeed the assembly of naturally occurring amyloid peptides into cytotoxic aggregates is linked to the disease pathogenesis. Although Aβ42 is a highly aggregating form of Aβ, the co-occurrence of shorter Aβ peptides might affect the aggregation potential of the Aβ pool. In this study we aimed to assess whether the structural behavior of human Aβ42 peptide inside the brain is influenced by the concomitant presence of N-terminal fragments produced by the proteolytic activity of glial cells. We show that the occurrence of the human C-terminal truncated 1–24 Aβ fragment impairs Aβ42 clearance through blood brain barrier and promotes the formation of Aβ42 aggregates even in the healthy brain. By showing that Aβ1-24 has seeding properties for aggregate formation in intracranially injected wild type mice, our study provide the proof-of-concept that peptides produced upon Aβ42 cleavage by activated glial cells may cause phenotypic defects even in the absence of genetic mutations associated with Alzheimer’s Disease, possibly contributing to the development of the sporadic form of the pathology.

Copper Coordination Features of Human Islet Amyloid Polypeptide: The Type 2 Diabetes Peptide

Human islet amyloid polypeptide (hIAPP) is the major component of amyloid deposits found in pancreatic β-cells of patients with type 2 diabetes (T2D). Copper ions have an inhibitory effect on the amyloid aggregation of hIAPP, and they may play a role in the etiology of T2D. However, deeper knowledge of the structural details of the copper-hIAPP interaction is required to understand the molecular mechanisms involved. Here, we performed a spectroscopic study of Cu(II) binding to hIAPP and several variants, using electron paramagnetic resonance (EPR), nuclear magnetic resonance (NMR), electronic absorption, and circular dichroism (CD) in the UV-vis region in combination with Born-Oppenheimer molecular dynamics (BOMD) and density functional theory geometry optimizations. We find that Cu(II) binds to the imidazole N1 of His18, the deprotonated amides of Ser19 and Ser20, and an oxygen-based ligand provided by Ser20, either via its hydroxyl group or its backbone carbonyl, while Asn22 might also play a role as an axial ligand. Ser20 plays a crucial role in stabilizing Cu(II) coordination toward the C-terminal, providing a potential link between the S20G mutation associated with early onset of T2D, its impact in Cu binding properties, and hIAPP amyloid aggregation. Our study defines the nature of the coordination environment in the Cu(II)-hIAPP complex, revealing that the amino acid residues involved in metal ion binding are also key residues for the formation of β-sheet structures and amyloid fibrils. Cu(II) binding to hIAPP may lead to the coexistence of more than one coordination mode, which in turn could favor different sets of Cu-induced conformational ensembles. Cu-induced hIAPP conformers would display a higher energetic barrier to form amyloid fibrils, hence explaining the inhibitory effect of Cu ions in hIAPP aggregation. Overall, this study provides further structural insights into the bioinorganic chemistry of T2D.

Exploring the aggregation free energy landscape of the amyloid-β protein (1-40)

A predictive coarse-grained protein force field [associative memory, water-mediated, structure, and energy model for molecular dynamics (AWSEM)-MD] is used to study the energy landscapes and relative stabilities of amyloid-β protein (1-40) in the monomer and all of its oligomeric forms up to an octamer. We find that an isolated monomer is mainly disordered with a short α-helix formed at the central hydrophobic core region (L17-D23). A less stable hairpin structure, however, becomes increasingly more stable in oligomers, where hydrogen bonds can form between neighboring monomers. We explore the structure and stability of both prefibrillar oligomers that consist of mainly antiparallel β-sheets and fibrillar oligomers with only parallel β-sheets. Prefibrillar oligomers are polymorphic but typically take on a cylindrin-like shape composed of mostly antiparallel β-strands. At the concentration of the simulation, the aggregation free energy landscape is nearly downhill. We use umbrella sampling along a structural progress coordinate for interconversion between prefibrillar and fibrillar forms to identify a conversion pathway between these forms. The fibrillar oligomer only becomes favored over its prefibrillar counterpart in the pentamer where an interconversion bottleneck appears. The structural characterization of the pathway along with statistical mechanical perturbation theory allow us to evaluate the effects of concentration on the free energy landscape of aggregation as well as the effects of the Dutch and Arctic mutations associated with early onset of Alzheimer's disease.

Binding of ACE-inhibitors to in vitro and patient-derived amyloid-β fibril models

Currently, no drugs exist that can prevent or reverse Alzheimer's disease, a neurodegenerative disease associated with the presence, in the brain, of plaques that are composed of β-amyloid (Aβ) peptides. Recent studies suggest that angiotensin-converting enzyme (ACE) inhibitors, a set of drugs used to treat hypertension, may inhibit amyloid formation in vitro. In the present study, we investigate through computer simulations the binding of ACE inhibitors to patient-derived Aβ fibrils and contrast it with that of ACE inhibitors binding to in vitro generated fibrils. The binding affinities of the ACE inhibitors are compared with that of Congo red, a dye that is used to identify amyloid structures and that is known to be a weak inhibitor of Aβ aggregation. We find that ACE inhibitors have a lower binding affinity to the patient-derived fibrils than to in vitro generated ones. For patient-derived fibrils, their binding affinities are even lower than that of Congo red. Our observations raise doubts on the hypothesis that these drugs inhibit fibril formation in Alzheimer patients by interacting directly with the amyloids.

Conformational-Sensitive Fast Photochemical Oxidation of Proteins and Mass Spectrometry Characterize Amyloid Beta 1-42 Aggregation

Preventing and treating Alzheimer's disease require understanding the aggregation of amyloid beta 1-42 (Aβ1-42) to give oligomers, protofibrils, and fibrils. Here we describe footprinting of Aβ1-42 by hydroxyl radical-based fast photochemical oxidation of proteins (FPOP) and mass spectrometry (MS) to monitor the time-course of Aβ1-42 aggregation. We resolved five distinct stages characterized by two sigmoidal behaviors, showing the time-dependent transitions of monomers-paranuclei-protofibrils-fibrillar aggregates. Kinetic modeling allows deciphering the amounts and interconversion of the dominant Aβ1-42 species. Moreover, the irreversible footprinting probe provides insights into the kinetics of oligomerization and subsequent fibrillar growth by allowing the conformational changes of Aβ1-42 at subregional and even amino-acid-residue levels to be revealed. The middle domain of Aβ1-42 plays a major role in aggregation, whereas the N-terminus retains most of its solvent-accessibility during aggregation, and the hydrophobic C-terminus is involved to an intermediate extent. This approach affords an in situ, real-time monitoring of the solvent accessibility of Aβ1-42 at various stages of oligomerization, and provides new insights on site-specific aggregation of Aβ1-42 for a sample state beyond the capabilities of most other biophysical methods.

An Atomistic View of Amyloidogenic Self-assembly: Structure and Dynamics of Heterogeneous Conformational States in the Pre-nucleation Phase

The formation of well-defined filamentous amyloid structures involves a polydisperse collection of oligomeric states for which relatively little is known in terms of structural organization. Here we use extensive, unbiased explicit solvent molecular dynamics (MD) simulations to investigate the structural and dynamical features of oligomeric aggregates formed by a number of highly amyloidogenic peptides at atomistic resolution on the μs time scale. A consensus approach has been adopted to analyse the simulations in multiple force fields, yielding an in-depth characterization of pre-fibrillar oligomers and their global and local structure properties. A collision cross section analysis revealed structurally heterogeneous aggregate ensembles for the individual oligomeric states that lack a single defined quaternary structure during the pre-nucleation phase. To gain insight into the conformational space sampled in early aggregates, we probed their substructure and found emerging β-sheet subunit layers and a multitude of ordered intermolecular β-structure motifs with growing aggregate size. Among those, anti-parallel out-of-register β-strands compatible with toxic β-barrel oligomers were particularly prevalent already in smaller aggregates and formed prior to ordered fibrillar structure elements. Notably, also distinct fibril-like conformations emerged in the oligomeric state and underscore the notion that pre-nucleated oligomers serve as a critical intermediate step on-pathway to fibrils.

Picosecond dissociation of amyloid fibrils with infrared laser: A nonequilibrium simulation study

Recently, mid-infrared free-electron laser technology has been developed to dissociate amyloid fibrils. Here, we present a theoretical framework for this type of experiment based on laser-induced nonequilibrium all-atom molecular dynamics simulations. We show that the fibril is destroyed due to the strong resonance between its amide I vibrational modes and the laser field. The effects of laser irradiation are determined by a balance between fibril formation and dissociation. While the overall rearrangements of the fibril finish over short time scales, the interaction between the peptides and the solvent continues over much longer times indicating that the waters play an important role in the dissociation process. Our results thus provide new insights into amyloid fibril dissociation by laser techniques and open up new venues to investigate the complex phenomena associated with amyloidogenesis.

Discovery of DNA dyes Hoechst 34580 and 33342 as good candidates for inhibiting amyloid beta formation: in silico and in vitro study

Combining Lipinski's rule with the docking and steered molecular dynamics simulations and using the PubChem data base of about 1.4 million compounds, we have obtained DNA dyes Hoechst 34580 and Hoechst 33342 as top-leads for the Alzheimer's disease. The binding properties of these ligands to amyloid beta (Aβ) fibril were thoroughly studied by in silico and in vitro experiments. Hoechst 34580 and Hoechst 33342 prefer to locate near hydrophobic regions with binding affinity mainly governed by the van der Waals interaction. By the Thioflavin T assay, it was found that the inhibition constant IC50 ≈ 0.86 and 0.68 μM for Hoechst 34580 and Hoechst 33342, respectively. This result qualitatively agrees with the binding free energy estimated using the molecular mechanic-Poisson Boltzmann surface area method and all-atom simulations with the AMBER-f99SB-ILDN force field and water model TIP3P. In addition, DNA dyes have the high capability to cross the blood brain barrier. Thus, both in silico and in vitro experiments have shown that Hoechst 34580 and 33342 are good candidates for treating the Alzheimer's disease by inhibiting Aβ formation.

Structural and Mechanical Properties of Amyloid Beta Fibrils: A Combined Experimental and Theoretical Approach

In this combined experimental (deep ultraviolet resonance Raman (DUVRR) spectroscopy and atomic force microscopy (AFM)) and theoretical (molecular dynamics (MD) simulations and stress-strain (SS)) study, the structural and mechanical properties of amyloid beta (Aβ40) fibrils have been investigated. The DUVRR spectroscopy and AFM experiments confirmed the formation of linear, unbranched and β-sheet rich fibrils. The fibrils (Aβ40)n, formed using n monomers, were equilibrated using all-atom MD simulations. The structural properties such as β-sheet character, twist, interstrand distance, and periodicity of these fibrils were found to be in agreement with experimental measurements. Furthermore, Young's modulus (Y) = 4.2 GPa computed using SS calculations was supported by measured values of 1.79 ± 0.41 and 3.2 ± 0.8 GPa provided by two separate AFM experiments. These results revealed size dependence of structural and material properties of amyloid fibrils and show the utility of such combined experimental and theoretical studies in the design of precisely engineered biomaterials.