Discussion
The first structure of P450 21A2 was that of a slightly modified bovine enzyme (T241R/L442A) with the substrate 17α-OH-progesterone (
8- Zhao B.
- Lei L.
- Kagawa N.
- Sundaramoorthy M.
- Banerjee S.
- Nagy L.D.
- Guengerich F.P.
- Waterman M.R.
Three-dimensional structure of steroid 21-hydroxylase (cytochrome P450 21A2) with two substrates reveals locations of disease-associated variants.
), and two molecules of this substrate were present, one bound at the active site and the other at a distal site adjacent to the F′ helix. We have not obtained a structure of the bovine enzyme with the substrate progesterone bound. The structures of human P450 21A2 bound with progesterone (
10- Pallan P.S.
- Wang C.
- Lei L.
- Yoshimoto F.K.
- Auchus R.J.
- Waterman M.R.
- Guengerich F.P.
- Egli M.
Human cytochrome P450 21A2, the major steroid 21-hydroxylase: structure of the enzyme-progesterone substrate complex and rate-limiting C–H bond cleavage.
) also revealed a substrate molecule at a distal site that matches the site of the second 17α-OH-progesterone in the structure of the complex with bovine 21A2. Most likely due to the lower resolution of the structure of the human 21A2 enzyme with 17α-OH-progesterone (this work; see
Fig. 1), the electron density at this distal site was not defined well enough to identify a second substrate molecule.
The availability of structural, activity, stability, and substrate binding data offers an opportunity to potentially gain a refined understanding of the phenotypes of diverse CAH variants and the different degrees to which SW, SV, and NC mutations affect the P450 21A2 enzyme. Among the variants we investigated in this work are rare mutations as well as common pseudogene-derived ones. The latter include I172N (SV), P31L (NC), and V282L (the most common NC variant).
6R. J. Auchus, personal communication.
Ile-172 (I172N variant, SV) is located in the middle of helix E, which is part of a helical bundle (C, E, I, and M helices;
Fig. 1A). Ile-172 is part of a relatively hydrophobic environment including Val-140, Leu-176, and Leu-434 and also a hydrophilic residue, Glu-438 (
Fig. 7A). Replacing isoleucine with asparagine affects these hydrophobic interactions to some extent but also allows hydrogen bonds between the amino group of Asn-172 and Glu-438 (Oϵ2; 3.06 Å) and between the side chain carbonyl of Asn-172 and the main-chain N–H of Val-140 (3.60 Å). I172N displayed a considerable loss in catalytic activity (∼2% residual activity with 17α-OH-progesterone;
Table 2) and a 7-fold reduction in heme incorporation compared with wild-type P450 21A2 (
Table 3).
Val-282 (V282L variant, NC) is located near the N-terminal end of the long I helix and surrounded by hydrophobic residues from the G (Leu-242) and H helices (Met-258, Met-261, and Leu-262) (
Fig. 7B). The increased size of the leucine relative to the valine side chain can be accommodated to some degree at this site, consistent with the residual activity and the NC phenotype (
Table 2). The activity we measured here is considerably lower than that reported earlier (20% for progesterone and 50% for 17α-OH-progesterone;
http://www.cypalleles.ki.se/cyp21.htm)
7Please note that the JBC is not responsible for the long-term archiving and maintenance of this site or any other third party hosted site.
(
22- Tusie-Luna M.T.
- Traktman P.
- White P.C.
Determination of functional effects of mutations in the steroid 21-hydroxylase gene (CYP21) using recombinant vaccinia virus.
). Among the variants for which we evaluated the stability change of the Fe
2+-CO spectrum based on heating, the V282L mutant showed the greatest reduction in
Tm relative to the wild-type protein (
Fig. 6E). This variant also displayed a significantly distorted CD spectrum even at 25 °C (
Fig. 5D), supporting the notion that even the addition of a single methyl group of the leucine side chain relative to valine is not well-tolerated. This variant also showed very poor affinity for substrates (
Table 4).
Gly-292 (G292C variant, SW) maps to the central region of the I helix and lies in close proximity of substrate and heme (
Fig. 7C). The cysteine mutation results in steric clashes with ring A of 17α-OH-progesterone and progesterone and potentially heme. The CD spectrum was very distorted even at 25 °C (
Fig. 5B). The complete loss of activity in this severe mutant is therefore not surprising (
Table 2), despite a heme incorporation that was only marginally affected (75% of wild-type enzyme;
Table 3). The effect of temperature was not great, surprisingly (
Figure 5,
Figure 6 and
supplemental Fig. S1). Binding of progesterone was very poor (
Table 4), but 17α-OH-progesterone bound reasonably well. Interestingly, progesterone binding (relative
KS) is dramatically reduced compared with wild-type 21A2 (> 10
4-fold;
Table 4), although the effect is less prominent in the case of 17α-OH-progesterone.
The G292S variant (SW) is less damaging than the cysteine mutant as indicated by activity data that show ∼100-fold reduction compared with wild-type P450 21A2 (
Table 2). The serine side chain is slightly smaller, and the γ-hydroxyl may be more easily accommodated, given its vicinity to more polar moieties, such as the 17α-OH group and the main-chain carbonyl groups of Asp-288 and Leu-289 (
Fig. 7D). Indeed, this variant binds progesterone about 25-fold more tightly than the P450 G292C variant (
Table 4). Neither the G292C nor the G292S variant exhibits a dramatically altered thermal stability. In fact, in the Fe
2+-CO spectra, the
Tm of the latter was only 2 °C below that of wild-type P450 21A2 (
Fig. 6D). However, just like the V282L variant, G292C showed a significantly changed CD spectrum relative to wild-type protein even at 25 °C (
Fig. 5B).
The Arg-357 (R357W) variant (SW) is located at the C terminus of the K helix. Methylene groups of the arginine side chain are surrounded by relatively hydrophobic side chains, including Trp-303, Pro-464, and Leu-472, whereas the guanidino moiety faces Asp-44, Gln-390, and Gln-463. The mutation to tryptophan results in a stacking interaction between the side chains of Trp-357 and Trp-303 (∼3.8-Å separation;
Fig. 7E). The variant displayed some residual catalytic activity (< 1% with progesterone or 17α-OH-progesterone;
Table 2) (but displayed the lowest level of heme incorporation among the investigated CAH variants, < 10% relative to wild-type P450 21A2;
Table 3).
Arg-409 (R409C variant, SW) is located in a highly polar environment in the loop region that connects helices L and M and is part of the highly conserved E
XXR motif (see
Fig. 4D in Ref.
11- Pallan P.S.
- Lei L.
- Wang C.
- Waterman M.R.
- Guengerich F.P.
- Egli M.
Research resource: correlating human cytochrome P450 21A2 crystal structure and phenotypes of mutations in congenital adrenal hyperplasia.
). The mutation to cysteine introduces a hydrophobic and less sterically demanding entity into this polar environment (
Fig. 7F). Thus, the cysteine side chain is unlikely to create any short contacts, and it appears that it can be accommodated between the side chains of Trp-406 and Asn-416. Given the intricate interactions that Arg-409 is engaged in (hydrogen-bonding with four neighboring residues and involving both side and main chain atoms of these;
Fig. 7F), none of which can be mimicked by cysteine, it is perhaps surprising that the R409C CAH variant still displays between ∼0.5 and 2.5% of the wild-type activity (
Table 2).
We had evaluated potential steric and electrostatic effects of 24 SW, SV, and NC CAH variants based on the previously determined crystal structure of the complex between human P450 21A2 and progesterone (
10- Pallan P.S.
- Wang C.
- Lei L.
- Yoshimoto F.K.
- Auchus R.J.
- Waterman M.R.
- Guengerich F.P.
- Egli M.
Human cytochrome P450 21A2, the major steroid 21-hydroxylase: structure of the enzyme-progesterone substrate complex and rate-limiting C–H bond cleavage.
,
11- Pallan P.S.
- Lei L.
- Wang C.
- Waterman M.R.
- Guengerich F.P.
- Egli M.
Research resource: correlating human cytochrome P450 21A2 crystal structure and phenotypes of mutations in congenital adrenal hyperplasia.
). Among them were the P31L (NC), P31Q (SW), G65E (SW), L108R (SW), T296N (SW), and W303R (SW) variants, for which we now report catalytic activities with the progesterone and 17α-OH-progesterone substrates (
Table 2). Consistent with the severity of the CAH phenotype, the P31L mutant does not hamper catalysis to the same extent as the P31Q mutant (the latter is basically inactive). The deviating magnitude of the activity reductions seen for the two variants are matched by those on heme incorporation (3-fold lower for P31Q;
Table 3), and the loss in stability at least for P31L relative to the wild-type protein is quite minor (
Fig. 6B). Pro-31, along with an adjacent proline, marks the site of a sharp turn in the direction of the protein chain and helps to anchor the N-terminal hydrophobic tail that serves as the attachment of P450 21A2 to the membrane. The leucine mutation is sterically compromising but less severe, from a structural perspective, than the glutamine mutation that affects both sterics and electrostatics/hydrogen bonding in comparison with proline.
Gly-65 (G65E variant, SW) is close to the surface and at the apex of the short loop that connects the β2 and β3 strands, facing the N-terminal residue of the F′ helix (Pro-214;
Fig. 1A). Mutation to glutamate introduces a polar residue into a patch of hydrophobic side chains from residues Leu-38, Leu-40, Leu-64, Leu-66, Val-212, and Ile-213. In addition to disturbing the apolar character of this location, the relatively long glutamate side chain creates steric challenges, and the significant loss of catalytic activity with this variant (∼0.1% of wild-type protein with both substrates;
Table 2) is in line with the anticipated structural consequences. The poor incorporation of heme (25% of wild-type 21A2;
Table 3) and the significantly reduced binding, particularly of progesterone (some 2,000-fold relative to wild type;
Table 4), are perhaps more surprising, given the location of Gly-65 close to the surface and at a considerable distance from active site and heme.
Leu-108 (L108R variant, SW) is located in a loop that precedes helix C (
Fig. 1A) and is positioned immediately adjacent to methyl and propionic acid substituents of heme. Introduction of the longer arginine with its positively charged guanidino moiety in place of leucine that sits in a partly hydrophobic region but also has His-120, Lys-121, and Asp-288 among its neighbors was expected to distort the conformational/electronic balance. Consistent with this assessment, the CD spectrum for the L108R variant displays significant changes relative to wild-type protein even at 25 °C (
Fig. 5H). Moreover, the complete loss of activity (
Table 2) is in line with the structural picture and the severity of the CAH phenotype. The activities based on
in vitro kinetic assays and using purified mutant protein reported here are lower than those reported earlier and based on experiments conducted with transiently transfected COS-1 cells (∼0.3% activity with both substrates relative to wild-type protein (
23- Soardi F.C.
- Barbaro M.
- Lau I.F.
- Lemos-Marini S.H.
- Baptista M.T.
- Guerra-Junior G.
- Wedell A.
- Lajic S.
- de Mello M.P.
Inhibition of CYP21A2 enzyme activity caused by novel missense mutations identified in Brazilian and Scandinavian patients.
)). Perhaps surprisingly, given the close vicinity of this mutation to heme, the fraction of heme measured relative to wild-type 21A2 still amounts to 40% (
Table 3).
In comparison with the above L108R variant, mutation of Thr-296 to asparagine, although classified as an SW CAH phenotype as well, affects the catalytic activity to a smaller degree: 0.05 and 4% of wild-type 21A2 activity with the progesterone and 17α-OH-progesterone substrates, respectively (
Table 2). Heme incorporation was also increased relative to the L108R variant (62 and 25%, respectively;
Table 3), and the thermal stability was virtually the same as for wild-type protein (
Fig. 6F; −1 °C). Thr-296 is located in the long I helix and sits above the side of the heme (5.2-Å distance between Oγ and Fe
3+) and in the vicinity of ring A of the substrate molecule (5.5-Å distance between Oγ and O17) (
Fig. 1D). In the modeled Asn-296 mutant, the above distances would be reduced to 4.7 and 3.1 Å (Nδ2…Fe
3+ and Nδ2…O17, respectively). Apparently, swapping threonine and asparagine does not entirely abolish activity (some 25-fold reduced compared with wild-type protein with 17α-OH-progesterone), an observation that is noteworthy, given how closely it occurs to substrate and heme and given the more drastic consequences of mutations farther away from the active site (
e.g. P31Q and L65R).
The consequences of the T296N mutation are also more limited than those seen for another SW CAH variant, W303R, Trp-303 being located farther along the I helix and thus more removed from heme and substrate than Thr-296. Mutation of Trp-303 to arginine destroys the catalytic activity with both substrates (
Table 2) and also affects heme incorporation to a very significant extent (11%;
Table 3). An inspection of the structures of the human P450 21A2–substrate complexes suggests that an arginine at position 303 is likely to create a repulsion with Arg-357 that sits at the end of the K helix (
Fig. 1A). In the model of the mutant, the guanidino moieties of the two arginines are positioned at 3.5 Å from one another. Interestingly, the R357W mutation is also classified as an SW CAH variant, whereby its activity is slightly higher than that of W303R (∼0.01% and 0.1% for progesterone and 17α-OH-progesterone, respectively;
Table 2). Moreover, the levels of heme incorporation were quite similar for the two variants, ∼10% (
Table 3). The most likely reasons for the detrimental effects on activity and heme incorporation of these two mutations are altered positions and mobility of the I helix that spans the entire core of the protein and runs by the heme/active site. These changes in conformation and dynamic behavior occur either as a consequence of replacement of the Trp-303/Arg-357 pair by Arg/Arg (W303R variant) or Trp/Trp (R357W variant) pairs (in the model of the latter, the two indole moieties are ideally stacked at 3.4 Å,
Fig. 7E).
Overall, we find that single amino acid substitutions can have remarkable effects. All of the variants studied here were identified in clinical practice (
11- Pallan P.S.
- Lei L.
- Wang C.
- Waterman M.R.
- Guengerich F.P.
- Egli M.
Research resource: correlating human cytochrome P450 21A2 crystal structure and phenotypes of mutations in congenital adrenal hyperplasia.
). The summaries of catalytic activities of P450 21A2 activities (
11- Pallan P.S.
- Lei L.
- Wang C.
- Waterman M.R.
- Guengerich F.P.
- Egli M.
Research resource: correlating human cytochrome P450 21A2 crystal structure and phenotypes of mutations in congenital adrenal hyperplasia.
) (
http://www.cypalleles.ki.se/cyp21.htm)
7 are rather crude and contain values measured with many different expression systems and cellular and other assays. We provide here (
Table 2 and
Fig. 1) and elsewhere (
11- Pallan P.S.
- Lei L.
- Wang C.
- Waterman M.R.
- Guengerich F.P.
- Egli M.
Research resource: correlating human cytochrome P450 21A2 crystal structure and phenotypes of mutations in congenital adrenal hyperplasia.
) a robust set of
kcat,
Km, and
kcat/
Km values obtained under precise conditions with UPLC-UV measurements.
Although it is popular to use crystallographic results to consider alterations of substrate binding to explain properties of variants, our results suggest that the losses of catalytic activity (
Table 2) are too large to be explained by attenuation of substrate binding (Type I) (
24- Schenkman J.B.
- Remmer H.
- Estabrook R.W.
Spectral studies of drug interaction with hepatic microsomal cytochrome P-450.
). However, the heme perturbation assay may not necessarily reflect productive substrate binding. Overall, the changes of single amino acids often have dramatic effects on structural properties of P450 21A2. Evidence for this is seen in the incorporation of heme, both total heme (
Table 3) and the content of P450
versus inactive cytochrome P420 (
Fig. 6 and
supplemental Fig. S1). Some of the loss of heme from certain mutants might have occurred during purification, but large differences were also seen in the bacterial cells (
Fig. 6 and
supplemental Fig. S2).
The CD spectra of some of the purified variants indicated considerable distortion (
e.g. G292C, V282L, and L108R) (
Fig. 5). However, other variants with catalytic efficiencies just as low did not show the losses of α-helicity (
e.g. W303R). The loss of α-helicity with increasing temperature was no worse for the defective mutants than for the wild-type enzyme (
supplemental Fig. S1). Overall, the CD results by themselves did not provide explanations for the defective nature of the variants.
The Fe
2+-CO results clearly showed less stability of all variants (
Fig. 6 and
supplemental Fig. S2). Interestingly, the temperature dependence of these was rather invariant (
i.e. they started with a certain fraction of cytochrome P420, but the degree to which the residual P450 broke down was similar (V282L was somewhat faster;
Fig. 6E)).
Overall, we can conclude that a variety of issues are associated with the different clinical loss-of-function variants. The single amino acid substitutions can produce dramatic losses of functional integrity (
Fig. 5), the ability to bind the heme prosthetic group (
Table 3,
Fig. 6, and
supplemental Fig. S2), and the ability to bind substrates (
Table 4). Collectively, these problems result in low catalytic efficiency (
Fig. 1 and
Table 2). Nevertheless, C–H bond breaking appears to be the rate-limiting step across several orders of magnitude of catalytic efficiency (
Figure 1,
Figure 4 and
Table 2).
The ability of small changes in structure to influence large changes in catalytic efficiency seems remarkable, but it is of use to consider the Eyring equation,
where
kB is the Boltzmann constant,
h is Planck's constant (
); a mutation yielding a ΔΔ
G of 1.3 kcal mol
−1 leads to a 10-fold change in rate, and a ΔΔ
G of 6.4 kcal mol
−1 leads to a 50,000-fold change in
kobs.
In an earlier study of P450 21A2, we had mapped many SW, SV, and NC mutations onto the 3D structure of the enzyme and found patterns of distribution characteristic of the three types of CAH variants (
11- Pallan P.S.
- Lei L.
- Wang C.
- Waterman M.R.
- Guengerich F.P.
- Egli M.
Research resource: correlating human cytochrome P450 21A2 crystal structure and phenotypes of mutations in congenital adrenal hyperplasia.
). Thus, SW-causing mutations that are more common than SV variants were spread throughout the protein, often involving hydrophobic amino acids and dotting the active-site area and heme-binding pocket. By comparison, SV- and NC-causing mutations concern regions farther removed from the center and, particularly in the case of the NC variants, commonly map to areas on or near the enzyme surface, where they could be expected to be less damaging to stability and activity. The comparison of activities of SW, SV, and NC variants, based on various assays, confirmed the greater loss of activity in enzymes exhibiting the SW phenotype relative to those classified as SV and particularly NC (
11- Pallan P.S.
- Lei L.
- Wang C.
- Waterman M.R.
- Guengerich F.P.
- Egli M.
Research resource: correlating human cytochrome P450 21A2 crystal structure and phenotypes of mutations in congenital adrenal hyperplasia.
). The detailed analysis of a dozen variants presented here, using structural, activity, folding stability, and spectroscopic assays to measure heme incorporation, demonstrates in some cases that steady-state kinetic assays using a purified P450 21A2 variant and cell-based measurements of activity can produce considerably different outcomes. Thus, the most common NC variant, V282L, displayed much lower activity in the kinetic assays than in earlier cell-based experiments. Not only did this mutation impact activity to a considerable degree, but it also showed poor substrate binding and the greatest loss in thermodynamic stability among all of the CAH variants tested. Similarly, mutations at the surface can be as damaging as or more damaging than those at the heart of the enzyme, as demonstrated by the P31Q (SW) and the G292S/G292C (SW) variants, respectively. The crystal structure provides some insight into the origins of the relative degree of damage caused by individual mutations. Overall, our in-depth investigation of selected types of CAH variants, relying on an array of experimental approaches, provides new insight into disease-causing mutations in the P450 21A2 enzyme, which is among the metabolic enzymes with the largest number of naturally occurring variations (≫ 100) that alter activity.