Antioxidant potential of glutathione

A theoretical study

Béla Fiser, Milán Szóri, Balázs Jójárt, Róbert Izsák, I. Csizmadia, B. Viskolcz

Research output: Contribution to journalArticle

27 Citations (Scopus)

Abstract

All possible X-H (where X can be C, N, O or S) bond dissociation energies (BDEs) of glutathione (γ-l-glutamyl-l-cysteinyl-glycine, GSH) and its fragments have been calculated by first principle methods, and the antioxidant potential of GSH was revealed to be higher than expected in earlier studies. Electron delocalization was found to have an important influence on the antioxidant potential. All structures were optimized and their harmonic vibrational frequencies were calculated in the gas phase at the B3LYP/6-31G(d) level of theory. Solvent effects were taken into account for optimizations at the same level of theory by applying the conductor-like polarizable continuum model (CPCM). Hydrogen cleavage from glutathione proved that the G3MP2B3 composite method provides results consistent with the experimental values for bond dissociation enthalpies (DH298) of S-H, O-H, C-H, and N-H bonds. In order to replace the G3MP2B3 energies with accurate single point calculations, six density functionals, namely, MPWKCIS, MPWKCIS1K, M06, TPSS1KCIS, TPSSh, and B3LYP, were tested against G3MP2B3 for obtaining accurate bond dissociation energies. The MPWKCIS1K/6-311++G(3df,2p)//B3LYP/6-31G(d) level of theory provides the best correlation with the G3MP2B3 method for BDEs in both phases, and therefore, it is recommended for similar calculations. Gas phase results showed that the O-H bond was the weakest, while in aqueous phase the N-H bond in the ammonium group proved to have the smallest BDE value in the studied system. In both cases, the cleavage of the X-H bond was followed by decarboxylation which was responsible for the energetic preference of these processes over the S-H dissociation, which was regarded as the most favorable one until now. The calculated BDE values showed that in aqueous phase the most preferred H-abstraction site is at the weakest N-H bond (BDEaq = 349.3 kJ mol-1) in the glutamine fragment near the α-carbon. In water, the formation of N-centered radicals compared to S-centered ones (BDEaq = 351.7 kJ mol-1) is more endothermic by 2.54 kJ mol-1, due to decarboxylation. Hydrogen dissociation energies from the α-carbons are also comparable in energy with those of the thiol hydrogen, within computational error. The higher stability of the radicals-except the S-centered ones-is due to various degrees of electron delocalization. In aqueous phase, four quasi-equivalent stable radical centers (the α-carbons, the N-centered radical of the NH2 group, and the S-centered radical) were found which provide the antioxidant behavior of glutathione.

Original languageEnglish
Pages (from-to)11269-11277
Number of pages9
JournalJournal of Physical Chemistry B
Volume115
Issue number38
DOIs
Publication statusPublished - Sep 29 2011

Fingerprint

glutathione
antioxidants
Antioxidants
Glutathione
dissociation
Hydrogen
decarboxylation
Carbon
energy
cysteinylglycine
cleavage
carbon
hydrogen
Gases
fragments
vapor phases
glutamine
Electrons
Vibrational spectra
glycine

ASJC Scopus subject areas

  • Surfaces, Coatings and Films
  • Physical and Theoretical Chemistry
  • Materials Chemistry

Cite this

Antioxidant potential of glutathione : A theoretical study. / Fiser, Béla; Szóri, Milán; Jójárt, Balázs; Izsák, Róbert; Csizmadia, I.; Viskolcz, B.

In: Journal of Physical Chemistry B, Vol. 115, No. 38, 29.09.2011, p. 11269-11277.

Research output: Contribution to journalArticle

Fiser, Béla ; Szóri, Milán ; Jójárt, Balázs ; Izsák, Róbert ; Csizmadia, I. ; Viskolcz, B. / Antioxidant potential of glutathione : A theoretical study. In: Journal of Physical Chemistry B. 2011 ; Vol. 115, No. 38. pp. 11269-11277.
@article{2da34a183ad9465c94663425e35b6380,
title = "Antioxidant potential of glutathione: A theoretical study",
abstract = "All possible X-H (where X can be C, N, O or S) bond dissociation energies (BDEs) of glutathione (γ-l-glutamyl-l-cysteinyl-glycine, GSH) and its fragments have been calculated by first principle methods, and the antioxidant potential of GSH was revealed to be higher than expected in earlier studies. Electron delocalization was found to have an important influence on the antioxidant potential. All structures were optimized and their harmonic vibrational frequencies were calculated in the gas phase at the B3LYP/6-31G(d) level of theory. Solvent effects were taken into account for optimizations at the same level of theory by applying the conductor-like polarizable continuum model (CPCM). Hydrogen cleavage from glutathione proved that the G3MP2B3 composite method provides results consistent with the experimental values for bond dissociation enthalpies (DH298) of S-H, O-H, C-H, and N-H bonds. In order to replace the G3MP2B3 energies with accurate single point calculations, six density functionals, namely, MPWKCIS, MPWKCIS1K, M06, TPSS1KCIS, TPSSh, and B3LYP, were tested against G3MP2B3 for obtaining accurate bond dissociation energies. The MPWKCIS1K/6-311++G(3df,2p)//B3LYP/6-31G(d) level of theory provides the best correlation with the G3MP2B3 method for BDEs in both phases, and therefore, it is recommended for similar calculations. Gas phase results showed that the O-H bond was the weakest, while in aqueous phase the N-H bond in the ammonium group proved to have the smallest BDE value in the studied system. In both cases, the cleavage of the X-H bond was followed by decarboxylation which was responsible for the energetic preference of these processes over the S-H dissociation, which was regarded as the most favorable one until now. The calculated BDE values showed that in aqueous phase the most preferred H-abstraction site is at the weakest N-H bond (BDEaq = 349.3 kJ mol-1) in the glutamine fragment near the α-carbon. In water, the formation of N-centered radicals compared to S-centered ones (BDEaq = 351.7 kJ mol-1) is more endothermic by 2.54 kJ mol-1, due to decarboxylation. Hydrogen dissociation energies from the α-carbons are also comparable in energy with those of the thiol hydrogen, within computational error. The higher stability of the radicals-except the S-centered ones-is due to various degrees of electron delocalization. In aqueous phase, four quasi-equivalent stable radical centers (the α-carbons, the N-centered radical of the NH2 group, and the S-centered radical) were found which provide the antioxidant behavior of glutathione.",
author = "B{\'e}la Fiser and Mil{\'a}n Sz{\'o}ri and Bal{\'a}zs J{\'o}j{\'a}rt and R{\'o}bert Izs{\'a}k and I. Csizmadia and B. Viskolcz",
year = "2011",
month = "9",
day = "29",
doi = "10.1021/jp2049525",
language = "English",
volume = "115",
pages = "11269--11277",
journal = "Journal of Physical Chemistry B Materials",
issn = "1520-6106",
publisher = "American Chemical Society",
number = "38",

}

TY - JOUR

T1 - Antioxidant potential of glutathione

T2 - A theoretical study

AU - Fiser, Béla

AU - Szóri, Milán

AU - Jójárt, Balázs

AU - Izsák, Róbert

AU - Csizmadia, I.

AU - Viskolcz, B.

PY - 2011/9/29

Y1 - 2011/9/29

N2 - All possible X-H (where X can be C, N, O or S) bond dissociation energies (BDEs) of glutathione (γ-l-glutamyl-l-cysteinyl-glycine, GSH) and its fragments have been calculated by first principle methods, and the antioxidant potential of GSH was revealed to be higher than expected in earlier studies. Electron delocalization was found to have an important influence on the antioxidant potential. All structures were optimized and their harmonic vibrational frequencies were calculated in the gas phase at the B3LYP/6-31G(d) level of theory. Solvent effects were taken into account for optimizations at the same level of theory by applying the conductor-like polarizable continuum model (CPCM). Hydrogen cleavage from glutathione proved that the G3MP2B3 composite method provides results consistent with the experimental values for bond dissociation enthalpies (DH298) of S-H, O-H, C-H, and N-H bonds. In order to replace the G3MP2B3 energies with accurate single point calculations, six density functionals, namely, MPWKCIS, MPWKCIS1K, M06, TPSS1KCIS, TPSSh, and B3LYP, were tested against G3MP2B3 for obtaining accurate bond dissociation energies. The MPWKCIS1K/6-311++G(3df,2p)//B3LYP/6-31G(d) level of theory provides the best correlation with the G3MP2B3 method for BDEs in both phases, and therefore, it is recommended for similar calculations. Gas phase results showed that the O-H bond was the weakest, while in aqueous phase the N-H bond in the ammonium group proved to have the smallest BDE value in the studied system. In both cases, the cleavage of the X-H bond was followed by decarboxylation which was responsible for the energetic preference of these processes over the S-H dissociation, which was regarded as the most favorable one until now. The calculated BDE values showed that in aqueous phase the most preferred H-abstraction site is at the weakest N-H bond (BDEaq = 349.3 kJ mol-1) in the glutamine fragment near the α-carbon. In water, the formation of N-centered radicals compared to S-centered ones (BDEaq = 351.7 kJ mol-1) is more endothermic by 2.54 kJ mol-1, due to decarboxylation. Hydrogen dissociation energies from the α-carbons are also comparable in energy with those of the thiol hydrogen, within computational error. The higher stability of the radicals-except the S-centered ones-is due to various degrees of electron delocalization. In aqueous phase, four quasi-equivalent stable radical centers (the α-carbons, the N-centered radical of the NH2 group, and the S-centered radical) were found which provide the antioxidant behavior of glutathione.

AB - All possible X-H (where X can be C, N, O or S) bond dissociation energies (BDEs) of glutathione (γ-l-glutamyl-l-cysteinyl-glycine, GSH) and its fragments have been calculated by first principle methods, and the antioxidant potential of GSH was revealed to be higher than expected in earlier studies. Electron delocalization was found to have an important influence on the antioxidant potential. All structures were optimized and their harmonic vibrational frequencies were calculated in the gas phase at the B3LYP/6-31G(d) level of theory. Solvent effects were taken into account for optimizations at the same level of theory by applying the conductor-like polarizable continuum model (CPCM). Hydrogen cleavage from glutathione proved that the G3MP2B3 composite method provides results consistent with the experimental values for bond dissociation enthalpies (DH298) of S-H, O-H, C-H, and N-H bonds. In order to replace the G3MP2B3 energies with accurate single point calculations, six density functionals, namely, MPWKCIS, MPWKCIS1K, M06, TPSS1KCIS, TPSSh, and B3LYP, were tested against G3MP2B3 for obtaining accurate bond dissociation energies. The MPWKCIS1K/6-311++G(3df,2p)//B3LYP/6-31G(d) level of theory provides the best correlation with the G3MP2B3 method for BDEs in both phases, and therefore, it is recommended for similar calculations. Gas phase results showed that the O-H bond was the weakest, while in aqueous phase the N-H bond in the ammonium group proved to have the smallest BDE value in the studied system. In both cases, the cleavage of the X-H bond was followed by decarboxylation which was responsible for the energetic preference of these processes over the S-H dissociation, which was regarded as the most favorable one until now. The calculated BDE values showed that in aqueous phase the most preferred H-abstraction site is at the weakest N-H bond (BDEaq = 349.3 kJ mol-1) in the glutamine fragment near the α-carbon. In water, the formation of N-centered radicals compared to S-centered ones (BDEaq = 351.7 kJ mol-1) is more endothermic by 2.54 kJ mol-1, due to decarboxylation. Hydrogen dissociation energies from the α-carbons are also comparable in energy with those of the thiol hydrogen, within computational error. The higher stability of the radicals-except the S-centered ones-is due to various degrees of electron delocalization. In aqueous phase, four quasi-equivalent stable radical centers (the α-carbons, the N-centered radical of the NH2 group, and the S-centered radical) were found which provide the antioxidant behavior of glutathione.

UR - http://www.scopus.com/inward/record.url?scp=80053111570&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=80053111570&partnerID=8YFLogxK

U2 - 10.1021/jp2049525

DO - 10.1021/jp2049525

M3 - Article

VL - 115

SP - 11269

EP - 11277

JO - Journal of Physical Chemistry B Materials

JF - Journal of Physical Chemistry B Materials

SN - 1520-6106

IS - 38

ER -