Impact of A2V Mutation and Histidine Tautomerism on Aβ42 Monomer Structures from Atomistic Simulations

Histidine tautomerism-mediated transthyretin amyloidogenesis: A molecular insight

Characterization of the conformational alterations involved in monomer misfolding is essential for elucidating the molecular basis of the initial stage of protein accumulation. Here, we report the first structural analyses of transthyretin (TTR) (26-57) fragments with two histidine tautomeric states (δ; N_δ1H and ε; N_ε2H) using replica-exchange molecular dynamics (REMD) simulations. Explaining the organizational properties and misfolding procedure is challenging because the δ and ε configurations can occur in the free neutral state. REMD revealed that β-sheet generation is favored for the δδ (16.8%) and εδ (6.7%) tautomeric isomers, showing frequent main-chain contacts between the stable regions near the head (N-terminus) and central (middle) part compared to the εε (4.8%) and δε (2.8%) isomers. The presence of smaller and wider local energy minima may be related to the structural stability and toxicity of δδ/εδ and εε/δε. Histidines31 and 56 were the parts of regular (such as β-strand) and nonregular (such as coil) secondary structures within the highly toxic TTR isomer. For TTR amyloidosis, focusing on hazardous isomeric forms with high sheet contents may be a potent treatment strategy. Overall, our findings support the tautomerism concept and aid in our comprehension of the basic tautomeric actions of neutral histidine throughout the misfolding process.

Exposure to the electric field: A potential way to block the aggregation of histidine tautomeric isomers of β-amyloid

Controlling protein misfolding and accumulation in neurodegeneration is a challenge in chemical neuroscience. The application of appropriate electric fields (EFs) can be a potential noninvasive therapy to treat neuro disorders. The effect of EFs of varying intensities and directions on the conformational dynamics of β-Amyloid40 (Aβ40) under histidine tautomerism has been investigated for the first time. Our findings suggest that peptides tend to align their dipole moments with the orientation of EF. Irrespective of the EF direction, the dipole moment magnitude is affected by the EF strength. With the conformational changes, the EF strength equal to 0.5 V/nm destroyed the β-sheet content of the δδδ isomer as a potentially toxic agent. The content of the alpha-helical structure which can be transformed into the β-sheet is reduced. The strength of the EF showed a significant influence on the reduction of the number of intra-protein hydrogen bonds especially when EF is equal to 0.5 V/nm which could facilitate destabilization of the structure of the peptides. Current findings provide quantitative insights into the tautomerization-mediated Aβ40 dynamic and conformational changes induced by the external EFs in aqueous solutions, which may provide beneficial information for use as a therapeutic technique.

Mechanistic insights into the conformational switch in profilin-1 subject to collective effects of mutation and histidine tautomerism

Mutations and histidine (His) tautomerism in profilin-1 (PFN1) are associated with amyotrophic lateral sclerosis (ALS). The conformational changes in PFN1 caused by the collective effects of G117V mutation and His tautomeric isomers εε, εδ, δε, and δδ were clarified using molecular dynamics (MD) simulations. The predominant structural variations were seen in α-helices, β-sheets, turns, and coils and the His tautomer's unique degree of disruption was seen in these conformations. The content of α-helices was 23.2 % in the εε and δδ isomers, but the observed α-helices in the isomers εδ and δε were 20.3 % and 21.7 % respectively. The percentage of β-sheet was found to be higher (34.1) in the εε isomer than in the εδ, δε, and δδ isomers, and the values were 30.4, 29.7, and 31.9, respectively. Intermolecular water dynamics analysis discloses that His 133 can form an intramolecular H-bond interaction (N^α-H---N^δ), confirming the experimental observations in the simulations of εε, δε, and δδ isomers of G117V PFN1 mutant. It was concluded that these solvent molecules are crucial for aggregation and must be considered in future research on the PFN1 associated with ALS. Overall, the study offers a thorough microscopic understanding of the pathogenic mechanisms behind conformational changes that cause aggregation illnesses like ALS.

Monitoring early-stage β-amyloid dimer aggregation by histidine site-specific two-dimensional infrared spectroscopy in a simulation study

Monitoring early-stage β-amyloid (Aβ) dimerization is a formidable challenge for understanding neurological diseases. We compared β-sheet formation and histidine site-specific two-dimensional infrared (2D IR) spectroscopic signatures of Aβ dimers with different histidine states (δ; N^δ1-H, ε; N^ε2-H, or π; both protonated). Molecular dynamics (MD) simulations revealed that β-sheet formation is favored for the δδδ:δδδ and πππ:πππ tautomeric isomers showing strong couplings and frequent contacts between the central hydrophobic core and C-terminus compared with the εεε:εεε isomer. Characteristic blue-shifts in the 2D IR central bands were observed upon monomer-dimer transformation. The εεε:εεε dimer exhibited larger frequency shifts than δδδ:δδδ and πππ:πππ implying that the red-shift may have a correlation with N^δ1-H(δ) protonation. Our results support the tautomerization/protonation hypothesis that attributes Aβ misfolding to histidine tautomers as a possible primary initiator for Aβ aggregation and facilitates the application of histidine site-specific 2D IR spectroscopy for studying early-stage Aβ self-assembly.

Histidine tautomerism dependent conformational transitions driven aggregation of profilin-1: Implications in amyotrophic lateral sclerosis

Aggregation of profilin-1 (PFN1) causes a fatal neurodegenerative disease, familial amyotrophic lateral sclerosis (fALS). Histidine (His) tautomerism has been linked to the formation of fibril aggregation causing neurodegenerative disease. Characterization of intermediate species that form during aggregation is crucial, however, this has proven very challenging for experimentalists due to their transient nature. Hence, molecular dynamics (MD) simulations have been performed on the His tautomeric isomers εε, εδ, δε, and δδ of PFN1 to explain the structural changes and to correlate them with its aggregation propensity. MD simulations show that His133 presumably plays a major role in the aggregation of PFN1 upon His tautomerism compared to His119. Further, the formation of a new 3₁₀-helix is observed in εε and δε but 3₁₀-helix is not observed in δδ and εδ isomers. In addition, our findings unveil that β-sheet dominating conformations are observed in His119(δ)-His133(δ) δδ isomer of PFN1 with significant antiparallel β-sheets between residues T15-G23, S29-A33, L63-L65, Q68-S76, F83-T89, T97-T105, and K107-K115, suggesting a novel aggregation mechanism possibly occur for the formation of PFN1 aggregates. Overall, these results propose that MD simulations of PFN1 His tautomers can provide a detailed microscopic understanding of the aggregation mechanisms which are hard to probe through experiments.

Reaction Mechanisms of Histidine and Carnosine with Hypochlorous Acid Along with Chlorination Reactivity of N-Chlorinated Intermediates: A Computational Study

Hypochlorous acid (HOCl) released from activated leukocytes not only plays a significant role in the human immune system but is also implicated in numerous diseases including atherosclerosis and some cancers due to its inappropriate production. Histidine (His) and carnosine (Car), as a respective mediator and protective agent of HOCl damage, have attracted considerable attention; however, their detailed reaction mechanisms are still unclear. In this study, using a His residue with two peptide bond groups (HisRes) as a model, the reaction mechanisms of HisRes and Car including N_εH and N_δH tautomers with HOCl along with the chlorination reactivity of N-chlorinated intermediates were investigated by quantum chemical methods. The obtained results indicate that in the imidazole side chain, the pyridine-like N is the most reactive site rather than the pyrrole-like N, and the kinetic order of all of the possible reaction sites in HisRes follows pyridine-like N > imidazole C_δ ≫ imidazole C_ε > pyrrole-like N, while that in Car is pyridine-like N ≫ imidazole C_δ ≫ amide N. As for N-chlorinated intermediates at imidazole, although the unprotonated form has a low chlorination reactivity as expected, it can still chlorinate tyrosine. Especially, the protonated form exhibits similar ability to HOCl, causing secondary damage in vivo. N-Chlorinated Car features higher internal chlorine migration ability than its intermolecular transchlorination, preventing further HOCl-induced damage. Additionally, a generally overlooked nucleophilic Cl^- shift is also found in N-chlorinated Car/HisRes, indicating that nucleophilic sites in biomolecules also need to be considered. The outcomes of this study are expected to expand our understanding of secondary damage and protective mechanisms involved in HOCl in humans.

Mechanistic Insights into the Polymorphic Associations and Cross-Seeding of Aβ and hIAPP in the Presence of Histidine Tautomerism: An All-Atom Molecular Dynamic Study

Hundreds of millions of people around the world have been affected by Type 2 diabetes (T2D) which is a metabolic disorder. Clinical research has revealed T2D as a possible risk factor for Alzheimer’s disease (AD) development (and vice versa). Amyloid-β (Aβ) and human islet amyloid polypeptide are the main pathological species in AD and T2D, respectively. However, the mechanisms by which these two amyloidogenic peptides co-aggregate are largely uninvestigated. Herein, for the first time, we present the cross-seeding between Amylin1-37 and Aβ40 considering the particular effect of the histidine tautomerism at atomic resolution applying the all-atom molecular dynamics (MD) simulations for heterodimeric complexes. The results via random seed MD simulations indicated that the Aβ40(δδδ) isomer in cross-talking with Islet(ε) and Islet(δ) isomers could retain or increase the β-sheet content in its structure that may make it more prone to further aggregation and exhibit higher toxicity. The other tautomeric isomers which initially did not have a β-sheet structure in their monomeric forms did not show any generated β-sheet, except for one seed of the Islet(ε) and Aβ40(εεε) heterodimers complex that displayed a small amount of formed β-sheet. This computational research may provide a different point of view to examine all possible parameters that may contribute to the development of AD and T2D and provide a better understanding of the pathological link between these two severe diseases.

Dynamics of Amyloid Formation from Simplified Representation to Atomistic Simulations

Amyloid fibril formation is an intrinsic property of short peptides, non-disease proteins, and proteins associated with neurodegenerative diseases. Aggregates of the Aβ and tau proteins, the α-synuclein protein, and the prion protein are observed in the brain of Alzheimer's, Parkinson's, and prion disease patients, respectively. Due to the transient short-range and long-range interactions of all species and their high aggregation propensities, the conformational ensemble of these devastating proteins, the exception being for the monomeric prion protein, remains elusive by standard structural biology methods in bulk solution and in lipid membranes. To overcome these limitations, an increasing number of simulations using different sampling methods and protein models have been performed. In this chapter, we first review our main contributions to the field of amyloid protein simulations aimed at understanding the early aggregation steps of short linear amyloid peptides, the conformational ensemble of the Aβ40/42 dimers in bulk solution, and the stability of Aβ aggregates in lipid membrane models. Then we focus on our studies on the interactions of amyloid peptides/inhibitors to prevent aggregation, and long amyloid sequences, including new results on a monomeric tau construct.

Structural analysis of dexrazoxane: Exploring tautomeric conformations

Structural analysis of dexrazoxane, as a cardioprotective agent, was done in this work by exploring formations of tautomeric conformations and investigating the corresponding effects. Density functional theory (DFT) calculations were performed to optimize the structures to evaluate their molecular and atomic descriptors. In addition to the original structure of dexrazoxane, eight tautomers were obtained with lower stability than the original compound. Movements of two hydrogen atoms in between nitrogen and oxygen atoms of heterocyclic ring put such significant effects. Moreover, electronic molecular orbital features showed effects of such tautomerism processes on distribution patterns and surfaces, in which evaluating the quadrupole coupling constants helped to show the role of atomic sites for resulting the features. As a consequence, the results indicated that the tautomeric formations could significantly change the features of dexrazoxane reminding the importance of carful medication of this drug for patients.

Theoretical Insights into Mutation and Histidine Tautomerism Effects on Tau Proteins

Research on misfolding of tau proteins will help to better understand the formation process of neurofibrillary tangles, a hallmark of Alzheimer's disease. Mutation and histidine tautomeric effects have been considered the two most important inherent factors for tau protein misfolding. In current research, replica-exchange molecular dynamics (REMD) were performed to characterize the structural properties of the key fragment R3 of tau protein under the collective effects of P332L mutation and histidine tautomerism. Simulation results suggest that though the content β-sheet of P332L R3 εδ isomer is slightly lower than that of the WT P332L R3 fragment, the total stable secondary structures including β-sheet and helix of P332L R3 isomers are generally more prevalent than those of wild type R3, which may be the reason that P332L R3 has a higher aggregation tendency. Further analysis showed that the hydrogen bond networks are affected by the mutation and histidine tautomerism. Furthermore, the interactions between N-terminus and C-terminus play a crucial role in β-hairpin formation in all isomers. The current study will contribute to revealing the collective effects of P332L and histidine tautomerism on the misfolding of tau proteins.

Unraveling the Histidine Tautomerism Effect on the Initial Stages of Prion Misfolding: New Insights from a Computational Perspective

The aggregation and structural conversion of normal prion peptide (PrP^C) into the pathogenic scrapie form (PrP^Sc), which can act as a seed to enhance prion amyloid fiber formation, is believed to be a crucial event in prionopathies. Previous research suggests that the prion monomer may play an important role in oligomer generation during disease pathogenesis. In the present study, extensive replica-exchange molecular dynamics (REMD) simulations were conducted to explore the conformational characteristics of the huPrP (125-160) monomer under the histidine tautomerism effect. Investigating the structural characteristics and fibrilization process is challenging because two histidine tautomers [N_ε2-H (ε) and N_δ1-H (δ)] can occur in the open neutral state. Molecular dynamics (MD) simulation outcomes have shown that the toxic εδ and δδ isomer (containing several and broader local minima) had the highest α-helix structures, with contents of 21.11% and 21.01%, respectively, and may have a strong influence on the organizational behavior of a monomeric prion. The amino acids aspartate 20 (D20)-asparagine 29 (N29) and isoleucine 15 (I15)-histidine 16 (H16), D20-arginine 27 (R27) as well as N29 formed α-helix with the highest probabilities in the δδ and εδ isomer, accordingly. On the basis of our findings, we propose the histidine tautomerization hypothesis as a new prion accumulation mechanism, which may exist to induce the formation of prion accumulates. Overall, our tautomerism hypothesis constitutes a promising perspective for enhancing understanding of prion disease pathobiology and may help in the design of a good inhibitor.

Role of the English (H6R) Mutation on the Structural Properties of Aβ40 and Aβ42 Owing to the Histidine Tautomeric Effect

As an intrinsic origin cause, histidine behaviors play a critical role in protein misfolding processes. Generally, the English (H6R) mutation will disrupt H6 interactions. However, the structural properties of Aβ40 H6R and Aβ42 H6R under the complex influence of a histidine tautomeric effect and an H6R mutation remain unclear. Therefore, we performed a replica exchange molecular dynamics simulation to unveil such structural properties. Our result showed that the H6R substitute could promote the generation of β-sheet structures in comparison to the wild type. Three β-strand structure properties were observed in Aβ40 (rδδ), Aβ42 (rεε), Aβ42 (rεδ), and Aβ42 (rδδ) with β-sheet contents of 47.5%, 37.2%, 46.9%, and 38.6%, respectively, and the dominant conformational properties of Aβ40 (rδδ), Aβ42 (rεε), Aβ42 (rεδ), and Aβ42 (rδδ) had top conformational states of 86.0%, 73.2%, 67.0%, and 56.5%, respectively. Further analysis confirmed that R6 had different mechanisms for controlling the conformational features in Aβ40 H6R and Aβ42 H6R. In the Aβ40 systems, H14 H-bond networks played a critical role in controlling the structural properties. However, in the Aβ42 systems, R6 was more important because it was directly involved in the β-strand formation and maintained the β-sheet between the N-terminus and the central hydrophobic core region. Our current study helps to elucidate the histidine tautomeric behaviors in H6R mutations, which will present opportunities to understand the correlation between with/without H6 and the Aβ40/Aβ42 H6R misfolding mechanisms.

Molecular insight into the early stage of amyloid-β(1-42) Homodimers aggregation influenced by histidine tautomerism

Aggregated amyloid β-peptide (Aβ) in small oligomeric forms inside the brain causes synaptic function disruption and the development of Alzheimer's disease (AD). Histidine is an important amino acid that may lead to structural changes. Aβ42 monomer chain includes 3 histidine residues that considering two ε and δ tautomers 8 isomers, including (εεε) and (εδδ) could be formed. Molecular dynamics simulation on homodimerization of (εεε) (the most common type of tautomers) and (εδδ) tautomers with different initial configurations using monomer chains from our previous work were performed to uncover the tautomeric behavior of histidine on Aβ42 aggregation in a physiological pH which is still largely unknown and impossible to observe experimentally. We found a higher propensity of forming β-sheet in (εδδ) homodimers and specifically in a greater amount from Aβ42 than from Aβ40. A smaller amount of β-sheet formation was observed for (εεε) homodimers compared with (εδδ). Additionally, interactions in (εδδ) homodimers may indicate the importance of the hydrophobic core and C-/N-terminals during oligomerization. Our findings indicate the important role of the tautomeric effect of histidine and (εδδ) homodimers at the early stage of Aβ aggregation.

Histidine Tautomeric Effect on the Key Fragment R3 of Tau Protein from Atomistic Simulations

Self-assembly of hyperphosphorylated tau proteins into neurofibrillary tangles (NFT) is a hallmark of Alzheimer's disease. Previous studies suggest that the tau monomer may play an important role in NFTs formation in two general categories: inert (M_i) monomer and seed-competent (M_s) monomer. In the current study, replica-exchange molecular dynamics (REMD) were performed to investigate the effect of histidine tautomerism on the structures of a key fragment (R3) of tau protein and the transformation between different conformations. Based on the simulation results, we propose the histidine tautomerism hypothesis for tau protein misfolding. Histidine tautomerism greatly expands the conformational library, which triggers the emergence of conformations with higher aggregation tendency. Moreover, the conversions existing in both isomers and conformations may cause protein misfolding to occur more readily.

The F19W mutation reduces the binding affinity of the transmembrane Aβ11-40 trimer to the membrane bilayer

Alzheimer's disease is linked to the aggregation of the amyloid-β protein (Aβ) of 40 or 42 amino acids. Lipid membranes are known to modulate the rate and mechanisms of the Aβ aggregation. Point mutations in Aβ can alter these rates and mechanisms. In particular, experiments show that F19 mutations influence the aggregation rate, but maintain the fibril structures. Here, we used molecular dynamics simulations to examine the effect of the F19W mutation in the 3Aβ11-40 trimer immersed in DPPC lipid bilayers submerged in aqueous solution. Substituting Phe by its closest (non-polar) aromatic amino acid Trp has a dramatic reduction in binding affinity to the phospholipid membrane (measured with respect to the solvated protein) compared to the wild type: the binding free energy of the protein-DPPC lipid bilayer increases by 40-50 kcal mol-1 over the wild-type. This is accompanied by conformational changes and loss of salt bridges, as well as a more complex free energy surface, all indicative of a more flexible and less stable mutated trimer. These results suggest that the impact of mutations can be assessed, at least partially, by evaluating the interaction of the mutated peptides with the lipid membranes.

Molecular mechanism of amyloidogenicity and neurotoxicity of a pro-aggregated tau mutant in the presence of histidine tautomerism via replica-exchange simulation

Considering ΔK280 tau mutation, δε isomer with highest sheet content may accelerate aggregation; generating small compounds to inhibit this would help tp prevent tauopathies.

Intrinsic Origin of Tau Protein Aggregation: Effects of Histidine Tautomerism on Tau267-312 Monomer

Histidine tautomerism is considered a crucial component that affects the constitutional and accumulation characteristics of the tau_267-312 monomer in the neutral condition, which are connected with the pathobiology of Alzheimer's disease (AD). Interpreting the organizational characteristics and accumulation procedure is a challenging task because two tautomeric conformations (the N^ε-H or N^δ-H tautomer) can occur in the open neutral condition. In the current work, replica-exchange molecular dynamics (REMD) simulations were performed to investigate the structural properties of the tau_267-312 monomer considering the histidine tautomeric effect. Based on the simulation outcomes, the histidine 268 (H268) (δ)-H299 (δ) (δδ) isomer had the highest β-sheet content with a value of 26.2%, which acquires a sheet-governing toxic conformer with the first abundant conformational state of 22.6%. In addition, δδ displayed notable antiparallel β-sheets between lysine 8 (K8)-asparagine 13 (N13) and valine 40 (V40)-tyrosine 44 (Y44) as well as between K32-H33 and V40-Y44 (β-meander supersecondary structure), indicating this tautomeric isomer may exist to stimulate tau oligomerization. Furthermore, H299 was found to play an essential role in the structural stabilization of the δδ isomer compared with H268. The present research will aid in obtaining insight into the organizational and accumulation properties of tau protein in the presence of histidine tautomerism to control AD.