Evidence for an induced-fit mechanism operating in pi class glutathione transferases

Observing How Glutathione and S-Hexyl Glutathione Bind to Glutathione S-Transferase from Rhipicephalus (Boophilus) microplus

Rhipicephalus (Boophilus) microplus is one of the most widespread ticks causing a massive loss to livestock production. The long-term use of acaracides rapidly develops acaracide resistance. In R. microplus, enhancing the metabolic activity of glutathione S-transferase (RmGST) is one of the mechanisms underlying acaracide resistance. RmGST catalyzes the conjugation of glutathione (GSH) to insecticides causing an easy-to-excrete conjugate. The active RmGST dimer contains two active sites (hydrophobic co-substrate binding site (H-site) and GSH binding site (G-site)) in each monomer. To preserve the insecticide efficacy, s-hexyl glutathione (GTX), a GST inhibitor, has been used as a synergist. To date, no molecular information on the RmGST-GSH/GTX complex is available. The insight is important for developing a novel RmGST inhibitor. Therefore, in this work, molecular dynamics simulations (MD) were performed to explore the binding of GTX and GSH to RmGST. GSH binds tighter and sits rigidly inside the G-site, while flexible GTX occupies both active sites. In GSH, the backbone mainly interacts with W8, R43, W46, K50, N59, L60, Q72, and S73, while its thiol group directs to Y7. In contrast, the aliphatic hexyl of GTX protrudes into the H-site and allows a flexible peptide core to form various interactions. Such high GTX flexibility and the protrusion of its hexyl moiety to the H-site suggest the dual role of GTX in preventing the conjugation reaction and the binding of acaracide. This insight can provide a better understanding of an important insecticide-resistance mechanism, which may in turn facilitate the development of novel approaches to tick control.

A structure-based mechanism of cisplatin resistance mediated by glutathione transferase P1-1

The resurgence of platinum-based chemotherapy in the last few years has renewed interest in the field, including clinical studies of cisplatin in combination with resistance modulators. Indeed, cisplatin is one of the most successful anticancer agents, effective against a wide range of solid tumors. However, its use is restricted by side effects and/or by intrinsic or acquired drug resistance. We propose here a new mechanism of cisplatin resistance mediated by glutathione transferase (GST) P1-1, as a cisplatin-binding protein. Our results show that cisplatin can be inactivated by this protein with the aid of 2 solvent-accessible and reactive cysteines. These findings may constitute the basis for the design and synthesis of new GST inhibitors able to circumvent cisplatin resistance.

Cisplatin [cis-diamminedichloroplatinum(II) (cis-DDP)] is one of the most successful anticancer agents effective against a wide range of solid tumors. However, its use is restricted by side effects and/or by intrinsic or acquired drug resistance. Here, we probed the role of glutathione transferase (GST) P1-1, an antiapoptotic protein often overexpressed in drug-resistant tumors, as a cis-DDP–binding protein. Our results show that cis-DDP is not a substrate for the glutathione (GSH) transferase activity of GST P1-1. Instead, GST P1-1 sequesters and inactivates cisplatin with the aid of 2 solvent-accessible cysteines, resulting in protein subunits cross-linking, while maintaining its GSH-conjugation activity. Furthermore, it is well known that GST P1-1 binding to the c-Jun N-terminal kinase (JNK) inhibits JNK phosphorylation, which is required for downstream apoptosis signaling. Thus, in turn, GST P1-1 overexpression and Pt-induced subunit cross-linking could modulate JNK apoptotic signaling, further confirming the role of GST P1-1 as an antiapoptotic protein.

Arabidopsis thaliana dehydroascorbate reductase 2: Conformational flexibility during catalysis

Dehydroascorbate reductase (DHAR) catalyzes the glutathione (GSH)-dependent reduction of dehydroascorbate and plays a direct role in regenerating ascorbic acid, an essential plant antioxidant vital for defense against oxidative stress. DHAR enzymes bear close structural homology to the glutathione transferase (GST) superfamily of enzymes and contain the same active site motif, but most GSTs do not exhibit DHAR activity. The presence of a cysteine at the active site is essential for the catalytic functioning of DHAR, as mutation of this cysteine abolishes the activity. Here we present the crystal structure of DHAR2 from Arabidopsis thaliana with GSH bound to the catalytic cysteine. This structure reveals localized conformational differences around the active site which distinguishes the GSH-bound DHAR2 structure from that of DHAR1. We also unraveled the enzymatic step in which DHAR releases oxidized glutathione (GSSG). To consolidate our structural and kinetic findings, we investigated potential conformational flexibility in DHAR2 by normal mode analysis and found that subdomain mobility could be linked to GSH binding or GSSG release.

Mice lacking three Loci encoding 14 glutathione transferase genes: a novel tool for assigning function to the GSTP, GSTM, and GSTT families

Glutathione S-transferases (GSTs) form a superfamily defined by their ability to catalyze the conjugation of glutathione with electrophilic substrates. These enzymes are proposed to play a critical role in protection of cellular components from damage mediated by reactive metabolites. Twenty-two cytosolic GSTs, grouped into seven families, are recognized in mice. This complexity hinders the assignment of function to a subset or family of these genes. We report generation of a mouse line in which the locus encoding three GST gene families is deleted. This includes the four Gstt genes spanning 65 kb on chromosome 10 and the seven Gstm genes found on a 150 kb segment of DNA chromosome 3. In addition, we delete two Gstp genes on chromosome 19 as well as a third related gene located 15 kb telomeric to Gstp1 and Gstp2, which we identify as a potential new member of this gene family. We show that, despite the loss of up to 75% of total GST activity in some tissues from these animals, the mice are healthy and fertile, with normal life expectancy. The normal development and health of these animals make them an appropriate model for defining the role of these families in redox homeostasis and metabolism of drugs and environmental pollutants.

S-nitrosation of glutathione transferase p1-1 is controlled by the conformation of a dynamic active site helix

S-Nitrosation is a post-translational modification of protein cysteine residues, which occurs in response to cellular oxidative stress. Although it is increasingly being linked to physiologically important processes, the molecular basis for protein regulation by this modification remains poorly understood. We used transient kinetic methods to determine a minimal mechanism for spontaneous S-nitrosoglutathione (GSNO)-mediated transnitrosation of human glutathione transferase (GST) P1-1, a major detoxification enzyme and key regulator of cell proliferation. Cys(47) of GSTP1-1 is S-nitrosated in two steps, with the chemical step limited by a pre-equilibrium between the open and closed conformations of helix α2 at the active site. Cys(101), in contrast, is S-nitrosated in a single step but is subject to negative cooperativity due to steric hindrance at the dimer interface. Despite the presence of a GSNO binding site at the active site of GSTP1-1, isothermal titration calorimetry as well as nitrosation experiments using S-nitrosocysteine demonstrate that GSNO binding does not precede S-nitrosation of GSTP1-1. Kinetics experiments using the cellular reductant glutathione show that Cys(101)-NO is substantially more resistant to denitrosation than Cys(47)-NO, suggesting a potential role for Cys(101) in long term nitric oxide storage or transfer. These results constitute the first report of the molecular mechanism of spontaneous protein transnitrosation, providing insight into the post-translational control of GSTP1-1 as well as the process of protein transnitrosation in general.

Identifying disordered regions in proteins by limited proteolysis

Limited proteolysis experiments can be successfully used to detect sites of disorder in otherwise folded globular proteins. The approach relies on the fact that the proteolysis of a polypeptide substrate requires its binding in an extended conformation at the protease's active site and thus an enhanced backbone flexibility or local unfolding of the site of proteolytic attack. A striking correlation was found between sites of limited proteolysis and sites of enhanced chain flexibility of the polypeptide chain, this last evaluated by the crystallographically determined B-factor. In numerous cases, it has been shown that limited proteolysis occurs at chain regions characterized by missing electron density and thus being disordered. Therefore, limited proteolysis is a simple and reliable experimental technique that can detect sites of disorder in proteins, thus complementing the results that can be obtained by the use of other physicochemical and computational approaches.

Crystal structures and kinetic studies of human Kappa class glutathione transferase provide insights into the catalytic mechanism

GSTs (glutathione transferases) are a family of enzymes that primarily catalyse nucleophilic addition of the thiol of GSH (reduced glutathione) to a variety of hydrophobic electrophiles in the cellular detoxification of cytotoxic and genotoxic compounds. GSTks (Kappa class GSTs) are a distinct class because of their unique cellular localization, function and structure. In the present paper we report the crystal structures of hGSTk (human GSTk) in apo-form and in complex with GTX (S-hexylglutathione) and steady-state kinetic studies, revealing insights into the catalytic mechanism of hGSTk and other GSTks. Substrate binding induces a conformational change of the active site from an ‘open’ conformation in the apo-form to a ‘closed’ conformation in the GTX-bound complex, facilitating formations of the G site (GSH-binding site) and the H site (hydrophobic substrate-binding site). The conserved Ser16 at the G site functions as the catalytic residue in the deprotonation of the thiol group and the conserved Asp69, Ser200, Asp201 and Arg202 form a network of interactions with γ-glutamyl carboxylate to stabilize the thiolate anion. The H site is a large hydrophobic pocket with conformational flexibility to allow the binding of different hydrophobic substrates. The kinetic mechanism of hGSTk conforms to a rapid equilibrium random sequential Bi Bi model.

Diuretic drug binding to human glutathione transferase P1-1: potential role of Cys-101 revealed in the double mutant C47S/Y108V

The diuretic drug ethacrynic acid (EA), both an inhibitor and substrate of pi class glutathione S-transferase (GST P1-1), has been tested in clinical trials as an adjuvant in chemotherapy. We recently studied the role of the active site residue Tyr-108 in binding EA to the enzyme and found that the analysis was complicated by covalent binding of this drug to the highly reactive Cys-47. Previous attempts to eliminate this binding by chemical modification yielded ambiguous results and therefore we decided here to produce a double mutant C47S/Y108V by site directed mutagenesis and further expression in Escherichia coli and the interaction of EA and its GSH conjugate (EASG) examined by calorimetric studies and X-ray diffraction. Surprisingly, in the absence of Cys-47, Cys-101 (located at the dimer interface) becomes a target for modification by EA, albeit at a lower conjugation rate than Cys-47. The Cys-47 → Ser mutation in the double mutant enzyme induces a positive cooperativity between the two subunits when ligands with affinity to G-site bind to enzyme. However, this mutation does not seem to affect the thermodynamic properties of ligand binding to the electrophilic binding site (H-site) and the thermal or chemical stability of this double mutant does not significantly affect the unfolding mechanism in either the absence or presence of ligand. Crystal structures of apo and an EASG complex are essentially identical with a few exceptions in the H-site and in the water network at the dimer interface.

Class Pi glutathione transferase unfolds via a dimeric and not monomeric intermediate: functional implications for an unstable monomer

Cytosolic class pi glutathione transferase P1-1 (GSTP1-1) is associated with drug resistance and proliferative pathways because of its catalytic detoxification properties and ability to bind and regulate protein kinases. The native wild-type protein is homodimeric, and whereas the dimeric structure is required for catalytic functionality, a monomeric and not dimeric form of class pi GST is reported to mediate its interaction with and inhibit the activity of the pro-apoptotic enzyme c-Jun N-terminal kinase (JNK) [Adler, V., et al. (1999) EMBO J. 18, 1321-1334]. Thus, the existence of a stable monomeric form of wild-type class pi GST appears to have physiological relevance. However, there are conflicting accounts of the subunit's intrinsic stability since it has been reported to be either unstable [Dirr, H., and Reinemer, P. (1991) Biochem. Biophys. Res. Commun. 180, 294-300] or stable [Aceto, A., et al. (1992) Biochem. J. 285, 241-245]. In this study, the conformational stability of GSTP1-1 was re-examined by equilibrium folding and unfolding kinetics experiments. The data do not demonstrate the existence of a stable monomer but that unfolding of hGSTP1-1 proceeds via an inactive, nativelike dimeric intermediate in which the highly dynamic helix 2 is unfolded. Furthermore, molecular modeling results indicate that a dimeric GSTP1-1 can bind JNK. According to the available evidence with regard to the stability of the monomeric and dimeric forms of GSTP1-1 and the modality of the GST-JNK interaction, formation of a complex between GSTP1-1 and JNK most likely involves the dimeric form of the GST and not its monomer as is commonly reported.

Buffer interference with protein dynamics: a case study on human liver fatty acid binding protein

Selection of suitable buffer types is often a crucial step for generating appropriate protein samples for NMR and x-ray crystallographic studies. Although the possible interaction between MES buffer (2-(N-morpholino)ethanesulfonic acid) and proteins has been discussed previously, the interaction is usually thought to have no significant effects on the structures of proteins. In this study, we demonstrate the direct, albeit weak, interaction between MES and human liver fatty acid binding protein (hLFABP). Rather than affecting the structure of hLFABP, we found that the dynamics of hLFABP, which were previously proposed to be relevant to its functions, were significantly affected by the binding of hLFABP with MES. Buffer interference with protein dynamics was also demonstrated with Bis-Tris buffer, which is quite different from MES and fatty acids in terms of their molecular structures and properties. This result, to our knowledge, is the first published report on buffer interference with protein dynamics on a microsecond to millisecond timescale and could represent a generic problem in the studies of functionally relevant protein dynamics. Although being a fortuity, our finding of buffer-induced changes in protein dynamics offers a clue to how hLFABP accommodates its ligands.

Nested Arg-specific bifunctional crosslinkers for MS-based structural analysis of proteins and protein assemblies

The combination of chemical probing and high-resolution mass spectrometry constitutes a powerful alternative for the structural elucidation of biomolecules possessing unfavorable size, solubility, and flexibility. We have developed nested Arg-specific bifunctional crosslinkers to obtain complementary information to typical Cys- and Lys-specific reagents available on the market. The structures of 1,4-phenyl-diglyoxal (PDG) and 4,4'-biphenyl-diglyoxal (BDG) include two identical 1,2-dicarbonyl functions capable of reacting with the guanido group of Arg residues in proteins, as well as the base-pairing face of guanine in nucleic acids. The reactive functions are separated by modular spacers consisting of one or two benzene rings, which confer greater rigidity to the crosslinker structure than it is afforded by typical aliphatic spacers. Analysis by electrospray ionization (ESI) Fourier transform ion cyclotron resonance (FTICR) mass spectrometry has shown that the probes provide both mono- and bifunctional products with model protein substrates, which are stabilized by the formation of diester derivatives in the presence of borate buffer. The identification of crosslinked sites was accomplished by employing complementary proteolytic procedures and peptide mapping by ESI-FTICR. The results showed excellent correlation with the solvent accessibility and structural context of susceptible residues, and highlighted the significance of possible dynamic effects in determining the outcome of crosslinking reactions. The application of nested reagents with different spacing has provided a new tool for experimentally recognizing flexible regions that may be involved in prominent dynamics in solution. The development of new bifunctional crosslinkers with diverse target specificity and different bridging spans is expected to facilitate the structure elucidation of progressively larger biomolecular assemblies by increasing the number and diversity of spatial constraints available for triangulating the position of crosslinked structures in the three dimensions.

Modulating catalytic activity by unnatural amino acid residues in a GSH-binding loop of GST P1-1

The loop following helix alpha2 in glutathione transferase P1-1 has two conserved residues, Cys48 and Tyr50, important for glutathione (GSH) binding and catalytic activity. Chemical modification of Cys48 thwarts the catalytic activity of the enzyme, and mutation of Tyr50 generally decreases the k(cat) value and the affinity for GSH in a differential manner. Cys48 and Tyr50 were targeted by site-specific mutations and chemical modifications in order to investigate how the alpha2 loop modulates GSH binding and catalysis. Mutation of Cys48 into Ala increased K(M)(GSH) 24-fold and decreased the binding energy of GSH by 1.5 kcal/mol. Furthermore, the protein stability against thermal inactivation and chemical denaturation decreased. The crystal structure of the Cys-free variant was determined, and its similarity to the wild-type structure suggests that the mutation of Cys48 increases the flexibility of the alpha2 loop rather than dislocating the GSH-interacting residues. On the other hand, replacement of Tyr50 with Cys, producing mutant Y50C, increased the Gibbs free energy of the catalyzed reaction by 4.8 kcal/mol, lowered the affinity for S-hexyl glutathione by 2.2 kcal/mol, and decreased the thermal stability. The targeted alkylation of Cys50 in Y50C increased the affinity for GSH and protein stability. Characterization of the most active alkylated variants, S-n-butyl-, S-n-pentyl-, and S-cyclobutylmethyl-Y50C, indicated that the affinity for GSH is restored by stabilizing the alpha2 loop through positioning of the key residue into the lock structure of the neighboring subunit. In addition, k(cat) can be further modulated by varying the structure of the key residue side chain, which impinges on the rate-limiting step of catalysis.

Calorimetric and structural studies of the nitric oxide carrier S-nitrosoglutathione bound to human glutathione transferase P1-1

The nitric oxide molecule (NO) is involved in many important physiological processes and seems to be stabilized by reduced thiol species, such as S-nitrosoglutathione (GSNO). GSNO binds strongly to glutathione transferases, a major superfamily of detoxifying enzymes. We have determined the crystal structure of GSNO bound to dimeric human glutathione transferase P1-1 (hGSTP1-1) at 1.4 A resolution. The GSNO ligand binds in the active site with the nitrosyl moiety involved in multiple interactions with the protein. Isothermal titration calorimetry and differential scanning calorimetry (DSC) have been used to characterize the interaction of GSNO with the enzyme. The binding of GSNO to wild-type hGSTP1-1 induces a negative cooperativity with a kinetic process concomitant to the binding process occurring at more physiological temperatures. GSNO inhibits wild-type enzyme competitively at lower temperatures but covalently at higher temperatures, presumably by S-nitrosylation of a sulfhydryl group. The C47S mutation removes the covalent modification potential of the enzyme by GSNO. These results are consistent with a model in which the flexible helix alpha2 of hGST P1-1 must move sufficiently to allow chemical modification of Cys47. In contrast to wild-type enzyme, the C47S mutation induces a positive cooperativity toward GSNO binding. The DSC results show that the thermal stability of the mutant is slightly higher than wild type, consistent with helix alpha2 forming new interactions with the other subunit. All these results suggest that Cys47 plays a key role in intersubunit cooperativity and that under certain pathological conditions S-nitrosylation of Cys47 by GSNO is a likely physiological scenario.

The structural roles of a conserved small hydrophobic core in the active site and an ionic bridge in domain I of Delta class glutathione S-transferase

GSTs (glutathione S-transferases; E.C.2.5.1.18) are a supergene family of dimeric multifunctional enzymes that have a major role in detoxification pathways. Using a GST from the mosquito Anopheles dirus (adGSTD4-4), we have characterized the enzymatic and physical properties of Leu-6, Thr-31, Leu-33, Ala-35, Glu-37, Lys-40 and Glu-42. These residues generate two motifs located in the N-terminal domain (domain I) that are functionally conserved across GST classes. The aim of this study was to understand the function of these two motifs. The first motif is a small hydrophobic core in the G-site (glutathione-binding site) wall, and the second motif contains an ionic bridge at the N-terminus of the alpha2 helix and is also part of the G-site. The mutations in the small hydrophobic core appear to have structural effects, as shown by the thermal stability, refolding rate and intrinsic fluorescence differences. In the Delta class GST, interactions form an ionic bridge motif located at the beginning of the alpha2 helix. The data suggest that electrostatic interactions in the alpha2 helix are involved in alpha-helix stabilization, and disruption of this ionic bridge interaction changes the movement of the alpha2-helix region, thereby modulating the interaction of the enzyme with substrates. These results show that the small hydrophobic core and ionic bridge have a major impact on structural stabilization, as well as being required to maintain structural conformation of the enzyme. These structural effects are also transmitted to the active site to influence substrate binding and specificity. Therefore changes in the conformation of the G-site wall in the active site appear to be capable of exerting influences on the tertiary structural organization of the whole GST protein.

Structure-based Prediction of bZIP Partnering Specificity

Predicting protein interaction specificity from sequence is an important goal in computational biology. We present a model for predicting the interaction preferences of coiled-coil peptides derived from bZIP transcription factors that performs very well when tested against experimental protein microarray data. We used only sequence information to build atomic-resolution structures for 1711 dimeric complexes, and evaluated these with a variety of functions based on physics, learned empirical weights or experimental coupling energies. A purely physical model, similar to those used for protein design studies, gave reasonable performance. The results were improved significantly when helix propensities were used in place of a structurally explicit model to represent the unfolded reference state. Further improvement resulted upon accounting for residue-residue interactions in competing states in a generic way. Purely physical structure-based methods had difficulty capturing core interactions accurately, especially those involving polar residues such as asparagine. When these terms were replaced with weights from a machine-learning approach, the resulting model was able to correctly order the stabilities of over 6000 pairs of complexes with greater than 90% accuracy. The final model is physically interpretable, and suggests specific pairs of residues that are important for bZIP interaction specificity. Our results illustrate the power and potential of structural modeling as a method for predicting protein interactions and highlight obstacles that must be overcome to reach quantitative accuracy using a de novo approach. Our method shows unprecedented performance in predicting protein-protein interaction specificity accurately using structural modeling and suggests that predicting coiled-coil interactions generally may be within reach.

Two glutathione S-transferase isoenzymes purified from Bulinus truncatus (Gastropoda: Planorbidae)

We purified and characterized two major glutathione S-transferase isoenzymes (GST2 and GST3) from snail Bulinus truncatus (Mollusca, Gastropoda, Planorbidae) tissue. The Km with respect to 1-chloro-2, 4-dinitrobenzene (CDNB) for both isoenzymes was increased as the pH decreased. Km of both isoenzymes with respect to glutathione (GSH) doubled when the pH was increased from 6.0 to 6.5. Acid inactivated GST2 and GST3 and the two enzymes were almost inactive at pH 3.5. However, they retain the full activity for at least 20 h when incubated at pH between 6.0 and 9.0. The optimum temperature was 45 degrees C for GST2 and 50 degrees C for GST3. The half lifetime at 50 degrees C was 70 min and 45 min for GST2 and GST3 isoenzymes, respectively. Addition of 5 mM GSH to the incubation buffer increased the half life of both isoenzymes more than fourfold. The activation energy for catalyzing the conjugation of CDNB was 1.826 and 3.435 kcal/mol for GST2 and GST3, respectively. I50 values for Cibacron blue, bromosulphophthalein, indocyanine green, hematin and ethacrynic acid were 0.76 microM, 47.9 microM, 7.59 microM, 0.03 microM and 0.79 microM for GST2, and 0.479 microM, 79.4 microM, 89.1 microM, 32.4 microM and 1.15 microM for GST3, respectively. Cibacron blue and indocyanine green were non-competitive inhibitors, while hematin was a mixed inhibitor. Bromosulphophthalein was found to be a competitive inhibitor for GST2 and a mixed inhibitor for GST3.

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.

Side-chain flexibility in protein-ligand binding: the minimal rotation hypothesis

The goal of this work is to learn from nature about the magnitudes of side-chain motions that occur when proteins bind small organic molecules, and model these motions to improve the prediction of protein-ligand complexes. Following analysis of protein side-chain motions upon ligand binding in 63 complexes, we tested the ability of the docking tool SLIDE to model these motions without being restricted to rotameric transitions or deciding which side chains should be considered as flexible. The model tested is that side-chain conformational changes involving more atoms or larger rotations are likely to be more costly and less prevalent than small motions due to energy barriers between rotamers and the potential of large motions to cause new steric clashes. Accordingly, SLIDE adjusts the protein and ligand side groups as little as necessary to achieve steric complementarity. We tested the hypothesis that small motions are sufficient to achieve good dockings using 63 ligands and the apo structures of 20 different proteins and compared SLIDE side-chain rotations to those experimentally observed. None of these proteins undergoes major main-chain conformational change upon ligand binding, ensuring that side-chain flexibility modeling is not required to compensate for main-chain motions. Although more frugal in the number of side-chain rotations performed, this model substantially mimics the experimentally observed motions. Most side chains do not shift to a new rotamer, and small motions are both necessary and sufficient to predict the correct binding orientation and most protein-ligand interactions for the 20 proteins analyzed.

Modeling analysis of GST (glutathione-S-transferases) from Wuchereria bancrofti and Brugia malayi

GST (glutathione S-transferases) are a family of detoxification enzymes that catalyze the conjugation of reduced GSH (glutathione) to xenobiotic (endogenous electrophilic) compounds. GST from Wb (Wuchereria bancrofti) and Bm (Brugia malayi) are significantly different from human GST in sequence and structure. Thus, Wb-GST and Bm-GST are potential chemotherapeutic targets for anti-filarial treatment. Comparison of modeled Wb and Bm GST with human GST show structural difference between them. Analysis of the active site residues for the binding of electrophilic co-substrates provides insight towards the design of parasite specific GST inhibitors.

Structure of the major cytosolic glutathione S-transferase from the parasitic nematode Onchocerca volvulus

Onchocerciasis is a debilitating parasitic disease caused by the filarial worm Onchocerca volvulus. Similar to other helminth parasites, O. volvulus is capable of evading the host's immune responses by a variety of defense mechanisms, including the detoxification activities of the glutathione S-transferases (GSTs). Additionally, in response to drug treatment, helminth GSTs are highly up-regulated, making them tempting targets both for chemotherapy and for vaccine development. We analyzed the three-dimensional x-ray structure of the major cytosolic GST from O. volvulus (Ov-GST2) in complex with its natural substrate glutathione and its competitive inhibitor S-hexylglutathione at 1.5 and 1.8 angstrom resolution, respectively. From the perspective of the biochemical classification, the Ov-GST2 seems to be related to pi-class GSTs. However, in comparison to other pi-class GSTs, in particular to the host's counterpart, the Ov-GST2 reveals significant and unusual differences in the sequence and overall structure. Major differences can be found in helix alpha-2, an important region for substrate recognition. Moreover, the binding site for the electrophilic co-substrate is spatially increased and more solvent-accessible. These structural alterations are responsible for different substrate specificities and will form the basis of parasite-specific structure-based drug design investigations.

Glutathione S-transferase Pi has at least three distinguishable xenobiotic substrate sites close to its glutathione-binding site

Benzyl isothiocyanate (BITC), present in cruciferous vegetables, is an efficient substrate of human glutathione S-transferase P1-1 (hGST P1-1). BITC also acts as an affinity label of hGST P1-1 in the absence of glutathione, yielding an enzyme inactive toward BITC as substrate. As monitored by using BITC as substrate, the dependence of k of inactivation (K(I)) of hGST P1-1 on [BITC] is hyperbolic, with K(I) = 66 +/- 7 microM. The enzyme incorporates 2 mol of BITC/mol of enzyme subunit upon complete inactivation. S-Methylglutathione and 8-anilino-1-naphthalene sulfonate (ANS) each yield partial protection against inactivation and decrease reagent incorporation, whereas S-(N-benzylthiocarbamoyl)glutathione or S-methylglutathione + ANS protects completely. Mapping of proteolytic digests of modified enzyme by using mass spectrometry reveals that Tyr(103) and Cys(47) are modified equally. S-Methylglutathione reduces modification of Cys(47), indicating this residue is at/near the glutathione binding region, whereas ANS decreases modification of Tyr(103), suggesting this residue is at/near the BITC substrate site, which is also near the binding site of ANS. The Y103F and Y103S mutant enzymes were generated, expressed, and purified. Both mutants handle substrate 1-chloro-2,4-dinitrobenzene normally; however, Y103S exhibits a 30-fold increase in K(m) for BITC and binds ANS poorly, whereas Y103F has a normal K(m) for BITC and K(d) for ANS. These results indicate that an aromatic residue at position 103 is essential for the binding of BITC and ANS. This study provides evidence for the existence of a novel xenobiotic substrate site in hGST P1-1, which can be occupied by benzyl isothiocyanate and is distinct from that of monobromobimane and 1-chloro-2,4 dinitrobenzene.

Monobromobimane occupies a distinct xenobiotic substrate site in glutathione S-transferase pi

Monobromobimane (mBBr), functions as a substrate of porcine glutathione S-transferase pi (GST pi): The enzyme catalyzes the reaction of mBBr with glutathione. S-(Hydroxyethyl)bimane, a nonreactive analog of monobromobimane, acts as a competitive inhibitor with respect to mBBr as substrate but does not affect the reaction of GST pi with another substrate, 1-chloro-2,4-dinitrobenzene (CDNB). In the absence of glutathione, monobromobimane inactivates GST pi at pH 7.0 and 25 degrees C as assayed using mBBr as substrate, with a lesser effect on the enzyme's use of CDNB as substrate. These results indicate that the sites occupied by CDNB and mBBr are not identical. Inactivation is proportional to the incorporation of 2 moles of bimane/mole of subunit. Modification of GST pi with mBBr does not interfere with its binding of 8-anilino-1-naphthalene sulfonate, indicating that this hydrophobic site is not the target of monobromobimane. S-Methylglutathione and S-(hydroxyethyl)bimane each yield partial protection against inactivation and decrease reagent incorporation, while glutathionyl-bimane protects completely against inactivation. Peptide analysis after trypsin digestion indicates that mBBr modifies Cys45 and Cys99 equally. Modification of Cys45 is reduced in the presence of S-methylglutathione, indicating that this residue is at or near the glutathione binding region. In contrast, modification of Cys99 is reduced in the presence of S-(hydroxyethyl)bimane, suggesting that this residue is at or near the mBBr xenobiotic substrate binding site. Modification of Cys99 can best be understood by reaction with monobromobimane while it is bound to its xenobiotic substrate site in an alternate orientation. These results support the concept that glutathione S-transferase accomplishes its ability to react with a diversity of substrates in part by harboring distinct xenobiotic substrate sites.

Functional role of the lock and key motif at the subunit interface of glutathione transferase p1-1

The glutathione transferases (GSTs) represent a superfamily of dimeric proteins. Each subunit has an active site, but there is no evidence for the existence of catalytically active monomers. The lock and key motif is responsible for a highly conserved hydrophobic interaction in the subunit interface of pi, mu, and alpha class glutathione transferases. The key residue, which is either Phe or Tyr (Tyr(50) in human GSTP1-1) in one subunit, is wedged into a hydrophobic pocket of the other subunit. To study how an essentially inactive subunit influences the activity of the neighboring subunit, we have generated the heterodimer composed of subunits from the fully active human wild-type GSTP1-1 and the nearly inactive mutant Y50A obtained by mutation of the key residue Tyr(50) to Ala. Although the key residue is located far from the catalytic center, the k(cat) value of mutant Y50A decreased about 1300-fold in comparison with the wild-type enzyme. The decrease of the k(cat) value of the heterodimer by about 27-fold rather than the expected 2-fold in comparison with the wild-type enzyme indicates that the two active sites of the dimeric enzyme work synergistically. Further evidence for cooperativity was found in the nonhyperbolic GSH saturation curves. A network of hydrogen-bonded water molecules, found in crystal structures of GSTP1-1, connects the two active sites and the main chain carbonyl group of Tyr(50), thereby offering a mechanism for communication between the two active sites. It is concluded that a subunit becomes catalytically competent by positioning the key residue of one subunit into the lock pocket of the other subunit, thereby stabilizing the loop following the helix alpha2, which interacts directly with GSH.

Thermodynamic description of the effect of the mutation Y49F on human glutathione transferase P1-1 in binding with glutathione and the inhibitor S-hexylglutathione

The thermodynamics of binding of both the substrate glutathione (GSH) and the competitive inhibitor S-hexylglutathione to the mutant Y49F of human glutathione S-transferase (hGST P1-1), a key residue at the dimer interface, has been investigated by isothermal titration calorimetry and fluorescence spectroscopy. Calorimetric measurements indicated that the binding of these ligands to both the Y49F mutant and wild-type enzyme is enthalpically favorable and entropically unfavorable over the temperature range studied. The affinity of these ligands for the Y49F mutant is lower than those for the wild-type enzyme due mainly to an entropy change. Therefore, the thermodynamic effect of this mutation is to decrease the entropy loss due to binding. Calorimetric titrations in several buffers with different ionization heat amounts indicate a release of protons when the mutant binds GSH, whereas protons are taken up in binding S-hexylglutathione at pH 6.5. This suggests that the thiol group of GSH releases protons to buffer media during binding and a group with low pKa (such as Asp98) is responsible for the uptake of protons. The temperature dependence of the free energy of binding, DeltaG0, is weak because of the enthalpy-entropy compensation caused by a large heat capacity change. The heat capacity change is -199.5 +/- 26.9 cal K-1 mol-1 for GSH binding and -333.6 +/- 28.8 cal K-1 mol-1 for S-hexylglutathione binding. The thermodynamic parameters are consistent with the mutation Tyr49 --> Phe, producing a slight conformational change in the active site.

Distilling the essential features of a protein surface for improving protein-ligand docking, scoring, and virtual screening

For the successful identification and docking of new ligands to a protein target by virtual screening, the essential features of the protein and ligand surfaces must be captured and distilled in an efficient representation. Since the running time for docking increases exponentially with the number of points representing the protein and each ligand candidate, it is important to place these points where the best interactions can be made between the protein and the ligand. This definition of favorable points of interaction can also guide protein structure-based ligand design, which typically focuses on which chemical groups provide the most energetically favorable contacts. In this paper, we present an alternative method of protein template and ligand interaction point design that identifies the most favorable points for making hydrophobic and hydrogen-bond interactions by using a knowledge base. The knowledge-based protein and ligand representations have been incorporated in version 2.0 of SLIDE and resulted in dockings closer to the crystal structure orientations when screening a set of 57 known thrombin and glutathione S-transferase (GST) ligands against the apo structures of these proteins. There was also improved scoring enrichment of the dockings, meaning better differentiation between the chemically diverse known ligands and a approximately 15,000-molecule dataset of randomly-chosen small organic molecules. This approach for identifying the most important points of interaction between proteins and their ligands can equally well be used in other docking and design techniques. While much recent effort has focused on improving scoring functions for protein-ligand docking, our results indicate that improving the representation of the chemistry of proteins and their ligands is another avenue that can lead to significant improvements in the identification, docking, and scoring of ligands.

Design of a monomeric human glutathione transferase GSTP1, a structurally stable but catalytically inactive protein

By the introduction of 10 site-specific mutations in the dimer interface of human glutathione transferase P1-1 (GSTP1-1), a stable monomeric protein variant, GSTP1, was obtained. The monomer had lost the catalytic activity but retained the affinity for a number of electrophilic compounds normally serving as substrates for GSTP1-1. Fluorescence and circular dichroism spectra of the monomer and wild-type proteins were similar, indicating that there are no large structural differences between the subunits of the respective proteins. The GSTs have potential as targets for in vitro evolution and redesign with the aim of developing proteins with novel properties. To this end, a monomeric GST variant may have distinct advantages.

An ensemble of theta class glutathione transferases with novel catalytic properties generated by stochastic recombination of fragments of two mammalian enzymes

The correlation between sequence diversity and enzymatic function was studied in a library of Theta class glutathione transferases (GSTs) obtained by stochastic recombination of fragments of cDNA encoding human GST T1-1 and rat GST T2-2. In all, 94 randomly picked clones were characterized with respect to sequence, expression level, and catalytic activity in the conjugation reactions between glutathione and six alternative electrophilic substrates. Out of these six different compounds, dichloromethane is a selective substrate for human GST T1-1, whereas 1-menaphthyl sulfate and 1-chloro-2,4-dinitrobenzene are substrates for rat GST T2-2. The other three substances serve as substrates for both enzymes. Through this broad characterization, we have identified enzyme variants that have acquired novel activity profiles that differ substantially from those of the original GSTs. In addition, the expression levels of many clones were improved in comparison to the parental enzyme. A library of mutants can thus display a distribution of properties from which highly divergent evolutionary pathways may emerge, resembling natural evolutionary processes. From the GST library, a clone was identified that, by the point mutation N49D in the rat GST T2-2 sequence, has a 1700% increased activity with 1-menaphthyl sulfate and a 60% decreased activity with 4-nitrophenethyl bromide. Through the N49D mutation, the ratio of these activities has thus been altered 40-fold. An extensive characterization of a population of stochastically mutated enzymes can accordingly be used to find variants with novel substrate-activity profiles and altered catalytic properties. Recursive recombination of selected sequences displaying optimized properties is a strategy for the engineering of proteins for medical and biochemical applications. Such sequential design is combinatorial protein chemistry based on remodeling of existing structural scaffolds and has similarities to evolutionary processes in nature.

Solution structure of polyglutamine tracts in GST-polyglutamine fusion proteins

Aggregation of expanded polyglutamine (polyQ) seems to be the cause of various genetic neurodegenerative diseases. Relatively little is known as yet about the polyQ structure and the mechanism that induces aggregation. We have characterised the solution structure of polyQ in a proteic context using a model system based on glutathione S-transferase fusion proteins. A wide range of biophysical techniques was applied. For the first time, nuclear magnetic resonance was used to observe directly and selectively the conformation of polyQ in the pathological range. We demonstrate that, in solution, polyQs are in a random coil conformation. However, under destabilising conditions, their aggregation behaviour is determined by the polyQ length.

The crystal structures of glutathione S-transferases isozymes 1-3 and 1-4 from Anopheles dirus species B

Glutathione S-transferases (GSTs) are dimeric proteins that play an important role in cellular detoxification. Four GSTs from the mosquito Anopheles dirus species B (Ad), an important malaria vector in South East Asia, are produced by alternate splicing of a single transcription product and were previously shown to have detoxifying activity towards pesticides such as DDT. We have determined the crystal structures for two of these alternatively spliced proteins, AdGST1-3 (complexed with glutathione) and AdGST1-4 (apo form), at 1.75 and 2.45 A resolution, respectively. These GST isozymes show differences from the related GST from the Australian sheep blowfly Lucilia cuprina; in particular, the presence of a C-terminal helix forming part of the active site. This helix causes the active site of the Anopheles GSTs to be enclosed. The glutathione-binding helix alpha2 and flanking residues are disordered in the AdGST1-4 (apo) structure, yet ordered in the AdGST1-3 (GSH-bound) structure, suggesting that insect GSTs operate with an induced fit mechanism similar to that found in the plant phi- and human pi-class GSTs. Despite the high overall sequence identities, the active site residues of AdGST1-4 and AdGST1-3 have different conformations.

Molecular basis of mechanotransduction in living cells

The simplest cell-like structure, the lipid bilayer vesicle, can respond to mechanical deformation by elastic membrane dilation/thinning and curvature changes. When a protein is inserted in the lipid bilayer, an energetic cost may arise because of hydrophobic mismatch between the protein and bilayer. Localized changes in bilayer thickness and curvature may compensate for this mismatch. The peptides alamethicin and gramicidin and the bacterial membrane protein MscL form mechanically gated (MG) channels when inserted in lipid bilayers. Their mechanosensitivity may arise because channel opening is associated with a change in the protein's membrane-occupied area, its hydrophobic mismatch with the bilayer, excluded water volume, or a combination of these effects. As a consequence, bilayer dilation/thinning or changes in local membrane curvature may shift the equilibrium between channel conformations. Recent evidence indicates that MG channels in specific animal cell types (e.g., Xenopus oocytes) are also gated directly by bilayer tension. However, animal cells lack the rigid cell wall that protects bacteria and plants cells from excessive expansion of their bilayer. Instead, a cortical cytoskeleton (CSK) provides a structural framework that allows the animal cell to maintain a stable excess membrane area (i.e., for its volume occupied by a sphere) in the form of membrane folds, ruffles, and microvilli. This excess membrane provides an immediate membrane reserve that may protect the bilayer from sudden changes in bilayer tension. Contractile elements within the CSK may locally slacken or tighten bilayer tension to regulate mechanosensitivity, whereas membrane blebbing and tight seal patch formation, by using up membrane reserves, may increase membrane mechanosensitivity. In specific cases, extracellular and/or CSK proteins (i.e., tethers) may transmit mechanical forces to the process (e.g., hair cell MG channels, MS intracellular Ca(2+) release, and transmitter release) without increasing tension in the lipid bilayer.

Tyrosine 50 at the subunit interface of dimeric human glutathione transferase P1-1 is a structural key residue for modulating protein stability and catalytic function

The dimer interface in human GSTP1-1 has been altered by site-directed mutagenesis of Tyr50. It is shown that the effects of some mutations of this single amino acid residue are as detrimental for enzyme function as mutations of Tyr8 in the active site. The dimeric structure is a common feature of the soluble glutathione transferases and the structural lock-and-key motif contributing to the subunit-subunit interface is well conserved in classes alpha, mu, and pi. The key residue Tyr50 in GSTP1-1 was replaced with 5 different amino acids with divergent properties and the mutant proteins expressed and characterized. Mutant Y50F is an improved variant, with higher thermal stability and higher catalytic efficiency than the wild-type enzyme. The other mutants studied are also dimeric proteins, but have lower stabilities and catalytic activities that are reduced by a factor of 10(2)-10(4) from the wild-type value. Mutants Y50L and Y50T are characterized by a markedly increased K(M) value for GSH, while the effect is mainly due to decreased k(cat) values for mutants Y50A and Y50R. In conclusion, residue 50 in the interface governs both structural stability and catalytic function.

Molecular dynamics simulations of human glutathione transferase P1-1: conformational fluctuations of the apo-structure

We have investigated by molecular dynamics simulations the conformational fluctuations of the monomer of human apo-glutathione transferase P1-1. After attainment of steady-state dynamics, the structural fluctuations involve mainly the protein segments that participate also in the holo-apo transition discussed in the accompanying article (Stella et al., 1999:37:1-9.). The most mobile region is the C-terminal segment of helix 2. In contrast, helices 1, 6, 7, and 8 constitute a relatively rigid protein core. An "essential dynamics" analysis of the simulation shows that the largest fluctuations involve specific regions of glutathione transferases. In such regions, atomic motions are correlated. Motions of helix 2 are accounted for by the second most prominent principal component, which reveals a fluctuation between two distinct conformations. The residues that constitute the H-site undergo a breathing motion, possibly relevant during the binding of hydrophobic cosubstrates. Based on our simulation, several experimental findings can be rationalized, including the viscosity-dependent reactivity of Cys 47 and Cys 101 as well as the selective proteolysis of the peptide bond between Lys 44 and Ala 45. We have also modeled the structural changes that lead to the formation of an intrachain disulfide bridge between cysteines 47 and 101 and to the inactivation of the enzyme. The resulting structure maintains essentially the native fold except for helix 2, which closes the G-site. Proteins 1999;37:10-19.

The ligandin (non-substrate) binding site of human Pi class glutathione transferase is located in the electrophile binding site (H-site)

Glutathione S -transferases (GSTs) play a pivotal role in the detoxification of foreign chemicals and toxic metabolites. They were originally termed ligandins because of their ability to bind large molecules (molecular masses >400 Da), possibly for storage and transport roles. The location of the ligandin site in mammalian GSTs is still uncertain despite numerous studies in recent years. Here we show by X-ray crystallography that the ligandin binding site in human pi class GST P1-1 occupies part of one of the substrate binding sites. This work has been extended to the determination of a number of enzyme complex crystal structures which show that very large ligands are readily accommodated into this substrate binding site and in all, but one case, causes no significant movement of protein side-chains. Some of these molecules make use of a hitherto undescribed binding site located in a surface pocket of the enzyme. This site is conserved in most, but not all, classes of GSTs suggesting it may play an important functional role.

Temperature adaptation of glutathione S-transferase P1-1. A case for homotropic regulation of substrate binding

Human glutathione S-transferase P1-1 (GST P1-1) is a homodimeric enzyme expressed in several organs as well as in the upper layers of epidermis, playing a role against carcinogenic and toxic compounds. A sophisticated mechanism of temperature adaptation has been developed by this enzyme. In fact, above 35 degrees C, glutathione (GSH) binding to GST P1-1 displays positive cooperativity, whereas negative cooperativity occurs below 25 degrees C. This binding mechanism minimizes changes of GSH affinity for GST P1-1 because of temperature fluctuation. This is a likely advantage for epithelial skin cells, which are naturally exposed to temperature variation and, incidentally, to carcinogenic compounds, always needing efficient detoxifying systems. As a whole, GST P1-1 represents the first enzyme which displays a temperature-dependent homotropic regulation of substrate (e.g. GSH) binding.

Mutations of Gly to Ala in human glutathione transferase P1-1 affect helix 2 (G-site) and induce positive cooperativity in the binding of glutathione

Previous kinetic studies on human glutathione transferase P1-1 have indicated that the motions of an irregular alpha-helix (helix 2) lining the glutathione (GSH) binding site are viscosity dependent and may modulate the affinity of GSH binding. The effect of single amino acid residue substitutions (Gly to Ala) in this region is investigated here by site-directed mutagenesis. Three mutants (Gly41Ala, Gly50Ala and Gly41Ala/Gly50Ala) were overexpressed in Escherichia coli, purified, and characterized by kinetic, structural, and spectroscopic studies. All these mutant enzymes show kcat values similar to that of the wild-type enzyme, while the [S]0.5 for GSH increases about eight-fold in the Gly41Ala mutant and more than 100-fold in the Gly41Ala/Gly50Ala double mutant. This change in affinity towards GSH is accompanied by an induced positive cooperativity as reflected by Hill coefficients of 1.4 (Gly41Ala) and 1.7 (Gly41Ala/Gly50Ala) upon substrate binding. Taken together, these data suggest that the region around helix 2 is markedly altered leading to the observed intersubunit communication. Molecular modeling of the Gly41Ala/Gly50Ala mutant and of the inactive oxidized form of the native enzyme provides a structural explanation of our results.