Investigation of the role of conserved residues Ser13, Asn48 and Pro49 in the catalytic mechanism of the tau class glutathione transferase from Glycine max

Plant glutathione transferases (GSTs) play a key role in the metabolism of various xenobiotics. In this report, the catalytic mechanism of the tau class GSTU4-4 isoenzyme from Glycine max (GmGSTU4-4) was investigated by site-directed mutagenesis and steady-state kinetic analysis. The catalytic properties of the wild-type enzyme and three mutants of strictly conserved residues (Ser13Ala, Asn48Ala and Pro49Ala) were studied in 1-chloro-2,4-dinitrobenzene (CDNB) conjugation reaction. The results showed that the mutations significantly affect substrate binding and specificity. The effect of Ser13Ala mutation on the catalytic efficiency of the enzyme could be explained by assuming the direct involvement of Ser13 to the reaction chemistry and the correct positioning of GSH and CDNB in the ternary catalytic complex. Asn48 and Pro49 were found to have a direct role on the structural integrity of the GSH-binding site (G-site). Moreover, mutation of Asn48 and Pro49 residues may bring about secondary effects altering the thermal stability and the catalytic activity (k_cat) of the enzyme without affecting the nature of the rate-limiting step of the catalytic reaction.     

Aromatic residues in the C-terminal region of glutathione transferase A1-1 influence rate-determining steps in the catalytic mechanism

Human glutathione transferase A1-1 (GST A1-1) has a flexible C-terminal segment that forms a helix (alpha9) closing the active site upon binding of glutathione and a small electrophilic substrate such as 1-chloro-2,4-dinitrobenzene (CDNB). In the absence of active-site ligands, the C-terminal segment is not fixed in one position and is not detectable in the crystal structure. A key residue in the alpha9-helix is Phe 220, which can interact with both the enzyme-bound glutathione and the second substrate, and possibly guide the reactants into the transition state. Mutation of Phe 220 into Ala and Thr was shown to reduce the catalytic efficiency of GST A1-1. The mutation of an additional residue, Phe 222, caused further decrease in activity. The presence of a viscosogen in the reaction medium decreased the kinetic parameters k(cat) and k(cat)/K(m) for the conjugation of CDNB catalyzed by wild-type GST A1-1, in agreement with the view that product release is rate limiting for the substrate-saturated enzyme. The mutations cause a decrease of the viscosity dependence of both kinetic parameters, indicating that the motion of the alpha9-helix is linked to catalysis in wild-type GST A1-1. The isomerization reaction with the alternative substrate Delta(5)-androstene-3,17-dione (AD) is affected in a similar manner by the viscosogens. The transition state energy of the isomerization reaction, like that of the CDNB conjugation, is lowered by Phe 220 as indicated by the effects of the mutations on k(cat)/K(m). The results demonstrate that Phe 220 and Phe 222, in the dynamic C-terminal segment, influence rate-determining steps in the catalytic mechanism of both the substitution and the isomerization reactions.     

Aromatic residues in the C-terminal region of glutathione transferase A1-1 influence rate-determining steps in the catalytic mechanism

Human glutathione transferase A1-1 (GST A1-1) has a flexible C-terminal segment that forms a helix (alpha 9) closing the active site upon binding of glutathione and a small electrophilic substrate such as 1-chloro-2,4-dinitrobenzene (CDNB). In the absence of active-site ligands, the C-terminal segment is not fixed in one position and is not detectable in the crystal structure. A key residue in the alpha 9-helix is Phe 220, which can interact with both the enzyme-bound glutathione and the second substrate, and possibly guide the reactants into the transition state. Mutation of Phe 220 into Ala and Thr was shown to reduce the catalytic efficiency of GST A1-1. The mutation of an additional residue, Phe 222, caused further decrease in activity. The presence of a viscosogen in the reaction medium decreased the kinetic parameters K(cat) and K(cat)/K(m) for the conjugation of CDNB catalyzed by wild-type GST A1-1, in agreement with the view that product release is rate limiting for the substrate-saturated enzyme. The mutations cause a decrease of the viscosity dependence of both kinetic parameters, indicating that the motion of the alpha 9-helix is linked to catalysis in wild-type GST A1-1. The isomerization reaction with the alternative substrate Delta(5)-androstene-3,17-dione (AD) is affected in a similar manner by the viscosogens. The transition state energy of the isomerization reaction, like that of the CDNB conjugation, is lowered by Phe 220 as indicated by the effects of the mutations on K(cat)/K(m). The results demonstrate that Phe 220 and Phe 222, in the dynamic C-terminal segment, influence rate-determining steps in the catalytic mechanism of both the substitution and the isomerization reactions.     

Naturally occurring Phe151Leu substitution near a conserved folding module lowers stability of glutathione transferase P1-1

Glutathione transferases (GSTs) are a family of enzymes that detoxify electrophilic compounds, such as carcinogens or drugs, by conjugating them to glutathione. The enzymes have contributed to the understanding of protein structure, due to large differences in amino acid sequence within the family, yet similar architecture and folding. Our objective was to conduct a systematic survey of GSTP1 polymorphisms and their function. Nearly all variants detected were known polymorphisms: IVS4+13C>A; Ile105Val; Ala114Val; and g.2596T>C (Ser185Ser). However, we also found a novel Phe151Leu substitution in an African-American subject (1 out of 111). Kinetic parameters for the conjugation reaction with 1-chloro-2,4-dinitrobenzene (CDNB) were determined for the novel variant enzyme purified via heterologous expression in Escherichia coli. Five substrates were used for measurement of specific activities, including isothiocyanate compounds that occur in cruciferous vegetables (benzylisothiocyanate, phenethylisothiocyanate, and sulforaphane). Such isothiocyanate substrates are potential cancer chemopreventive agents that are conjugated by GSTs. No major change in kinetic parameters was observed. However, the half-life at 50 degrees C of the Leu 151 enzyme was reduced to 12 min, as compared to 28 min for the Phe 151 enzyme. Residue 151 is located at the N-terminus of helix alpha6 in GST motif II, surrounded by hydrophobic residues, and near the conserved "hydrophobic staple" and N-capping box motifs. These local structural elements aid in formation of helix alpha6 and promote proper folding and protein stability. Analysis of the three-dimensional structure showed that substitution of Phe 151 with Leu produces a hydrophobic cavity in the GSTP1 core, thereby destabilizing its structure. Phe151Leu represents one of the first-described allelic variations in a protein folding motif.     

Synthesis and characterization of the oxygen and desthio analogues of glutathione as dead-end inhibitors of glutathione S-transferase

The oxygen analogue, gamma-L-Glu-L-SerGly (GOH) and desthio analogue, gamma-L-Glu-L-AlaGly (GH) have been synthesized by a simple three step procedure involving active ester coupling of N-t-BOC-alpha-(4-nitrophenyl)-L-glutamate to L-SerGly and L-AlaGly, respectively. The two peptides are excellent dead-end inhibitors of isozymes 3-3 and 4-4 of rat liver glutathione S-transferase. At low fixed concentrations of 1-chloro-2,4-dinitrobenzene (CDNB) GOH and GH are linear competitive inhibitors of isozyme 3-3 vs glutathione with KI values of 13.0 and 116 microM, respectively. Both peptides are non-competitive (mixed-type) inhibitors vs CDNB when glutathione is the fixed substrate. Similar results are obtained with both peptides and isozyme 4-4. The results rule out ordered or ping-pong kinetic mechanisms where the electrophile adds first.     

Aromatic residues in the C-terminal region of glutathione transferase A1-1 influence rate-determining steps in the catalytic mechanism [Biochimica et Biophysica Acta 1597 (2002) 157-163]

Human glutathione transferase A1-1 (GST A1-1) has a flexible C-terminal segment that forms a helix (α9) closing the active site upon binding of glutathione and a small electrophilic substrate such as 1-chloro-2,4-dinitrobenzene (CDNB). In the absence of active-site ligands, the C-terminal segment is not fixed in one position and is not detectable in the crystal structure. A key residue in the α9-helix is Phe 220, which can interact with both the enzyme-bound glutathione and the second substrate, and possibly guide the reactants into the transition state. Mutation of Phe 220 into Ala and Thr was shown to reduce the catalytic efficiency of GST A1-1. The mutation of an additional residue, Phe 222, caused further decrease in activity. The presence of a viscosogen in the reaction medium decreased the kinetic parameters k_cat and k_cat/K_m for the conjugation of CDNB catalyzed by wild-type GST A1-1, in agreement with the view that product release is rate limiting for the substrate-saturated enzyme. The mutations cause a decrease of the viscosity dependence of both kinetic parameters, indicating that the motion of the α9-helix is linked to catalysis in wild-type GST A1-1. The isomerization reaction with the alternative substrate Δ⁵-androstene-3,17-dione (AD) is affected in a similar manner by the viscosogens. The transition state energy of the isomerization reaction, like that of the CDNB conjugation, is lowered by Phe 220 as indicated by the effects of the mutations on k_cat/K_m. The results demonstrate that Phe 220 and Phe 222, in the dynamic C-terminal segment, influence rate-determining steps in the catalytic mechanism of both the substitution and the isomerization reactions.     

Enzymatic properties of mutant Escherichia coli tryptophan synthase alpha-subunits

Thirty-nine mutant tryptophan synthase alpha subunits have been purified and analyzed (in the presence of the beta 2-subunit) for their enzymatic (kcat, Km) behavior in the reactions catalyzed by the alpha 2.beta 2 complex, the fully constituted form of this enzyme. The mutant alpha subunits, obtained by in vitro random, saturation mutagenesis of the encoding trpA gene, contain single amino acid substitutions at sites within the first 121 residues of the alpha polypeptide. Four categories of altered residues have been tentatively assigned roles in the catalytic functions of this enzyme: 1) catalytic residues (Glu49 and Asp60); 2) residues involved in substrate binding or orientation (Phe22, Thr63, Gln65, Tyr102, and Leu105); 3) residues involved in alpha.beta subunit interactions (Gly51, Pro53, Asp56, Asp60, Pro62, Ala67, Phe72, Thr77, Pro78, Tyr102, Asn104, Leu105, and Asn108); and 4) residues with no apparent catalytic roles. Catalytic residue alterations result in no detectable activity in the alpha-subunit specific reactions. Substrate binding/orientation roles are detected enzymatically primarily as rate defects; alterations only at Tyr102 result in apparent Km effects. alpha.beta interaction roles are detected as rate defects in all tryptophan synthase reactions plus Km increases for the alpha-subunit substrate, indole-3-glycerol phosphate, only when L-serine is present at the beta 2-subunit active site. A substitution at only one site, Asn104, appears to be unique in its potential effect on intersubunit channeling of indole, the product of the alpha-subunit specific reaction, to the beta 2-subunit active site.     

Sigma-class glutathione transferase from Xenopus laevis: molecular cloning, expression, and site-directed mutagenesis

The structural gene for glutathione transferase (XlGSTS1-1) in the amphibia Xenopus laevis has been cloned from an embryo library and its nucleotide sequence has been determined. Open reading frame analysis indicated that xlgsts1 gene encodes the smallest protein of sigma class GST so far identified as being composed of only 194 amino acid residues. The recombinant XlGSTS1-1 shows a narrow range of substrate specificity as well as a significantly lower 1-chloro-2,4-dinitrobenzene conjugation capacity than that of squid sigma class GST. To compare the structural and functional differences between the squid and amphibian enzymes, several site-specific mutations were introduced in XlGSTS1-1, i.e., Ser100Asn, Phe102Tyr, Trp143Leu, Phe146Leu, and Trp148Cys. The results obtained indicate that Trp143 and Trp148 are more important determinants for the structural stability of XlGSTS1-1 rather than for its substrate specificity.     

Tyrosine-7 is an essential residue for the catalytic activity of human class PI glutathione S-transferase: chemical modification and site-directed mutagenesis studies.

The glutathione (GSH)-conjugating activity of human class Pi glutathione S-transferase (GST pi) toward 1-chloro-2,4-dinitrobenzene (CDNB) was significantly lowered by reaction with N-acetylimidazole, an O-acetylating reagent for tyrosine residues. Further, the replacement of Tyr7 in GST pi, which is conserved in all cytosolic GSTs, with phenylalanine by site-directed mutagenesis also lowered the activities toward CDNB and ethacrynic acid. The Km values of the mutant for both GSH and CDNB were almost equivalent to those of the wild type, while the Vmax of the former was about 55-fold smaller than that of the latter. Therefore, Tyr7 is considered to be an essential residue for the catalytic activity of GST pi.     

Naturally occurring Phe151Leu substitution near a conserved folding module lowers stability of glutathione transferase P1–1

Glutathione transferases (GSTs) are a family of enzymes that detoxify electrophilic compounds, such as carcinogens or drugs, by conjugating them to glutathione. The enzymes have contributed to the understanding of protein structure, due to large differences in amino acid sequence within the family, yet similar architecture and folding. Our objective was to conduct a systematic survey of GSTP1 polymorphisms and their function. Nearly all variants detected were known polymorphisms: IVS4+13C>A; Ile105Val; Ala114Val; and g.2596T>C (Ser185Ser). However, we also found a novel Phe151Leu substitution in an African-American subject (1 out of 111). Kinetic parameters for the conjugation reaction with 1-chloro-2,4-dinitrobenzene (CDNB) were determined for the novel variant enzyme purified via heterologous expression in Escherichia coli. Five substrates were used for measurement of specific activities, including isothiocyanate compounds that occur in cruciferous vegetables (benzylisothiocyanate, phenethylisothiocyanate, and sulforaphane). Such isothiocyanate substrates are potential cancer chemopreventive agents that are conjugated by GSTs. No major change in kinetic parameters was observed. However, the half-life at 50 °C of the Leu 151 enzyme was reduced to 12 min, as compared to 28 min for the Phe 151 enzyme. Residue 151 is located at the N-terminus of helix α6 in GST motif II, surrounded by hydrophobic residues, and near the conserved “hydrophobic staple” and N-capping box motifs. These local structural elements aid in formation of helix α6 and promote proper folding and protein stability. Analysis of the three-dimensional structure showed that substitution of Phe 151 with Leu produces a hydrophobic cavity in the GSTP1 core, thereby destabilizing its structure. Phe151Leu represents one of the first-described allelic variations in a protein folding motif.     

Structure-function analysis of Arg-Gly-Asp helix motifs in alpha v beta 6 integrin ligands

Data relating to the structural basis of ligand recognition by integrins are limited. Here we describe the physical requirements for high affinity binding of ligands to alpha v beta6. By combining a series of structural analyses with functional testing, we show that 20-mer peptide ligands, derived from high affinity ligands of alpha v beta6 (foot-and-mouth-disease virus, latency associated peptide), have a common structure comprising an Arg-Gly-Asp motif at the tip of a hairpin turn followed immediately by a C-terminal helix. This arrangement allows two conserved Leu/Ile residues at Asp(+1) and Asp(+4) to be presented on the outside face of the helix enabling a potential hydrophobic interaction with the alpha v beta6 integrin, in addition to the Arg-Gly-Asp interaction. The extent of the helix determines peptide affinity for alpha v beta6 and potency as an alpha v beta6 antagonist. A major role of this C-terminal helix is likely to be the correct positioning of the Asp(+1) and Asp(+4) residues. These data suggest an explanation for several biological functions of alpha v beta6 and provide a structural platform for design of alpha v beta6 antagonists.     

Activation of rat glutathione transferases in class mu by active oxygen species

The activities of rat glutathione transferases (GSTs) 3-3, 3-4, 4-4 in Class mu towards 1-chloro-2,4-dinitrobenzene (CDNB) but not 1,2-dichloro-4-nitrobenzene were increased up to 5-fold during preincubation with 0.4 mM xanthine and xanthine oxidase in 50 mM potassium phosphate, pH 7.8, containing 0.1 mM EDTA. The activated GST 3-4, purified by S-hexylglutathione affinity chromatography after the treatment, had a higher specific activity (130 units/mg) than that of the nontreated (35 units/mg), the Km and Vmax values for glutathione or CDNB also were increased. Other rat GSTs in Class alpha and pi were inactivated by the same treatment. In the presence of superoxide dismutase, the activation of GST 3-4 did not occur.     

A conserved N-capping motif contributes significantly to the stabilization and dynamics of the C-terminal region of class Alpha glutathione S-transferases

Helix 9, the major structural element in the C-terminal region of class Alpha glutathione transferases, forms part of the active site of these enzymes where its dynamic properties modulate both catalytic and ligandin functions. A conserved aspartic acid N-capping motif for helix 9 was identified by sequence alignments of the C-terminal regions of class Alpha glutathione S-transferases (GSTs) and an analysis by the helix-coil algorithm AGADIR. The contribution of the N-capping motif to the stability and dynamics of the region was investigated by replacing the N-cap residue Asp-209 with a glycine in human glutathione S-transferase A1-1 (hGST A1-1) and in a peptide corresponding to its C-terminal region. Far-UV circular dichroism and AGADIR analyses indicate that, in the absence of tertiary interactions, the wild-type peptide displays a low intrinsic tendency to form a helix and that this tendency is reduced significantly by the Asp-to-Gly mutation. Disruption of the N-capping motif of helix 9 in hGST A1-1 alters the conformational dynamics of the C-terminal region and, consequently, the features of the H-site to which hydrophobic substrates (e.g. 1-chloro-2,4-dinitrobenzene (CDNB)) and nonsubstrates (e.g. 8-anilino-1-naphthalene sulfonate (ANS)) bind. Isothermal calorimetric and fluorescence data for complex formation between ANS and protein suggest that the D209G-induced perturbation in the C-terminal region prevents normal ligand-induced localization of the region at the active site, resulting in a less hydrophobic and more solvent-exposed H-site. Therefore, the catalytic efficiency of the enzyme with CDNB is diminished due to a lowered affinity for the electrophilic substrate and a lower stabilization of the transition state.     

Microtubule disassembly and morphologic alterations induced by 1-chloro-2,4-dinitrobenzene, a substrate for glutathione S-transferase

1-chloro-2,4-dinitrobenzene (CDNB), a potent substrate for glutathione S-transferase, is known to rapidly deplete cellular glutathione (GSH) via conjugate formation. Treatment of quiescent 3T3 cells with 5 uM CDNB results in disassembly of microtubules (MT) within 1 hr as revealed by indirect immunofluorescence microscopy. In addition, CDNB treatment also induces dramatic morphologic alterations similar to those mediated by colchicine. Furthermore, taxol prevents both MT disassembly and morphologic changes normally occurring in CDNB as well as colchicine-treated cells. The mechanism of CDNB-mediated MT disassembly in vivo and its possible relationship to cellular GSH metabolism are under current studies.     

Induction of anionic glutathione transferases in rat liver by propylthiouracil

The treatment of rats with propylthiouracil (PTU) resulted in an induction of not only cationic but also anionic glutathione (GSH) transferases. DE-52 column chromatography of the anionic GSH transferases divided into five main peaks (I-V). Peaks I-IV had a high activity toward 1-chloro-2,4-dinitrobenzene (CDNB). Peak V, which is a new form of the enzyme, showed high activity to ethacrynic acid. PTU induced peaks I and V, whereas phenobarbital increased the activity of only peak I.     

The action of the glutathione transferase substrate, 1-chloro-2,4-dinitrobenzene on synaptosomal glutathione content and the release of hydrogen peroxide

We studied the action of the glutathione transferase substrate, 1-chloro-2,4-dinitrobenzene (CDNB) on the synaptosomal production of H2O2. We found that CDNB (30-40 microM) readily depletes the cytosolic glutathione but is almost without effect on the mitochondrial fraction. The depletion of the cytosolic glutathione induced by CDNB affords the detection in the extracellular space of H2O2 produced intrasynaptosomally upon increasing the cytosolic Ca2+ concentration that is otherwise destroyed by glutathione peroxidase. Higher concentrations of CDNB induce a H2O2 production which is not related to the glutathione content. This H2O2 is of mitochondrial origin and requires that NAD be reduced. The primary product of the mitochondrial CD-NB-dependent oxygen reduction is at least in part the superoxide anion.     

Inhibition of glutathione transferase pi from human placenta by 1-chloro-2,4-dinitrobenzene occurs because of covalent reaction with cysteine 47.

Human placenta glutathione transferase pi is irreversibly inhibited when incubated with 1-chloro-2,4-dinitrobenzene (CDNB) in the absence of the cosubstrate glutathione. The enzyme is protected against CDNB inactivation by the presence of S-methylglutathione and glutathione. The kinetics of inactivation is pseudo-first-order with k(obs) = 0.04 min-1 when 44 microM enzyme is incubated in presence of 1 mM CDNB at pH 6.5. The inhibition is saturable with respect to the CDNB concentration and the enzyme-CDNB complex shows a K(i) = 2.7 mM. Concomitant to the inhibition process is formation of an absorption band at 340 nm. After trypsin digestion and HPLC analysis, the CDNB-reacted enzyme gives a single peptide absorbing at 340 nm. Automated Edman degradation of this peptide indicates cysteine 47 to be the residue alkylated by CDNB.     

Structural changes in the C-terminus of Ca2+-bound rat S100B (beta beta) upon binding to a peptide derived from the C-terminal regulatory domain of p53.

S100B(beta beta) is a dimeric Ca2+-binding protein that interacts with p53, inhibits its phosphorylation by protein kinase C (PKC) and promotes disassembly of the p53 tetramer. Likewise, a 22 residue peptide derived from the C-terminal regulatory domain of p53 has been shown to interact with S100B(beta beta) in a Ca2+-dependent manner and inhibits its phosphorylation by PKC. Hence, structural studies of Ca2+-loaded S100B(beta beta) bound to the p53 peptide were initiated to characterize this interaction. Analysis of nuclear Overhauser effect (NOE) correlations, amide proton exchange rates, 3J(NH-H alpha) coupling constants, and chemical shift index data show that, like apo- and Ca2+-bound S100B(beta beta), S100B remains a dimer in the p53 peptide complex, and each subunit has four helices (helix 1, Glu2-Arg20; helix 2, Lys29-Asn38; helix 3, Gln50-Asp61; helix 4, Phe70-Phe87), four loops (loop 1, Glu21-His25; loop 2, Glu39-Glu49; loop 3, Glu62-Gly66; loop 4, Phe88-Glu91), and two beta-strands (beta-strand 1, Lys26-Lys28; beta-strand 2, Glu67-Asp69), which forms a short antiparallel beta-sheet. However, in the presence of the p53 peptide helix 4 is longer by five residues than in apo- or Ca2+-bound S100B(beta beta). Furthermore, the amide proton exchange rates in helix 3 (K55, V56, E58, T59, L60, D61) are significantly slower than those of Ca2+-bound S100B(beta beta). Together, these observations plus intermolecular NOE correlations between the p53 peptide and S100B(beta beta) support the notion that the p53 peptide binds in a region of S100B(beta beta), which includes residues in helix 2, helix 3, loop 2, and the C-terminal loop, and that binding of the p53 peptide interacts with and induces the extension of helix 4.     

一个新的产黄青霉谷胱甘肽转移酶基因的克隆和鉴定

以产黄青霉(Penicillium chrysogenum Thom)cDNA为模板,克隆得到一个新的谷胱甘肽转移酶基因PcgstB,其开放阅读框长651bp,编码216个氨基酸的蛋白质。与已知序列进行BLASTp比较显示,该蛋白具有保守的GST结构域,与烟曲霉GstB的序列一致性最高,达65%。将PcgstB与原核表达载体pTrc99A连接得到表达质粒pTrc-gstB,转化大肠杆菌DH5α,经IPTG诱导后获得以可溶形式表达的重组PcGstB蛋白。以1-chloro-2,4-dinitrobenzene(CDNB)为底物检测,确认该蛋白具有GST活性。     

Molecular determinants of xenobiotic metabolism: QM/MM simulation of the conversion of 1-chloro-2,4-dinitrobenzene catalyzed by M1-1 glutathione S-transferase

Modeling methods allow the identification and analysis of determinants of reactivity and specificity in enzymes. The reaction between glutathione and 1-chloro-2,4-dinitrobenzene (CDNB) is widely used as a standard activity assay for glutathione S-transferases (GSTs). It is important to understand the causes of differences between catalytic GST isoenzymes and the effects of mutations and genetic polymorphisms. Quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations have been performed here to investigate the addition of the glutathione anion to CDNB in the wild-type M1-1 GST isoenzyme from rat and in three single point mutant (Tyr6Phe, Tyr115Phe, and Met108Ala) M1-1 GST enzymes. We have developed a specifically parameterized QM/MM method (AM1-SRP/CHARMM22) to model this reaction by fitting to experimental heats of formation and ionization potentials. Free energy profiles were obtained from molecular dynamics simulations of the reaction using umbrella sampling and weighted histogram analysis techniques. The reaction in solution has also been simulated and is compared to the enzymatic reaction. The free energies are in excellent agreement with experimental results. Overall the results of the present study show that QM/MM reaction pathway analysis provides detailed insight into the chemistry of GST and can be used to obtain mechanistic insight into the effects of specific mutations on this catalytic process.