Tautomeric Effect of Histidine on the Monomeric Structure of Amyloid β-Peptide(1–40)

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.

Delving into the crucial role of the initial structure in the dynamic and self-assembly of amyloid beta

Alzheimer's disease involves the accumulation of amyloid beta (Aβ) monomers that form oligomers and fibrils in the brain. Studying the Aβ monomer is critical for understanding Aβ assembly and peptide behavior and has implications for drug design. Choosing a starting structure with a higher aggregation tendency for cost-effective MD studies and drug design is crucial. Previous studies have utilized distinct initial conformations, leading to varying results. Hence, this study was conducted to compare different initial conformations using the same MD simulation protocol to investigate the behavior and oligomerization propensity of different starting structures of Aβ during 1μs. The behavior of the monomers and their self-assembly systems were studied thoroughly, and the results revealed that highly helical Aβ monomers which used as starting structures retain high helix content during the simulation, and their tautomerization states did not cause significant changes in the structure. On the other hand, the Aβ extended and S-shaped monomers displayed the fingerprints of the fibril structure, which is believed to be more favorable for self-assembly. Self-assembly behaviors were seen for three S-shaped and three Aβ extended peptides. However, both conformations did not show stable β-sheet intermolecular interaction. For the Aβ16-22 monomer as a fragment of the Aβ that can assemble into fibrils, the impacts of capping and uncapping on the initial structure were also investigated. The results displayed that capped and uncapped structures can form oligomers with β-sheet at termini. However, in the capped state, β-sheet interactions were more stable and remained relatively longer than uncapped.

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.

Molecular Simulation of the Binding of Amyloid Beta to Apolipoprotein A-I in High-Density Lipoproteins

Disrupted clearance of amyloid beta (Aβ) from the brain enhances its aggregation and formation of amyloid plaques in Alzheimer's disease. The most abundant protein constituent of circulating high-density lipoprotein (HDL) particles, apoA-I, readily crosses the blood-brain barrier from periphery circulation, exhibits low-micromolar binding affinity for soluble, neurotoxic forms of Aβ, and modulates Aβ aggregation and toxicity in vitro. Its highly conserved N-terminal sequence, 42LNLKLLD48 ('LN'), has been proposed as a binding region for Aβ. However, high-resolution structural characterisation of the mechanism of HDL-Aβ interaction is very difficult to attain. Molecular dynamics simulations were conducted to investigate for the first time the interaction of Aβ and the 'LN' segment of apoA-I. Favourable binding of Aβ by HDLs was found to be driven by hydrophobic and hydrogen-bonding interactions predominantly between the 'LN' segment of apoA-I and Aβ. Preferential binding of Aβ may proceed in small, protein-rich HDLs whereby solvent-exposed hydrophobic 'LN' segments of apoA-I interact specifically with Aβ, stabilising it on the HDL surface in a possibly non-amyloidogenic conformation, facilitating effective Aβ clearance. These findings rationalise the potentially therapeutic role of HDLs in reducing Aβ aggregation and toxicity, and of peptide mimics of the apoA-I interacting region in blocking Aβ aggregation.

Survey of the Aβ-peptide structural diversity: molecular dynamics approaches

The review deals with the application of Molecular Dynamics (MD) to the structure modeling of beta-amyloids (Aβ), currently classified as intrinsically disordered proteins (IDPs). In this review, we strive to relate the main advances in this area but specifically focus on the approaches and methodology. All relevant papers on the Aβ modeling are cited in the Tables in Supplementary Data, including a concise description of the applied approaches, sorted according to the types of the studied systems: modeling of the monomeric Aβ and Aβ aggregates. Similar sections focused according to the type of modeled object are present in the review. In the final part of the review, novel methods of general IDP modeling not confined to Aβ are described.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-024-01253-y.

Edge Substitution Effects of Histidine Tautomerization Behaviors on the Structural Properties and Aggregation Properties of Aβ(1-42) Mature Fibril

Histidine behaviors play critical roles in folding and misfolding processes due to the changes in net charge and the various N/N-H orientations on imidazole rings. However, the effect of histidine tautomerization (HIE (Nε-H, ε) and HID (Nδ-H, δ) states) behaviors on the edge chain of Aβ mature fibrils remains inadequately understood, which is critical for finding a strategy to disturb fibril elongation and growth. In the current study, eight independent molecular dynamics simulations were conducted to investigate such impacts on the structural and aggregation properties. Our results from three different binding models revealed that the binding contributions of edge substitution effects are primarily located between chains 1 and 2. Histidine states significantly influence the secondary structure of each domain. Further analysis confirmed that the C1_H6//C1_E11 intrachain interaction is essential in maintaining the internal stability of chain 1, while the C1_H13//C2_H13 and C1_H14//C2_H13 interchain interactions are critical in maintaining the interchain stability of the fibril structure. Our subsequent analysis revealed that the current edge substitution leads to the loss of the C1_H13//C1_E11 intrachain and C1_H13//C2_H14 interchain interactions. The N-terminal regularity was significantly directly influenced by histidine states, particularly by the residue of C1_H13. Our study provides valuable insights into the effect of histidine behaviors on the edge chain of Aβ mature fibril, advancing our understanding of the histidine behavior hypothesis in misfolding diseases.

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.

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.

New Insights into the Structural and Binding Properties on Aβ Mature Fibrils Due to Histidine Protonation Behaviors

Histidine tautomeric behaviors have been considered origin factors for controlling the structure and aggregation properties of misfolding peptides. Except for tautomeric behaviors, histidine protonation behaviors definitely have the same capacities due to the net charge changes and the various N/N-H orientations on imidazole rings. However, such phenomena are still unknown. In the current study, Aβ mature fibrils substituted with various protonation states were performed by molecular dynamics simulations to investigate the structure and binding properties. Our results show that all kinds of protonation states can increase the ΔG1 stability and decrease ΔG2 and ΔG3 stabilities. A significantly higher averaged β-sheet content was detected in (εεp), (εpp), and (ppp) fibrils in one, two, and three protonation stages, respectively. Impressively, we found that the substituted fibril with specific protonated states can control the N-terminus structural properties. Further analysis confirmed that H6 and H13 are more important than H14 since the H-bond donor and receptor cooperate among C1/C3/C8_H6, C1/C3/C8_H13, and C1/C3/C8_E11. Furthermore, the mechanism of protonation behaviors was discussed. The current study is helpful for understanding the histidine protonation behaviors on one, two, and three protonation stages, which provides new horizons for exploring the origin of protein folding and misfolding.

Structural properties of Aβ (1-40) peptide in protonation stage of one, two, and three: New insights from the histidine protonation behaviors

Structural properties and aggregation tendency can be significantly influenced by histidine behaviors (histidine on N^δ-H state is defined as δ, likewise, N^ε-H: ε and both N^δ-H and N^ε-H: p). In current study, we investigated structural properties of Aβ(1-40) peptide during protonation evolution stage of one, two, and three by total 19 independent replica exchange molecular dynamics simulations using implicit solvent. Our results show that any kind of protonated state will promote β-sheet structure formation in comparison with deprotonated (εεε). With increase in number of protonation, the lowest β-sheet content increased. The highest averaged β-sheet structure content was detected in (δpδ) (46.0 %), (εpp) (36.8 %), and (ppp) (16.0 %) in each protonation stage. With three β-strand structures, (δpδ) shows more stable features and high hydrophobic properties. Further analysis confirmed that H13 and H14 are more important than H6. Specifically, H13 and H14 have a synergistic effect for structural formations by controlling H-bond networks in H13(p) with V39/V40 and H14(p/δ) with G37/G38. Finally, the Pearson correlation coefficient results confirmed that experimental result (ref. 44) is corresponding to our (εpp) system. Our current study will be conducive to understanding the effects of the histidine behaviors, it provides new insights for exploration protein folding and misfolding processes.

Mechanistic insight into the tautomerization of histidine initiated by water-catalyzed N-H and C-H cleavages

The N-H and C-H activation is of great significance in organic chemistry and chemical industry fields, especially, in the utilization of petroleum raw materials. High NδH (tautomer of natural histidine) content would increase Alzheimer's disease risk. To inhibit this and improve the activation of N-H and C-H bonds, the isomerization mechanism from NδH to NεH of histidine-containing dipeptide catalyzed by water cluster was explored. The results discovered that water cluster assists this reaction by reducing the activation energies from 68.20 to 9.60 kcal mol^-1, and its size not only affects the reaction rate but also determines the reaction pathway in a degree. Moreover, water cluster, taken as a potential green catalyst, is more effective on the reactions involving N-H and C-H bond cleavages than reported common toxic organometallic compounds and has different catalytic mechanisms. This work also provides some theoretical guidance for the modulation of Alzheimer's disease induced by histidine isomerization.

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.

On Aromaticity of the Aromatic α-Amino Acids and Tuning of the NICS Indices to Find the Aromaticity Order

The NICS aromaticity indices of the rings in flexible phenylalanine (Phe), tryptophan (Trp), tyrosine (Tyr), and histidine (His) chiral molecules were analyzed. These molecules have several dozens of conformers, and their rings are slightly non-planar. Therefore, the population-averaged NICS_pav index was defined, and the NICS scans had to be performed with respect to planes found by the least-squares routine. A rule differentiating an obverse and a reverse ring face in aromatic amino acids was formulated. The NICS scan minima corresponding to the obverse and reverse face were unequal, which prompted us to use the term ring face aromaticity/ring face tropicity. It appeared that for Phe, Trp, Tyr, and His, the reverse face has always had higher ring face aromaticity/ring face tropicity than the obverse one. Despite the NICS modifications, uncertainty about the amino acid aromaticity order remained. This motivated us to use the integral INICS index newly proposed by Stanger as well. Then, the following sequence was obtained: Trp(phenyl) > Phe > Trp(pyrrole) > His > Tyr. The juxtaposition of the INICS indices of amino acids with that of some model rings revealed a fair transferability of the values. Finally, analysis of the substituent effect on INICS demonstrated that the aromaticity of Tyr is the lowest due to the strength of the OH group π-electron-donating effect able to perturb enough the ring charge distribution and its magnetic aromaticity. The NICS calculations were executed using the ARONICS program written within the project.

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.

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.

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.

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.

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.

Small-Molecule Targeted Aβ42 Aggregate Degradation: Negatively Charged Small Molecules Are More Promising than the Neutral Ones

Heavy evidence has confirmed that Aβ42 oligomers are the most neurotoxic aggregates and play a critical role in the occurrence and development of Alzheimer's disease by causing functional neuron death, cognitive damage, and dementia. Disordered Aβ42 oligomers are challenging therapeutic targets, and no drug is currently in clinical use that modifies the properties of their monomeric states. Here, a negatively charged molecule (ER), rather than the neutral TS1 one, is identified by a molecular dynamics simulation method to be more capable of binding and sequestering the intrinsically disordered amyloid-β peptide Aβ42 in its soluble pentameric state as well as its monomeric components. Results reveal that the ERs interact with Aβ and inhibit the primary nucleation pathways in its aggregation process in entropic expansion mechanism for both Aβ42 and Aβ40 oligomers but with opposite characteristics of hydrophobic surface area (HSA). The interaction between Aβ42 oligomer and either charged ER or neutral TS1/TS0 characterizes decreased HSA, and the decrease in ER-involved case is highly visible, consistent with the observations from in silico and in vitro studies. By contrast, the presence of these inhibitors causes the HSA of Aβ40 oligomer to change undetectably and there is even a bit of increase in the histidine isomerized Aβ40 oligomer. The HSA distinction between Aβ42 and Aβ40 oligomer is possibly derived from the different effects of M35-inhibitor interaction, which is analogous to the effect of M35 oxidation. In comparison with the neutral TS1/TS0 inhibitors, ER is more prone to bind the residues located in the central (β1) and C-terminal (β2) regions of Aβ42 peptide, two key nucleation regions for Aβ intramolecular folding, intermolecular aggregation, and assembly. Notably, ER can strongly bind the charged residues, such as K16, K28, D23, to greatly disturb the potential stabilizer (e.g., salt-bridge, etc.) in metastable Aβ42 oligomers and protofibrils. These results illustrate the strategy of overcoming Alzheimer's disease from inhibiting its early stage Aβ aggregation with two kinds of small molecules to alter their behavior for therapeutic purposes and strongly recommend paying more attention to the engineering and development of negatively charged inhibitors, the long-term underappreciated ones, targeting the early stage Aβ aggregates.

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.

An α-helix mimetic oligopyridylamide, ADH-31, modulates Aβ42 monomer aggregation and destabilizes protofibril structures: insights from molecular dynamics simulations

Alzheimer's disease (AD), an epidemic growing worldwide due to no effective medical aid available in the market, is a neurological disorder. AD is known to be directly associated with the toxicity of amyloid-β (Aβ) aggregates. In search of potent inhibitors of Aβ aggregation, Hamilton and co-workers reported an α-helix mimetic, ADH-31, which acts as a powerful antagonist of Aβ42 aggregation. To identify the key interactions between protein-ligand complexes and to gain insights into the inhibitory mechanism of ADH-31 against Aβ42 aggregation, molecular dynamics (MD) simulations were performed in the present study. The MD simulations highlighted that ADH-31 showed distinct binding capabilities with residues spanning from the N-terminal to the central hydrophobic core (CHC) region of Aβ42 and restricted the conformational transition of the helix-rich structure of Aβ42 into another form of secondary structures (coil/turn/β-sheet). Hydrophobic contacts, hydrogen bonding and π-π interaction contribute to the strong binding between ADH-31 and Aβ42 monomer. The Dictionary of Secondary Structure of Proteins (DSSP) analysis highlighted that the probability of helical content increases from 38.5% to 50.2% and the turn content reduces from 14.7% to 6.2% with almost complete loss of the β-sheet structure (4.5% to 0%) in the Aβ42 monomer + ADH-31 complex. The per-residue binding free energy analysis demonstrated that Arg5, Tyr10, His14, Gln15, Lys16, Val18, Phe19 and Lys28 residues of Aβ42 are responsible for the favourable binding free energy in Aβ42 monomer + ADH-31 complex, which is consistent with the 2D HSQC NMR of the Aβ42 monomer that depicted a change in the chemical shift of residues spanning from Glu11 to Phe20 in the presence of ADH-31. The MD simulations highlighted the prevention of sampling of amyloidogenic β-strand conformations in Aβ42 trimer in the presence of ADH-31 as well as the ability of ADH-31 to destabilize Aβ42 trimer and protofibril structures. The lower binding affinity between Aβ42 trimer chains in the presence of ADH-31 highlights the destabilization of the Aβ42 trimer structure. Overall, MD results highlighted that ADH-31 inhibited Aβ42 aggregation by constraining Aβ peptides into helical conformation and destabilized Aβ42 trimer as well as protofibril structures. The present study provides a theoretical insight into the atomic level details of the inhibitory mechanism of ADH-31 against Aβ42 aggregation as well as protofibril destabilization and could be implemented in the structure-based drug design of potent therapeutic agents for AD.

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.

Tautomeric Effect of Histidine on β-Sheet Formation of Amyloid Beta 1-40: 2D-IR Simulations

Histidine state (protonated or δ or ε tautomer) has been considered the origin of abnormal misfolding and aggregation of β-amyloid (Aβ). Our previous studies reported that the δδδ isomer of Aβ (1-40) has a greater propensity for β-sheet conformation compared to other isomers. However, direct proof of the tautomeric effect has not been reported. In this context, we calculated histidine site-specific two-dimensional infrared spectroscopy of the δδδ, εεε, and πππ (all protonated histidine) systems within the framework of classical molecular dynamics simulations aiming at connecting our previous results with the current experimental observations. Our results showed that β-sheet formation is favored for the δδδ and πππ tautomers compared with the εεε tautomer, consistent with our previous studies. This result was further supported by contact map analyses and the strength of dipole coupling between the amide-I bonds of each residue. The two-dimensional infrared diagonal trace for each tautomer included three distinctive spectrally resolvable peaks near 1680, 1686, and 1693 cm^-1, as was also observed for histidine dipeptides. However, the peak positions at His6, His13, and His14 did not show a consensus trend with the histidine or protonation state but were instead affected by the presence of surrounding hydrogen bonds. Our study provides a deeper insight into the influence of tautomerism and protonation of histidine residues in Aβ (1-40) on amyloid misfolding and provides a connection between our previous simulations and experimental observations.

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

The self-assembly of amyloid-β (Aβ) peptides into senile plaques in the brain is a hallmark of Alzheimer's disease (AD) pathology. Mutation and histidine tautomerism are considered intrinsic origins in the accumulation of Aβ. As a first step toward understanding the impact of A2V mutation and histidine tautomerism on the Aβ42 isoform, we performed replica-exchange molecular dynamics (REMD) simulations to investigate the effects of histidine tautomerism on the structural properties of A2V Aβ42 peptides. There are generally more β-sheet and less α-helix secondary structures in A2V Aβ42 monomers than in WT Aβ42, implying a higher aggregation tendency in A2V Aβ42, which is consistent with previous studies. The current research will help develop the histidine tautomerism hypothesis of misfolded protein aggregation and eventually elucidate the pathogenesis of AD.

Interactions of a multifunctional di-triazole derivative with Alzheimer's Aβ42monomer and Aβ42protofibril: a systematic molecular dynamics study

The molecular dynamics simulations results highlighted that the multi-target-directed ligand6nstabilizes the native α-helix conformation of the Aβ₄₂monomer and induces a sizable destabilization in the Aβ₄₂protofibril structure.

Structural and Binding Properties on Aβ Mature Fibrils Due to the Histidine Tautomeric Effect

Many studies have focused on histidine behaviors in misfolding diseases. However, histidine behaviors on mature fibrils are still unknown. In the current study, we investigated mature fibrils with various histidine states to understand the structural properties of the histidine tautomeric effect on mature fibrils. Our results show that substituting chain 1 with different histidine states affects Aβ structural properties in A2, D7-G9, H14-Q15, S26-N27, and G33-G37 regions. The binding free energies with substituted fibrils were influenced not only along the axial direction, but also between duplex fibrils. Our results suggest that substituted (εδδ) preferentially disturbed the stability among the current mature fibrils. Further, H-bonded network differences indicate that twisted morphologies in mature fibrils are derived from the position and orientation of the imidazole ring in histidines. Our current study helps to elucidate histidine behaviors on mature fibrils, which will present opportunities to understand the misfolding mechanisms.

Intrinsic Origin of Amyloid Aggregation: Collective Effects of the Mutation and Tautomerism of Histidine

Mutation is considered an important factor in the accumulation of amyloid-β (Aβ), which is a hallmark of Alzheimer's disease (AD). A2V Aβ40 shows a higher aggregation tendency; however, the existing knowledge is not sufficient to explain the mechanism. We performed replica-exchange molecular dynamics simulations (REMD) to investigate the structural properties of A2V Aβ40 monomers and consider the tautomerism of histidine. The collective effects of the mutation and tautomerism leads A2V Aβ40 to much higher β-sheet and lower α-helix contents than WT Aβ40, which may explain the enhanced aggregation kinetics of A2V Aβ40 with respect to WT Aβ40. The current research provides new insights on understanding the pathology of AD.

Tautomerization Effect of Histidines on Oligomer Aggregation of β-Amyloid(1-40/42) during the Early Stage: Tautomerism Hypothesis for Misfolding Protein Aggregation

As the intrinsic origin of the hypothesis for β-amyloid (Aβ) from Alzheimer's disease, histidine behaviors were found to play a crucial role in Aβ aggregation. To investigate the histidine behaviors during the early stage of aggregation, Aβ40/42 pentamers with different histidine isomer states were simulated at the atomic level. Results show that five Aβ40 (δδδ) and Aβ42 (εδδ) monomers can rapidly decrease the aggregation threshold, promote stable pentamer formation, and increase pentamer contents by 51.8% and 56.7%, respectively, as compared with the values of their wild-type (εεε) counterparts. Additionally, pentamers of Aβ40 (δδδ) and Aβ42 (εδδ) have different aggregation pathways and disassembly species, Tr+D and Te+M, during the growth of the pentamer. This work discloses the significance of histidine tautomerization in Aβ aggregation, implying a potential way to control Aβ aggregation and develop the assembly inhibitors.

Molecular interactions with redox sites and salt bridges modulate the anti-aggregatory effect of flavonoid, tannin and cardenolide moieties against amyloid-beta (1-42) in silico

In this study, the interactions of flavonoid, tannin and cardenolide moieties as well as their known metabolites were docked against the apolar NMR structure of the aggregatory amyloid-beta fragment (Aβ_1-42). Results showed that the catechin moiety favorably bound Aβ_1-42 peptide at Asp23, Asn27, Ser26 and Glu22 residues, with chalcone similarly binding the middle region of the peptide. Remarkably, hippuric and ferulic acids exhibited hydrophobic interactions with Aβ_1-42 at the latter portion of the peptide, possibly blocking the salt bridges formed by Glu22-Lys28 which stabilizes Phe19-Gly25, as well as the β-sheet Leu34-Gly38 that are known to exist in peptide aggregation. Meanwhile, the metabolites of hydrolyzable tannins, such as urolithin A and gallic acid, exhibited H-bonding interactions with different residues of Aβ_1-42, including Asp1, Asp23 and hydrophobic interactions by gallic acid planar ring to the Hsd6 residue. The coverage was lessened in pyrogallol, suggesting that gallic acid loses its efficacy when further metabolized. Lastly, the different binding poses of the cardenolide moiety interacted with Hsp6 (protonated His) and Tyr10 via hydrophobic interactions. Due to these interactions, the large polycyclic moiety of the ligand would also block further interactions with Hsd6 (prototropic tautomer of His), Asp7, Ser8 and Gly9 that are integral to His6-His13-His14, Arg5-Asp7and Leu34-Gly38 β-sheets, salt bridges in Glu22-Lys28 and turn conformation Phe19-Gly25. Together, these data suggest that the known metabolites of anthocyanins and hydrolyzable tannins contribute the most effective anti-aggregatory interactions with Aβ_1-42, with an unexpected role for cardiac glycosides such as the cardenolie moiety. These bring to light the important role of metabolism in vivo, and suggests further investigation on the effects of these metabolites when concentrated in vivo.

Tautomeric Effect of Histidine on the Monomeric Structure of Amyloid β-Peptide(1-42)

Tautomeric state of histidine is one of the factors that influence the structural and aggregation properties of amyloid β (Aβ)-peptide in neutral state. It is worth it to uncover the monomeric properties of Aβ(1-42) peptide in comparison with Aβ(1-40) peptide. Our replica-exchange molecular dynamics simulations results show that the sheet content of each tautomeric isomer in Aβ(1-42) monomer is slightly higher than that in Aβ(1-40) monomer except His6(δ)-His13(δ)-His14(δ) (δδδ) isomer, implying higher aggregation tendency in Aβ(1-42), which is in agreement with previous experimental and theoretical studies. Further analysis indicates that (εεε), (εδε), (εδδ), and (δδε) isomers prefer sheet conformation although they are in nondominating states. Particularly, it is confirmed that antiparallel β-sheets of (εδδ) were formed at K16-E22 (22.0-43.9%), N27-A30 except G29 (21.9-40.2%), and M35-I41 except G37 (24.1-43.4%). Furthermore, (εδδ) may be the easiest one to overcome structural transformation due to nonobstructing interactions between K16 and/or L17 and histidine residues. The current study will help to understand the tautomeric effect of Aβ(1-42) peptide to overcome Alzheimer's disease.