Structure of the human lung cytochrome P450 2A13

The human lung cytochrome P450 2A13 (CYP2A13) activates the nicotine-derived procarcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) into DNA-altering compounds that cause lung cancer. Another cytochrome P450, CYP2A6, is also present in human lung, but at much lower levels. Although these two enzymes are 93.5% identical, CYP2A13 metabolizes NNK with much lower K(m) values than does CYP2A6. To investigate the structural differences between these two enzymes the structure of CYP2A13 was determined to 2.35A by x-ray crystallography and compared with structures of CYP2A6. As expected, the overall CYP2A13 and CYP2A6 structures are very similar with an average root mean square deviation of 0.5A for the Calpha atoms. Like CYP2A6, the CYP2A13 active site cavity is small and highly hydrophobic with a cluster of Phe residues composing the active site roof. Active site residue Asn(297) is positioned to hydrogen bond with an adventitious ligand, identified as indole. Amino acid differences between CYP2A6 and CYP2A13 at positions 117, 300, 301, and 208 relate to different orientations of the ligand plane in the two protein structures and may underlie the significant variations observed in binding and catalysis of many CYP2A ligands. In addition, docking studies suggest that residues 365 and 366 may also contribute to differences in NNK metabolism.     

Reverse type I inhibitor of Mycobacteriumtuberculosis CYP125A1

Cytochrome P450 CYP125A1 of Mycobacteriumtuberculosis, a potential therapeutic target for tuberculosis in humans, initiates degradation of the aliphatic chain of host cholesterol and is essential for establishing M. tuberculosis infection in a mouse model of disease. We explored the interactions of CYP125A1 with a reverse type I inhibitor by X-ray structure analysis and UV-vis spectroscopy. Compound LP10 (α-[(4-methylcyclohexyl)carbonyl amino]-N-4-pyridinyl-1H-indole-3-propanamide), previously identified as a potent type II inhibitor of Trypanosomacruzi CYP51, shifts CYP125A1 to a water-coordinated low-spin state upon binding with low micromolar affinity. When LP10 is present in the active site, the crystal structure and spectral characteristics both demonstrate changes in lipophilic and electronic properties favoring coordination of the iron axial water ligand. These results provide an insight into the structural requirements for developing selective CYP125A1 inhibitors.     

Nicotine and 4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone Binding and Access Channel in Human Cytochrome P450 2A6 and 2A13 Enzymes

Cytochromes P450 (CYP) from the 2A subfamily are known for their roles in the metabolism of nicotine, the addictive agent in tobacco, and activation of the tobacco procarcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK). Although both the hepatic CYP2A6 and respiratory CYP2A13 enzymes metabolize these compounds, CYP2A13 does so with much higher catalytic efficiency, but the structural basis for this has been unclear. X-ray structures of nicotine complexes with CYP2A13 (2.5 Å) and CYP2A6 (2.3 Å) yield a structural rationale for the preferential binding of nicotine to CYP2A13. Additional structures of CYP2A13 with NNK reveal either a single NNK molecule in the active site with orientations corresponding to metabolites known to form DNA adducts and initiate lung cancer (2.35 Å) or with two molecules of NNK bound (2.1 Å): one in the active site and one in a more distal staging site. Finally, in contrast to prior CYP2A structures with enclosed active sites, CYP2A13 conformations were solved that adopt both open and intermediate conformations resulting from an ∼2.5 Å movement of the F to G helices. This channel occurs in the same region where the second, distal NNK molecule is bound, suggesting that the channel may be used for ligand entry and/or exit from the active site. Altogether these structures provide multiple new snapshots of CYP2A13 conformations that assist in understanding the binding and activation of an important human carcinogen, as well as critical comparisons in the binding of nicotine, one of the most widely used and highly addictive drugs in human use.     

Structural comparison of cytochromes P450 2A6, 2A13, and 2E1 with pilocarpine

Human xenobiotic-metabolizing cytochrome P450 (CYP) enzymes can each bind and monooxygenate a diverse set of substrates, including drugs, often producing a variety of metabolites. Additionally, a single ligand can interact with multiple CYP enzymes, but often the protein structural similarities and differences that mediate such overlapping selectivity are not well understood. Even though the CYP superfamily has a highly canonical global protein fold, there are large variations in the active site size, topology, and conformational flexibility. We have determined how a related set of three human CYP enzymes bind and interact with a common inhibitor, the muscarinic receptor agonist drug pilocarpine. Pilocarpine binds and inhibits the hepatic CYP2A6 and respiratory CYP2A13 enzymes much more efficiently than the hepatic CYP2E1 enzyme. To elucidate key residues involved in pilocarpine binding, crystal structures of CYP2A6 (2.4 ?), CYP2A13 (3.0 ?), CYP2E1 (2.35 ?), and the CYP2A6 mutant enzyme, CYP2A6?I208S/I300F/G301A/S369G (2.1 ?) have been determined with pilocarpine in the active site. In all four structures, pilocarpine coordinates to the heme iron, but comparisons reveal how individual residues lining the active sites of these three distinct human enzymes interact differently with the inhibitor pilocarpine.     

Structural Analysis of the Streptomyces avermitilis CYP107W1-Oligomycin A Complex and Role of the Tryptophan 178 Residue

CYP107W1 from Streptomyces avermitilis is a cytochrome P450 enzyme involved in the biosynthesis of macrolide oligomycin A. A previous study reported that CYP107W1 regioselectively hydroxylated C12 of oligomycin C to produce oligomycin A, and the crystal structure of ligand free CYP107W1 was determined. Here, we analyzed the structural properties of the CYP107W1-oligomycin A complex and characterized the functional role of the Trp178 residue in CYP107W1. The crystal structure of the CYP107W1 complex with oligomycin A was determined at a resolution of 2.6 Å. Oligomycin A is bound in the substrate access channel on the upper side of the prosthetic heme mainly by hydrophobic interactions. In particular, the Trp178 residue in the active site intercalates into the large macrolide ring, thereby guiding the substrate into the correct binding orientation for a productive P450 reaction. A Trp178 to Gly mutation resulted in the distortion of binding titration spectra with oligomycin A, whereas binding spectra with azoles were not affected. The Gly178 mutant’s catalytic turnover number for the 12-hydroxylation reaction of oligomycin C was highly reduced. These results indicate that Trp178, located in the open pocket of the active site, may be a critical residue for the productive binding conformation of large macrolide substrates.     

First structures of an active bacterial tyrosinase reveal copper plasticity

Tyrosinase is a member of the type 3 copper enzyme family that is involved in the production of melanin in a wide range of organisms. The crystal structures of a tyrosinase from Bacillus megaterium were determined at a resolution of 2.0-2.3 Å. The enzyme crystallized as a dimer in the asymmetric unit and was shown to be active in crystal. The overall monomeric structure is similar to that of the monomer of the previously determined tyrosinase from Streptomyces castaneoglobisporus, but it does not contain an accessory Cu-binding “caddie” protein. Two Cu(II) ions, serving as the major cofactors within the active site, are coordinated by six conserved histidine residues. However, determination of structures under different conditions shows varying occupancies and positions of the copper ions. This apparent mobility in copper binding modes indicates that there is a pathway by which copper is accumulated or lost by the enzyme. Additionally, we suggest that residues R209 and V218, situated in a second shell of residues surrounding the active site, play a role in substrate binding orientation based on their flexibility and position. The determination of a structure with the inhibitor kojic acid, the first tyrosinase structure with a bound ligand, revealed additional residues involved in the positioning of substrates in the active site. Comparison of wild-type structures with the structure of the site-specific variant R209H, which possesses a higher monophenolase/diphenolase activity ratio, lends further support to a previously suggested mechanism by which monophenolic substrates dock mainly to CuA.     

N-Myristoyltransferase from Leishmania donovani: Structural and Functional Characterisation of a Potential Drug Target for Visceral Leishmaniasis

N-Myristoyltransferase (NMT) catalyses the attachment of the 14-carbon saturated fatty acid, myristate, to the amino-terminal glycine residue of a subset of eukaryotic proteins that function in multiple cellular processes, including vesicular protein trafficking and signal transduction. In these pathways, N-myristoylation facilitates association of substrate proteins with membranes or the hydrophobic domains of other partner peptides. NMT function is essential for viability in all cell types tested to date, demonstrating that this enzyme has potential as a target for drug development. Here, we provide genetic evidence that NMT is likely to be essential for viability in insect stages of the pathogenic protozoan parasite, Leishmania donovani, causative agent of the tropical infectious disease, visceral leishmaniasis. The open reading frame of L. donovaniNMT has been amplified and used to overproduce active recombinant enzyme in Escherichia coli, as demonstrated by gel mobility shift assays of ligand binding and peptide-myristoylation activity in scintillation proximity assays. The purified protein has been crystallized in complex with the non-hydrolysable substrate analogue S-(2-oxo)pentadecyl-CoA, and its structure was solved by molecular replacement at 1.4 Å resolution. The structure has as its defining feature a 14-stranded twisted β-sheet on which helices are packed so as to form an extended and curved substrate-binding groove running across two protein lobes. The fatty acyl-CoA is largely buried in the N-terminal lobe, its binding leading to the loosening of a flap, which in unliganded NMT structures, occludes the protein substrate binding site in the carboxy-terminal lobe. These studies validate L. donovani NMT as a potential target for development of new therapeutic agents against visceral leishmaniasis.     

In-silico and In-vitro based studies of Streptomyces peucetius CYP107N3 for oleic acid epoxidation

Certain members of the cytochromes P450 superfamily metabolize polyunsaturated long-chain fatty acids to several classes of oxygenated metabolites. An approach based on in silico analysis predicted that Streptomyces peucetius CYP107N3 might be a fatty acid-metabolizing enzyme, showing high homology with epoxidase enzymes. Homology modeling and docking studies of CYP107N3 showed that oleic acid can fit directly into the active site pocket of the double bond of oleic acid within optimum distance of 4.6 Å from the Fe. In order to confirm the epoxidation activity proposed by in silico analysis, a gene coding CYP107N3 was expressed in Escherichia coli. The purified CYP107N3 was shown to catalyze C₉-C₁₀ epoxidation of oleic acid in vitro to 9,10-epoxy stearic acid confirmed by ESI-MS, HPLC-MS and GC-MS spectral analysis. [BMB Reports 2012; 45(12): 736-741]     

Identification of mitochondrial cytochrome P450 induced in response to polycyclic aromatic hydrocarbons in the mummichog (Fundulus heteroclitus)

Increasing evidence suggests that polycyclic aromatic hydrocarbons (PAHs) such as benzo[a]pyrene (BaP) are localized to the mitochondria. Because the toxic effects of many PAHs are the result of metabolism by cytochrome P4501A (CYP1A), it is important to investigate whether active forms of these enzymes can be identified in the mitochondria. In this study, we identified mitochondrial P450s with a monoclonal antibody against scup (Stenotomus chrysops) CYP1A in the isolated mitochondrial fraction of the liver from adult male mummichog (Fundulus heteroclitus) livers. The size of the protein in the mitochondria was similar to that of microsomal CYP1A. Fish dosed with 10 mg/kg BaP had increased EROD activity in the mitochondrial fraction compared to controls. In mummichog larvae dosed with 100 µg/L BaP and 100 µg/L benzo[k]fluoranthene, CYP1A protein levels as well as enzyme activity were elevated. However, fish from a PAH-polluted Superfund site (Elizabeth River, Portsmouth VA) showed recalcitrant mitochondrial CYP1A protein levels and enzyme activity in a similar manner to microsomal CYP1A.     

Ligand binding to cytochrome P450 3A4 in phospholipid bilayer nanodiscs: the effect of model membranes

The membrane-bound protein cytochrome P450 3A4 (CYP3A4) is a major drug-metabolizing enzyme. Most studies of ligand binding by CYP3A4 are currently carried out in solution, in the absence of a model membrane. Therefore, there is little information concerning the membrane effects on CYP3A4 ligand binding behavior. Phospholipid bilayer Nanodiscs are a novel model membrane system derived from high density lipoprotein particles, whose stability, monodispersity, and consistency are ensured by their self-assembly. We explore the energetics of four ligands (6-(p-toluidino)-2-naphthalenesulfonic acid (TNS), alpha-naphthoflavone (ANF), miconazole, and bromocriptine) binding to CYP3A4 incorporated into Nanodiscs. Ligand binding to Nanodiscs was monitored by a combination of environment-sensitive ligand fluorescence and ligand-induced shifts in the fluorescence of tryptophan residues present in the scaffold proteins of Nanodiscs; binding to the CYP3A4 active site was monitored by ligand-induced shifts in the heme Soret band absorbance. The dissociation constants for binding to the active site in CYP3A4-Nanodiscs were 4.0 microm for TNS, 5.8 microm for ANF, 0.45 microm for miconazole, and 0.45 microm for bromocriptine. These values are for CYP3A4 incorporated into a lipid bilayer and are therefore presumably more biologically relevant that those measured using CYP3A4 in solution. In some cases, affinity measurements using CYP3A4 in Nanodiscs differ significantly from solution values. We also studied the equilibrium between ligand binding to CYP3A4 and to the membrane. TNS showed no marked preference for either environment; ANF preferentially bound to the membrane, and miconazole and bromocriptine preferentially bound to the CYP3A4 active site.     

Crystal structure of CYP199A2, a para-substituted benzoic acid oxidizing cytochrome P450 from Rhodopseudomonas palustris

CYP199A2, a cytochrome P450 enzyme from Rhodopseudomonas palustris, oxidatively demethylates 4-methoxybenzoic acid to 4-hydroxybenzoic acid. 4-Ethylbenzoic acid is converted to a mixture of predominantly 4-(1-hydroxyethyl)-benzoic acid and 4-vinylbenzoic acid, the latter being a rare example of CC bond dehydrogenation of an unbranched alkyl group. The crystal structure of CYP199A2 has been determined at 2.0-Å resolution. The enzyme has the common P450 fold, but the B′ helix is missing and the G helix is broken into two (G and G′) by a kink at Pro204. Helices G and G′ are bent back from the extended BC loop and the I helix to open up a clearly defined substrate access channel. Channel openings in this region of the P450 fold are rare in bacterial P450 enzymes but more common in eukaryotic P450 enzymes. The channel is hydrophobic except for the basic residue Arg246 at the entrance, which probably plays a role in the specificity of this enzyme for charged benzoates over neutral phenols and benzenes. The substrate binding pocket is hydrophobic, with Ser97 and Ser247 being the only polar residues. Computer docking of 4-ethylbenzoic acid into the active site suggests that the substrate carboxylate oxygens interact with Ser97 and Ser247, and the β-methyl group is located over the heme iron by Phe185, the side chain of which is only 6.35 Å above the iron in the native structure. This binding orientation is consistent with the observed product profile of exclusive attack at the para substituent. Putidaredoxin of the CYP101A1 system from Pseudomonas putida supports substrate oxidation by CYP199A2 at ∼6% of the activity of the physiological ferredoxin. Comparison of the heme proximal faces of CYP199A2 and CYP101A1 suggests that charge reversal surrounding the surface residue Leu369 in CYP199A2 may be a significant factor in this low cross-activity.     

Crystallographic and Functional Characterization of the Fluorodifen-inducible Glutathione Transferase from Glycine max Reveals an Active Site Topography Suited for Diphenylether Herbicides and a Novel L-site

Glutathione transferases (GSTs) from the tau class (GSTU) are unique to plants and have important roles in stress tolerance and the detoxification of herbicides in crops and weeds. A fluorodifen-induced GST isoezyme (GmGSTU4-4) belonging to the tau class was purified from Glycine max by affinity chromatography. This isoenzyme was cloned and expressed in Escherichia coli, and its structural and catalytic properties were investigated. The structure of GmGSTU4-4 was determined at 1.75 Å resolution in complex with S-(p-nitrobenzyl)-glutathione (Nb-GSH). The enzyme adopts the canonical GST fold but with a number of functionally important differences. Compared with other plant GSTs, the three-dimensional structure of GmGSTU4-4 primarily shows structural differences in the hydrphobic substrate binding site, the linker segment and the C-terminal region. The X-ray structure identifies key amino acid residues in the hydrophobic binding site (H-site) and provides insights into the substrate specificity and catalytic mechanism of the enzyme. The isoenzyme was highly active in conjugating the diphenylether herbicide fluorodifen. A possible reaction pathway involving the conjugation of glutathione with fluorodifen is described based on site-directed mutagenesis and molecular modeling studies. A serine residue (Ser13) is present in the active site, at a position that would allow it to stabilise the thiolate anion of glutathione and enhance its nucleophilicity. Tyr107 and Arg111 present in the active site are important structural moieties that modulate the catalytic efficiency and specificity of the enzyme, and participate in k_cat regulation by affecting the rate-limiting step of the catalytic reaction. A hitherto undescribed ligand-binding site (L-site) located in a surface pocket of the enzyme was also found. This site is formed by conserved residues, suggesting it may have an important functional role in the transfer and delivery of bound ligands, presumably to specific protein receptors.     

Structure–function analysis of cytochromes P450 2B

In the last 4 years, breakthroughs were made in the field of P450 2B (CYP2B) structure–function through determination of one ligand-free and two inhibitor-bound X-ray crystal structures of CYP2B4, which revealed many of the structural features required for binding ligands of different size and shape. Large conformational changes of several plastic regions of CYP2B4 can dramatically reshape the active site of the enzyme to fit the size and shape of the bound ligand without perturbing the overall P450 fold. Solution biophysical studies using isothermal titration calorimetry (ITC) have revealed the large difference in the thermodynamic parameters of CYP2B4 in binding inhibitors of different ring chemistry and side chains. Other studies have revealed that the effects of site-specific mutations on steady-state kinetic parameters and mechanism-based inactivation are often substrate dependent. These findings agree with the structural data that the enzymes adopt different conformations to bind various ligands. Thus, the substrate specificity of an individual enzyme is determined not only by active site residues but also non-active site residues that modulate conformational changes that are important for substrate access and rearrangement of the active site to accommodate the bound substrate.     

Structural Insights into the Substrate Specificity of (S)-Ureidoglycolate Amidohydrolase and Its Comparison with Allantoate Amidohydrolase

In plants, the ureide pathway is a metabolic route that converts the ring nitrogen atoms of purine into ammonia via sequential enzymatic reactions, playing an important role in nitrogen recovery. In the final step of the pathway, (S)-ureidoglycolate amidohydrolase (UAH) catalyzes the conversion of (S)-ureidoglycolate into glyoxylate and releases two molecules of ammonia as by-products. UAH is homologous in structure and sequence with allantoate amidohydrolase (AAH), an upstream enzyme in the pathway with a similar function as that of an amidase but with a different substrate. Both enzymes exhibit strict substrate specificity and catalyze reactions in a concerted manner, resulting in purine degradation. Here, we report three crystal structures of Arabidopsis thaliana UAH (bound with substrate, reaction intermediate, and product) and a structure of Escherichia coli AAH complexed with allantoate. Structural analyses of UAH revealed a distinct binding mode for each ligand in a bimetal reaction center with the active site in a closed conformation. The ligand directly participates in the coordination shell of two metal ions and is stabilized by the surrounding residues. In contrast, AAH, which exhibits a substrate-binding site similar to that of UAH, requires a larger active site due to the additional ureido group in allantoate. Structural analyses and mutagenesis revealed that both enzymes undergo an open-to-closed conformational transition in response to ligand binding and that the active-site size and the interaction environment in UAH and AAH are determinants of the substrate specificities of these two structurally homologous enzymes.     

Scanning chimeragenesis: the approach used to change the substrate selectivity of fatty acid monooxygenase CYP102A1 to that of terpene ω-hydroxylase CYP4C7

CYP102A1 is a highly active, water-soluble, bacterial monooxygenase enzyme that contains both substrate-binding heme and diflavin reductase subunits, both in a single polypeptide. Recently we developed a procedure which uses the known structure of the substrate-bound heme domain of CYP102A1 and its sequence homology with a cytochrome P450 of unknown structure, both of which react with a common substrate but produce different products, to create recombinant enzymes which have substrate selectivity different from that of CYP102A1, and produce the product of the enzyme of unknown structure. Insect CYP4C7, a terpene hydroxylase from the cockroach, was chosen as the cytochrome P450 of unknown structure, and farnesol was chosen as the substrate. CYP102A1 oxidizes farnesol to three products (2,3-epoxyfarnesol, 10,11-epoxyfarnesol, and 9-hydroxyfarnesol), whereas CYP4C7 produces 12-hydroxyfarnesol as the major product. In earlier work it was found that the chimera C(78-82,F87L) showed a change in substrate selectivity from fatty acids to farnesol, and was approximately sixfold more active than wild-type CYP102A1 (Chen et al. in J Biol Inorg Chem 13:813–824, 2008), but neither it nor any other earlier chimera produced 12-hydroxyfarnesol. In this work we added amino acid residues 327–332, to create six new full-length, functional chimeric proteins. Four of these, the most active of which was C(78-82,F87L,328-330), produce 12-hydroxyfarnesol as the major product, with approximately twofold increase in turnover number as compared with wild-type CYP102A1 toward farnesol. Methylfarnesoate was metabolized to 12-hydroxymethylfarnesoate (70%) and 10,11-epoxymethylfarnesoate (juvenile hormone III) (30%). The latter is metabolized to 65% 12-hydroxy-10,11-epoxymethylfarnesoate and 35% 15-hydroxy-10,11-epoxymethylfarnesoate. Substitution of residues 328–330, APA, by VPL was crucial to accomplishing this change in product.     

Structure of microsomal cytochrome P450 2B4 complexed with the antifungal drug bifonazole: insight into P450 conformational plasticity and membrane interaction

To better understand ligand-induced structural transitions in cytochrome P450 2B4, protein-ligand interactions were investigated using a bulky inhibitor. Bifonazole, a broad spectrum antifungal agent, inhibits monooxygenase activity and induces a type II binding spectrum in 2B4dH(H226Y), a modified enzyme previously crystallized in the presence of 4-(4-chlorophenyl)imidazole (CPI). Isothermal titration calorimetry and tryptophan fluorescence quenching indicate no significant burial of protein apolar surface nor altered accessibility of Trp-121 upon bifonazole binding, in contrast to recent results with CPI. A 2.3 A crystal structure of 2B4-bifonazole reveals a novel open conformation with ligand bound in the active site, which is significantly different from either the U-shaped cleft of ligand-free 2B4 or the small active site pocket of 2B4-CPI. The O-shaped active site cleft of 2B4-bifonazole is widely open in the middle but narrow at the top. A bifonazole molecule occupies the bottom of the active site cleft, where helix I is bent approximately 15 degrees to accommodate the bulky ligand. The structure also defines unanticipated interactions between helix C residues and bifonazole, suggesting an important role of helix C in azole recognition by mammalian P450s. Comparison of the ligand-free 2B4 structure, the 2B4-CPI structure, and the 2B4-bifonazole structure identifies structurally plastic regions that undergo correlated conformational changes in response to ligand binding. The most plastic regions are putative membrane-binding motifs involved in substrate access or substrate binding. The results allow us to model the membrane-associated state of P450 and provide insight into how lipophilic substrates access the buried active site.     

Metabolism of anabolic steroids by recombinant human cytochrome P450 enzymes: Gas chromatographic–mass spectrometric determination of metabolites

Metabolism of steroid hormones with anabolic properties was studied in vitro using human recombinant CYP3A4, CYP2C9 and 2B6 enzymes. The enzyme formats used for CYP3A4 and CYP2C9 were insect cell microsomes expressing human CYP enzymes and purified recombinant human CYP enzymes in a reconstituted system. CYP3A4 enzyme formats incubated with anabolic steroids, testosterone, 17α-methyltestosterone, metandienone, boldenone and 4-chloro-1,2-dehydro-17α-methyltestosterone, produced 6β-hydroxyl metabolites identified as trimethylsilyl (TMS)-ethers by a gas chromatography–mass spectrometry (GC–MS) method. When the same formats of CYP2C9 were incubated with the anabolic steroids, no 6β-hydroxyl metabolites were formed. Human lymphoblast cell microsomes expressing human CYP2B6 incubated with the steroids investigated produced traces of 6β-hydroxyl metabolites with testosterone and 17α-methyltestosterone only. We suggest that the electronic effects of the 3-keto-4-ene structural moiety contribute to the selectivity within the active site of CYP3A4 enzyme resulting in selective 6β-hydroxylation.     

Stereoselective Formation and Metabolism of 4-Hydroxy-Retinoic Acid Enantiomers by Cytochrome P450 Enzymes

All-trans-retinoic acid (atRA), the major active metabolite of vitamin A, plays a role in many biological processes, including maintenance of epithelia, immunity, and fertility and regulation of apoptosis and cell differentiation. atRA is metabolized mainly by CYP26A1, but other P450 enzymes such as CYP2C8 and CYP3As also contribute to atRA 4-hydroxylation. Although the primary metabolite of atRA, 4-OH-RA, possesses a chiral center, the stereochemical course of atRA 4-hydroxylation has not been studied previously. (4S)- and (4R)-OH-RA enantiomers were synthesized and separated by chiral column HPLC. CYP26A1 was found to form predominantly (4S)-OH-RA. This stereoselectivity was rationalized via docking of atRA in the active site of a CYP26A1 homology model. The docked structure showed a well defined niche for atRA within the active site and a specific orientation of the β-ionone ring above the plane of the heme consistent with stereoselective abstraction of the hydrogen atom from the pro-(S)-position. In contrast to CYP26A1, CYP3A4 formed the 4-OH-RA enantiomers in a 1:1 ratio and CYP3A5 preferentially formed (4R)-OH-RA. Interestingly, CYP3A7 and CYP2C8 preferentially formed (4S)-OH-RA from atRA. Both (4S)- and (4R)-OH-RA were substrates of CYP26A1 but (4S)-OH-RA was cleared 3-fold faster than (4R)-OH-RA. In addition, 4-oxo-RA was formed from (4R)-OH-RA but not from (4S)-OH-RA by CYP26A1. Overall, these findings show that (4S)-OH-RA is preferred over (4R)-OH-RA by the enzymes regulating atRA homeostasis. The stereoselectivity observed in CYP26A1 function will aid in better understanding of the active site features of the enzyme and the disposition of biologically active retinoids.