|
|
||||||||
-hydroxysteroid dehydrogenase reveals an "induced-fit" mechanism and a conserved basic motif involved in the binding of androgen
Oncology and Molecular Endocrinology Research Center, Laval University Medical Center (CHUL) and Laval University, Québec, G1V 4G2, Canada
Reprint requests to: Rock Breton, Centre de Recherche en Endocrinologie Moléculaire et Oncologique, Centre Hospitalier de lUniversité Laval, 2705, boul. Laurier, Ste-Foy, QC, G1V 4G2, Canada; e-mail: rock.breton{at}crchul.ulaval.ca; fax: (418) 654-2761.
(RECEIVED January 11, 2005; FINAL REVISION March 17, 2005; ACCEPTED March 18, 2005)
| Abstract |
|---|
|
|
|---|
-hydroxysteroid dehydrogenase (h3
HSD3, AKR1C2) plays a crucial role in the regulation of the intracellular concentrations of testosterone and 5
-dihydrotestosterone (5
-DHT), two steroids directly linked to the etiology and the progression of many prostate diseases and cancer. This enzyme also binds many structurally different molecules such as 4-hydroxynonenal, polycyclic aromatic hydrocarbons, and indanone. To understand the mechanism underlying the plasticity of its substrate-binding site, we solved the binary complex structure of h3
HSD3-NADP(H) at 1.9 Å resolution. During the refinement process, we found acetate and citrate molecules deeply engulfed in the steroid-binding cavity. Superimposition of this structure with the h3
HSD3-NADP(H)-testosterone/acetate ternary complex structure reveals that one of themobile loops forming the binding cavity operates a slight contraction movement against the citrate molecule while the side chains of many residues undergo numerous conformational changes, probably to create an optimal binding site for the citrate. These structural changes, which altogether cause a reduction of the substrate-binding cavity volume (from 776 Å3 in the presence of testosterone/acetate to 704 Å3 in the acetate/citratecomplex), are reminiscent of the "induced-fit" mechanism previously proposed for the aldose reductase, another member of the AKR superfamily. We also found that the replacement of residues Arg301 and Arg304, localized near the steroid-binding cavity, significantly affects the 3
HSD activity of this enzyme toward 5
-DHT and completely abolishes its 17
HSD activity on 4-dione. All these results have thus been used to reevaluate the binding mode of this enzyme for androgens. Keywords: aldo-keto reductase; hydroxysteroid dehydrogenase; crystal structure; induced-fit mechanism
Abbreviations: h3
HSD3, human 3
-hydroxysteroid dehydrogenase type 3 h20
HSD, human 20
-hydroxysteroid dehydrogenase h3
HSD1, human 3
-hydroxysteroid dehydrogenase type 1 17
HSD, 17
-hydroxysteroid dehydrogenase NADP(H), reduce form of nicotinamide adenine dinucleotide phosphate GST, glutathione sulfotransferase PDB, Protein Data Bank AKR, aldo-keto reductase progesterone, 5-pregnen- 3,20-dione 4-dione, 4-androsten-3,20-dione PAH, polycyclic aromatic hydrocarbon EDTA, ethylene diamine tetraacetic acid.
Article and publication are at http://www.proteinscience.org/cgi/doi/10.1110/ps.051353205.
| Introduction |
|---|
|
|
|---|
HSD (h3
HSD3; AKR1C2; EC 1.1.1.213
[EC]
) is able to exert its activity on structurally different steroids, such as androgens (C19-steroids) and progestins (C21-steroids). In addition to steroids (Dufort et al. 1996, 2001; Penning et al. 2000), this enzyme is also able to bind/transform numerous molecules with various shapes, such as acetate (Nahoum et al. 2001), prostaglandins (Lovering et al. 2004), polycyclic aromatic hydrocarbons (PAHs) (Palackal et al. 2001, 2002), 1-indanol (Matsuura et al. 1997), 9, 10-phenanthrenequinone (Penning et al. 1984), and 4-hydroxy-2- nonenal (Burczynski et al. 2001) (Fig. 1B
HSD3 must adapt the structure of its binding cavity in order to form an optimal binding site for each of them. This capacity to alter the orientation of numerous side chains surrounding the ligand-binding site, which was first described for the yeast hexokinase, is named the "induced-fit" mechanism (DelaFuente and Sols 1970; DelaFuente et al. 1970). The induced-fit mechanism is also known for members of the AKR superfamily and has been well studied, notably in the case of the aldose reductase, in a strategy for the design of inhibitors (Klebe et al. 2004; Sotriffer et al. 2004).
|
HSDs. Rat 3
HSD structure (r3
HSD, AKR1C9) has been solved in ternary complex with testosterone (the product of the 17
-reduction of 4-dione) (Hoog et al. 1994; Bennett et al. 1996, 1997) (Fig. 1A
HSD3 has been determined in ternary complex with testosterone/acetate (Nahoum et al. 2001) and ursodeoxycholate (Jin et al. 2001), a bile acid possessing a steroid-like architecture. These enzymes are folded as triose-phosphate isomerase barrels (
/
)8 (Hoog et al. 1994). They maintain the cofactor in an extended conformation (Bennett et al. 1996) and stabilize the steroids through numerous contacts, mainly with hydrophobic residues belonging to three large loops named A, B, and C. Biochemical studies performed on r3
HSD demonstrated that the reaction catalyzed by this enzyme is triggered by the 4-Pro-R hydride transfer from the C-4 of the nicotinamide head of the cofactor to the carbon atom of the ketone group, thus reducing the steroid substrate. The reduction reaction is completed by a proton transfer from Tyr55 residue, which acts as a proton donor, to the oxygen atom through a "pushpull" mechanism (Schlegel et al. 1998). Recent studies have suggested that the orbital overlapping produced by an optimal orientation of the reacting atoms could play a major quantitative role in the catalytic power of these enzymes (Mesecar et al. 1997; Heredia et al. 2003). They determined that this optimal orientation is obtained when the carbon atom of the target ketone group is maintained at about 3.0 Å of the C-4 of the pyridine head of the cofactor and when a line joining these two atoms makes an angle of about 100° with the pyridine head (Fig. 1C
Considering these recent findingsthe pronounced induced-fit adaptations shown for at least one member of the AKR enzyme family and the precise localization of the ketone group to be reduced inside the active siteand taking advantage of ternary complex crystal structures we have determined so far for many AKR1C enzymes (Nahoum et al. 2001; Couture et al. 2003, 2004), we decided to revisit all these structures, especially those for which the steroid substrate was found in the steroid-binding cavity but at a nonproductive position (Nahoum et al. 2001; Couture et al. 2003, 2004). We wished to investigate some pending issues: (1) the mechanism by which h3
HSD3 binds/transforms numerous structurally different molecules, which implies an ability to discriminate various ligands from different steroid hormone classes; (2) the orientation adopted by a steroid in a productive complex; and (3) the geometry of the residues that surround the ligand in the productive complexes.
Here we report the crystal structure of h3
HSD3 in binary complex with NADP(H), a complex in which a citrate and an acetate molecule, added during the crystallization step, were found deeply bound inside the steroid-binding cavity. Comparison with the testosterone/acetate ternary complex structure reveals that loop A undergoes a slight contraction toward the citrate. Moreover, numerous residues found in the cavity, namely Val54, Val128, Ile129, Trp227, Leu306, and Leu308, alter the orientation of their side chains, narrowing the cavity. All these structural changes are reminiscent of the induced-fit mechanism previously proposed for the aldose reductase (Sotriffer et al. 2004). We also found that, in order to respect the orbital steering constraint, the optimal position of the ketone of a steroid during the reduction reaction is somewhat similar to that of the acetate molecule found in the active site. Manual positioning of 4-dione and DHT in an ideal productive ternary complex revealed that the nucleus of the steroid must be oriented toward the C-terminal domain (loop C) of the protein. To verify this, we proceeded to a leucine-scanning mutagenesis study on many residues of loop C, and found that the replacement of residue Arg301 or Arg304, located 10 Å from the catalytic Tyr55, significantly affects the activity of this enzyme. These results suggest that residues involved in the formation of a mature binding cavity in h3
HSD3 differ according to the nature of the steroid substrate, the position of the ketone group on the steroid nucleus, and the activity exerted by this enzyme.
| Results and Discussion |
|---|
|
|
|---|
HSD3 binary complex
HSD3/NADPH binary complex was refined to a crystallographic R-factor of 17.4% (Rfree=19.9%) at 1.9 Å of resolution. No major differences were found when both monomers were compared (root-mean-square deviation, calculated for 327 C
atoms is 0.36 Å). The model obtained was of good quality since more than 99.7% of the dihedral angles of the main chain were found in the allowed region of the Ramachandran plot (Table 1
/
)8-barrel], in which eight
strands forming the cylindrical core of the barrel are surrounded by eight helices. The first 17 residues are folded as a two-strands anti-parallel beta-sheet that closes the bottom of the barrel. Finally, three large loops, named A (Phe118 to Asp143), B (Ser217 to Asp238), and C (Leu299 to Tyr323) are situated in a manner that creates a funnel-shaped pocket at the carboxy end of the
-strands, the top of the barrel (Fig. 2A
|
|
HSD3. Because no cofactor was added to the protein preparation during the purification and crystallization steps, we suppose that the NADP(H) found in all our crystal structures came from the Escherichia coli cells used to overproduce the h3
HSD3 protein. The fact that the endogenous bacterial NADP(H) was kept by h3
HSD3 through all the purification steps is probably an indication of the high affinity of this enzyme for its cofactor.
Steroid-binding site
During the refinement process, two unexpected positive peaks appeared inside the steroid-binding cavity of monomer A. Citrate and acetate molecules, added to the mother liquor during the crystallization step, were fitted into the FoFc difference Fourier map and refined with the rest of the model (Fig. 2B
). These two molecules were surrounded by residues Tyr24, Ala25, Val54, Tyr55, Trp86, His117, Ile129, Asn167, Tyr216, His222, Glu224, Trp227, Leu306, Leu308, Ile310, and Phe311 (Fig. 2C
). Together, these residues form a narrow cavity where the acetate and the citrate are tightly bound. At the bottom of the cavity, the acetate molecule occupies the catalytic binding site of the h3
HSD3 enzyme with the C=O of its carboxyl group mimicking the ketone group of a steroid substrate. Indeed, one of the oxygen atoms of this carboxylate group is involved in a hydrogen bond network that includes the OH group of Tyr55 (2.7 Å), the N
atom of His117 (2.7 Å), in addition to the N
of Lys84 (located at 3.1 Å from the tyrosine hydroxyl group), while the carbon atom is located at 3.3 Å from the C-4 of the nicotinamide head of the cofactor, the hydride donor atom. The presence of an acetate molecule well stabilized inside the catalytic binding site of an enzyme of the AKR superfamily has been reported before (Nahoum et al. 2001; Lovering et al. 2004), but this is the first time it is seen orientated so as to interact with the catalytic residues (Tyr55 and His117) and, at the same time, participate in the ligand stabilization, here the citrate molecule. Indeed, one of the oxygen atoms of the acetate molecule carboxylate group is equidistant (2.5 Å) from two oxygen atoms (O2
and O7
) of the citrate molecule with which it could share a hydrogen atom (Fig. 2C
). In addition, two other oxygen atoms of the citrate molecule are involved in hydrogen bonds with residues of the h3
HSD3 enzyme: O2
with the side chain (N
2) of residue His222 (2.6 Å) and O5
with the hydroxyl group of Tyr24 (2.7 Å) and with the nitrogen atom of the indole ring of Trp227 (2.8 Å).
The numerous potential hydrogen bonds that stabilize the acetate and citrate in the substrate-binding cavity of h3
HSD3 led us to think that these two molecules could probably interfere in the binding of the steroid, which normally forms only one hydrogen bond with the enzyme, and thus could be potential inhibitors of h3
HSD3. We thus tested the capacity of the acetate and citrate to inhibit the 3
-reduction of DHT and found that none of them, alone or together, could measurably inhibit the h3
HSD3 activity (data not shown). This result suggests that weaker and less specific forces such as van der Waals bonds or hydrophobic interactions, which are well known to mediate the binding of steroids and other structurally different ligands by h3
HSD3, are probably much more important than anticipated, provided of course that the plasticity of the substrate pocket is sufficient to insure the steric complementarity between the enzyme and its ligands.
Citrate bound in a steroid-binding cavity suggests a ligand induced-fit mechanism
The side chains of residues Tyr55 and His117 form an anionic binding site that could bind molecules bearing a negative charge. This, combined with the fact that the substrate-binding site is funnel shaped and largely open on the protein surface, probably explains why human 3
HSD3 can bind and transform numerous molecules that vary in length and shape. Indeed, it has been shown that this enzyme can bind acetate (Nahoum et al. 2001), citrate (this paper), indanol (Matsuura et al. 1997), steroids (Dufort et al. 1996; Penning et al. 2000), and PAHs (Burczynski et al. 1998). To understand the mechanism underlying this plasticity we have superposed and compared our model (h3
HSD3/NADP(H)/Citrate/Acetate) with another h3
HSD3 complex structure (h3
HSD3/NADP(H)/Testosterone/Acetate). We first observed that a short segment of loop A, between residues Pro124 and Phe139, makes a movement against the citrate molecule that slightly closes the top of the binding cavity (Fig. 2A
). Loop A is one of the three loops involved in the formation of the testosterone-binding cavity in rat and human 3
HSD enzymes. It has been previously suggested for HSDs of the AKR superfamily that this loop plays a key role in the recognition, binding, and transformation of their substrates since a chimera rat 3
HSD having the loop A found in human 20
HSD sequence acquires 17
- and 20
HSD activities on testosterone and progesterone, respectively (Ma and Penning 1999). Among conformational changes observed in this segment, the backbone of residues Val128 and Ile129 moves 0.8 Å toward the citrate, contributing to a reduction of the cavity opening. These two residues are fairly important since they directly interact with and participate in the binding of testosterone and 20
-OHProg by h3
HSD3 and h20
HSD, respectively (Nahoum et al. 2001; Couture et al. 2003). In addition, previous studies on h20
HSD enzymes revealed that the side chain of another residue of loop A, Glu127, rotates upon binding of the steroid and slightly closes the top of the cavity, preserving the hydrophobic environment of the steroid-binding site (Couture et al. 2003). All these observations support the importance of the role played by loop A in the formation of the ligand-binding cavity and in the plasticity of this site. In comparison with loop A, the main chains of loops B and C do not undergo a significant structural shift (Fig. 2A
) but close inspection of the super-imposed binding site regions reveals that some residues that belong to these loops (Trp227, Leu306, and Leu308) or which delineate the cavity (Val54) make important conformational changes and become closer to the citrate molecule (Fig. 2D
).
Compared with the model of h3
HSD3 in ternary complex with testosterone/acetate, the side chain of residue Val54 in the present model undergoes a rotation of 180° around its C
, placing its C
atom toward the binding cavity (Fig. 2D
), probably to stabilize the citrate molecule. Interestingly, the importance of this residue for the selectivity of HSDs of the AKR superfamily has been previously demonstrated by site-directed mutagenesis experiments (Matsuura et al. 1997; Couture et al. 2003, 2004). For example, Matsuura et al. (1997) demonstrated that the replacement of residue Val54 by a leucine (residue found in the h20
HSD sequence) confers to the mutated h3
HSD3 almost the same catalytic properties as h20
HSD. Now, our structural results tend to demonstrate that residue Val54, through its capacity to adopt different conformations, is also important for the plasticity of the ligand-binding cavity of h3
HSD3.
We also observed that the indole ring of residue Trp227 moves toward the citrate molecule, where it is stabilized through a hydrogen bond between the O5
and the nitrogen atom of its indole ring (Fig. 2D
). In the ternary complex structure of h3
HSD3 with testosterone, this residue makes a
-stacking interaction with cycles B and C of testosterone and greatly contributes to the formation of a mature steroid-binding cavity (Nahoum et al. 2001). In r3
HSD, the replacement of Trp227 has been shown to affect the capacity of this enzyme to bind and transform testosterone and progesterone (Heredia et al. 2004) and, more importantly, its replacement can affect the activity on decalone and 9,10-phenanthrenequinone (Fig. 1B
). Combined with the h3
HSD3/NADP(H)/Acetate/Citrate ternary complex structure study, these results suggest that upon the binding of a non-steroid-shaped ligand, this residue is of great importance for the optimal formation of a mature cavity and for the plasticity of this enzyme.
Finally, the side chains of residues Leu306 and Leu308 also rotate 80° around their C
toward the citrate molecule (Fig. 2D
). This "leucine-flip" thus seems to play a role in the plasticity of this cavity. Similarly, a recent study on aldose reductase has pointed out the Leu300, a residue found in loop C (also called specificity pocket), as an important factor in the modulation of the shape of this cavity. It appears that the conformation adopted by the side chain of this residue with the backbone shift observed for residues Ala299 and Leu300 are directly related to the shape of the bound molecule (Sotriffer et al. 2004). Our finding indicating that the conformational change undergone by Leu306 and Leu308 reduced the size of the binding cavity of h3
HSD3 upon binding of the citrate is thus in total agreement with the model proposed for the aldose reductase, even if the contraction movement observed here does not affect the backbone position of these two residues.
To ascertain that the size of the h3
HSD3 substrate cavity really changes upon binding of its structurally different substrates, and thus confirm our first visual finding, we decided to precisely measure the volume of this surface groove in the presence of the acetate/citrate (monomer A) and compare it with that of the same enzyme in complex with testosterone/acetate (Nahoum et al. 2001) (RCSB PDB entry, 1J96
[PDB]
) and ursodeoxycholate (Jin et al. 2001) (RCSB PDB entry, 1IHI
[PDB]
). We thus found that the substrate cavity volume was directly related to the size of the bound ligand. Indeed, the cavity volume for the acetate/citrate (Fig. 3A
) was 704 Å3 compared to 776 Å3 for testosterone/acetate (Fig. 3B
), in agreement with the size of these molecules (236 Å3 vs. 333 Å3). Interestingly, the cavity volume in the ursodeoxycholate complex (734 Å3) is smaller than the testosterone cavity, although ursodeoxycholate (362 Å3) is slightly bigger (Fig. 3C
). This apparent contradiction could be easily explained by the difference in the binding mode adopted by these two steroidal substrates. The ursodeoxycholate is bound in such a manner that its long and thin 17
-(1-methyl-3- carboxypropyl) arm is positioned close to the catalytic site of the enzyme, pushing almost half of the steroidal nucleus (cycle A and a large part of cycle B) out of the cavity (Fig. 3D
). In this position, the ursodeoxycholate is in fact a "smaller" ligand than testosterone (248 Å3 vs. 277 Å3), which is almost totally inside the cavity.
|
HSD3 and probably to many other members of the AKR1C enzymes subfamily.
What is the real productive ternary complex?
Residues found in loops A, B, and C are of great importance in the plasticity of the steroid-binding cavity of h3
HSD3, a characteristic that allows this enzyme to bind molecules with various shapes and structures. In addition, these loops must optimally orient the incoming molecules to promote their transformation. Unfortunately, each HSD member of the AKR superfamily for which the ternary complex structure has been solved shows its steroid substrate in a nonproductive position, i.e., in a position (or in an orientation) that does not permit the oxido-reduction reaction to take place (Bennett et al. 1997; Nahoum et al. 2001; Couture et al. 2004). According to recent reports, HSDs of the AKR superfamily would transform steroids through an orbital steering mechanism where the optimal distance and angle between the carbon atom of the substrate ketone group and C4 of the nicotinamide would be around 3.0 Å and 100°, respectively (Heredia et al. 2003). Close inspection of the h3
HSD3 structure with the citrate and acetate molecules bound in the substrate cavity (monomer A) reveals that this optimal position corresponds exactly to the carbon atom of the carboxylate group of the acetate molecule.
To better understand how h3
HSD3 can bind androgen in a productive complex, we manually modeled steroid molecules (4-dione and DHT) in a position that respects the orbital steering constraints. When these two structurally similar steroids were modeled with their reactive group in an optimal position (the C3 for the DHT and C17 for the 4-dione), their steroid nuclei were found in different positions (Fig. 4
). For DHT, the nucleus of the steroid was roughly directed toward loop C and clashed with the indole ring of Trp227, a residue belonging to loop B and that is highly conserved among AKR members. Interestingly, we had previously found that Trp227, with many other residues of loop B, can undergo significant conformational changes to promote the binding of the steroid in the crystal structure of rb20
HSD in ternary complex with testosterone (Couture et al. 2004). In fact, we have found that Trp227 with some other residues of loop B that point toward the steroid-binding cavity in the rb20
HSD binary complex, probably thus maintaining the hydrophobic character of the steroid cavity, are at the surface of the protein and exposed to the solvent upon binding of testosterone. This relaxation movement results in an enlargement of the steroid-binding cavity and in the mobilization of some residues belonging to the C-terminal part of the protein (loop C) to create the mature steroid-binding site (Couture et al. 2004). Crystal structure determination of r3
HSD in ternary complex with testosterone also reveals that the indole ring of Trp227 can rotate by 180° around its C
, a movement that drastically changes the shape of the cavity and allows DHT to bind in a productive complex.
|
-reduction, the body of the steroid occupied a space delineated by the side chain of residues Leu306 and Leu308 with its C3-ketone group lying down near the side chain of residue Arg304. When modeled in a productive ternary complex, DHT and 4-dione do not occupy the same position in the binding cavity but both steroids interact with residues of the C-terminal end of the protein (Fig. 4
HSD greatly impairs the activity of this enzyme toward DHT and 4-dione (Couture et al. 2004). Likewise, the exchange of loop C sequence between h3
HSD1 and h20
HSD gives to the chimerical enzyme 20
HSD3
HSD1 (h20
HSD with 39 C-terminal residues of h3
HSD type 1) new properties and transforms it into a potent reductase enzyme for C19-steroids (Matsuura et al. 1998). Together, all these observations strongly support the hypothesis that the C-terminal end of these enzymes is of paramount importance in the discrimination of the incoming steroid substrates.
A C-terminal basic motif, a concerted way to bind androgen?
We thus found that, in order to respect the orbital steering constraints, the steroid should be maintained approximately parallel to the pyridine head of the cofactor with the body of the steroid pointing toward the C-terminal part of the protein. To evaluate the role of residues of this domain in steroid discrimination and binding, we undertook a systematic leucine scanning mutagenesis study. We found that the replacement of two residues, Arg301 and Arg304, leads to a drastic effect on the capacity of h3
HSD3 to transform androgen. Indeed, when the kinetic constants for the reduction of DHT were compared with the wild-type h3
HSD3 (Table 2
), we found that Km for the 3
-reduction of DHT was barely affected, whereas kcat was 60- and 227-fold lower for the Arg301Leu and Arg304Leu mutants, respectively. For the 17
-reduction of 4-dione, we found that the replacement of Arg301 affects the Km and the kcat values by 3.8- and 7.3-fold, respectively. More strikingly, we found that the 17
-activity on 4-dione was totally abolished for the Arg304Leu mutant enzyme. No activity was detectable even by using radio-activity, the most sensitive technique to detect the transformation of a steroid (Luu-The et al. 1995).
|
HSD3, the Kd values for NADPH were almost identical for the Arg301Leu and Arg304Leu mutant enzymes, demonstrating that the replacement of this basic motif does not impair the binding of the cofactor (Table 3
-reduction activity of h3
HSD3, whereas its role in the 3
-reduction of DHT is less important. More importantly, they also demonstrate that the residues involved in the formation of a mature cavity for 4-dione are not necessarily the same as those having the same role for a structurally similar C19-steroid, DHT.
|
|
HSD severely affected the Km value of this enzyme for the 20
-reduction of progesterone (Couture et al. 2003). Moreover, site-directed mutagenesis studies on human aldehyde reductase, another representative member of the AKR superfamily, have demonstrated that Arg311 is likely involved in the binding of aldehydes containing a carboxyl group and being longer than 78 Å (Barski et al. 1996). A superposition of the C-terminal loop sequence for HSDs of the AKR superfamily, the less conserved region among the AKR1C subfamily, reveals that Arg301 and Arg304 are conserved in 10 HSDs out of 12 (Fig. 5
HSD3 possesses these structural characteristics since a basic motif, the Arg301 and Arg304 residues acting as the proton donors, is at one end of the steroid molecule, whereas the catalytic Tyr55, the proton acceptor/donor residue, is at the other end.
|
HSD3 enzyme suggest that, for triggering of the reaction, 4-dione is likely maintained parallel to the pyridine head of the cofactor in order to optimally orient its C17-ketone group, with respect to the 4-proR hydride of the cofactor. In this position, the nucleus of the steroid lies down toward the C loop of the enzyme where the guanido group of Arg304 interacts with the steroid, either directly by forming an ionic bond with the C3-ketone group or indirectly by the formation of a positively charged environment essential for perfect positioning of the steroid substrate. It thus appears that, together, these structural characteristics are essential for the formation of the mature steroid-binding cavity and are the basis of the specificity of this enzyme for a steroid, a subgroup of steroids as well as other molecules.
These observations also have important physiological relevance since, in the prostate tissues, h3
HSD3 is able to inactivate the most potent androgen, DHT, into a weak androgen, 3
-diol, while this enzyme is also involved in the formation of testosterone, a potent androgen, through its 17
-reduction activity on 4-dione. If the residues involved in the binding of DHT and 4-dione differ, it is tempting to speculate that it would be feasible to synthesize an inhibitor specific for the 17
-reduction of 4-dione that does not affect, or very weakly affects, the 3
-inactivation of DHT. In this regard, it will be essential to take into account that some residues of the h3
HSD3 binding cavity involved in the binding and maintaining of the steroid substrates can undergo significant conformational changes, in part dictated by the shape of the steroid itself. Similar to the "ligand induced-fit mechanism" reported for the aldose reductase, this binding site adaptation is of major importance in the efficiency of inhibitors targeted against this enzyme.
| Materials and methods |
|---|
|
|
|---|
Site-directed mutagenesis
Arg301Leu and Arg304Leu mutations were created using Quick- Change Site-Directed Mutagenesis Kit (Stratagene) with these forward primers: Arg301Leu, AAAGCCATAGATGGCC TAAACCTAAATGTGCGATATTTGACCCTTGATAT and Arg304Leu, CATAGATGGCCTAAACAGAAATGTGCTAT ATTTGACCCTTGATATTTTTGATGGCC. To confirm the presence of the desired mutation and to ascertain that no other mutations had occurred, the complete coding region of the mutated cDNA was sequenced with a T7 sequencing kit (Amersham Bioscience).
Overexpression and purification of wild-type and mutant h3
HSD3
Briefly, the cDNA encoding h3
HSD3 was subcloned in a pGEX vector (Amersham Bioscience), overproduced in E. coli BL21(DE3)pLysS as a fusion protein with glutathione-S transferase (GST) and purified by a combination of affinity, anionic, and size-exclusion chromatographies, yielding about 5 mg of a highly purified enzyme preparation per 100 mL of cell culture. Mutated h3
HSD3s were overproduced and purified to homogeneity as described for the wild-type enzyme.
Determination of kinetic parameters
Reduction reactions on [14C]-steroids were performed as previously described, with minor modifications (Luu-The et al. 1995; Couture et al. 2002). Kinetic parameters (Km and kcat/Km) were calculated using the EnzFitter program (Biosoft). Because some HSDs of the AKR superfamily have been shown to lose their activity when they are purified to homogeneity (Dufort et al. 1999), a kinetic characterization was performed with DHT and 4-dione, two natural ligands of this enzyme, to verify that the purification steps did not affect the enzyme activity. Within the experimental error, Km values obtained for the 3
-reduction of the DHT (5.7 ± 0.8 µM) and for the 17
-reduction of 4-dione (10.5 ± 1.8 µM) (Table 2
) were very similar to those published earlier (1.4 ± 0.66 µM and 9.73 ± 1.8 µM [Penning et al. 2000; Dufort et al. 2001]), indicating that our h3
HSD3 preparation was suitable for subsequent experiments.
Reduction of 1-indanone was measured at room temperature in a 50-µL system containing 100 µM 1-indanone in 2% (v/v) ethanol, 200 µM NADPH, and 100 mM KH2PO4K2HPO4 (pH 7.5). Reactions were started when 10 µg of homogenous recombinant type 3 h3
HSDs were added to the system in which 1-indanone transformation was monitored with a Beckman DU-7400 spectrophotometer by measuring the rate of change of absorbance of NADPH (
= 6270 cm1 M1) in a 1-cm light path.
Fluorescence
Measurements of the fluorescence of wild-type, Arg301Leu, and Arg304Leu enzymes were conducted on an SLM 8000 fluorometer. Scans of the fluorescence emission spectra were performed at 25°C in a 200-µL microcell containing 2 µM of the purified enzyme, 50 mM Na2HPO4/NaH2PO4 buffer (pH 7.5), 20% glycerol, 1 mM EDTA, and variable concentrations of NADPH (04000 nM). Tryptophan fluorescence was detected using 290 nm as the excitation wavelength, while the emission was scanned from 300 to 600 nm (at 120 nm/min). Emission and excitation bandpasses were both set to 4 nm.
Crystallization
A homogenous preparation of h3
HSD3 enzyme freshly purified was used for crystallization experiments using the hanging-drop vapor diffusion method at 4°C. Crystals were obtained by mixing equal volumes of h3
HSD3 concentrated to 14 mg/mL and of mother liquor (0.1 M NaCitrate at pH 6.5, 0.1 M NH4Ac, and 24%30% PEG4000). Crystal used for X-ray diffraction experiments was directly soaked into the cryopreserving solution (20% of ethylene glycol added to the mother liquor) and flash-cooled at 100 K in a nitrogen gas stream.
Structure determination and refinement
X-ray diffraction oscillation images of h3
HSD3/NADPH binary complex were recorded on an R-AXIS IIc image-plate detector mounted on a Rigaku RU-200 rotating-anode generator equipped with focusing mirrors (Rigaku MSC). Data were integrated and scaled using the December 2000 version of the XDS program (Kabsch 1993) and the crystal structure was directly determined by performing a rigid body refinement using the Refmac program (Murshudov et al. 1997). The h3
HSD3 structure (PDB code 1J96
[PDB]
; Nahoum et al. 2001) was used as a model and molecular replacement solution was calculated with two molecules per asymmetric unit. Crystallographic statistics are shown in Table 1
. The refinement procedure was performed using CNS (Brunger et al. 1998) and Refmac (Murshudov et al. 1997) programs. The initial model issued from rigid body refinement was submitted to a cycle of simulated annealing at 3000 K followed by energy minimization cycles. After electronic density maps calculations and manual rebuilding using O (Jones et al. 1991), NADPH was added to the model. Atomic models were then refined by simple energy minimization followed by isotropic B-factors refinement (restrained, individual B-factor refinement) and corrected by manual rebuilding. Missing parts of the model, sulfate ion,
-mercaptoethanol, ethylene glycol, citric acid, acetate, and water molecules were progressively added during the refinement procedure. Finally, the quality of the model was verified with PROCHECK (Table 1
) (Laskowski et al. 1993).
The substrate cavity volumes of h3
HSD3 in the presence of its different ligands were measured with the program VOIDOO (Kleywegt 1994) using program default parameters. The volume values obtained this way were finally corrected using the equation given by Chakravarty et al. (2002).
Accession numbers
Atomic coordinates of this binary complex have been deposited in the Protein Data Bank with the accession code 1XJB.
| Footnotes |
|---|
2 Synchrotron Soleil, LOrme des Merisiers, Saint-Aubin BP 48, 91192 Gif-sur-Yvette CEDEX, France. ![]()
| Acknowledgments |
|---|
| References |
|---|
|
|
|---|
Bennett, M.J., Schlegel, B.P., Jez, J.M., Penning, T.M., and Lewis, M. 1996. Structure of 3
-hydroxysteroid/dihydrodiol dehydrogenase complexed with NADP+. Biochemistry 35: 1070210711.[CrossRef][Medline]
Bennett, M.J., Albert, R.H., Jez, J.M., Ma, H., Penning, T.M., and Lewis, M. 1997. Steroid recognition and regulation of hormone action: Crystal structure of testosterone and NADP+ bound to 3
-hydroxysteroid/dihydrodiol dehydrogenase. Structure 5: 799812.[Medline]
Brunger, A.T., Adams, P.D., Clore, G.M., Delano, W.L., Gros, P., Grosse-Kunstleve, R.W., Jiang, J.-S., Kuszewski, J., Nilges, M., Pannu, N.S., et al. 1998. Crystallography & NMR system: A new software system for macromolecular structure determination. Acta Crystallogr. D Biol. Crystallogr. 54: 905921.[CrossRef][Medline]
Burczynski, M.E., Harvey, R.G., and Penning, T.M. 1998. Expression and characterization of four recombinant human dihydrodiol dehydrogenase isoforms: Oxidation of trans-7, 8-dihydroxy-7,8-dihydrobenzo[a]pyrene to the activated o-quinone metabolite benzo[a]pyrene-7,8-dione. Biochemistry 37: 67816790.[CrossRef][Medline]
Burczynski, M.E., Sridhar, G.R., Palackal, N.T., and Penning, T.M. 2001. The reactive oxygen species- and Michael acceptor-inducible human aldo-keto reductase AKR1C1 reduces the
,
-unsaturated aldehyde 4-hydroxy- 2-nonenal to 1,4-dihydroxy-2-nonene. J. Biol. Chem. 276: 28902897.
Chakravarty, S., Bhinge, A., and Varadarajan, R. 2002. A procedure for detection and quantitation of cavity volumes proteins. Application to measure the strength of the hydrophobic driving force in protein folding. J. Biol. Chem. 277: 3134531353.
Couture, J.F., Cantin, L., Legrand, P., Luu-The, V., Labrie, F., and Breton, R. 2002. Expression, crystallization and preliminary X-ray analysis of human and rabbit 20
-hydroxysteroid dehydrogenases in complex with NADP(H) and various steroid substrates. Acta Crystallogr. D Biol. Crystallogr. 58: 135139.[CrossRef][Medline]
Couture, J.F., Legrand, P., Cantin, L., Luu-The, V., Labrie, F., and Breton, R. 2003. Human 20
-hydroxysteroid dehydrogenase: Crystallographic and site-directed mutagenesis studies lead to the identification of an alternative binding site for C21-steroids. J. Mol. Biol. 331: 593604.[CrossRef][Medline]
Couture, J.F., Legrand, P., Cantin, L., Labrie, F., Luu-The, V., and Breton, R. 2004. Loop relaxation, a mechanism that explains the reduced specificity of rabbit 20
-hydroxysteroid dehydrogenase, a member of the aldo-keto reductase superfamily. J. Mol. Biol. 339: 89102.[CrossRef][Medline]
DelaFuente, G. and Sols, A. 1970. The kinetics of yeast hexokinase in the light of the induced fit involved in the binding of its sugar substrate. Eur. J. Biochem. 16: 234239.[Medline]
DelaFuente, G., Lagunas, R., and Sols, A. 1970. Induced fit in yeast hexokinase. Eur. J. Biochem. 16: 226233.[Medline]
Dufort, I., Soucy, P., Labrie, F., and Luu-The, V. 1996. Molecular cloning of human type 3 3
-hydroxysteroid dehydrogenase that differs from 20
-hydroxysteroid dehydrogenase by seven amino acids. Biochem. Biophys. Res. Commun. 228: 474479.[CrossRef][Medline]
Dufort, I., Rheault, P., Huang, X.F., Soucy, P., and Luu-The, V. 1999. Characteristics of a highly labile human type 5 17
-hydroxysteroid dehydrogenase. Endocrinology 140: 568574.
Dufort, I., Labrie, F., and Luu-The, V. 2001. Human types 1 and 3 3
-hydroxysteroid dehydrogenases: Differential lability and tissue distribution. J. Clin. Endocrinol. Metab. 86: 841846.
Fenn, T.D., Ringe, D., and Petsko, G.A. 2003. POVScript+: A program for model and data visualization using persistence of vision ray-tracing. J. Appl. Crystallogr. 36: 944947.[CrossRef]
Heredia, V.V., Kruger, R.G., and Penning, T.M. 2003. Steroid-binding site residues dictate optimal substrate positioning in rat 3
-hydroxysteroid dehydrogenase (3
HSD or AKR1C9). Chem. Biol. Interact. 143144: 393400.
Heredia, V.V., Cooper, W.C., Kruger, R.G., Jin, Y., and Penning, T.M. 2004. Alanine scanning mutagenesis of the testosterone binding site of rat 3
-hydroxysteroid dehydrogenase demonstrates contact residues influence the rate-determining step. Biochemistry 43: 58325841.[CrossRef][Medline]
Hoog, S.S., Pawlowski, J.E., Alzari, P.M., Penning, T.M., and Lewis, M. 1994. Three-dimensional structure of rat liver 3
-hydroxysteroid/dihydrodiol dehydrogenase: A member of the aldo-keto reductase superfamily. Proc. Natl. Acad. Sci. 91: 25172521.
Jin, Y., Stayrook, S.E., Albert, R.H., Palackal, N.T., Penning, T.M., and Lewis, M. 2001. Crystal structure of human type III 3
-hydroxysteroid dehydrogenase/bile acid binding protein complexed with NADP(+) and ursodeoxycholate. Biochemistry 40: 1016110168.[CrossRef][Medline]
Jones, T.A., Zou, J.Y., Cowan, S.W., and Kjeldgaard, M. 1991. Improved methods for the building of protein models in electron density maps and the location of errors in these models. Acta Crystallogr. A 47: 110119.
Kabsch, W. 1993. Automatic processing of rotation diffraction data from crystals of initially unknown symmetry and cell constants. J. Appl. Crystallogr. 26: 795800.[CrossRef]
Kaplan, W. and Littlejohn, T.G. 2001. Swiss-PDB Viewer (Deep View). Brief Bioinform. 2: 195197.
Kavanagh, K.L., Klimacek, M., Nidetzky, B., and Wilson, D.K. 2002. The structure of apo and holo forms of xylose reductase, a dimeric aldo-keto reductase from Candida tenuis. Biochemistry 41: 87858795.[CrossRef][Medline]
Khurana, S., Powers, D.B., Anderson, S., and Blaber, M. 1998. Crystal structure of 2,5-diketo-D-gluconic acid reductase A complexed with NADPH at 2.1-A resolution. Proc. Natl. Acad. Sci. 95: 67686773.
Klebe, G., Kramer, O., and Sotriffer, C. 2004. Strategies for the design of inhibitors of aldose reductase, an enzyme showing pronounced induced-fit adaptations. Cell. Mol. Life Sci. 61: 783793.[CrossRef][Medline]
Kleywegt, G.J. 1994. Detection, delineation, measurement and display of cavities in macromolecular structures. Acta Crystallogr. D Biol. Crystallogr. 50: 178185.[CrossRef][Medline]
Kozma, E., Brown, E., Ellis, E.M., and Lapthorn, A.J. 2002. The crystal structure of rat liver AKR7A1. A dimeric member of the aldo-keto reductase superfamily. J. Biol. Chem. 277: 1628516293.
Kraulis, P.J. 1991. MOLSCRIPT: A program to produce both detailed and schematic plots of protein structures. J. Appl. Crystallogr. 24: 946950.[CrossRef]
Laskowski, R.A., MacArthur, M.W., Moss, D.S., and Thornton, J.M. 1993. PROCHECK: A program to check the stereochemical quality of protein structures. J. Appl. Crystallogr. 26: 283291.[CrossRef]
Lovering, A.L., Ride, J.P., Bunce, C.M., Desmond, J.C., Cummings, S.M., and White, S.A. 2004. Crystal structures of prostaglandin D(2) 11-ketoreductase (AKR1C3) in complex with the nonsteroidal anti-inflammatory drugs flufenamic acid and indomethacin. Cancer Res. 64: 18021810.
Luu-The, V., Zhang, Y., Poirier, D., and Labrie, F. 1995. Characteristics of humantypes1, 2 and 317
-hydroxysteroid dehydrogenase activities: Oxidation/reduction and inhibition. J. Steroid Biochem. Mol. Biol. 55: 581587.[CrossRef][Medline]
Ma, H. and Penning, T.M. 1999. Conversion of mammalian 3
-hydroxysteroid dehydrogenase to 20
-hydroxysteroid dehydrogenase using loop chimeras: Changing specificity from androgens to progestins. Proc. Natl. Acad. Sci. 96: 1116111166.
Matsuura, K., Deyashiki, Y., Sato, K., Ishida, N., Miwa, G., and Hara, A. 1997. Identification of amino acid residues responsible for differences in substrate specificity and inhibitor sensitivity between two human liver dihydrodiol dehydrogenase isoenzymes by site-directed mutagenesis. Biochem. J. 323 (Pt. 1): 6164.
Matsuura, K., Hara, A., Deyashiki, Y., Iwasa, H., Kume, T., Ishikura, S., Shiraishi, H., and Katagiri, Y. 1998. Roles of the C-terminal domains of human dihydrodiol dehydrogenase isoforms in the binding of substrates and modulators: Probing with chimaeric enzymes. Biochem. J. 336 (Pt 2): 429436.
Mesecar, A.D., Stoddard, B.L., and Koshland Jr., D.E. 1997. Orbital steering in the catalytic power of enzymes: Small structural changes with large catalytic consequences. Science 277: 202206.
Murshudov, G.N., Vagin, A.A., and Dodson, E.J. 1997. Refinement of macromolecular structures by the maximum-likelihood method. Acta Crystallogr. D 53: 240255.[CrossRef][Medline]
Nahoum, V., Gangloff, A., Legrand, P., Zhu, D.W., Cantin, L., Zhorov, B.S., Luu-The, V., Labrie, F., Breton, R., and Lin, S.X. 2001. Structure of the human 3
-hydroxysteroid dehydrogenase type 3 in complex with testosterone and NADP at 1.25-A resolution. J. Biol. Chem. 276: 4209142098.
Nahoum, V., Gangloff, A., Shi, R., and Lin, S.X. 2003. How estrogen-specific proteins discriminate estrogens from androgens: A common steroid binding site architecture. FASEB J. 17: 13341336.
Palackal, N.T., Burczynski, M.E., Harvey, R.G., and Penning, T.M. 2001. The ubiquitous aldehyde reductase (AKR1A1) oxidizes proximate carcinogen trans-dihydrodiols to o-quinones: Potential role in polycyclic aromatic hydrocarbon activation. Biochemistry 40: 1090110910.[CrossRef][Medline]