Identification of residues in glutathione transferase capable of driving functional diversification in evolution. A novel approach to protein redesign

Integrated structural and evolutionary analysis reveals common mechanisms underlying adaptive evolution in mammals

Understanding the molecular basis of adaptation to the environment is a central question in evolutionary biology, yet linking detected signatures of positive selection to molecular mechanisms remains challenging. Here we demonstrate that combining sequence-based phylogenetic methods with structural information assists in making such mechanistic interpretations on a genomic scale. Our integrative analysis shows that positively selected sites tend to colocalize on protein structures and that positively selected clusters are found in functionally important regions of proteins, indicating that positive selection can contravene the well-known principle of evolutionary conservation of functionally important regions. This unexpected finding, along with our discovery that positive selection acts on structural clusters, opens previously unexplored strategies for the development of better models of protein evolution. Remarkably, proteins where we detect the strongest evidence of clustering belong to just two functional groups: Components of immune response and metabolic enzymes. This gives a coherent picture of pathogens and xenobiotics as important drivers of adaptive evolution of mammals.

Molecular Evolution of the Glutathione S-Transferase Family in the Bemisia tabaci Species Complex

The glutathione S-transferase (GST) family plays an important role in the adaptation of herbivorous insects to new host plants and other environmental constrains. The family codes for enzymes that neutralize reactive oxygen species and phytotoxins through the conjugation of reduced glutathione. Here, we studied the molecular evolution of the GST family in Bemisia tabaci, a complex of >35 sibling species, differing in their geographic and host ranges. We tested if some enzymes evolved different functionality, by comparing their sequences in six species, representing five of the six major genetic clades in the complex. Comparisons of the nonsynonymous to synonymous substitution ratios detected positive selection events in 11 codons of 5 cytosolic GSTs. Ten of them are located in the periphery of the GST dimer, suggesting a putative involvement in interactions with other proteins. Modeling the tertiary structure of orthologous enzymes, identified additional 19 mutations in 9 GSTs, likely affecting the enzymes' functionality. Most of the mutation events were found in the environmentally responsive classes Delta and Sigma, indicating a slightly different delta/sigma tool box in each species. At a broader genomic perspective, our analyses indicated a significant expansion of the Delta GST class in B. tabaci and a general association between the diet breadth of hemipteran species and their total number of GST genes. We raise the possibility that at least some of the identified changes improve the fitness of the B. tabaci species carrying them, leading to their better adaptation to specific environments.

Structure-based design and application of an engineered glutathione transferase for the development of an optical biosensor for pesticides determination

In the present work, a structure-based design approach was used for the generation of a novel variant of synthetic glutathione transferase (PvGmGSTU) with higher sensitivity towards pesticides. Molecular modelling studies revealed Phe117 as a key residue that contributes to the formation of the hydrophobic binding site (H-site) and modulates the affinity of the enzyme towards xenobiotic compounds. Site-saturation mutagenesis of position Phe117 created a library of PvGmGSTU variants with altered kinetic and binding properties. Screening of the library against twenty-five different pesticides, showed that the mutant enzyme Phe117Ile displays 3-fold higher catalytic efficiency and exhibits increased affinity towards α-endosulfan, compared to the wild-type enzyme. Based on these catalytic features the mutant enzyme Phe117Ile was explored for the development of an optical biosensor for α-endosulfan. The enzyme was entrapped in alkosixylane sol-gel system in the presence of two pH indicators (bromocresol purple and phenol red). The sensing signal was based on the inhibition of the sol-gel entrapped GST, with subsequent decrease of released [H+] by the catalytic reaction, measured by sol-gel entrapped indicators. The assay response at 562 nm was linear in the range pH = 4-7. Linear calibration curves were obtained for α-endosulfan in the range of 0-30 μΜ. The reproducibility of the assay response, expressed by relative standard deviation, was in the order of 4.1% (N = 28). The method was successfully applied to the determination of α-endosulfan in real water samples without sample preparation steps.

Rapid birth-death evolution and positive selection in detoxification-type glutathione S-transferases in mammals

Glutathione S-Transferases (GSTs) are phase II detoxification enzymes that may have evolved in response to changes of environmental substrates. GST genes formed a multigene family and in mammals, there are six classes known as Alpha, Mu, Omega, Pi, Theta, and Zeta. Recent studies in phase I detoxification system specifically the cytochrome P450s provided a general explanation on why genes from a common origin such as those in a multigene family have both phylogenetically stable and unstable genes. Genes that participate in core functions of organisms such as development and physiology are stable whereas genes that play a role in detoxification are unstable and evolve in a process known as birth-death evolution, which is characterised by frequent gene gains and losses. The generality of the birth-death model at explaining the evolution of detoxification enzymes beyond the phase I enzyme has not been comprehensively explored. This work utilized 383 Gst genes and 300 pseudogenes across 22 mammalian species to study gene gains and losses. GSTs vary greatly in their phylogenetic stability despite their overall sequence similarity. Stable Gst genes from Omega and Zeta classes do not show fluctuation in gene numbers from human to opossum. These genes play a role in biosynthesis related functions. Unstable genes that include Alpha, Mu, Pi and Theta undergo frequent gene gain and loss in a process known as birth-death evolution. Gene members of these four classes are well known for their roles in detoxification. Our positive selection screen identified five positively selected sites in mouse GSTA3. Previous studies showed two of these sites (108H and 208E) were biochemically tested as important residues that conferred catalytic activity against the toxic aflatoxin B1-8,9-epoxide. The functional significance against aflatoxin of the remaining three positively selected sites warrant further investigation.

Exploring sequence-function space of a poplar glutathione transferase using designed information-rich gene variants

Exploring the vicinity around a locus of a protein in sequence space may identify homologs with enhanced properties, which could become valuable in biotechnical and other applications. A rational approach to this pursuit is the use of ‘infologs’, i.e. synthetic sequences with specific substitutions capturing maximal sequence information derived from the evolutionary history of the protein family. Ninety-five such infolog genes of poplar glutathione transferase were synthesized and expressed in Escherichia coli, and the catalytic activities of the proteins determined with alternative substrates. Sequence–activity relationships derived from the infologs were used to design a second set of 47 infologs in which 90% of the members exceeded wild-type properties. Two mutants, C2 (V55I/E95D/D108E/A160V) and G5 (F13L/C70A/G122E), were further functionally characterized. The activities of the infologs with the alternative substrates 1-chloro-2,4-dinitrobenzene and phenethyl isothiocyanate, subject to different chemistries, were positively correlated, indicating that the examined mutations were affecting the overall catalytic competence without major shift in substrate discrimination. By contrast, the enhanced protein expressivity observed in many of the mutants were not similarly correlated with the activities. In conclusion, small libraries of well-defined infologs can be used to systematically explore sequence space to optimize proteins in multidimensional functional space.

Recent advances in protein engineering and biotechnological applications of glutathione transferases

Glutathione transferases (GSTs, EC 2.5.1.18) are a widespread family of enzymes that play a central role in the detoxification, metabolism, and transport or sequestration of endogenous or xenobiotic compounds. During the last two decades, delineation of the important structural and catalytic features of GSTs has laid the groundwork for engineering GSTs, involving both rational and random approaches, aiming to create new variants with new or altered properties. These approaches have expanded the usefulness of native GSTs, not only for understanding the fundamentals of molecular detoxification mechanisms, but also for the development medical, analytical, environmental, and agricultural applications. This review article attempts to summarize successful examples and current developments on GST engineering, highlighting in parallel the recent knowledge gained on their phylogenetic relationships, structural/catalytic features, and biotechnological applications.

Adaptive molecular evolution of the two-pore channel 1 gene TPC1 in the karst-adapted genus Primulina (Gesneriaceae)

BACKGROUND AND AIMS

Limestone karst areas possess high floral diversity and endemism. The genus Primulina, which contributes to the unique calcicole flora, has high species richness and exhibit specific soil-based habitat associations that are mainly distributed on calcareous karst soils. The adaptive molecular evolutionary mechanism of the genus to karst calcium-rich environments is still not well understood. The Ca²⁺-permeable channel TPC1 was used in this study to test whether its gene is involved in the local adaptation of Primulina to karst high-calcium soil environments.

METHODS

Specific amplification and sequencing primers were designed and used to amplify the full-length coding sequences of TPC1 from cDNA of 76 Primulina species. The sequence alignment without recombination and the corresponding reconstructed phylogeny tree were used in molecular evolutionary analyses at the nucleic acid level and amino acid level, respectively. Finally, the identified sites under positive selection were labelled on the predicted secondary structure of TPC1.

KEY RESULTS

Seventy-six full-length coding sequences of Primulina TPC1 were obtained. The length of the sequences varied between 2220 and 2286 bp and the insertion/deletion was located at the 5' end of the sequences. No signal of substitution saturation was detected in the sequences, while significant recombination breakpoints were detected. The molecular evolutionary analyses showed that TPC1 was dominated by purifying selection and the selective pressures were not significantly different among species lineages. However, significant signals of positive selection were detected at both TPC1 codon level and amino acid level, and five sites under positive selective pressure were identified by at least three different methods.

CONCLUSIONS

The Ca²⁺-permeable channel TPC1 may be involved in the local adaptation of Primulina to karst Ca²⁺-rich environments. Different species lineages suffered similar selective pressure associated with calcium in karst environments, and episodic diversifying selection at a few sites may play a major role in the molecular evolution of Primulina TPC1.

Enhancing the pharmaceutical properties of protein drugs by ancestral sequence reconstruction

Optimization of a protein's pharmaceutical properties is usually carried out by rational design and/or directed evolution. Here we test an alternative approach based on ancestral sequence reconstruction. Using available genomic sequence data on coagulation factor VIII and predictive models of molecular evolution, we engineer protein variants with improved activity, stability, and biosynthesis potential and reduced inhibition by anti-drug antibodies. In principle, this approach can be applied to any protein drug based on a conserved gene sequence.

Drosophila GSTs display outstanding catalytic efficiencies with the environmental pollutants 2,4,6-trinitrotoluene and 2,4-dinitrotoluene

The nitroaromatic explosive 2,4,6-trinitrotoluene (TNT) and the related 2,4-dinitrotoluene (DNT) are toxic environmental pollutants. The biotransformation and detoxication of these persistent compounds in higher organisms are of great significance from a health perspective as well as for the biotechnological challenge of bioremediation of contaminated soil. We demonstrate that different human glutathione transferases (GSTs) and GSTs from the fruit fly Drosophila melanogaster are catalysts of the biotransformation of TNT and DNT. The human GSTs had significant but modest catalytic activities with both DNT and TNT. However, D. melanogaster GSTE6 and GSTE7 displayed outstanding high activities with both substrates.

•

The explosive TNT is a carcinogenic environmental pollutant spread world-wide.

•

TNT and the related DNT can be detoxified by conjugation with cellular glutathione.

•

Previously studied plant glutathione transferases display modest activity with TNT.

•

We found that human GSTs from four classes have low activity with TNT and DNT.

•

By contrast Drosophila GSTE6 and GSTE7 displayed outstanding TNT and DNT activities.

Proteomic profiling of gill GSTs in Mytilus galloprovincialis from the North of Portugal and Galicia evidences variations at protein isoform level with a possible relation with water quality

Glutathione transferases (GSTs) are key for xenobiotic detoxification at the molecular level across phyla. These enzymes are therefore likely to be part of the defence mechanisms used by marine organisms, such as mussels, that thrive in highly polluted environments. Taking this hypothesis into account, we used proteomics to characterize the profile of GSTs from the gills of marine mussel Mytilus galloprovincialis in order to discriminate natural mussel populations exposed to different levels of pollution. Samples were collected between Cabo Home (Spain) and Matosinhos (Portugal) covering a north-south transect of approximately 122 Km of the Atlantic Ocean along the Western Coast of the Iberian Peninsula. GSTs from mussel gills were extracted and purified by affinity chromatography with glutathione as the binding substrate to the solid medium. We studied the abundance of GST isoforms by two-dimensional gel electrophoresis and matrix assisted laser desorption/ionization-time of flight mass spectrometry and assessed total activity. Eleven putative individual GSTs from classes Mu, Pi and Sigma were identified by proteomics. Few variations were observed in total GST activity of post-mitochondrial samples between sampling sites, with animals from Matosinhos (polluted site) showing highest GST activity and Cabo Home (clean site) showing lowest. This contrasts with the increased number of differences in the individual GST isoforms. Each mussel population showed unique GST proteomic profiles. Based on the results we conclude that proteomics surpasses the conventional GST enzymatic activity method to discriminate natural mussel populations and has potential application in environmental monitoring. It is reasonable to suggest that the GST proteomic profiles observed may reflect differences in contamination levels.

Structural and functional evolution of positively selected sites in pine glutathione S-transferase enzyme family

Phylogenetic analyses have identified positive selection as an important driver of protein evolution, both structural and functional. However, the lack of appropriate combined functional and structural assays has generally hindered attempts to elucidate patterns of positively selected sites and their effects on enzyme activity and substrate specificity. In this study we investigated the evolutionary divergence of the glutathione S-transferase (GST) family in Pinus tabuliformis, a pine that is widely distributed from northern to central China, including cold temperate and drought-stressed regions. GSTs play important roles in plant stress tolerance and detoxification. We cloned 44 GST genes from P. tabuliformis and found that 26 of the 44 belong to the largest (Tau) class of GSTs and are differentially expressed across tissues and developmental stages. Substitution models identified five positively selected sites in the Tau GSTs. To examine the functional significance of these positively selected sites, we applied protein structural modeling and site-directed mutagenesis. We found that four of the five positively selected sites significantly affect the enzyme activity and specificity; thus their variation broadens the GST family substrate spectrum. In addition, positive selection has mainly acted on secondary substrate binding sites or sites close to (but not directly at) the primary substrate binding site; thus their variation enables the acquisition of new catalytic functions without compromising the protein primary biochemical properties. Our study sheds light on selective aspects of the functional and structural divergence of the GST family in pine and other organisms.

Systematic identification and characterization of glutathione S-transferases in cynomolgus macaque

Glutathione S-transferases (GSTs) are essential drug-metabolizing enzymes, involved in conjugation of various endogenous and exogenous substrates. Cynomolgus macaque is an important primate species in drug metabolism studies; however, cynomolgus GSTs have not been fully characterized. In this study the cDNAs of 12 GSTs (GSTA3-A5, GSTK1, GSTM2-M4, GSTO2, GSTP1, GSTS1, and GSTT1/2) were isolated from cynomolgus macaque and rhesus macaque liver. Cynomolgus GSTM1 cDNA was not amplified and only an aberrantly spliced GSTM1 transcript was isolated from rhesus macaque. Amino acid sequences of these 12 GSTs shared high sequence identities (93-98%) and were clustered into the same clades as the human orthologs in the phylogenetic tree. The 12 GSTs had exon-intron structures similar to the human orthologs, and exhibited distinct tissue expression patterns. GSTA3, GSTA5, and GSTM3/O2 were expressed predominantly in adrenal gland, jejunum, and testis, respectively, whereas the other GSTs showed universal expression patterns in the 10 tissues analyzed. Comparison of expression levels showed that GSTA1, GSTK1, GSTA3, and GSTM3 were most abundantly expressed in liver/jejunum, kidney, adrenal gland, and testis, respectively. Metabolic assays of proteins expressed heterologously in Escherichia coli, showed that all 12 GSTs and 5 previously identified GSTs, GSTA1/2, GSTM5, GSTO1, and GSTZ1, catalyzed the conjugation of GST substrate(s) 1-chloro-2,4-dinitrobenzene and/or 1,2-epoxy-3-(p-nitrophenoxy)propane, indicating that these 17 GSTs are functional drug-metabolizing enzymes. These results suggest that the 12 GST genes examined in this study are expressed and encoded functional enzymes in cynomolgus macaque.

Characterization and functional analysis of four glutathione S-transferases from the migratory locust, Locusta migratoria

Glutathione S-transferases (GSTs) play an important role in detoxification of xenobiotics in both prokaryotic and eukaryotic cells. In this study, four GSTs (LmGSTd1, LmGSTs5, LmGSTt1, and LmGSTu1) representing different classes were identified from the migratory locust, Locusta migratoria. These four proteins were heterologously expressed in Escherichia coli as soluble fusion proteins, purified by Ni(2+)-nitrilotriacetic acid agarose column and biochemically characterized. LmGSTd1, LmGSTs5, and LmGSTu1 showed high activities with 1-chloro-2, 4-dinitrobenzene (CDNB), detectable activity with p-nitro-benzyl chloride (p-NBC) and 1, 2-dichloro-4-nitrobenzene (DCNB), whereas LmGSTt1 showed high activity with p-NBC and detectable activity with CDNB. The optimal pH of the locust GSTs ranged between 7.0 to 9.0. Ethacrynic acid and reactive blue effectively inhibited all four GSTs. LmGSTs5 was most sensitive to heavy metals (Cu(2+) and Cd(2+)). The maximum expression of the four GSTs was observed in Malpighian tubules and fat bodies as evaluated by western blot. The nymph mortalities after carbaryl treatment increased by 28 and 12% after LmGSTs5 and LmGSTu1 were silenced, respectively. The nymph mortalities after malathion and chlorpyrifos treatments increased by 26 and 18% after LmGSTs5 and LmGSTu1 were silenced, respectively. These results suggest that sigma GSTs in L. migratoria play a significant role in carbaryl detoxification, whereas some of other GSTs may also involve in the detoxification of carbaryl and chlorpyrifos.

cDNA cloning and mRNA expression of four glutathione S-transferase (GST) genes from Mytilus galloprovincialis

Glutathione S-transferases (GSTs) are phase II enzymes involved in the regulation of redox homeostasis and innate immune responses against bacterial infection, which also play important roles in the detoxification of xenobiotics. In this study, we reported four genes of the GST family (named MgGSTα, MgGSTS1, MgGSTS2, and MgGSTS3, respectively) from Mytilus galloprovincialis. MgGSTα, MgGSTS1, MgGSTS2, and MgGSTS3 consisted of open reading frame (ORF) of 648 bp, 612 bp, 621 bp and 609 bp respectively, which encoded proteins of 216, 204, 207 and 203 amino acids residues, respectively. Sequence analysis showed that the predicted protein sequence of MgGSTs contained the conserved domain of the GST_N and GST_C. Alignment analysis indicated that the MgGSTs were divided into two types, one was of alpha GST, and the others were of sigma class. Tissue distribution study revealed that MgGSTα, MgGSTS2, MgGSTS3 transcripts were highly expressed in hemocytes, while MgGSTS1 mRNA was most abundantly expressed in hepatopancreas. After bacterial challenge, the expression level of these MgGSTs in hemocytes were all significantly higher than that of the control group. These results suggested that MgGSTs might play important roles in the modulation of immune response in M. galloprovincialis.

Glycan evolution in response to collaboration, conflict, and constraint

Glycans, oligo- and polysaccharides secreted or attached to proteins and lipids, cover the surfaces of all cells and have a regulatory capacity and structural diversity beyond any other class of biological molecule. Glycans may have evolved these properties because they mediate cellular interactions and often face pressure to evolve new functions rapidly. We approach this idea two ways. First, we discuss evolutionary innovation. Glycan synthesis, regulation, and mode of chemical interaction influence the spectrum of new forms presented to evolution. Second, we describe the evolutionary conflicts that arise when alleles and individuals interact. Glycan regulation and diversity are integral to these biological negotiations. Glycans are tasked with such an amazing diversity of functions that no study of cellular interaction can begin without considering them. We propose that glycans predominate the cell surface because their physical and chemical properties allow the rapid innovation required of molecules on the frontlines of evolutionary conflict.

Effects of atrazine and chlorpyrifos on activity and transcription of glutathione S-transferase in common carp (Cyprinus carpio L.)

Glutathione S-transferase isoenzymes (GSTs) play a critical role in detoxification pathways. Here we report the tissue distribution of four antioxidant GSTs gene in common carp, and their expression profiles. We also investigated the GSTs activity in different tissues after exposure to the agricultural chemicals atrazine (ATR), chlorpyrifos (CPF), and their mixture. Relative changes in the mRNA abundance of the GST isoforms were examined by real time PCR in liver, brain, kidney and gill of common carp. After exposure and recovery, we observed a statistically significant decrease in the GSTs activity in animals exposed to high concentrations of ATR (428 μg/L), CPF (116 μg/L), and their mixture (113 μg/L). At basal levels of tissue expression, four GSTs transcript were detected in liver, brain, kidney, and gill. High expression levels were found in all examined tissues. Transcription of some GST isoforms, GST kappa (GSTK), GST theta (GSTT) and GST rho (GSTR), decreased after exposure to CPF and ATR for the entire experimental period in both the kidney and gill. However, increased transcription of GST mu (GSTM) was observed in the kidney or gill 20 d after exposure to ATR or CPF, respectively. Transcription of both GSTT and GSTR was inhibited for the entire experimental period in the brain, kidney and gill of animals exposed to the ATR/CPF mixture, but transcription of GSTM was induced in the liver after 40 d of exposure. In summary, changes in the GSTs activity and their transcription varied within each organ and among organs of common carp after exposure to ATR, CPF, and their mixture.

Mutational analysis of human glutathione transferase A2-2 identifies structural elements supporting high activity with the prodrug azathioprine

Glutathione transferase (GST) A2-2 is the human enzyme displaying the highest catalytic activity with the prodrug azathioprine (Aza). The reaction releases pharmacologically active 6-mercaptopurine by displacing the imidazole moiety from the Aza molecule. The GST-catalyzed reaction is of medical significance, since high rates of Aza activation may lead to adverse side effects in treated patients. The present study involves structure-activity relationships in GST A2-2 variants. Chimeric GSTs were previously generated by DNA shuffling and two peptide segments, one N-terminal and one C-terminal, were identified as primary determinants of Aza activity. The segments contain several residues of the substrate-binding H-site and their significance for supporting high Aza activity was investigated. Substitution of the corresponding two small regions in the low-activity human GST A3-3 or rat GST A3-3 by the human GST A2-2 segments generated chimeras with ∼10-fold enhanced Aza activity. The H-site residues Met208 and Leu213 in the C-terminal segment of GST A2-2 were mutated to produce a library with all possible residue combinations. At a calculated 93% library coverage, all of the 1880 mutants examined showed wild-type or decreased Aza activity, even though some retained activities with alternative substrates, further emphasizing the importance of this region for the targeted activity.

Heterologous expression and characterization of a sigma glutathione S-transferase involved in carbaryl detoxification from oriental migratory locust, Locusta migratoria manilensis (Meyen)

Glutathione S-transferases (GSTs) play a major role in detoxification of xenobiotics and resistance to insecticides in insects. In the present study, a sigma-class GST gene (LmGSTs3) was identified from the locust, Locusta migratoria manilensis. Its full-length cDNA sequence is 828 bp containing an open reading frame (ORF) of 612 bp that encodes 204 amino acid residues. The predicted protein molecular mass and pI are 23.4 kDa and 7.62, respectively. Recombinant LmGSTs3 was heterologously expressed in Escherichia coli as a soluble fusion protein. Its optimal activity was observed at pH 8.0. Incubation for 30 min at temperatures below 40 °C scarcely affected activity. The LmGSTs3 at pH values between 4.0 and 11.0 retained more than 80% of its original activity. Ethacrynic acid and cibacron blue were very effective inhibitors of LmGSTs3 with I50-values 1.7 and 3.7 μM, respectively. In response to heavy metal (CuSO4, CdCl2) exposure there was a concentration-dependent and time-dependent decrease in activity. The nymph mortalities after carbaryl treatment increased 38.7% after LmGSTs3 were silenced. These results suggest that LmGSTs3 may be involved in carbaryl detoxification in L. migratoria manilensis.

Molecular evolution of candidate sour taste receptor gene PKD1L3 in mammals

The PKD1L3 gene encodes an ion channel protein that can interact with the PKD2L1 protein to form a candidate sour taste receptor. In the present study, we have analyzed the evolutionary patterns of PKD1L3 genes from 10 mammalian species. The results showed that PKD1L3 genes have evolved under a dominant purifying selection force. However, for some branches and sites, PKD1L3 genes were detected to have been operated by positive selection. Moreover, some of these positive evolutionary sites are likely to participate in acid stimulus recognition. In rodents, PKD1L3 genes evolved more rapidly than other mammalian lineages. Combined with other functional research reports, our results suggest that rodents may not be the most appropriate model for functional research on the PKD1L3 gene.

Cys-X scanning for expansion of active-site residues and modulation of catalytic functions in a glutathione transferase

We propose Cys-X scanning as a semisynthetic approach to engineer the functional properties of recombinant proteins. As in the case of Ala scanning, key residues in the primary structure are identified, and one of them is replaced by Cys via site-directed mutagenesis. The thiol of the residue introduced is subsequently modified by alternative chemical reagents to yield diverse Cys-X mutants of the protein. This chemical approach is orthogonal to Ala or Cys scanning and allows the expansion of the repertoire of amino acid side chains far beyond those present in natural proteins. In its present application, we have introduced Cys-X residues in human glutathione transferase (GST) M2-2, replacing Met-212 in the substrate-binding site. To achieve selectivity of the modifications, the Cys residues in the wild-type enzyme were replaced by Ala. A suite of simple substitutions resulted in a set of homologous Met derivatives ranging from normethionine to S-heptyl-cysteine. The chemical modifications were validated by HPLC and mass spectrometry. The derivatized mutant enzymes were assayed with alternative GST substrates representing diverse chemical reactions: aromatic substitution, epoxide opening, transnitrosylation, and addition to an ortho-quinone. The Cys substitutions had different effects on the alternative substrates and differentially enhanced or suppressed catalytic activities depending on both the Cys-X substitution and the substrate assayed. As a consequence, the enzyme specificity profile could be changed among the alternative substrates. The procedure lends itself to large-scale production of Cys-X modified protein variants.

Ancient and recent adaptive evolution of primate non-homologous end joining genes

In human cells, DNA double-strand breaks are repaired primarily by the non-homologous end joining (NHEJ) pathway. Given their critical nature, we expected NHEJ proteins to be evolutionarily conserved, with relatively little sequence change over time. Here, we report that while critical domains of these proteins are conserved as expected, the sequence of NHEJ proteins has also been shaped by recurrent positive selection, leading to rapid sequence evolution in other protein domains. In order to characterize the molecular evolution of the human NHEJ pathway, we generated large simian primate sequence datasets for NHEJ genes. Codon-based models of gene evolution yielded statistical support for the recurrent positive selection of five NHEJ genes during primate evolution: XRCC4, NBS1, Artemis, POLλ, and CtIP. Analysis of human polymorphism data using the composite of multiple signals (CMS) test revealed that XRCC4 has also been subjected to positive selection in modern humans. Crystal structures are available for XRCC4, Nbs1, and Polλ; and residues under positive selection fall exclusively on the surfaces of these proteins. Despite the positive selection of such residues, biochemical experiments with variants of one positively selected site in Nbs1 confirm that functions necessary for DNA repair and checkpoint signaling have been conserved. However, many viruses interact with the proteins of the NHEJ pathway as part of their infectious lifecycle. We propose that an ongoing evolutionary arms race between viruses and NHEJ genes may be driving the surprisingly rapid evolution of these critical genes.

Because all cells experience DNA damage, they must also have mechanisms for repairing DNA. When the proteins that repair DNA malfunction, mutation and disease often result. Based on their fundamental importance, DNA repair proteins would be expected to be well preserved over evolutionary time in order to ensure optimal DNA repair function. However, a previous genome-wide study of molecular evolution in Saccharomyces yeast identified the non-homologous end joining (NHEJ) DNA repair pathway as one of the two most rapidly evolving pathways in the yeast genome. In order to analyze the evolution of this pathway in humans, we have generated large evolutionary sequence sets of NHEJ genes from our primate relatives. Similar to the scenario in yeast, several genes in this pathway are evolving rapidly in primate genomes and in modern human populations. Thus, complex and seemingly opposite selective forces are shaping the evolution of these important DNA repair genes. The finding that NHEJ genes are rapidly evolving in species groups as diverse as yeasts and primates indicates a systematic perturbation of the NHEJ pathway, one that is potentially important to human health.

Molecular evolution and the role of oxidative stress in the expansion and functional diversification of cytosolic glutathione transferases

Background

Cytosolic glutathione transferases (cGST) are a large group of ubiquitous enzymes involved in detoxification and are well known for their undesired side effects during chemotherapy. In this work we have performed thorough phylogenetic analyses to understand the various aspects of the evolution and functional diversification of cGSTs. Furthermore, we assessed plausible correlations between gene duplication and substrate specificity of gene paralogs in humans and selected species, notably in mammalian enzymes and their natural substrates.

Results

We present a molecular phylogeny of cytosolic GSTs that shows that several classes of cGSTs are more ubiquitous and thus have an older ancestry than previously thought. Furthermore, we found that positive selection is implicated in the diversification of cGSTs. The number of duplicate genes per class is generally higher for groups of enzymes that metabolize products of oxidative damage.

Conclusions

1) Protection against oxidative stress seems to be the major driver of positive selection in mammalian cGSTs, explaining the overall expansion pattern of this subfamily;

2) Given the functional redundancy of GSTs that metabolize xenobiotic chemicals, we would expect the loss of gene duplicates, but by contrast we observed a gene expansion of this family, which likely has been favored by: i) the diversification of endogenous substrates; ii) differential tissue expression; and iii) increased specificity for a particular molecule;

3) The increased availability of sequence data from diversified taxa is likely to continue to improve our understanding of the early origin of the different cGST classes.

Residue 234 is a master switch of the alternative-substrate activity profile of human and rodent theta class glutathione transferase T1-1

BACKGROUND

The Theta class glutathione transferase GST T1-1 is a ubiquitously occurring detoxication enzyme. The rat and mouse enzymes have high catalytic activities with numerous electrophilic compounds, but the homologous human GST T1-1 has comparatively low activity with the same substrates. A major structural determinant of substrate recognition is the H-site, which binds the electrophile in proximity to the nucleophilic sulfur of the second substrate glutathione. The H-site is formed by several segments of amino acid residues located in separate regions of the primary structure. The C-terminal helix of the protein serves as a lid over the active site, and contributes several residues to the H-site.

METHODS

Site-directed mutagenesis of the H-site in GST T1-1 was used to create the mouse Arg234Trp for comparison with the human Trp234Arg mutant and the wild-type rat, mouse, and human enzymes. The kinetic properties were investigated with an array of alternative electrophilic substrates to establish substrate selectivity profiles for the different GST T1-1 variants.

RESULTS

The characteristic activity profile of the rat and mouse enzymes is dependent on Arg in position 234, whereas the human enzyme features Trp. Reciprocal mutations of residue 234 between the rodent and human enzymes transform the substrate-selectivity profiles from one to the other.

CONCLUSIONS

H-site residue 234 has a key role in governing the activity and substrate selectivity profile of GST T1-1.

GENERAL SIGNIFICANCE

The functional divergence between human and rodent Theta class GST demonstrates that a single point mutation can enable or suppress enzyme activities with different substrates.

Minor modifications of the C-terminal helix reschedule the favored chemical reactions catalyzed by theta class glutathione transferase T1-1

Adaptive responses to novel toxic challenges provide selective advantages to organisms in evolution. Glutathione transferases (GSTs) play a pivotal role in the cellular defense because they are main contributors to the inactivation of genotoxic compounds of exogenous as well as of endogenous origins. GSTs are promiscuous enzymes catalyzing a variety of chemical reactions with numerous alternative substrates. Despite broad substrate acceptance, individual GSTs display pronounced selectivities such that only a limited number of substrates are transformed with high catalytic efficiency. The present study shows that minor structural changes in the C-terminal helix of mouse GST T1-1 induce major changes in the substrate-activity profile of the enzyme to favor novel chemical reactions and to suppress other reactions catalyzed by the parental enzyme.

Glutathione S-transferase as a biomarker in the Antarctic bivalve Laternula elliptica after exposure to the polychlorinated biphenyl mixture Aroclor 1254

Glutathione S-transferases (GSTs) are a family of multifunctional enzymes involved in cellular detoxification that catalyze the attachment of electrophilic substrates to glutathione. Two classes of GSTs related to the rho and sigma classes of enzymes in Antarctic bivalves have been cloned from Laternula elliptica. The full-length cDNA of rho class GST (leGSTr) is 1530bp in length and contains an open reading frame (ORF) of 672bp encoding 223 amino acid residues. The deduced amino acid sequences of this gene have 41% and 40% identity to rho class GSTs from Ctenopharyngodon idella and Pleuronectes platessa, respectively. The sigma class GST (leGSTs) cDNA, however, is 1127bp in length and contains an ORF of 696bp encoding 231 amino acid residues. The deduced amino acid sequences share only 22% identity with sigma class GST from Xenopus laevis. The transcriptional expression of leGSTr, leGSTs, and leGSTp cloned in our previous study were examined using real-time polymerase chain reaction in response to exposure to a polychlorinated biphenyl (PCB) mixture. The expressions of these three GST transcripts were rapidly upregulated, although they showed different expression levels and patterns within each isoform. Moreover, leGSTs was the most upregulated in the gill and digestive gland in response to PCB exposure. The recombinant GSTs were highly expressed in transformed Escherichia coli, and their kinetic properties were studied with various substrates. As a result, the three classes of GSTs were found to have diverse biological functions and were responsible for different enzymatic features.

Residue 234 in glutathione transferase T1-1 plays a pivotal role in the catalytic activity and the selectivity against alternative substrates

GST (glutathione transferase) T1-1 plays an important role in the biotransformation of halogenated alkanes, which are used in large quantities as solvents and occur as environmental pollutants. Many reactions that are catalysed by GST T1-1 qualify as detoxification processes, but some reactions with dihalogenated alkanes lead to reactive products more toxic than the substrates. Murine GST T1-1 is particularly active with dichloromethane, which may explain the high carcinogenicity of dichloromethane in the mouse. Human GST T1-1 activity is considerably lower with halogenated hydrocarbons and some related substrates. Human GST T1-1 is polymorphic with a frequent null phenotype, suggesting that it is advantageous, under some circumstances, to lack the functional enzyme, which catalyses GSH conjugations that may cause bioactivation. The present study shows that amino acid residue 234 is a determinant of the differences in catalytic efficiency between the human and the rodent enzymes. The replacement of Trp234 in human GST T1-1 by arginine, found in the rodent enzyme, enhanced the alkyltransferase activity by an order of magnitude with a series of homologous iodoalkanes and some typical GST substrates. The specific activity of the alternative mutant Trp234-->Lys was lower than for the parental human GST T1-1 with many substrates, showing that a positive charge is not sufficient for increased activity. The enhanced activity of Trp234-->Arg with alkylating agents was dependent on the substrate tested, whereas no increase of the peroxidase activity with cumene hydroperoxide was noted. Residue 234 therefore is also involved in the control of the substrate selectivity of GST T1-1.

Gene cloning of a sigma class glutathione S-transferase from abalone (Haliotis diversicolor) and expression analysis upon bacterial challenge

Glutathione S-transferases (GSTs) are a multigene family of xenobiotic metabolizing phase II detoxification enzymes which take part in many pathological and physiological processes, and which can potentially be used as indicators and biomarkers for cancer diagnoses and organic or inorganic pollutant exposure. In this study, a full-length cDNA of a sigma class GST (abGSTsigma) (GenBank accession number EF546619) from variously colored abalone (Haliotis diversicolor) was identified. It was 1328bp containing an open reading frame of 624bp, encoding 208 amino acid residues with a predicted protein molecular weight of 23.67kDa and an estimated pI of 5.67. Sequence analysis showed that the predicted protein sequence of abGSTsigma cDNA contained the conserved domain of the GST_N_Sigma_like (PSSM: cd03039) and GST_C_Sigma_like (PSSM: cd03192). Alignment analysis demonstrated that the abGSTsigma of H. diversicolor was in a branch position with other known class sigma GSTs from different organisms. The abGSTsigma mRNA was distributed in multiple tissues tested and was highly demonstrated in the gill and mantle of normal abalones. In bacteria-challenged abalone, the abGSTsigma gene was significantly expressed in the hemocytes, gill, mantle and digestive gland and the total GSTs enzyme and SOD were also induced in the four tissues. The increased activities of SOD and GSTs can result in the elimination of reactive oxygen species (ROS) indicating antioxidant activities involved. The preliminary work revealed that the sigma class glutathione S-transferase gene abGSTsigma, a phase II detoxification enzyme, had a positive response to bacterial challenge, and that will lead to an insightful study on elucidating the interactions between immune responses and biotransformation exerted by abGSTsigma.

Divergent evolution and molecular adaptation in the Drosophila odorant-binding protein family: inferences from sequence variation at the OS-E and OS-F genes

BACKGROUND

The Drosophila Odorant-Binding Protein (Obp) genes constitute a multigene family with moderate gene number variation across species. The OS-E and OS-F genes are the two phylogenetically closest members of this family in the D. melanogaster genome. In this species, these genes are arranged in the same genomic cluster and likely arose by tandem gene duplication, the major mechanism proposed for the origin of new members in this olfactory-system family.

RESULTS

We have analyzed the genomic cluster encompassing OS-E and OS-F genes (Obp83 genomic region) to determine the role of the functional divergence and molecular adaptation on the Obp family size evolution. We compared nucleotide and amino acid variation across 18 Drosophila and 4 mosquito species applying a phylogenetic-based maximum likelihood approach complemented with information of the OBP three-dimensional structure and function. We show that, in spite the OS-E and OS-F genes are currently subject to similar and strong selective constraints, they likely underwent divergent evolution. Positive selection was likely involved in the functional diversification of new copies in the early stages after the gene duplication event; moreover, it might have shaped nucleotide variation of the OS-E gene concomitantly with the loss of functionally related members. Besides, molecular adaptation likely affecting the functional OBP conformational changes was supported by the analysis of the evolution of physicochemical properties of the OS-E protein and the location of the putative positive selected amino acids on the OBP three-dimensional structure.

CONCLUSION

Our results support that positive selection was likely involved in the functional differentiation of new copies of the OBP multigene family in the early stages after their birth by gene duplication; likewise, it might shape variation of some members of the family concomitantly with the loss of functionally related genes. Thus, the stochastic gene gain/loss process coupled with the impact of natural selection would influence the observed OBP family size.

Molecular evolution of bovine Toll-like receptor 2 suggests substitutions of functional relevance

Background

There is accumulating evidence that polymorphism in Toll-like receptor (TLR) genes might be associated with disease resistance or susceptibility traits in livestock. Polymorphic sites affecting TLR function should exhibit signatures of positive selection, identified as a high ratio of non-synonymous to synonymous nucleotide substitutions (ω). Phylogeny based models of codon substitution based on estimates of ω for each amino acid position can therefore offer a valuable tool to predict sites of functional relevance. We have used this approach to identify such polymorphic sites within the bovine TLR2 genes from ten Bos indicus and Bos taurus cattle breeds. By analysing TLR2 gene phylogeny in a set of mammalian species and a subset of ruminant species we have estimated the selective pressure on individual sites and domains and identified polymorphisms at sites of putative functional importance.

Results

The ω were highest in the mammalian TLR2 domains thought to be responsible for ligand binding and lowest in regions responsible for heterodimerisation with other TLR-related molecules. Several positively-selected sites were detected in or around ligand-binding domains. However a comparison of the ruminant subset of TLR2 sequences with the whole mammalian set of sequences revealed that there has been less selective pressure among ruminants than in mammals as a whole. This suggests that there have been functional changes during ruminant evolution. Twenty newly-discovered non-synonymous polymorphic sites were identified in cattle. Three of them were localised at positions shaped by positive selection in the ruminant dataset (Leu227Phe, His305Pro, His326Gln) and in domains involved in the recognition of ligands. His326Gln is of particular interest as it consists of an exchange of differentially-charged amino acids at a position which has previously been shown to be crucial for ligand binding in human TLR2.

Conclusion

Within bovine TLR2, polymorphisms at amino acid positions 227, 305 and 326 map to functionally important sites of TLR2 and should be considered as candidate SNPs for immune related traits in cattle. A final proof of their functional relevance requires further studies to determine their functional effect on the immune response after stimulation with relevant ligands and/or their association with immune related traits in animals.

Emergence of a novel highly specific and catalytically efficient enzyme from a naturally promiscuous glutathione transferase

Redesign of glutathione transferases (GSTs) has led to enzymes with remarkably enhanced catalytic properties. Exchange of substrate-binding residues in GST A1-1 created a GST A4-4 mimic, called GIMFhelix, with >300-fold improved activity with nonenal and suppressed activity with other substrates. In the present investigation GIMFhelix was compared with the naturally-evolved GSTs A1-1 and A4-4 by determining catalytic efficiencies with nine alternative substrates. The enzymes can be represented by vectors in multidimensional substrate-activity space, and the vectors of GIMFhelix and GST A1-1, expressed in kcat/Km values for the alternative substrates, are essentially orthogonal. By contrast, the vectors of GIMFhelix and GST A4-4 have approximately similar lengths and directions. The broad substrate acceptance of GST A1-1 contrasts with the high selectivity of GST A4-4 and GIMFhelix for alkenal substrates. Multivariate analysis demonstrated that among the diverse substrates used, nonenal, cumene hydroperoxide, and androstenedione are major determinants in the portrayal of the three enzyme variants. These GST substrates represent diverse chemistries of naturally occurring substrates undergoing Michael addition, hydroperoxide reduction, and steroid double-bond isomerization, respectively. In terms of function, GIMFhelix is a novel enzyme compared to its progenitor GST A1-1 in spite of 94% amino-acid sequence identity between the enzymes. The redesign of GST A1-1 into GIMFhelix therefore serves as an illustration of divergent evolution leading to novel enzymes by minor structural modifications in the active site. Notwithstanding low sequence identity (60%), GIMFhelix is functionally an isoenzyme of GST A4-4.

Molecular cloning and characterization of three sigma glutathione S-transferases from disk abalone (Haliotis discus discus)

Three novel glutathione S-transferase (GSTs) cDNAs were cloned from a disk abalone (Haliotis dicus discus) cDNA library. Multiple alignment and phylogenetic analysis of three GSTs revealed that their closest relationship is with insect sigma GSTs. Recombinant GSTs were over-expressed in Escherichia coli as soluble fusion proteins. HdGSTS1 and HdGSTS2 were active towards 1-chloro-2,4-dinitrobenzene and ethacrynic acid, whereas HdGSTS3 appeared to be a non-enzymatic GST. Two active GSTs had similar optimum conditions for enzymatic reaction at pH 8.0 and temperature of approximately 30 degrees C. Molecular modeling analysis of three GSTs implicates their diverse active sites as being responsible for their different enzymatic features. Three sigma GSTs had significantly different expression patterns and levels of expression in abalone tissues, indicating their different functions. After 48 h-exposure to three model marine pollutants, only HdGSTS1 exhibited a proper inducibility, exhibiting its good biomarker potential for organic contaminants in marine environment. In contrast, the other two sigma GSTs revealed a minor role in the response of pollutants exposure.

Proteomic identification, cDNA cloning and enzymatic activity of glutathione S-transferases from the generalist marine gastropod, Cyphoma gibbosum

Glutathione S-transferases (GST) were characterized from the digestive gland of Cyphoma gibbosum (Mollusca; Gastropoda), to investigate the possible role of these detoxification enzymes in conferring resistance to allelochemicals present in its gorgonian coral diet. We identified the collection of expressed cytosolic Cyphoma GST classes using a proteomic approach involving affinity chromatography, HPLC and nano-spray liquid chromatography-tandem mass spectrometry (LC-MS/MS). Two major GST subunits were identified as putative mu-class GSTs; while one minor GST subunit was identified as a putative theta-class GST, apparently the first theta-class GST identified from a mollusc. Two Cyphoma GST cDNAs (CgGSTM1 and CgGSTM2) were isolated by RT-PCR using primers derived from peptide sequences. Phylogenetic analyses established both cDNAs as mu-class GSTs and revealed a mollusc-specific subclass of the GST-mu clade. These results provide new insights into metazoan GST diversity and the biochemical mechanisms used by marine organisms to cope with their chemically defended prey.

Role of invariant tyrosines in a crustacean mu-class glutathione S-transferase from shrimp Litopenaeus vannamei: site-directed mutagenesis of Y7 and Y116

Y6 and Y115 are key amino acids involved in enzyme-substrate interactions in mu-class glutathione S-transferase (GST). They provide electrophilic assistance and stabilize substrates through their hydroxyl groups. Two site-directed mutants (Y7F and Y116F) and the wild-type shrimp GSTs were expressed in Escherichia coli, and the steady-state kinetic parameters were determined using CDNB as the second substrate. The mutants were modeled based on a crystal structure of a mu-class GST to obtain further insights about the changes at the active site. The Y116F mutant had an increase in kcat contrary to Y7F compared to the wild type. Molecular modeling showed that the shrimp GST has a H108 residue that may contribute to compensate and lead to a less deleterious change when conserved tyrosine residues are mutated. This work indicates that shrimp GST is a useful model to understand the catalysis mechanisms in this critical enzyme.

Molecular evolution of glutathione S-transferases in the genus Drosophila

As classical phase II detoxification enzymes, glutathione S-transferases (GSTs) have been implicated in insecticide resistance and may have evolved in response to toxins in the niche-defining feeding substrates of Drosophila species. We have annotated the GST genes of the 12 Drosophila species with recently sequenced genomes and analyzed their molecular evolution. Gene copy number variation is attributable mainly to unequal crossing-over events in the large delta and epsilon clusters. Within these gene clusters there are also GST genes with slowly diverging orthologs. This implies that they have their own unique functions or have spatial/temporal expression patterns that impose significant selective constraints. Searches for positively selected sites within the GSTs identified G171K in GSTD1, a protein that has previously been shown to be capable of metabolizing the insecticide DDT. We find that the same radical substitution (G171K) in the substrate-binding domain has occurred at least three times in the Drosophila radiation. Homology-modeling places site 171 distant from the active site but adjacent to an alternative DDT-binding site. We propose that the parallel evolution observed at this site is an adaptive response to an environmental toxin and that sequencing of historical alleles suggests that this toxin was not a synthetic insecticide.

Conceptual bases for quantifying the role of the environment on gene evolution: the participation of positive selection and neutral evolution

To demonstrate that a given change in the environment has contributed to the emergence of a given genotypic and phenotypic shift during the course of evolution, one should ask to what extent such shifts would have occurred without environmental change. Of course, such tests are rarely practical but phenotypic novelties can still be correlated to genomic shifts in response to environmental changes if enough information is available. We surveyed and re-evaluated the published data in order to estimate the role of environmental changes on the course of species and genomic evolution. Only a few published examples clearly demonstrate a causal link between a given environmental change and the fixation of a genomic variant resulting in functional modification (gain, loss or alteration of function). Many others suggested a link between a given phenotypic shift and a given environmental change but failed to identify the underlying genomic determinant(s) and/or the associated functional consequence(s). The proportion of genotypic and phenotypic variation that is fixed concomitantly with environmental changes is often considered adaptive and hence, the result of positive selection, even though alternative causes, such as genetic drift, are rarely investigated. Therefore, the second aim herein is to review evidence for the mechanisms leading to fixation.

Approaches for identifying targets of positive selection

Despite significant advancements in both empirical and theoretical population genetics throughout the past century, fundamental questions about the evolutionary forces that shape genomic diversity remain unresolved. Perhaps foremost among these are the strength and frequency of adaptive evolution. To quantify these parameters, statistical tools are needed that are capable of effectively identifying targets of positive selection throughout the genome in an unbiased manner, and functional approaches are needed that are capable of connecting these identified genotypes with the resulting adaptively significant phenotypes. Here we review recent advancements in both statistical and empirical methodology, and discuss important challenges and opportunities that remain as researchers continue to uncouple the relative importance of stochastic and deterministic factors in the evolution of natural populations.

Engineering the enantioselectivity of glutathione transferase by combined active-site mutations and chemical modifications

Based on the crystal structure of human glutathione transferase M1-1, cysteine residues were introduced in the substrate-binding site of a Cys-free mutant of the enzyme, which were subsequently alkylated with 1-iodoalkanes. By different combinations of site-specific mutations and chemical modifications of the enzyme the enantioselectivity in the conjugation of glutathione with the epoxide-containing substrates 1-phenylpropylene oxide and styrene-7,8-oxide were enhanced up to 9- and 10-fold. The results also demonstrate that the enantioselectivity can be diminished, or even reversed, by suitable modifications, which can be valuable under some conditions. The redesign of the active-site structure for enhanced or diminished enantioselectivities have divergent requirements for different epoxides, calling for a combinatorial approach involving alternative mutations and chemical modifications to optimize the enantioselectivity for a targeted substrate. This approach outlines a general method of great potential for fine-tuning substrate specificity and tailoring stereoselectivity of recombinant enzymes.

Importance of a hypervariable active-site residue in human Mu class glutathione transferases catalyzing the bioactivation of chemotherapeutic thiopurine prodrugs

Glutathione transferases (GSTs) catalyze the bioactivation of the thiopurine prodrugs azathioprine, cis-6-(2-acetylvinylthio)purine (cAVTP) and trans-6-(2-acetylvinylthio)guanine (tAVTG), thereby releasing the antimetabolites 6-mercaptopurine and 6-thioguanine. In the GST Mu class, GST M1-1 has the highest catalytic efficiency, whereas GST M2-2 and other enzymes are less active. In the evolution of Mu class GSTs, residue 210 appears hypervariable and has particular functional significance. We demonstrate that the catalytic activity of GST M1-1 with cAVTP or tAVTG is successively diminished when wild-type Ser-210 is mutated into Ala followed by Thr. Conversely, mutating wild-type Thr-210 in GST M2-2 into Ala and Ser enhanced the corresponding activities. Comparisons were also made with GST M2-2 distinguished by Gly or Pro in position 210, as well as wild-type GSTs M4-4 and M5-5. The results suggest that the hydroxyl group of Ser in position 210 stabilizes the transition state of the GST-catalyzed reaction. The low activity of GSTs containing Thr in position 210 is probably due to steric hindrance caused by the beta-methyl group of the side chain. The ratios of the different catalytic efficiencies were translated into differences in the Gibbs free energies of transition state stabilization. The effects of the mutations were qualitatively parallel for the alternative substrates, but vary significantly in magnitude. From the evolutionary perspective the data show that a point mutation can alternatively enhance or attenuate the activity with a particular substrate and illustrate the functional plasticity of GSTs.

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.

Multivariate-activity mining for molecular quasi-species in a glutathione transferase mutant library

A library of recombinant glutathione transferases (GSTs) generated by shuffling of DNA encoding human GST M1-1 and GST M2-2 was screened with eight alternative substrates, and the activities were subjected to multivariate analysis. Assays were made in lysates of bacteria in which the GST variants had been expressed. The primary data showed clustering of the activities in eight-dimensional substrate-activity space. For an incisive analysis, the rows of the data matrix, corresponding to the different enzyme variants, were individually scaled to unit length, thus accounting for different expression levels of the enzymes. The columns representing the activities with alternative substrates were subsequently individually normalized to unit variance and a zero mean. By this standardization, the data were adjusted to comparable orders of magnitude. Three molecular quasi-species were recognized by multivariate K-means and principal component analyses. Two of them encompassed the parental GST M1-1 and GST M2-2. A third one diverged functionally by displaying enhanced activities with some substrates and suppressed activities with signature substrates for GST M1-1 and GST M2-2. A fourth cluster contained mutants with impaired functions and was not regarded as a quasi-species. Sequence analysis of representatives of the mutant clusters demonstrated that the majority of the variants in the diverging novel quasi-species were structurally similar to the M1-like GSTs, but distinguished themselves from GST M1-1 by a Ser to Thr substitution in the active site. The data show that multivariate analysis of functional profiles can identify small structural changes influencing the evolution of enzymes with novel substrate-activity profiles.