Structural determinants of glutathione transferases with azathioprine activity identified by DNA shuffling of alpha class members

A library of alpha class glutathione transferases (GSTs), composed of chimeric enzymes derived from human (A1-1, A2-2 and A3-3), bovine (A1-1) and rat (A2-2 and A3-3) cDNA sequences was constructed by the method of DNA shuffling. The GST variants were screened in bacterial lysates for activity with the immunosuppressive agent azathioprine, a prodrug that is transformed into its active form, 6-mercaptopurine, by reaction with the tripeptide glutathione catalyzed by GSTs. Important structural determinants for activity with azathioprine were recognized by means of primary structure analysis and activities of purified enzymes chosen from the screening. The amino acid sequences could be divided into 23 exchangeable segments on the basis of the primary structures of 45 chosen clones. Segments 2, 20, 21, and 22 were identified as primary determinants of the azathioprine activity representing two of the regions forming the substrate-binding H-site. Segments 21 and 22 are situated in the C-terminal helix characterizing alpha class GSTs, which is instrumental in their catalytic function. The study demonstrates the power of DNA shuffling in identifying segments of primary structure that are important for catalytic activity with a targeted substrate. GSTs in combination with azathioprine have potential as selectable markers for use in gene therapy. Knowledge of activity-determining segments in the structure is valuable in the protein engineering of glutathione transferase for enhanced or suppressed activity.     

Polymorphisms of glutathione S-transferases M1, T1, P1 and A1 genes in the Tunisian population: an intra and interethnic comparative approach

Genetic polymorphisms in glutathione S-transferases (GSTs) genes might influence the detoxification activities of the enzymes predisposing individuals to cancer risk. Owing to the presence of these genetic variants, inter-individual and ethnic differences in GSTs detoxification capacity have been observed in various populations. Therefore, the present study was performed to determine the prevalence GSTM1 0/0, GSTT1 0/0, GSTP1 Ile(105)Val, and GSTA1 A/B polymorphisms in 154 healthy individuals from South Tunisia, and to compare them with those observed in North and Centre Tunisian populations and other ethnic groups. GSTM1 and GSTT1 polymorphisms were analyzed by a Multiplex-PCR approach, whereas GSTP1 and GSTA1 polymorphisms were examined by PCR-RFLP. The frequencies of GSTM10/0 and GSTT1 0/0 genotypes were 53.9% and 27.9%, respectively. The genotype distribution of GSTP1 was 47.4% (Ile/Ile), 40.9% (Ile/Val), and 11.7% (Val/Val). For GSTA1, the genotype distribution was 24.7% (A/A), 53.9% (A/B), and 21.4% (B/B). The combined genotypes distribution of GSTM1, GSTT1, GSTP1 and GSTA1 polymorphisms showed that thirty one of the 36 possible genotypes were present in our population; eight of them have a frequency greater than 5%. To the best of our knowledge, this is the first report of GSTs polymorphisms in South Tunisian population. Our findings demonstrate the impact of ethnicity and reveal a characteristic pattern for Tunisian population. The molecular studies in these enzymes provide basis for further epidemiological investigations in the population where these functional polymorphisms alter therapeutic response and act as susceptibility markers for various clinical conditions.     

Metal binding ability of glutathione transferases conserved between two animal species,the vanadium-rich ascidian Ascidia sydneiensis samea and the schistosome Schistosoma japonicum

Glutathione transferases (GSTs) are multifunctional enzymes found in many organisms. We recently identified vanadium-binding GSTs, designated AsGSTs, from the vanadium-rich ascidian, Ascidia sydneiensis samea. In this study, the metal-selectivity of AsGST-I was investigated. Immobilized metal ion affinity chromatography (IMAC) analysis revealed that AsGST-I binds to V(IV), Fe(III), and Cu(II) with high affinity in the following order Cu(II) > V(IV) > Fe(III), and to Co(II), Ni(II), and Zn(II) with low affinity. The GST activity of AsGST-I was inhibited dose-dependently by not V(IV) but Cu(II). A competition experiment demonstrated that the binding of V(IV) to AsGST-I was not inhibited by Cu(II). These results suggest that AsGST-I has high V(IV)-selectivity, which can confer the specific vanadium accumulation of ascidians. Because there are few reports on the metal-binding ability of GSTs, we performed the same analysis on SjGST (GST from the schistosome, Schistosoma japonicum). SjGST also demonstrated metal-binding ability although the binding pattern differed from that of AsGST-I. The GST activity of SjGST was inhibited by Cu(II) only, as that of AsGST-I. Our results indicate a possibility that metal-binding abilities of GSTs are conserved among organisms, at least animals, which is suggestive of a new role for these enzymes in metal homeostasis or detoxification.     

Structural Basis for Featuring of Steroid Isomerase Activity in Alpha Class Glutathione Transferases

Glutathione transferases (GSTs) are abundant enzymes catalyzing the conjugation of hydrophobic toxic substrates with glutathione. In addition to detoxication, human GST A3-3 displays prominent steroid double-bond isomerase activity; e.g. transforming Δ⁵-androstene-3-17-dione into Δ⁴-androstene-3-17-dione (AD). This chemical transformation is a crucial step in the biosynthesis of steroids, such as testosterone and progesterone. In contrast to GST A3-3, the homologous GST A2-2 does not show significant steroid isomerase activity. We have solved the 3D structures of human GSTs A2-2 and A3-3 in complex with AD. In the GST A3-3 crystal structure, AD was bound in an orientation suitable for the glutathione (GSH)-mediated catalysis to occur. In GST A2-2, however, AD was bound in a completely different orientation with its reactive double bond distant from the GSH-binding site. The structures illustrate how a few amino acid substitutions in the active site spectacularly alter the binding mode of the steroid substrate in relation to the conserved catalytic groups and an essentially fixed polypeptide chain conformation. Furthermore, AD did not bind to the GST A2-2-GSH complex. Altogether, these results provide a first-time structural insight into the steroid isomerase activity of any GST and explain the 5000-fold difference in catalytic efficiency between GSTs A2-2 and A3-3. More generally, the structures illustrate how dramatic diversification of functional properties can arise via minimal structural alterations. We suggest a novel structure-based mechanism of the steroid isomerization reaction.     

A Novel Quasi-Species of Glutathione Transferase with High Activity towards Naturally Occurring Isothiocyanates Evolves from Promiscuous Low-Activity Variants

Glutathione transferases (GSTs) are known as promiscuous enzymes capable of catalyzing the conjugation of glutathione with a broad range of electrophilic substrates. A previous study based on recombinant chimeras derived from human GST M1-1 and GST M2-2 demonstrated the formation of a subset of F1 generation GSTs, which had lost high activity with substrates distinguishing parental enzymes. In the present study, the members of this subset were recombined by DNA shuffling to produce an F2 generation of GSTs. Screening of 930 bacterial clones demonstrated that 83% of recombinant enzyme variants were active with at least one of three alternative substrates: phenethyl isothiocyanate (PEITC), 1-chloro-2,4-dinitrobenzene, or p-nitrophenyl acetate. The majority had similar low activity as the parental GSTs in the F1 generation. However, 17 novel enzymes displayed high activity with PEITC. Half of these enzymes were similar to GST M1-1, which also has high activity with the same substrate, and all of these GSTs featured Tyr116/Ser210 in the active site. This group of F2 variants apparently had reverted to the GST M1-1 type. A second group of F2 variants with high PEITC activity was characterized by His116 in the active site. This category represented a new variety of GSTs, which demonstrated higher selectivity for isothiocyanate substrates than the GST M1-1 type. The different groups of GSTs can be considered as distinct molecular quasi-species, each of which comprises variant amino acid sequences. The quasi-species are structurally distinguished by active-site residues that govern their substrate selectivities. Clearly, minimal alterations of the active site can generate enzymes with highly distinctive functional properties.     

Comparative genomics identifies new alpha class genes within the avian glutathione S-transferase gene cluster

Glutathione S-transferases (GSTs: EC2.5.1.18) are a superfamily of multifunctional dimeric enzymes that catalyze the conjugation of glutathione (GSH) to electrophilic chemicals. In most animals and in humans, GSTs are the principal enzymes responsible for detoxifying the mycotoxin aflatoxin B₁ (AFB₁) and GST dysfunction is a known risk factor for susceptibility towards AFB₁. Turkeys are one of the most susceptible animals known to AFB₁, which is a common contaminant of poultry feeds. The extreme susceptibility of turkeys is associated with hepatic GSTs unable to detoxify the highly reactive and electrophilic metabolite exo-AFB₁-8,9-epoxide (AFBO). In this study, comparative genomic approaches were used to amplify and identify the α-class tGST genes (tGSTA1.1, tGSTA1.2, tGSTA1.3, tGSTA2, tGSTA3 and tGSTA4) from turkey liver. The conserved GST domains and four α-class signature motifs in turkey GSTs (with the exception of tGSTA1.1 which lacked one motif) confirm the presence of hepatic α-class GSTs in the turkey. Four signature motifs and conserved residues found in α-class tGSTs are (1) xMExxxWLLAAAGVE, (2) YGKDxKERAxIDMYVxG, (3) PVxEKVLKxHGxxxL and (4) PxIKKFLXPGSxxKPxxx. A BAC clone containing the α-class GST gene cluster was isolated and sequenced. The turkey α-class GTS genes genetically map to chromosome MGA2 with synteny between turkey and human α-class GSTs and flanking genes. This study identifies the α-class tGST gene cluster and genetic markers (SNPs, single nucleotide polymorphisms) that can be used to further examine AFB₁ susceptibility and resistance in turkeys. Functional characterization of heterologously expressed proteins from these genes is currently underway.     

Identification and cloning of a key insecticide-metabolizing glutathione S-transferase (MdGST-6A) from a hyper insecticide-resistant strain of the housefly Musca domestica

Strains of the housefly, Musca domestica, highly resistant to organophosphate (OP) and other insecticides are known because they overproduce glutathione S-transferases (GSTs). Previous work has shown that overproduction in these strains involved numerous isozymes with glutathione conjugating activities (Pesticide Biochem. Physiol., 25 (1986) 169; Mol. General Genetics, 227 (1991) 355; J. Biol. Chem., 267 (1992) 1840; Mol. General Genetics, 245 (1994) 236; J. Mol. Evol., 43 (1996) 236). The current work describes the purification and identification of a M. domestica GST isozyme (pI 7.1) broadly specific for substrates from a housefly strain, Cornell-HR, that is highly resistant against OP-insecticides, and the isolation of two new MdGST genes using the antibody made against it. This isozyme, which was identified from amongst more than 20 isoelectric forms of GSTs of the same subunit size, was highly active for conjugating GSH to the model substrate 3,4-dichloronitrobenzne (DCNB). When expressed in Escherichia coli, one of the cloned GSTs, MdGST-6A, produces an enzyme that conjugates glutathione to the insecticides methyl parathion and lindane. On indication that it was the most active isozyme toward several xenobiotics among several MdGSTs tested, we advance the notion that MdGST-6A probably plays an important role in M. domestica Cornell-HR's resistance towards OP-insecticides. MdGST-6A and a second closely related one found in this work, MdGST-6B, are members of the traditional insect class I family (theta-class) and share the greatest homologies with a cluster of Drosophila GSTs on locus 55. In addition to having the unusually broad substrate specificity, the sequence of the new group of enzymes reveals that it has a highly diverged hydrophobic motif in its active site as compared to other class I GSTs from insects.     

Phytodetoxification of the environmental pollutant and explosive 2,4,6-trinitrotoluene

Our recent study highlights the role of 2 glutathione transferases (GSTs) in the detoxification of the environmental pollutant, 2,4,6-trinitrotoluene (TNT) in Arabidopsis thaliana. TNT is toxic and highly resistant to biodegradation in the environment, raising both health and environmental concerns. Two GSTs, GST-U24 and GST-U25, are upregulated in response to TNT treatment, and expressed predominantly in the root tissues; the site of TNT location following uptake. Plants overexpressing GST-U24 and GST-U25 exhibited significantly enhanced ability to withstand and detoxify TNT, and remove TNT from contaminated soil. Analysis of the catalytic activities of these 2 enzymes revealed that they form 3 TNT-glutathionyl products. Of particular interest is 2-glutathionyl-4,6-dinitrotoluene as this represents a potentially favorable step toward subsequent degradation and mineralization of TNT. We demonstrate how GSTs fit into what is already known about pathways for TNT detoxification, and discuss the short and longer-term fate of TNT conjugates in planta.     

Catalytic properties of glutathione-binding residues in a tau class glutathione transferase (PtGSTU1) from Pinus tabulaeformis

Glutathione transferases (GSTs) play important roles in stress tolerance and detoxification in plants. However, there is extremely little information on the molecular characteristics of GSTs in gymnosperms. In a previous study, we cloned a tau class GST (PtGSTU1) from a gymnosperm (Pinus tabulaeformis) for the first time. Based on the N-terminal amino acid sequence identity to the available crystal structures of plant tau GSTs, Ser13, Lys40, Ile54, Glu66 and Ser67 of PtGSTU1 were proposed as glutathione-binding (G-site) residues. The importance of Ser13 as a G-site residue was investigated previously. The functions of Lys40, Ile54, Glu66 and Ser67 of PtGSTU1 are examined in this study through site-directed mutagenesis. Enzyme assays and thermal stability measurements on the purified recombinant PtGSTU1 showed that substitution at each of these sites significantly affects the enzyme's substrate specificity and affinity for GSH, and these residues are essential for maintaining the stability of PtGSTU1. The results of protein expression and refolding analyses suggest that Ile54 is involved in the protein folding process. The findings demonstrate that the aforementioned residues are critical components of active sites that contribute to the enzyme's catalytic activity and structural stability.     

Emergence of novel enzyme quasi-species depends on the substrate matrix

Current research on enzyme evolution has shown that many enzymes are promiscuous and have activities with alternative substrates. Mutagenesis tends to relax substrate selectivity, and evolving enzymes can be regarded (summed over evolutionary time) as clusters of enzyme variants, or “quasi-species,” tested against a “substrate matrix” defined by all chemical substances to which the evolvants are exposed.In this investigation, the importance of the substrate matrix for identification of evolvable clusters of enzymes was evaluated by random sampling of variants from a library of glutathione transferase (GST) mutants. The variant GSTs were created by DNA shuffling of homologous Alpha class sequences. The substrate matrix was an array of alternative substrates used under defined experimental conditions. The measured enzyme activities produced a rectangular matrix, in which the rows can be projected as enzyme vectors in substrate-activity space and, reciprocally, the columns can be projected as alternative substrate vectors in enzyme-activity space. Multivariate analysis of the catalytic activities demonstrated that the enzyme vectors formed two primary clusters or functional “molecular quasi-species.” These quasi-species serve as the raw material from which more specialized enzymes eventually could evolve. The substrate vectors similarly formed two major groups. Identification of separate quasi-species of GSTs in a mutant library was critically dependent on the nature of the substrate matrix. When substrates from just one of the two groups were used, only one cluster of enzymes could be recognized. On the other hand, expansion of the substrate matrix to include additional substrates showed the presence of a third quasi-species among the GST variants already analyzed. Thus, the portrayal of the functional quasi-species is intimately linked to the effective substrate matrix. In natural evolution, the substrates actually encountered therefore play a pivotal role in determining whether latent catalytic abilities become manifest in novel enzymes.     

The Impact of Nitric Oxide Toxicity on the Evolution of the Glutathione Transferase Superfamily: A PROPOSAL FOR AN EVOLUTIONARY DRIVING FORCE*

Glutathione transferases (GSTs) are protection enzymes capable of conjugating glutathione (GSH) to toxic compounds. During evolution an important catalytic cysteine residue involved in GSH activation was replaced by serine or, more recently, by tyrosine. The utility of these replacements represents an enigma because they yield no improvements in the affinity toward GSH or in its reactivity. Here we show that these changes better protect the cell from nitric oxide (NO) insults. In fact the dinitrosyl·diglutathionyl·iron complex (DNDGIC), which is formed spontaneously when NO enters the cell, is highly toxic when free in solution but completely harmless when bound to GSTs. By examining 42 different GSTs we discovered that only the more recently evolved Tyr-based GSTs display enough affinity for DNDGIC (K_D < 10⁻⁹ m) to sequester the complex efficiently. Ser-based GSTs and Cys-based GSTs show affinities 10²–10⁴ times lower, not sufficient for this purpose. The NO sensitivity of bacteria that express only Cys-based GSTs could be related to the low or null affinity of their GSTs for DNDGIC. GSTs with the highest affinity (Tyr-based GSTs) are also over-represented in the perinuclear region of mammalian cells, possibly for nucleus protection. On the basis of these results we propose that GST evolution in higher organisms could be linked to the defense against NO.     

Cooperative kinetics of the recombinant glutathione transferase of Taenia solium and characterization of the enzyme

Glutathione transferases (GSTs) are essential enzymes in many organisms due their diverse functions and, in helminths they are the main detoxification system. For Taenia solium, two cytosolic GSTs with molecular masses of 25.5 and 26.5 kDa (Ts26GST) have been found. Ts26GST was cloned to be studied in its recombinant form (recTs26GST). Although the primary structure is related to the mu class, the kinetic parameters for CDNB (V_max = 51.5 μmol min⁻¹ mg⁻¹; K_m = 1.06 mM; k_cat = 22.2 s⁻¹) are related with some alpha GSTs. The substrate and inhibitor class markers reaffirmed these bimodal characteristics. Inhibition studies with anthelminthics indicate that recTs26GST is sensitive to mebendazole, displaying a non competitive inhibition pattern suggesting that at least two molecules are binding to recTs26GST. On the other hand, the kinetic curves for CDNB and GSH showed a positive cooperativity that was corroborated using fluorometric assays. Those assays indicate that CDNB binding is highly influenced by GSH, probably by modulation of the Ts26GST conformational ensamble.     

Heterologous expression and functional characterization of avian mu-class glutathione S-transferases

Hepatic glutathione S-transferases (GSTs: EC2.5.1.1.8) catalyze the detoxification of reactive electrophilic compounds, many of which are toxic and carcinogenic intermediates, via conjugation with the endogenous tripeptide glutathione (GSH). Glutathione S-transferase (GST)-mediated detoxification is a critical determinant of species susceptibility to the toxic and carcinogenic mycotoxin aflatoxin B₁ (AFB₁), which in resistant animals efficiently detoxifies the toxic intermediate produced by hepatic cytochrome P450 bioactivation, the exo-AFB₁-8,9-epoxide (AFBO). Domestic turkeys (Meleagris gallopavo) are one of the most sensitive animals known to AFB₁, a condition associated with a deficiency of hepatic GST-mediated detoxification of AFBO. We have recently shown that unlike their domestic counterparts, wild turkeys (Meleagris gallopavo silvestris), which are relatively resistant, express hepatic GST-mediated detoxification activity toward AFBO. Because of the importance of GSTs in species susceptibility, and to explore possible GST classes involved in AFB₁ detoxification, we amplified, cloned, expressed and functionally characterized the hepatic mu-class GSTs tGSTM3 (GenBank accession no. JF340152), tGSTM4 (JF340153) from domestic turkeys, and a GSTM4 variant (ewGSTM4, JF340154) from Eastern wild turkeys. Predicted molecular masses of tGSTM3 and two tGSTM4 variants were 25.6 and 25.8 kDa, respectively. Multiple sequence comparisons revealed four GSTM motifs and the mu-loop in both proteins. tGSTM4 has 89% amino acid sequence identity to chicken GSTM2, while tGSTM3 has 73% sequence identity to human GSTM3 (hGSTM3). Specific activities of Escherichia coli-expressed tGSTM3 toward 1-chloro-2,4-dinitrobenzene (CDNB) and peroxidase activity toward cumene hydroperoxide were five-fold greater than tGSTM4 while tGSTM4 possessed more than three-fold greater activity toward 1,2-dichloro-4-nitrobenzene (DCNB). The two enzymes displayed equal activity toward ethacrynic acid (ECA). However, none of the GSTM proteins had AFBO detoxification capability, in contrast to recombinant alpha-class GSTs shown in our recent study to possess this important activity. In total, our data indicate that although turkey hepatic GSTMs may contribute to xenobiotic detoxification, they probably play no role in detoxification of AFBO in the liver.     

Detection of anti-oxidant enzymatic activities and purification of glutathione transferases from Angiostrongylus cantonensis

There are several anti-oxidant enzyme families that play pivotal roles in facilitating the survival of parasites. Glutathione transferases (GSTs) are members of the anti-oxidant family that can detoxify a broad range of exogenous or endogenous compounds including reactive oxidative species. GSTs have been studied as vaccine candidates, immunodiagnostic markers and as treatment targets. Helminths of the genus Angiostrongylus live inside arteries of vertebrates and two main species are associated with accidental human infections: Angiostrongylus costaricensis adult worms live inside the mesenteric arteries and larvae of Angiostrongylus cantonensis become trapped in the central nervous system vasculature. Since the interactions between angiostrongylid nematodes and their vertebrate hosts are poorly understood, this study characterized the anti-oxidant enzymatic activities of A. cantonensis from female worms by collecting excreted and secreted (ES) and total extract (TE) molecules. Catalase (CAT) and superoxide dismutase (SOD) activities were found both in the ES and TE while glutathione peroxidase (GPX) and GST were found only in the TE. GSTs were purified by glutathione agarose affinity column (AcGST) and the pool of eluted GSTs was analyzed by mass spectrometry (LC-MS/MS) and de novo sequencing (Masslynx software). Sequences from two peptides (AcGSTpep1 and AcGSTpep2) present high identity to the N-terminal and C-terminal from sigma class GSTs of nematodes. It is known that these GST enzymes are associated with host immune regulation. Furthermore, understanding the role of parasite-derived anti-oxidant molecules is important in understanding host-parasite interactions.     

Plasmodium vivax: molecular cloning, expression and characterization of glutathione S-transferase

Malaria parasite glutathione S-transferases (GSTs) are postulated to be essential for parasite survival by protecting the parasite against oxidative stress and buffering the detoxification of heme-binding compounds; therefore, GSTs are considered potential targets for drug development. In this study, we identified a Plasmodium vivax gene encoding GST (PvGST) and characterized the biochemical properties of the recombinant enzyme. The PvGST contained 618 bp that encoded 205 amino acids and shared a significant degree of sequence identity with GSTs from other Plasmodium species. The recombinant homodimeric enzyme had an approximate molecular mass of 50kDa and exhibited GSH-conjugating and GSH-peroxidase activities towards various model substrates. The optimal pH for recombinant PvGST (rPvGST) activity was pH 8.0, and the enzyme was moderately unstable at 37 degrees C. The K(m) values of rPvGST with respect to GSH and CDNB were 0.17+/-0.09 and 2.1+/-0.4mM, respectively. The significant sequence homology and similar biochemical properties of PvGST and Plasmodium falciparum GST (PfGST) indicate that they may have similar molecular structures. This information may be useful for the design of specific inhibitors for plasmodial GSTs as potential antimalarial drugs.     

Identification and cloning of a key insecticide-metabolizing glutathione S-transferase (MdGST-6A) from a hyper insecticide-resistant strain of the housefly Musca domestica

Strains of the housefly, Musca domestica, highly resistant to organophosphate (OP) and other insecticides are known because they overproduce glutathione S-transferases (GSTs). Previous work has shown that overproduction in these strains involved numerous isozymes with glutathione conjugating activities (Pesticide Biochem. Physiol., 25 (1986) 169; Mol. General Genetics, 227 (1991) 355; J. Biol. Chem., 267 (1992) 1840; Mol. General Genetics, 245 (1994) 236; J. Mol. Evol., 43 (1996) 236). The current work describes the purification and identification of a M. domestica GST isozyme (pI 7.1) broadly specific for substrates from a housefly strain, Cornell-HR, that is highly resistant against OP-insecticides, and the isolation of two new MdGST genes using the antibody made against it. This isozyme, which was identified from amongst more than 20 isoelectric forms of GSTs of the same subunit size, was highly active for conjugating GSH to the model substrate 3,4-dichloronitrobenzne (DCNB). When expressed in Escherichia coli, one of the cloned GSTs, MdGST-6A, produces an enzyme that conjugates glutathione to the insecticides methyl parathion and lindane. On indication that it was the most active isozyme toward several xenobiotics among several MdGSTs tested, we advance the notion that MdGST-6A probably plays an important role in M. domestica Cornell-HR's resistance towards OP-insecticides. MdGST-6A and a second closely related one found in this work, MdGST-6B, are members of the traditional insect class I family (theta-class) and share the greatest homologies with a cluster of Drosophila GSTs on locus 55. In addition to having the unusually broad substrate specificity, the sequence of the new group of enzymes reveals that it has a highly diverged hydrophobic motif in its active site as compared to other class I GSTs from insects.     

Stereochemical Features of Glutathione-dependent Enzymes in the Sphingobium sp. Strain SYK-6 β-Aryl Etherase Pathway

Glutathione-dependent enzymes play important protective, repair, or metabolic roles in cells. In particular, enzymes in the glutathione S-transferase (GST) superfamily function in stress responses, defense systems, or xenobiotic detoxification. Here, we identify novel features of bacterial GSTs that cleave β-aryl ether bonds typically found in plant lignin. Our data reveal several original features of the reaction cycle of these GSTs, including stereospecific substrate recognition and stereoselective formation of β-S-thioether linkages. Products of recombinant GSTs (LigE, LigP, and LigF) are β-S-glutathionyl-α-keto-thioethers that are degraded by a β-S-thioetherase (LigG). All three Lig GSTs produced the ketone product (β-S-glutathionyl-α-veratrylethanone) from an achiral side chain-truncated model substrate (β-guaiacyl-α-veratrylethanone). However, when β-etherase assays were conducted with a racemic model substrate, β-guaiacyl-α-veratrylglycerone, LigE- or LigP-catalyzed reactions yielded only one of two potential product (β-S-glutathionyl-α-veratrylglycerone) epimers, whereas the other diastereomer (differing in configuration at the β-position (i.e. its β-epimer)) was produced only in the LigF-catalyzed reaction. Thus, β-etherase catalysis causes stereochemical inversion of the chiral center, converting a β(R)-substrate to a β(S)-product (LigE and LigP), and a β(S)-substrate to a β(R)-product (LigF). Further, LigG catalyzed glutathione-dependent β-S-thioether cleavage with β-S-glutathionyl-α-veratrylethanone and with β(R)-configured β-S-glutathionyl-α-veratrylglycerone but exhibited no or significantly reduced β-S-thioether-cleaving activity with the β(S)-epimer, demonstrating that LigG is a stereospecific β-thioetherase. We therefore propose that multiple Lig enzymes are needed in this β-aryl etherase pathway in order to cleave the racemic β-ether linkages that are present in the backbone of the lignin polymer.