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

Insights into amyloid-β misfolding: The impact of histidine tautomerism and Au(111) surfaces through MD simulations and 2DIR spectroscopy

The structural and dynamic properties of histidine tautomeric isomers (δδδ and εεε) in the amyloid-β (Aβ40) peptide are thoroughly examined, both in their isolated state and in interaction with the Au(111) surface, to explore their behavior in diverse environments. Utilizing molecular dynamics simulations and advanced 2D spectroscopy techniques, we reveal that the Au(111) surface significantly modulates the peptides' conformational flexibility, hydrogen bonding patterns, and secondary structures. Notably, the presence of gold leads to enhanced stability and a reduction in β-sheet formation, especially for the δδδ isomer, favoring the formation of coils and α-helices instead. The analysis, which includes root mean square fluctuation (RMSF), free energy landscape (FEL), root mean square deviation (RMSD), and contact maps, highlights the transformative impact of Au(111) on the peptides' secondary structures and dipole coupling behaviors. Spectroscopic insights, such as frequency shifts and two-dimensional infrared spectroscopy (2DIR), further illustrate the interaction-induced spectral changes and conformational adjustments. These findings underscore the critical role of metal surface interactions in modulating peptide aggregation and stability, offering new insights into Alzheimer's disease mechanisms and potential therapeutic approaches.

Interpreting the Histidine-Containing Small Peptides on Tau Protein Tautomerism: A Theoretical Perspective

Exploring the nature of histidine residue tautomerization via a systematic conformational study is essential for understanding the pathology and toxicity of several neurodegenerative diseases, as well as for their diagnosis and treatment. Herein, density functional theory (DFT) calculations were used to determine the Tau protein's histidine-containing dipeptide (Lys-His, His-Gln, and His-Val) and tripeptide (Lys-His-Gln and Lys-His-Val) isomeric conformations via intramolecular hydrogen bond interactions, with particular attention to the influence of N-H group isomeric forms on their properties. The calculated infrared (IR) spectroscopy of the N-H stretch region of each isomer and nuclear magnetic resonance (NMR) shielding of the imidazole ring carbon atoms (13C1, 13C2, and 13C3) were investigated. The results show that both the IR spectrum of the N-H group and the NMR shielding of 13C nuclei on the imidazole ring can be used to identify the histidine-containing dipeptide and tripeptide tautomeric isomers. Systematically analyzing the hydrogen bonding interactions, the atomic charge distribution, the potential energy distribution, and the HOMO-LUMO transitions of each isomer further verified the above conclusions. This study provides theoretical evidence for the conformation identification of the histidine-containing dipeptide and tripeptide isomers on Tau protein.

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.

Sequence-Dependent Conformational Properties of PGGG Motif in Tau Repeats: Insights from Molecular Dynamics Simulations of Narrow Pick Filament

Tauopathies are a class of neurodegenerative diseases correlated with the presence of pathological Tau fibrils as a diagnostic marker. The microtubule-binding repeat region of Tau protein, which includes R1, R2, R3, and R4 repeats, constitutes the core of these fibrils. Each repeat consists of a semiconserved C-terminal hexapeptide flanked by KxGS and PGGG motifs. Previous studies have shown the influence of these peptides on protein aggregation, yet their repeat-specific properties are less explored. Using molecular dynamics, we probed the sequence-specific influence of the C-terminal hexapeptide (²⁶⁴ENLKHQ²⁶⁹) in determining the compact local conformation of the R1 repeat of the narrow Pick filament (NPF) with a homologous E264G mutation. In addition to that, we also studied the influence of ²⁶²S phosphorylation on this conformation as the phosphorylation is proposed to alleviate the pathogenesis of Pick's disease. Interestingly, we determined that E264G mutation induces a conformational shift of ²⁷⁰PGGG²⁷³ from a turn to a random coil. This conformational dependence is experimentally verified with the R1R3-E264G mutant construct, which displayed accelerated aggregation compared with the R1R3 wild-type construct. A significant delay in aggregation of the R1R3-G326E mutant further demonstrates the importance of ³²⁶G in determining the conformation of the R3 repeat. Thus, we conclude that the conformational properties of the PGGG motif in Tau repeats are strongly dependent on the repeat-specific sequence of the C-terminal hexapeptide.

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.

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.

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.

Copper (Cu2+) ion-induced misfolding of tau protein R3 peptide revealed by enhanced molecular dynamics simulation

Tau misfolding plays a significant role in some neurodegenerative diseases such as Alzheimer's disease (AD). It is intrinsically disordered and highly soluble under normal physiological conditions. While the protein will aggregate to form paired helical filaments (PHFs) under copper homeostasis at pathological conditions, which is the main substance of neurofibrillary tangles (NFTs) in the brain of AD patients. However, the molecular mechanism underlying the copper (Cu2+) ion-induced tau misfolding is not fully understood. In this study, using the 1/2 third repeat fragment (R3 peptide) of tau protein (residues 318-335: VTSKCGSLGNIHHKPGGG) as a model, a Gaussian accelerated molecular dynamics (GaMD) method followed by efficient trajectory analysis was carried out to investigate the influences of Cu2+ on the tau about the protein fold and the free energy landscape along the simulation. The two-dimensional potential of mean force (PMF) profiles obtained from reweighting of the GaMD simulations as well as clustering analysis revealed the Cu2+ ion induced α-helix fold R3 peptide located at the low-energy wells of free energy map, which is in agreement with the reported experimental result. In contrast, there is no α-helix fold of R3 peptide that appeared during the GaMD simulation without Cu2+ ion existing. Furthermore, the definition of secondary structure of protein (DSSP) analysis indicated that the R3 peptide with Cu2+ ion forms a stable structure of the helix (Lys321-His330 interval of the peptide) at between 400 and 500 ns. Therefore, the structures and free energy profiles from GaMD simulations proposed that Cu2+ triggers the aggregation of R3 peptide into toxic PHFs through a stable α-helix fold form.

Histidine Tautomerism Driving Human Islet Amyloid Polypeptide Aggregation in the Early Stages of Diabetes Mellitus Progression: Insight at the Atomistic Level

Early oligomerization of human islet amyloid polypeptide (hIAPP), which is accountable for β-cell death, has been implicated in the progression of type 2 diabetes mellitus. Some researches have shown the connection between hIAPP and Alzheimer's disease as well. However, the mechanism of peptide accumulation and associated cytotoxicity remains unclear. Due to the unique properties and significant role of histidine in protein sequences, here for the first time, the tautomeric effect of histidine at the early stages of amylin misfolding was investigated via molecular dynamics simulations. Considering Tau and Pi tautomeric forms of histidine (Tau and Pi tautomers are denoted as ϵ and δ, respectively), simulations were performed on two possible isomers of amylin. Our analysis revealed a higher probability of transient α-helix generation in the δ isomer in monomeric form. In dimeric forms, the δδ and δϵ conformations showed an elevated amount of α-helix and lower coil in comparison to the ϵϵ dimer. Due to the significant role of α-helix in membrane disruption and transition to β-sheet structure, these results may imply a noticeable contribution of the δ isomer and the δδ and δϵ dimers rather than ϵ and ϵϵ conformations in the early stages of diabetes initiation. Our results may aid in elucidating the hIAPP self-association process in the etiology of amyloidosis.

Misfolding and Self-Assembly Dynamics of Microtubule-Binding Repeats of the Alzheimer-Related Protein Tau

Pathological aggregation of intrinsically disordered tau protein, driven by the interactions between microtubule-binding (MTB) domains, is associated with Alzheimer's disease. The MTB domain contains either three or four repeats with sequence similarities. Compared to amyloid β, many aspects of the misfolding and aggregation mechanisms of tau are largely unknown. In this study, we systematically investigated the dynamics of monomer misfolding and dimerization of each MTB repeat using atomistic discrete molecular dynamic simulations. Our results revealed that all the four repeat monomers (R1-R4) were very dynamic, featuring frequent conformational conversion and lacking stable conformations. While R1, R2, and R4 monomers occasionally adopted partially helical conformations, R3 monomers frequently formed β-sheets. In dimerization simulations, R3 displayed the strongest aggregation propensity with high β-sheet contents, while R1 was the least prone to aggregation. The R2 and R4 dimers contained both helix and β-sheet structures. The β-sheets in R4 assemblies were dominant with β-hairpin conformation. In R2 and R3 dimers, intermolecular β-sheets were mainly driven by residues around the paired helical filament (PHF) regions. Residues around the PHF6* in R2 and PHF6 in R3 had significantly higher intermolecular contacts than other regions, suggesting that these residues play a key role in the amyloid aggregation of tau. Our results on the structural ensembles and early aggregation dynamics of each tau MTB repeat will help understand the nucleation and fibrillization of tau.

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.