Exploring the Alzheimer amyloid-β peptide conformational ensemble: A review of molecular dynamics approaches

Alzheimer's disease is one of the most common dementia among elderly worldwide. There is no therapeutic drugs until now to treat effectively this disease. One main reason is due to the poorly understood mechanism of Aβ peptide aggregation, which plays a crucial role in the development of Alzheimer's disease. It remains challenging to experimentally or theoretically characterize the secondary and tertiary structures of the Aβ monomer because of its high flexibility and aggregation propensity, and its conformations that lead to the aggregation are not fully identified. In this review, we highlight various structural ensembles of Aβ peptide revealed and characterized by computational approaches in order to find converging structures of Aβ monomer. Understanding how Aβ peptide forms transiently stable structures prior to aggregation will contribute to the design of new therapeutic molecules against the Alzheimer's disease.     

Staphylococcal display for combinatorial protein engineering of a head-to-tail affibody dimer binding the Alzheimer amyloid-β peptide

We have previously generated an affibody molecule for the disease-associated amyloid beta (Aβ) peptide, which has been shown to inhibit the formation of various Aβ aggregates and revert the neurotoxicity of Aβ in a fruit fly model of Alzheimer's disease. In this study, we have investigated a new bacterial display system for combinatorial protein engineering of the Aβ-binder as a head-to-tail dimeric construct for future optimization efforts, e.g. affinity maturation. Using the bacterial display platform, we have: (i) demonstrated functional expression of the dimeric binder on the cell surface, (ii) determined the affinity and investigated the pH sensitivity of the interaction, (iii) demonstrated the importance of an intramolecular disulfide bond through selections from a cell-displayed combinatorial library, as well as (iv) investigated the effects from rational truncation of the N-terminal part of the affibody molecule on surface expression level and Aβ binding. Overall, the detailed engineering and characterization of this promising Aβ-specific affibody molecule have yielded valuable insights concerning its unusual binding mechanism. The results also demonstrated that our bacterial display system is a suitable technology for future protein engineering and characterization efforts of homo- or heterodimeric affinity proteins.     

Carbon nanotube prevents the secondary structure formation of amyloid-β trimers: an all-atom molecular dynamics study

Abstract

Increasing evidence shows that the formation of misfolded aggregates amyloid-β (Aβ) peptide is associated with the Alzheimer’s disease (AD). Recent experiments reveal a significant correlation of the amount of trimer species bound to neurons with increasing neuro-toxicity. Our previous computational study demonstrated that carbon nanotubes (CNT) can inhibit effectively the formation of β-sheet-rich oligomers of Aβ(16-22) – a hydrophobic peptide essential for Aβ fibrillization. However, the influence of CNT on the oligomers formed by full-length Aβ remains elusive. In this study, we have investigated the conformational dynamics of Aβ(1-42) trimer, built from an NMR structure of α-helical monomer, in the absence and presence of a single-walled carbon nanotube (SWCNT). Our simulations demonstrate that SWCNT can significantly hinder the trimerisation and prevents the secondary structure formation of Aβ(1-42) peptide. Hydrophobic and aromatic stacking interactions between SWCNT and Aβ play important roles in the secondary structure formation of the Aβ trimer. This study reveals a complete picture of the detailed preventable mechanism of full-length Aβ(1-42) by SWCNT, providing theoretical insights into the development of drug candidates of AD.     

Prediction of the binding mode between GSK3β and a peptide derived from GSKIP using molecular dynamics simulation

GSK3β plays an important role in many physiological functions; dysregulated GSK3β is involved in human diseases such as diabetes, cancer, and Alzheimer's disease. This study uses MD simulations to determine the interaction between GSK3β and a peptide derived from GSKIP, a novel GSK3β interacting protein. Results show that GSKIPtide is inlaid in a binding pocket consisting of an α-helix and an extended loop near the carboxy-terminal end. This binding pocket is hydrophobic, and is responsible for the protein-protein interaction of two other GSK3β interacting proteins: FRAT and Axin. The GSKIPtide binding mode is closer to that of AxinGID (in the Axin-GSK3-interacting domain). The single-point mutations of V267G and Y288F in GSK3β differentiate the binding modes between GSK3 and GSKIPtide, AxinGID, and FRATide. The V2677G mutation of GSK3β reduces the GSKIPtide binding affinity by 70% and abolishes the binding affinity with AxinGID, but has no effect on FRATide. However, GSK3β Y288F completely abolishes the FRATide binding without affecting GSKIPtide or AxinGID binding. An analysis of the GSK3β-GSKIPtide complex structure and the X-ray crystal structures of GSK3β-FRATide and GSK3β-AxinGID complexes suggests that the hydroxyl group of Y288 is crucial to maintaining a hydrogen bond network in GSK3β-FRATide. The hydrophobic side chain of V267 maintains the integrity of helix-helix ridge-groove hydrophobic interaction for GSK3β-GSKIPtide and GSK3β-AxinGID. This study simulates these two mutant systems to provide atomic-level evidence of the aforementioned experimental results and validate the wild-type complex structure prediction.     

Computational study of the binding of CuII to Alzheimer’s amyloid-β peptide: Do Aβ42 and Aβ40 bind copper in identical fashion?

One of the many hypotheses on the pathogenesis of Alzheimer’s disease is that the amyloid-β peptide (Aβ) binds Cu^II and can catalytically generate H₂O₂, leading to oxidative damage in brain tissues. For a molecular level understanding of such catalysis it is critical to know the structure of the Aβ–Cu^II complex precisely. Unfortunately, no high-resolution structure is available to date and there is considerable debate over the copper coordination environment with no clear consensus on which residues are directly bound to Cu^II. Considering all plausible isomers of the copper-bound Aβ42 and Aβ40 using a combination of density functional theory and classical molecular dynamics methods, we report an atomic resolution structure for each possible complex. We evaluated the relative energies of these isomeric structures and surprisingly found that Aβ42 and Aβ40 display very different binding modes, suggesting that shorter peptides that are truncated at the C-terminus may not be realistic models for understanding the chemistry of the most neurotoxic peptide, Aβ42. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

On the Origin of the Stronger Binding of PIB over Thioflavin T to Protofibrils of the Alzheimer Amyloid-β Peptide: A Molecular Dynamics Study

Pittsburgh compound B (PIB) is a neutral derivative of the fluorescent dye Thioflavin T (ThT), which displays enhanced hydrophobicity and binding affinity to amyloid fibrils. We present molecular dynamics simulations of binding of PIB and ThT to a common cross-β-subunit of the Alzheimer Amyloid-β peptide (Aβ). Our simulations of binding to Aβ_9-40 protofibrils show that PIB, like ThT, selectively binds to the hydrophobic or aromatic surface grooves on the β-sheet surface along the fibril axis. The lack of two methyl groups and charge in PIB not only improves its hydrophobicity but also leads to a deeper insertion of PIB compared to ThT into the surface grooves. This significantly increases the steric, aromatic, and hydrophobic interactions, and hence leads to stronger binding. Simulations on protofibrils consisting of the more-toxic Aβ_17-42 revealed an additional binding mode in which PIB and ThT insert into the channel that forms in the loop region of the protofibril, sandwiched between two sheet layers. Our simulations indicate that the rotation between the two ring parts of the dyes is significantly more restricted when the dyes are bound to the surface of the cross-β-subunits or to the channel inside the Aβ_17-42 cross-β-subunit, compared with free solution. The specific conformations of the dyes are influenced by small chemical modifications (ThT versus PIB) and by the environment in which the dye is placed.     

Binding of congo red to amyloid protofibrils of the Alzheimer aβ(9-40) Peptide probed by molecular dynamics simulations

Congo red (CR) is a commonly used histological amyloid dye and a weak amyloid inhibitor. There is currently no experimentally available structure of CR bound to an amyloid fibril and the binding modes, and the mechanisms governing its inhibitory and optical properties are poorly understood. In this work, we present the first, to our knowledge, atomistically detailed picture of CR binding to protofibrils of the Alzheimer Aβ_9–40 peptide. We identify three major binding modes, with the primary mode residing in the grooves formed by the β-sheets, and observe a restriction of the torsional rotation of the CR molecule upon binding. Our simulations reveal a novel, to our knowledge, electrostatic steering mechanism that plays an important role in the initial recognition and binding of CR to the positively charged surface residues of the fibril. Our simulations provide new, to our knowledge, insights into the striking spectrophotometric and inhibitory properties of CR. In particular, we show that birefringence upon CR binding is due to the anisotropic orientation of the CR dipoles resulting from the spatial ordering of these molecules in the grooves along the fibril axis. The fluorescent enhancement of the bound CR, in turn, is associated with the torsional restriction of this molecule upon binding.     

pH-dependence of the specific binding of Cu(II) and Zn(II) ions to the amyloid-β peptide

Metal ions like Cu(II) and Zn(II) are accumulated in Alzheimer's disease amyloid plaques. The amyloid-β (Aβ) peptide involved in the disease interacts with these metal ions at neutral pH via ligands provided by the N-terminal histidines and the N-terminus. The present study uses high-resolution NMR spectroscopy to monitor the residue-specific interactions of Cu(II) and Zn(II) with (15)N- and (13)C,(15)N-labeled Aβ(1-40) peptides at varying pH levels. At pH 7.4 both ions bind to the specific ligands, competing with one another. At pH 5.5 Cu(II) retains its specific histidine ligands, while Zn(II) seems to lack residue-specific interactions. The low pH mimics acidosis which is linked to inflammatory processes in vivo. The results suggest that the cell toxic effects of redox active Cu(II) binding to Aβ may be reversed by the protective activity of non-redox active Zn(II) binding to the same major binding site under non-acidic conditions. Under acidic conditions, the protective effect of Zn(II) may be decreased or changed, since Zn(II) is less able to compete with Cu(II) for the specific binding site on the Aβ peptide under these conditions.     

Soluble aggregates of the amyloid-β peptide are trapped by serum albumin to enhance amyloid-β activation of endothelial cells

Background  

Self-assembly of the amyloid-β peptide (Aβ) has been implicated in the pathogenesis of Alzheimer's disease (AD). As a result, synthetic molecules capable of inhibiting Aβ self-assembly could serve as therapeutic agents and endogenous molecules that modulate Aβ self-assembly may influence disease progression. However, increasing evidence implicating a principal pathogenic role for small soluble Aβ aggregates warns that inhibition at intermediate stages of Aβ self-assembly may prove detrimental. Here, we explore the inhibition of Aβ_1–40 self-assembly by serum albumin, the most abundant plasma protein, and the influence of this inhibition on Aβ_1–40 activation of endothelial cells for monocyte adhesion.     

Binding of histamine to the H1 receptor—a molecular dynamics study

Binding of histamine to the G-protein coupled histamine H₁ receptor plays an important role in the context of allergic reactions; however, no crystal structure of the resulting complex is available yet. To deduce the histamine binding site, we performed unbiased molecular dynamics (MD) simulations on a microsecond time scale, which allowed to monitor one binding event, in which particularly the residues of the extracellular loop 2 were involved in the initial recognition process. The final histamine binding pose in the orthosteric pocket is characterized by interactions with Asp107^3.32, Tyr108^3.33, Thr194^5.43, Asn198^5.46, Trp428^6.48, Tyr431^6.51, Phe432^6.52, and Phe435^6.55, which is in agreement with existing mutational data. The conformational stability of the obtained complex structure was subsequently confirmed in 2 μs equilibrium MD simulations, and a metadynamics simulation proved that the detected binding site represents an energy minimum. A complementary investigation of a D107A mutant, which has experimentally been shown to abolish ligand binding, revealed that this exchange results in a significantly weaker interaction and enhanced ligand dynamics. This finding underlines the importance of the electrostatic interaction between the histamine ammonium group and the side chain of Asp107^3.32 for histamine binding.     

Cryogenic solid state NMR studies of fibrils of the Alzheimer’s disease amyloid-β peptide: perspectives for DNP

Dynamic Nuclear Polarization solid-state NMR holds the potential to enable a dramatic increase in sensitivity by exploiting the large magnetic moment of the electron. However, applications to biological solids are hampered in uniformly isotopically enriched biomacromolecules due to line broadening which yields a limited spectral resolution at cryogenic temperatures. We show here that high magnetic fields allow to overcome the broadening of resonance lines often experienced at liquid nitrogen temperatures. For a fibril sample of the Alzheimer’s disease β-amyloid peptide, we find similar line widths at low temperature and at room temperature. The presented results open new perspectives for structural investigations in the solid-state.     

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

Abstract

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

Communicated by Ramaswamy H. Sarma     

Critical test of bead–spring model to resolve the scaling laws of polymer melts: a molecular dynamics study

Abstract

To examine the intrinsic nature of the bead–spring Kremer–Grest (KG) model, long-time molecular dynamics simulations are performed. Certain scaling laws for representative polymer properties are compared with theoretical predictions. The results for static properties satisfy the expected static Gaussian nature, irrespective of the chain length. In contrast, the results for the dynamic properties of short chains show a clear discrepancy from theoretical predictions that assume ideal chain motion. This is clear evidence that the Gaussian nature of the dynamics of short chains is not necessarily established for the actual KG model, despite it being designed to have Gaussian characteristics by virtue of its stochastic equations of motion. This intrinsic nature of the KG model should be considered carefully when using this model for applications that involve relatively short chains.     

Specific binding of a β-cyclodextrin dimer to the amyloid β peptide modulates the peptide aggregation process

Alzheimer's disease involves progressive neuronal loss. Linked to the disease is the amyloid β (Aβ) peptide, a 38-43-amino acid peptide found in extracellular amyloid plaques in the brain. Cyclodextrins are nontoxic, cone-shaped oligosaccharides with a hydrophilic exterior and a hydrophobic cavity making them suitable hosts for aromatic guest molecules in water. β-Cyclodextrin consists of seven α-d-glucopyranoside units and has been shown to reduce the level of fibrillation and neurotoxicity of Aβ. We have studied the interaction between Aβ and a β-cyclodextrin dimer, consisting of two β-cyclodextrin monomers connected by a flexible linker. The β-cyclodextrin monomer has been found to interact with Aβ(1-40) at sites Y10, F19, and/or F20 with a dissociation constant (K(D)) of 3.9 ± 2.0 mM. Here (1)H-(15)N and (1)H-(13)C heteronuclear single-quantum correlation nuclear magnetic resonance (NMR) spectra show that in addition, the β-cyclodextrin monomer and dimer bind to the histidines. NMR translational diffusion experiments reveal the increased affinity of the β-cyclodextrin dimer (apparent K(D) of 1.1 ± 0.5 mM) for Aβ(1-40) compared to that of the β-cyclodextrin monomer. Kinetic aggregation experiments based on thioflavin T fluorescence indicate that the dimer at 0.05-5 mM decreases the lag time of Aβ aggregation, while a concentration of 10 mM increases the lag time. The β-cyclodextrin monomer at a high concentration decreases the lag time of the aggregation. We conclude that cyclodextrin monomers and dimers have specific, modulating effects on the Aβ(1-40) aggregation process. Transmission electron microscopy shows that the regular fibrillar aggregates formed by Aβ(1-40) alone are replaced by a major fraction of amorphous aggregates in the presence of the β-cyclodextrin dimer.     

Modulation of aggregate size- and shape-distributions of the amyloid-β peptide by a designed β-sheet breaker

A peptide with 42 amino acid residues (Aβ42) plays a key role in the pathogenesis of the Alzheimer’s disease. It is highly prone to self aggregation leading to the formation of fibrils which are deposited in amyloid plaques in the brain of diseased individuals. In our study we established a method to analyze the aggregation behavior of the Aβ peptide with a combination of sedimentation velocity centrifugation and enhanced data evaluation software as implemented in the software package UltraScan. Important information which becomes accessible by this methodology is the s-value distribution and concomitantly also the shape-distribution of the Aβ peptide aggregates generated by self-association. With this method we characterized the aggregation modifying effect of a designed β-sheet breaker molecule. This compound is built from three head-to-tail connected aminopyrazole moieties and represents a derivative of the already described Tripyrazole. By addition of this compound to a solution of the Aβ42 peptide the maximum of the s-value distribution was clearly shifted to smaller s-values as compared to solutions where only the vehicle DMSO was added. This shift to smaller s-values was stable for at least 7 days. The information about size- and shape-distributions present in aggregated Aβ42 solutions was confirmed by transmission electron microscopy and by measurement of amyloid formation by thioflavin T fluorescence.     

Acetylation and phosphorylation processes modulate Tau’s binding to microtubules: A molecular dynamics study

The microtubule-associated protein Tau has its normal function impaired when undergoing post-translational modifications. In this work, molecular modelling techniques were used to infer the effects of acetylation and phosphorylation in Tau's overall conformation, electrostatics, and interactions, but mostly in Tau's ability to bind microtubules. Reported harmful Lys sites were mutated by its acetylated form, generating eight different acetylated Tau (aTau) analogues. Similarly, phosphorylation sites found in normal brains and in Alzheimer’s lesioned brains were considered to design phosphorylated Tau (pTau) analogues. All these designed variants were evaluated in intracellular fluid and near a microtubule (MT) model. Our in silico findings demonstrated that the electrostatic changes, due to the absence of positive Lys’ charges in acetylation cases, or the increasingly negative charge in the phosphorylated forms, hamper the association to the MT tubulins in most cases. Post-translational modifications also pose very distinct conformations to the ones described for native Tau, which hinders the microtubule-binding region (MTBR) and turns difficult the expected binding. Our study elucidates important molecular processes behind Tau abnormal function which can inspire novel therapeutics to address Alzheimer’s disease.     

Insight into the stability of cross-β amyloid fibril from VEALYL short peptide with molecular dynamics simulation

Amyloid fibrils are found in many fatal neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, type II diabetes, and prion disease. The VEALYL short peptide from insulin has been confirmed to aggregate amyloid-like fibrils. However, the aggregation mechanism of amyloid fibril is poorly understood. Here, we utilized molecular dynamics simulation to analyse the stability of VEALYL hexamer. The statistical results indicate that hydrophobic residues play key roles in stabilizing VEALYL hexamer. Single point and two linkage mutants confirmed that Val1, Leu4, and Tyr5 of VEALYL are key residues. The consistency of the results for the VEALYL oligomer suggests that the intermediate states might be trimer (3-0) and pentamer(3-2). These results can help us to obtain an insight into the aggregation mechanism of amyloid fibril. These methods can be used to study the stability of amyloid fibril from other short peptides.     

Binding studies of truncated variants of the Aβ peptide to the V-domain of the RAGE receptor reveal Aβ residues responsible for binding

Alzheimer's disease (AD) symptoms correlate with the concentration of soluble, although not necessarily monomeric forms of Aβ peptide in the brain parenchyma. The RAGE receptor has been implicated as the protein responsible for active transport of Aβ from blood circulation to the brain. In murine models of AD, inhibition of the Aβ:RAGE interaction decreases the levels of Aβ in the brain. Inhibition of the Aβ:RAGE interaction would be a promising alternative for the therapy of AD. Rational design of an Aβ:RAGE interaction blocker requires detailed knowledge of the structure of the complex. However, the binding domain of RAGE is natively unfolded in physiological conditions, which severely hampers the application of classic methods of protein structure analysis to the design of an antagonist. Here, alternative methods are used to characterize the structural properties of the RAGE-ligand binding domain and to monitor the binding of a series of truncated variants of Aβ. Using intrinsic RAGE tryptophan fluorescence and mass spectrometry of non-covalent protein-ligand complexes we have identified shorter versions of Aβ that bind to the RAGE V-domain. We have also shown in cell culture experiments that a selected shortened version of Aβ effectively inhibits full-length Aβ, RAGE-mediated, cell uptake. Thus, a truncated version of Aβ capable of blocking its receptor-mediated internalization was established, revealing the binding code and providing the lead compound in the process of drug design.