|
Originally published In Press as doi:10.1074/jbc.M703414200 on September 11, 2007
J. Biol. Chem., Vol. 282, Issue 46, 33484-33493, November 16, 2007
Elucidation of a Complete Kinetic Mechanism for a Mammalian Hydroxysteroid Dehydrogenase (HSD) and Identification of All Enzyme Forms on the Reaction CoordinateTHE EXAMPLE OF RAT LIVER 3 -HSD (AKR1C9)*
William C. Cooper 1,
Yi Jin 1, and
Trevor M. Penning 2
From the
Center of Excellence in Environmental Toxicology, Department of Pharmacology, University of Pennsylvania School of Medicine, Philadelphia, Pennsylvania 19104-6084
Received for publication, April 24, 2007
, and in revised form, September 7, 2007.
 |
ABSTRACT
|
|---|
Hydroxysteroid dehydrogenases (HSDs) are essential for the biosynthesis and mechanism of action of all steroid hormones. We report the complete kinetic mechanism of a mammalian HSD using rat 3 -HSD of the aldo-keto reductase superfamily (AKR1C9) with the substrate pairs androstane-3,17-dione and NADPH (reduction) and androsterone and NADP+ (oxidation). Steady-state, transient state kinetics, and kinetic isotope effects reconciled the ordered bi-bi mechanism, which contained 9 enzyme forms and permitted the estimation of 16 kinetic constants. In both reactions, loose association of the NADP(H) was followed by two conformational changes, which increased cofactor affinity by >86-fold. For androstane-3,17-dione reduction, the release of NADP+ controlled kcat, whereas the chemical event also contributed to this term. kcat was insensitive to [2H]NADPH, whereas Dkcat/Km and the Dklim (ratio of the maximum rates of single turnover) were 1.06 and 2.06, respectively. Under multiple turnover conditions partial burst kinetics were observed. For androsterone oxidation, the rate of NADPH release dominated kcat, whereas the rates of the chemical event and the release of androstane-3,17-dione were 50-fold greater. Under multiple turnover conditions full burst kinetics were observed. Although the internal equilibrium constant favored oxidation, the overall Keq favored reduction. The kinetic Haldane and free energy diagram confirmed that Keq was governed by ligand binding terms that favored the reduction reactants. Thus, HSDs in the aldo-keto reductase superfamily thermodynamically favor ketosteroid reduction.
 |
INTRODUCTION
|
|---|
Mammalian hydroxysteroid dehydrogenases (HSDs)3 play pivotal roles in steroid hormone biosynthesis and metabolism (1, 2). In target tissues, they regulate the amount of steroid hormone available for their cognate nuclear receptor (3–7). This is achieved by the positional and stereospecific oxidoreduction of potent steroid hormones to their corresponding inactive metabolites or vice versa. Often these interconversions are catalyzed by pairs of HSDs that function preferentially as reductases or oxidases (5–8). There are two superfamilies of HSDs, the soluble aldo-keto reductase (AKRs) (9) and the mainly membrane-associated short-chain dehydrogenase reductases (SDRs) (10). Several HSDs have become drug targets to manipulate steroid hormone action in target tissues (11–13).
Four human HSDs from the AKR superfamily, AKR1C1–AKR1C4, can reduce 3-, 17-, and 20-ketosteroids in different ratios (14, 15). They catalyze the 4-pro-R hydride transfer from NADPH to the acceptor carbonyl on the steroid substrate in an ordered bi-bi reaction where NAD(P)(H) binds first and leaves last (16, 17). Rat liver 3 -HSD (AKR1C9), which shares 69% sequence identity with its human homologs, represents the best system to relate structure to function for HSDs, because crystal structures exist for the apoenzyme (E) (18), the E·NADP+ binary complex (19), and the E·NADP+·testosterone ternary complex (where testosterone is a competitive inhibitor) (20). These "snapshots" of the enzyme along the reaction pathway reveal that significant conformational changes occur upon the binding of each ligand and provide a structural basis for the ordered bi-bi mechanism (Fig. 1).
Like all AKRs, AKR1C9 adopts an ( / )8-barrel structure that has three large associated loops (loop A–C). In the binary complex, the cofactor is woven through the -helices by 22 contacts with the enzyme. Superimposition of the structures of the apoenzyme and the binary complex showed that the tight binding of the cofactor is associated with significant movement of loops A and B. Interestingly, the 1- 1 loop and part of loop B, and the C-terminal tail were disordered in the binary complex, but become ordered in the AKR1C9·NADP+·testosterone complex to form a mature steroid binding cavity. The influence of loop dynamics on catalysis remains unclear and raises the possibility that additional enzyme forms can exist on the reaction pathway. A complete detailed kinetic mechanism, including microscopic rate constants for conformational changes of all enzyme forms involved in ligand binding, and the chemistry of transformation have yet to be described for any HSD in real-time. Knowledge of the number of enzyme forms that exist is a pre-requisite to identify all the steps that could be intercepted by HSD inhibitors.
The hydride transfer reaction catalyzed by HSDs is reversible in vitro. The seemingly directional preference of an individual HSD, either an AKR or an SDR, seen in transfection studies is believed to be governed by the cofactor ratio and cofactor preference. In the cytosol, the relationships are [NADPH] >> [NADP+] and [NAD+] >> [NADH] (21, 22). Thus enzymes that have high affinity for NADPH will act as reductases, and those that preferentially utilize NADH/NAD will act as oxidases. AKRs have nanomolar affinity for the NADPH/NADP+ pair. We found that the NADPH-dependent reduction of 5 -DHT catalyzed by human type 3 3 -HSD (AKR1C2) occurred unimpeded in the presence of 1 mM NAD+ but that the NAD+-dependent oxidation of 3 -androstanediol oxidation catalyzed by AKR1C2 was potently inhibited by low micromolar concentrations of NADPH, suggesting that the reaction is unidirectional under normal cellular redox conditions (5). Recently, transfection studies with HSDs showed that a functional pseudo-equilibrium was reached such that the rates of the forward and reverse reactions were identical (23, 24). Thus, directionality and/or reversibility of an HSD are also related to the equilibrium constant of the reaction, which has generally not been reported for mammalian HSDs.
In this study we elucidate the complete kinetic mechanism for a mammalian HSD at physiological pH using AKR1C9. We find that the reaction sequence is governed by 9 enzyme forms, and 16 kinetic constants were estimated. We also find that the internal equilibrium constant for the reaction favors oxidation but this is not reflected in the overall Keq, which favors reduction. The kinetic Haldane and the energy diagram confirmed that Keq was governed by ligand binding terms that favored reactants in the reduction direction. Our findings are discussed within a physiological context.
 |
EXPERIMENTAL PROCEDURES
|
|---|
Materials—Cofactors were purchased from Roche Diagnostics. Steroids were purchased from Steraloids. Deuterated compounds for the synthesis of deuterated NADPH were purchased from Cambridge Isotope Laboratories. NADP+-specific alcohol dehydrogenase from Thermoanaerobium brockii was purchased from Sigma. Aldehyde dehydrogenase was purchased from Roche Biochemicals. The expression of pET16b-3 -HSD in Escherichia coli DE3 cells and purification of recombinant AKR1C9 have been described previously (25). A specific activity of 1.6 µmol of androsterone oxidized/min/mg of protein was determined. All other reagents were of ACS grade or better. Unless stated otherwise, experiments were carried out at 25 °C in 10 mM potassium phosphate buffer (pH 7.0) and 4% acetonitrile.
Steady-state Kinetic Analysis—Initial velocities were measured on a Hitachi F-4500 fluorescence spectrophotometer by monitoring NADPH depletion or production by the loss or gain of emission at 450 nm upon excitation at 340 nm.
Double reciprocal plots of families of lines in which the initial velocities were plotted against substrate concentrations, androstanedione 0.07–7 µM and NADPH 0.07–8 µM for reduction, and androsterone 0.1–30 µM and NADP+ 0.39–27 µM for oxidation, were constructed (26). Data were fitted to the minimal rate equation for an ordered bi-bi reaction, Equation 1 (also see Scheme 1),
 | (Eq. 1) | with the Cleland FORTRAN program and analyzed for a sequential mechanism. Patterns from plots of families of lines (see "Results") using the SEQUEN program were judged to be the best fit according to the criteria of Cleland (27).
The individual rate constants described in Scheme 1 were obtained using the measured steady-state kinetic parameters based on the definitions of these parameters in terms of rate constants for an ordered bi-bi reaction where the isomerization of the central complex is included in the reaction sequence (28). kAb, kAr, kQr, and kQb, which describe cofactor binding and release, were directly calculated. However, the remaining rate constants could only be derived by iteratively fitting six independent algebraic equations for each of the six steady-state parameters that were measured (i.e. kcatRed, KmNADPH, KmAndrostanedione, kcatOxi, KmNADP+, and KmAndrosterone). Initial estimates and constraints of rate constants were made based on results of transient kinetic experiments. (see "Transient State Kinetic Analysis"). Values of the rate constants were adjusted until all the fits converged simultaneously, i.e. the calculated values of all six steady-state parameters matched the measured values.
The synthesis of 4-pro-R-[2H]NADPH (NADPD) was performed as described by Hermes et al. (29) and Viola et al. (30) with minor modifications. Analysis of the NADPD product and assessment of deuterium incorporation was determined at 2.9 ppm by [1H]NMR. The [2H] content at the 4-pro-R position was estimated to be >99%. For estimates of primary KIEs, initial velocities were performed as described previously by substituting NADPD in the assays. The nomenclature of Northrop (31) was used where Dkcat is the ratio of kcat determined in the presence of NADPH versus NADPD, and Dkcat/Km is the ratio of kcat/Km determined in the presence of NADPH versus NADPD.
Measurement of Equilibrium Constant—The equilibrium constant of the reaction was determined experimentally using a method described by Talalay and Levy (32). In this method, the equilibrium was approached from both the reduction and oxidation directions. Experiments were performed at different cofactor and steroid concentrations. In a typical experiment, the initial reaction mixture contained 4.1 µg of AKR1C9, 22.6 µM androstanedione, 34.6 µM androsterone, and 23 µM NADPH (in 100 mM phosphate buffer, pH 7.0, 5% methanol at 25 °C). When the equilibrium was reached, 8.5 nmol of NADP+ was added to force the oxidative reaction and establish a new equilibrium. The changes in absorbance at 340 nm were monitored to calculate the concentrations of reactants and products at each equilibrium end-point permitting the calculation of Keq.
Transient State Kinetic Analysis—Transient state kinetic assays were performed on an Applied Photophysics SX.18MV-R stopped-flow reaction analyzer instrument fitted with a 20-µl flow cell and a dead time of 1 ms. All concentrations are given as final.
Kinetics of Binary Complex Formation—The binding of NADPH to AKR1C9 was monitored in real-time using stopped-flow spectroscopy by exciting the protein at 290 nm and capturing the emission of the energy transfer band at 450 nm through a 10-nm band-pass filter. The energy transfer band is produced as a result of the interaction of NADPH with Trp-86 (33). Conversely, to monitor the binding of NADP+, the decrease in protein fluorescence emission at 330 nm was recorded upon quenching the excited tryptophanyl residues at 290 nm. In a typical experiment, AKR1C9 (0.25 µM) was mixed with NADPH (2.0–11 µM). For NADP+ binding, enzyme (0.70 µM) was mixed with NADP+ (4.4–36 µM). The time courses of the fluorescence increase or decrease at each concentration of ligand were best fitted to double exponential equation (Equation 2),
 | (Eq. 2) | where Ft is the fluorescence at time t, F is the amplitude of the fluorescence increase or decrease, kobs1 and kobs2 are the apparent first-order rate constants for the fast and slow phases, respectively, and Feq is the fluorescence at equilibrium. Data were analyzed based on a three-step model for cofactor binding (Scheme 2). Estimates of the kinetic constants K1 NADPH and k3–k5 were obtained by the method of Stone and Le Bonniec (34). K1 NADPH is the apparent dissociation constant of the initial loose complex. The plot of kobs2 versus [NADPH] provided the estimate of k5. Values of KdNADPH were determined by steady-state fluorescence titration (25), and estimates of K1 and k3–k5 were substituted into Equation 3 to solve for k6.
 | (Eq. 3) | Similarly, estimates of kinetic constants for the binding of NADP+ (k13–k16 and K1NADP+) were obtained.
Ligand Chase Experiments—Enzyme (0.4 µM) was incubated with NADPH (12 µM) in one syringe and mixed with NADP+ (40 µM) from a second syringe. The decrease of the energy transfer band corresponding to dissociation of NADPH from the binary complex was monitored at 450 nm. Transients were best fitted to a bi-exponential (Equation 2).
Kinetics of Ternary Complex Formation—A solution containing enzyme (0.4 µM) and NADPH (12 µM) from one syringe was mixed with increasing concentrations of testosterone (3.60–6.15 µM) from a second syringe. Decreases in fluorescence at 450 nm were monitored with excitation at 295 nm. Fluorescence traces were marked by a single exponential decay. Five transient traces were averaged for each steroid concentration, and values for kobs were fitted to Equation 4 to solve for kon and koff (35).
 | (Eq. 4) |
Transient State Turnover Experiments—For single turnover experiments, AKR1C9 (0.5 µM) was incubated with substoichiometric amounts of NADPH (0.3 µM) in one syringe and mixed with a varying excess of androstanedione (4–15.0 µM) added from a second syringe in the stopped-flow instrument. Similarly for single turnover of androsterone oxidation, a solution of AKR1C9 (0.8 µM) and NADP+ (0.6 µM) from one syringe were mixed with an excess of androsterone ranging from 7.0 to 30 µM from a second syringe. Reactions were monitored fluorimetrically (excitation at 340 nm and emission at 450 nm). Data were fitted to a single exponential equation. Plots of kobs versus steroid concentration were fitted to the hyperbolic equation (Equation 5), where A is the steroid concentration, klim is the maximum rate of single turnover, and Kd is a quasi dissociation constant for steroid substrate (34).
 | (Eq. 5) | Values of klim provided the lower limits for the k+p and k-p and initial constraints on kBb, kBr, kPr, and kPb were set based on values of Kd during steady-state fitting analysis.
To study multiple turnover of androstanedione reduction, a solution of AKR1C9 (1.5 µM) and NADPH (25 µM) was mixed with androstanedione (13.5 µM). A complete concentration series with androstanedione was not possible due to steroid solubility constraints. Kinetic transients of NADPH depletion at 450 nm were collected. Experiments were repeated in the absorbance mode. Similarly, multiple turnover experiments of androsterone oxidation were conducted in which AKR1C9 (1 µM) and NADP+ (40 µM) from one syringe were mixed with increasing concentrations of androsterone (10–35 µM) added from a second syringe. Transient traces were fitted to an equation (Equation 6) that contained both burst and steady-state terms,
 | (Eq. 6) | where Ft equals fluorescence at time t, Amp equals the amplitude of the burst, kobs is the apparent rate constant of the burst in units of s-1, vss is the steady-state rate in units of F/ t, and c is the fluorescence at equilibrium (35). The maximum value of kobs yielded kburst, and the maximum value of vss was, subsequently, mathematically transformed into micromolar/second and divided by enzyme concentration to yield kss (s-1).
Transient state single and multiple turnover primary KIEs were determined by substituting NADPD in single turnover and multiple turnover experiments. Dklim is the ratio of klim determined in the presence of NADPH versus NADPD, and Dklim/Kd is the ratio of klim/Kd determined in the presence of NADPH versus NADPD. klim and Kd are given by Equation 5.
 |
RESULTS
|
|---|
Dissection of the complete kinetic mechanism for AKR1C9 was conducted using androstanedione and NADPH, and androsterone and NADP+, as substrate pairs at physiological pH. Androstanedione and androsterone were chosen, because these steroids have been standards with which all other assays on this enzyme in our laboratory have been compared. NADP(H) cofactors were chosen, because bound NADP+ was present in the crystal structures of the binary and ternary complexes. Thus, the existence of kinetically inferred complexes could be related to structural changes in the enzyme.
Steady-state Initial Velocity Studies—Initial velocity studies for androstanedione reduction and androsterone oxidation using NADP(H) were performed. Double reciprocal plots of the initial velocity data (see supplemental data) were characterized by a fan of lines that converged to an intersecting point to the left of the origin and above or on the x-axis. This is indicative of a sequential mechanism defined by the requirement of the cofactor and substrate to bind to the enzyme to form a ternary complex before any product can be produced (28). The sequential mechanism was previously shown to be an ordered bi-bi reaction as seen in Scheme 1. Fits of family of lines were generated through the Cleland SEQUEN program to yield steady-state kinetic parameters given in Table 1.
View this table:
[in this window]
[in a new window]
|
TABLE 1 Steady-state kinetic parameters for androstanedione reduction and androsterone oxidation catalyzed by AKR1C9
|
|
It was found that kcat = 0.82 s-1 for the oxidation direction was 2-fold greater than the reduction direction, kcat = 0.42 s-1, whereas KmNADPH was smaller than KmNADP+ by > 18-fold. Ki values for the cofactors were determined to be 0.17 µM and 1.51 µM for NADPH and NADP+, respectively, and were in reasonable agreement with the Kd values of 0.14 µM and 0.32 µM determined directly by fluorescence titration (25). The internal consistency of the kinetic parameters in Table 1 was checked with the kinetic Haldane relationship for an ordered bi-bi mechanism (28) by Equation 7.
 | (Eq. 7) | The steady-state parameters in Table 1 gave a value of 22.7 for this expression where the forward direction is reduction. This is in good agreement with a value of 20.9 ± 1.8 determined experimentally using a method described by Talalay and Levy (32). The kinetic Haldane revealed that the Keq is dominated by submicromolar Ki and Km values for NADPH and androstanedione, respectively.
Rate Constants Obtained in the Steady State—The individual rate constants described in Scheme 1 were calculated and derived from the measured steady-state parameters (Table 2). In the reduction direction these constants showed that the binding of NADPH (kAb) was slow and did not occur at the diffusion limit. The binding of androstanedione (kBb) to the binary complex occurred 5-times faster than kAb, and the rate of conversion of the central complex (k+p) was 6.6 s-1. The release of the steroid product androsterone (kPr) occurred at a comparable rate of 10.5 s-1. The main contribution to the overall kcat of 0.42 s-1 came from the release of the NADP+ product, kQr = 0.65 s-1.
In the oxidation direction, macroscopic rate constants showed that the binding of NADP+ (kQb) was 10-fold slower than the binding of NADPH (kAb), and the binding of androsterone (kPb) was 100-fold slower than the binding of androstanedione (kBb). However, the conversion of the central complex (k-p of 43.2 s-1) was >50-fold faster than the kcat of oxidation (0.82 s-1). The release of steroid product was also fast (kBr of 55.2 s-1), and the reaction was rate-controlled by the release of NADPH (kAr of 0.86 s-1). Because these rate constants may reflect more than one event, the reaction was further dissected using transient kinetic approaches.

View larger version (16K):
[in this window]
[in a new window]
|
FIGURE 2. Fluorescence kinetic transients for the binding of NADP+ to AKR1C9. A, a representative fluorescence transient trace for the binding of NADP+ (2.5 µM) to AKR1C9 (0.25 µM) fitted to a mono-exponential function; B, the same trace in A fitted to a bi-exponential function; The residuals are shown under A and B; C, concentration dependence of kobs (s-1) from the fast phase () and slow phase ( ) of the bi-exponential fittings.
|
|
Transient Kinetics of Cofactor Binding—The binding of NADP(H) to AKR1C9 is associated with the quenching of the intrinsic tryptophan fluorescence of the protein and yields a kinetic fluorescent kinetic transient that can be measured by stopped-flow kinetics (25). This kinetic transient was abolished by the R276M mutant and therefore assigned to cofactor anchoring by the guanidinium group of Arg-276. A two-step binding model was proposed that involved rapid formation of a loose complex followed by conformational changes to yield a tight binding complex that occurred as a result of the anchoring event. In re-examining this kinetic transient we found that these fluorescence kinetic transients can be best fitted to a bi-exponential decay rather than to a mono-exponential decay based on the shape and size of the residuals (Fig. 2). The differences in the apparent rate constants kobs1 and kobs2 for the fast and slow phases were 10-fold, and the amplitudes of the slow phase contributed to 25–30% of the total fluorescence change. Both kobs1 and kobs2 showed saturation kinetics versus the ligand concentration, thus a three-step binding model was proposed, which invoked the formation of a loose enzyme cofactor complex followed by two conformational changes (Scheme 2).
The kinetic rate constants defined in Scheme 2 were estimated using the transient data and are listed in Table 3. The dissociation constants for the initial loose complexes were determined to be 12 µM and 35 µM for E·NADPH and E·NADP+, respectively. Compared with the apparent Kd values determined experimentally by fluorescence titration, the sequential conformational change steps resulted in 86- and 110-fold increases in affinity for NADPH and NADP+ binding, respectively.
View this table:
[in this window]
[in a new window]
|
TABLE 3 The complete set of experimentally derived kinetic constants for the reduction of androstanedione by NADPH and oxidation of androsterone by NADP+
|
|
Ligand exchange experiments were performed on the premixed enzyme and NADPH solution by chasing the bound NADPH with excess NADP+ and measuring the decay in amplitude at 450 nm of the energy transfer band (see supplemental data). The kinetic traces were best fitted to a bi-exponential equation consistent with the three-step model for the release of cofactor.
Kinetics of Ketosteroid Binding—Because it is not possible to directly measure the kinetics of the binding of 3-ketosteroid steroid substrates we conducted experiments with testosterone, a competitive inhibitor and related 3-ketosteroid. Testosterone is also the ligand present in the crystal structure of the ternary complex. Quenching of the energy transfer band upon testosterone binding to the enzyme NADPH binary complex was marked by a decrease in fluorescence at 450 nm in the steady state and yielded a Kd value of 0.5 µM obtained by fluorescence titration.4 Stopped-flow experiments were then performed to determine whether steroid binding was accompanied with a fluorescence kinetic transient. Each transient trace was fitted to a single exponential equation (Fig. 3). At high testosterone concentrations the amplitude decayed too rapidly and the observable rate constant could not be discerned. As a result, saturation kinetics was not detected. At lower steroid concentrations, the concentration dependence on kobs was linear. The microscopic association and dissociation rate constants were estimated to be 95.7 ± 4.8 µM-1 s-1 and 14.1 ± 1.9 s-1, respectively, and gave a Kd of 0.15 µM.

View larger version (29K):
[in this window]
[in a new window]
|
FIGURE 3. Transient kinetics of testosterone binding. A, a representative fluorescence transient trace for the quenching of the energy transfer band at 450 nm upon mixing a preincubated solution of AKR1C9 (0.4 µM) and saturating NADPH (12 µM) with testosterone (3.6 µM); B, the dependence of kobs on varying testosterone concentrations.
|
|

View larger version (27K):
[in this window]
[in a new window]
|
FIGURE 4. Reduction of androstanedione by NADPH or NADPD and oxidation of androsterone by NADP+ under single turnover conditions. A, a representative fluorescence transient trace at 450 nm for the single turnover of androstanedione reduction obtained by mixing a solution of AKR1C9 (0.5 µM) and NADPH (lower trace) or NADPD (upper trace) (0.3 µM) with androstanedione (10 µM). The transient was fitted to a mono-exponential function to obtain kobs; B, the dependence of kobs on androstanedione concentration with NADPH () and NADPD ( ); C, a representative transient kinetic trace for the single turnover of androsterone oxidation obtained by mixing a solution of enzyme (0.8 µM) and cofactor (0.6 µM) with androsterone (20 µM). The transient was fitted to a mono-exponential function to obtain kobs; D, dependence of kobs on substrate concentration for androsterone oxidation.
|
|
Single Turnover Experiments—Transient kinetic traces of a decrease in NADPH fluorescence over time were recorded for androstanedione reduction catalyzed under single turnover conditions (Fig. 4). Rate constant kobs values were obtained from fitting transients to a single exponential equation. As predicted, the amplitude of the signal was directly proportional to the concentration of NADPH oxidized. Subsequently, kobs values were plotted against [androstanedione]. The plot was characterized by a hyperbola whereby: 1) the plateau yielded a klim value that estimates a maximal rate of single turnover and 2) substrate concentration at half the maximum klim yields an apparent Kd (Fig. 4B). All klim and Kd values (from transient as well as steady-state experiments) are given in Table 4. klim for androstanedione reduction approximated a value of 4.5 s-1 that was 10-fold >kcat (0.42 s-1) and, therefore, suggested that a step other than chemistry was rate-limiting. The klim value was in good agreement with the rate of conversion of the central complex k+p = 6.6 s-1 observed in the steady-state solution. Under single turnover conditions a klim of 35.7 s-1 was estimated for androsterone oxidation, which was much greater than kcat (0.82 s-1) and suggested that in the oxidation direction another step was also rate-determining. The klim value was consistent with the rate of conversion of the central complex k-p = 43.2 s-1 obtained in the steady-state solution (Fig. 4, C and D).
Multiple Turnover Reduction Experiments—Transient state multiple turnover experiments were also performed to see whether androstanedione reduction was accompanied by burst-phase kinetics. Transient traces were characterized by a small but yet discernable burst that was not stoichiometric with enzyme concentration (Fig. 5A). To validate this observation multiple turnover experiments for androstanedione reduction were also performed in the absorbance mode. A transient trace obtained over 5 s was characterized by a burst phase of distinguishable amplitude followed by a steady-state phase as illustrated in Fig. 5B. Fitting the traces to Equation 6 gave an apparent kburst of 5.6 s-1 and kss of 0.41 s-1, which agreed well with a klim of 4.5 s-1 and a kcat of 0.42 s-1, respectively.
Transient-state multiple turnover experiments were also performed to determine whether androsterone oxidation was accompanied by burst-phase kinetics. Androsterone oxidation was characterized by a pronounced burst that yielded a kburst (55.8 s-1) that was similar to klim (35.7 s-1) (Fig. 5C). The amplitude of the burst coincided with the concentration of the enzyme. The kss (0.77 s-1) from androsterone oxidation under multiple turnovers was in accord with kcat (0.82 s-1).
Primary KIEs on Steady-state and Transient State Kinetic Parameters—To confirm that the chemical step was not rate-determining for androstanedione reduction, primary KIE values were performed using NADPD. There was no significant Dkcat nor Dkcat/Km, and values of 1.08 and 1.06 were obtained (Table 4). By contrast under single turnover conditions, Dklim = 2.06 was observed (Fig. 4). Under multiple turnover conditions a Dkburst equal to 2.07 was noted and was reduced to 0.93 for Dkss (Fig. 5A). The data were internally consistent with steps other than chemistry being rate-determining in the reduction of androstanedione.
 |
DISCUSSION
|
|---|
A complete kinetic mechanism for a mammalian HSD has been presented (see Table 3 and Scheme 3). The mechanism contains 9 enzyme forms, and 16 individual kinetic constants were assigned. The steady-state kinetic parameters for the mechanism were reconciled with rate constants obtained from the transient-state kinetic measurements. This mechanism can be used as a guide for studies on other HSDs.
Rate-determining Steps in AKR1C9 Catalysis—Our study focused on the NADP(H)-dependent interconversion of androstanedione and androsterone catalyzed by AKR1C9 using steady-state and transient-state kinetics. The contributions of different steps to rate determination and hence kcat can be dissected with Equations 8, 9, 10, 11. In the reduction direction, the release of NADP+ (0.65 s-1) largely (64%) controls the kcat of 0.42 s-1. However, the net rate for the chemical reduction of androstanedione (kchemRed of 1.3 s-1) also contributed significantly (33%), and this was due to the size of k-p for the oxidation reaction in the steady state. The small difference between kchemRed and the net rate of the subsequent product release steps (0.61 s-1) resulted in the partial burst seen under multiple turnover conditions. In contrast, the kcat of 0.82 s-1 for the oxidation direction was dominated (96%) by the release of NADPH (0.85 s-1). The 30-fold greater kchemOxi (24.2 s-1 caused the full burst observed under multiple turnover conditions for the oxidation reaction.

View larger version (18K):
[in this window]
[in a new window]
|
FIGURE 5. Reduction of androstanedione by NADPH or NADPD and oxidation of androsterone by NADP+ under multiple turnover conditions. A, a premixed solution of AKR1C9 (1.5 µM) and NADPD (upper trace) or NADPH (lower trace) (25 µM) was mixed with androstanedione (13.5 µM). Product formation corresponding to the decrease in fluorescence emission at 450 nm was monitored; B, a representative transient trace obtained in the absorbance mode where a solution of enzyme (2 µM) and NADPH (30 µM) were mixed with androstanedione (16 µM). The decrease in absorbance at 340 nm was monitored. C, a representative transient kinetic trace for androsterone oxidation under multiple turnover conditions where AKR1C9 (1 µM) and NADP+ (40 µM) were mixed with androsterone (15 µM).
|
|
 | (Eq. 8) |
 | (Eq. 9) |
 | (Eq. 10) |
 | (Eq. 11) | We have previously shown that the NADPH reduction of 5 -DHT catalyzed by AKR1C9 had a distinctively different kinetic profile (36). In this reaction, there was no burst-phase kinetics in multiple turnover experiments, and a presence of significant primary KIEs provided evidence that chemistry was a major rate-determining step. Thus AKR1C9 can catalyze identical reactions on two highly related 3-ketosteroids (5 -DHT and androstanedione) with different rate-determining profiles. This is likely due to the different positioning of the steroid at the active site and changes in proximity between the reactants. We have previously shown using alanine-scanning mutagenesis that subtle changes in the binding of steroids to AKR1C9 can alter the rate-determining step of the reaction (37).
Loop Dynamics—Fluorescence kinetic transients were observed for the binding of NADP(H) and the competitive inhibitor testosterone, and these may report the conformational changes that are observed upon ligand binding in the crystal structures of AKR1C9 and its complexes. The binding of cofactor to AKR1C9 was dissected into a three-step process that can be related to conformational changes. The formation of the initial loose complex E·NADP(H) was followed by two subsequent conformational changes that result in a tight binding complex. The first conformational change may be attributed to the formation of an electrostatic linkage between Arg-276 and the 2'-phosphate of AMP, because the fast phase of the kinetic transient can be eliminated by the R276M mutant (25). In the second conformational change step that yields the tight E**·NADP(H) complex, the on rate is 20 fold higher than the off rate. This step may be related to the formation of the cofactor tunnel that involves the residues on loop 1- 1 and loop B, the nicotinamide ring being sandwiched between Trp-86 and Tyr-216, and the formation of hydrogen bonds with the C3 carboxamide group that enable the head group to orient correctly. The lower off rates suggest that the disruption of the tunnel and binding of the head group is slow.
To progress from the binary complex to the ternary complex, structural data imply that additional ordering of the loops occurs (19). Kinetic transients were captured possibly reflecting these changes for the binding of testosterone to AKR1C9. The microscopic association and dissociation rate constants for testosterone binding were estimated to be 95.7 µM-1s-1 and 14.1 s-1 and approached those for androstanedione obtained in the steady-state solution, which were 55.2 µM-1s-1 and 27.8 s-1. Significantly, these values differ from the association and dissociation rate constants for androsterone obtained from steady-state analysis, which were 0.26 µM-1s-1 and 10.5 s-1, respectively. This reflects a 80-fold decrease in the Kd value for androstanedione. Considering the structural differences between androstanedione (substrate) and androsterone (product) in the reduction direction, it is apparent that the enzyme has evolved to favor 3-ketosteroid binding, i.e. the binding of reductive substrates.
Equilibrium Constants and Free Energy Diagram—The concept that HSDs work in pairs to act as either reductases or oxidases and thereby regulate ligand access to steroid hormone receptors is based to a large extent on transfection assays using radiotracers. In transfection paradigms HSDs are forced to use the prevailing concentrations of reduced or oxidized cofactor that exist in the recipient cell line. By contrast HSDs in vitro can be bi-directional providing that the NAD(P)(H) concentration is manipulated to drive the reaction in the desired direction. There have been few attempts to reconcile these disparate findings with the Keq of these reactions. In this study we measured a Keq of 20.9 for the NADPH-dependent reduction of androstanedione, which is in close agreement with that determined by the kinetic Haldane, which gave a value of 22.7. These values also agree with the Keq of 14.4 for the NADPH-dependent reduction of 5 -DHT catalyzed by AKR1C9 (36) and the Keq of 12.0 for the NADH-dependent reduction of androstanedione catalyzed by AKR1C9 (17). Thus, irrespective of the cofactor pair, AKR1C9 behaves preferentially as a reductase.

View larger version (18K):
[in this window]
[in a new window]
|
FIGURE 6. Free energy diagram for the NADPH-dependent reduction of androstanedione catalyzed by AKR1C9. Changes of the Gibbs free energy associated with each step described in Scheme 1 for the NADP(H)-dependent interconversion of androstanedione and androsterone were calculated according to Equations 12 and 13.
|
|
Our kinetic mechanism allows us to depict the HSD reaction in a free-energy diagram (Fig. 6). The difference in free energy for each step G and the activation free energy G were calculated using Equations 12 and 13,
 | (Eq. 12) |
 | (Eq. 13) | where kf and kr are the forward and reverse rate constants of the step, respectively, R is the gas constant, T is the absolute temperature, and kB and h are the Boltzmann and Planck constants, respectively. The free-energy profile was similar for both directions. In the reduction direction, the binding of NADPH and androstanedione are associated with free energy changes ( 1 and 2) of -9.2 and -7.8 kcal/mol, respectively. However, the free-energy change for the chemical conversion ( 3) is only +1.1 kcal/mol, while the energy required to reach the transition state for the reaction ( 3) was +16.3 kcal/mole. The subsequent gain in free energy for the product releases steps were 4 of +6.0 kcal/mol and 5 of +8.1 kcal/mol.
The overall free-energy change ( total) for the reaction is -1.8 kcal/mol, indicating that the reduction is favored thermodynamically. However, the internal Keq based on single turnover klimRed:klimOxi (0.13) or the steady-state solution k+p:k-p (0.15) indicates that the oxidation reaction is favored. By contrast, the Keq for the entire action is 22.7 and is 170 times larger than the internal Keq. The reason for this large change can be accounted for by ligand binding terms, which favor NADPH and androstanedione. In comparison to other members of the AKR superfamily, the Keq of AKR1C9 is low. Thus aldose reductase, a prototypical AKR, has a Keq of 150 in favor of reduction (38). The Keq values and internal Keq values for other mammalian HSDs, including the SDRs, have yet to be reported.
Physiological Significance—In humans four highly related AKR1C isoforms exist (AKR1C1–AKR1C4). We have reported a kinetic analysis for the NADPH-dependent reduction of 5 -DHT catalyzed by AKR1C2, which allows calculation of a Keq of 8.0 for its reduction reaction (39). When combined with the data reported here, it is clear that all AKR1C isoforms will preferentially act as reductases. Thus the human enzymes will function as 3-, 17-, and 20-ketoreductases. The unexpected finding of this study is the Keq does not appear to be governed by cofactor, because the same Keq is observed irrespective of whether reduction is coupled to NADPH/NADP+ or NADH/NAD+. Rather, the Keq is governed by favorable binding terms for the reduction reactants. Also, the Keq is not overwhelmingly in favor of reduction. The Keq values reported are rather modest and indicate that under some conditions the oxidative reaction may be observed. Two observations suggest that this can happen. Transfection of AKR1C9 into HEK293 cells showed that a pseudo-equilibrium was reached in which the product to substrate ratio was 97:3 (24). At this equilibrium the forward and reverse reaction rates were identical. Second, in studies on the metabolism of 3 -hydroxysteroids in human hepatoma cells we noted that, although the favored metabolic route was formation of steroid conjugates, a small amount of epimerization (<5%) to the 3 -hydroxysteroid metabolites was observed that could only be explained by oxidation back to the 3-ketosteroid (40).
 |
CONCLUSIONS
|
|---|
AKR1C9, an HSD, is a member of the steroid hormone-transforming AKRs that regulate ligand occupancy of steroid receptors. Understanding the number of enzyme forms and their effects on rate determination of steroid transformation identifies steps for interception by inhibitors and regulation of these reactions. This study revealed that in the reaction sequence there are nine enzyme forms and six contain NAD(P)H. It also showed that, although the internal Keq of AKR1C9 clearly favors oxidation, which is unique among AKRs, the overall Keq favors reduction. This difference in AKRs is achieved by slow product release steps that control the reaction in both directions. The Keq that favors reduction is dominated by ligand binding terms for the reduction reactants.
 |
FOOTNOTES
|
|---|
* This work was supported by National Institutes of Health Grants R01-DK47015 and P30 ES013508 (to T. M. P.). 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. 
The on-line version of this article (available at http://www.jbc.org) contains supplemental Figs. S1 and S2. 
1 Both authors contributed equally to this work. 
2 To whom correspondence should be addressed: Dept. of Pharmacology, University of Pennsylvania, 130C John Morgan Bldg., 3620 Hamilton Walk, Philadelphia, PA 19104-6084. Tel.: 215-573-2236; Fax: 215-573-2236; E-mail: penning{at}pharm.med.upenn.edu.
3 The abbreviations used are: HSD, hydroxysteroid dehydrogenase; AKR, aldoketo reductase; AKR1C2, human type 3 3 -HSD; AKR1C9, rat liver 3 -HSD; androstanedione, androstane-3,17-dione; androsterone, 17 -hydroxy-5 -androstane-3-one; 5 -DHT, 5 -dihydrotestosterone; KIE, kinetic isotope effect; SDR, short-chain dehydrogenase reductase; NADPD, [2H]NADPH. 
4 W. C. Cooper, unpublished results. 
 |
ACKNOWLEDGMENTS
|
|---|
We thank Jason E. Drury and Carol A. Shultz for providing purified AKR1C9 and Vladi V. Heredia for helpful discussions.
 |
REFERENCES
|
|---|
- Penning, T. M. (1997) Endocr. Rev. 18, 281-305[Abstract/Free Full Text]
- Payne, A. H., and Hales, D. B. (2004) Endocr. Rev. 25, 947-970[Abstract/Free Full Text]
- Funder, J. W., Pearce, P. T., Smith, R., and Smith, A. I. (1988) Science 242, 583-585[Abstract/Free Full Text]
- Draper, N., and Stewart, P. M. (2005) J. Endocrinol. 186, 251-271[Abstract/Free Full Text]
- Ri
ner, T. L., Lin, H.-K., Peehl, D. M., Steckelbroeck, S., Bauman, D. R., and Penning, T. M. (2003) Endocrinology 144, 2922-2932[Abstract/Free Full Text] - Bauman, D. R., Steckelbroeck, S., Williams, M. V., Peehl, D. M., and Penning, T. M. (2006) Mol. Endocrinol. 20, 444-458[Abstract/Free Full Text]
- Penning, T. M., Bauman, D. R., Jin, Y., and Ri
ner, T. L. (2007) Mol. Cell. Endocrinol. 265-266, 77-82[CrossRef][Medline]
[Order article via Infotrieve] - Liden, M., Tryggvason, K., and Eriksson, U. (2003) Mol. Aspects of Med. 24, 403-409[CrossRef]
- Jez, J. M., Bennett, M. J., Schlegel, B. P., Lewis, M., and Penning, T. M. (1997) Biochem. J. 326, 625-636[Medline]
[Order article via Infotrieve]
- Jörnvall, H., Persson, B., Krook, M., Atrian, S., Gonzàlez-Duarte, R., Jeffery, J., and Ghosh, D. (1995) Biochemistry 34, 6003-6013[CrossRef][Medline]
[Order article via Infotrieve]
- Penning, T. M. (2003) Hum. Reprod. Update 9, 193-205[Abstract/Free Full Text]
- Ogg, D., Elleby, B., Norstrom, C., Stefansson, K., Abrahmsen, L., Oppermann, U., and Svensson, S. (2005) J. Biol. Chem. 280, 3789-3794[Abstract/Free Full Text]
- Hosfield, D. J., Wu, Y., Skene, R. J., Hilgers, M., Jennings, A., Snell, G. P., and Aertgeerts, K. (2005) J. Biol. Chem. 280, 4639-4648[Abstract/Free Full Text]
- Penning, T. M., Burczynski, M. E., Jez, J. M., Hung, C.-F., Lin, H.-K., Ma, H., Moore, M., Palackal, N., and Ratnam, K. (2000) Biochem. J. 351, 67-77[CrossRef][Medline]
[Order article via Infotrieve]
- Steckelbroeck, S., Jin, Y., Gopishetty, S., Oyesanmi, B., and Penning, T. M. (2004) J. Biol. Chem. 279, 10784-10795[Abstract/Free Full Text]
- Ricigliano, J. W., and Penning, T. M. (1990) Biochem. J. 269, 749-755[Medline]
[Order article via Infotrieve]
- Askonas, L. J., Ricigliano, J. W., and Penning, T. M. (1991) Biochem. J. 278, 835-841[Medline]
[Order article via Infotrieve]
- Hoog, S. S., Pawlowski, J. E., Alzari, P. M., Penning, T. M., and Lewis, M. (1994) Proc. Natl. Acad. Sci. U. S. A. 91, 2517-2521[Abstract/Free Full Text]
- Bennett, M. J., Schlegel, B. P., Jez, J. M., Penning, T. M., and Lewis, M. (1996) Biochemistry 35, 10702-10711[CrossRef][Medline]
[Order article via Infotrieve]
- Bennett, M. J., Albert, R. H., Jez, J. M., Ma, H., Penning, T. M., and Lewis, M. (1997) Steroids 5, 799-812
- Veech, R. L., Eggleston, I. V., and Krebs, H. A. (1969) Biochem. J. 115, 609-619[Medline]
[Order article via Infotrieve]
- Williamson, D. H., Lund, P., and Krebs, H. A. (1967) Biochem. J. 103, 514-527[Medline]
[Order article via Infotrieve]
- Khun, N., Sharma, K. K., Andersson, S., and Auchus, R. J. (2004) Arch. Biochem. Biophys. 429, 50-59[CrossRef][Medline]
[Order article via Infotrieve]
- Papari-Zareei, M., Brandmaier, A., and Auchus, R. J. (2006) Endocrinology 147, 1591-1597[Abstract/Free Full Text]
- Ratnam, K., Ma, H., and Penning, T. M. (1999) Biochemistry 38, 7856-7864[CrossRef][Medline]
[Order article via Infotrieve]
- Cleland, W. W., (1979) Methods Enzymol. 63, 103-138[Medline]
[Order article via Infotrieve]
- Cleland, W. W. (1977) Adv. in Enzymol. 45, 273-387
- Segel, I. H. (1993) Enzyme Kinetics: Behavior and Analysis of Rapid Equilibrium and Steady-State Enzyme Systems, Wiley & Sons, New York, pp. 274-574
- Hermes, J. D., Morrical, S. W., O'Leary, M. H., and Cleland, W. W. (1984) Biochemistry 23, 5479-5488[CrossRef][Medline]
[Order article via Infotrieve]
- Viola, R. E., Cook, P. F., and Cleland, W. W. (1979) Anal. Biochem. 96, 334-340[CrossRef][Medline]
[Order article via Infotrieve]
- Northrop, D. B. (1982) Methods Enzymol. 87, 607-625[Medline]
[Order article via Infotrieve]
- Talalay, P., and Levy, H. R. (1959) Ciba Foundation Study Group No. 2: Steric Course of Microbiological Reactions, pp. 64-66, Little, Brown & Co., Boston
- Jez, J. M., Schlegel, B. P., and Penning, T. M. (1996) J. Biol. Chem. 271, 30190-30198[Abstract/Free Full Text]
- Stone, S. R., and Le Bonniec, B. F. (1997) J. Mol. Biol. 265, 344-362[CrossRef][Medline]
[Order article via Infotrieve]
- Johnson, K. A. (1992) The Enzymes 20, 1-62
- Heredia, V. V., and Penning, T. M. (2004) Biochemistry 43, 12028-12037[CrossRef][Medline]
[Order article via Infotrieve]
- Heredia, V. V., Cooper, W. C., Kruger, R. G., Jin, Y., and Penning, T. M. (2004) Biochemistry 43, 5832-5841[CrossRef][Medline]
[Order article via Infotrieve]
- Grimshaw, C. E., Bohren, K. M., Lai, C.-J., and Gabbay, K. H. (1995) Biochemistry 34, 14356-14365[CrossRef][Medline]
[Order article via Infotrieve]
- Jin, Y., and Penning, T. M. (2006) Biochemistry 45, 13054-13063[CrossRef][Medline]
[Order article via Infotrieve]
- Steckelbroeck, S., Oyesanmi, B., Jin, Y., Lee, S. H., Kloosterboer, H. J., and Penning, T. M. (2006) J. Pharmacol. Exp. Ther. 316, 1300-1309[Abstract/Free Full Text]

CiteULike Complore Connotea Del.icio.us Digg Reddit Technorati What's this?
This article has been cited by other articles:

|
 |

|
 |
 
Y. Jin, L. Duan, S. H. Lee, H. J. Kloosterboer, I. A. Blair, and T. M. Penning
Human Cytosolic Hydroxysteroid Dehydrogenases of the Aldo-ketoreductase Superfamily Catalyze Reduction of Conjugated Steroids: IMPLICATIONS FOR PHASE I AND PHASE II STEROID HORMONE METABOLISM
J. Biol. Chem.,
April 10, 2009;
284(15):
10013 - 10022.
[Abstract]
[Full Text]
[PDF]
|
 |
|
Copyright © 2007 by the American Society for Biochemistry and Molecular Biology.
|
Advertisement
Advertisement
|