Structural analysis of human alpha-class glutathione transferase A1-1 in the apo-form and in complexes with ethacrynic acid and its glutathione conjugate

Integrating Electron Paramagnetic Resonance Spectroscopy and Computational Modeling to Measure Protein Structure and Dynamics

Electron paramagnetic resonance (EPR) has become a powerful probe of conformational heterogeneity and dynamics of biomolecules. In this Review, we discuss different computational modeling techniques that enrich the interpretation of EPR measurements of dynamics or distance restraints. A variety of spin labels are surveyed to provide a background for the discussion of modeling tools. Molecular dynamics (MD) simulations of models containing spin labels provide dynamical properties of biomolecules and their labels. These simulations can be used to predict EPR spectra, sample stable conformations and sample rotameric preferences of label sidechains. For molecular motions longer than milliseconds, enhanced sampling strategies and de novo prediction software incorporating or validated by EPR measurements are able to efficiently refine or predict protein conformations, respectively. To sample large-amplitude conformational transition, a coarse-grained or an atomistic weighted ensemble (WE) strategy can be guided with EPR insights. Looking forward, we anticipate an integrative strategy for efficient sampling of alternate conformations by de novo predictions, followed by validations by systematic EPR measurements and MD simulations. Continuous pathways between alternate states can be further sampled by WE-MD including all intermediate states.

Changing human hair fibre colour and shape from the follicle

INTRODUCTION

Natural hair curvature and colour are genetically determined human traits, that we intentionally change by applying thermal and chemical treatments to the fibre. Presently, those cosmetic methodologies act externally and their recurrent use is quite detrimental to hair fibre quality and even to our health.

OBJECTIVES

This work represents a disruptive concept to modify natural hair colour and curvature. We aim to model the fibre phenotype as it is actively produced in the follicle through the topical delivery of specific bioactive molecules to the scalp.

METHODS

Transcriptome differences between curly and straight hairs were identified by microarray. In scalp samples, the most variable transcripts were mapped by in situ hybridization. Then, by using appropriate cellular models, we screened a chemical library of 1200 generic drugs, searching for molecules that could lead to changes in either fibre colour or curvature. A pilot-scale, single-centre, investigator-initiated, prospective, blind, bilateral (split-scalp) placebo-controlled clinical study with the intervention of cosmetics was conducted to obtain a proof of concept (RNEC n.92938).

RESULTS

We found 85 genes transcribed significantly different between curly and straight hair, not previously associated with this human trait. Next, we mapped some of the most variable genes to the inner root sheath of follicles, reinforcing the role of this cell layer in fibre shape moulding. From the drug library screening, we selected 3 and 4 hits as modulators of melanin synthesis and gene transcription, respectively, to be further tested in 33 volunteers. The intentional specific hair change occurred: 8 of 14 volunteers exhibited colour changes, and 16 of 19 volunteers presented curvature modifications, by the end of the study.

CONCLUSION

The promising results obtained are the first step towards future cosmetics, complementary or alternative to current methodologies, taking hair styling to a new level: changing hair from the inside out.

Direct observation of negative cooperativity in a detoxification enzyme at the atomic level by Electron Paramagnetic Resonance spectroscopy and simulation

The catalytic activity of human glutathione S-transferase A1-1 (hGSTA1-1), a homodimeric detoxification enzyme, is dependent on the conformational dynamics of a key C-terminal helix α9 in each monomer. However, the structural details of how the two monomers interact upon binding of substrates is not well understood and the structure of the ligand-free state of the hGSTA1-1 homodimer has not been resolved. Here, we used a combination of electron paramagnetic resonance (EPR) distance measurements and weighted ensemble (WE) simulations to characterize the conformational ensemble of the ligand-free state at the atomic level. EPR measurements reveal a broad distance distribution between a pair of Cu(II) labels in the ligand-free state that gradually shifts and narrows as a function of increasing ligand concentration. These shifts suggest changes in the relative positioning of the two α9 helices upon ligand binding. WE simulations generated unbiased pathways for the seconds-timescale transition between alternate states of the enzyme, leading to the generation of atomically detailed structures of the ligand-free state. Notably, the simulations provide direct observations of negative cooperativity between the monomers of hGSTA1-1, which involve the mutually exclusive docking of α9 in each monomer as a lid over the active site. We identify key interactions between residues that lead to this negative cooperativity. Negative cooperativity may be essential for interaction of hGSTA1-1 with a wide variety of toxic substrates and their subsequent neutralization. More broadly, this work demonstrates the power of integrating EPR distances with WE rare-events sampling strategy to gain mechanistic information on protein function at the atomic level.

Ethacrynic Acid: A Promising Candidate for Drug Repurposing as an Anticancer Agent

Ethacrynic acid (ECA) is a diuretic that inhibits Na-K-2Cl cotransporter (NKCC2) present in the thick ascending loop of Henle and muculo dens and is clinically used for the treatment of edema caused by excessive body fluid. However, its clinical use is limited due to its low bioavailability and side effects, such as liver damage and hearing loss at high doses. Despite this, ECA has recently emerged as a potential anticancer agent through the approach of drug repositioning, with a novel mechanism of action. ECA has been shown to regulate cancer hallmark processes such as proliferation, apoptosis, migration and invasion, angiogenesis, inflammation, energy metabolism, and the increase of inhibitory growth factors through various mechanisms. Additionally, ECA has been used as a scaffold for synthesizing a new material, and various derivatives have been synthesized. This review explores the potential of ECA and its derivatives as anticancer agents, both alone and in combination with adjuvants, by examining their effects on ten hallmarks of cancer and neuronal contribution to cancer. Furthermore, we investigated the trend of synthesis research of a series of ECA derivatives to improve the bioavailability of ECA. This review highlights the importance of ECA research and its potential to provide a cost-effective alternative to new drug discovery and development for cancer treatment.

A Key Role in Catalysis and Enzyme Thermostability of a Conserved Helix H5 Motif of Human Glutathione Transferase A1-1

Glutathione transferases (GSTs) are promiscuous enzymes whose main function is the detoxification of electrophilic compounds. These enzymes are characterized by structural modularity that underpins their exploitation as dynamic scaffolds for engineering enzyme variants, with customized catalytic and structural properties. In the present work, multiple sequence alignment of the alpha class GSTs allowed the identification of three conserved residues (E137, K141, and S142) at α-helix 5 (H5). A motif-directed redesign of the human glutathione transferase A1-1 (hGSTA1-1) was performed through site-directed mutagenesis at these sites, creating two single- and two double-point mutants (E137H, K141H, K141H/S142H, and E137H/K141H). The results showed that all the enzyme variants displayed enhanced catalytic activity compared to the wild-type enzyme hGSTA1-1, while the double mutant hGSTA1-K141H/S142H also showed improved thermal stability. X-ray crystallographic analysis revealed the molecular basis of the effects of double mutations on enzyme stability and catalysis. The biochemical and structural analysis presented here will contribute to a deeper understanding of the structure and function of alpha class GSTs.

Interactions Between Odorants and Glutathione Transferases in the Human Olfactory Cleft

Xenobiotic metabolizing enzymes and other proteins, including odorant-binding proteins located in the nasal epithelium and mucus, participate in a series of processes modulating the concentration of odorants in the environment of olfactory receptors (ORs) and finely impact odor perception. These enzymes and transporters are thought to participate in odorant degradation or transport. Odorant biotransformation results in 1) changes in the odorant quantity up to their clearance and the termination of signaling and 2) the formation of new odorant stimuli (metabolites). Enzymes, such as cytochrome P450 and glutathione transferases (GSTs), have been proposed to participate in odorant clearance in insects and mammals as odorant metabolizing enzymes. This study aims to explore the function of GSTs in human olfaction. Using immunohistochemical methods, GSTs were found to be localized in human tissues surrounding the olfactory epithelium. Then, the activity of 2 members of the GST family toward odorants was measured using heterologously expressed enzymes. The interactions/reactions with odorants were further characterized using a combination of enzymatic techniques. Furthermore, the structure of the complex between human GSTA1 and the glutathione conjugate of an odorant was determined by X-ray crystallography. Our results strongly suggest the role of human GSTs in the modulation of odorant availability to ORs in the peripheral olfactory process.

Discovery of a new non-substrate inhibitor of the 26.5 kDa glutathione transferase from Taenia solium by virtual screening

The inappropriate use of anthelmintics, such as praziquantel and albendazole, has generated resistance and the need to develop new drugs. Glutathione transferases, GSTs, are bisubstrate dimeric enzymes that constitute the main detoxification mechanism against electrophiles, drugs and oxidative damage in Taenia solium. Therefore, GSTs are important targets for the development of new anthelmintics. In this work, we reported a successful virtual screen aimed at the identification of novel inhibitors of a 26.5 kDa GST from T. solium (TsGST26). We found that a compound, i7, able to inhibit selectively TsGST26 concerning human GSTs, showing a non-competitive inhibition mechanism towards substrate glutathione with a K_i (GSH) of 55.7 μM and mixed inhibition towards the electrophilic substrate 1-chloro-2,4-dinitrobenzene with a K_i (CDNB) of 8.64 μM. These results are in agreement with those of docking simulations, which showed i7 binds a site adjacent to the electrophilic site and furthest from the glutathione site.

Structural and functional divergence of GDAP1 from the glutathione S-transferase superfamily

Mutations in ganglioside-induced differentiation-associated protein 1 (GDAP1) alter mitochondrial morphology and result in several subtypes of the inherited peripheral neuropathy Charcot-Marie-Tooth disease; however, the mechanism by which GDAP1 functions has remained elusive. GDAP1 contains primary sequence homology to the GST superfamily; however, the question of whether GDAP1 is an active GST has not been clearly resolved. Here, we present biochemical evidence, suggesting that GDAP1 has lost the ability to bind glutathione without a loss of substrate binding activity. We have revealed that the α-loop, located within the H-site motif is the primary determinant for substrate binding. Using structural data of GDAP1, we have found that critical residues and configurations in the G-site which canonically interact with glutathione are altered in GDAP1, rendering it incapable of binding glutathione. Last, we have found that the overexpression of GDAP1 in HeLa cells results in a mitochondrial phenotype which is distinct from oxidative stress-induced mitochondrial fragmentation. This phenotype is dependent on the presence of the transmembrane domain, as well as a unique hydrophobic domain that is not found in canonical GSTs. Together, we data point toward a non-enzymatic role for GDAP1, such as a sensor or receptor.

Mechanisms of promiscuity among drug metabolizing enzymes and drug transporters

Detoxication, or 'drug-metabolizing', enzymes and drug transporters exhibit remarkable substrate promiscuity and catalytic promiscuity. In contrast to substrate-specific enzymes that participate in defined metabolic pathways, individual detoxication enzymes must cope with substrates of vast structural diversity, including previously unencountered environmental toxins. Presumably, evolution selects for a balance of 'adequate' kcat /KM values for a wide range of substrates, rather than optimizing kcat /KM for any individual substrate. However, the structural, energetic, and metabolic properties that achieve this balance, and hence optimize detoxication, are not well understood. Two features of detoxication enzymes that are frequently cited as contributions to promiscuity include the exploitation of highly reactive versatile cofactors, or cosubstrates, and a high degree of flexibility within the protein structure. This review examines these intuitive mechanisms in detail and clarifies the contributions of the classic ligand binding models 'induced fit' (IF) and 'conformational selection' (CS) to substrate promiscuity. The available literature data for drug metabolizing enzymes and transporters suggest that IF is exploited by these promiscuous detoxication enzymes, as it is with substrate-specific enzymes, but the detoxication enzymes uniquely exploit 'IFs' to retain a wide range of substrates at their active sites. In contrast, whereas CS provides no catalytic advantage to substrate-specific enzymes, promiscuous enzymes may uniquely exploit it to recruit a wide range of substrates. The combination of CS and IF, for recruitment and retention of substrates, can potentially optimize the promiscuity of drug metabolizing enzymes and drug transporters.

Molecular characterization of kappa class glutathione S-transferase from the disk abalone (Haliotis discus discus) and changes in expression following immune and stress challenges

Glutathione S-transferase (GST; EC 2.5.1.18) isoenzymes represent a complex group of proteins that are involved in phase II detoxification in several organisms. In this study, GST kappa (GSTκ) from the disk abalone (Haliotis discus discus; AbGSTκ) was characterized at both the transcriptional and functional levels to determine its potential capacity to perform as a detoxification agent under conditions of different stress. The predicted AbGSTκ protein consists of 227 amino acids, with a predicted molecular weight of 25.6 kDa and a theoretical isoelectric point (pI) of 7.78. In silico analysis reveals that AbGSTκ is a disulfide bond formation protein A (DsbA), consisting of a thioredoxin domain, GSH binding sites (G-sites), and a catalytic residue. In contrast, no hydrophobic ligand binding site (H-site), or signal peptides, were detected. AbGSTκ showed the highest sequence identity with the orthologue from pufferfish (Takifugu obscurus) (60.0%). In a phylogenetic tree, AbGSTκ clustered closely together with other fish GSTκs, and was evolutionarily distanced from other cytosolic GSTs. The predicted three-dimensional structure clearly demonstrates that the dimer adopts a butterfly-like shape. A tissue distribution analysis revealed that GSTκ was highly expressed in the digestive tract, suggesting it has detoxification ability. Depending on the tissue and time, AbGSTκ showed different expression patterns, and levels of expression, following challenge of the abalone with immune stimulants. Enzyme kinetics of the purified recombinant proteins demonstrated its conjugating ability using 1-Chloro-2,4-dinitrobenzene (CDNB) and glutathione (GSH) as substrates, and suggested it has a low affinity for both substrates. The optimum temperature and pH for the rAbGSTκ GSH: CDNB conjugating activity were found to be 35 °C and pH 8, respectively indicating that the abalone is well adapted to a wide range of environmental conditions. Cibacron blue (100 μM) was capable of completely inhibiting rAbGSTκ (100%) with an IC₅₀ (half maximal inhibitory concentration) of 0.05 μM. A disk diffusion assay revealed that rAbGSTκ could significantly protect cells from H₂O₂, CdCl_2, and ZnCl₂. Altogether, this current study suggests that AbGSTκ is involved in detoxification and immunological host defense mechanisms and allows abalones to overcome stresses in order for them to have an increased chance of survival.

Activity-Based Probes for Isoenzyme- and Site-Specific Functional Characterization of Glutathione S-Transferases

Glutathione S-transferases (GSTs) comprise a diverse family of phase II drug metabolizing enzymes whose shared function is the conjugation of reduced glutathione (GSH) to endo- and xenobiotics. Although the conglomerate activity of these enzymes can be measured, the isoform-specific contribution to the metabolism of xenobiotics in complex biological samples has not been possible. We have developed two activity-based probes (ABPs) that characterize active GSTs in mammalian tissues. The GST active site is composed of a GSH binding "G site" and a substrate binding "H site". Therefore, we developed (1) a GSH-based photoaffinity probe (GSTABP-G) to target the "G site", and (2) an ABP designed to mimic a substrate molecule and have "H site" activity (GSTABP-H). The GSTABP-G features a photoreactive moiety for UV-induced covalent binding to GSTs and GSH-binding enzymes. The GSTABP-H is a derivative of a known mechanism-based GST inhibitor that binds within the active site and inhibits GST activity. Validation of probe targets and "G" and "H" site specificity was carried out using a series of competition experiments in the liver. Herein, we present robust tools for the characterization of enzyme- and active site-specific GST activity in mammalian model systems.

Identification of Potent Chloride Intracellular Channel Protein 1 Inhibitors from Traditional Chinese Medicine through Structure-Based Virtual Screening and Molecular Dynamics Analysis

Chloride intracellular channel 1 (CLIC1) is involved in the development of most aggressive human tumors, including gastric, colon, lung, liver, and glioblastoma cancers. It has become an attractive new therapeutic target for several types of cancer. In this work, we aim to identify natural products as potent CLIC1 inhibitors from Traditional Chinese Medicine (TCM) database using structure-based virtual screening and molecular dynamics (MD) simulation. First, structure-based docking was employed to screen the refined TCM database and the top 500 TCM compounds were obtained and reranked by X-Score. Then, 30 potent hits were achieved from the top 500 TCM compounds using cluster and ligand-protein interaction analysis. Finally, MD simulation was employed to validate the stability of interactions between each hit and CLIC1 protein from docking simulation, and Molecular Mechanics/Generalized Born Surface Area (MM-GBSA) analysis was used to refine the virtual hits. Six TCM compounds with top MM-GBSA scores and ideal-binding models were confirmed as the final hits. Our study provides information about the interaction between TCM compounds and CLIC1 protein, which may be helpful for further experimental investigations. In addition, the top 6 natural products structural scaffolds could serve as building blocks in designing drug-like molecules for CLIC1 inhibition.

Directed evolution of glutathione transferases towards a selective glutathione-binding site and improved oxidative stability

BACKGROUND

Glutathione transferases (GSTs) are a family of detoxification enzymes that catalyze the conjugation of glutathione (GSH) to electrophilic compounds.

METHODS

A library of alpha class GSTs was constructed by DNA shuffling using the DNA encoding the human glutathione transferase A1-1 (hGSTA1-1) and the rat glutathione transferase A1-1 (rGSTA1-1).

RESULTS

Activity screening of the library allowed the selection of a chimeric enzyme variant (GSTD4) that displayed high affinity towards GSH and GSH-Sepharose affinity adsorbent, higher k_cat/K_m and improved thermal stability, compared to the parent enzymes. The crystal structures of the GSTD4 enzyme in free form and in complex with GSH were determined to 1.6Å and 2.3Å resolution, respectively. Analysis of the GSTD4 structure showed subtle conformational changes in the GSH-binding site and in electron-sharing network that may contribute to the increased GSH affinity. The shuffled variant GSTD4 was further optimized for improved oxidative stability employing site-saturation mutagenesis. The Cys112Ser mutation confers optimal oxidative stability and kinetic properties in the GSTD4 enzyme.

CONCLUSIONS

DNA shuffling allowed the creation of a chimeric enzyme variant with improved properties, compared to the parent enzymes. X-ray crystallography shed light on how recombination of a specific segment from homologous GSTA1-1 together with point mutations gives rise to a new functionally competent enzyme with improved binding, catalytic properties and stability.

GENERAL SIGNIFICANCE

Such an engineered GST would be useful in biotechnology as affinity tool in affinity chromatography as well as a biocatalytic matrix for the construction of biochips or enzyme biosensors.

Members of the chloride intracellular ion channel protein family demonstrate glutaredoxin-like enzymatic activity

The Chloride Intracellular Ion Channel (CLIC) family consists of six evolutionarily conserved proteins in humans. Members of this family are unusual, existing as both monomeric soluble proteins and as integral membrane proteins where they function as chloride selective ion channels, however no function has previously been assigned to their soluble form. Structural studies have shown that in the soluble form, CLIC proteins adopt a glutathione S-transferase (GST) fold, however, they have an active site with a conserved glutaredoxin monothiol motif, similar to the omega class GSTs. We demonstrate that CLIC proteins have glutaredoxin-like glutathione-dependent oxidoreductase enzymatic activity. CLICs 1, 2 and 4 demonstrate typical glutaredoxin-like activity using 2-hydroxyethyl disulfide as a substrate. Mutagenesis experiments identify cysteine 24 as the catalytic cysteine residue in CLIC1, which is consistent with its structure. CLIC1 was shown to reduce sodium selenite and dehydroascorbate in a glutathione-dependent manner. Previous electrophysiological studies have shown that the drugs IAA-94 and A9C specifically block CLIC channel activity. These same compounds inhibit CLIC1 oxidoreductase activity. This work for the first time assigns a functional activity to the soluble form of the CLIC proteins. Our results demonstrate that the soluble form of the CLIC proteins has an enzymatic activity that is distinct from the channel activity of their integral membrane form. This CLIC enzymatic activity may be important for protecting the intracellular environment against oxidation. It is also likely that this enzymatic activity regulates the CLIC ion channel function.

Structure of Escherichia coli Grx2 in complex with glutathione: a dual-function hybrid of glutaredoxin and glutathione S-transferase

The structure of glutaredoxin 2 (Grx2) from Escherichia coli co-crystallized with glutathione (GSH) was solved at 1.60 Å resolution. The structure of a mutant with the active-site residues Cys9 and Cys12 changed to serine crystallized in the absence of glutathione was solved to 2.4 Å resolution. Grx2 has an N-terminal domain characteristic of glutaredoxins, and the overall structure is congruent with the structure of glutathione S-transferases (GSTs). Purified Grx2 exhibited GST activity. Grx2, which is the physiological electron donor for arsenate reduction by E. coli ArsC, was docked with ArsC. The docked structure could be fitted with GSH bridging the active sites of the two proteins. It is proposed that Grx2 is a novel Grx/GST hybrid that functions in two steps of the ArsC catalytic cycle: as a GST it catalyzes glutathionylation of the ArsC-As(V) intermediate and as a glutaredoxin it catalyzes deglutathionylation of the ArsC-As(III)-SG intermediate.

The interaction of the chemotherapeutic drug chlorambucil with human glutathione transferase A1-1: kinetic and structural analysis

Glutathione transferases (GSTs) are enzymes that contribute to cellular detoxification by catalysing the nucleophilic attack of glutathione (GSH) on the electrophilic centre of a number of xenobiotic compounds, including several chemotherapeutic drugs. In the present work we investigated the interaction of the chemotherapeutic drug chlorambucil (CBL) with human GSTA1-1 (hGSTA1-1) using kinetic analysis, protein crystallography and molecular dynamics. In the presence of GSH, CBL behaves as an efficient substrate for hGSTA1-1. The rate-limiting step of the catalytic reaction between CBL and GSH is viscosity-dependent and kinetic data suggest that product release is rate-limiting. The crystal structure of the hGSTA1-1/CBL-GSH complex was solved at 2.1 Å resolution by molecular replacement. CBL is bound at the H-site attached to the thiol group of GSH, is partially ordered and exposed to the solvent, making specific interactions with the enzyme. Molecular dynamics simulations based on the crystal structure indicated high mobility of the CBL moiety and stabilization of the C-terminal helix due to the presence of the adduct. In the absence of GSH, CBL is shown to be an alkylating irreversible inhibitor for hGSTA1-1. Inactivation of the enzyme by CBL followed a biphasic pseudo-first-order saturation kinetics with approximately 1 mol of CBL per mol of dimeric enzyme being incorporated. Structural analysis suggested that the modifying residue is Cys112 which is located at the entrance of the H-site. The results are indicative of a structural communication between the subunits on the basis of mutually exclusive modification of Cys112, indicating that the two enzyme active sites are presumably coordinated.

Differential scanning fluorometry signatures as indicators of enzyme inhibitor mode of action: case study of glutathione S-transferase

Differential scanning fluorometry (DSF), also referred to as fluorescence thermal shift, is emerging as a convenient method to evaluate the stabilizing effect of small molecules on proteins of interest. However, its use in the mechanism of action studies has received far less attention. Herein, the ability of DSF to report on inhibitor mode of action was evaluated using glutathione S-transferase (GST) as a model enzyme that utilizes two distinct substrates and is known to be subject to a range of inhibition modes. Detailed investigation of the propensity of small molecule inhibitors to protect GST from thermal denaturation revealed that compounds with different inhibition modes displayed distinct thermal shift signatures when tested in the presence or absence of the enzyme's native co-substrate glutathione (GSH). Glutathione-competitive inhibitors produced dose-dependent thermal shift trendlines that converged at high compound concentrations. Inhibitors acting via the formation of glutathione conjugates induced a very pronounced stabilizing effect toward the protein only when GSH was present. Lastly, compounds known to act as noncompetitive inhibitors exhibited parallel concentration-dependent trends. Similar effects were observed with human GST isozymes A1-1 and M1-1. The results illustrate the potential of DSF as a tool to differentiate diverse classes of inhibitors based on simple analysis of co-substrate dependency of protein stabilization.

Ensemble perspective for catalytic promiscuity: calorimetric analysis of the active site conformational landscape of a detoxification enzyme

Enzymological paradigms have shifted recently to acknowledge the biological importance of catalytic promiscuity. However, catalytic promiscuity is a poorly understood property, and no thermodynamic treatment has described the conformational landscape of promiscuous versus substrate-specific enzymes. Here, two structurally similar glutathione transferase (GST, glutathione S-transferase) isoforms with high specificity or high promiscuity are compared. Differential scanning calorimetry (DSC) indicates a reversible low temperature transition for the promiscuous GSTA1-1 that is not observed with substrate-specific GSTA4-4. This transition is assigned to rearrangement of the C terminus at the active site of GSTA1-1 based on the effects of ligands and mutations. Near-UV and far-UV circular dichroism indicate that this transition is due to repacking of tertiary contacts with the remainder of the subunit, rather than "unfolding" of the C terminus per se. Analysis of the DSC data using a modified Landau theory indicates that the local conformational landscape of the active site of GSTA1-1 is smooth, with barrierless transitions between states. The partition function of the C-terminal states is a broad unimodal distribution at all temperatures within this DSC transition. In contrast, the remainder of the GSTA1-1 subunit and the GSTA4-4 protein exhibit folded and unfolded macrostates with a significant energy barrier separating them. Their partition function includes a sharp unimodal distribution of states only at temperatures that yield either folded or unfolded macrostates. At intermediate temperatures the partition function includes a bimodal distribution. The barrierless rearrangement of the GSTA1-1 active site within a local smooth energy landscape suggests a thermodynamic basis for catalytic promiscuity.

Glutathione transferases: a structural perspective

The glutathione transferases (GSTs) are one of the most important families of detoxifying enzymes in nature. The classic activity of the GSTs is conjugation of compounds with electrophilic centers to the tripeptide glutathione (GSH), but many other activities are now associated with GSTs, including steroid and leukotriene biosynthesis, peroxide degradation, double-bond cis-trans isomerization, dehydroascorbate reduction, Michael addition, and noncatalytic "ligandin" activity (ligand binding and transport). Since the first GST structure was determined in 1991, there has been an explosion in structural data across GSTs of all three families: the cytosolic GSTs, the mitochondrial GSTs, and the membrane-associated proteins in eicosanoid and glutathione metabolism (MAPEG family). In this review, the major insights into GST structure and function will be discussed.

Interactions of glutathione transferases with 4-hydroxynonenal

Electrophilic products of lipid peroxidation are important contributors to the progression of several pathological states. The prototypical α,β-unsaturated aldehyde, 4-hydroxynonenal (HNE), triggers cellular events associated with oxidative stress, which can be curtailed by the glutathione-dependent elimination of HNE. The glutathione transferases (GSTs) are a major determinate of the intracellular concentration of HNE and can influence susceptibility to toxic effects, particularly when HNE and GST levels are altered in disease states. In this article, we provide a brief summary of the cellular effects of HNE, followed by a review of its GST-catalyzed detoxification, with an emphasis on the structural attributes that play an important role in the interactions with alpha-class GSTs. Some of the key determining characteristics that impart high alkenal activity reside in the unique C-terminal interactions of the GSTA4-4 enzyme. Studies encompassing both kinetic and structural analyses of related isoforms will be highlighted, with additional attention to stereochemical aspects that demonstrate the capacity of GSTA4-4 to detoxify both enantiomers of the biologically relevant racemic mixture while generating a select set of diastereomeric products with subsequent implications. A summary of the literature that examines the interplay between GSTs and HNE in model systems relevant to oxidative stress will also be discussed to demonstrate the magnitude of importance of GSTs in the overall detoxification scheme.

Identification of ligand binding sites of proteins using the Gaussian Network Model

The nonlocal nature of the protein-ligand binding problem is investigated via the Gaussian Network Model with which the residues lying along interaction pathways in a protein and the residues at the binding site are predicted. The predictions of the binding site residues are verified by using several benchmark systems where the topology of the unbound protein and the bound protein-ligand complex are known. Predictions are made on the unbound protein. Agreement of results with the bound complexes indicates that the information for binding resides in the unbound protein. Cliques that consist of three or more residues that are far apart along the primary structure but are in contact in the folded structure are shown to be important determinants of the binding problem. Comparison with known structures shows that the predictive capability of the method is significant.

Substrate specificity combined with stereopromiscuity in glutathione transferase A4-4-dependent metabolism of 4-hydroxynonenal

Conjugation to glutathione (GSH) by glutathione transferase A4-4 (GSTA4-4) is a major route of elimination for the lipid peroxidation product 4-hydroxynonenal (HNE), a toxic compound that contributes to numerous diseases. Both enantiomers of HNE are presumed to be toxic, and GSTA4-4 has negligible stereoselectivity toward them, despite its high catalytic chemospecificity for alkenals. In contrast to the highly flexible, and substrate promiscuous, GSTA1-1 isoform that has poor catalytic efficiency with HNE, GSTA4-4 has been postulated to be a rigid template that is preorganized for HNE metabolism. However, the combination of high substrate chemoselectivity and low substrate stereoselectivity is intriguing. The mechanism by which GSTA4-4 achieves this combination is important, because it must metabolize both enantiomers of HNE to efficiently detoxify the biologically formed mixture. The crystal structures of GSTA4-4 and an engineered variant of GSTA1-1 with high catalytic efficiency toward HNE, cocrystallized with a GSH-HNE conjugate analogue, demonstrate that GSTA4-4 undergoes no enantiospecific induced fit; instead, the active site residue Arg15 is ideally located to interact with the 4-hydroxyl group of either HNE enantiomer. The results reveal an evolutionary strategy for achieving biologically useful stereopromiscuity toward a toxic racemate, concomitant with high catalytic efficiency and substrate specificity toward an endogenously formed toxin.

Influence of the H-site residue 108 on human glutathione transferase P1-1 ligand binding: structure-thermodynamic relationships and thermal stability

The effect of the Y108V mutation of human glutathione S-transferase P1-1 (hGST P1-1) on the binding of the diuretic drug ethacrynic acid (EA) and its glutathione conjugate (EASG) was investigated by calorimetric, spectrofluorimetric, and crystallographic studies. The mutation Tyr 108 --> Val resulted in a 3D-structure very similar to the wild type (wt) enzyme, where both the hydrophobic ligand binding site (H-site) and glutathione binding site (G-site) are unchanged except for the mutation itself. However, due to a slight increase in the hydrophobicity of the H-site, as a consequence of the mutation, an increase in the entropy was observed. The Y108V mutation does not affect the affinity of EASG for the enzyme, which has a higher affinity (K(d) approximately 0.5 microM) when compared with those of the parent compounds, K(d) (EA) approximately 13 microM, K(d) (GSH) approximately 25 microM. The EA moiety of the conjugate binds in the H-site of Y108V mutant in a fashion completely different to those observed in the crystal structures of the EA or EASG wt complex structures. We further demonstrate that the Delta C(p) values of binding can also be correlated with the potential stacking interactions between ligand and residues located in the binding sites as predicted from crystal structures. Moreover, the mutation does not significantly affect the global stability of the enzyme. Our results demonstrate that calorimetric measurements maybe useful in determining the preference of binding (the binding mode) for a drug to a specific site of the enzyme, even in the absence of structural information.

Structural analysis of a glutathione transferase A1-1 mutant tailored for high catalytic efficiency with toxic alkenals

The specificity of human glutathione transferase (GST) A1-1 is drastically altered to favor alkenal substrates in the GIMFhelix mutant designed to mimic first-sphere interactions utilized by GSTA4-4. This redesign serves as a model for improving our understanding of the structural determinants that contribute to the distinct specificities of alpha class GSTs. Herein we report the first crystal structures of GIMFhelix, both in complex with GSH and in apo form at 1.98 and 2.38 A resolution. In contrast to the preorganized hydrophobic binding pocket that accommodates alkenals in GSTA4-4, GSTA1-1 includes a dynamic alpha9 helix that undergoes a ligand-dependent localization to complete the active site. Comparisons of the GIMFhelix structures with previously reported structures show a striking similarity with the GSTA4-4 active site obtained within an essentially GSTA1-1 scaffold and reveal the alpha9 helix assumes a similar localized structure regardless of active site occupancy in a manner resembling that of GSTA4-4. However, we cannot fully account for all the structural elements important in GSTA4-4 within the mutant's active site. The contribution of Phe10 to the Tyr212-Phe10-Phe220 network prevents complete C-terminal closure and demonstrates that the presence of Phe10 within the context of a GSTA4-4-like active site may ultimately hinder Phe220, a key C-terminal residue, from effectively contributing to the active site. In total, these results illustrate the remaining structural differences presumably reflected in the previously reported catalytic efficiencies of GIMFhelix and GSTA4-4 and emphasize the F10P mutation as being necessary to completely accomplish the transformation to a highly specific GST from the more promiscuous GSTA1-1 enzyme.

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.

Structure-Based Design of Anticancer Prodrug PABA/NO

Glutathione S-transferase (GST) is a superfamily of detoxification enzymes, represented by GSTα, GSTμ, GSTπ, etc. GSTα is the predominant isoform of GST in human liver, playing important roles for our well being. GSTπ is overexpressed in many forms of cancer, thus presenting an opportunity for selective targeting of cancer cells. Our structure-based design of prodrugs intended to release cytotoxic levels of nitric oxide in GSTπ-overexpressing cancer cells yielded PABA/NO, which exhibited anticancer activity both in vitro and in vivo with a potency similar to that of cisplatin. Here, we present the details on structural modification, molecular modeling, and enzymatic characterization for the design of PABA/NO. The design was efficient because it was on the basis of the reaction mechanism and the structures of related GST isozymes at both the ground state and the transition state. The ground-state structures outlined the shape and property of the substrate-binding site in different isozymes, and the structural information at the transition-state indicated distinct conformations of the Meisenheimer complex of prodrugs in the active site of different isozymes, providing guidance for the modifications of the molecular structure of the prodrug molecules. Two key alterations of a GSTα-selective compound led to the GSTπ-selective PABA/NO.

The stereochemical course of 4-hydroxy-2-nonenal metabolism by glutathione S-transferases

4-Hydroxy-2-nonenal (HNE) is a toxic aldehyde generated during lipid peroxidation and has been implicated in a variety of pathological states associated with oxidative stress. Glutathione S-transferase (GST) A4-4 is recognized as one of the predominant enzymes responsible for the metabolism of HNE. However, substrate and product stereoselectivity remain to be fully explored. The results from a product formation assay indicate that hGSTA4-4 exhibits a modest preference for the biotransformation of S-HNE in the presence of both enantiomers. Liquid chromatography mass spectrometry analyses using the racemic and enantioisomeric HNE substrates explicitly demonstrate that hGSTA4-4 conjugates glutathione to both HNE enantiomers in a completely stereoselective manner that is not maintained in the spontaneous reaction. Compared with other hGST isoforms, hGSTA4-4 shows the highest degree of stereoselectivity. NMR experiments in combination with simulated annealing structure determinations enabled the determination of stereochemical configurations for the GSHNE diastereomers and are consistent with an hGSTA4-4-catalyzed nucleophilic attack that produces only the S-configuration at the site of conjugation, regardless of substrate chirality. In total these results indicate that hGSTA4-4 exhibits an intriguing combination of low substrate stereoselectivity with strict product stereoselectivity. This behavior allows for the detoxification of both HNE enantiomers while generating only a select set of GSHNE diastereomers with potential stereochemical implications concerning their effects and fates in biological tissues.

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.

Functional Promiscuity Correlates with Conformational Heterogeneity in A-class Glutathione S-Transferases

The structurally related glutathione S-transferase isoforms GSTA1-1 and GSTA4-4 differ greatly in their relative catalytic promiscuity. GSTA1-1 is a highly promiscuous detoxification enzyme. In contrast, GSTA4-4 exhibits selectivity for congeners of the lipid peroxidation product 4-hydroxynonenal. The contribution of protein dynamics to promiscuity has not been studied. Therefore, hydrogen/deuterium exchange mass spectrometry (H/DX) and fluorescence lifetime distribution analysis were performed with glutathione S-transferases A1-1 and A4-4. Differences in local dynamics of the C-terminal helix were evident as expected on the basis of previous studies. However, H/DX demonstrated significantly greater solvent accessibility throughout most of the GSTA1-1 sequence compared with GSTA4-4. A Phe-111/Tyr-217 aromatic-aromatic interaction in A4-4, which is not present in A1-1, was hypothesized to increase core packing. "Swap" mutants that eliminate this interaction from A4-4 or incorporate it into A1-1 yield H/DX behavior that is intermediate between the wild type templates. In addition, the single Trp-21 residue of each isoform was exploited to probe the conformational heterogeneity at the intrasubunit domain-domain interface. Excited state fluorescence lifetime distribution analysis indicates that this core residue is more conformationally heterogeneous in GSTA1-1 than in GSTA4-4, and this correlates with greater stability toward urea denaturation for GSTA4-4. The fluorescence distribution and urea sensitivity of the mutant proteins were intermediate between the wild type templates. The results suggest that the differences in protein dynamics of these homologs are global. The results suggest also the possible importance of extensive conformational plasticity to achieve high levels of functional promiscuity, possibly at the cost of stability.

Bivalent inhibitors of glutathione S-transferase: the effect of spacer length on isozyme selectivity

Glutathione S-transferases (GSTs) are cytosolic enzymes that catalyze the conjugation of glutathione with a variety of exogenous and endogenous electrophiles. High affinity, isozyme-specific inhibitors of GST are required for use as pharmacological tools as well as potential therapeutics. The design of selective inhibitors is hindered due to the broad substrate binding capabilities of the GST enzymes. GSTs are dimeric enzymes, and therefore offer a unique discriminator for achieving inhibitor selectivity: the distance between binding sites on each monomer unit as a function of its quaternary organization. Bivalent analogs of the non-selective GST inhibitor ethacrynic acid were prepared, and selectivity for the GST A1-1 isozyme over GST P1-1 (IC50 values of 13.7 vs 1022 nM, respectively) was achieved through the optimization of the spacer length between the ethacrynic acid ligand domains.

Chemical modification at subunit 1 of rat kidney Alpha class glutathione transferase with 2,3,5,6-tetrachloro-1,4-benzoquinone: Close structural connectivity between glutathione conjugation activity and non-substrate ligand binding

2, 3, 5, 6-Tetrachloro-1, 4-benzoquinone (TCBQ) is a metabolite of pentachlorophenol known to react with cysteines of glutathione transferases (GSTs). TCBQ treatment of rat kidney rGSTA1-2 and rGSTA1-1 abolishes 70-80% conjugation of glutathione (GSH) to 1-chloro-2, 4-dinitrobenzene and results in strongly correlated quenching of intrinsic fluorescence of Trp-20 (R>0.96). rGSTA2-2 is only inhibited by 25%. Approximately 70% (rGSTA1-1) and 60% (rGSTA1-2) conjugation activity is abolished at TCBQ: GST stoichiometries near 1:1. The inactivation follows a Kitz/Wilson model with K(D) of 4.77+/-2.5microM for TCBQ and k(3) for inactivation of 0.036+/-0.01min(-1). A single tryptic peptide labelled with TCBQ was isolated from kidney rGSTA1-2 containing Cys-17 which we identify as the site of modification. Treatment with more than stoichiometric amounts of TCBQ modified other residues but resulted in only modest further inhibition of catalysis. We interpret these findings in terms of localised steric effects on the relatively rigid alpha-helix 1 adjacent to the catalytic site of subunit 1 possibly affecting the Alpha class-specific alpha-helix 9 which acts as a "lid" on the hydrophobic part of the active site. Homology modelling of rGSTA1-1 modified at Cys-17 of one subunit revealed only modest structural perturbations in the second subunit and tends to exclude global structural effects.

The intersubunit lock-and-key motif in human glutathione transferase A1-1: role of the key residues Met51 and Phe52 in function and dimer stability

The dimeric structure of certain cytosolic GSTs (glutathione S-transferases) is stabilized by a hydrophobic lock-and-key motif at their subunit interface. In hGSTA1-1 (human class Alpha GST with two type-1 subunits), the key consists of two residues, Met51 and Phe52, that fit into a hydrophobic cavity (lock) in the adjacent subunit. SEC (size-exclusion chromatography)-HPLC, far-UV CD and tryptophan fluorescence of the M51A and M51A/F52S mutants indicated the non-disruptive nature of these mutations on the global structure. While the M51A mutant retained 80% of wild-type activity, the activity of the M51A/F52S was markedly diminished, indicating the importance of Phe52 in maintaining the correct conformation at the active site. The M51A and M51A/F52S mutations altered the binding of ANS (8-anilinonaphthalene-l-sulphonic acid) at the H-site by destabilizing helix 9 in the C-terminal region. Data from urea unfolding studies show that the dimer is destabilized by both mutations and that the dimer dissociates to aggregation-prone monomers at low urea concentrations before global unfolding. Although not essential for the assembly of the dimeric structure of hGSTA1-1, both Met51 and Phe52 in the intersubunit lock-and-key motif play important structural roles in maintaining the catalytic and ligandin functions and stability of the GST dimer.

Structural basis of the suppressed catalytic activity of wild-type human glutathione transferase T1-1 compared to its W234R mutant

The crystal structures of wild-type human theta class glutathione-S-transferase (GST) T1-1 and its W234R mutant, where Trp234 was replaced by Arg, were solved both in the presence and absence of S-hexyl-glutathione. The W234R mutant was of interest due to its previously observed enhanced catalytic activity compared to the wild-type enzyme. GST T1-1 from rat and mouse naturally contain Arg in position 234, with correspondingly high catalytic efficiency. The overall structure of GST T1-1 is similar to that of GST T2-2, as expected from their 53% sequence identity at the protein level. Wild-type GST T1-1 has the side-chain of Trp234 occupying a significant portion of the active site. This bulky residue prevents efficient binding of both glutathione and hydrophobic substrates through steric hindrance. The wild-type GST T1-1 crystal structure, obtained from co-crystallization experiments with glutathione and its derivatives, showed no electron density for the glutathione ligand. However, the structure of GST T1-1 mutant W234R showed clear electron density for S-hexyl-glutathione after co-crystallization. In contrast to Trp234 in the wild-type structure, the side-chain of Arg234 in the mutant does not occupy any part of the substrate-binding site. Instead, Arg234 is pointing in a different direction and, in addition, interacts with the carboxylate group of glutathione. These findings explain our earlier observation that the W234R mutant has a markedly improved catalytic activity with most substrates tested to date compared to the wild-type enzyme. GST T1-1 catalyzes detoxication reactions as well as reactions that result in toxic products, and our findings therefore suggest that humans have gained an evolutionary advantage by a partially disabled active site.

Regio- and enantioselectivities in epoxide conjugations are modulated by residue 210 in Mu class glutathione transferases

The homologous human glutathione transferases (GSTs) M1-1 and M2-2 have similar catalytic activities with many electrophilic substrates, but differ strikingly in their conjugation of epoxides with glutathione. Residue 210, Thr in GST M2-2 and Ser in GST M1-1, is a key active-site component in determining the activity profile with epoxide substrates. This residue is hypervariable in Mu class GSTs, suggesting that it has special significance in the evolution of new functions. The present study shows that minor modifications of this residue can have major consequences for the enzyme-catalyzed epoxide conjugations. In general, a Ser at position 210 gives the highest catalytic efficiency, but the relatively high activity with an Ala placed on this position demonstrates that a hydroxyl group is not required. In contrast, a Thr residue suppresses the activity with epoxides by several orders of magnitude without major effects on the activity with alternative GST substrates. Residue 210 influences both the regio- and enantioselectivity with chiral and prochiral epoxides of stilbene and styrene and influences the distribution of isomeric glutathione conjugates. Thus, residue 210 contributes to both stereoselective recognition of the substrates and to partitioning of the isomeric reactants to the alternative transition states leading to separate chiral products.

Communication between the two active sites of glutathione S-transferase A1-1, probed using wild-type-mutant heterodimers

To study the communication between the two active sites of dimeric glutathione S-transferase A1-1, we used heterodimers containing one wild-type (WT) active site and one active site with a single mutation at either Tyr9, Arg15, or Arg131. Tyr9 and Arg15 are part of the active site of the same subunit, while Arg131 contributes to the active site of the opposite subunit. The V(max) values of Tyr9 and Arg15 mutant enzymes were less than 2% that of WT, indicating their importance in catalysis. In contrast, V(max) values of Arg131 mutant enzymes were about 50-90% of that of WT enzyme while K(m)(GSH) values were approximately 3-8 times that of WT, suggesting that Arg131 plays a role in glutathione binding. The mutant enzyme (with a His(6) tag) and the WT enzyme (without a His(6) tag) were used to construct heterodimers (WT-Y9F, WT-Y9T, WT-R15Q, WT-R131M, WT-R131Q, and WT-R131E) by incubation of a mixture of wild-type and mutant enzyme at pH 7.5 in buffer containing 1,6-hexanediol, followed by dialysis against buffer lacking the organic solvent. The resultant heterodimers were separated from the wild-type and mutant homodimers using chromatography on nickel-nitrilotriacetic acid agarose. The V(max) values of all heterodimers were lower than expected for independent active sites. Our experiments demonstrate that mutation of an amino acid residue in one active site affects the activity in the other active site. Modeling studies show that key amino acid residues and water molecules connect the two active sites. This connectivity is responsible for the cross-talk between the active sites.