The Anomalous pKa of Tyr-9 in Glutathione S-Transferase A1-1 Catalyzes Product Release*
Catherine A. Ibarra
,
Pramit Chowdhury
,
Jacob W. Petrich
and
William M. Atkins
¶
From the
Department of Medicinal Chemistry, University of Washington, Seattle, Washington 98195-7610,
Department of Chemistry, Iowa State University, Ames, Iowa 50011
Received for publication, February 13, 2003
, and in revised form, March 11, 2003.
 |
ABSTRACT
|
---|
The pKa of the catalytic Tyr-9 in glutathione S-transferase (GST) A1-1 is lowered from 10.3 to
8.1 in the apoenzyme and
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.
 |
INTRODUCTION
|
---|
Glutathione S-transferases (GSTs)1 comprise a superfamily of enzymes that catalyze the nucleophilic conjugation of tri-peptide (
-Glu-Cys-Gly) glutathione (GSH) with a wide variety of electrophilic compounds, including xenobiotics, products of oxidative stress, and endogenous toxins (1, 2). To date, several cytosolic GSTs in mammals have been identified on the basis of sequence similarity and substrate specificity: Alpha (A), Mu (M), Pi (P), Sigma (S), Theta (T), Zeta (Z), Kappa (K), and Omega (O) (3, 4, 5). In addition to their detoxification role, GSTs have been implicated in the development of drug resistance toward anticancer therapeutics (6, 7, 8).
A structural feature common to nearly all cytosolic GSTs is a conserved hydrogen bond between the hydroxyl group of catalytic Tyr or Ser and the thiol group of GSH. This interaction contributes to the lowering of the GSH pKa from
9.3 in solution to 6.57.4 at the active site of various GSTs. In the A-class GSTs, the catalytic Tyr has an unusually low pKa (9, 10). The spectroscopic pKa values of Tyr-9 in the human and rat A1-1 isoforms are 8.1 and 8.3, respectively, compared with 10.3 for Tyr in solution. In the GST A4-4 isoform, the pKa of the corresponding residue is even lower (11). Upon discovery of the anomalous pKa of Tyr-9, a general base mechanism was considered for the A-class GSTs, wherein Tyr abstracts a proton from GSH prior to attack of the resulting GS- at the electrophiles (9). However, solvent deuterium studies indicate that the catalytic Tyr in GSTM1-1, which has a "normal" pKa, does not undergo a general base catalysis (12). Subsequently, kinetic studies on hGSTA1-1 substituted with fluorotyrosine-9 analogs suggested that the pKa of Tyr-9 in the native enzyme is not sufficiently low to support general base catalysis (13). Thus, a functional role of the anomalous pKa of catalytic Tyr-9 in A-class GST has not been identified.
A structural feature unique to GSTA1-1 is the presence of a dynamic C-terminal helix (residues 208222) that covers the active site upon binding of ligand at the H-site or both sub-sites of the enzyme (Fig. 1). As described previously (14, 15, 16) for hGSTA1-1, wild-type rGSTA1-1, and several mutants, the binding of product GSH conjugate (GS-EA) follows a biphasic mechanism, which includes a bimolecular docking step (k1 and k-1) and a slower unimolecular isomerization step (k2 and k-2) that corresponds to the ordering and closure of the C terminus (Fig. 1). Indirect evidence from several laboratories suggests an intimate link between the ionization state of Tyr-9 and the dynamics of the C terminus (14, 15, 17, 18). One possible implication is that ligands will preferentially bind to or dissociate from the ionized state of the enzyme. Although substitution of Tyr-9 with Phe has demonstrated the importance of this residue in ligand binding, there is no direct evidence that has addressed whether the ionization state of Tyr-9 alters ligand binding kinetics. Similarly, no studies have probed the mechanism of "communication" between Tyr-9 and the C terminus. The objectives of the present studies were to determine directly whether the local conformational status of the C terminus is in fact controlled by the ionization properties of Tyr-9, and whether the extent and rate of ligand binding and release are dependent on the ionization state of Tyr-9. Our results indicate that the rates of the C-terminal disorder-to-order transition are faster when Tyr-9 is ionized, and equilibria for C-terminal opening and release of ligand are modestly shifted toward the dissociated state. Therefore, these results provide a clear, albeit modest, functional role for the unusual pKa of Tyr-9.

View larger version (31K):
[in this window]
[in a new window]
|
FIG. 1. Location of the C terminus in GSTA1-1 and ligand binding mechanism. The dimeric structure of [GSTA1-1·GS-EA] is shown with the C-terminal helix colored blue. The GS portion of the conjugate is green, and the ethacrynic acid portion is red. With product GSH conjugates bound, this helix is ordered and covers the active site. In the absence of ligand, this helix becomes crystallographically invisible or disordered. The binding of GSH conjugate GS-EA follows a biphasic mechanism, with the slow step corresponding to closure of the C terminus. The protein structure is from Ref. 10. See text for details.
|
|
 |
MATERIALS AND METHODS
|
---|
Chemicals and InstrumentationGSH, EA, MES, HEPES, CAPSO, and
were obtained from Sigma. GS-EA was synthesized and its purity analyzed as described previously (14). Stopped-flow experiments were performed with BioLogic SFM/QFM fluorimeter. Steady-state fluorescence was performed with Aminco SLM-8100 spectrofluorimeter.
Time-correlated Single Photon Counting Anisotropy DecayA coherent 701 rhodamine 6 G dye laser is pumped with about 1 watt of 532 nm radiation from an Antares 76-s CW mode-locked Nd:YAG laser. The 701 dye laser is cavity-dumped at 3.8 MHz. The pulses have an autocorrelation of about 10 ps frequency width at half-maximum. Excitation of tryptophan at 295 nm was achieved by focusing the dye laser pulses with a 5-cm lens onto a crystal of potassium diphosphate. Fluorescence is collected at right angles through a polarizer mounted at 54.7° to the excitation polarization and then passed through an ISA H-10 monochromator with a 160nm bandpass or through cut-off filters. A Hamamatsu 2809u microchannel plate, amplified by a Minicircuits ZHL-1042J, and a FFD 100 EG & G photodiode provide the start and stop signals, respectively. Constant-fraction discrimination of these signals is performed by a Tennelec TC 455 and time-to-amplitude conversion by an ORTEC 457. Data are stored in a Norland 5500 multichannel analyzer before transfer to and analysis with a PC. The instrument function of this system has a frequency width at half-maximum of 100150 ps. The polarized fluorescence traces used to obtain fluorescence anisotropy decay parameters are collected to a maximum of 16,000 counts in the peak channel. Samples were 5 µM protein, at the indicated pH, in 35 mM MES, 35 mM Tris, 35 mM CAPSO, with either no ligand or 50 µM S-hexyl-GSH.
Global Data FittingTime-resolved fluorescence data, acquired at various wavelengths, were subjected to a global fitting procedure found in Spectra SolveTM. A local
2 is calculated for each local fit. After all curves are fit (typically 611), a global
2 is calculated. The global parameters are varied, and the whole process is repeated to calculate a new global
2. This process is repeated until a minimum is reached for
2. The quality of the fit is determined by visual inspection and residuals.
Site-directed MutagenesisThe construction of W21F/F222W rGSTA1-1 has been described previously (14). The triple mutant Y9F/W21F/F222W was constructed by subcloning the available W21F/F222W insert, which was acquired from BglII and SalI digests, into a Y9F/W21F plasmid that had been linearized with a similar set of restriction enzymes. Originally, the single and double Y9F mutants were generated by PCR-based amplification of a fragment spanning the EcoRI and BglII restriction sites contained in the linearized pKKGTB34 or pKKGTB34-W21F plasmid (19). All mutants employed in the present studies were validated by DNA sequencing.
Protein Expression and PurificationThe expression and purification of wild-type rGSTA1-1 and its various mutants have been described previously (20). The purity of each protein was assessed by SDS-PAGE. Enzymatic activity of each mutant protein was assessed as described previously (21).
Determination of Conformational Transition in GSTA1-1Emission spectra were monitored at excitation wavelength of 295 nm, with slit widths of both monochromators set to 4 nm. In contrast to the previously employed tyrosine titration method (9), the emission properties of tryptophan were exploited to determine the nature of the surrounding environment of this residue. The pH dependence of Trp emission was monitored in 2 µM GST prepared in mixed buffer of 35 mM MES, 35 mM Tris, 35 mM CAPSO at 25 °C. The pKa values were determined by fitting the spectral center of mass, emission peak intensity, relative intensity at 305 nm, and relative intensity at 345 nm versus pH to an equation describing single residue ionization via Enzfitter. Fluorescence titrations were also performed in the presence of 200 µM S-hexylglutathione.
Equilibrium Binding StudiesThe equilibrium dissociation constants (KD) for the Michael addition product GS-EA were determined by the decrease in the fluorescence of W21F/F222W, Y9F/W21F/F222W, and
209222 upon ligand binding. During the titration experiments, serial dilutions of equilibrated [GST·GS-EA] complex, consisting initially of 0.5 µM GST and 100 µM GS-EA in mixed buffer system at pH 6.5 and 9, were made. The steady-state emission of each protein at 25 °C was monitored at an excitation wavelength of 295 nm, using slit widths of 4 nm. The observed fluorescence intensity at each dilution was corrected for the emission intensity of the reference buffer. The mole fraction of GS-EA bound to total protein was determined from a ratio of change in protein intensity when a ligand is bound and emission intensity of free enzyme. This fraction is further corrected for the dilution factor. KD values were then derived from Equation 1 via DeltaGraph 4,
 | (Eq. 1) |
where
is the mole fraction of GS-EA bound to protein, and P and C are the upper plateau and lower plateau, respectively, of the
versus [GS-EA]free plot.
Ligand Binding and Dissociation KineticsStopped-flow fluorescence studies were performed to determine the rates of ligand binding and dissociation as described in previous studies (14, 15, 16, 20). The rate constants for the binding reaction were measured by monitoring the decrease in protein fluorescence after rapid mixing of an equal volume of 2 µM GST and 6.25100 µM GS-EA at 15 °C and at various pH values. The multicomponent buffer used in these pH-dependent studies was 35 mM MES, 35 mM Tris, 35 mM CAPSO to ensure constant ionic strength. Kinetic data were analyzed and fitted to a single (Equation 2) or double (Equation 3) exponential decay equation,
 | (Eq. 2) |
 | (Eq. 3) |
where a1 and a2 are the amplitudes of the two observed rate constants kobs1 and kobs2; t is the time; and C is the offset. In cases where two rate constants were experimentally observed, the rapid association and dissociation rates, k1 and k-1(calc), were fitted to Equation 4 that describes the linear dependence of kobs1 on GS-EA concentration,
 | (Eq. 4) |
It should be noted that the k-1(calc) values derived from the binding experiment and the above equation are different from k-1(exp) values that were directly determined from the dissociation experiments. The association and dissociation rates of the slower step of the binding reaction, k2 and k-1(calc), were determined from a hyperbolic dependence (Equation 5) of kobs2 on GS-EA concentration,
 | (Eq. 5) |
Alternatively, a single observed rate, kobs, could be plotted against GS-EA concentration to yield a hyperbolic curve fit (Equation 6),
 | (Eq. 6) |
The dissociation of product from the [GST·GS-EA] complex was measured by monitoring the increase in fluorescence intensity after rapid mixing of equal volume of 2 µM complex and 4 mM of the trapping agent, GS-sulfonate (
). The dissociation rate constants, k-1(exp) and k-2(exp), were measured from a single (Equation 1) or double (Equation 2) exponential fit to the raw data.
Kinetic SimulationsSimulations were performed with Kinecyte, using the following expressions (Equation 7) for the time-dependent concentration of relevant species. The specific rate constants and initial concentrations used are indicated under "Results."
 | (Eq. 7) |
 | (Eq. 8) |
 | (Eq. 9) |
 | (Eq. 10) |
 |
RESULTS
|
---|
Characterization of C-terminal Dynamics by FluorescenceOur previous studies on the anomalous pKa value of catalytic Tyr-9 in rGSTA1-1 suggest that the ionization state of the tyrosyl hydroxyl group influences the dynamics of the C terminus (20, 18). In human and rat isoforms, the only tryptophan residue in the native GSTA1-1, Trp-21, near the domain-domain interface serves as a probe of global microconformational changes that occur upon ligand binding (22, 23). However, to investigate directly the localized structural changes at the C terminus, a Trp residue was engineered at the ultimate position of rGSTA1-1. The kinetics of ligand binding and dissociation for this mutant, W21F/F222W, are comparable with the true wild-type, indicating that global changes in protein conformation are coupled to local conformational changes at the C terminus (16). Thus, the W21F/F222W serves as a faithful mimic of ligand binding to the wild-type.
In order to determine the relationship between ionization state of Tyr-9 and C-terminal dynamics, here the W21F/F222W mutant was titrated between pH 7 and 10 to monitor the fluorescence emission properties of the terminal Trp. The emission of Trp, which absorbs at 275295 nm and emits at 340353 nm, is very sensitive to its local environment and is often used as a direct probe of conformational changes in proteins (23). The changes in the emission spectra of Trp-222 in the site-directed mutants reflect the polarity of the surrounding environment. Fig. 2 summarizes the typical emission spectra for the W21F/F222W and Y9F/W21F/F222W at several pH values. The emission intensity decreases, and the spectrum is red-shifted up to 7 nm with increasing pH in the W21F/F222W mutant. In contrast, the emission intensity of the Y9F/W21F/F222W mutant remained relatively constant, except at the very high pH values, at which the maximal spectral red shift was
2 nm. Clearly, W21F/F222W exhibits a large pH-dependent change in environment only when Tyr-9 is present. From the fluorescence titrations employed, the pKa values for this transition could be determined by fitting the spectral peak intensity, relative intensity at 305 nm, or relative intensity at 345 nm against pH via single ionization model. The average pKa values recovered from all curve fittings based on intensity change are 7.9 and 10.5 for W21F/F222W and Y9F/W21F/F222W, respectively (Table I).

View larger version (37K):
[in this window]
[in a new window]
|
FIG. 2. Emission of Trp-222 in the C terminus of rGSTA1-1 mutants. The fluorescence emission spectra of W21F/F222W (A) and Y9F/W21F/F222W (C) are shown at variable pH values. Although the emission of W21F/F222W was very sensitive to pH in the range of 78, the emission of Y9F/W21F/F222W was constant except at higher pH values. The pH dependence of emission peak intensity (B) and spectral center of mass (D) in the W21F/F222W (squares) and Y9F/W21F/F222W (circles) reflect the changes in polarity of the active site. Note that the additional ordinate axis on the right applies to the Y9F/W21F/F222W mutant. The recovered pKa values (Table I) are for the spectral transition based on peak intensity.
|
|
A more sensitive measurement of polarity at the active site of the enzyme is the pH dependence of spectral center of mass, which closely resembles the trend observed in the pH dependence of peak intensity (Fig. 2). However, the spectral center of mass in W21F/F222W does not reach a plateau after pH 9.5. In fact, the center of mass continues to decrease even after pH 10. Presumably, at high pH, tyrosine residues in the protein become ionized and contribute to the additional spectral red shift as observed in the control Y9F/W21F/F222W mutant. Whereas tyrosine or tyrosinate would not contribute significantly to the change in intensity at the emission maximum of Trp-222, their spectral signals would contribute to the overall spectral center of mass. Therefore, an accurate pKa value for the conformational transition based on Trp-222 parameters could not be obtained. Although the spectral center of mass data could not be fitted to an equation yielding a pKa, they clearly demonstrate a change in Trp polarity in the pH range of 7.58.5, which does not occur in the triple mutant. In addition, the overall pKa recovered from the intensity data of W21F/F222W most likely reflects the pKa of active site Tyr-9. Although the catalytic Tyr in the triple mutant had been removed, the apparent pKa of 10.5 is probably due to seven other tyrosines in the protein that become ionized at very high pH.
The binding of S-hexylglutathione (S-hexyl-GSH) to GSTA1-1 has been shown to induce closure of the C terminus, with Tyr-9 pKa (
9) remaining lower than the pKa of free tyrosine in solution (13, 24, 25). Here the emission spectra of W21F/F222W and Y9F/W21F/F222W were also obtained in the presence of saturating concentrations of S-hexyl-GSH. The S-hexyl-GSH was used because it provides no interference with the emission in fluorescence measurements. If the pH-dependent changes observed with W21F/F222W correlate with ionization of Tyr-9 in the apoenzyme, and if ionization is a major determinant of helix dynamics, then spectral changes would still be observed in the presence of S-hexyl-GSH but at a higher pH. The recovered pKa values for the conformational transition (Table I) indicate an increase of 1.4 pKa unit for the ligand-bound double mutant compared with the ligand-free double mutant, in excellent agreement with previously observed ligand-dependent increase in pKa of Tyr-9. The increase in Tyr-9 pKa of the conjugate-bound W21F/F222W suggests that, at equilibrium, only
3% of the complex remains unprotonated at pH 7.5. Also, the triple mutant Y9F/W21F/F222W was titrated with ligand at each pH, in order to determine whether the observed pH dependence was related to Tyr-9 ionization. As expected for a mutant with no ionizable group at residue 9, the bound conjugate had little effect on the apparent pKa for the conformational transition of the triple mutant (Table I).
Together, these spectroscopic studies directly demonstrate for the first time that the local conformational equilibria of the C terminus are pH-dependent in addition to being ligand-dependent. Moreover, this pH dependence exhibits an apparent pKa equivalent to the pKa of Tyr-9 (
8) in the apo form and pKa of
9 for the ligand-bound form, and this dependence is eliminated upon removal of the Tyr-9 phenolic hydroxyl.
Time-resolved AnisotropyThe observed pH dependence in steady-state spectral properties of Trp-222 could be due to a pH-dependent shift in the macroscopic equilibrium between closed and open states of the C terminus. Alternatively, there could be a pH-dependent change in the characteristics of the "open" conformation or the "closed," without a significant change in the relative amount of each. In order to probe these possibilities, time-resolved anisotropy experiments were performed. The anisotropy would be expected to increase upon opening the C terminus to fully expose the C-terminal Trp-222. The anisotropy decay of W21F/F222W was measured for the ligand-free protein at two limiting pH values, 6.5 and 9, in the presence and absence of S-hexyl-GSH, as probes of the open and closed states. Identical experiments were done with the triple mutant in order to assess the contribution of Tyr-9 to any observed effects. The recovered rotational correlation times, and the order parameters, S2, are summarized in Table II. The data were analyzed globally according to a biexponential decay with a fixed rotational correlation time for the overall protein motion (20 ns). This was based on the assumption that pH or ligand would not change the overall structure sufficiently to alter global protein hydrodynamics. The recovered rotational relaxation time for the whole protein, 20 ns, is consistent with a protein of 48 kDa. In addition, for each sample a fast anisotropy decay component was detected with correlation times ranging from 0.8 to 1.8 ns, clearly reflecting local motion of the C terminus. The results clearly demonstrate the expected increase in the order of Trp-222 upon binding of S-hexyl-GSH, at both pH values. This provides a benchmark for the magnitude of the change in S2 that is expected for closure of the C terminus. Furthermore, comparison of the double and triple mutants demonstrates that the Tyr-9 hydroxyl group is critical for maintaining the conformational integrity of the C terminus and the ligand-dependent switch. However, the moderate pH dependence observed for the ligand-free double mutant is observed also with the triple mutant. Based on these data, the ionization of Tyr-9 does not shift significantly the equilibrium between open and closed states, in contrast to the ligand, which does. Rather, the open state, the closed state, or both are different in the ionized versus the un-ionized forms, yielding different steady-state emission spectra. However, these differences do not correspond to changes in the fast local motion of the C-terminal residue expected for closed versus open conformations.
View this table:
[in this window]
[in a new window]
|
TABLE II
Anisotropy decay of Trp-222
1,r and 2,r are the rotational correlation times with fractional "zero time" limiting anisotropy values of r1 (0) and r2 (0), respectively. The anisotropy at time t is r(t) = r1(0)e(-t/r1) + r2(0)e(-t/r2), 1 and 2 are the rotational correlation times of the individual components, and S2 is the order parameter, defined as S2 = r2(0)/(r1(0) + r2(0)).
|
|
Equilibrium Binding AffinityEquilibrium binding studies were explored as an initial strategy to determine whether the pH-dependent microconformations of the C terminus differentially affected the binding of GSH conjugates. For several substrates, product release is rate-limiting for GSTA1-1 turnover, so the C-terminal conformational status is an important aspect of catalysis. These studies were performed at pH 6.5 and 9 at which the active site Tyr-9 in bound rGSTA1-1 (pKa
9.3) is completely protonated and significantly unprotonated (ionized), respectively. The ligand employed in the equilibrium studies was GS-EA, which is a conjugate formed from ethacrynic acid and GSH via Michael addition, and for which we have previously characterized the binding mechanism in detail. Moreover, crystal structures demonstrate that the C-terminal helix is closed and ordered with GS-EA bound. To prevent the proteins from unfolding and the GSH conjugate from hydrolyzing to its parent compounds at highly basic conditions, equilibrium binding studies beyond pH 9 were not pursued. Throughout our studies, the initial stock of GS-EA was kept on ice until use to avoid hydrolysis, which was monitored in a separate experiment via high performance liquid chromatography (data not shown).
The maximal emission intensities of W21F/F222W and Y9F/W21F/F222W at pH 6.5 were measured at 333 and 330 nm, respectively, at varying ratios of GS-EA:protein. In the presence of GS-EA at low pH, these emission maxima blue-shifted to 321 nm, as shown previously (25). At pH 9, only the W21F/F222W mutant exhibited red-shifted emission spectra. Equilibrium ligand binding studies were also performed in the C-terminal truncation mutant,
209222, for comparison. Typical binding isotherms, from which KD values were determined, are shown on Fig. 3. The recovered KD values are reported in Table III. At pH 6.5, the W21F/F222W mutant yielded a KD of 3.3 µM, which is comparable with a previously determined KD of wild-type protein (14). The KD for GS-EA significantly increased by 2.5-fold at pH 9, indicating lower affinity at equilibrium when Tyr-9 is ionized. In previous studies, the C-terminal truncation mutant,
209222, has been shown to have Tyr-9 pKa that is slightly lower than the wild-type (20). Moreover, with an open, solvent-accessible active site, GS-EA was found to bind to the
209222 mutant (KD
6.4 µM) with 3-fold increased KD compared with wild-type (14). In the present studies, the measured KD for GS-EA in this mutant is 8.4 µM at pH 6.5, which is 2.5-fold greater than W21F/F222W. Thus, physical truncation of the C terminus or increased pH lowers the ligand affinity to a similar extent. As expected for a mutant lacking the C terminus to restrict the access of solvent and ligand at the active site, the
209222 mutant did not show a pH dependence of KD (Table III). Similarly, the Y9F/W21F/F222W mutant showed no apparent pH dependence, as KD values at pH 6.5 and 9 remained essentially unchanged. The lack of pH dependence of KD for GS-EA, along with spectroscopic pKa of 10.5 for the conformational transition, suggests that ionization of Tyr-9 contributes to the ligand affinity through some type of communication with the C terminus. This is amplified by the observed pH dependence of KD in the W21F/F222W mutant. However, when access to the active site is no longer limited by the C terminus, such as in the
209222 mutant, facile ligand binding and release follow. As a result, KD is no longer dependent on pH in the absence of the C terminus.

View larger version (17K):
[in this window]
[in a new window]
|
FIG. 3. Equilibrium binding of GS-EA to rGSTA1-1 mutants at pH 6.5. The fluorescence of W21F/F222W (squares), Y9F/W21F/F222W (circles), and C-terminal truncation mutant (triangles) in complex with GS-EA was monitored at an excitation wavelength of 295 nm. The KD values for GS-EA (Table III) were determined from nonlinear regression fits to mole fraction of GS-EA bound to protein versus concentration of free GS-EA via Equation 1 ("Materials and Methods").
|
|
Kinetics of Ligand Binding and DissociationA further direct approach to examine the functional role of ionized Tyr-9 in the dynamics of the C terminus is a study of the pH dependence of individual rate constants that determine ligand binding. Stopped-flow fluorescence was employed to measure the individual rate constants that compose the binding mechanism of GSTA1-1 (Fig. 1). As described previously, the binding of GS-EA to rGSTA1-1 follows a two-step process, which includes a bimolecular docking step and a slow isomerization step that corresponds to the ordering and closure of the C terminus. The binding and dissociation of GS-EA were monitored by the decrease and increase in protein fluorescence, respectively, as we have demonstrated previously (20, 14, 15, 16). The rate constants for the association and dissociation of GS-EA were determined by fitting to appropriate equations as described under "Materials and Methods."
As we have shown previously, the raw data for wild-type fit best to double exponential equation (Equation 3), whereas the raw data for Y9F fit best to single exponential equation (Equation 2). Notably, the approach to equilibrium, described by k2 and k-2, is dramatically slower in the Y9F-containing mutant. The docking step could not be experimentally observed, presumably due to a fast k1. Fig. 4 depicts the dependence of observed rate constants (kobs) on GS-EA concentration in wild-type and Y9F at pH 6 and 9. Although only a single rate constant could be obtained from the raw data, the concentration dependence of kobs in the Y9F mutant yielded a hyperbolic curve, consistent with a two-step mechanism. Moreover, the rates of ligand binding were significantly increased at higher pH in the wild type. In contrast, ligand binding at the extreme pH values, 6 versus 9, is essentially identical in the Y9F mutant. As with the spectral changes for Trp-222 and the equilibrium ligand binding titration, the pH dependence observed in the double mutant is lost in the triple mutant. Apparently, Tyr-9 controls this pH dependence.

View larger version (13K):
[in this window]
[in a new window]
|
FIG. 4. Kinetics of GS-EA binding to wild-type rGSTA1-1 (A and B) and Y9F mutant (C). The stopped-flow fluorescence data for GS-EA binding in wild-type and Y9F fit best to double exponential Equation 3 and single exponential Equation 2, respectively. The rate constants kobs1 and kobs2 for wild-type showed linear dependence (Equation 4) and hyperbolic dependence (Equation 5), respectively, on GS-EA concentration (R2 = 0.99 and R2 = 0.92). Meanwhile, the single rate constant for Y9F mutant showed hyperbolic dependence (Equation 6) on GS-EA concentration (R2 = 0.98). Note that at pH 6 (squares) and pH 9 (circles), the rate constants for GS-EA binding to wild-type rGSTA1-1 showed apparent dependence on pH, although the rate constants for the Y9F mutant were essentially identical at the extreme pH values. Kinetic parameters are reported in Table IV.
|
|
Upon discovery of the pH dependence of GSH-conjugate binding and dissociation rates for the wild type, other stopped-flow experiments were undertaken to determine whether apparent pKa values could be extracted from the recovered rate constants as a function of pH. The rate constants for the docking and isomerization steps of the binding reaction (Table IV) showed a pH dependence and fit well to a single ionization model, with recovered pKa values ranging from 8.2 to 8.8 (Fig. 5). From this set of rate constants, calculated equilibrium dissociation constants (KD(calc)) could be obtained. Of course, the calculated KD values are pH-dependent, with an apparent pKa of 7.7 (Fig. 6). Apparently, all of the recovered pKa values (Table IV) reflect the ionization of Tyr-9, further suggesting that ligand docking and dissociation are favored in the open conformation of the enzyme. Due to difficulties in obtaining significant amounts of the Y9F proteins, which bind poorly to GSH affinity column, combined with the need for high protein concentrations in stopped-flow experiments, the kinetics of ligand binding and dissociation were determined at only three pH values. In contrast to the wild type, the Y9F mutant indicated no apparent pH dependence of rate constants (Table IV). The apparent calculated KD values for GS-EA in Y9F were not significantly different at any given pH (Fig. 6). Clearly, there is no indication that binding and dissociation of GS-EA in Y9F mutant depends on the ionization of a protein residue in this pH range.

View larger version (19K):
[in this window]
[in a new window]
|
FIG. 5. pH dependence of observed rate constants for GS-EA binding (A and C) and dissociation (B and D) in rGSTA1-1. Solid lines are nonlinear regression fits to a single ionization model (Enzfitter). The recovered pKa values from each fit (Table III) corresponded to the pKa of Tyr-9 in the ionized, open form of the enzyme.
|
|

View larger version (12K):
[in this window]
[in a new window]
|
FIG. 6. pH dependence of equilibrium ligand binding. The calculated KD for GS-EA with wild-type (A) as a function of pH resulted in pKa of 7.7, which closely reflects the ionization of Tyr-9 in the open conformational state. In contrast, Y9F (B) showed no apparent pH dependence of ligand dissociation. Solid lines represent fits to single ionization model.
|
|
It is important to note that the recovered microscopic rate constants at pH 9 are lower limits for the true rate constants when Tyr-9 is fully ionized. At pH 9, the Tyr-9 is
50% ionized, based on the experimentally measured pKa with GSH-conjugate bound. Presumably, if higher pH values were experimentally accessible, higher rates would have been observed, reflecting the true rates for Tyr-9 in the completely unprotonated form.
Kinetic SimulationsBecause the extent of ionization is low for the ligand-bound species, it is reasonable to consider whether the ionized state could have any effect on the rate of ligand binding and dissociation. What magnitude of rate change could be obtained with such a low extent of ionization? In order to relate the pH dependence of the experimentally determined microscopic rate constants to steady-state turnover of the enzyme, kinetic models were constructed, and simulations were performed. The kinetic scheme shown in Fig. 1 was used as a template but was expanded to include a parallel pathway of ligand binding and C-terminal ordering, via the state with unprotonated Tyr-9. The expanded model is shown in Fig. 7 and is oriented such that the reaction starts with the bound state, with the C terminus closed, [Tyr-9-OH·L]closed. The C terminus opens and closes, for the protonated state, with rate constants k-2 and k2, as above. For the unprotonated state the analogous rate constants are k-2' and k2'. Similarly, the ligand egress and association rate constants are k-1 and k1 for the protonated state and k-1' and k1' for the unprotonated state. The rate constants in Table IV for the wild-type enzyme at pH 6 were used directly to model the protonated pathway. The rate equations describing the concentration of each species are provided under "Materials and Methods." For the unprotonated pathway, each of the rate constants in Table IV for pH 9 was doubled, in order to mimic the expected increase in rate constant at full ionization; based on the pKa of 9 in the ligand-bound state, the experimental values are lower limits, as described above. The equilibria defined by each pKa value were modeled with a ratio of rate constants that yielded the appropriate ratio of protonated to unprotonated species, at equilibrium pH 7.4, i.e. 3% unprotonated for the closed, ligand-bound species ([Tyr-9-O-·L]closed; pKa = 9.0), 9% unprotonated for the open ligand bound species ([Tyr-9-O-·L]open; pKa = 8.5), and 24% unprotonated for the unbound species ([Tyr-9-O-]; pKa = 8.0). Starting with 1 µM ligand bound, closed, protonated species ([Tyr-9-OH·L]closed; upper left), in silico stopped-flow simulations were run to compare the time-dependent population of each species in the model and to simulate the experimental observations. The intent of these simulations is not to quantitatively reproduce the experimental results, but rather to qualitatively model the potential effect of Tyr-9 ionization on the rates of ligand dissociation. Unlike the actual experiments, the in silico stopped-flow experiment contains no trapping agent (
). Still, the time scale for approach to equilibrium in silico, with a half-time of
3 s, is in good agreement with the present and past experimental results (14). Fig. 7A shows the simulated results for decay of the starting species, [Tyr-9-OH·L]closed, for the simple model with no Tyr-9 ionization (dashed curve), and with the parallel pathway available through the unprotonated species (solid curve). As expected, the decay is faster for the later case. The half-life for the process is 2.4-fold shorter when the unprotonated pathway is available. The results in Fig. 7B demonstrate the time-dependent build up and decay of the protonated, open, ligand-bound intermediate [Tyr-9-OH·L]open. At its peak concentration, there is slightly less of this intermediate when the parallel unprotonated, pathway is available (solid curve) than when it is excluded (dashed curve). This is expected, because the parallel pathway competitively removes some of the available enzyme. Also, the rate of decay of this intermediate is faster when the parallel pathway is available. This is also expected as this intermediate has two pathways through which to decay, rather than one. Fig. 7C shows the build up of the ligand free open state for the protonated species. This behavior is interesting in as much as the rate of accumulation of ligand free [Tyr-9-OH] is faster at early times when the parallel pathway through the unprotonated states is available (solid curve). However, it is slower at later times and reaches a lower equilibrium concentration, when the parallel pathway is available. This behavior makes sense and validates the accuracy of the model. The equilibrium concentration of [Tyr-9-OH] is reduced by 24% because of population of the unprotonated species, [Tyr-9-O-]. When the simulation is run to 20 s, the difference in concentration is exactly 24% (not shown). The fact that the rate of appearance of unbound [Tyr-9-OH] is faster at early times reflects the faster release of ligand when the parallel pathway is available. In fact, Fig. 7D demonstrates this directly. Fig. 7D is the rate of appearance of free ligand. When the parallel pathway is available (solid curve), the free ligand accumulates faster. The half-life for approach to equilibrium is 2.4-fold faster (0.9 versus 2.1 s) when the unprotonated species contribute to the flux. These simulated results match reasonably well with the difference in experimental kobs at pH 6 versus pH 9, which differs by 4-fold. Based on these simulations, it is clear that even a minor extent of ionization of the initial ligand bound, closed, state can allow flux through the faster, unprotonated pathway with a modest increase in the rate of ligand dissociation.

View larger version (21K):
[in this window]
[in a new window]
|
FIG. 7. Parallel pathways for the ligand dissociation from ionized and un-ionized forms of Tyr-9. Top panels are the simulated time-dependent changes as follows: A, [Tyr-9-OH·L]closed; B, [Tyr-9-OH·L]open; C, ligand-free[Tyr-9-OH·L]free; and D, free ligand. The scheme below summarizes all possible states and rate constants. Although only 3% of the Tyr-9 is ionized in the ligand-bound state (left) at pH 7.5, the rates of C-terminal transitions (k2/k-2 and k2'/k-2') and product egress (k1/k-1 versus k1'/k-1) are faster for the ionized state. In particular, k-2' is at least 10-fold faster than k-2, and k-1' is at least 4.5-fold greater than k-1. In addition the equilibria defined by K2/k-2 versus k2'/k-', and by k1/k-1 versus k1'/k-' lie further toward the open and ligand-free states for the ionized form. See text for details. The rate constants use for simulation are as follows: k1 = 3.08 M- s-1; k-1 = 24.7 s-1; k2 = 16.4 s-1; k-2 = 0.63 s-1; k1' = 16.8 M-1 s-1; k-1' = 212 s-1; k2' = 108.6 s-1; k-2' = 13.36 s-1.
|
|
 |
DISCUSSION
|
---|
As already noted, there are many indirect data that indicate a functional role for Tyr-9 in ligand binding, which is mediated through the C terminus. For example, a thermodynamic coupling between Tyr-9 ionization and C-terminal dynamics is indicated by our previous pressure studies (9) in which we initially demonstrated a pressure-dependent ionization of Tyr-9, which was sensitive to C-terminal mutations. Also, the rates of ligand binding and dissociation are dramatically affected in the Y9F mutant compared with wild-type (15, 16, 17, see below), further implicating a role for Tyr-9 in C-terminal dynamics. However, there has been no direct evidence regarding the functional significance of coupling between C-terminal conformational status and Tyr-9 ionization.
Although various approaches have been used to address the pH dependence of ligand binding to GSTs, the focus in those studies had been on the ionization state of GSH and how its binding and release were subsequently affected (26, 27, 28, 29). In those cases, the rates of ligand binding and release were dependent on whether or not the GSH at the active site of several GST isoforms was unprotonated. Apparently, the rates of GSH binding and release were faster at lower pH, where the ligand is mostly protonated. Importantly, the available crystal structures indicate that the C terminus of GSTA1-1 does not close upon binding of GSH in the absence of an H-site ligand (30), so any pH dependence observed in those studies would not address C-terminal dynamics.
In contrast, the emphasis of our studies is on the ionization state of active site Tyr-9 and how it controls the rates of binding and release of GSH conjugates, for which crystallographic models clearly indicate ligand-induced closure of the C terminus. Because release of the product conjugate is rate-limiting for several substrates, these studies reveal new information about kinetically relevant steps of the reaction cycle. No previous studies have directly investigated the role of Tyr-9 pKa in the product dissociation mechanism.
To examine directly the relationship between the ionization state of Tyr-9 and the conformational status of the C terminus, here we have exploited four approaches as follows: 1) the pH dependence of conformational dynamics of the C terminus, as reported by an engineered Trp-222; 2) pH-dependent equilibrium binding studies; 3) stopped-flow studies in which binding and dissociation of ligand were directly measured at variable pH; and 4) kinetic simulations that mimic the dissociation of a GSH conjugate, based on experimentally derived parameters.
On the basis of emission properties of Trp-222, the local environments of the double and triple mutants were indeed different from each other. The decrease in the emission intensity and the red shift in the spectrum of W21F/F222W from pH 7 to 9.5 are indicative of polar, solvent-accessible environment at higher pH, or some type of conformational change. This pH-dependent emission profile was lost in the triple mutant containing Y9F. Moreover, the dependence of Trp-222 spectral intensity on pH yielded an apparent pKa of 7.9 in the W21F/F222W, which is remarkably close to the pKa for the ionization of Tyr-9 at the active site. In contrast, the pH dependence of spectral intensity in the Y9F/W21F/F222W resulted in a pKa of 10.5, which likely reflects the ionization of other tyrosines in the protein at the extreme pH values. In fact, no appreciable Tyr ionization of analogous mutant (Y9F:W21F) was observed at pH <11 (9). Presumably, the C terminus of Y9F/W21F/F222W mutant is significantly less solvated than W21F/F222W at pH values above 8.0, suggesting that the C terminus in the former likely adopts a different conformation. This conclusion is further supported by similar studies in the presence of ligand known to induce closure of the C terminus and alter the pKa of Tyr-9. The increase in Tyr-9 pKa of W21F/F222W, as a result of S-hexyl-GSH binding, is not only consistent with our previous studies but also lends support to our hypothesis that the ionization of Tyr-9 also occurs to a limited extent even with GSH conjugate bound. Indeed, the pH-dependent spectral properties of W21F/F222W in complex with S-hexyl-GSH are shifted to pKa 9.1, in agreement with the Tyr-9 pKa in this complex. Importantly, these are the first spectroscopic data that correlate local conformation of the C terminus with ionization of Tyr-9.
One possible scheme consistent with the steady-state spectral properties of the two mutants is that the C terminus is open and disordered when Tyr-9 is ionized, and it is closed and static when Tyr-9 is protonated, corresponding to the two macroconformations observed in crystal structures. However, the anisotropy decay results are not consistent with this scheme. If the open and closed states of the C terminus corresponded directly with the ionized and neutral forms of Tyr-9, then a much larger pH dependence of the fluorescence anisotropy order parameter would be expected. Obviously, the local dynamics of the ultimate residue in the C terminus are more sensitive to addition of ligand than to ionization of Tyr-9 (Table II). Whereas the anisotropy results nicely confirm the ligand-dependent disorder-to-order transition observed in the crystal and results from NMR, the results also indicate that a simple two-state model, with only open and closed states, is insufficient. The anisotropy data also confirm the importance of the Tyr-9 hydroxyl group, in as much as the ligand-dependent effects are lost upon substitution of Tyr-9 with Phe.
In order to accommodate all of the data, we propose that there is a conformational distribution of open states and a distribution of closed states. This is supported by the crystal structures with different ligands bound. Different GSH conjugates, or the inhibitor ethacrynic acid in the absence of GSH, yield readily definable electron density for the C terminus corresponding to the closed macroconformation. However, the exact location of the C terminus is different in each case, defining several possible microconformations within the closed state. Presumably, there are more microconformations sampled in the open state than in the closed, but these experiments do not address that issue. Together, all of the fluorescence data suggest that, whereas the ionization state of Tyr-9 is not the major determinant of the distribution between open and close states (based on anisotropy), it does perturb the local conformational distribution of each state (based on emission properties and kinetic behavior). Alternatively, it is formally possible that even when the majority of the helix is in the closed state, the ultimate Trp residue remains isotropic and indistinguishable from the open state.
Regardless of these details, the equilibrium and pre-steady-state binding experiments indicate a modest, but clear, functional advantage for the ionized form of Tyr-9. To the extent that product release is rate-limiting for many substrates, retention of product conjugates is a kinetic "problem" for GSTA1-1. Furthermore, product inhibition is commonly observed as a result of high affinity of GSH conjugates for the enzyme. The equilibrium binding experiments reveal a modest, 34-fold, decrease in affinity of ligand as the Tyr-9 becomes unprotonated. Also, all of the individual rate constants describing the biphasic binding and dissociation reactions are increased upon deprotonation of Tyr-9. In particular, the microscopic rate constant k-2, which corresponds to the first step of rate-limiting product dissociation, and which we propose corresponds to opening of the helix, is increased by 10-fold upon changing the pH from 6 to 9. Furthermore, although it is experimentally challenging to determine a residue pKa by stopped-flow, the variable pH experiments yield pKa values remarkably close to that of Tyr-9, and this pH dependence is lost in the triple mutant.
It is instructive to compare these results with a recent study by Armstrong and co-workers (31). They observe a similar modest,
3-fold, increase in the kobs for product release of a mutant relative to wild-type GSTM1-1, based on burst phase analysis of substrate turnover. The increase in the rate of product release is correlated with a disruption of an intramolecular hydrogen bonding network that restricts segmental motion. Upon disruption of these interactions, product release is accelerated modestly. The similar magnitude in rate increase for product release upon alteration of conformational dynamics of the two proteins provides a benchmark for the extent to which diffusive steps are controlled by enzymes; perhaps, very large rate enhancements are inaccessible by fine-tuning ligand binding and dissociation when they are coupled to conformational changes.
The new information obtained from our studies clarifies a functional role for Tyr-9 ionization in GSTA1-1 catalysis. To the extent that the Tyr-9 pKa remains low (
9.0) even with a GSH conjugate bound, a small fraction will be unprotonated even at physiological pH. Similarly, our previous results suggest a persistent disordered state, wherein even at saturating ligand and pH 7.4 the C-terminal helix is not completely closed and ordered. Thus, the equilibria involved yield small populations of deprotonated or disordered helix even at saturating concentration of conjugate. Apparently, several microconformational states are nearly isoenergetic. Our results indicate that several microconformations compose the closed state, and deprotonation of Tyr-9 alters this ensemble in favor of microconformations that promote faster C-terminal opening and release of product. As long as the microconformational states corresponding to the protonated and unprotonated Tyr-9 are in rapid equilibrium, even a small population of the unprotonated form could provide a kinetically favored, parallel, pathway for product release, as depicted in Fig. 7. Based on this, the ionized form of Tyr-9 provides a modest catalytic advantage by lowering the kinetic barrier to C-terminal opening and product release, and by decreasing the extent of product inhibition.
 |
FOOTNOTES
|
---|
* This work was supported by National Institutes of Health Grant GM62284 and the Department of Medicinal Chemistry at the University of Washington. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. 
¶ To whom correspondence should be addressed: Dept. of Medicinal Chemistry, Box 357610, University of Washington, Seattle, WA 98195-7610. Tel.: 206-685-0379l; Fax: 206-685-3252; E-mail: winky{at}u.washington.edu.
1 The abbreviations used are: GST, glutathione S-transferase; GS-, thiolate anion of glutathione;
, glutathione sulfonate; MES, 2-(N-morpholino)ethanesulfonic acid; CAPSO, 3-(cyclohexylamino)-2-hydroxyl-1-propanesulfonic acid; GS-EA, Michael adduct of ethacrynic acid and glutathione; rGST, GST derived from rats; hGST, GST derived from humans. 
 |
ACKNOWLEDGMENTS
|
---|
We thank Dr. Dharma Prabaharan for cloning the Y9F/W21F/F222W mutant employed in this study.
 |
REFERENCES
|
---|
- Mannervik, B., and Danielson, U. H. (1988) CRC Crit. Rev. Biochem. 23, 283337[Medline]
[Order article via Infotrieve]
- Hayes, J. D., and Pulford, D. J. (1995) Crit. Rev. Biochem. Mol. Biol. 30, 445600[Abstract]
- Armstrong, R. N. (1997) Chem. Res. Toxicol. 10, 218[CrossRef][Medline]
[Order article via Infotrieve]
- Board, P. G., Baker, R. T., Chelvanayagan, G., and Jermiin, L. S. (1997) Biochem. J. 328, 929935[Medline]
[Order article via Infotrieve]
- Board, P. G., Coggan, M., Chelvanayagam, G., Easteal, S., Jermiin, L. S., Schulte, G. K., Danley, D. E., Hoth, L. R., Griffor, M. C., Kamath, A. V., Rosner, M. H., Chrunyk, B. A., Perregaux, D. E., Gabel, C. A., Geoghegan, K. F., and Pandit, J. (2000) J. Biol. Chem. 275, 2479824806[Abstract/Free Full Text]
- Coles, B., and Ketterer, B. (1990) CRC Crit. Rev. Biochem. Mol. Biol. 25, 4767
- Tew, K. D. (1994) Cancer Res. 54, 43134320[Abstract]
- O'Brien, M., Kruh, G. D., and Tew, K. D. (2000) J. Pharmacol. Exp. Ther. 294, 480487[Abstract/Free Full Text]
- Atkins, W. M., Wang, R. W., Bird, A. W., Newton, D. J., and Lu, A. Y. H. (1993) J. Biol. Chem. 268, 1918819191[Abstract/Free Full Text]
- Bjornestedt, A., Stenberg, G., Wildersten, M., Board, P. G., Sinning, I., Jones, T. A., and Mannervik, B. (1995) J. Mol. Biol. 247, 765773[CrossRef][Medline]
[Order article via Infotrieve]
- Nilsson, L. O., Gustafsson, A., and Mannervik, B. (2000) Proc. Natl. Acad. Sci. U. S. A. 97, 94089412[Abstract/Free Full Text]
- Parsons, J. F., and Armstrong, R. N. (1996) J. Am. Chem. Soc. 118, 22952296[CrossRef]
- Thorson, J. S., Shin, I., Chapman, E., Stenberg, G., Mannervik, B., and Shultz, P. G. (1998) J. Am. Chem. Soc. 120, 451452[CrossRef]
- Nieslanik, B. S., Dabrowski, M. J., Lyon, R. P., and Atkins, W. M. (1999) Biochemistry 38, 69716980[CrossRef][Medline]
[Order article via Infotrieve]
- Nieslanik, B. S., and Atkins, W. M. (2000) J. Biol. Chem. 275, 17447174451[Abstract/Free Full Text]
- Nieslanik, B. S., Ibarra, C., and Atkins, W. M. (2001) Biochemistry 40, 35363543[CrossRef][Medline]
[Order article via Infotrieve]
- Allardyce, C. S., McDonagh, P. D., Lian, L.-Y., Wolf, C. R., and Roberts, G. C. K. (1999) Biochem. J. 343, 525531[CrossRef][Medline]
[Order article via Infotrieve]
- Atkins, W. M., Dietze, E. C., and Ibarra, C. (1997) Protein Sci. 6, 873881[Abstract/Free Full Text]
- Wang, R. W., Pickett, C. B., and Lu, A. Y. H. (1989) Arch. Biochem. Biophys. 269, 536543[Medline]
[Order article via Infotrieve]
- Ibarra, C., Nieslanik, B. S., and Atkins, W. M. (2001) Biochemistry 40, 1061410624[CrossRef][Medline]
[Order article via Infotrieve]
- Habig, W. H., Pabst, M. J., and Jakoby, W. B. (1974) J. Biol. Chem. 249, 71307139[Abstract/Free Full Text]
- Dirr, H. W., and Wallace, L. A. (1999) Biochemistry 38, 1563115640[CrossRef][Medline]
[Order article via Infotrieve]
- Lakowicz, J. R. (1999) Principles of Fluorescence Spectroscopy, 2nd Ed., pp. 156180, Kluwer Acadamic/Plenum Publishers, New York
- Wang, R. W., Bird, A. W., Newton, D. J., Lu, A. Y. H., and Atkins, W. M. (1993) Protein Sci. 2, 20852094[Abstract/Free Full Text]
- Adman, E. T., Le Trong, I., Stenkamp, R. E., Nieslanik, B. S., Dietze, E. C., Tai, G., Ibarra, C., and Atkins, W. M. (2001) Proteins Struct. Funct. Genet. 42, 192200[CrossRef][Medline]
[Order article via Infotrieve]
- Le Trong, I., Stenkamp, R. E., Ibarra, C., Atkins, W. M., and Adman, E. T. (2002) Proteins Struct. Funct. Genet. 48, 618627[CrossRef][Medline]
[Order article via Infotrieve]
- Gustafsson, A., Etahadieh, M., Jemth, P., and Mannervik, B. (1999) Biochemistry 38, 1626816275[CrossRef][Medline]
[Order article via Infotrieve]
- Caccuri, A. M., Antonini, G., Board, P. G., Parker, M. W., Nicotra, M., Lo Bello, M., Federici, G., and Ricci, G. (1999) Biochem. J. 344, 419425[CrossRef][Medline]
[Order article via Infotrieve]
- Caccuri, A. M., Lo Bello, M., Nuccetelli, M., Nicotra, M., Rossi, P., Antonini, G., Federici, G., and Ricci, G. (1998) Biochemistry 37, 30283034[CrossRef][Medline]
[Order article via Infotrieve]
- Widersten, M., Bjornestedt, R., and Mannervik, B. (1996) Biochemistry 35, 77317742[CrossRef][Medline]
[Order article via Infotrieve]
- Codreanu, S. G., Ladner, J. E., Xiao, G., Stourman, N. V., Hachey, D. L., Gilliland, G. L., and Armstrong, R. N. (2002) Biochemistry 41, 1516115172[CrossRef][Medline]
[Order article via Infotrieve]