Abstract
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.