First-sphere and second-sphere electrostatic effects in the active site of a class mu gluthathione transferase

Switching a HO···π Interaction to a Nonconventional OH···π Hydrogen Bond: A Completed Crystallographic Puzzle

In this article, we present crystallographic and spectroscopic evidence of a tunable system wherein a HO···π interaction switches incrementally to a nonconventional OH···π hydrogen bonding (HB) interaction. More specifically, we report the synthesis of substituted forms of model system 1 to study the effects of aryl ring electronic density on the qualitative characteristics of OH···π hydrogen bonds therein. The OH stretch in experimental infrared data, in agreement with density-functional theory (DFT) calculations, shows continuous red-shifts as the adjacent ring becomes more electron rich. For example, the OH stretch of an amino-substituted analogue is red-shifted by roughly 50 cm^-1 compared to the same stretch in the CF₃ analogue, indicating a significantly stronger HB interaction in the former. Moreover, DFT calculations (ωB97XD/6-311+G**) predict that increasing electronic density on the adjacent top ring reduces the aryl π-OH σ* energy gap with a concomitant enhancement of the OH n-π* energy gap. Consequently, a dominant π-σ* interaction in the amino substituted analogue locks the system in the in-form while a favorable n-π* interaction reverses the orientation of the oxygen-bound hydrogen in its protonated form. Additionally, the ¹H NMR data of various analogues reveal that stronger OH···π interactions in systems with electron-rich aromatic rings slow exchange of the alcoholic proton, thereby revealing coupling with the geminal proton. Finally, X-ray crystallographic analyses of a spectrum of analogues clearly visualize the three distinct stages of "switch"-starting with exclusive HO···π, to partitioned HO···π/OH···π, and finally to achieving exclusive OH···π forms. Furthermore, the crystal structure of the amino analogue reveals an interesting feature in which an extended HB network, involving two conventional (NH···O) and two nonconventional (OH···π) HBs, dimerizes and anchors the molecule in the unit cell.

The mercapturic acid pathway

The mercapturic acid pathway is a major route for the biotransformation of xenobiotic and endobiotic electrophilic compounds and their metabolites. Mercapturic acids (N-acetyl-l-cysteine S-conjugates) are formed by the sequential action of the glutathione transferases, γ-glutamyltransferases, dipeptidases, and cysteine S-conjugate N-acetyltransferase to yield glutathione S-conjugates, l-cysteinylglycine S-conjugates, l-cysteine S-conjugates, and mercapturic acids; these metabolites constitute a "mercapturomic" profile. Aminoacylases catalyze the hydrolysis of mercapturic acids to form cysteine S-conjugates. Several renal transport systems facilitate the urinary elimination of mercapturic acids; urinary mercapturic acids may serve as biomarkers for exposure to chemicals. Although mercapturic acid formation and elimination is a detoxication reaction, l-cysteine S-conjugates may undergo bioactivation by cysteine S-conjugate β-lyase. Moreover, some l-cysteine S-conjugates, particularly l-cysteinyl-leukotrienes, exert significant pathophysiological effects. Finally, some enzymes of the mercapturic acid pathway are described as the so-called "moonlighting proteins," catalytic proteins that exert multiple biochemical or biophysical functions apart from catalysis.

Learning protein constitutive motifs from sequence data

Statistical analysis of evolutionary-related protein sequences provides information about their structure, function, and history. We show that Restricted Boltzmann Machines (RBM), designed to learn complex high-dimensional data and their statistical features, can efficiently model protein families from sequence information. We here apply RBM to 20 protein families, and present detailed results for two short protein domains (Kunitz and WW), one long chaperone protein (Hsp70), and synthetic lattice proteins for benchmarking. The features inferred by the RBM are biologically interpretable: they are related to structure (residue-residue tertiary contacts, extended secondary motifs (α-helixes and β-sheets) and intrinsically disordered regions), to function (activity and ligand specificity), or to phylogenetic identity. In addition, we use RBM to design new protein sequences with putative properties by composing and 'turning up' or 'turning down' the different modes at will. Our work therefore shows that RBM are versatile and practical tools that can be used to unveil and exploit the genotype-phenotype relationship for protein families.

Crystal structure of a class-mu glutathione S-transferase from whiteleg shrimp Litopenaeus vannamei: structural changes in the xenobiotic binding H-site may alter the spectra of molecules bound

Glutathione S-transferases (GSTs) are dimeric proteins that play a key role in phase II cellular detoxification. Here, the first crystal structure of a GST class-mu from marine crustacean shrimp Litopenaeus vannamei is reported at a resolution of 2.0 Å. The coordinates reported here have the lowest sequence identity with previously reported GSTs class-mu deposited at the Protein Data Bank (PDB), although they have subtle conformational differences. One key feature of GST class-mu from L. vannamei is the active site crevice markedly reduced when it is compared with other GSTs class-mu. This finding together with the chemical change of residues into the cavity (F112 and Y210) points to a particular specialization in which smallest xenobiotics with nonstandard chemical characteristics can be bound to the H-site. This suggests that marine organisms have evolved structural strategies to provide efficient selectivity toward xenobiotics to be disposed of by the phase II detoxification process.

Molecular and biochemical characterization of dehydroascorbate reductase from a stress adapted C4 plant, pearl millet [Pennisetum glaucum (L.) R. Br]

KEY MESSAGE

PgDHAR was isolated from Pennisetum glaucum. PgDHAR responded to abiotic stress and exhibited enzyme activity at broad ranges of temperature, pH and substrate concentrations suggesting its role in stress tolerance.

ABSTRACT

Dehydroascorbate reductase (EC 1.8.5.1) is a crucial enzyme actively involved in the recycling of ascorbate redox pool in the cellular environment. In this study, the full-length cDNA coding for DHAR polypeptide and its corresponding gene was isolated from Pennisetum glaucum (PgDHAR). PgDHAR encodes a polypeptide of 213 amino acids with a predicted molecular mass of 23.4 kDa and shares 80-75 % sequence homology with DHAR from other plants. The heterologously expressed recombinant PgDHAR protein exhibited activity in a wide range of substrate concentrations. The recombinant PgDHAR is thermostable and retains its activity over a broad pH range. Furthermore, transcript level of PgDHAR is quantitatively up-regulated in response to temperature. On the whole, PgDHAR alone or in combination with other genes of ascorbate-glutathione cycle can be used for the development of stress tolerant as well as nutritionally improved food crop with enhanced ascorbic acid content.

Chemical proteomics with sulfonyl fluoride probes reveals selective labeling of functional tyrosines in glutathione transferases

Chemical probes have great potential for identifying functional residues in proteins in crude proteomes. Here we studied labeling sites of chemical probes based on sulfonyl fluorides (SFs) on plant and animal proteomes. Besides serine proteases and many other proteins, SF-based probes label Tyr residues in glutathione transferases (GSTs). The labeled GSTs represent four different GST classes that share less than 30% sequence identity. The targeted Tyr residues are located at similar positions in the promiscuous substrate binding site and are essential for GST function. The high selectivity of SF-based probes for functional Tyr residues in GSTs illustrates how these probes can be used for functional studies of GSTs and other proteins in crude proteomes.

Low-molecular-weight thiols in thiol-disulfide exchange

SIGNIFICANCE

Oxidative stress is widely invoked in inflammation, aging, and complex diseases. To avoid unwanted oxidations, the redox environment of cellular compartments needs to be tightly controlled. The complementary action of oxidoreductases and of high concentrations of low-molecular-weight (LMW) nonprotein thiols plays an essential role in maintaining the redox potential of the cell in balance.

RECENT ADVANCES

While LMW thiols are central players in an extensive range of redox regulation/metabolism processes, not all organisms use the same thiol cofactors to this effect, as evidenced by the recent discovery of mycothiol (MSH) and bacillithiol (BSH) among different gram-positive bacteria.

CRITICAL ISSUES

LMW thiol-disulfide exchange processes and their cellular implications are often oversimplified, as only the biology of the free thiols and their symmetrical disulfides is considered. In bacteria under oxidative stress, especially where concentrations of different LMW thiols are comparable [e.g., BSH, coenzyme A (CoA), and cysteine (Cys) in many low-G+C gram-positive bacteria (Firmicutes)], mixed disulfides (e.g., CoASSB and CySSCoA) must surely be major thiol-redox metabolites that need to be taken into consideration.

FUTURE DIRECTIONS

There are many microorganisms whose LMW thiol-redox buffers have not yet been identified (either bioinformatically or experimentally). Many elements of BSH and MSH redox biochemistry remain to be explored. The fundamental biophysical properties, thiol pK(a) and redox potential, have not yet been determined, and the protein interactome in which the biothiols MSH and BSH are involved needs further exploration.

Understanding the pK(a) of redox cysteines: the key role of hydrogen bonding

Many cellular functions involve cysteine chemistry via thiol-disulfide exchange pathways. The nucleophilic cysteines of the enzymes involved are activated as thiolate. A thiolate is much more reactive than a neutral thiol. Therefore, determining and understanding the pK(a)s of functional cysteines are important aspects of biochemistry and molecular biology with direct implications for redox signaling. Here, we describe the experimental and theoretical methods to determine cysteine pK(a) values, and we examine the factors that control these pK(a)s. Drawing largely on experience gained with the thioredoxin superfamily, we examine the roles of solvation, charge-charge, helix macrodipole, and hydrogen bonding interactions as pK(a)-modulating factors. The contributions of these factors in influencing cysteine pK(a)s and the associated chemistry, including the relevance for the reaction kinetics and thermodynamics, are discussed. This analysis highlights the critical role of direct hydrogen bonding to the cysteine sulfur as a key factor modulating the equilibrium between thiol S-H and thiolate S(-). This role is easily understood intuitively and provides a framework for biochemical functional insights.

Structures of yeast glutathione-S-transferase Gtt2 reveal a new catalytic type of GST family

Glutathione-S-transferases (GSTs) are ubiquitous detoxification enzymes that catalyse the conjugation of electrophilic substrates to glutathione. Here, we present the crystal structures of Gtt2, a GST of Saccharomyces cerevisiae, in apo and two ligand-bound forms, at 2.23 A, 2.20 A and 2.10 A, respectively. Although Gtt2 has the overall structure of a GST, the absence of the classic catalytic essential residues--tyrosine, serine and cysteine--distinguishes it from all other cytosolic GSTs of known structure. Site-directed mutagenesis in combination with activity assays showed that instead of the classic catalytic residues, a water molecule stabilized by Ser129 and His123 acts as the deprotonator of the glutathione sulphur atom. Furthermore, only glycine and alanine are allowed at the amino-terminus of helix-alpha1 because of stereo-hindrance. Taken together, these results show that yeast Gtt2 is a novel atypical type of cytosolic GST.

Closely related antibody receptors exploit fundamentally different strategies for steroid recognition

Molecular recognition by the adaptive immune system relies on specific high-affinity antibody receptors that are generated from a restricted set of starting sequences through homologous recombination and somatic mutation. The steroid binding antibody DB3 and the catalytic Diels-Alderase antibody 1E9 derive from the same germ line sequences but exhibit very distinct specificities and functions. However, mutation of only two of the 36 sequence differences in the variable domains, Leu(H47)Trp and Arg(H100)Trp, converts 1E9 into a high-affinity steroid receptor with a ligand recognition profile similar to DB3. To understand how these changes switch binding specificity and function, we determined the crystal structures of the 1E9 Leu(H47)Trp/Arg(H100)Trp double mutant (1E9dm) as an unliganded Fab at 2.05 A resolution and in complex with two configurationally distinct steroids at 2.40 and 2.85 A. Surprisingly, despite the functional mimicry of DB3, 1E9dm employs a distinct steroid binding mechanism. Extensive structural rearrangements occur in the combining site, where residue H47 acts as a specificity switch and H100 adapts to different ligands. Unlike DB3, 1E9dm does not use alternative binding pockets or different sets of hydrogen-bonding interactions to bind configurationally distinct steroids. Rather, the different steroids are inserted more deeply into the 1E9dm combining site, creating more hydrophobic contacts that energetically compensate for the lack of hydrogen bonds. These findings demonstrate how subtle mutations within an existing molecular scaffold can dramatically modulate the function of immune receptors by inducing unanticipated, but compensating, mechanisms of ligand interaction.

Characterization of the binding of 8-anilinonaphthalene sulfonate to rat class Mu GST M1-1

Molecular docking and ANS-displacement experiments indicated that 8-anilinonaphthalene sulfonate (ANS) binds the hydrophobic site (H-site) in the active site of dimeric class Mu rGST M1-1. The naphthalene moiety provides most of the van der Waals contacts at the ANS-binding interface while the anilino group is able to sample different rotamers. The energetics of ANS binding were studied by isothermal titration calorimetry (ITC) over the temperature range of 5-30 degrees C. Binding is both enthalpically and entropically driven and displays a stoichiometry of one ANS molecule per subunit (or H-site). ANS binding is linked to the uptake of 0.5 protons at pH 6.5. Enthalpy of binding depends linearly upon temperature yielding a DeltaC(p) of -80+/-4 cal K(-1) mol(-1) indicating the burial of solvent-exposed nonpolar surface area upon ANS-protein complex formation. While ion-pair interactions between the sulfonate moiety of ANS and protein cationic groups may be significant for other ANS-binding proteins, the binding of ANS to rGST M1-1 is primarily hydrophobic in origin. The binding properties are compared with those of other GSTs and ANS-binding proteins.

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.

Marine glutathione S-transferases

The aquatic environment is generally affected by the presence of environmental xenobiotic compounds. One of the major xenobiotic detoxifying enzymes is glutathione S-transferase (GST), which belongs to a family of multifunctional enzymes involved in catalyzing nucleophilic attack of the sulfur atom of glutathione (gamma-glutamyl-cysteinylglycine) to an electrophilic group on metabolic products or xenobiotic compounds. Because of the unique nature of the aquatic environment and the possible pollution therein, the biochemical evolution in terms of the nature of GSTs could by uniquely expressed. The full complement of GSTs has not been studied in marine organisms, as very few aquatic GSTs have been fully characterized. The focus of this article is to present an overview of the GST superfamily and their critical role in the survival of organisms in the marine environment, emphasizing the critical roles of GSTs in the detoxification of marine organisms and the unique characteristics of their GSTs compared to those from non-marine organisms.

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.

Probing the Mechanism of GSH Activation in Schistosoma haematobium Glutathione-S-transferase by Site-directed Mutagenesis and X-ray Crystallography

During turnover, the catalytic tyrosine residue (Tyr10) of the sigma class Schistosoma haematobium wild-type glutathione-S-transferase is expected to switch alternately in and out of the reduced glutathione-binding site (G-site). The Tyrout10 conformer forms a pi-cation interaction with the guanidinium group of Arg21. As in other similar glutathione-S-transferases, the catalytic Tyr has a low pKa of 7.2. In order to investigate the catalytic role of Tyr10, and the structural and functional roles of Arg21, we carried out structural studies on two Arg21 mutants (R21L and R21Q) and a Tyr10 mutant, Y10F. Our crystallographic data for the two Arg21 mutants indicate that only the Tyrout10 conformation is populated, thereby excluding a role of Arg21 in the stabilisation of the out conformation. However, Arg21 was confirmed to be catalytically important and essential for the low pKa of Tyr10. Upon comparison with structural data generated for reduced glutathione-bound and inhibitor-bound wild-type enzymes, it was observed that the orientations of Tyr10 and Arg35 are concerted and that, upon ligand binding, minor rearrangements occur within conserved residues in the active site loop. These rearrangements are coupled to quaternary rigid-body movements at the dimer interface and alterations in the localisation and structural order of the C-terminal domain.

Contribution of the mu loop to the structure and function of rat glutathione transferase M1-1

The "mu loop," an 11-residue loop spanning amino acid residues 33-43, is a characteristic structural feature of the mu class of glutathione transferases. To assess the contribution of the mu loop to the structure and function of rat GST M1-1, amino acid residues 35-44 (35GDAPDYDRSQ44) were excised by deletion mutagenesis, resulting in the "Deletion Enzyme." Kinetic studies reveal that the Km values of the Deletion Enzyme are markedly increased compared with those of the wild-type enzyme: 32-fold for 1-chloro-2,4-dinitrobenzene, 99-fold for glutathione, and 880-fold for monobromobimane, while the Vmax value for each substrate is increased only modestly. Results from experiments probing the structure of the Deletion Enzyme, in comparison with that of the wild-type enzyme, suggest that the secondary and quaternary structures have not been appreciably perturbed. Thermostability studies indicate that the Deletion Enzyme is as stable as the wild-type enzyme at 4 degrees C and 10 degrees C, but it rapidly loses activity at 25 degrees C, unlike the wild-type enzyme. In the temperature range of 4 degrees C through 25 degrees C, the loss of activity of the Deletion Enzyme is not the result of a change in its structure, as determined by circular dichroism spectroscopy and sedimentation equilibrium centrifugation. Collectively, these results indicate that the mu loop is not essential for GST M1-1 to maintain its structure nor is it required for the enzyme to retain some catalytic activity. However, it is an important determinant of the enzyme's affinity for its substrates.

Are AM1 ligand-protein binding enthalpies good enough for use in the rational design of new drugs?

We have examined the performance of semiempirical quantum mechanical methods in solving the problem of accurately predicting protein-ligand binding energies and geometries. Firstly, AM1 and PM3 geometries and binding enthalpies between small molecules that simulate typical ligand-protein interactions were compared with high level quantum mechanical techniques that include electronic correlation (e.g., MP2 or B3LYP). Species studied include alkanes, aromatic systems, molecules including groups with hypervalent sulfur or with donor or acceptor hydrogen bonding capability, as well as ammonium or carboxylate ions. B3LYP/6-311+G(2d,p) binding energies correlated very well with the BSSE corrected MP2/6-31G(d) values. AM1 binding enthalpies also showed good correlation with MP2 values, and their systematic deviation is acceptable when enthalpies are used for the comparison of interaction energies between ligands and a target. PM3 otherwise gave erratic energy differences in comparison to the B3LYP or MP2 approaches. As one would expect, the geometries of the binding complexes showed the known limitations of the semiempirical and DFT methods. AM1 calculations were subsequently applied to a test set consisting of "real" protein active site-ligand complexes. Preliminary results indicate that AM1 could be a valuable tool for the design of new drugs using proteins as templates. This approach also has a reasonable computational cost. The ligand-protein X-ray structures were reasonably reproduced by AM1 calculations and the corresponding AM1 binding enthalpies are in agreement with the results from the "small molecules" test set.

Insights into the Catalytic Mechanism of Glutathione S-Transferase: The Lesson from Schistosoma haematobium

Glutathione S-transferases (GSTs) are involved in detoxification of xenobiotic compounds and in the biosynthesis of important metabolites. All GSTs activate glutathione (GSH) to GS(-); in many GSTs, this is accomplished by a Tyr at H-bonding distance from the sulfur of GSH. The high-resolution structure of GST from Schistosoma haematobium revealed that the catalytic Tyr occupies two alternative positions, one external, involving a pi-cation interaction with the conserved Arg21, and the other inside the GSH binding site. The interaction with Arg21 lowers the pK(a) of the catalytic Tyr10, as required for catalysis. Examination of several other GST structures revealed the presence of an external pocket that may accommodate the catalytic Tyr, and suggested that the change in conformation and acidic properties of the catalytic Tyr may be shared by other GSTs. Arginine and two other residues of the external pocket constitute a conserved structural motif, clearly identified by sequence comparison.

Tertiary Interactions Stabilise the C-terminal Region of Human Glutathione Transferase A1-1: a Crystallographic and Calorimetric Study

The C-terminal region in class Alpha glutathione transferase A1-1 (GSTA1-1), which forms an amphipathic alpha-helix (helix 9), is known to contribute to the catalytic and non-substrate ligand-binding functions of the enzyme. The region in the apo protein is proposed to be disordered which, upon ligand binding at the active-site, becomes structured and localised. Because Ile219 plays a pivotal role in the stability and localisation of the region, the role of tertiary interactions mediated by Ile219 in determining the conformation and dynamics of the C-terminal region were studied. Ligand-binding microcalorimetric and X-ray structural data were obtained to characterise ligand binding at the active-site and the associated localisation of the C-terminal region. In the crystal structure of the I219A hGSTA1-1.S-hexylglutathione complex, the C-terminal region of one chain is mobile and not observed (unresolved electron density), whereas the corresponding region of the other chain is localised and structured as a result of crystal packing interactions. In solution, the mutant C-terminal region of both chains in the complex is mobile and delocalised resulting in a hydrated, less hydrophobic active-site and a reduction in the affinity of the protein for S-hexylglutathione. Complete dehydration of the active-site, important for maintaining the highly reactive thiolate form of glutathione, requires the binding of ligands and the subsequent localisation of the C-terminal region. Thermodynamic data demonstrate that the mobile C-terminal region in apo hGSTA1-1 is structured and does not undergo ligand-induced folding. Its close proximity to the surface of the wild-type protein is indicated by the concurrence between the observed heat capacity change of complex formation and the type and amount of surface area that becomes buried at the ligand-protein interface when the C-terminal region in the apo protein assumes the same localised structure as that observed in the wild-type complex.

A mu-class glutathione S-transferase from the marine shrimp Litopenaeus vannamei: molecular cloning and active-site structural modeling

A cDNA clone coding for a mu-class glutathione S-transferase (GST) was isolated from a hepatopancreas cDNA library from the shrimp Litopenaeus vannamei. The deduced amino acid sequence (215 amino acids) has >50% identity to rodents and other mammals mu-class GSTs. Using RT-PCR, the shrimp GST transcript was detected in hepatopancreas, hemocytes, gills, and muscle, but not in pleopods. The shrimp GST sequence was computer modeled and found to fit the classical two-domain GST structure. Domain I, containing the glutathione (GSH) binding site, is more conserved compared to the flexible C-terminal domain II. Residue Q208 appears to be a key to substrate specificity by comparison with mammalian GST mutants. This position is commonly occupied by serine or threonine in mammalian mu-class GSTs, and shrimp Q208 may affect the affinity to substrates like aminochrome or 1,3-dimethyl-2-cyano-1-nitrosoguanidine. This is the first report of molecular cloning and structural modeling of a crustacean GST and provides new insights into the nature of the detoxification response on marine invertebrates.

Structure-function relationships in glutathione and its analogues

The results are presented of measurements of protonation constants (potentiometry and NMR), UV spectroscopic properties and redox potentials of GSH and its five analogues, which are modified at the C-terminal glycine residue (gammaGlu-Cys-X, X = Gly, Gly-NH2, Gly-OEt, Ala, Glu, Ser). Strong linear correlations were found between various properties of the thiol and other functions of these peptides. These results allow discussion of the relationships between the structures and properties in glutathione and its analogues, and provide a novel chemical background for the issue of control of GSH reactivity.

Exploration of in vitro pro-drug activation and futile cycling by glutathione S-transferases: thiol ester hydrolysis and inhibitor maturation

In addition to glutathione (GSH) conjugating activity, glutathione S-transferases (GSTs) catalyze "reverse" reactions, such as the hydrolysis of GSH thiol esters. Reverse reactions are of interest as potential tumor-directed pro-drug activation strategies and as mechanisms for tissue redistribution of carboxylate-containing drugs. However, the mechanism and specificity of GST-mediated GSH thiol ester hydrolysis are uncharacterized. Here, the GSH thiol esters of ethacrynic acid (E-SG) and several nonsteroidal antiinflammatory agents have been tested as substrates with human GSTs. The catalytic hydrolysis of these thiol esters appears to be a general property of GSTs. The hydrolysis of the thiol ester of E-SG was studied further with GSTA1-1 and GSTP1-1, as a model pro-drug with several possible fates for the hydrolysis products: competitive inhibition, covalent enzyme adduction, and sequential metabolism. In contrast to hydrolysis rates, significant isoform-dependent differences in the subsequent fate of the products ethacrynic acid and GSH were observed. At low [E-SG], only the GSTP1-1 efficiently catalyzed sequential metabolism, via a dissociative mechanism.

Structure, function and evolution of glutathione transferases: implications for classification of non-mammalian members of an ancient enzyme superfamily

The glutathione transferases (GSTs; also known as glutathione S-transferases) are major phase II detoxification enzymes found mainly in the cytosol. In addition to their role in catalysing the conjugation of electrophilic substrates to glutathione (GSH), these enzymes also carry out a range of other functions. They have peroxidase and isomerase activities, they can inhibit the Jun N-terminal kinase (thus protecting cells against H(2)O(2)-induced cell death), and they are able to bind non-catalytically a wide range of endogenous and exogenous ligands. Cytosolic GSTs of mammals have been particularly well characterized, and were originally classified into Alpha, Mu, Pi and Theta classes on the basis of a combination of criteria such as substrate/inhibitor specificity, primary and tertiary structure similarities and immunological identity. Non-mammalian GSTs have been much less well characterized, but have provided a disproportionately large number of three-dimensional structures, thus extending our structure-function knowledge of the superfamily as a whole. Moreover, several novel classes identified in non-mammalian species have been subsequently identified in mammals, sometimes carrying out functions not previously associated with GSTs. These studies have revealed that the GSTs comprise a widespread and highly versatile superfamily which show similarities to non-GST stress-related proteins. Independent classification systems have arisen for groups of organisms such as plants and insects. This review surveys the classification of GSTs in non-mammalian sources, such as bacteria, fungi, plants, insects and helminths, and attempts to relate them to the more mainstream classification system for mammalian enzymes. The implications of this classification with regard to the evolution of GSTs are discussed.

Inactivation of glutathione S-transferases by nitric oxide-derived oxidants: exploring a role for tyrosine nitration

Reactive intermediates derived from nitric oxide ((*)NO) are thought to play a contributing role in disease states associated with inflammation and infection. We show here that glutathione S-transferases (GSTs), principal enzymes responsible for detoxification of endogenous and exogenous electrophiles, are susceptible to inactivation by reactive nitrogen species (RNS). Treatment of isolated GSTs or rat liver homogenates with either peroxynitrite, the myeloperoxidase/hydrogen peroxide/nitrite system, or tetranitromethane, resulted in loss of GST activity with a concomitant increase in the formation of protein-associated 3-nitrotyrosine (NO(2)Tyr). This inactivation was only partially (<25%) reversible by dithiothreitol, and exposure of GSTs to hydrogen peroxide or S-nitrosoglutathione was only partially inhibitory (<25%) and did not result in protein nitration. Thus, irreversible modifications such as tyrosine nitration may have contributed to GST inactivation by RNS. Since all GSTs contain a critical, highly conserved, active-site tyrosine residue, we postulated that this Tyr residue might present a primary target for nitration by RNS, thus leading to enzyme inactivation. To directly investigate this possibility, we analyzed purified mouse liver GST-mu, following nitration by several RNS, by trypsin digestion, HPLC separation, and matrix-assisted laser desorption/ionization-time of flight analysis, to determine the degree of tyrosine nitration of individual Tyr residues. Indeed, nitration was found to occur preferentially on several tyrosine residues located in and around the GST active site. However, RNS concentrations that resulted in near complete GST inactivation only caused up to 25% nitration of even preferentially targeted tyrosine residues. Hence, nitration of active-site tyrosine residues may contribute to GST inactivation by RNS, but is unlikely to fully account for enzyme inactivation. Overall, our studies illustrate a potential mechanism by which RNS may promote (oxidative) injury by environmental pollutants in association with inflammation.

The C-terminus of glutathione S-transferase A1-1 is required for entropically-driven ligand binding

Binding of a hydrophobic glutathione product conjugate to rGST A1-1 proceeds via a two-step mechanism, including rapid ligand docking, followed by a slow isomerization to the final [GST.ligand] complex, which involves the localization of the flexible C-terminal helix. These kinetically resolved steps have been observed previously by stopped-flow fluorescence with the wild-type rGST A1-1, which contains a native Trp-21 approximately 20 A from the ligand binding site at the intrasubunit domain-domain interface. To confirm this binding mechanism, as well as elucidate the effects of truncation of the C-terminus, we have further characterized the binding and dissociation of the glutathione-ethacrynic acid product conjugate (GS-EA) to wild-type, F222W:W21F, and Delta209-222 rGST A1-1 and wild-type hGST A1-1. Although modest kinetic differences were observed between the hGST A1-1 and rGST A1-1, stopped-flow binding studies with GS-EA verified that the two-step mechanism of ligand binding is not unique to the GST A1-1 isoform from rat. An F222W:W21F rGST A1-1 double mutant provides a direct fluorescence probe of changes in the environment of the C-terminal residue. The observation of two relaxation times during ligand binding and dissociation to F222W:W21F suggests that the C-terminus has an intermediate conformation following ligand docking, which is distinct from its conformation in the apoenzyme or localized helical state. For the wild-type, Delta209-222, and F222W:W21F proteins, variable-temperature stopped-flow experiments were performed and activation parameters calculated for the individual steps of the binding reaction. Activation parameters for the binding reaction coordinate illustrate that the C-terminus provides a significant entropic contribution to ligand binding, which is completely realized within the initial docking step of the binding mechanism. In contrast, the slow isomerization step is enthalpically driven. The partitioning of entropic and enthalpic components of binding energy was confirmed by isothermal titration calorimetry with wild-type and Delta209-222 rGST A1-1.

A highly acidic tyrosine 9 and a normally titrating tyrosine 212 contribute to the catalytic mechanism of human glutathione transferase A4-4

Human glutathione transferase A4-4 is an enzyme catalyzing the detoxication of intracellularly produced electrophiles such as 4-hydroxynonenal and other alkenal products of lipid peroxidation. Two tyrosines in the active site of the enzyme have been studied with help of UV difference spectroscopy and site-directed mutagenesis. The titration curve of GST A4-4 shows a pK(a) of 6.7 attributable to tyrosine 9, which in the Y212F mutant was shifted to pK(a) 7.1. In both cases the pK(a) was independent of the absence or presence of GSH. Thus, the active-site tyrosine 9 of this isoenzyme is more than one unit more acidic than the corresponding tyrosine of other Alpha class glutathione transferases. The tyrosines remaining in the Y9F mutant titrate like free tyrosine with pK(a) values > or = 10. A mechanism involving a tyrosine-9-bound water molecule acting as a proton shuttle is proposed for the Michael additions catalyzed by GST A4-4.

Tyr115, gln165 and trp209 contribute to the 1, 2-epoxy-3-(p-nitrophenoxy)propane-conjugating activity of glutathione S-transferase cGSTM1-1

We investigated the epoxidase activity of a class mu glutathione S-transferase (cGSTM1-1), using 1,2-epoxy-3-(p-nitrophenoxy)propane (EPNP) as substrate. Trp209 on the C-terminal tail, Arg107 on the alpha4 helix, Asp161 and Gln165 on the alpha6 helix of cGSTM1-1 were selected for mutagenesis and kinetic studies. A hydrophobic side-chain at residue 209 is needed for the epoxidase activity of cGSTM1-1. Replacing Trp209 with histidine, isoleucine or proline resulted in a fivefold to 28-fold decrease in the k(cat)(app) of the enzyme, while a modest 25 % decrease in the k(cat)(app) was observed for the W209F mutant. The rGSTM1-1 enzyme has serine at the correponding position. The k(cat)(app) of the S209W mutant is 2. 5-fold higher than that of the wild-type rGSTM1-1. A charged residue is needed at position 107 of cGSTM1-1. The K(m)(app)(GSH) of the R107L mutant is 38-fold lower than that of the wild-type enzyme. On the contrary, the R107E mutant has a K(m)(app)(GSH) and a k(cat)(app) that are 11-fold and 35 % lower than those of the wild-type cGSTM1-1. The substitutions of Gln165 with Glu or Leu have minimal effect on the affinity of the mutants towards GSH or EPNP. However, a discernible reduction in k(cat)(app) was observed. Asp161 is involved in maintaining the structural integrity of the enzyme. The K(m)(app)(GSH) of the D161L mutant is 616-fold higher than that of the wild-type enzyme. In the hydrogen/deuterium exchange experiments, this mutant has the highest level of deuteration among all the proteins tested. We also elucidated the structure of cGSTM1-1 co-crystallized with the glutathionyl-conjugated 1, 2-epoxy-3-(p-nitrophenoxy)propane (EPNP) at 2.8 A resolution. The product found in the active site was 1-hydroxy-2-(S-glutathionyl)-3-(p-nitrophenoxy)propane, instead of the conventional 2-hydroxy isomer. The EPNP moiety orients towards Arg107 and Gln165 in dimer AB, and protrudes into a hydrophobic region formed by the loop connecting beta1 and alpha1 and part of the C-terminal tail in dimer CD. The phenoxyl ring forms strong ring stacking with the Trp209 side-chain in dimer CD. We hypothesize that these two conformations represent the EPNP moiety close to the initial and final stages of the reaction mechanism, respectively.

Noncovalent associations of glutathione S-transferase and ligands: a study using electrospray quadrupole/time-of-flight mass spectrometry

Human glutathione S-transferase A1-1 was observed predominantly as dimeric ions (51 kDa) during electrospray mass spectrometric analysis from aqueous solution at pH 7.4, in keeping with the known dimeric structure in solution. When analyses were performed on solutions of the enzyme containing glutathione (GSH), noncovalent adducts of protein dimer and one or two ligand molecules were observed; each mass increment, which exceeded the mass of GSH alone, was provisionally interpreted to indicate concomitant association of two water molecules per bound GSH. Noncovalent adducts of ligand and protein dimer were similarly observed for oxidized glutathione and for two glutathione inhibitors, both incorporating substituted thiol structures. In these instances, the mass increments exactly matched the ligand masses, suggesting that the apparent concomitant binding of water was associated with the presence in the ligand of a free thiol group. Collisionally activated decomposition during tandem mass spectrometry analyses of noncovalent adducts incorporating protein dimer and ligands yielded initially the denuded dimer; at higher collision energies the monomer and a protein fragment were formed.

Ab initio calculations on hidden modulators of theta class glutathione transferase activity

The glutathione transferases decrease the pKa of glutathione, allowing its deprotonation and the formation of the more reactive thiolate anion. The thiolate is maintained in the active site through a weak conventional hydrogen bond first sphere interaction donated by a Tyr hydroxyl in the Alpha, Mu, Pi, and Sigma glutathione transferase classes that can be modified by other second sphere or indirect thiolate contacts. However, the Theta and Delta class isoforms use a Ser hydroxyl for stabilizing the GSH thiolate, and as such, have a different chemical system compared with that of the Tyr possessed by other classes. We have used high level ab initio methods to investigate this interaction by using a simple methanol methanethiol system as a model. The hydrogen bond strength of this initial first sphere interaction was calculated to be less than that of the Tyr interaction. A putative second sphere interaction exists in the Theta and Delta class structures between Cys or Ser-14 and Ser-11 in the mammalian Theta subclass 1 and 2, respectively. The effect of this interaction on the first sphere interaction has also been investigated and found to significantly increase the energy of the bond.

Substrate- and pH-dependent contribution of oxyanion binding site to the catalysis of prolyl oligopeptidase, a paradigm of the serine oligopeptidase family

Prolyl oligopeptidase, an enzyme implicated in memory disorders, is a member of a new serine peptidase family. Crystallographic studies (Fülöp et al., 1998) revealed a novel oxyanion binding site containing a tyrosine residue, Tyr473. To study the importance of Tyr473 OH, we have produced prolyl oligopeptidase and its Tyr473Phe variant in Escherichia coli. The specificity rate constant, k(cat)/Km, for the modified enzyme decreased by a factor of 8-40 with highly specific substrates, Z-Gly-Pro-Nap, and a fluorogenic octapeptide. With these compounds, the decline in k(cat) was partly compensated for by reduction in Km, a difference from the extensively studied subtilisin. With the less specific suc-Gly-Pro-Nap, the Km value, which approximates Ks, was not significantly changed, resulting in greater diminution (approximately 500-fold) in k(cat)/Km. The second-order rate constant for the reaction with Z-Pro-prolinal, a slow tight-binding transition-state analogue inhibitor, and the Ki values for a slow substrate and two product-like inhibitors were not significantly affected by the Tyr473 OH group. The mechanism of transition-state stabilization was markedly dependent upon the nature of substrate and varied with pH as the enzyme interconverted between its two catalytically competent forms.

N5-(L-1-carboxyethyl)-L-ornithine synthase: physical and spectral characterization of the enzyme and its unusual low pKa fluorescent tyrosine residues

N5-(L-1-carboxyethyl)-L-ornithine synthase [E.C. 1.5.1.24] (CEOS) from Lactococcus lactis has been cloned, expressed, and purified from Escherichia coli in quantities sufficient for characterization by biophysical methods. The NADPH-dependent enzyme is a homotetramer (Mr approximately equal to 140,000) and in the native state is stabilized by noncovalent interactions between the monomers. The far-ultraviolet circular dichroism spectrum shows that the folding pattern of the enzyme is typical of the alpha,beta family of proteins. CEOS contains one tryptophan (Trp) and 19 tyrosines (Tyr) per monomer, and the fluorescence spectrum of the protein shows emission from both Trp and Tyr residues. Relative to N-acetyltyrosinamide, the Tyr quantum yield of the native enzyme is about 0.5. All 19 Tyr residues are titratable and, of these, two exhibit the uncommonly low pKa of approximately 8.5, 11 have pKa approximately 10.75, and the remaining six titrate with pKa approximately 11.3. The two residues with pKa approximately 8.5 contribute approximately 40% of the total tyrosine emission, implying a relative quantum yield >1, probably indicating Tyr-Tyr energy transfer. In the presence of NADPH, Tyr fluorescence is reduced by 40%, and Trp fluorescence is quenched completely. The latter result suggests that the single Trp residue is either at the active site, or in proximity to the sequence GSGNVA, that constitutes the beta alphabeta fold of the nucleotide-binding domain. Chymotrypsin specifically cleaves native CEOS after Phe255. Although inactivated by this single-site cleavage of the subunit, the enzyme retains the capacity to bind NADPH and tetramer stability is maintained. Possible roles in catalysis for the chymotrypsin sensitive loop and for the low pKa Tyr residues are discussed.

Affinity labeling of pig lung glutathione S-transferase pi by 4-(fluorosulfonyl)benzoic acid

The compound 4-(fluorosulfonyl)benzoic acid (4-FSB) functions as an affinity label of the dimeric pig lung pi class glutathione S-transferase yielding a completely inactive enzyme. Protection against inactivation is provided by glutathione-based ligands, suggesting that the reaction target is near or part of the glutathione binding site. Radioactive 4-FSB is incorporated to the extent of 1 mol per mole of enzyme subunit. Peptide mapping revealed that 4-FSB reacts with two tyrosine residues in the ratio 69% Tyr7 and 31% Tyr106. The ratio is not changed by the addition of ligands. The results suggest that only one of the tyrosine residues can be labeled in the active site of a given subunit; i.e., reactions with Tyr7 and Tyr106 are mutually exclusive. We propose that the difference in labeling of these tyrosine residues is related to their pKa values, with Tyr7 exhibiting the lower pKa. The modified enzyme no longer binds to a S-hexylglutathione-agarose affinity column, even when only one of the active sites contains 4-FSB; these results may reflect interaction between the subunits. We conclude that Tyr7 and Tyr106 of the pig lung class pi glutathione S-transferase are important for function and are located at or close to the substrate binding site of the enzyme.

Conformational changes in the crystal structure of rat glutathione transferase M1-1 with global substitution of 3-fluorotyrosine for tyrosine

The structure of the tetradeca-(3-fluorotyrosyl) M1-1 GSH transferase (3-FTyr GSH transferase), a protein in which tyrosine residues are globally substituted by 3-fluorotyrosines has been determined at 2.2 A resolution. This variant was produced to study the effect on the enzymatic mechanism and the structure was undertaken to assess how the presence of the 3-fluorotyrosyl residue influences the protein conformation and hence its function. Although fluorinated amino acid residues have frequently been used in biochemical and NMR investigations of proteins, no structure of a protein that has been globally substituted with a fluorinated amino acid has previously been reported. Thus, this structure represents the first crystal structure of such a protein containing a library of 14 (28 crystallographically distinct) microenvironments from which the nature of the interactions of fluorine atoms with the rest of the protein can be evaluated. Numerous conformational changes are observed in the protein structure as a result of substitution of 3-fluorotyrosine for tyrosine. The results of the comparison of the crystal structure of the fluorinated protein with the native enzyme reveal that conformational changes are observed for most of the 3-fluorotyrosines. The largest differences are seen for residues where the fluorine, the OH, or both are directly involved in interactions with other regions of the protein or with a symmetry-related molecule. The fluorine atoms of the 3-fluorotyrosine interact primarily through hydrogen bonds with other residues and water molecules. In several cases, the conformation of a 3-fluorotyrosine is different in one of the monomers of the enzyme from that observed in the other, including different hydrogen-bonding patterns. Altered conformations can be related to differences in the crystal packing interactions of the two monomers in the asymmetric unit. The fluorine atom on the active-site Tyr6 is located near the S atom of the thioether product (9R,10R)-9-(S-glutathionyl)-10-hydroxy-9,10-dihydrophenanthrene and creates a different pattern of interactions between 3-fluorotyrosine 6 and the S atom. Studies of these interactions help explain why 3-FTyr GSH transferase exhibits spectral and kinetic properties distinct from the native GSH transferase.

The three-dimensional structure of an avian class-mu glutathione S-transferase, cGSTM1-1 at 1.94 A resolution

Glutathione S-transferase cGSTM1-1, an avian class-mu enzyme with high sequence identity with rGSTM3-3, was expressed heterologously in Escherichia coli. The three-dimensional structure of this protein that co-crystallized with an inhibitor, S-hexylglutathione, was determined by the molecular replacement method and refined to 1.94 A resolution. The three-dimensional structure and the folding topology of the dimeric cGSTM1-1 closely resembles those of other class-mu GSTs. The bound inhibitor, S-hexylglutathione, orients in disparate directions in the two subunits. The combined space occupied by the hexyl moiety of the inhibitors overlaps with that reported for rGSTM1-1 co-crystallized with (9 S,10 S)-9-(S-glutathionyl)-10-hydroxy-9,10-dihydrophenanthrene. Conformational differences at a flexible loop (residue 35 to 40) were also observed between the crystal structures of cGSTM1-1 and rGSTM1-1.cGSTM1-1 has the highest epoxidase activity among all the class-mu enzymes reported. Tyr115, has been identified as a residue that participates in the epoxidase activity of class-mu glutathione S-transferase and is conserved in cGSTM1-1. The epoxidase and trans-4-phenyl-3-buten-2-one conjugating activity of cGSTM1-1 are decreased drastically but not abolished by replacing Tyr115 with phenylalanine. The specificity constant of the cGSTM1-1(Y115F) mutant, with 1-chloro-2,4-dinitrobenzene as substrate, is 15-fold higher than that of the wild-type enzyme.

Structures of class pi glutathione S-transferase from human placenta in complex with substrate, transition-state analogue and inhibitor

BACKGROUND

Glutathione S-transferases (GSTs) are detoxification enzymes, found in all aerobic organisms, which catalyse the conjugation of glutathione with a wide range of hydrophobic electrophilic substrates, thereby protecting the cell from serious damage caused by electrophilic compounds. GSTs are classified into five distinct classes (alpha, mu, pi, sigma and theta) by their substrate specificity and primary structure. Human GSTs are of interest because tumour cells show increased levels of expression of single classes of GSTs, which leads to drug resistance. Structural differences between classes of GST can therefore be utilised to develop new anti-cancer drugs. Many mutational and structural studies have been carried out on the mu and alpha classes of GST to elucidate the reaction mechanism, whereas knowledge about the pi class is still limited.

RESULTS

We have solved the structures of the pi class GST hP1-1 in complex with its substrate, glutathione, a transition-state complex, the Meisenheimer complex, and an inhibitor, S-(rho-bromobenzyl)-glutathione, and refined them to resolutions of 1.8 A, 2.0 A and 1.9 A, respectively. All ligand molecules are well-defined in the electron density. In all three structures, an additionally bound N-morpholino-ethansulfonic acid molecule from the buffer solution was found.

CONCLUSIONS

In the structure of the GST-glutathione complex, two conserved water molecules are observed, one of which hydrogen bonds directly to the sulphur atom of glutathione and the other forms hydrogen bonds with residues around the glutathione-binding site. These water molecules are absent from the structure of the Meisenheimer complex bound to GST, implicating that deprotonation of the cysteine occurs during formation of the ternary complex which involves expulsion of the inner bound water molecule. The comparison of our structures with known mu class GST structures show differences in the location of the electrophile-binding site (H-site), explaining the different substrate specificities of the two classes. Fluorescence measurements are in agreement with the position of the N-morpholino-ethansulfonic acid, close to Trp28, identifying a possible ligandin-substrate binding site.

Characterising non-covalent interactions with the Cambridge Structural Database

This review describes how the CSD can be used to study non-covalent interactions. Several different types of information may be obtained. First, the relative frequencies of various interactions can be studied; for example, we have shown that the terminal oxygen atoms of phosphate groups accept hydrogen bonds far more often than the linkage oxygens. Secondly, information can be obtained about the geometries of nonbonded contacts; for example, hydrogen bonds to P-O groups rarely form along the extension of the P-O bond, whereas short contacts between oxygen and carbon-bound iodine show a strong preference for linear C-I ... O angles. Thirdly, the CSD can be searched for novel interactions which may be exploited in inhibitor design; for example, the I ... O contacts just mentioned, and N-H ... pi hydrogen bonds. Finally, the CSD can suggest synthetic targets for medicinal chemistry; for example, molecules containing delocalised electron deficient groups such as trimethylammonium, pyridinium, thaizolium and dinitrophenyl have a good chance of binding to an active-site tryptophan. Although the CSD contains small-molecule crystal structures, not protein-ligand complexes, there is considerable evidence that the contacts seen in the two types of structures are similar. We have illustrated this a number of times in the present review and additional evidence has been given previously by Klebe. The major advantages of the CSD are its size, diversity and experimental accuracy. For these reasons, it is a useful tool for modellers engaged in rational inhibitor design.