Analysis of selective regions in the active sites of human cytochromes P450, 2C8, 2C9, 2C18, and 2C19 homology models using GRID/CPCA.

This study demonstrates a selectivity analysis using the GRID/CPCA strategy on four human cytochrome P450 2C homology models (CYP2C8, 2C9, 2C18, and 2C19). Although the four enzymes share more than 80% amino acid sequence identity, the substrate specificity differs. To investigate the selectivity of the enzymes and the amino acids that determine the specificity of each CYP2C enzyme, a selectivity analysis was made using GRID/CPCA. In the GRID calculations 10 probes were used covering hydrophobic, steric, and hydrogen bond acceptor and donor interactions. The selectivity analysis showed that the most important determinants of selectivity among the CYP2C models are the geometrical features of the active sites and the hydrophobic interactions. The selectivity analysis singled out CYP2C8 as the most different of the four CYP2C enzymes with amino acids with distinct properties in positions 114, 205, and 476 (Ser, Phe, and Ile, respectively) compared to the other enzymes. An inverse pharmacophore model for CYP2C9 was constructed from the selective regions, and the model agreed with the docking of diclofenac where the properties of the ligand overlapped with the pharmacophoric points in the model.

Competitive CYP2C9 inhibitors: enzyme inhibition studies, protein homology modeling, and three-dimensional quantitative structure-activity relationship analysis.

This study describes the generation of a three-dimensional quantitative structure activity relationship (3D-QSAR) model for 29 structurally diverse, competitive CYP2C9 inhibitors defined experimentally from an initial data set of 73 compounds. In parallel, a homology model for CYP2C9 using the rabbit CYP2C5 coordinates was built. For molecules with a known interaction mode with CYP2C9, this homology model, in combination with the docking program GOLD, was used to select conformers to use in the 3D-QSAR analysis. The remaining molecules were docked, and the GRID interaction energies for all conformers proposed by GOLD were calculated. This was followed by a principal component analysis (PCA) of the GRID energies for all conformers of all compounds. Based on the similarity in the PCA plot to the inhibitors with a known interaction mode, the conformer to be used in the 3D-QSAR analysis was selected. The compounds were randomly divided into two groups, the training data set (n = 21) to build the model and the external validation set (n = 8). The PLS (partial least-squares) analysis of the interaction energies against the K(i) values generated a model with r(2) = 0.947 and a cross-validation of q(2) = 0.730. The model was able to predict the entire external data set within 0.5 log units of the experimental K(i) values. The amino acids in the active site showed complementary features to the grid interaction energies in the 3D-QSAR model and were also in agreement with mutagenesis studies.

The active-site residue tyr-175 in human glyoxalase II contributes to binding of glutathione derivatives.

Tyrosine-175 located in the active site of human glyoxalase II was replaced by phenylalanine in order to study the contribution of this residue to catalysis. The mutation had a marginal effect on the k(cat) value determined using S-D-lactoylglutathione as substrate. However, the Y175F mutant had an 8-fold higher K(m) value than the wild-type enzyme. The competitive inhibitor S-(N-hydroxy-N-bromophenylcarbamoyl)glutathione had a 30-fold higher K(i) value towards the mutant, than that of the wild-type. Pre-equilibrium fluorescence studies with the inhibitor showed that this was due to a significantly increased off-rate for the mutant enzyme. The phenolic hydroxyl group of tyrosine-175 is within hydrogen bonding distance of the amide nitrogen of the glycine in the glutathione moiety and the present study shows that this interaction makes a significant contribution to the binding of the active-site ligand.

Arginines 97 and 108 in CYP2C9 are important determinants of the catalytic function.

Human cytochrome P450 2C9 (CYP2C9) is one of the major drug metabolising enzymes which exhibits a broad substrate specificity. The B-C loop is located in the active-site but has been difficult to model, owing to its diverse and flexible structure. To elucidate the function of the B-C loop we used homology modelling based on the Cyp102 structure in combination with functional studies of mutants using diclofenac as a model substrate for CYP2C9. The study shows the importance of the conserved arginine in position 97 and the arginine in position 108 for the catalytic function. The R97A mutant had a 13-fold higher K(m) value while the V(max) was in the same order as the wild type. The R108 mutant had a 100-fold lower activity with diclofenac compared to the wild-type enzyme. The other six mutants (S95A, F100A, L102A, E104A, R105A, and N107A) had kinetic parameters similar to the CYP2C9 wild-type. Our homology model based on the CYP102 structure as template indicates that R97, L102, and R105 are directed into the active site, whereas R108 is not. The change in catalytic function when arginine 97 was replaced with alanine and the orientation of this amino acid in our homology model indicates its importance for substrate interaction.

The human glutathione transferase P1-1 specific inhibitor TER 117 designed for overcoming cytostatic-drug resistance is also a strong inhibitor of glyoxalase I.

gamma-L-Glutamyl-S-(benzyl)-L-cysteinyl-R-(-)-phenylglycine (TER 117) has previously been developed for selective inhibition of human glutathione S-transferase P1-1 (GST P1-1) based on the postulated contribution of this isoenzyme to the development of drug resistance in cancer cells. In the present investigation, the inhibitory effect of TER 117 on the human glyoxalase system was studied. Although designed as an inhibitor specific for GST P1-1, TER 117 also competitively inhibits glyoxalase I (K(I) = 0.56 microM). In contrast, no inhibition of glyoxalase II was detected. Reduced glyoxalase activity is expected to raise intracellular levels of toxic 2-oxoaldehydes otherwise eliminated by glyoxalase I. The resulting toxicity would accompany the potentiation of cytostatic drugs, caused by inhibition of the detoxication effected by GST P1-1. TER 117 was designed for efficient inhibition of the most abundant form GST P1-1/Ile105. Therefore, the inhibitory effect of TER 117 on a second allelic variant GST P1-1/Val105 was also studied. TER 117 was shown to competitively inhibit both GST P1-1 variants. The apparent K(I) values at glutathione concentrations relevant to the intracellular milieu were in the micromolar range for both enzyme forms. Extrapolation to free enzyme produced K(I) values of approximately 0.1 microM for both isoenzymes, reflecting the high affinity of GST P1-1 for the inhibitor. Thus, the allelic variation in position 105 of GST P1-1 does not affect the inhibitory potency of TER 117. The inhibitory effects of TER 117 on GST P1-1 and glyoxalase I activities may act in synergy in the cell and improve the effectiveness of chemotherapy.

Reaction mechanism of glyoxalase I explored by an X-ray crystallographic analysis of the human enzyme in complex with a transition state analogue.

The structures of human glyoxalase I in complexes with S-(N-hydroxy-N-p-iodophenylcarbamoyl)glutathione (HIPC-GSH) and S-p-nitrobenzyloxycarbonylglutathione (NBC-GSH) have been determined at 2.0 and 1.72 A resolution, respectively. HIPC-GSH is a transition state analogue mimicking the enediolate intermediate that forms along the reaction pathway of glyoxalase I. In the structure, the hydroxycarbamoyl function is directly coordinated to the active site zinc ion. In contrast, the equivalent group in the NBC-GSH complex is approximately 6 A from the metal in a conformation that may resemble the product complex with S-D-lactoylglutathione. In this complex, two water molecules occupy the liganding positions at the zinc ion occupied by the hydroxycarbamoyl function in the enediolate analogue complex. Coordination of the transition state analogue to the metal enables a loop to close down over the active site, relative to its position in the product-like structure, allowing the glycine residue of the glutathione moiety to hydrogen bond with the protein. The structure of the complex with the enediolate analogue supports an "inner sphere mechanism" in which the GSH-methylglyoxal thiohemiacetal substrate is converted to product via a cis-enediolate intermediate. The zinc ion is envisioned to play an electrophilic role in catalysis by directly coordinating this intermediate. In addition, the carboxyl of Glu 172 is proposed to be displaced from the inner coordination sphere of the metal ion during substrate binding, thus allowing this group to facilitate proton transfer between the adjacent carbon atoms of the substrate. This proposal is supported by the observation that in the complex with the enediolate analogue the carboxyl group of Glu 172 is 3.3 A from the metal and is in an ideal position for reprotonation of the transition state intermediate. In contrast, Glu 172 is directly coordinated to the zinc ion in the complexes with S-benzylglutathione and with NBC-GSH.

Crystal structure of human glyoxalase II and its complex with a glutathione thiolester substrate analogue.

BACKGROUND: Glyoxalase II, the second of two enzymes in the glyoxalase system, is a thiolesterase that catalyses the hydrolysis of S-D-lactoylglutathione to form glutathione and D-lactic acid. RESULTS: The structure of human glyoxalase II was solved initially by single isomorphous replacement with anomalous scattering and refined at a resolution of 1.9 A. The enzyme consists of two domains. The first domain folds into a four-layered beta sandwich, similar to that seen in the metallo-beta-lactamases. The second domain is predominantly alpha-helical. The active site contains a binuclear zinc-binding site and a substrate-binding site extending over the domain interface. The model contains acetate and cacodylate in the active site. A second complex was derived from crystals soaked in a solution containing the slow substrate, S-(N-hydroxy-N-bromophenylcarbamoyl)glutathione. This complex was refined at a resolution of 1.45 A. It contains the added ligand in one molecule of the asymmetric unit and glutathione in the other. CONCLUSIONS: The arrangement of ligands around the zinc ions includes a water molecule, presumably in the form of a hydroxide ion, coordinated to both metal ions. This hydroxide ion is situated 2.9 A from the carbonyl carbon of the substrate in such a position that it could act as the nucleophile during catalysis. The reaction mechanism may also have implications for the action of metallo-beta-lactamases.

Unfolding and refolding of human glyoxalase II and its single-tryptophan mutants.

Here the structure of human glyoxalase II has been investigated by studying unfolding at equilibrium and refolding. Human glyoxalase II contains two tryptophan residues situated at the N-terminal (Trp57) and C-terminal (Trp199) regions of the molecule. Trp57 is a non-conserved residue located within a "zinc binding motif" (T/SHXHX57DH) which is strictly conserved in all known glyoxalase II sequences as well as in metal-dependent beta-lactamase and arylsulfatase. Site-directed mutagenesis has been used to construct single-tryptophan mutants in order to characterize better the guanidine-induced unfolding intermediates. The denaturation at equilibrium of wild-type glyoxalase II, as followed by activity, intrinsic fluorescence and CD, is multiphasic, suggesting that different regions of varying structural stability characterize the native structure of glyoxalase II. At intermediate denaturant concentration (1.2 M guanidine) a molten globule state is attained. The reactivation of the denatured wild-type enzyme occurs only in the presence of Zn(II) ions. The results show that Zn(II) is essential for the maintenance of the native structure of glyoxalase II and that its binding to the apoenzyme occurs during an essential step of refolding. The comparison of unfolding fluorescence transitions of single-trypthophan mutants with that of wild-type enzyme indicates that the strictly conserved "zinc binding motif" is located in a flexible region of the active site in which Zn(II) participates in catalysis.

Human glutathione transferase A4-4 crystal structures and mutagenesis reveal the basis of high catalytic efficiency with toxic lipid peroxidation products.

The oxidation of lipids and cell membranes generates cytotoxic compounds implicated in the etiology of aging, cancer, atherosclerosis, neurodegenerative diseases, and other illnesses. Glutathione transferase (GST) A4-4 is a key component in the defense against the products of this oxidative stress because, unlike other Alpha class GSTs, GST A4-4 shows high catalytic activity with lipid peroxidation products such as 4-hydroxynon-2-enal (HNE). The crystal structure of human apo GST A4-4 unexpectedly possesses an ordered C-terminal alpha-helix, despite the absence of any ligand. The structure of human GST A4-4 in complex with the inhibitor S-(2-iodobenzyl) glutathione reveals key features of the electrophilic substrate-binding pocket which confer specificity toward HNE. Three structural modules form the binding site for electrophilic substrates and thereby govern substrate selectivity: the beta1-alpha1 loop, the end of the alpha4 helix, and the C-terminal alpha9 helix. A few residue changes in GST A4-4 result in alpha9 taking over a predominant role in ligand specificity from the N-terminal loop region important for GST A1-1. Thus, the C-terminal helix alpha9 in GST A4-4 provides pre-existing ligand complementarity rather than acting as a flexible cap as observed in other GST structures. Hydrophobic residues in the alpha9 helix, differing from those in the closely related GST A1-1, delineate a hydrophobic specificity canyon for the binding of lipid peroxidation products. The role of residue Tyr212 as a key catalytic residue, suggested by the crystal structure of the inhibitor complex, is confirmed by mutagenesis results. Tyr212 is positioned to interact with the aldehyde group of the substrate and polarize it for reaction. Tyr212 also coopts part of the binding cleft ordinarily formed by the N-terminal substrate recognition region in the homologous enzyme GST A1-1 to reveal an evolutionary swapping of function between different recognition elements. A structural model of catalysis is presented based on these results.

Involvement of an active-site Zn2+ ligand in the catalytic mechanism of human glyoxalase I.

The Zn2+ ligands glutamate 99 and glutamate 172 in the active site of human glyoxalase I were replaced, each in turn, by glutamines by site-directed mutagenesis to elucidate their potential significance for the catalytic properties of the enzyme. To compensate for the loss of the charged amino acid residue, another of the metal ligands, glutamine 33, was simultaneously mutated into glutamate. The double mutants and the single mutants Q33E, E99Q, and E172Q were expressed in Escherichia coli, purified on an S-hexylglutathione matrix, and characterized. Metal analysis demonstrated that mutant Q33E/E172Q contained 1.0 mol of zinc/mol of enzyme subunit, whereas mutant Q33E/E99Q contained only 0.3 mol of zinc/mol of subunit. No catalytic activity could be detected with the double mutant Q33E/E172Q (<10(-8) of the wild-type activity). The second double mutant Q33E/E99Q had 1.5% of the specific activity of the wild-type enzyme, whereas the values for mutants Q33E and E99Q were 1.3 and 0. 1%, respectively; the E172Q mutant had less than 10(-5) times the specific activity of the wild-type. The crystal structure of the catalytically inactive double mutant Q33E/E172Q demonstrated that Zn2+ was bound without any gross changes or perturbations. The results suggest that the metal ligand glutamate 172 is directly involved in the catalytic mechanism of the enzyme, presumably serving as the base that abstracts a proton from the hemithioacetal substrate.

An evolutionary approach to the design of glutathione-linked enzymes.

Studies of protein structure provide information about principles of protein design that have come into play in natural evolution. This information can be exploited in the redesign of enzymes for novel functions. The glutathione-binding domain of glutathione transferases has similarities with structures in other glutathione-linked proteins, such as glutathione peroxidases and thioredoxin (glutaredoxin), suggesting divergent evolution from a common ancestral protein fold. In contrast, the binding site for glutathione in human glyoxalase I is located at the interface between the two identical subunits of the protein. Comparison with the homologous, but monomeric, yeast glyoxalase I suggests that new domains have originated through gene duplications, and that the oligomeric structure of the mammalian glyoxalase I has arisen by 'domain swapping'. Recombinant DNA techniques are being used for the redesign of glutathione-linked proteins in attempts to create binding proteins with novel functions and catalysts with tailored specificities. Enzymes with desired properties are selected from libraries of variant structures by use of phage display and functional assays.

Human glutathione transferase A4-4: an alpha class enzyme with high catalytic efficiency in the conjugation of 4-hydroxynonenal and other genotoxic products of lipid peroxidation.

A sequence encoding a novel glutathione transferase, GST A4-4, has been identified in a human fetal brain cDNA library. The protein has been produced in Escherichia coli after optimization of the codon usage for high-level heterologous expression. The dimeric protein has a subunit molecular mass of 25704 Da based on the deduced amino acid composition. Human GST A4-4 is a member of the Alpha class but shows only 53% amino acid sequence identity with the major liver enzyme GST A1-1. High catalytic efficiency with 4-hydroxyalkenals and other cytotoxic and mutagenic products of radical reactions and lipid peroxidation is a significant feature of GST A4-4. The kcat/Km values for 4-hydroxynonenal and 4-hydroxydecenal are > 3 x 10(6) M-1. s-1, several orders of magnitude higher than the values for conventional GST substrates. 4-Hydroxynonenal and other reactive electrophiles produced by oxidative metabolism have been linked to aging, atherosclerosis, cataract formation, Parkinson's disease and Alzheimer's disease, as well as other degenerative human conditions, suggesting that human GST A4-4 fulfills an important protective role and that variations in its expression may have significant pathophysiological consequences.

Crystal structure of human glyoxalase I--evidence for gene duplication and 3D domain swapping.

The zinc metalloenzyme glyoxalase I catalyses the glutathione-dependent inactivation of toxic methylglyoxal. The structure of the dimeric human enzyme in complex with S-benzyl-glutathione has been determined by multiple isomorphous replacement (MIR) and refined at 2.2 A resolution. Each monomer consists of two domains. Despite only low sequence homology between them, these domains are structurally equivalent and appear to have arisen by a gene duplication. On the other hand, there is no structural homology to the 'glutathione binding domain' found in other glutathione-linked proteins. 3D domain swapping of the N- and C-terminal domains has resulted in the active site being situated in the dimer interface, with the inhibitor and essential zinc ion interacting with side chains from both subunits. Two structurally equivalent residues from each domain contribute to a square pyramidal coordination of the zinc ion, rarely seen in zinc enzymes. Comparison of glyoxalase I with other known structures shows the enzyme to belong to a new structural family which includes the Fe2+-dependent dihydroxybiphenyl dioxygenase and the bleomycin resistance protein. This structural family appears to allow members to form with or without domain swapping.

Molecular cloning and characterization of the thiolesterase glyoxalase II from Arabidopsis thaliana.

cDNA encoding glyoxalase II from Arabidopsis thaliana has been cloned and sequenced. The isolated 894 bp segment included a sequence of 774 bp encoding a protein with a calculated molecular mass of 28,791 Da. The amino acid sequence deduced from the A. thaliana cDNA showed 54% identity with that of the human enzyme. Searches in databanks identified seven additional DNA sequences from different species with high similarity to glyoxalase II. Certain limited regions, one rich in histidine residues, shared 100% identity. A 29 kDa protein with an isoelectric point of 6.2 was obtained by heterologous expression of the A. thaliana cDNA in Escherichia coli. Homogeneous enzyme was obtained by affinity purification and its catalytic parameters with thiolesters of glutathione were similar to those for human glyoxalase II. The structural and functional similarities between glyoxalase II from A. thaliana and from human tissues suggest a common evolutionary origin.

Mutagenesis of residue 157 in the active site of human glyoxalase I.

Met-157 in the active site of human glyoxalase I was changed by site-directed mutagenesis into alanine, glutamine or histidine in order to evaluate its possible role in catalysis. The glyoxalase I mutants were expressed in Escherichia coli and purified on an S-hexylglutathione affinity gel. The physicochemical properties of the mutant proteins were similar to those of the wild-type enzyme. The glutamine mutant exhibited the same high specific activity as wild-type glyoxalase I, whereas the alanine and histidine mutants had approx. 20% of wild-type activity. The kcat/Km values of the mutant glyoxalase I determined with the hemithioacetal adduct of glutathione and methylglyoxal were reduced to between 10 and 40% of the wild-type value. This reduction was due to lower kcat values for the alanine and histidine mutants and a twofold increase in the Km value for the glutamine mutant. With the hemithioacetal of glutathione and phenylglyoxal, the kinetic parameters of the mutants were also of the same magnitude as those of wild-type glyoxalase I. Studies with the competitive inhibitors S-hexyl- and S-benzyl-glutathione revealed that the affinity was reduced to 7-11% of the wild-type affinity for the glutamine and alanine mutants and to 30-40% for the histidine mutant, as measured by a comparison of Ki values. The results show that Met-157 has no direct role in catalysis, but is rather involved in forming the substrate-binding site of human glyoxalase I. The high activity of the glutamine mutant suggests that a structurally equivalent glutamine residue in the N-terminal half of Saccharomyces cerevisiae glyoxalase I may be part of a catalytically competent active site.

Optimized heterologous expression of the human zinc enzyme glyoxalase I.

DNA coding for human glyoxalase I was isolated from a HeLa cell cDNA library by means of PCR. The deduced amino acid sequence differs form previously isolated sequences in that a glutamic acid replaces an alanine in position 111. This variant cDNA may represent the more acidic isoform of glyoxalase I originally identified at the protein level. An expression clone was constructed for high-level production of glyoxalase I in Escherichia coli. For optimal yield of the recombinant protein, silent random mutations were introduced in the cDNA coding region. Antisera against human glyoxalase I were used to select a high-level expression clone. This clone afforded 60 mg of purified enzyme per litre of culture medium. Addition of a zinc salt to the culture medium was essential to obtain an active enzyme and a stoicheiometric metal content. The functional characterization of the recombinant enzyme included determination of kinetic constants for methylglyoxal, phenylglyoxal and p-phenylphenylglyoxal, as well as inhibition studies. The kinetic properties of recombinant glyoxalase I were indistinguishable from those of the enzyme purified from human tissues.

Molecular cloning, heterologous expression, and characterization of human glyoxalase II.

A clone encoding glyoxalase II has been isolated from a human adult liver cDNA library. The sequence of 1011 base pairs consists of a full-length coding region of 780 base pairs, corresponding to a protein with a calculated molecular mass of 28,861 daltons. Identities (50-60%) were found to partial 5' and 3' cDNA sequences from Arabidopsis thaliana as well as within a limited region of glutathione transferase I cDNA from corn. A vector was constructed for heterologous expression of glyoxalase II in Escherichia coli. For optimal yield of enzyme, silent random mutations were introduced in the 5' coding region of the cDNA. A yield of 25 mg of glyoxalase II per liter of culture medium was obtained after affinity purification with immobilized glutathione. The recombinant enzyme had full catalytic activity and kinetic parameters indistinguishable from those of the native enzyme purified from human erythrocytes.

Similar toxic effect of 1,3-bis(2-chloroethyl)-1-nitrosourea on lymphocytes from human subjects differing in the expression of glutathione transferase M1-1.

Sixteen healthy donors were investigated for the presence or absence of glutathione transferase (GST) M1-1 in lymphocytes by immunodetection with polyclonal antibodies against human GST M1-1. Nine out of 16 individuals (56%) were categorized as GST M1-1 positive. Phytohaemagglutinin stimulated lymphocytes from GST M1-1 positive and negative donors were treated with 1,3-bis(2-chloroethyl)-1-nitrosourea (BCNU) and compared regarding inhibition of [3H]thymidine incorporation as a measure of cytotoxicity. No significant differences in the effect of BCNU were observed between the two groups, indicating that GST M1-1 is not an important resistance factor for BCNU.

Co-variation of glutathione transferase expression and cytostatic drug resistance in HeLa cells: establishment of class Mu glutathione transferase M3-3 as the dominating isoenzyme.

Qualitative and quantitative analyses of glutathione, glutathione transferases (GSTs) and other glutathione-linked enzymes in HeLa cells have been made in order to study their significance in cellular resistance to electrophilic cytotoxic agents. The cytosolic concentrations of three GSTs, GST M1-1 (53 +/- 9 ng/mg of cytosolic protein), GST P1-1 (11 +/- 3 ng/mg) and GST A1-1 (1.1 +/- 0.4 ng/mg) were quantified by isoenzyme-specific enzyme-linked immunoassays. Electrophoretic analysis and immunoblotting demonstrated another component, GST M3-3, which was identified by amino acid sequence analysis. GST M3-3 was quantified (1550 +/- 250 ng/mg) by slot-blot immunoanalysis and was the most abundant GST in HeLa cells. An additional cytosolic 13 kDa protein with high affinity for immobilized glutathione or S-hexyglutathione was found to be identical with a macrophage migration-inhibitory factor, previously identified as a lymphokine. Cells grown in roller bottles (HR) rather than in ordinary culture flasks contain a significantly lower concentration of all the GSTs and were found to be more sensitive to the cytostatic agents doxorubicin (2.3-fold), cisplatin (1.7-fold) and melphalan (1.4-fold). The cytosolic concentrations of glutathione reductase and glyoxalase I were also lower in HR cells, whereas the total glutathione concentration was unchanged and the glutathione peroxidase activity was increased. The results indicate that GSTs contribute to the cellular resistance phenotype.