Omega-hydroxylation of farnesol by mammalian cytochromes p450

Studies have shown that mammalian cytochromes p450 participate in the metabolism of terpenes, yet their role in the biotransformation of farnesol, an endogenous 15-carbon isoprenol, is unknown. In this report, [(14)C]-farnesol was transformed to more polar metabolites by NADPH-supplemented mammalian microsomes. In experiments with microsomes isolated from acetone-treated animals, the production of one polar metabolite was induced, suggesting catalysis by CYP2E1. The metabolite was identified as (2E, 6E, 10E)-12-hydroxyfarnesol. In studies with purified CYP2E1, 12-hydroxyfarnesol was obtained as the major product of farnesol metabolism. Among a series of available human p450 enzymes, only CYP2C19 also produced 12-hydroxyfarnesol. However, in individual human microsomes, CYP2E1 was calculated to contribute up to 62% toward total 12-hydroxyfarnesol production, suggesting CYP2E1 as the major catalyst. Mammalian cells expressing CYP2E1 demonstrated further farnesol metabolism to alpha,omega-prenyl dicarboxylic acids. Since such acids were identified in animal urine, the data suggest that CYP2E1 could be an important regulator of farnesol homeostasis in vivo. In addition, CYP2E1-dependent 12-hydroxyfarnesol formation was inhibited by pharmacological alcohol levels. Given that farnesol is a signaling molecule implicated in the regulation of tissue and cell processes, the biological activity of ethanol may be mediated in part by interaction with CYP2E1-dependent farnesol metabolism.     

Molecular definitions of fatty acid hydroxylases in Arabidopsis thaliana

Towards defining the function of Arabidopsis thaliana fatty acid hydroxylases, five members of the CYP86A subfamily have been heterologously expressed in baculovirus-infected Sf9 cells and tested for their ability to bind a range of fatty acids including unsubstituted (lauric acid (C12:0) and oleic acid (C18:1)) and oxygenated (9,10-epoxystearic acid and 9,10-dihydroxystearic acid). Comparison between these five P450s at constant P450 content over a range of concentrations for individual fatty acids indicates that binding of different fatty acids to CYP86A2 always results in a higher proportion of high spin state heme than binding titrations conducted with CYP86A1 or CYP86A4. In comparison to these three, CYP86A7 and CYP86A8 produce extremely low proportions of high spin state heme even with the most effectively bound fatty acids. In addition to their previously demonstrated lauric acid hydroxylase activities, all CYP86A proteins are capable of hydroxylating oleic acid but not oxygenated 9,10-epoxystearic acid. Homology models have been built for these five enzymes that metabolize unsubstituted fatty acids and sometimes bind oxygenated fatty acids. Comparison of the substrate binding modes and predicted substrate access channels indicate that all use channel pw2a consistent with the crystal structures and models of other fatty acid-metabolizing P450s in bacteria and mammals. Among these P450s, those that bind internally oxygenated fatty acids contain polar residues in their substrate binding cavity that help stabilize these charged/polar groups within their largely hydrophobic catalytic site.     

Progress in cytochrome P450 active site modeling

Models capable of predicting the possible involvement of cytochromes P450 in the metabolism of drugs or drug candidates are important tools in drug discovery and development. Ideally, functional information would be obtained from crystal structures of all the cytochromes P450 of interest. Initially, only crystal structures of distantly related bacterial cytochromes P450 were available-comparative modeling techniques were used to bridge the gap and produce structural models of human cytochromes P450, and thereby obtain some useful functional information. A significant step forward in the reliability of these models came four years ago with the first crystal structure of a mammalian cytochrome P450, rabbit CYP2C5, followed by the structures of two human enzymes, CYP2C8 and CYP2C9, and a second rabbit enzyme, CYP2B4. The evolution of a CYP2D6 model, leading to the validation of the model as an in silico tool for predicting binding and metabolism, is presented as a case study.     

Contribution of cytosolic cysteine residues to the gating properties of the Kir2.1 inward rectifier

The topological model proposed for the Kir2.1 inward rectifier predicts that seven of the channel 13 cysteine residues are distributed along the N- and C-terminus regions, with some of the residues comprised within highly conserved domains involved in channel gating. To determine if cytosolic cysteine residues contribute to the gating properties of Kir2.1, each of the N- and C-terminus cysteines was mutated into either a polar (S, D, N), an aliphatic (A,V, L), or an aromatic (W) residue. Our patch-clamp measurements show that with the exception of C76 and C311, the mutation of individual cytosolic cysteine to serine (S) did not significantly affect the single-channel conductance nor the channel open probability. However, mutating C76 to a charged or polar residue resulted either in an absence of channel activity or a decrease in open probability. In turn, the mutations C311S (polar), C311R (charged), and to a lesser degree C311A (aliphatic) led to an increase of the channel mean closed time due to the appearance of long closed time intervals (T(c) >or= 500 ms) and to a reduction of the reactivation by ATP of rundown Kir2.1 channels. These changes could be correlated with a weakening of the interaction between Kir2.1 and PIP(2), with C311R and C311S being more potent at modulating the Kir2.1-PIP(2) interaction than C311A. The present work supports, therefore, molecular models whereby the gating properties of Kir2.1 depend on the presence of nonpolar or neutral residues at positions 76 and 311, with C311 modulating the interaction between Kir2.1 and PIP(2).     

Massive endocytosis driven by lipidic forces originating in the outer plasmalemmal monolayer: a new approach to membrane recycling and lipid domains

The pore-lining amino acids of ion channel proteins reside on the interface between a polar (the pore) and a nonpolar environment (the rest of the protein). The structural dynamics of this region, which physically controls ionic flow, are essential components of channel gating. Using large-conductance, Ca(2+)-dependent K(+) (BK) channels, we devised a systematic charge-substitution method to probe conformational changes in the pore region during channel gating. We identified a deep-pore residue (314 in hSlo1) as a marker of structural dynamics. We manipulated the charge states of this residue by substituting amino acids with different valence and pKa, and by adjusting intracellular pH. We found that the charged states of the 314 residues stabilized an open state of the BK channel. With models based on known structures of related channels, we postulate a dynamic rearrangement of the deep-pore region during BK channel opening/closing, which involves a change of the degree of pore exposure for 314.     

Exploring the structure characteristics and major channels of cytochrome P450 2A6, 2A13, and 2E1 with pilocarpine

The majority of cytochromes P450 play a critical role in metabolism of endogenous and exogenous substrates, some of its products are carcinogens. Therefore, inhibition of P450 enzymes activity can promote the detoxification and elimination of chemical carcinogens. In this study, molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed to explore the structure features and channel dynamics of three P450 isoforms 2A6, 2A13, and 2E1 bound with the common inhibitor pilocarpine. The binding free energy results combined with the PMF calculations give a reasonable ranking of binding affinity, which are consistent with the experimental data. Our results uncover how a sequence divergence of different CYP2 enzymes causes individual variations in major channel selections. On the basis of channel bottleneck and energy decomposition analysis, we propose a gating mechanism of their respective major channels in three enzymes, which may be attributed to a reversal of Phe209 in CYP2A6/2A13, as well as the rotation of Phe116 and Phe298 in CYP2E1. The hydrophobic residues not only make strong hydrophobic interactions with inhibitor, but also act as gatekeeper to regulate the opening of channel. The present study provides important insights into the structure–function relationships of three cytochrome P450s and the molecular basis for development of potent inhibitors.     

Modelling human cytochromes P450 involved in drug metabolism from the CYP2C5 crystallographic template

A historical background to homology modelling of human P450s involved in drug metabolism is outlined, showing that the progress in crystallographic studies of bacterial forms of enzyme and, latterly, determination of a mammalian P450 crystal structure, has enabled the production of increasingly satisfactory models of human P450 enzymes. The methodology for the generation of P450 models by homology with crystallographic template structures is summarized, and recent results of CYP2C5-constructed models of P450s are described. These indicate that selective substrates are able to fit within the putative active sites of each enzyme, where key contacts with complementary amino acid residues are largely consistent with the results of site-directed mutagenesis experiments and metabolic studies. Consequently, the CYP2C5 crystal structure can be regarded at the current paradigm for homology modelling of the drug metabolizing P450s, especially those from the CYP2 family.     

Substrate recognition sites in cytochrome P450 family 2 (CYP2) proteins inferred from comparative analyses of amino acid and coding nucleotide sequences.

The substrate recognition regions in cytochrome P450 family 2 (CYP2) proteins were inferred by group-to-group alignment of CYP2 sequences and those of bacterial P450s, including Pseudomonas putida P450 101A (P450cam), whose substrate-binding residues have been definitely identified by x-ray crystallography of a substrate-bound form (Poulos T. L., Finzel, B. C., and Howard, A. J. (1987) J. Mol. Biol. 195, 687-700). The six putative substrate recognition sites, SRSs, thus identified are dispersively located along the primary structure and constitute about 16% of the total residues. All the reported point mutations and chimeric fragments that significantly affect the substrate specificities of the parental CYP2 enzymes fell within or overlapped some of the six SRSs. Analysis of nucleotide substitution patterns in closely related members in four subfamilies, CYP2A, 2B, 2C, and 2D, consistently indicated that the SRSs have accumulated more nonsynonymous (amino acid-changing) substitutions than the rest of the sequence. This observation supports the idea that diversification of duplicate genes of drug-metabolizing P450s occurs primarily in substrate recognition regions to cope with an increasing number of foreign compounds.     

Homology model of 1alpha,25-dihydroxyvitamin D3 24-hydroxylase cytochrome P450 24A1 (CYP24A1): active site architecture and ligand binding

Homology models of cytochrome P450 24A1 (CYP24A1) were constructed using three human P450 structures, CYP2C8, CYP2C9 and CYP3A4 as templates for the model building. Using molecular operating environment (MOE) software the lowest energy CYP24A1 model was then assessed for stereochemical quality and side chain environment. Further active site optimisation of the CYP24A1 model built using the CYP3A4 template was performed by molecular dynamics to generate a final CYP24A1 model. The natural substrate, 1,25-dihydroxyvitamin D(3) (calcitriol) and the CYP24 inhibitor (R)-N-(2-(1H-imidazol-1-yl)-2-phenylethyl)-4'-chlorobiphenyl-4-carboxamide ((R)-VID-400) were docked into the model allowing further validation of the active site architecture. Using the docking studies structurally and functionally important residues were identified with subsequent characterisation of secondary structure.     

Homology modeling and substrate binding study of human CYP2C9 enzyme

The CYP2C subfamily of human liver P450 isozymes is of major importance in drug metabolism. The most abundant 2C isozyme, CYP2C9, regioselectively hydroxylates a wide variety of substrates. A major obstacle to understanding this specificity in human CYP2C9 is the absence of a 3D structure. A 3D model of CYP2C9 was built, assessed, and used to characterize explicit enzyme-substrate complexes using methods previously developed in our laboratory. The 3D model was assessed by determining its stability to unconstrained molecular dynamics and by comparison of specific properties with those of known protein structures. The CYP2C9 model was then used to characterize explicit enzyme complexes with three structurally and chemically diverse substrates: (S)-naproxen, phenytoin, and progesterone. Each substrate was found to bind to the enzyme with a favorable interaction energy and to remain in the binding site during unconstrained molecular dynamics. Moreover, the mode of binding of each substrate led to calculated preferred hydroxylation sites consistent with experiment. Binding-site residues identified for the models included Arg 105 and Arg97 as key cationic residues, as well as Asn 202, Asp 293, Pro 101, Leu 102, Gly 296, and Phe 476. Site-specific mutations are proposed for further integrated computational and experimental study.     

On the human CYP2C9*13 variant activity reduction: a molecular dynamics simulation and docking study

Cytochrome P450 2C9 (CYP2C9) plays a key role in the metabolism of clinical drugs. CYP2C9 is a genetically polymorphic enzyme and some of its allelic variants have less activity compared to the wild-type form. Drugs with a narrow therapeutic index may cause serious toxicity to the individuals who carry such allele. CYP2C9*13, firstly identified by some of the present authors in a Chinese poor metabolizer of lornoxicam, is characterized by mutation encoding Leu90Pro substitution. Kinetic experiments show that CYP2C9*13 has less catalytic activity in elimination of diclofenac and lornoxicam in vitro. In order to explore the structure-activity relationship of CYP2C9*13, the three-dimensional structure models of the substrate-free CYP2C9*1 and its variant CYP2C9*13 are constructed on the basis of the X-ray crystal structure of human CYP2C9*1 (PDB code 1R9O) by molecular dynamics simulations. The structure change caused by Leu90Pro replacement is revealed and used to explain the dramatic decrease of the enzymatic activity in clearance of the two CYP2C9 substrates: diclofenac and lornoxicam. The trans configuration of the bond between Pro90 and Asp89 in CYP2C9*13 is firstly identified. The backbone of residues 106-108 in CYP2C9*13 turns over and their side chains block the entrance for substrates accessing so that the entrance of *13 shrinks greatly than that in the wild-type, which is believed to be the dominant mechanism of the catalytic activity reduction. Consequent docking study which is consistent with the results of the kinetic experiments by Guo et al. identifies the most important residues for enzyme-substrate complexes.     

Homology model of human retinoic acid metabolising enzyme cytochrome P450 26A1 (CYP26A1): active site architecture and ligand binding

Homology models of cytochrome P450 RA1 (CYP26A1) were constructed using three human P450 structures, CYP2C8, CYP2C9 and CYP3A4 as templates for the model building. Using MOE software the lowest energy CYP26A1 model was then assessed for stereochemical quality and side chain environment. Further active site optimisation of the CYP26A1 model built using the CYP3A4 template was performed by molecular dynamics to generate a final CYP26A1 model. The natural substrate, all-trans-retinoic acid (atRA), and inhibitor R 15866, were docked into the model allowing further validation of the active site architecture. Using the docking studies structurally and functionally important residues were identified with subsequent characterisation of secondary structure. Multiple hydrophobic interactions, including the side chains of TRP112, PHE299, PHE222, PHE84, PHE374 and PRO371, are important for binding of atRA and R115866. Additional hydrogen bonding interactions were noted as follows: atRA-- C==O of the atRA carboxylate group and ARG86; R115866--benzothiazole nitrogen and the backbone NH of SER115.     

Three-dimensional model of cytochrome P450 human aromatase

A three-dimensional (3-D) structure of human aromatase (CYP 19) was modeled on the basis of the crystal structure of rabbit CYP2C5, the first solved X-ray structure of an eukaryotic cytochrome P450 and was evaluated by docking S-fadrozole and the steroidal competitive inhibitor (19R)-10-thiiranylestr-4-ene-3,17-dione, into the enzyme active site. According to a previous pharmacophoric hypothesis described in the literature, the cyano group of S-fadrozole partially mimics the steroid backbone C(17) carbonyl group of (19R)-10-thiiranylestr-4-ene-3,17-dione, and was oriented in a favorable position for H-bonding with the newly identified positively charged residues Lys 119 and Arg435. In addition, this model is consistent with the recent combined mutagenesis/modeling studies already published concerning the roles ofAsp309 and His480 in the aromatization of the steroid A ring.     

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.     

Homology modeling, molecular dynamics simulations, and analysis of CYP119, a P450 enzyme from extreme acidothermophilic archaeon Sulfolobus solfataricus

The recent characterization of a thermophilic and barophilic CYP119 from Sulfolobus solfataricus offers a new opportunity to identify the origin of its stability by comparing it with mesophilic P450s with known structures. Since the three-dimensional structure of CYP119 is not yet available, homology modeling techniques were used to build model structures for this enzyme. The overall quality and stability of the models were assessed using three protein analysis programs and by monitoring structural stability during 1 ns of molecular dynamics simulations at 300 and 390 K. The results show the CYP119 models to be of good quality. Possible origins of the thermo- and barostability of CYP119 were then investigated by examining the amino acid compositions and the three-dimensional structure of CYP119 compared with the five mesophilic templates. Three possible factors were identified that could contribute to the enhanced stability of CYP119. The first was the higher relative population of salt bridges and the presence of a few unique salt bridges found in CYP119 that were absent in all five template CYP450s. The second factor was a decreased population of Ala and an increased population of Ile found in the interior of CYP119, which are likely to improve packing in CYP119. The third factor was a more extensive aromatic cluster seen in CYP119 which was not found in all five template P450s. In addition, the model CYP119 three-dimensional structures were also used to determine key properties related to its function. Specifically, binding site residues and surface residues for redox partner interactions were identified. These residues identified together with those residues found that might contribute to the increased stability are suggested for future mutagenesis studies. The results obtained from these experimental studies shall then provide further validation of the suggested origins of stability and the structure-function relationships derived from the model structures of this enzyme.     

The polar T1 interface is linked to conformational changes that open the voltage-gated potassium channel

Kv voltage-gated potassium channels share a cytoplasmic assembly domain, T1. Recent mutagenesis of two T1 C-terminal loop residues implicates T1 in channel gating. However, structural alterations of these mutants leave open the question concerning direct involvement of T1 in gating. We find in mammalian Kv1.2 that gating depends critically on residues at complementary T1 surfaces in an unusually polar interface. An isosteric mutation in this interface causes surprisingly little structural alteration while stabilizing the closed channel and increasing the stability of T1 tetramers. Replacing T1 with a tetrameric coiled-coil destabilizes the closed channel. Together, these data suggest that structural changes involving the buried polar T1 surfaces play a key role in the conformational changes leading to channel opening.     

Homology modeling of the three membrane proteins of the dhurrin metabolon: catalytic sites, membrane surface association and protein-protein interactions

Formation of metabolons (macromolecular enzyme complexes) facilitates the channelling of substrates in biosynthetic pathways. Metabolon formation is a dynamic process in which transient structures mediated by weak protein-protein interactions are formed. In Sorghum, the cyanogenic glucoside dhurrin is derived from l-tyrosine in a pathway involving the two cytochromes P450 (CYPs) CYP79A1 and CYP71E1, a glucosyltransferase (UGT85B1), and the redox partner NADPH-dependent cytochrome P450 reductase (CPR). Experimental evidence suggests that the enzymes of this pathway form a metabolon. Homology modeling of the three membrane bound proteins was carried out using the Sybyl software and available relevant crystal structures. Residues involved in tight positioning of the substrates and intermediates in the active sites of CYP79A1 and CYP71E1 were identified. In both CYPs, hydrophobic surface domains close to the N-terminal trans-membrane anchor and between the F′ and G helices were identified as involved in membrane anchoring. The proximal surface of both CYPs showed positively charged patches complementary to a negatively charged bulge on CPR carrying the FMN domain. A patch of surface exposed, positively charged amino acid residues positioned on the opposite face of the membrane anchor was identified in CYP71E1 and might be involved in binding UGT85B1 via a hypervariable negatively charged loop in this protein.     

Identification of Val117 and Arg372 as critical amino acid residues for the activity difference between human CYP2A6 and CYP2A13 in coumarin 7-hydroxylation

Human cytochrome P450 (CYP) 2A6 and 2A13 play an important role in catalyzing the metabolism of many environmental chemicals including coumarin, nicotine, and several tobacco-specific carcinogens. Both CYP2A6 and CYP2A13 proteins are composed of 494 amino acid residues. Although CYP2A13 shares a 93.5% identity with CYP2A6 in the amino acid sequence, it is only about one-tenth as active as CYP2A6 in catalyzing coumarin 7-hydroxylation. To identify the key amino acid residues that account for such a remarkable difference, we generated a series of CYP2A6 and CYP2A13 mutants by site-directed mutagenesis/heterologous expression and compared their coumarin 7-hydroxylation activities. In CYP2A6, the amino acid residues at position 117 and 372 are valine (Val) and arginine (Arg), respectively; whereas in CYP2A13, they are alanine (Ala) and histidine (His). Kinetic analysis revealed that the catalytic efficiency (Vmax/Km) of the CYP2A6 Val(117)--> Ala and Arg(372)--> His mutants was drastically reduced (0.41 and 0.64 versus 3.23 for the wild-type CYP2A6 protein). In contrast, the catalytic efficiency of the CYP2A13 Ala(117) --> Val and His(372) --> Arg mutants was greatly increased (2.65 and 2.60 versus 0.31 for wild-type CYP2A13 protein). These results clearly demonstrate that the Val at position 117 and Arg at position 372 are critical amino acid residues for coumarin 7-hydroxylation. Based on the crystal structure of CYP2C5, we have generated the homology models of CYP2A6 and CYP2A13 and docked the substrate coumarin to the active site. Together with the kinetic characterization, our structural modeling provides explanations for the amino acid substitution results and the insights of detailed enzyme-substrate interactions.     

Phylogenetic and Functional Analysis of the Vertebrate Cytochrome P450 2 Family

Cytochrome P450 (CYP) proteins compose a highly diverse superfamily found in all domains of life. These proteins are enzymes involved in metabolism of endogenous and exogenous compounds. In vertebrates, the CYP2 family is one of the largest, most diverse and plays an important role in mammalian drug metabolism. However, there are more than 20 vertebrate CYP2 subfamilies with uncertain evolution and fairly discrete subfamily composition within vertebrate classes, hindering extrapolation of knowledge across subfamilies. To better understand CYP2 diversity, a phylogenetic analysis of 196 CYP2 protein sequences from 16 species was performed using a maximum likelihood approach and Bayesian inference. The analyses included the CYP2 compliment from human, fugu, zebrafish, stickleback, medaka, cow, and dog genomes. Additional sequences were included from rabbit, marsupial, platypus, chicken, frog, and salmonid species. Three CYP2 sequences from the tunicate Ciona intestinalis were utilized as the outgroup. Results indicate a single ancestral vertebrate CYP2 gene and monophyly of all CYP2 subfamilies. Two subfamilies (CYP2R and CYP2U) pre-date vertebrate diversification, allowing direct comparison across vertebrate classes, while all other subfamilies originated during vertebrate diversification, often within specific vertebrate lineages. Analysis of site-specific evolution indicates that some substrate recognition sites (SRS) previously proposed for CYP genes do not have elevated rates of evolution, suggesting that these regions of the protein are not necessarily important in recognition of CYP2 substrates. Type II functional divergence analysis identified multiple residues in the active site of CYP2F, CYP2A, and CYP2B proteins that have undergone radical biochemical changes and may be functionally important.     

Molecular characterization and chromosomal localization of cytochrome P450 genes involved in the biosynthesis of cyclic hydroxamic acids in hexaploid wheat

The cyclic hydroxamic acids, 2,4-dihydroxy-1,4-benzoxazin-3-one (DIBOA) and 2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one (DIMBOA), are defensive secondary metabolites found in gramineous plants including wheat, maize and rye. cDNAs for five cytochromes P450 (P450s) involved in DIBOA biosynthesis (CYP71C6, CYP71C7v2, CYP71C8v2, CYP71C9v1 and CYP71C9v2) were isolated from seedlings of hexaploid wheat [( Triticum aestivum L. cv. Chinese Spring (2n=6x=42, genomes AABBDD)] by RT-PCR and screening of a cDNA library. CYP71C9v1 and CYP71C9v2 are 97% identical to each other in amino acid and nucleotide sequences. The cloned P450 species showed 76-79% identity at the amino acid level to the corresponding maize P450 species CYP71C1-C4, which are also required for DIBOA biosynthesis. The wheat P450 cDNAs were heterologously expressed in the yeast ( Saccharomyces cerevisiae) strain AH22. Microsome fractions from yeast cells expressing these P450 species catalyzed the same reactions as their maize orthologs. The chromosomes carrying the cyp71C6- C9v1 orthologs were identified by Southern hybridization using aneuploid lines of Chinese Spring wheat. The cyp71C9v1 orthologs were located on the chromosomes of wheat homoeologous group-4. The orthologs of the other P450 genes, cyp71C7v2, cyp71C6 and cyp71C8v2, were located on group-5 chromosomes. The same P450 genes were also present in the three ancestral diploid species of hexaploid wheat, T. monococcum (AA), Aegilops speltoides [BB (approximately SS)] and Ae. squarrosa (DD).