Structures and free energy landscapes of aqueous zinc(II)-bound amyloid-β(1-40) and zinc(II)-bound amyloid-β(1-42) with dynamics

Molecular Insights into Cu/Zn Metal Response to the Amyloid β-Peptide (1-42)

Aβ1-40 peptide and Aβ1-42 peptide are the building units of beta-amyloid plaques present in Alzheimer's disease (AD)-affected brain. The binding affinity of various divalent metal ions such as Cu and Zn present in AD-affected brain with different amino acids available in Aβ-peptide became the focus to explore their role in soluble neurotoxic oligomer formation. Cu2+ metal ions are known to enhance the neurotoxicity of the Aβ1-42 peptide by catalyzing the formation of soluble neurotoxic oligomers. The competitive preference of both Cu2+ and Zn2+ simultaneously to interact with the Aβ-peptide is unknown. The divalent Cu and Zn ions were inserted in explicit aqueous Aβ1-42 peptide configurations to get insights into the binding competence of these metal ions with peptides using classical molecular dynamics (MD) simulations. The metal-ion interactions reveal that competitive binding preferences of various peptide sites become metal-ion-specific and differ significantly. For Cu2+, interactions are found to be more significant with respect to those of Asp-7, His-6, Glu-11, and His-14. Asp-1, Glu-3, Asp-7, His-6, Glu-11, and His-13 amino acid residues show higher affinity toward Zn2+ ions. MD simulations show notable variation in the solvent-accessible surface area in the hydrophobic region of the peptide. Infinitesimal mobility was obtained for Zn2+ compared to Cu2+ in an aqueous solution and Cu2+ diffusivity deviated significantly at different time scales, proving its labile features in aqueous Aβ1-42 peptides.

Atomistic molecular simulations of Aβ-Zn conformational ensembles

The amyloid-forming Aβ peptide is able to interact with metal cations to form very stable complexes that influence fibril formation and contribute to the onset of Alzheimer's disease. Multiple structures of peptides derived from Aβ in complex with different metals have been resolved experimentally to provide an atomic-level description of the metal-protein interactions. However, Aβ is intrinsically disordered, and hence more amenable to an ensemble description. Molecular dynamics simulations can now reach the timescales needed to generate ensembles for these type of complexes. However, this requires accurate force fields both for the protein and the protein-metal interactions. Here we use state-of-the-art methods to generate force field parameters for the Zn(II) cations in a set of complexes formed by different Aβ variants and combine them with the Amber99SB*-ILDN optimized force field. Upon comparison of NMR experiments with the simulation results, further optimized with a Bayesian/Maximum entropy approach, we provide an accurate description of the molecular ensembles for most Aβ-metal complexes. We find that the resulting conformational ensembles are more heterogeneous than the NMR models deposited in the Protein Data Bank.

Rules governing metal coordination in Aβ-Zn(II) complex models from quantum mechanical calculations

Transition metals directly contribute to the neurotoxicity of the aggregates of the amyloid-forming Aβ peptide. The understanding and rationalization of the coordination modes of metals to Aβ amyloid is, therefore, of paramount importance to understand the capacity of a given metal to promote peptide aggregation. Experimentally, multiple Aβ-metal structures have been resolved, which exhibit different modes of coordination in both the monomeric and oligomeric forms of Aβ. Although Zn(II) metalloproteins are very abundant and often involve cysteine residues in the first coordination shell, in the case of Aβ-Zn(II), though, Zn(II) is coordinated by glutamic/aspartic acid and/or histidine residues exclusively, making for an interesting case study. Here we present a systematic analysis of the underlying chemistry on Aβ-Zn(II) coordination, where relative stabilities of different coordination arrangements indicate that a mixture of Glu/Asp and His residues is favored. A detailed comparison between different coordination shell geometries shows that tetrahedral coordination is generally favored in the aqueous phase. Our calculations show an interplay between dative covalent interactions and electrostatics which explains the observed trends. Multiple structures deposited in the Protein Data Bank support our findings, suggesting that the trends found in our work may be transferable to other Zn(II) metalloproteins with this type of coordination.

Concentration Effect, Structural Properties, and Driving Force on Aβ28 Dimerization with and without Zn2+ Cooperation: Learning from Replica Exchange Sampling

Zn²⁺ is a very important factor in promoting the formation of amyloid beta (Aβ) aggregates and amyloid plaques. The Zn²⁺ -bound Aβ species generate amorphous or low molecular-weight oligomers. However, it is a lack of studies to approach the starting structural features (dimerization) in Aβ nucleation processes with and without Zn²⁺ , which is the key point in understanding Zn²⁺ -induced nucleation mechanisms. To better understand the effect of concentration, structural properties, and the driving force, 14 independent replica exchange molecular dynamics simulations were performed in Aβ₂₈ dimerization with and without Zn²⁺ (zAβ₂₈ ) cooperation. Our scanning results show that the aggregation propensity is easier in Aβ₂₈ -Aβ₂₈ and Aβ₂₈ -zAβ₂₈ systems than zAβ₂₈ -zAβ₂₈ system. In binding property, the Aβ₂₈ -Aβ₂₈ model (-61.5 kcal mol^-1 ) is stronger than zAβ₂₈ -zAβ₂₈ (-26.6 kcal mol^-1 ) and Aβ₂₈ -zAβ₂₈ (-7.24 kcal mol^-1 ) models. Further analysis confirmed that H13 and H14 residues play specific roles in the three systems. The key point is the orientation of N atom of the imidazole ring in histidine residues. Furthermore, we discovered different driving forces for each system. Our current study contributes to the understanding of how the Aβ₂₈ dimer interacts with Zn²⁺ , which could lead to new insights into Zn²⁺ -induced nucleation mechanisms.

Zinc Induced Aβ16 Aggregation Modeled by Molecular Dynamics

It is widely accepted that the addition of zinc leads to the formation of neurotoxic nonfibrillar aggregates of beta-amyloid peptides Aβ₄₀ and Aβ₄₂ and at the same time destabilizes amyloid fibrils. However, the mechanism of the effect of zinc on beta-amyloid is not fully understood. In this study, a fast zinc-induced aggregation of Aβ₁₆ (as compared to a system without zinc) via the formation of Aβ₁₆ dimers with one zinc ion coordinated in the metal-binding site ₁₁EVHH₁₄, followed by their polymerization, has been studied by molecular dynamics. The best aggregation was shown by the system composed of Aβ₁₆ dimers bound by one zinc ion, with no additional zinc in solution. The presence of Aβ₁₆ dimers was a major condition, sufficient for fast aggregation into larger complexes. It has been shown that the addition of zinc to a system with already formed dimers does not substantially affect the characteristics and rate of aggregation. At the same time, an excessive concentration of zinc at the early stages of the formation of conglomerates can negatively affect aggregation, since in systems where zinc ions occupied the ₁₁EVHH₁₄ coordination center and the His6 residue of every Aβ₁₆ monomer, the aggregation proceeded more slowly and the resulting complexes were not as large as in the zinc-free Aβ system. Thus, this study has shown that the formation of Aβ₁₆ dimers bound through zinc ions at the ₁₁EVHH₁₄ sites of the peptides plays an important role in the formation of neurotoxic non-fibrillar aggregates of beta-amyloid peptide Aβ₁₆. The best energetically favorable structure has been obtained for the complex of two Aβ₁₆ dimers with two zinc ions.

Defining the Neuropathological Aggresome across in Silico, in Vitro, and ex Vivo Experiments

The loss of proteostasis over the life course is associated with a wide range of debilitating degenerative diseases and is a central hallmark of human aging. When left unchecked, proteins that are intrinsically disordered can pathologically aggregate into highly ordered fibrils, plaques, and tangles (termed amyloids), which are associated with countless disorders such as Alzheimer's disease, Parkinson's disease, type II diabetes, cancer, and even certain viral infections. However, despite significant advances in protein folding and solution biophysics techniques, determining the molecular cause of these conditions in humans has remained elusive. This has been due, in part, to recent discoveries showing that soluble protein oligomers, not insoluble fibrils or plaques, drive the majority of pathological processes. This has subsequently led researchers to focus instead on heterogeneous and often promiscuous protein oligomers. Unfortunately, significant gaps remain in how to prepare, model, experimentally corroborate, and extract amyloid oligomers relevant to human disease in a systematic manner. This Review will report on each of these techniques and their successes and shortcomings in an attempt to standardize comparisons between protein oligomers across disciplines, especially in the context of neurodegeneration. By standardizing multiple techniques and identifying their common overlap, a clearer picture of the soluble neuropathological aggresome can be constructed and used as a baseline for studying human disease and aging.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

Transition metal complexes as imaging or therapeutic agents for neurodegenerative diseases

Neurodegenerative diseases are the result of neurodegeneration, which is the process of losing neuronal functions gradually due to the irreversible damage and death of neurons. Metal complexes have attracted intense interest over recent decades as probes or inhibitors of biomolecules.

Identification of a Novel Multifunctional Ligand for Simultaneous Inhibition of Amyloid-Beta (Aβ42) and Chelation of Zinc Metal Ion

Zinc binding to β-amyloid structure could promote amyloid-β aggregation, as well as reactive oxygen species (ROS) production, as suggested in many experimental and theoretical studies. Therefore, the introduction of multifunctional drugs capable of chelating zinc metal ion and inhibiting Aβ aggregation is a promising strategy in the development of AD treatment. The present study has evaluated the efficacy of a new bifunctional peptide drug using molecular docking and molecular dynamics (MD) simulations. This drug comprises two different domains, an inhibitor domain, obtained from the C-terminal hydrophobic region of Aβ, and a Zn²⁺ chelating domain, derived from rapeseed meal, merge with a linker. The multifunctionality of the ligand was evaluated using a comprehensive set of MD simulations spanning up to 3.2 μs including Aβ relaxation, ligand-Zn²⁺ bilateral interaction, and, more importantly, ligand-Zn²⁺-Aβ₄₂ trilateral interactions. Analysis of the results strongly indicated that the bifunctional ligand can chelate zinc metal ion and avoid Aβ aggregation simultaneously. The present study illustrated that the proposed ligand has considerable hydrophobic interactions and hydrogen bonding with monomeric Aβ in the presence of zinc metal ion. Therefore, in light of these considerable interactions and contacts, the α-helical structure of Aβ has been enhanced, while the β-sheet formation is prevented and the α-helix native structure is protected. Furthermore, the analysis of interactions between Aβ and ligand-zinc complex revealed that the zinc metal ion is coordinated to Met13, the ending residue of the ligand and merely one residue in Aβ. The results have proven the previous experimental and theoretical findings in the literature about Aβ interactions with zinc metal ion and also Aβ interactions with the first domain of the proposed ligand. Moreover, the current research has evaluated the chelation using MD simulation and linear interaction energy (LIE) methods, and the result has been satisfactorily verified with previous experimental and theoretical (DFT) studies.

Intrinsically disordered proteins in various hypotheses on the pathogenesis of Alzheimer's and Parkinson's diseases

Amyloid-β (Aβ) and α-synuclein (αS) are two intrinsically disordered proteins (IDPs) at the centers of the pathogenesis of Alzheimer's and Parkinson's diseases, respectively. Different hypotheses have been proposed for explanation of the molecular mechanisms of the pathogenesis of these two diseases, with these two IDPs being involved in many of these hypotheses. Currently, we do not know, which of these hypothesis is more accurate. Experiments face challenges due to the rapid conformational changes, fast aggregation processes, solvent and paramagnetic effects in studying these two IDPs in detail. Furthermore, pathological modifications impact their structures and energetics. Theoretical studies using computational chemistry and computational biology have been utilized to understand the structures and energetics of Aβ and αS. In this chapter, we introduce Aβ and αS in light of various hypotheses, and discuss different experimental and theoretical techniques that are used to study these two proteins along with their weaknesses and strengths. We suggest that a promising solution for studying Aβ and αS at the center of varying hypotheses could be provided by developing new techniques that link quantum mechanics, statistical mechanics, thermodynamics, bioinformatics to machine learning. Such new developments could also lead to development in experimental techniques.

Effects of Metal Ions on Aβ42 Peptide Conformations from Molecular Simulation Studies

In this study, we investigate the conformational characteristics of full-length Aβ₄₂ peptide monomers in the presence of Na⁺ and Zn²⁺ metal ions using atomistic molecular dynamics (MD) simulations with an aim to explore the possible driving forces behind enhanced aggregation rates of the peptides in the presence of salts. The calculations reveal that the presence of metal ions shifts the conformational equilibrium more toward the compact ordered Aβ structures. Such compact ordered structures stabilized by distant nonlocal contacts between two crucial hydrophobic segments, hp1 and hp2, primarily through two important hydrophobic aromatic residues, Phe-19 and Phe-20, are expected to trigger the aggregation process at a faster rate by populating and stabilizing the aggregation prone structures. Formation of a significant number of such distant contacts in the presence of Na⁺ ions has also been found to result in breaking of the N-terminal helix. On the contrary, binding of Zn²⁺ ion to Aβ peptide is highly specific, which stabilizes the N-terminal helix instead of breaking it. This explains why the aggregation rate of Aβ peptides is higher in the presence of divalent Zn²⁺ ions than monovalent Na⁺ ions. Relatively higher overall stability of the most populated Aβ peptide monomers in the presence of Zn²⁺ ions has been found to be associated with specific Zn²⁺-Aβ binding and significant free energy gain.

Transition Metal Ion Interactions with Disordered Amyloid-β Peptides in the Pathogenesis of Alzheimer's Disease: Insights from Computational Chemistry Studies

Monomers and oligomers of the amyloid-β peptide aggregate to form the fibrils found in the brains of Alzheimer's disease patients. These monomers and oligomers are largely disordered and can interact with transition metal ions, affecting the mechanism and kinetics of amyloid-β aggregation. Due to the disordered nature of amyloid-β, its rapid aggregation, as well as solvent and paramagnetic effects, experimental studies face challenges in the characterization of transition metal ions bound to amyloid-β monomers and oligomers. The details of the coordination chemistry between transition metals and amyloid-β obtained from experiments remain debated. Furthermore, the impact of transition metal ion binding on the monomeric or oligomeric amyloid-β structures and dynamics are still poorly understood. Computational chemistry studies can serve as an important complement to experimental studies and can provide additional knowledge on the binding between amyloid-β and transition metal ions. Many research groups conducted first-principles calculations, ab initio molecular dynamics simulations, quantum mechanics/classical mechanics simulations, and classical molecular dynamics simulations for studying the interplay between transition metal ions and amyloid-β monomers and oligomers. This review summarizes the current understanding of transition metal interactions with amyloid-β obtained from computational chemistry studies. We also emphasize the current view of the coordination chemistry between transition metal ions and amyloid-β. This information represents an important foundation for future metal ion chelator and drug design studies aiming to combat Alzheimer's disease.

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

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

Insights into the Molecular Mechanisms of Alzheimer's and Parkinson's Diseases with Molecular Simulations: Understanding the Roles of Artificial and Pathological Missense Mutations in Intrinsically Disordered Proteins Related to Pathology

Amyloid-β and α-synuclein are intrinsically disordered proteins (IDPs), which are at the center of Alzheimer's and Parkinson's disease pathologies, respectively. These IDPs are extremely flexible and do not adopt stable structures. Furthermore, both amyloid-β and α-synuclein can form toxic oligomers, amyloid fibrils and other type of aggregates in Alzheimer's and Parkinson's diseases. Experimentalists face challenges in investigating the structures and thermodynamic properties of these IDPs in their monomeric and oligomeric forms due to the rapid conformational changes, fast aggregation processes and strong solvent effects. Classical molecular dynamics simulations complement experiments and provide structural information at the atomic level with dynamics without facing the same experimental limitations. Artificial missense mutations are employed experimentally and computationally for providing insights into the structure-function relationships of amyloid-β and α-synuclein in relation to the pathologies of Alzheimer's and Parkinson's diseases. Furthermore, there are several natural genetic variations that play a role in the pathogenesis of familial cases of Alzheimer's and Parkinson's diseases, which are related to specific genetic defects inherited in dominant or recessive patterns. The present review summarizes the current understanding of monomeric and oligomeric forms of amyloid-β and α-synuclein, as well as the impacts of artificial and pathological missense mutations on the structural ensembles of these IDPs using molecular dynamics simulations. We also emphasize the recent investigations on residual secondary structure formation in dynamic conformational ensembles of amyloid-β and α-synuclein, such as β-structure linked to the oligomerization and fibrillation mechanisms related to the pathologies of Alzheimer's and Parkinson's diseases. This information represents an important foundation for the successful and efficient drug design studies.

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.

TMPyP Inhibits Amyloid-β Aggregation and Alleviates Amyloid-Induced Cytotoxicity

The aggregation or misfolding of amyloid-β (Aβ) is a major pathological hallmark of Alzheimer's disease (AD). The regulation of Aβ aggregation is thought to be an effective strategy for AD treatment. The capability of a water-soluble porphyrin, 5,10,15,20-tetrakis(N-methyl-4-pyridyl)porphyrin (TMPyP), to inhibit Aβ aggregation and to lower Aβ-induced toxicity was demonstrated. As evidenced by surface plasmon resonance and circular dichroism, TMPyP can not only disrupt Aβ aggregation but also disassemble the preformed Aβ aggregates. The atomic force microscopy imaging proves that TMPyP inhibits the formation of both oligomers and fibrils. Molecular dynamic simulations provide an insight into the interaction between TMPyP and Aβ at the molecular level. The half-maximal inhibitory concentrations of TMPyP acting on the oligomers and fibrils were determined to be 0.6 and 0.43 μM, respectively. As a member of porphyrin family, TMPyP is of rather low cytotoxicity, and the cytotoxicity of the Aβ aggregates was also relieved upon coincubation with TMPyP. The excellent performance of TMPyP thus makes it a potential drug candidate for AD therapy.

Molecular Dynamics Study on the Inhibition Mechanisms of Drugs CQ1-3 for Alzheimer Amyloid-β40 Aggregation Induced by Cu(2.)

The aggregation of amyloid-β (Aβ) peptide induced by Cu(2+) is a key factor in development of Alzheimer's disease (AD), and metal ion chelation therapy enables treatment of AD. Three CQi (i = 1, 2, and 3 with R = H, Cl, and NO2, respectively) drugs had been verified experimentally to be much stronger inhibitors than the pioneer clioquinol (CQ) in both disaggregation of Aβ40 aggregate and reduction of toxicity induced by Cu(2+) binding at low pH. Due to the multiple morphologies of Cu(2+)-Aβ40 complexes produced at different pH states, we performed a series of molecular dynamics simulations to explain the structural changes and morphology characteristics as well as intrinsic disaggregation mechanisms of three Cu(2+)-Aβ40 models in the presence of any of the three CQi drugs at both low and high pH states. Three inhibition mechanisms for CQi were proposed as "insertion", "semi-insertion", and "surface" mechanisms, based on the morphologies of CQi-model x (CQi-x, x = 1, 2, and 3) and the strengths of binding between CQi and the corresponding model x. The insertion mechanism was characterized by the morphology with binding strength of more than 100 kJ/mol and by CQi being inserted or embedded into the hydrophobic cavity of model x. In those CQi-x morphologies with lower binding strength, CQi only attaches on the surface or inserts partly into Aβ peptide. Given the evidence that the binding strength is correlated positively with the effectiveness of drug to inhibit Aβ aggregation and thus to reduce toxicity, the data of binding strength presented here can provide a reference for one to screen drugs. From the point of view of binding strength, CQ2 is the best drug. Because of the special role of Asp23 in both Aβ aggregation and stabilizing the Aβ fibril, the generation of a H-bond between CQ3 and Asp23 of the Aβ40 peptide is believed to be responsible for CQ3 having the strongest disaggregation capacity. Therefore, besides strong binding, stronger propensity to H-bond with Asp23 would be another key factor to be taken seriously into account in drug screens. Meanwhile, the structural characteristics of drug CQi itself are also worthy of attention. First, the increasing polarity from CQ1 and CQ2 to CQ3 in turn results in increasing probability and strength of the interaction between the drug and the N-terminal (NT) region of Aβ40, which obviously inhibits Aβ peptide aggregation induced by Cu(2+) binding. Second, both the benzothiazole ring and phenol ring of CQi can overcome the activation energy barrier (∼16 kJ/mol) to rotate flexibly around the intramolecular C7-N14 bond to achieve the maximum match and interaction with the ambient Aβ40 residues. Such a structural feature of CQi paves the new way for ones in selection and modification of a drug.

Different zinc(II) complex species and binding modes at Aβ N-terminus drive distinct long range cross-talks in the Aβ monomers

The present study addresses the reconstruction of the free-energy landscapes of amyloid-beta1-42 (Aβ42) coordinated respectively with one and two zinc ions, to scrutinize whether different Aβ-zinc complex species, i.e., mononuclear and dinuclear metal complexes, induce different Aβ conformation features. We found a subtle switch of intramolecular interactions, depending both on the zinc coordination environment and on the peptide to zinc stoichiometric ratio. On the one side, hairpin-like structures are predominant in mononuclear complexes, where a salt-bridge that involves Lys28-Glu22 and Lys16-Asp23 is stabilized. On the other side, elongated conformations are instead stabilized in the dinuclear zinc complexes. Experimental studies of atomic force microscopy as well as of zinc-Aβ complex species distribution diagrams provide evidence that the theoretical calculations can be rationalized in terms of the correlation between the increased amount of amorphous aggregates and the Aβ/Zn(2+) ratio.

Coupling of Zinc-Binding and Secondary Structure in Nonfibrillar Aβ40 Peptide Oligomerization

Nonfibrillar neurotoxic amyloid β (Aβ) oligomer structures are typically rich in β-sheets, which could be promoted by metal ions like Zn(2+). Here, using molecular dynamics (MD) simulations, we systematically examined combinations of Aβ40 peptide conformations and Zn(2+) binding modes to probe the effects of secondary structure on Aβ dimerization energies and kinetics. We found that random conformations do not contribute to dimerization either thermodynamically or kinetically. Zn(2+) couples with preformed secondary structures (α-helix and β-hairpin) to speed dimerization and stabilize the resulting dimer. Partial α-helices increase the dimerization speed, and dimers with α-helix rich conformations have the lowest energy. When Zn(2+) coordinates with residues D1, H6, H13, and H14, Aβ40 β-hairpin monomers have the fastest dimerization speed. Dimers with experimentally observed zinc coordination (E11, H6, H13, and H14) form with slower rate but have lower energy. Zn(2+) cannot stabilize fibril-like β-arch dimers. However, Zn(2+)-bound β-arch tetramers have the lowest energy. Collectively, zinc-stabilized β-hairpin oligomers could be important in the nucleation-polymerization of cross-β structures. Our results are consistent with experimental findings that α-helix to β-structural transition should accompany Aβ aggregation in the presence of zinc ions and that Zn(2+) stabilizes nonfibrillar Aβ oligomers and, thus, inhibits formation of less toxic Aβ fibrils.

Aβ42 and Aβ40: similarities and differences

The abnormal accumulation of amyloid-β (Aβ) peptide in the brain is one of the most important hallmarks of Alzheimer's disease. Aβ is an aggregation-prone and toxic polypeptide with 39-43 residues, derived from the amyloid precursor protein proteolysis process. According to the amyloid hypothesis, abnormal accumulation of Aβ in the brain is the primary influence driving Alzheimer's disease pathologies. Among all kinds of Aβ isoforms, Aβ40 and Aβ42 are believed to be the most important ones. Although these two kinds of Aβ differ only in two amino acid residues, recent studies show that they differ significantly in their metabolism, physiological functions, toxicities, and aggregation mechanism. In this review, we mainly summarize the similarities and differences between Aβ42 and Aβ40, recent studies on selective inhibitors as well as probes will also be mentioned.

Arginine and disordered amyloid-β peptide structures: molecular level insights into the toxicity in Alzheimer's disease

Recent studies present that the single arginine (R) residue in the sequence of Aβ42 adopts abundant β-sheet structure and forms stable salt bridges with various residues. Furthermore, experiments proposed that R stimulates the Aβ assembly and arginine (R) to alanine (A) mutation (R5A) decreases both aggregate formation tendency and the degree of its toxicity. However, the exact roles of R and R5A mutation in the structures of Aβ42 are poorly understood. Extensive molecular dynamics simulations along with thermodynamic calculations present that R5A mutation impacts the structures and free energy landscapes of the aqueous Aβ42 peptide. The β-sheet structure almost disappears in the Ala21-Ala30 region but is more abundant in parts of the central hydrophobic core and C-terminal regions of Aβ42 upon R5A mutation. More abundant α-helix is adopted in parts of the N-terminal and mid-domain regions and less prominent α-helix formation occurs in the central hydrophobic core region of Aβ42 upon R5A mutation. Interestingly, intramolecular interactions between N- and C-terminal or mid-domain regions disappear upon R5A mutation. The structures of Aβ42 are thermodynamically less stable and retain reduced compactness upon R5A mutation. R5A mutant-type structure stability increases with more prominent central hydrophobic core and mid-domain or C-terminal region interactions. Based on our results reported in this work, small organic molecules and antibodies that avoid β-sheet formation in the Ala21-Ala30 region and hinder the intramolecular interactions occurring between the N-terminal and mid-domain or C-terminal regions of Aβ42 may help to reduce Aβ42 toxicity in Alzheimer's disease.

Effects of Zn2+ binding on the structural and dynamic properties of amyloid β peptide associated with Alzheimer's disease: Asp1 or Glu11?

Extensive experimental and computational studies have suggested that multiple Zn(2+) binding modes in amyloid β (Aβ) peptides could exist simultaneously. However, consistent results have not been obtained for the effects of Zn(2+) binding on Aβ structure, dynamics, and kinetics in particular. Some key questions such as why it is so difficult to distinguish the polymorphic states of metal ions binding to Aβ and what the underlying rationale is, necessitate elucidation. In this work, two 3N1O Zn(2+) binding modes were constructed with three histidines (His(6), His(13), and His(14)), and Asp(1)/Glu(11) of Aβ40 coordinated to Zn(2+). Results from molecular dynamics simulations reveal that the conformational ensembles of different Zn(2+)-Aβ40 complexes are nonoverlapping. The formation of turn structure and, especially, the salt bridge between Glu(22)/Asp(23) and Lys(28) is dependent on specific Zn(2+) binding mode. Agreement with available NMR observations of secondary and tertiary structures could be better achieved if the two simulation results are considered together. The free energy landscape constructed by combining both conformations of Aβ40 indicates that transitions between distinct Aβ40 conformations thar are ready for Zn(2+) binding could be possible in aqueous solution. Markov state model analyses reveal the complex network of conformational space of Aβ40 modeulated by Zn(2+) binding, suggesting various misfolding pathways. The binding free energies evaluated using a combination of quantum mechanics calculations and the MM/3D-RISM method suggest that Glu(11) is the preferred oxygen ligand of Zn(2+). However, such preference is dependent on the relative populations of different conformations with specific Zn(2+) binding modes, and therefore could be shifted when experimental or simulation conditions are altered.

Structures and free energy landscapes of the A53T mutant-type α-synuclein protein and impact of A53T mutation on the structures of the wild-type α-synuclein protein with dynamics

The A53T genetic missense mutation of the wild-type α-synuclein (αS) protein was initially identified in Greek and Italian families with familial Parkinson's disease. Detailed understanding of the structures and the changes induced in the wild-type αS structure by the A53T mutation, as well as establishing the direct relationships between the rapid conformational changes and free energy landscapes of these intrinsically disordered fibrillogenic proteins, helps to enhance our fundamental knowledge and to gain insights into the pathogenic mechanism of Parkinson's disease. We employed extensive parallel tempering molecular dynamics simulations along with thermodynamic calculations to determine the secondary and tertiary structural properties as well as the conformational free energy surfaces of the wild-type and A53T mutant-type αS proteins in an aqueous solution medium using both implicit and explicit water models. The confined aqueous volume effect in the simulations of disordered proteins using an explicit model for water is addressed for a model disordered protein. We also assessed the stabilities of the residual secondary structure component interconversions in αS based on free energy calculations at the atomic level with dynamics using our recently developed theoretical strategy. To the best of our knowledge, this study presents the first detailed comparison of the structural properties linked directly to the conformational free energy landscapes of the monomeric wild-type and A53T mutant-type α-synuclein proteins in an aqueous solution environment. Results demonstrate that the β-sheet structure is significantly more altered than the helical structure upon A53T mutation of the monomeric wild-type αS protein in aqueous solution. The β-sheet content close to the mutation site in the N-terminal region is more abundant while the non-amyloid-β component (NAC) and C-terminal regions show a decrease in β-sheet abundance upon A53T mutation. Obtained results utilizing our new theoretical strategy show that the residual secondary structure conversion stabilities resulting in α-helix formation are not significantly affected by the mutation. Interestingly, the residual secondary structure conversion stabilities show that secondary structure conversions resulting in β-sheet formation are influenced by the A53T mutation and the most stable residual transition yielding β-sheet occurs directly from the coil structure. Long-range interactions detected between the NAC region and the N- or C-terminal regions of the wild-type αS disappear upon A53T mutation. The A53T mutant-type αS structures are thermodynamically more stable than those of the wild-type αS protein structures in aqueous solution. Overall, the higher propensity of the A53T mutant-type αS protein to aggregate in comparison to the wild-type αS protein is related to the increased β-sheet formation and lack of strong intramolecular long-range interactions in the N-terminal region in comparison to its wild-type form. The specific residual secondary structure component stabilities reported herein provide information helpful for designing and synthesizing small organic molecules that can block the β-sheet forming residues, which are reactive toward aggregation.

Structures of the E46K mutant-type α-synuclein protein and impact of E46K mutation on the structures of the wild-type α-synuclein protein

The E46K genetic missense mutation of the wild-type α-synuclein protein was recently identified in a family of Spanish origin with hereditary Parkinson's disease. Detailed understanding of the structures of the monomeric E46K mutant-type α-synuclein protein as well as the impact of the E46K missense mutation on the conformations and free energy landscapes of the wild-type α-synuclein are required for gaining insights into the pathogenic mechanism of Parkinson's disease. In this study, we use extensive parallel tempering molecular dynamics simulations along with thermodynamic calculations to assess the secondary and tertiary structural properties as well as the conformational preferences of the monomeric wild-type and E46K mutant-type α-synuclein proteins in an aqueous solution environment. We also present the residual secondary structure component conversion stabilities with dynamics using a theoretical strategy, which we most recently developed. To the best of our knowledge, this study presents the first detailed comparison of the structural and thermodynamic properties of the wild-type and E46K mutant-type α-synuclein proteins in an aqueous solution environment at the atomic level with dynamics. We find that the E46K mutation results not only in local but also in long-range changes in the structural properties of the wild-type α-synuclein protein. The mutation site shows a significant decrease in helical content as well as a large increase in β-sheet structure formation upon E46K mutation. In addition, the β-sheet content of the C-terminal region increases significantly in the E46K mutant-type αS in comparison to the wild-type αS. Our theoretical strategy developed to assess the thermodynamic preference of secondary structure transitions indicates that this shift in secondary structure is the result of a decrease in the thermodynamic preference of turn to helix conversions while the coil to β-sheet preference increases for these residues. Long-range intramolecular protein interactions of the C-terminal with the N-terminal and NAC regions increase upon E46K mutation, resulting in more compact structures for the E46K mutant-type rather than wild-type αS. However, the E46K mutant-type αS structures are less stable than the wild-type αS. Overall, our results show that the E46K mutant-type αS has a higher propensity to aggregate than the wild-type αS and that the N-terminal and C-terminal regions are reactive toward fibrillization and aggregation upon E46K mutation and we explain the associated reasons based on the structural properties herein. Small molecules or drugs that can block the specific residues forming abundant β-sheet structure, which we report here, might help to reduce the reactivity of these intrinsically disordered fibrillogenic proteins toward aggregation and their toxicity.

Structures and free energy landscapes of the wild-type and A30P mutant-type α-synuclein proteins with dynamics

The genetic missense A30P mutation of the wild-type α-synuclein protein results in the replacement of the 30th amino acid residue from alanine (Ala) to proline (Pro) and was initially found in the members of a German family who developed Parkinson's disease. Even though the structures of these proteins have been measured before, detailed understanding about the structures and their relationships with free energy landscapes is lacking, which is of interest to provide insights into the pathogenic mechanism of Parkinson's disease. We report the secondary and tertiary structures and conformational free energy landscapes of the wild-type and A30P mutant-type α-synuclein proteins in an aqueous solution environment via extensive parallel tempering molecular dynamics simulations along with thermodynamic calculations. In addition, we present the residual secondary structure component transition stabilities at the atomic level with dynamics in terms of free energy change calculations using a new strategy that we reported most recently. Our studies yield new interesting results; for instance, we find that the A30P mutation has local as well as long-range effects on the structural properties of the wild-type α-synuclein protein. The helical content at Ala18-Gly31 is less prominent in comparison to the wild-type α-synuclein protein. The β-sheet structure abundance decreases in the N-terminal region upon A30P mutation of the wild-type α-synuclein, whereas the NAC and C-terminal regions possess larger tendencies for β-sheet structure formation. Long-range intramolecular protein interactions are less abundant upon A30P mutation, especially between the NAC and C-terminal regions, which is linked to the less compact and less stable structures of the A30P mutant-type rather than the wild-type α-synuclein protein. Results including the usage of our new strategy for secondary structure transition stabilities show that the A30P mutant-type α-synuclein tendency toward aggregation is higher than the wild-type α-synuclein but we also find that the C-terminal and NAC regions of the A30P mutant-type α-synuclein are reactive toward fibrillzation and aggregation based on atomic level studies with dynamics in an aqueous solution environment. Therefore, we propose that small molecules or drugs blocking the specific residues, which we report herein, located in the NAC- and C-terminal regions of the A30P mutant-type α-synuclein protein might help to reduce the toxicity of the A30P mutant-type α-synuclein protein.

A New Quantum Calibrated Force Field for Zinc-Protein Complex

A quantum calibrated polarizable-charge transfer force field (QPCT) has been proposed to accurately describe the interaction dynamics of zinc-protein complexes. The parameters of the QPCT force field were calibrated by quantum chemistry calculation and capture the polarization and charge transfer effect. QPCTs are validated by molecular dynamic simulation of the hydration shell of the zinc ion, five proteins containing the most common zinc-binding sites (ZnCys2His2, ZnCys3His1, ZnCys4, Zn2Cys6), as well as protein-ligand binding energy in zinc protein MMP3. The calculated results show excellent agreement with the experimental measurement and with results from QM/MM simulation, demonstrating that QPCT is accurate enough to maintain the correct structural integrity of the zinc binding pocket and provide accurate interaction dynamics of the zinc-residue complex. The current approach can also be extended to the study of interaction dynamics of other metal-containing proteins by recalibrating the corresponding parameters to the specific complexes.

The structures of the E22Δ mutant-type amyloid-β alloforms and the impact of E22Δ mutation on the structures of the wild-type amyloid-β alloforms

Structural differences between the intrinsically disordered fibrillogenic wild-type Aβ40 and Aβ42 peptides are linked to Alzheimer's disease. Recently, the E22Δ genetic missense mutation was detected in patients exhibiting Alzheimer's-disease type dementia. However, detailed knowledge about the E22Δ mutant-type Aβ40 and Aβ42 alloform structures as well as the differences from the wild-type Aβ40 and Aβ42 alloform structures is currently lacking. In this study, we present the structures of the E22Δ mutant-type Aβ40 and Aβ42 alloforms as well as the impact of E22Δ mutation on the wild-type Aβ40 and Aβ42 alloform structures. For this purpose, we performed extensive microsecond-time scale parallel tempering molecular dynamics simulations coupled with thermodynamic calculations. For studying the residual secondary structure component transition stabilities, we developed and applied a new theoretical strategy in our studies. We find that the E22Δ mutant-type Aβ40 might have a higher tendency toward aggregation due to more abundant β-sheet formation in the C-terminal region in comparison to the E22Δ mutant-type Aβ42 peptide. More abundant α-helix is formed in the mid-domain regions of the E22Δ mutant-type Aβ alloforms rather than in their wild-type forms. The turn structure at Ala21-Ala30 of the wild-type Aβ, which has been linked to the aggregation process, is less abundant upon E22Δ mutation of both Aβ alloforms. Intramolecular interactions between the N-terminal and central hydrophobic core (CHC), N- and C-terminal, and CHC and C-terminal regions are less abundant or disappear in the E22Δ mutant-type Aβ alloform structures. The thermodynamic trends indicate that the wild-type Aβ42 tends to aggregate more than the wild-type Aβ40 peptide, which is in agreement with experiments. However, this trend is vice versa for the E22Δ mutant-type alloforms. The structural properties of the E22Δ mutant-type Aβ40 and Aβ42 peptides reported herein may prove useful for the development of new drugs to block the formation of toxic E22Δ mutant-type oligomers by either stabilizing helical or destabilizing β-sheet structure in the C-terminal region of these two mutant alloforms.