Abstract
Non-reactive, comparative (2 × 1.2 μs) molecular dynamics simulations were carried out to characterize the interactions between glutathione (GSH, host molecule) and hydroxyl radical (OH(•), guest molecule). From this analysis, two distinct steps were identified in the recognition process of hydroxyl radical by glutathione: catching and steering, based on the interactions between the host-guest molecules. Over 78% of all interactions are related to the catching mechanism via complex formation between anionic carboxyl groups and the OH radical, hence both terminal residues of GSH serve as recognition sites. The glycine residue has an additional role in the recognition of OH radical, namely the steering. The flexibility of the Gly residue enables the formation of further interactions of other parts of glutathione (e.g. thiol, α- and β-carbons) with the lone electron pair of the hydroxyl radical. Moreover, quantum chemical calculations were carried out on selected GSH/OH(•) complexes and on appropriate GSH conformers to describe the energy profile of the recognition process. The relative enthalpy and the free energy changes of the radical recognition of the strongest complexes varied from -42.4 to -27.8 kJ/mol and from -21.3 to 9.8 kJ/mol, respectively. These complexes, containing two or more intermolecular interactions, would be the starting configurations for the hydrogen atom migration to quench the hydroxyl radical via different reaction channels.