1.3-A resolution structure of human glutathione S-transferase with S-hexyl glutathione bound reveals possible extended ligandin binding site

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.

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.

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.

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.

Myricetin as a Potential Adjuvant in Chemotherapy: Studies on the Inhibition of Human Glutathione Transferase A1–1

Glutathione transferases (GSTs) are a family of Phase II detoxification enzymes that are involved in the development of multi-drug resistance (MDR) phenomena toward chemotherapeutic agents. GST inhibitors are considered candidate compounds able to chemomodulate and reverse MDR. The natural flavonoid myricetin (MYR) has been shown to exhibit a wide range of pharmacological functions, including antitumor activity. In the present work, the interaction of MYR with human glutathione transferase A1–1 (hGSTA1–1) was investigated by kinetics inhibition analysis and molecular modeling studies. The results showed that MYR binds with high affinity to hGSTA1–1 (IC50 2.1 ± 0.2 μΜ). It functions as a non-competitive inhibitor towards the electrophile substrate 1-chloro−2,4-dinitrobenzene (CDNB) and as a competitive inhibitor towards glutathione (GSH). Chemical modification studies with the irreversible inhibitor phenethyl isothiocyanate (PEITC), in combination with in silico molecular docking studies allowed the prediction of the MYR binding site. MYR appears to bind at a distinct location, partially overlapping the GSH binding site (G-site). The results of the present study show that MYR is a potent inhibitor of hGSTA1–1 that can be further exploited towards the development of natural, safe, and effective GST-targeted cancer chemosensitizers.

An optimal acquisition scheme for Q-band EPR distance measurements using Cu2+-based protein labels

Recent advances in site-directed Cu²⁺ labeling of proteins and nucleic acids have added an attractive new methodology to measure the structure-function relationship in biomolecules. Despite the promise, accessing the higher sensitivity of Q-band Double Electron Electron Resonance (DEER) has been challenging for Cu²⁺ labels designed for proteins. Q-band DEER experiments on this label typically require many measurements at different magnetic fields, since the pulses can excite only a few orientations at a given magnetic field. Herein, we analyze such orientational effects through simulations and show that three DEER measurements, at strategically selected magnetic fields, are generally sufficient to acquire an orientational-averaged DEER time trace for this spin label at Q-band. The modeling results are experimentally verified on Cu²⁺ labeled human glutathione S-transferase (hGSTA1-1). The DEER distance distribution measured at the Q-band shows good agreement with the distance distribution sampled by molecular dynamics (MD) simulations and X-band experiments. The concordance of MD sampled distances and experimentally measured distances adds growing evidence that MD simulations can accurately predict distances for the Cu²⁺ labels, which remains a key bottleneck for the commonly used nitroxide label. In all, this minimal collection scheme reduces data collection time by as much as six-fold and is generally applicable to many octahedrally coordinated Cu²⁺ systems. Furthermore, the concepts presented here may be applied to other metals and pulsed EPR experiments.

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.

The Interaction of Human Glutathione Transferase GSTA1-1 with Reactive Dyes

Human glutathione transferase A1-1 (hGSTA1-1) contributes to developing resistance to anticancer drugs and, therefore, is promising in terms of drug-design targets for coping with this phenomenon. In the present study, the interaction of anthraquinone and diazo dichlorotriazine dyes (DCTD) with hGSTA1-1 was investigated. The anthraquinone dye Procion blue MX-R (PBMX-R) appeared to interact with higher affinity and was selected for further study. The enzyme was specifically and irreversibly inactivated by PBMX-R, following a biphasic pseudo-first-order saturation kinetics, with approximately 1 mol of inhibitor per mol of the dimeric enzyme being incorporated. Molecular modeling and protein chemistry data suggested that the modified residue is the Cys112, which is located at the entrance of the solvent channel at the subunits interface. The results suggest that negative cooperativity exists upon PBMX-R binding, indicating a structural communication between the two subunits. Kinetic inhibition analysis showed that the dye is a competitive inhibitor towards glutathione (GSH) and mixed-type inhibitor towards 1-chloro-2,4-dinitrobenzene (CDNB). The present study results suggest that PBMX-R is a useful probe suitable for assessing by kinetic means the drugability of the enzyme in future drug-design efforts.

Comparative analyses and structural insights of new class glutathione transferases in Cryptosporidium species

Cryptosporidiosis, caused by protozoan parasites of the genus Cryptosporidium, is estimated to rank as a leading cause in the global burden of neglected zoonotic parasitic diseases. This diarrheal disease is the second leading cause of death in children under 5 years of age. Based on the C. parvum transcriptome data, glutathione transferase (GST) has been suggested as a drug target against this pathogen. GSTs are diverse multifunctional proteins involved in cellular defense and detoxification in organisms and help pathogens to alleviate chemical and environmental stress. In this study, we performed genome-wide data mining, identification, classification and in silico structural analysis of GSTs in fifteen Cryptosporidium species. The study revealed the presence three GSTs in each of the Cryptosporidium species analyzed in the study. Based on the percentage identity and comprehensive comparative phylogenetic analysis, we assigned Cryptosporidium species GSTs to three new GST classes, named Vega (ϑ), Gamma (γ) and Psi (ψ). The study also revealed an atypical thioredoxin-like fold in the C. parvum GST1 of the Vega class, whereas C. parvum GST2 of the Gamma class and C. melagridis GST3 of the Psi class has a typical thioredoxin-like fold in the N-terminal region. This study reports the first comparative analysis of GSTs in Cryptosporidium species.

The mercapturic acid pathway

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

Sequence and Structure Properties Uncover the Natural Classification of Protein Complexes Formed by Intrinsically Disordered Proteins via Mutual Synergistic Folding

Intrinsically disordered proteins mediate crucial biological functions through their interactions with other proteins. Mutual synergistic folding (MSF) occurs when all interacting proteins are disordered, folding into a stable structure in the course of the complex formation. In these cases, the folding and binding processes occur in parallel, lending the resulting structures uniquely heterogeneous features. Currently there are no dedicated classification approaches that take into account the particular biological and biophysical properties of MSF complexes. Here, we present a scalable clustering-based classification scheme, built on redundancy-filtered features that describe the sequence and structure properties of the complexes and the role of the interaction, which is directly responsible for structure formation. Using this approach, we define six major types of MSF complexes, corresponding to biologically meaningful groups. Hence, the presented method also shows that differences in binding strength, subcellular localization, and regulation are encoded in the sequence and structural properties of proteins. While current protein structure classification methods can also handle complex structures, we show that the developed scheme is fundamentally different, and since it takes into account defining features of MSF complexes, it serves as a better representation of structures arising through this specific interaction mode.

ESR Resolves the C Terminus Structure of the Ligand-free Human Glutathione S-Transferase A1-1

Nitroxide- and Cu²⁺-based electron spin resonance (ESR) are combined to provide insight into the conformational states of the functionally important α-helix of the human glutathione S-transferase A1. Distance measurements on various spin-labeled dimeric human glutathione S-transferase A1-1 all result in bimodal distance distributions, indicating that the C-terminus exists in two distinct conformations in solution, one of which closely matches that found in the crystal structure of the ligand-bound enzyme. These measurements permit the generation of a model of the unliganded conformation. Room temperature ESR indicates that the second conformation has high mobility, potentially enabling the enzyme's high degree of substrate promiscuity. This model is then validated using computational modeling and further Cu²⁺-based ESR distance measurements. Cu²⁺-based ESR also provides evidence that the secondary structure of the second conformation is of helical nature. Addition of S-hexyl glutathione results in a shift in relative populations, favoring the state that is similar to the previously known structure of the ligand-bound enzyme.

Wucherria bancrofti glutathione S-Transferase: Insights into the 2.3 Å resolution X-ray structure and function, a therapeutic target for human lymphatic filariasis

The notoriety of parasitic nematode survival is directly related to chronic pathogenicity, which is evident in human lymphatic filariasis. It is a disease of poverty which causes severe disability affecting more than 120 million people worldwide. These nematodes down-regulate host immune system through a myriad of strategies that includes secretion of antioxidant and detoxification enzymes like glutathione-S-transferases (GSTs). Earlier studies have shown Wuchereria bancrofti GST to be a potential therapeutic target. Parasite GSTs catalyse the conjugation of glutathione to xenobiotic and other endogenous electrophiles and are essential for their long-term survival in lymph tissues. Hence, the crystal structure of WbGST along with its cofactor GSH at 2.3 Å resolution was determined. Structural comparisons against host GST reveal distinct differences in the substrate binding sites. The parasite xenobiotic binding site is more substrate/solvent accessible. The structure also suggests the presence of putative non-catalytic binding sites that may permit sequestration of endogenous and exogenous ligands. The structure of WbGST also provides a case for the role of the π-cation interaction in stabilizing catalytic Tyr compared to stabilization interactions described for other GSTs. Hence, the obtained information regarding crucial differences in the active sites will support future design of parasite specific inhibitors. Further, the study also evaluates the inhibition of WbGST and its variants by antifilarial diethylcarbamazine through kinetic assays.

The Cu2+-nitrilotriacetic acid complex improves loading of α-helical double histidine site for precise distance measurements by pulsed ESR

Site-directed spin labeling using two strategically placed natural histidine residues allows for the rigid attachment of paramagnetic Cu²⁺. This double histidine (dHis) motif enables extremely precise, narrow distance distributions resolved by Cu²⁺-based pulsed ESR. Furthermore, the distance measurements are easily relatable to the protein backbone-structure. The Cu²⁺ ion has, till now, been introduced as a complex with the chelating agent iminodiacetic acid (IDA) to prevent unspecific binding. Recently, this method was found to have two limiting concerns that include poor selectivity towards α-helices and incomplete Cu²⁺-IDA complexation. Herein, we introduce an alternative method of dHis-Cu²⁺ loading using the nitrilotriacetic acid (NTA)-Cu²⁺ complex. We find that the Cu²⁺-NTA complex shows a four-fold increase in selectivity toward α-helical dHis sites. Furthermore, we show that 100% Cu²⁺-NTA complexation is achievable, enabling precise dHis loading and resulting in no free Cu²⁺ in solution. We analyze the optimum dHis loading conditions using both continuous wave and pulsed ESR. We implement these findings to show increased sensitivity of the Double Electron-Electron Resonance (DEER) experiment in two different protein systems. The DEER signal is increased within the immunoglobulin binding domain of protein G (called GB1). We measure distances between a dHis site on an α-helix and dHis site either on a mid-strand or a non-hydrogen bonded edge-strand β-sheet. Finally, the DEER signal is increased twofold within two α-helix dHis sites in the enzymatic dimer glutathione S-transferase exemplifying the enhanced α-helical selectivity of Cu²⁺-NTA.

Dual Lifetimes for Complexes between Glutathione-S-transferase (hGSTA1-1) and Product-like Ligands Detected by Single-Molecule Fluorescence Imaging

Single-molecule fluorescence techniques were used to characterize the binding of products and inhibitors to human glutathione S-transferase A1-1 (hGSTA1-1). The identification of at least two different bound states for the wild-type enzyme suggests that there are at least two conformations of the protein, consistent with the model that ligand binding promotes closure of the carboxy-terminal helix over the active site. Ligand induced changes in ensemble fluorescence energy transfer support this proposed structural change. The more predominant state in the ensemble of single molecules shows a significantly faster off-rate, suggesting that the carboxy-terminal helix is delocalized in this state, permitting faster exit of the bound ligand. A point mutation (I219A), which is known to interfere with the association of the carboxy-terminal helix with the enzyme, shows increased rates of interconversion between the open and closed state. Kinematic traces of fluorescence from single molecules show that a single molecule readily samples a number of different conformations, each with a characteristic off-rate.

Energetics of Glutathione Binding to Human Eukaryotic Elongation Factor 1 Gamma: Isothermal Titration Calorimetry and Molecular Dynamics Studies

The energetics of ligand binding to human eukaryotic elongation factor 1 gamma (heEF1γ) was investigated using reduced glutathione (GSH), oxidised glutathione (GSSG), glutathione sulfonate and S-hexylglutathione as ligands. The experiments were conducted using isothermal titration calorimetry, and the findings were supported using computational studies. The data show that the binding of these ligands to heEF1γ is enthalpically favourable and entropically driven (except for the binding of GSSG). The full length heEF1γ binds GSSG with lower affinity (K _d = 115 μM), with more hydrogen-bond contacts (ΔH = -73.8 kJ/mol) and unfavourable entropy (-TΔS = 51.7 kJ/mol) compared to the glutathione transferase-like N-terminus domain of heEF1γ, which did not show preference to any specific ligand. Computational free binding energy calculations from the 10 ligand poses show that GSSG and GSH consistently bind heEF1γ, and that both ligands bind at the same site with a folded bioactive conformation. This study reveals the possibility that heEF1γ is a glutathione-binding protein.

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.

Comparison of epsilon- and delta-class glutathione S-transferases: the crystal structures of the glutathione S-transferases DmGSTE6 and DmGSTE7 from Drosophila melanogaster

Cytosolic glutathione transferases (GSTs) comprise a large family of enzymes with canonical structures that diverge functionally and structurally among mammals, invertebrates and plants. Whereas mammalian GSTs have been characterized extensively with regard to their structure and function, invertebrate GSTs remain relatively unstudied. The invertebrate GSTs do, however, represent potentially important drug targets for infectious diseases and agricultural applications. In addition, it is essential to fully understand the structure and function of invertebrate GSTs, which play important roles in basic biological processes. Invertebrates harbor delta- and epsilon-class GSTs, which are not found in other organisms. Drosophila melanogaster GSTs (DmGSTs) are likely to contribute to detoxication or antioxidative stress during development, but they have not been fully characterized. Here, the structures of two epsilon-class GSTs from Drosophila, DmGSTE6 and DmGSTE7, are reported at 2.1 and 1.5 Å resolution, respectively, and are compared with other GSTs to identify structural features that might correlate with their biological functions. The structures of DmGSTE6 and DmGSTE7 are remarkably similar; the structures do not reveal obvious sources of the minor functional differences that have been observed. The main structural difference between the epsilon- and delta-class GSTs is the longer helix (A8) at the C-termini of the epsilon-class enzymes.

Enzymatic detoxication, conformational selection, and the role of molten globule active sites

The role of conformational ensembles in enzymatic reactions remains unclear. Discussion concerning "induced fit" versus "conformational selection" has, however, ignored detoxication enzymes, which exhibit catalytic promiscuity. These enzymes dominate drug metabolism and determine drug-drug interactions. The detoxication enzyme glutathione transferase A1-1 (GSTA1-1), exploits a molten globule-like active site to achieve remarkable catalytic promiscuity wherein the substrate-free conformational ensemble is broad with barrierless transitions between states. A quantitative index of catalytic promiscuity is used to compare engineered variants of GSTA1-1 and the catalytic promiscuity correlates strongly with characteristics of the thermodynamic partition function, for the substrate-free enzymes. Access to chemically disparate transition states is encoded by the substrate-free conformational ensemble. Pre-steady state catalytic data confirm an extension of the conformational selection model, wherein different substrates select different starting conformations. The kinetic liability of the conformational breadth is minimized by a smooth landscape. We propose that "local" molten globule behavior optimizes detoxication enzymes.

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.

Cloning, expression and analysis of the olfactory glutathione S-transferases in coho salmon

The glutathione S-transferases (GSTs) provide cellular protection by detoxifying xenobiotics, maintaining redox status, and modulating secondary messengers, all of which are critical to maintaining olfaction in salmonids. Here, we characterized the major coho salmon olfactory GSTs (OlfGSTs), namely omega, pi, and rho subclasses. OlfGST omega contained an open reading frame of 720bp and encoded a protein of 239 amino acids. OlfGST pi and OlfGST rho contained open reading frames of 627 and 681nt, respectively, and encoded proteins of 208 and 226 amino acids. Whole-protein mass spectrometry yielded molecular weights of 29,950, 23,354, and 26,655Da, respectively, for the GST omega, pi, and rho subunits. Homology modeling using four protein-structure prediction algorithms suggest that the active sites in all three OlfGST isoforms resembled counterparts in other species. The olfactory GSTs conjugated prototypical GST substrates, but only OlfGST rho catalyzed the demethylation of the pesticide methyl parathion. OlfGST pi and rho exhibited thiol oxidoreductase activity toward 2-hydroxyethyl disulfide (2-HEDS) and conjugated 4-hydroxynonenal (HNE), a toxic aldehyde with neurodegenerative properties. The kinetic parameters for OlfGST pi conjugation of HNE were K(M)=0.16 ± 0.06mM and V(max)=0.5 ± 0.1μmolmin⁻¹mg⁻¹, whereas OlfGST rho was more efficient at catalyzing HNE conjugation (K(M)=0.022 ± 0.008 mM and V(max)=0.47 ± 0.05μmolmin⁻¹mg⁻¹). Our findings indicate that the peripheral olfactory system of coho expresses GST isoforms that detoxify certain electrophiles and pesticides and that help maintain redox status and signal transduction.

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.

A structural basis for cellular uptake of GST-fold proteins

It has recently emerged that glutathione transferase enzymes (GSTs) and other structurally related molecules can be translocated from the external medium into many different cell types. In this study we aim to explore in detail, the structural features that govern cell translocation and by dissecting the human GST enzyme GSTM2-2 we quantatively demonstrate that the α-helical C-terminal domain (GST-C) is responsible for this property. Attempts to further examine the constituent helices within GST-C resulted in a reduction in cell translocation efficiency, indicating that the intrinsic GST-C domain structure is necessary for maximal cell translocation capacity. In particular, it was noted that the α-6 helix of GST-C plays a stabilising role in the fold of this domain. By destabilising the conformation of GST-C, an increase in cell translocation efficiency of up to ∼2-fold was observed. The structural stability profiles of these protein constructs have been investigated by circular dichroism and differential scanning fluorimetry measurements and found to impact upon their cell translocation efficiency. These experiments suggest that the globular, helical domain in the ‘GST-fold’ structural motif plays a role in influencing cellular uptake, and that changes that affect the conformational stability of GST-C can significantly influence cell translocation efficiency.

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.

The role of a topologically conserved isoleucine in glutathione transferase structure, stability and function

The common fold shared by members of the glutathione-transferase (GST) family has a topologically conserved isoleucine residue at the N-terminus of helix 3 which is involved in the packing of helix 3 against two beta-strands in domain 1. The role of the isoleucine residue in the structure, function and stability of GST was investigated by replacing the Ile71 residue in human GSTA1-1 by alanine or valine. The X-ray structures of the I71A and I71V mutants resolved at 1.75 and 2.51 A, respectively, revealed that the mutations do not alter the overall structure of the protein compared with the wild type. Urea-induced equilibrium unfolding studies using circular dichroism and tryptophan fluorescence suggest that the mutation of Ile71 to alanine or valine reduces the stability of the protein. A functional assay with 1-chloro-2,4-dinitrobenzene shows that the mutation does not significantly alter the function of the protein relative to the wild type. Overall, the results suggest that conservation of the topologically conserved Ile71 maintains the structural stability of the protein but does not play a significant role in catalysis and substrate binding.

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.

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.

Methods for Predicting Human Drug Metabolism

Drug metabolism information is a necessary component of drug discovery and development. The key issues in drug metabolism include identifying: the enzyme(s) involved, the site(s) of metabolism, the resulting metabolite(s), and the rate of metabolism. Methods for predicting human drug metabolism from in vitro and computational methodologies and determining relationships between the structure and metabolic activity of molecules are also critically important for understanding potential drug interactions and toxicity. There are numerous experimental and computational approaches that have been developed in order to predict human metabolism which have their own limitations. It is apparent that few of the computational tools for metabolism prediction alone provide the major integrated functions needed to assist in drug discovery. Similarly the different in vitro methods for human drug metabolism themselves have implicit limitations. The utilization of these methods for pharmaceutical and other applications as well as their integration is discussed as it is likely that hybrid methods will provide the most success.

Differential gene expression profiling of vocal fold polyps and Reinke's edema by complementary DNA microarray

OBJECTIVES

Our purpose was to determine whether complementary DNA (cDNA) microarray analysis (MA) can establish distinct gene expression profiles for 2 phenotypically similar vocal fold lesions: Reinke's edema (RE) and polyps. Established transcript profiles can provide insight into the molecular and cellular processes involved in these diseases.

METHODS

Eleven RE specimens and 17 polyps were analyzed with MA for 8,745 genes. Further MA profiling was attempted within each lesion group to identify molecular markers for reflux exposure and smoking. Prediction analysis was used to predict lesion classification for 2 unclassified samples. A real-time polymerase chain reaction was performed to corroborate MA transcript levels for selected significant genes.

RESULTS

Sixty-five genes were found to differentiate RE and polyps (p = .0088). For RE, 19 genes were differentiated for reflux exposure (p = .016). No genes were found to differentiate smokers from nonsmokers. For polyps, no genes were found to differentiate for reflux (p = .16) and smoking (p = .565). Categorization of unclassified lesions was possible with a minimum of 13 genes.

CONCLUSIONS

We demonstrate the feasibility of benign lesion classification based on MA. Microarray analysis is useful not only for improving diagnosis and classification of such lesions, but also for potentially generating prognostic indicators and targets for therapy.

An intersubunit lock-and-key 'clasp' motif in the dimer interface of Delta class glutathione transferase

Structural investigations of a GST (glutathione transferase), adGSTD4-4, from the malaria vector Anopheles dirus show a novel lock-and-key 'Clasp' motif in the dimer interface of the Delta class enzyme. This motif also appears to be highly conserved across several insect GST classes, but differs from a previously reported mammalian lock-and-key motif. The aromatic 'key' residue not only inserts into a hydrophobic pocket, the 'lock', of the neighbouring subunit, but also acts as part of the 'lock' for the other subunit 'key'. The 'key' residues from both subunits show aromatic ring stacking with each other in a pi-pi interaction, generating a 'Clasp' in the middle of the subunit interface. Enzyme catalytic and structural characterizations revealed that single amino acid replacements in this 'Clasp' motif impacted on catalytic efficiencies, substrate selectivity and stability. Substitutions to the 'key' residue create strong positive co-operativity for glutathione binding, with a Hill coefficient approaching 2. The lock-and-key motif in general and especially the 'Clasp' motif with the pi-pi interaction appear to play a pivotal role in subunit communication between active sites, as well as in stabilizing the quaternary structure. Evidence of allosteric effects suggests an important role for this particular intersubunit architecture in regulating catalytic activity through conformational transitions of subunits. The observation of co-operativity in the mutants also implies that glutathione ligand binding and dimerization are linked. Quaternary structural changes of all mutants suggest that subunit assembly or dimerization basically manipulates subunit communication.

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.

Estradiol metabolites as isoform-specific inhibitors of human glutathione S-transferases

Numerous studies have suggested that the lifetime dose of unopposed estrogen is a significant risk factor for breast and uterine cancer. Estradiol (E2) plays a putative role as a tumor promoter through interaction with estrogen receptors but can also be metabolized to redox active and/or mutagenic semiquinones and quinones. Similarly, equine estrogens (components of certain hormone replacement therapy preparations) are converted to quinone metabolites. The use of hormone replacement therapy has also been associated with increased breast and endometrial cancer risk. Recently, metabolites of certain equine estrogens have been shown to inhibit human glutathione S-transferases (hGSTs). Since E2 and equine estrogens share similarities in other biological interactions, we have investigated the inhibitory capacity of endogenously formed E2 metabolites toward various hGSTs. The quinone metabolite of 2-hydroxy-17-beta-estradiol (2-OH-E2) was synthesized, and inhibition of hGST-mediated biotransformation of model substrates was assessed. Inhibition of purified recombinant hGSTM1-1 and hGSTA1-1 occurred in a concentration-dependent manner with IC50-values of approximately 250 and 350 nM, respectively. hGSTs M2-2, P1-1 and T1-1 were significantly less sensitive to inhibition. Specific glutathione-conjugates of the estrogen quinone also potently inhibited hGSTM1-1 and hGSTA1-1. Mass spectrometry data indicate that the inhibition was not mediated via covalent adduction. Although we have demonstrated hGST inhibition via E2 metabolites, our findings indicate that the isoform specificity and potency of GST inhibition by endogenous E2 metabolites is different than that of equine estrogen metabolites.

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

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

A conserved N-capping motif contributes significantly to the stabilization and dynamics of the C-terminal region of class Alpha glutathione S-transferases

Helix 9, the major structural element in the C-terminal region of class Alpha glutathione transferases, forms part of the active site of these enzymes where its dynamic properties modulate both catalytic and ligandin functions. A conserved aspartic acid N-capping motif for helix 9 was identified by sequence alignments of the C-terminal regions of class Alpha glutathione S-transferases (GSTs) and an analysis by the helix-coil algorithm AGADIR. The contribution of the N-capping motif to the stability and dynamics of the region was investigated by replacing the N-cap residue Asp-209 with a glycine in human glutathione S-transferase A1-1 (hGST A1-1) and in a peptide corresponding to its C-terminal region. Far-UV circular dichroism and AGADIR analyses indicate that, in the absence of tertiary interactions, the wild-type peptide displays a low intrinsic tendency to form a helix and that this tendency is reduced significantly by the Asp-to-Gly mutation. Disruption of the N-capping motif of helix 9 in hGST A1-1 alters the conformational dynamics of the C-terminal region and, consequently, the features of the H-site to which hydrophobic substrates (e.g. 1-chloro-2,4-dinitrobenzene (CDNB)) and nonsubstrates (e.g. 8-anilino-1-naphthalene sulfonate (ANS)) bind. Isothermal calorimetric and fluorescence data for complex formation between ANS and protein suggest that the D209G-induced perturbation in the C-terminal region prevents normal ligand-induced localization of the region at the active site, resulting in a less hydrophobic and more solvent-exposed H-site. Therefore, the catalytic efficiency of the enzyme with CDNB is diminished due to a lowered affinity for the electrophilic substrate and a lower stabilization of the transition state.

Characterization of the ligandin site of maize glutathione S-transferase I

Cytosolic GSTs (glutathione S-transferases) are a major reserve of high-capacity binding proteins and exhibit ligand-binding properties for a large variety of compounds. In the present study, the binding of two non-substrate anthraquinone dyes VBAR (Vilmafix Blue A-R) and CB3GA (Cibacron Blue 3GA) to maize (Zea mays) GST I was investigated. The results showed that the enzyme was specifically and irreversible inactivated by VBAR with a K(d) of 35.5+/-2.2 microM and a k(3) of 0.47 min(-1). Proteolytic cleavage of the VBAR-modified enzyme and subsequent separation of peptides gave only one modified peptide. Sequencing of the modified peptide revealed the target site of VBAR reaction to be Lys(41). CB3GA binds reversibly to GST I and behaves as a competitive inhibitor towards CDNB (1-chloro-2,4-dinitrobenzene) and glutathione. CB3GA binding to GST I is accompanied by a characteristic spectral change in the absorption at positive maximum (670 nm) which exhibited a hyperbolic dependence on dye concentration with a K(d) of 12.1+/-0.5 microM. Site-directed mutagenesis of selected residues (Trp(12), Phe(35), Lys(41), Asn(49), Gln(53), Ser(67) and Ile(118)) was employed, and the mutated enzymes were assessed for CB3GA binding. These results, together with molecular-modelling studies, established that the ligandin-binding site of GST I is located mainly in the hydrophobic binding site. The ability of VBAR to specifically inactivate GST I was exploited further to demonstrate the specific binding of several plant hormones and flavonoids to GST I. The inactivation of other GST isoenzymes by VBAR was also investigated, and it was concluded that VBAR may have wide applicability as an affinity label for probing structure-function relationships of GST isoenzymes.

Characterization of bromosulphophthalein binding to human glutathione S-transferase A1-1: thermodynamics and inhibition kinetics

In addition to their catalytic functions, GSTs (glutathione S-transferases) bind a wide variety of structurally diverse non-substrate ligands. This ligandin function is known to result in the inhibition of catalytic function. The interaction between hGSTA1-1 (human class Alpha GST with two type 1 subunits) and a non-substrate anionic ligand, BSP (bromosulphophthalein), was studied by isothermal titration calorimetry and inhibition kinetics. The binding isotherm is biphasic, best described by a set of two independent sites: a high-affinity site and a low-affinity site(s). The binding stoichiometries for these sites are 1 and 3 molecules of BSP respectively. BSP binds to the high-affinity site 80 times more tightly (K(d)=0.12 microM) than it does to the low-affinity site(s) (K(d)=9.1 microM). Binding at these sites is enthalpically and entropically favourable, with no linkage to protonation events. Temperature- and salt-dependent studies indicate the significance of hydrophobic interactions in the binding of BSP, and that the low-affinity site(s) displays low specificity towards the anion. Binding of BSP results in the release of ordered water molecules at these hydrophobic sites, which more than offsets unfavourable entropic changes during binding. BSP inhibition studies show that the binding of BSP to its high-affinity site does not inhibit hGSTA1-1. This site, located near Trp-20, may be related to the buffer-binding site observed in GSTP1-1. The low-affinity-binding site(s) for BSP is most probably located at or near the active site of hGSTA1-1. Binding to this site(s) results in non-competitive inhibition with respect to CDNB (1-chloro-2,4-dinitrobenzene) (K(i)(BSP)=16.8+/-1.9 microM). Given the properties of the H site and the relatively small size of the electrophilic substrate CDNB, it is plausible that the active site of the enzyme can simultaneously accommodate both BSP and CDNB. This would explain the non-competitive behaviour of certain inhibitors that bind the active site (e.g. BSP).

The anomalous pKa of Tyr-9 in glutathione S-transferase A1-1 catalyzes product release

The pKa of the catalytic Tyr-9 in glutathione S-transferase (GST) A1-1 is lowered from 10.3 to approximately 8.1 in the apoenzyme and approximately 9.0 with a GSH conjugate bound at the active site. However, a clear functional role for the unusual Tyr-9 pKa has not been elucidated. GSTA1-1 also includes a dynamic C terminus that undergoes a ligand-dependent disorder-to-order transition. Previous studies suggest a functional link between Tyr-9 ionization and C-terminal dynamics. Here we directly probe the role of Tyr-9 ionization in ligand binding and C-terminal conformation. An engineered mutant of rGSTA1-1, W21F/F222W, which contains a single Trp at the C terminus, was used as a fluorescent reporter of pH-dependent C-terminal dynamics. This mutant exhibited a pH-dependent change in Trp-222 emission properties consistent with changes in C-terminal solvation or conformation. The apparent pKa values for the conformational transition were 7.9 +/- 0.1 and 9.3 +/- 0.1 for the apoenzyme and ligand-bound enzyme, respectively, in excellent agreement with the pKa for Tyr-9 in these states. The Y9F/W21F/F222W mutant, however, exhibited no such pH-dependent changes. Time-resolved fluorescence anisotropy studies revealed a ligand-dependent, Tyr-9-dependent, change in the order parameter of Trp-222. However, no pH dependence was observed. In equilibrium and pre-steady-state ligand binding studies, product conjugate had a decreased equilibrium binding affinity (KD), concomitant with increased binding and dissociation rates, at higher pH values. Furthermore, the recovered pKa values for the pH-dependent microscopic rate constants ranged from 7.7 to 8.4, also in agreement with the pKa of Tyr-9. In contrast, the Y9F/W21F/F222W mutant had no pH-dependent transition in KD or rate constants for ligand binding or dissociation. The combined results indicate that the macroscopic populations of "open" and "closed" states of the C terminus are not determined solely by the ionization state of Tyr-9. However, the rates of transition between these states are faster for the ionized Tyr-9. The ionized Tyr-9 states provide a parallel pathway for product dissociation, which is kinetically and thermodynamically favored. In silico kinetic models further support the functional role for the parallel dissociation pathway provided by ionized Tyr-9.