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

N-domain of angiotensin-converting enzyme hydrolyzes human and rat amyloid-β(1-16) peptides as arginine specific endopeptidase potentially enhancing risk of Alzheimer's disease

Alzheimer’s disease (AD) is a multifactorial neurodegenerative disorder. Amyloid-β (Aβ) aggregation is likely to be the major cause of AD. In contrast to humans and other mammals, that share the same Aβ sequence, rats and mice are invulnerable to AD-like neurodegenerative pathologies, and Aβ of these rodents (ratAβ) has three amino acid substitutions in the metal-binding domain 1-16 (MBD). Angiotensin-converting enzyme (ACE) cleaves Aβ-derived peptide substrates, however, there are contradictions concerning the localization of the cleavage sites within Aβ and the roles of each of the two ACE catalytically active domains in the hydrolysis. In the current study by using mass spectrometry and molecular modelling we have tested a set of peptides corresponding to MBDs of Aβ and ratAβ to get insights on the interactions between ACE and these Aβ species. It has been shown that the N-domain of ACE (N-ACE) acts as an arginine specific endopeptidase on the Aβ and ratAβ MBDs with C-amidated termini, thus assuming that full-length Aβ and ratAβ can be hydrolyzed by N-ACE in the same endopeptidase mode. Taken together with the recent data on the molecular mechanism of zinc-dependent oligomerization of Aβ, our results suggest a modulating role of N-ACE in AD pathogenesis.

Ferrocene-based anilides: synthesis, structural characterization and inhibition of butyrylcholinesterase

Twenty-three compounds in two series of ferrocene-based anilides, with the general formula C₅H₅-Fe-C₅H₄-C₆H₄-NH-CO-C₆H₄-R (where R = H, F, Cl, CH₃ and OCH₃), have been synthesized and found to inhibit butyrylcholinesterase.

A computational study of the Fenton reaction in different pH ranges

The mechanism of the Fenton reaction is pH dependent and four distinct reactive species have been identified and found to display quite different oxidation reactivities.

Dihydrochalcone molecules destabilize Alzheimer's amyloid-β protofibrils through binding to the protofibril cavity

Dihydrochalcone molecules destabilize Aβ_17–42protofibrils by disrupting the N-terminal β1 region and the turn region through binding to the protofibril cavity.

Investigation of the interaction of amyloid β peptide (11–42) oligomers with a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane using molecular dynamics simulation

The mechanism of pore formation in model neural cell membranes by β amyloid (Aβ) peptides was investigated using molecular dynamics simulation which indicated that Aβ oligomers of size equal or greater than 3 has a higher tendency for pore formation than monomers and that cholesterol tends to retard Aβ binding and insertion into the membrane.

Fluorescence quenching by lipid encased nanoparticles shows that amyloid-β has a preferred orientation in the membrane

Short range plasmonic fields around a nanoparticle can modulate fluorescence or Raman processes.

Bexarotene cannot reduce amyloid beta plaques through inhibition of production of amyloid beta peptides:in silicoandin vitrostudy

Recently, it has been reported that anti-cancer drug bexarotene can remarkably destroy amyloid beta (Aβ) plaques in mouse models suggesting therapeutic potential for Alzheimer's disease.

Effects of force fields on the conformational and dynamic properties of amyloid β(1-40) dimer explored by replica exchange molecular dynamics simulations

The conformational space and structural ensembles of amyloid beta (Aβ) peptides and their oligomers in solution are inherently disordered and proven to be challenging to study. Optimum force field selection for molecular dynamics (MD) simulations and the biophysical relevance of results are still unknown. We compared the conformational space of the Aβ(1-40) dimers by 300 ns replica exchange MD simulations at physiological temperature (310 K) using: the AMBER-ff99sb-ILDN, AMBER-ff99sb*-ILDN, AMBER-ff99sb-NMR, and CHARMM22* force fields. Statistical comparisons of simulation results to experimental data and previously published simulations utilizing the CHARMM22* and CHARMM36 force fields were performed. All force fields yield sampled ensembles of conformations with collision cross sectional areas for the dimer that are statistically significantly larger than experimental results. All force fields, with the exception of AMBER-ff99sb-ILDN (8.8 ± 6.4%) and CHARMM36 (2.7 ± 4.2%), tend to overestimate the α-helical content compared to experimental CD (5.3 ± 5.2%). Using the AMBER-ff99sb-NMR force field resulted in the greatest degree of variance (41.3 ± 12.9%). Except for the AMBER-ff99sb-NMR force field, the others tended to under estimate the expected amount of β-sheet and over estimate the amount of turn/bend/random coil conformations. All force fields, with the exception AMBER-ff99sb-NMR, reproduce a theoretically expected β-sheet-turn-β-sheet conformational motif, however, only the CHARMM22* and CHARMM36 force fields yield results compatible with collapse of the central and C-terminal hydrophobic cores from residues 17-21 and 30-36. Although analyses of essential subspace sampling showed only minor variations between force fields, secondary structures of lowest energy conformers are different.

Resolving the Atomistic Modes of Anle138b Inhibitory Action on Peptide Oligomer Formation

The diphenyl-pyrazole compound anle138b is a known inhibitor of oligomeric aggregate formation in vitro and in vivo. Therefore, anle138b is considered a promising drug candidate to beneficially interfere with neurodegenerative processes causing devastating pathologies in humans. The atomistic details of the aggregation inhibition mechanism, however, are to date unknown since the ensemble of small nonfibrillar aggregates is structurally heterogeneous and inaccessible to direct structural characterization. Here, we set out to elucidate anle138b's mode of action using all-atom molecular dynamics simulations on the multi-microsecond time scale. By comparing simulations of dimeric to tetrameric aggregates from fragments of four amyloidogenic proteins (Aβ, hTau40, hIAPP, and Sup35N) in the presence and absence of anle138b, we show that the compound reduces the overall number of intermolecular hydrogen bonds, disfavors the sampling of the aggregated state, and remodels the conformational distributions within the small oligomeric peptide aggregates. Most notably, anle138b preferentially interacts with the disordered structure ensemble via its pyrazole moiety, thereby effectively blocking interpeptide main chain interactions and impeding the spontaneous formation of ordered β-sheet structures, in particular those with out-of-register antiparallel β-strands. The structurally very similar compound anle234b was previously identified as inactive by in vitro experiments. Here, we show that anle234b has no significant effect on the aggregation process in terms of reducing the β-structure content. Moreover, we demonstrate that the hydrogen bonding capabilities are autoinhibited due to steric effects imposed by the molecular geometry of anle234b and thereby indirectly confirm the proposed inhibitory mechanism of anle138b. We anticipate that the prominent binding of anle138b to partially disordered and dynamical aggregate structures is a generic basis for anle138b's ability to suppress toxic oligomer formation in a wide range of amyloidogenic peptides and proteins.

Chemistry of mammalian metallothioneins and their interaction with amyloidogenic peptides and proteins

Cu and Zn ions are essential in most living beings. Their metabolism is critical for health and mis-metabolism can be lethal. In the last two decades, a large body of evidence has reported the role of copper, zinc and iron, and oxidative stress in several neurodegenerative diseases like Alzheimer's, Parkinson's, prion diseases, etc. To what extent this mis-metabolism is causative or a consequence of these diseases is still a matter of research. In this context metallothioneins (MTs) appear to play a central gate-keeper role in controlling aberrant metal-protein interactions. MTs are small proteins that can bind high amounts of Zn(ii) and Cu(i) ions in metal-cluster arrangements via their cysteine thiolates. Moreover, MTs are well known antioxidants. The present tutorial outlines the chemistry underlying the interconnection between copper(i/ii) and zinc(ii) coordination to amyloidogenic proteins and MTs, and their redox properties in generation and/or silencing reactive oxygen species (overproduced in oxidative stress) and other reactants. These studies have revealed the coordination chemistry involved in neurodegenerative diseases and the interactions between MTs and amyloidogenic protein metal-complexes (like amyloid-β, α-synuclein and prion-protein). Overall, the protective role of MTs in neurodegenerative processes is emerging, serving as a foundation for exploring MT chemistry as inspiration for therapeutic approaches.

Thermodynamic properties of amyloid fibrils in equilibrium

In this manuscript we use a two-dimensional coarse-grained model to study how amyloid fibrils grow towards an equilibrium state where they coexist with proteins dissolved in a solution. Free-energies to dissociate proteins from fibrils are estimated from the residual concentration of dissolved proteins. Consistent with experiments, the concentration of proteins in solution affects the growth rate of fibrils but not their equilibrium state. Also, studies of the temperature dependence of the equilibrium state can be used to estimate thermodynamic quantities, e.g., heat capacity and entropy.

Initial Structural Models of the Aβ42 Dimer from Replica Exchange Molecular Dynamics Simulations

Experimental characterization of the molecular structure of small amyloid (A)β oligomers that are currently considered as toxic agents in Alzheimer's disease is a formidably difficult task due to their transient nature and tendency to aggregate. Such structural information is of importance because it can help in developing diagnostics and an effective therapy for the disease. In this study, molecular simulations and protein-protein docking are employed to explore a possible connection between the structure of Aβ monomers and the properties of the intermonomer interface in the Aβ42 dimer. A structurally diverse ensemble of conformations of the monomer was sampled in microsecond timescale implicit solvent replica exchange molecular dynamics simulations. Representative structures with different solvent exposure of hydrophobic residues and secondary structure content were selected to build structural models of the dimer. Analysis of these models reveals that formation of an intramonomer salt bridge (SB) between Asp23 and Lys28 residues can prevent the building of a hydrophobic interface between the central hydrophobic clusters (CHCs) of monomers upon dimerization. This structural feature of the Aβ42 dimer is related to the difference in packing of hydrophobic residues in monomers with the Asp23-Lys28 SB in on and off states, in particular, to a lower propensity to form hydrophobic contacts between the CHC domain and C-terminal residues in monomers with a formed SB. These findings could have important implications for understanding the difference between aggregation pathways of Aβ monomers leading to neurotoxic oligomers or inert fibrillar structures.

Comparing the Aggregation Free Energy Landscapes of Amyloid Beta(1-42) and Amyloid Beta(1-40)

Using a predictive coarse-grained protein force field, we compute and compare the free energy landscapes and relative stabilities of amyloid-β protein (1-42) and amyloid-β protein (1-40) in their monomeric and oligomeric forms up to the octamer. At the same concentration, the aggregation free energy profile of Aβ42 is more downhill, with a computed solubility that is about 10 times smaller than that of Aβ40. At a concentration of 40 μM, the clear free energy barrier between the pre-fibrillar tetramer form and the fibrillar pentamer in the Aβ40 aggregation landscape disappears for Aβ42, suggesting that the Aβ42 tetramer has a more diverse structural range. To further compare the landscapes, we develop a cluster analysis based on the structural similarity between configurations and use it to construct an oligomerization map that captures the paths of easy interconversion between different but structurally similar states of oligomers for both species. A taxonomy of the oligomer species based on β-sheet stacking topologies is proposed. The comparison of the two oligomerization maps highlights several key differences in the landscapes that can be attributed to the two additional C-terminal residues that Aβ40 lacks. In general, the two terminal residues strongly stabilize the oligomeric structures for Aβ42 relative to Aβ40, and greatly facilitate the conversion from pre-fibrillar trimers to fibrillar tetramers.

How accurate are your simulations? Effects of confined aqueous volume and AMBER FF99SB and CHARMM22/CMAP force field parameters on structural ensembles of intrinsically disordered proteins: Amyloid-β42 in water

Amyloid-β₄₂ (Aβ₄₂) is an intrinsically disordered peptide intimately related to the pathogenesis of several neurodegenerative diseases. Molecular dynamics (MD) simulations are extensively utilized in the characterization of the structures and conformational dynamics of intrinsically disordered proteins (IDPs) including Aβ₄₂, with AMBER and CHARMM parameters being commonly used in these studies. Recently, comparison of the effects of force field parameters on the Aβ₄₂ structures has started to gain significant attention. In this study, the structures of Aβ₄₂ are simulated using AMBER FF99SB and CHARMM22/CMAP parameters via replica exchange MD simulations utilizing a widely used clustering algorithm. These analyses show that the structural properties (extent and positioning of the elements of secondary and tertiary structure), radius of gyration values, number and position of salt bridges are extremely dependent on the chosen force field parameters notably with the usage of clustering algorithms. For example, predicted secondary structure elements, which are of the great importance for better understanding of the molecular mechanisms of neurodegenerative diseases, deviate enormously in models generated using currently available force field parameters for proteins. Based on the derived models, chemical shift values are calculated and compared to the experimentally determined data. This comparison revealed that although both force field parameters yield results in agreement with experiments, the obtained structural properties were rather different using a clustering algorithm. In other words, these results show that the predicted structures depend heavily on the force field parameters. Importantly, since none of the force field parameters currently utilized in MD studies were developed specifically taking into account the disordered nature of IDPs, these findings clearly indicate that new force field parameters have to be developed for IDPs considering their rapid flexibility and dynamics with high amplitude. Furthermore, molecular simulations of IDPs are typically conducted using one water volume. We show that the confined aqueous volume impacts the predicted structural properties of Aβ₄₂ in water. Although up to date, confined aqueous volume effects have been ignored in the MD simulations of IDPs in water, our data indicate that these effects have to be taken into account in predicting the structural and thermodynamic properties of disordered proteins in solution.

Binding of protofibrillar Aβ trimers to lipid bilayer surface enhances Aβ structural stability and causes membrane thinning

Alzheimer's disease, a common neurodegenerative disease, is characterized by the aggregation of amyloid-β (Aβ) peptides. The interactions of Aβ with membranes cause changes in membrane morphology and ion permeation, which are responsible for its neurotoxicity and can accelerate fibril growth. However, the Aβ-lipid interactions and how these induce membrane perturbation and disruption at the atomic level and the consequences for the Aβ organization are not entirely understood. Here, we perform multiple atomistic molecular dynamics simulations on three protofibrillar Aβ9-40 trimers. Our simulations show that, regardless of the morphologies and the initial orientations of the three different protofibrillar Aβ9-40 trimers, the N-terminal β-sheet of all trimers preferentially binds to the membrane surface. The POPG lipid bilayers enhance the structural stability of protofibrillar Aβ trimers by stabilizing inter-peptide β-sheets and D23-K28 salt-bridges. The interaction causes local membrane thinning. We found that the trimer structure related to Alzheimer's disease brain tissue () is the most stable both in water solution and at membrane surface, and displays slightly stronger membrane perturbation capability. These results provide mechanistic insights into the membrane-enhanced structural stability of protofibrillar Aβ oligomers and the first step of Aβ-induced membrane disruption at the atomic level.

Direct Evidence of Intrinsic Blue Fluorescence from Oligomeric Interfaces of Human Serum Albumin

The molecular origin behind the concentration-dependent intrinsic blue fluorescence of human serum albumin (HSA) is not known yet. This unusual blue fluorescence is believed to be a characteristic feature of amyloid-like fibrils of protein/peptide and originates due to the delocalization of peptide bond electrons through the extended hydrogen bond networks of cross-β-sheet structure. Herein, by combining the results of spectroscopy, size exclusion chromatography, native gel electrophoresis, and confocal microscopy, we have shown that the intrinsic blue fluorescence of HSA exclusively originates from oligomeric interfaces devoid of any amyloid-like fibrillar structure. Our study suggests that this low energy fluorescence band is not due to any particular residue/sequence, but rather it is a common feature of self-assembled peptide bonds. The present findings of intrinsic blue fluorescence from oligomeric interfaces pave the way for future applications of this unique visual phenomenon for early stage detection of various protein aggregation related human diseases.

Using chirality to probe the conformational dynamics and assembly of intrinsically disordered amyloid proteins

Intrinsically disordered protein (IDP) conformers occupy large regions of conformational space and display relatively flat energy surfaces. Amyloid-forming IDPs, unlike natively folded proteins, have folding trajectories that frequently involve movements up shallow energy gradients prior to the “downhill” folding leading to fibril formation. We suggest that structural perturbations caused by chiral inversions of amino acid side-chains may be especially valuable in elucidating these pathways of IDP folding. Chiral inversions are subtle in that they do not change side-chain size, flexibility, hydropathy, charge, or polarizability. They allow focus to be placed solely on the question of how changes in amino acid side-chain orientation, and the resultant alterations in peptide backbone structure, affect a peptide’s conformational landscape (Ramachandran space). If specific inversions affect folding and assembly, then the sites involved likely are important in mediating these processes. We suggest here a “focused chiral mutant library” approach for the unbiased study of amyloid-forming IDPs.

Amyloid β Fibril Elongation by Monomers Involves Disorder at the Tip

The growth of amyloid fibrils from Aβ_1-42 peptide, one of the key pathogenic players in Alzheimer's disease, is believed to follow a nucleation-elongation mechanism. Fibril elongation is often described as a "dock-lock" procedure, where a disordered monomer adsorbs to an existing fibril in a relatively fast process (docking), followed by a slower conformational transition toward the ordered state of the template (locking). Here, we use molecular dynamics simulations of an ordered pentamer of Aβ42 at fully atomistic resolution, which includes solvent, to characterize the elongation process. We construct a Markov state model from an ensemble of short trajectories generated by an advanced sampling algorithm that efficiently diversifies a subset of the system without any bias forces. This subset corresponds to selected dihedral angles of the peptide chain at the fibril tip favored to be the fast growing one experimentally. From the network model, we extract distinct locking pathways covering time scales in the high microsecond regime. Slow steps are associated with the exchange of hydrophobic contacts, between nonnative and native intermolecular contacts as well as between intra- and intermolecular ones. The N-terminal segments, which are disordered in fibrils and typically considered inert, are able to shield the lateral interfaces of the pentamer. We conclude by discussing our findings in the context of a refined dock-lock model of Aβ fibril elongation, which involves structural disorder for more than one monomer at the growing tip.

Bexarotene Does Not Clear Amyloid Beta Plaques but Delays Fibril Growth: Molecular Mechanisms

In 2012, it was reported that anticancer drug bexarotene reduced amyloid plaque and improved mental functioning in a small sample of mice engineered to exhibit Alzheimer's like symptoms. It has been suggested that bexarotene stimulates expression of apolipoprotein E (ApoE) leading to intracellular clearance of amyloid beta (Aβ). However, the effect of bexarotene on clearance of plaques has not been seen in some mouse models. Two interesting questions include whether bexarotene can destroy Aβ fibrils via direct interaction with them and how this compound impacts the lag phase in the fibril growth process. By the Thioflavin T fluorescence assay and atomic force microscopy, we have shown that bexarotene prolongs the lag phase, but it does not degrade Aβ fibrils. The impotence of bexarotene in destroying fibrils means that this compound is weakly bound to Aβ. On the other hand, the weak binding would prevent bexarotene from prolonging the lag phase. Thus, our two main in vitro observations seem to contradict each other. In order to settle this problem at the atomic level, we have performed all-atom molecular dynamics simulations in explicit water. We have demonstrated that bexarotene is not capable to reduce amyloid deposits due to weak binding to Aβ fibrils. However, it delays the self-assembly through reduction of the β-content of Aβ monomers at high enough ligand concentrations. Bexarotene is the first compound which displays such an unusual behavior. We have also shown that bexarotene has a low binding propensity to Aβ monomer and dimer.

Multi-scale simulations of biological systems using the OPEP coarse-grained model

Biomolecules are complex machines that are optimized by evolution to properly fulfill or contribute to a variety of biochemical tasks in the cellular environment. Computer simulations based on quantum mechanics and atomistic force fields have been proven to be a powerful microscope for obtaining valuable insights into many biological, physical, and chemical processes. Many interesting phenomena involve, however, a time scale and a number of degrees of freedom, notably if crowding is considered, that cannot be explored at an atomistic resolution. To bridge the gap between reality and simulation, many different advanced computational techniques and coarse-grained (CG) models have been developed. Here, we report some applications of the CG OPEP protein model to amyloid fibril formation, the response of catch-bond proteins to two types of fluid flow, and interactive simulations to fold peptides with well-defined 3D structures or with intrinsic disorder.

Biocompatible and blood-brain barrier permeable carbon dots for inhibition of Aβ fibrillation and toxicity, and BACE1 activity

Amyloid-β peptide (Aβ) fibrillation is pathologically associated with Alzheimer's disease (AD), and this has resulted in the development of an Aβ inhibitor which is essential for the treatment of AD. However, the design of potent agents which can target upstream secretases, inhibit Aβ toxicity and aggregation, as well as cross the blood-brain barrier remains challenging. In, this research carbon dots for AD treatment were investigated in vitro using experimental and computational methods for the first time. The results presented here demonstrate a novel strategy for the discovery of novel antiamyloidogenic agents for AD treatments.

Major Reaction Coordinates Linking Transient Amyloid-β Oligomers to Fibrils Measured at Atomic Level

The structural underpinnings for the higher toxicity of the oligomeric intermediates of amyloidogenic peptides, compared to the mature fibrils, remain unknown at present. The transient nature and heterogeneity of the oligomers make it difficult to follow their structure. Here, using vibrational and solid-state nuclear magnetic resonance spectroscopy, and molecular dynamics simulations, we show that freely aggregating Aβ₄₀ oligomers in physiological solutions have an intramolecular antiparallel configuration that is distinct from the intermolecular parallel β-sheet structure observed in mature fibrils. The intramolecular hydrogen-bonding network flips nearly 90°, and the two β-strands of each monomeric unit move apart, to give rise to the well-known intermolecular in-register parallel β-sheet structure in the mature fibrils. Solid-state nuclear magnetic resonance distance measurements capture the interstrand separation within monomer units during the transition from the oligomer to the fibril form. We further find that the D23-K28 salt-bridge, a major feature of the Aβ₄₀ fibrils and a focal point of mutations linked to early onset Alzheimer's disease, is not detectable in the small oligomers. Molecular dynamics simulations capture the correlation between changes in the D23-K28 distance and the flipping of the monomer secondary structure between antiparallel and parallel β-sheet architectures. Overall, we propose interstrand separation and salt-bridge formation as key reaction coordinates describing the structural transition of the small Aβ₄₀ oligomers to fibrils.

Membrane-Accelerated Amyloid-β Aggregation and Formation of Cross-β Sheets

Amyloid- β aggregates play a causative role in Alzheimer's disease. These aggregates are a product of the physical environment provided by the basic neuronal membrane, composed of a lipid bilayer. The intrinsic properties of the lipid bilayer allow amyloid- β peptides to nucleate and form well-ordered cross- β sheets within the membrane. Here, we correlate the aggregation of the hydrophobic fragment of the amyloid- β protein, A β 25 - 35 , with the hydrophobicity, fluidity, and charge density of a lipid bilayer. We summarize recent biophysical studies of model membranes and relate these to the process of aggregation in physiological systems.

Protocol for fast screening of multi-target drug candidates: Application to Alzheimer's disease

The treatment of many diseases may require drugs that are capable to attack multiple targets simultaneously. Obviously, the virtual screening of multi-target drug candidates is much more time consuming compared to the single-target case. This, in particular, concerns the last step of virtual screening where the binding free energy is computed by conventional molecular dynamics simulation. To overcome this difficulty we propose a simple protocol which is relied on the fast steered molecular dynamics simulation and on available experimental data on binding affinity of reference ligand to a given target. Namely, first we compute non-equilibrium works generated during pulling ligands from the binding site using the steered molecular dynamics method. Then as top leads we choose only those compounds that have the non-equilibrium work larger than that of a reference compound for which the binding free energy has been already known from experiment. Despite many efforts no cures for AD (Alzheimer's disease) have been found. One of possible reasons for this failure is that drug candidates were developed for a single target, while there are exist many possible pathways to AD. Applying our new protocol to five targets including amyloid beta fibril, peroxisome proliferator-activated receptor γ, retinoic X receptor α, β- and γ-secretases, we have found two potential drugs (CID 16040294 and CID 9998128) for AD from the large PubChem database. We have also shown that these two ligands can interfere with the activity of popular Acetylcholinesterase target through strong binding towards it.

Mechanism of Initiation, Association, and Formation of Amyloid Fibrils Modeled with the N-Terminal Peptide Fragment, IKYLEFIS, of Myoglobin G-Helix

Some peptides and proteins undergo self-aggregation under certain conditions, leading to amyloid fibrils formation, which is related to many disease conditions. It is important to understand such amyloid fibrils formation to provide mechanistic detail that governs the process. A predominantly α-helical myoglobin has been reported recently to readily form amyloid fibrils at a higher temperature, similar to its G-helix segment. Here, we have investigated the mechanism of amyloid fibrils formation by performing multiple long molecular dynamics simulations (27 μs) on the N-terminal segment of the G-helix of myoglobin. These simulations resulted in the formation of a single-layered tetrameric β-sheet with mixed parallel and antiparallel β-strands and this is the most common event irrespective of many different starting structures. Formation of the single-layered tetrameric β-sheet takes place following three distinctive pathways. The process of fibril initiation is dependent on temperature. Further, this study provides mechanistic insights into the formation of multilayered fibrilar structure, which could be applicable to a wider variety of peptides or proteins to understand the amyloidogenesis.

Probing oligomerization of amyloid beta peptide in silico

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

Structural studies of amyloid-β peptides: Unlocking the mechanism of aggregation and the associated toxicity

Alzheimer's disease (AD) is one of the most prevalent neurodegenerative diseases worldwide. Formation of amyloid plaques consisting of amyloid-β peptides (Aβ) is one of the hallmarks of AD. Several lines of evidence have shown a correlation between the Aβ aggregation and the disease development. Extensive research has been conducted with the aim to reveal the structures of the neurotoxic Aβ aggregates. However, the exact structure of pathological aggregates and mechanism of the disease still remains elusive due to complexity of the occurring processes and instability of various disease-relevant Aβ species. In this article we review up-to-date structural knowledge about amyloid-β peptides, focusing on data acquired using solution and solid state NMR techniques. Furthermore, we discuss implications from these structural studies on the mechanisms of aggregation and neurotoxicity.

Nanopore Sensing of Protein Folding

Single-molecule studies of protein folding hold keys to unveiling protein folding pathways and elusive intermediate folding states-attractive pharmaceutical targets. Although conventional single-molecule approaches can detect folding intermediates, they presently lack throughput and require elaborate labeling. Here, we theoretically show that measurements of ionic current through a nanopore containing a protein can report on the protein's folding state. Our all-atom molecular dynamics (MD) simulations show that the unfolding of a protein lowers the nanopore ionic current, an effect that originates from the reduction of ion mobility in proximity to a protein. Using a theoretical model, we show that the average change in ionic current produced by a folding-unfolding transition is detectable despite the orientational and conformational heterogeneity of the folded and unfolded states. By analyzing millisecond-long all-atom MD simulations of multiple protein transitions, we show that a nanopore ionic current recording can detect folding-unfolding transitions in real time and report on the structure of folding intermediates.

Why Is Research on Amyloid-β Failing to Give New Drugs for Alzheimer's Disease?

The two hallmarks of Alzheimer's disease (AD) are the presence of neurofibrillary tangles (NFT) made of aggregates of the hyperphosphorylated tau protein and of amyloid plaques composed of amyloid-β (Aβ) peptides, primarily Aβ1-40 and Aβ1-42. Targeting the production, aggregation, and toxicity of Aβ with small molecule drugs or antibodies is an active area of AD research due to the general acceptance of the amyloid cascade hypothesis, but thus far all drugs targeting Aβ have failed. From a review of the recent literature and our own experience based on in vitro, in silico, and in vivo studies, we present some reasons to explain this repetitive failure.

Molecular Dynamics Simulations of Amyloid β-Peptide (1-42): Tetramer Formation and Membrane Interactions

The aggregation cascade and peptide-membrane interactions of the amyloid β-peptide (Aβ) have been implicated as toxic events in the development and progression of Alzheimer's disease. Aβ42 forms oligomers and ultimately plaques, and it has been hypothesized that these oligomeric species are the main toxic species contributing to neuronal cell death. To better understand oligomerization events and subsequent oligomer-membrane interactions of Aβ42, we performed atomistic molecular-dynamics (MD) simulations to characterize both interpeptide interactions and perturbation of model membranes by the peptides. MD simulations were utilized to first show the formation of a tetramer unit by four separate Aβ42 peptides. Aβ42 tetramers adopted an oblate ellipsoid shape and showed a significant increase in β-strand formation in the final tetramer unit relative to the monomers, indicative of on-pathway events for fibril formation. The Aβ42 tetramer unit that formed in the initial simulations was used in subsequent MD simulations in the presence of a pure POPC or cholesterol-rich raft model membrane. Tetramer-membrane simulations resulted in elongation of the tetramer in the presence of both model membranes, with tetramer-raft interactions giving rise to the rearrangement of key hydrophobic regions in the tetramer and the formation of a more rod-like structure indicative of a fibril-seeding aggregate. Membrane perturbation by the tetramer was manifested in the form of more ordered, rigid membranes, with the pure POPC being affected to a greater extent than the raft membrane. These results provide critical atomistic insight into the aggregation pathway of Aβ42 and a putative toxic mechanism in the pathogenesis of Alzheimer's disease.

Light-induced self-assembly of bi-color CdTe quantum dots allows the discrimination of multiple proteins

We have found that the addition of proteins can greatly influence the light-induced self-assembly (LISA) behavior of bi-color thioglycolic acid (TGA)-capped CdTe Quantum Dots (QDs) and thus cause significant changes of their fluorescence (FL) signals (color and intensity), according to which a dual-channel FL sensor can be established for simultaneous discrimination of multiple proteins. The sensor is successfully used for the identification of ten native proteins and ten thermally denatured proteins and eight native proteins artificially added in human urine, respectively, during which process principal component analysis (PCA) is utilized to differentiate the targets based on their corresponding FL change patterns. This assay has provided a visual and simple method for the discrimination of various analytes, which may have great potential in the study of conformational changes of biomacromolecules and the analysis of real biological fluids.

A Binuclear Zinc Interaction Fold Discovered in the Homodimer of Alzheimer's Amyloid-β Fragment with Taiwanese Mutation D7H

Zinc-induced oligomerization of amyloid-β peptide (Aβ) produces potentially pathogenic agents of Alzheimer's disease. Mutations and modifications in the metal binding domain 1-16 of Aβ peptide crucially affect its zinc-induced oligomerization by changing intermolecular zinc mediated interface. The 3D structure of this interface appearing in a range of Aβ species is a prospective drug target for disease modifying therapy. Using NMR spectroscopy, EXAFS spectroscopy, mass spectrometry, and isothermal titration calorimetry the interaction of zinc ions with Aβ fragments 1-7 and 1-10 carrying familial Taiwanese mutation D7H was studied. Zinc ions induce formation of a stable homodimer formed by the two peptide chains fastened by two zinc ions and stacking interactions of imidazole rings. A binuclear zinc interaction fold in the dimer structure was discovered. It can be used for designing zinc-regulated proteins and zinc-mediated self-assembling peptides.

Study of structural stability and damaging effect on membrane for four Aβ42 dimers

Increasing evidence shows that Aβ oligomers are key pathogenic molecules in Alzheimer’s disease. Among Aβ oligomers, dimer is the smallest aggregate and toxic unit. Therefore, understanding its structural and dynamic properties is quite useful to prevent the formation and toxicity of the Aβ oligomers. In this study, we performed molecular dynamic simulations on four Aβ42 dimers, 2NCb, CNNC, NCNC and NCCN, within the hydrated DPPC membrane. Four Aβ42 dimers differ in the arrangements of two Aβ42 peptides. This study aims to investigate the impact of aggregation pattern of two Aβ peptides on the structural stability of the Aβ42 dimer and its disruption to the biological membrane. The MD results demonstrate that the NCCN, CNNC and NCNC have the larger structural fluctuation at the N-terminus of Aβ42 peptide, where the β-strand structure converts into the coil structure. The loss of the N-terminal β-strand further impairs the aggregate ability of Aβ42 dimer. In addition, inserting Aβ42 dimer into the membrane can considerably decrease the average APL of DPPC membrane. Moreover this decrease effect is largely dependent on the distance to the location of Aβ42 dimer and its secondary structure forms. Based on the results, the 2NCb is considered as a stable dimeric unit for aggregating the larger Aβ42 oligomer, and has a potent ability to disrupt the membrane.

Probing supramolecular protein assembly using covalently attached fluorescent molecular rotors

Changes in microscopic viscosity and macromolecular crowding accompany the transition of proteins from their monomeric forms into highly organised fibrillar states. Previously, we have demonstrated that viscosity sensitive fluorophores termed 'molecular rotors', when freely mixed with monomers of interest, are able to report on changes in microrheology accompanying amyloid formation, and measured an increase in rigidity of approximately three orders of magnitude during aggregation of lysozyme and insulin. Here we extend this strategy by covalently attaching molecular rotors to several proteins capable of assembly into fibrils, namely lysozyme, fibrinogen and amyloid-β peptide (Aβ(1-42)). We demonstrate that upon covalent attachment the molecular rotors can successfully probe supramolecular assembly in vitro. Importantly, our new strategy has wider applications in cellulo and in vivo, since covalently attached molecular rotors can be successfully delivered in situ and will colocalise with the aggregating protein, for example inside live cells. This important advantage allowed us to follow the microscopic viscosity changes accompanying blood clotting and during Aβ(1-42) aggregation in live SH-SY5Y cells. Our results demonstrate that covalently attached molecular rotors are a widely applicable tool to study supramolecular protein assembly and can reveal microrheological features of aggregating protein systems both in vitro and in cellulo not observable through classical fluorescent probes operating in light switch mode.

High-Resolution Structures of the Amyloid-β 1-42 Dimers from the Comparison of Four Atomistic Force Fields

The dimer of the amyloid-β peptide Aβ of 42 residues is the smallest toxic species in Alzheimer's disease, but its equilibrium structures are unknown. Here we determined the equilibrium ensembles generated by the four atomistic OPLS-AA, CHARMM22*, AMBER99sb-ildn, and AMBERsb14 force fields with the TIP3P water model. On the basis of 144 μs replica exchange molecular dynamics simulations (with 750 ns per replica), we find that the four force fields lead to random coil ensembles with calculated cross-collision sections, hydrodynamics properties, and small-angle X-ray scattering profiles independent of the force field. There are, however, marked differences in secondary structure, with the AMBERsb14 and CHARMM22* ensembles overestimating the CD-derived helix content, and the OPLS-AA and AMBER99sb-ildn secondary structure contents in agreement with CD data. Also the intramolecular beta-hairpin content spanning residues 17-21 and 30-36 varies between 1.5% and 13%. Overall, there are significant differences in tertiary and quaternary conformations among all force fields, and the key finding, irrespective of the force field, is that the dimer is stabilized by nonspecific interactions, explaining therefore its possible transient binding to multiple cellular partners and, in part, its toxicity.

Influence of Hyperglycemic Conditions on Self-Association of the Alzheimer's Amyloid β (Aβ1-42) Peptide

Clinical studies have identified a correlation between type-2 diabetes mellitus and cognitive decrements en route to the onset of Alzheimer's disease (AD). Recent studies have established that post-translational modifications of the amyloid β (Aβ) peptide occur under hyperglycemic conditions; particularly, the process of glycation exacerbates its neurotoxicity and accelerates AD progression. In view of the assertion that macromolecular crowding has an altering effect on protein self-assembly, it is crucial to characterize the effects of hyperglycemic conditions via crowding on Aβ self-assembly. Toward this purpose, fully atomistic molecular dynamics simulations were performed to study the effects of glucose crowding on Aβ dimerization, which is the smallest known neurotoxic species. The dimers formed in the glucose-crowded environment were found to have weaker associations as compared to that of those formed in water. Binding free energy calculations show that the reduced binding strength of the dimers can be mainly attributed to the overall weakening of the dispersion interactions correlated with substantial loss of interpeptide contacts in the hydrophobic patches of the Aβ units. Analysis to discern the differential solvation pattern in the glucose-crowded and pure water systems revealed that glucose molecules cluster around the protein, at a distance of 5-7 Å, which traps the water molecules in close association with the protein surface. This preferential exclusion of glucose molecules and resulting hydration of the Aβ peptides has a screening effect on the hydrophobic interactions, which in turn diminishes the binding strength of the resulting dimers. Our results imply that physical effects attributed to crowded hyperglycemic environments are incapable of solely promoting Aβ self-assembly, indicating that further mechanistic studies are required to provide insights into the self-assembly of post-translationally modified Aβ peptides, known to possess aggravated toxicity, under these conditions.

Aβ42 pentamers/hexamers are the smallest detectable oligomers in solution

Amyloid β (Aβ) oligomers may play a decisive role in Alzheimer's disease related neurodegeneration, but their structural properties are poorly understood. In this report, sedimentation velocity centrifugation, small angle neutron scattering (SANS) and molecular modelling were used to identify the small oligomeric species formed by the 42 amino acid residue long isoform of Aβ (Aβ42) in solution, characterized by a sedimentation coefficient of 2.56 S, and a radius of gyration between 2 and 4 nm. The measured sedimentation coefficient is in close agreement with the sedimentation coefficient calculated for Aβ42 hexamers using MD simulations at µM concentration. To the best of our knowledge this is the first report detailing the Aβ42 oligomeric species by SANS measurements. Our results demonstrate that the smallest detectable species in solution are penta- to hexamers. No evidences for the presence of dimers, trimers or tetramers were found, although the existence of those Aβ42 oligomers at measurable quantities had been reported frequently.