Molecular modelling of steroidogenic cytochromes P450 from families CYP11, CYP17, CYP19 and CYP21 based on the CYP102 crystal structure

Polymorphism in the Androgen Biosynthesis Gene (CYP17), a Risk for Prostate Cancer: A Meta-Analysis

Gene polymorphism is one of the few factors that increases the risk of prostate cancer. T to C substitution in the 5’ promoter region of the CYP17 gene is hypothesized to increase the rate of gene transcription, increase androgen production, and thereby increase the risk of prostate cancer. Nevertheless, the inconsistencies originating from studies on CYP17 polymorphism and prostate cancer prompted this meta-analysis, to decipher the association between CYP17 polymorphism and prostate cancer. Most case-control studies addressing CYP17 polymorphism and prostate cancer were exhaustively searched from Web of Science, Google Scholar, and PubMed. The various genotype distributions as well as the minor allele distributions were retrieved. Pooled odds ratios (ORs) with their 95% CI and estimates of the Hardy–Weinberg Equilibrium were calculated. Analyses were performed using the RevMan v.5.3 software and SPSS v.21. There was high-pooled heterogeneity (I² = 87.0%, OR = .42, CI [.39, .45], and p < .001) among the A2 versus A1 allele. With the per-allele model (A2 versus A1), ethnicity was a major risk factor to prostate cancer, with Asians recording the highest risk (OR = 12.61, 95% CI [8.77, 18.12]). From the genotype models, A1/A1 versus A2/A2 (OR = 3.02, 95% CI [2.65, 3.44]) and A1/A2 versus A2/A2 (OR = 4.39, 95% CI [3.86, 5.00]) were all significantly associated with prostate cancer. Although some genotype models were associated with the risk of prostate cancer, we should be mindful when interpreting the results of this study because of the limited number of studies and the small sample size used.

Ile351, Leu355 and Ile461 residues are essential for catalytic activity of bovine cytochrome P450scc (CYP11A1)

Cytochrome P450scc (CYP11A1) is a mammalian mitochondrial enzyme which catalyzes cholesterol side chain cleavage to form pregnenolone. Along with cholesterol, some other steroids including sterols with a branched side chain like β-sitosterol are the substrates for the enzyme, but the activity towards β-sitosterol is rather low. Modification of the catalytic site conformation could provide more effective β-sitosterol bioconversion by the enzyme. This study was aimed to find out the amino acid residues substitution of which could modify the conformation of the active site providing possible higher enzyme activity towards β-sitosterol. After structural and bioinformatics analysis three amino acid residues I351, L355, I461 were chosen. Molecular dynamics simulations of P450scc evidenced the stability of the wild type, double (I351A/L355A) and triple (I351A/L355A/I461A) mutants. Mutant variants of cDNA encoding P450scc with the single, double and triple mutations were obtained by site-directed mutagenesis. However, the experimental data indicate that the introduced single mutations Ile351A, Leu355A and Ile461A dramatically decrease the target catalytic activity of CYP11A1, and no activity was observed for double and triple mutants obtained. Therefore, isoleucine residues 351 and 461, and leucine residue 355 are important for the cytochrome P450scc functioning towards sterols both with unbranched (cholesterol) and branched (sitosterol) side chains.

Molecular identification of steroidogenesis-related genes in scallops and their potential roles in gametogenesis

Sex steroids are crucial for controlling gametogenesis and germ cell maturation in vertebrates. It has been proposed that Yesso scallop (Mizuhopecten yessoensis) has the same sex steroids as those animals, but the scallop biosynthetic pathway is unclear. In this study, we characterized several steroidogenesis-related genes in M. yessoensis and proposed a putative biosynthetic pathway for sex steroids that is similar to that of vertebrates. Specifically, we identified several steroidogenesis-related gene sequences that encode steroid metabolizing enzymes: StAR-related lipid transfer (START) protein, 17α-hydroxylase, 17,20-lyase (cyp17a), 17β-hydroxysteroid dehydrogenase (hsd17b), and 3β-hydroxysteroid dehydrogenase (hsd3b). We sampled adult scallops throughout their reproductive phase to compare their degree of maturation with their intensity of mRNA expression. Semi-quantitative RT-PCR analysis revealed a ubiquitous expression of transcripts for steroid metabolizing enzymes (i.e., star, cyp17a, hsd17b, and hsd3b) in peripheral and gonadal tissues. Real-time PCR analysis revealed a high level of expression of star3 and cyp17a genes in gonadal tissues at the early stage of cell differentiation in scallops. Interestingly, mRNA expression of hsd3b and hsd17b genes showed a synchronous pattern related to degree of gonad maturity. These results indicate that both hsd3b and hsd17b genes are likely involved in steroidogenesis in scallops. We therefore believe that these steroid-metabolizing enzymes allow scallops to endogenously produce sex steroids to regulate reproductive events.

Catalytically relevant electrostatic interactions of cytochrome P450c17 (CYP17A1) and cytochrome b5

Two acidic residues, Glu-48 and Glu-49, of cytochrome b5 (b5) are essential for stimulating the 17,20-lyase activity of cytochrome P450c17 (CYP17A1). Substitution of Ala, Gly, Cys, or Gln for these two glutamic acid residues abrogated all capacity to stimulate 17,20-lyase activity. Mutations E49D and E48D/E49D retained 23 and 38% of wild-type activity, respectively. Using the zero-length cross-linker ethyl-3-(3-dimethylaminopropyl)carbodiimide, we obtained cross-linked heterodimers of b5 and CYP17A1, wild-type, or mutations R347K and R358K. In sharp contrast, the b5 double mutation E48G/E49G did not form cross-linked complexes with wild-type CYP17A1. Mass spectrometric analysis of the CYP17A1-b5 complexes identified two cross-linked peptide pairs as follows: CYP17A1-WT: (84)EVLIKK(89)-b5: (53)EQAGGDATENFEDVGHSTDAR(73) and CYP17A1-R347K: (341)TPTISDKNR(349)-b5: (40)FLEEHPGGEEVLR(52). Using these two sites of interaction and Glu-48/Glu-49 in b5 as constraints, protein docking calculations based on the crystal structures of the two proteins yielded a structural model of the CYP17A1-b5 complex. The appositional surfaces include Lys-88, Arg-347, and Arg-358/Arg-449 of CYP17A1, which interact with Glu-61, Glu-42, and Glu-48/Glu-49 of b5, respectively. Our data reveal the structural basis of the electrostatic interactions between these two proteins, which is critical for 17,20-lyase activity and androgen biosynthesis.

Evolutionary origin of the mitochondrial cholesterol transport machinery reveals a universal mechanism of steroid hormone biosynthesis in animals

Steroidogenesis begins with the transport of cholesterol from intracellular stores into mitochondria via a series of protein-protein interactions involving cytosolic and mitochondrial proteins located at both the outer and inner mitochondrial membranes. In adrenal glands and gonads, this process is accelerated by hormones, leading to the production of high levels of steroids that control tissue development and function. A hormone-induced multiprotein complex, the transduceosome, was recently identified, and is composed of cytosolic and outer mitochondrial membrane proteins that control the rate of cholesterol entry into the outer mitochondrial membrane. More recent studies unveiled the steroidogenic metabolon, a bioactive, multimeric protein complex that spans the outer-inner mitochondrial membranes and is responsible for hormone-induced import, segregation, targeting, and metabolism of cholesterol by cytochrome P450 family 11 subfamily A polypeptide 1 (CYP11A1) in the inner mitochondrial membrane. The availability of genome information allowed us to systematically explore the evolutionary origin of the proteins involved in the mitochondrial cholesterol transport machinery (transduceosome, steroidogenic metabolon, and signaling proteins), trace the original archetype, and predict their biological functions by molecular phylogenetic and functional divergence analyses, protein homology modeling and molecular docking. Although most members of these complexes have a history of gene duplication and functional divergence during evolution, phylogenomic analysis revealed that all vertebrates have the same functional complex members, suggesting a common mechanism in the first step of steroidogenesis. An archetype of the complex was found in invertebrates. The data presented herein suggest that the cholesterol transport machinery is responsible for steroidogenesis among all vertebrates and is evolutionarily conserved throughout the entire animal kingdom.

CYP17 T27C polymorphism and prostate cancer risk: a meta-analysis based on 31 studies

Objective

The cytochrome P450 17α-hydroxylase (CYP17) plays a vital role in androgen biosynthesis. A T-to-C polymorphism in the 5′ promoter region of CYP17 has been implicated as a risk factor for prostate cancer, but the results of individual studies are inconclusive or controversial. To derive a more precise estimation of the relationship, we performed an updated meta-analysis from 31 studies based on 27 publications.

Methods

A comprehensive search was conducted to examine all the eligible studies of CYP17 polymorphism and prostate cancer risk. We used odds ratios (ORs) with 95% confidence intervals (CIs) to assess the strength of the association.

Results

Overall, individuals with CC/CT genotype were not associated with prostate cancer risk (CC vs. TT: OR = 1.03, 95% CI = 0.86-1.24, P = 0.72, P_{heterogeneity} < 0.0001; CT vs. TT: OR = 0.99, 95% CI = 0.87-1.12, P = 0.88, P_{heterogeneity} = 0.0006). In the stratified analysis by ethnicity, there was a significantly increased risk of prostate cancer among individuals of African descent under the recessive model (OR = 1.56, 95% CI = 1.01-2.39, P = 0.04, P_{heterogeneity} = 0.65).

Conclusion

This meta-analysis suggested that CYP17 polymorphism might be associated with prostate cancer risk among individuals of African descent.

An analysis of CYP19, CYP21 and ER genotypes in Polish Holstein-Friesian cows with regard to the selected reproductive traits

The aim of this study was to relate polymorphic variants of CYP19, CYP21 and ER1 genes to reproductive traits in 472 Polish Holstein-Friesian cows. High frequencies of one of the homozygous genotypes were found. The ER1/SnaBIAA homozygotes were not identified. In the first and third lactation, an average calving-to-conception interval (CLVC) in cows of ER1/SnaBIGG genotype was significantly shorter (P ≤ 0.05) than in heterozygous cows. In the cows of ER1/BglIGG genotype, significantly shorter CLVC (P ≤ 0.05) was observed compared to heterozygotes in the first lactation, whereas in the third lactation, CLVC in homozygous cows was significantly longer (P ≤ 0.05) than in heterozygous ones. It was also found that homozygous cows were characterized by significantly longer calving interval (CLVI; P ≤ 0.05) compared to heterozygotes in the third lactation. Longer CLVCs in CYP19AA cows were found, compared to heterozygotes, and this difference was significant in the first and third lactation (P ≤ 0.05). Similarly, the average CLVIs were longer in CYP19AA homozygotes than in heterozygous cows; however, significance was proven only in the third lactation (P ≤ 0.05). Description of the molecular mechanisms regulating reproduction, and thus identification of the individuals of genotypes with optimal potential may facilitate the employment of selected reproductive model by a breeder.

Systematic genetic and genomic analysis of cytochrome P450 enzyme activities in human liver

Liver cytochrome P450s (P450s) play critical roles in drug metabolism, toxicology, and metabolic processes. Despite rapid progress in the understanding of these enzymes, a systematic investigation of the full spectrum of functionality of individual P450s, the interrelationship or networks connecting them, and the genetic control of each gene/enzyme is lacking. To this end, we genotyped, expression-profiled, and measured P450 activities of 466 human liver samples and applied a systems biology approach via the integration of genetics, gene expression, and enzyme activity measurements. We found that most P450s were positively correlated among themselves and were highly correlated with known regulators as well as thousands of other genes enriched for pathways relevant to the metabolism of drugs, fatty acids, amino acids, and steroids. Genome-wide association analyses between genetic polymorphisms and P450 expression or enzyme activities revealed sets of SNPs associated with P450 traits, and suggested the existence of both cis-regulation of P450 expression (especially for CYP2D6) and more complex trans-regulation of P450 activity. Several novel SNPs associated with CYP2D6 expression and enzyme activity were validated in an independent human cohort. By constructing a weighted coexpression network and a Bayesian regulatory network, we defined the human liver transcriptional network structure, uncovered subnetworks representative of the P450 regulatory system, and identified novel candidate regulatory genes, namely, EHHADH, SLC10A1, and AKR1D1. The P450 subnetworks were then validated using gene signatures responsive to ligands of known P450 regulators in mouse and rat. This systematic survey provides a comprehensive view of the functionality, genetic control, and interactions of P450s.

Molecular modeling on inhibitor complexes and active-site dynamics of cytochrome P450 C17, a target for prostate cancer therapy

A molecular model for the P450 enzyme cytochrome P450 C17 (CYP17) is presented based on sequence alignments of multiple template structures and homology modeling. This enzyme plays a central role in the biosynthesis of testosterone and is emerging as a major target in prostate cancer, with the recently developed inhibitor abiraterone currently in advanced clinical trials. The model is described in detail, together with its validation, by providing structural explanations to available site-directed mutagenesis data. The CYP17 molecule in this model is in the form of a triangular prism, with an edge of approximately 55 A and a thickness of approximately 37 A. It is predominantly helical, comprising 13 alpha helices interspersed by six 3(10) helices and 11 beta-sheets. Multinanosecond molecular dynamics simulations in explicit solvent have been carried out, and principal components analysis has been used to reveal the details of dynamics around the active site. Coarse-grained methods have also been used to verify low-frequency motions, which have been correlated with active-site gating. The work also describes the results of docking synthetic inhibitors, including the drug abiraterone and the natural substrate pregnenolone, in the CYP17 active site together with molecular dynamics simulations on the complexes.

A single amino acid residue, Ala 105, confers 16alpha-hydroxylase activity to human cytochrome P450 17alpha-hydroxylase/17,20 lyase

In adrenal steroidogenesis, CYP17 catalyses the 17alpha-hydroxylation of pregnenolone and progesterone and the subsequent 17,20-lyase reaction, yielding adrenal androgens. The enzyme exhibits distinctly different selectivities towards these substrates in various species. CYP17 has also been shown to exhibit 16alpha-hydroxylase activity towards progesterone in some species, with only human and chimp CYP17 catalysing the biosynthesis of substantial amounts of 16-OHprogesterone. The 16alpha-hydroxylase activity was investigated by introducing an Ala105Leu substitution into human CYP17. The converse mutation, Leu105Ala was introduced into the baboon, goat and pig enzymes. Wt human CYP17 converted approximately 30% progesterone to 16-OHprogesterone while the Ala105Leu mutant converted negligible amounts to 16-OHprogesterone ( approximately 9%), comparable to wt CYP17 of the other three species when expressed in COS-1 cells. The ratio of 17-hydroxylated products to 16-OHprogesterone of human CYP17 was 2.7 and that of the mutant human construct 10.5. Similar ratios were observed for human and goat CYP17 with the corresponding Ala or Leu residues. Although the Leu105Ala mutation of both baboon and pig CYP17 exhibited the same trend regarding the ratios, the rate of progesterone conversion was reduced. Coexpression with cytochrome b(5) significantly decreased the ratio of 17-hydroxylated products to 16-OHprogesterone in the Leu105 constructs, while effects were negligible with Ala at this position. Homology models show that Ala105 faces towards the active pocket in the predicted B'-C domain of CYP17. The smaller residue allows more flexibility of movement in the active pocket than Leu, presenting both the C16 and C17 of progesterone to the iron-oxy complex.

CYP17 inhibition as a hormonal strategy for prostate cancer

Androgen receptor (AR) signaling has a key role in the pathogenesis of prostate cancer. AR gene amplification, AR overexpression, and activating mutations in the AR occur more frequently as castration-resistant prostate cancer (CRPC) evolves, with intratumoral androgen levels remaining sufficient for AR activation despite castration. The source of these androgens might be either adrenal or intratumoral. AR signaling, therefore, remains a valid treatment target for patients with CRPC. CYP17 is a key enzyme for androgen biosynthesis. The imidazole antifungal agent ketoconazole weakly and nonspecifically inhibits CYP17, but remains unlicensed for this indication. Chemists at the Cancer Research UK Centre for Cancer Therapeutics have designed a novel, selective, irreversible inhibitor of CYP17 called abiraterone, which is more than 20 times more potent than ketoconazole. Abiraterone acetate, a prodrug, has undergone phase I assessment, and is rapidly progressing from phase II to phase III trials, in view of its high level of antitumor activity. This agent is safe and well tolerated, and activity profiles suggest that approximately 50% of CRPC remains AR-ligand driven. Other CYP17 inhibitors with alternative mechanisms of action, for example VN/124-1, are in preclinical development. The rationale for and implications of CYP17 inhibition and the CYP17-targeting agents in development are discussed in this Review.

Comparison of the 17α-Hydroxylase/C17,20-Lyase Activities of Porcine, Guinea Pig and Bovine P450c17 Using Purified Recombinant Fusion Proteins Containing P450c17 Linked to NADPH-P450 Reductase

The cDNAs for cytochrome P450c17 (P450c17) of three species, pig, guinea pig, and cow, representing three families of mammals (suidae, procaviidae, and bovidae, respectively) were each engineered into an expression plasmid (pCWori+). The P450c17 domain of the coding sequence was connected to a truncated form of rat NADPH-P450 reductase by a linker sequence encoding two amino acids (SerThr). These fusion proteins were expressed in E. coli and purified for use in enzymatic assays to determine similarities and differences in 17 alpha-hydroxylase and lyase activities. The fusion proteins were found to catalyze both the 17 alpha-hydroxylation of progesterone (P4) and pregnenolone (P5) to 17 alpha-hydroxylated P4 and P5 (17 alpha-OH P4 and 17 alpha-OH P5) followed by the C17,20-lyase reaction for the conversion of these C(21)-17 alpha-hydroxylated steroids to C(19)-steroids (the C17,20-lyase reaction). These in vitro studies show that (a) porcine P450c17 possesses cytochrome b(5) (b(5))-stimulated C17,20-lyase activity that converts 17 alpha OH-P4 to androstenedione (AD) but also converts 17 alpha-OHP5 to dehydroepiandrosterone (DHEA); (b) guinea pig P450c17 possesses a b(5)-stimulated C17,20-lyase activity that converts 17 alpha-OH P4 to AD but does not convert 17 alpha-OH P5 to DHEA., and (c) bovine P450c17 possesses a b(5)-stimulated C17,20-lyase activity that converts 17 alpha-OH P5 to DHEA but does not convert 17 alpha-OH P4 to AD. Thus, the P450c17 of each species differs in its ability to catalyze in vitro the conversion of C(21)-steroids to C(19)-steroids. In addition, each P450c17 is capable of catalyzing additional hydroxylation reactions leading to low levels of 2 alpha-, 6 beta-, 16- and 21-hydroxy-metabolites. Porcine P450c17 also catalyzes the b(5)-dependent synthesis of andien-beta (androsta-5,16-dien-3beta-ol) from P5. When the amino acid sequences of the three P450c17s were aligned there was an approximate 50% variation in the alignment identity (227 differences in the sequences of 509 amino acids). Alignment did not permit the assignment of specific amino acids or domains to the observed differences in enzymatic activities.

Genetic polymorphisms in CYP17, CYP3A4, CYP19A1, SRD5A2, IGF-1, and IGFBP-3 and prostate cancer risk in African-American men: the Flint Men's Health Study

BACKGROUND

Association studies have examined the significance of several candidate genes based on biological pathways relevant to prostate carcinogenesis, including both the androgen and insulin-like growth factor pathways. Clinical and epidemiologic evidence suggest that androgens, specifically testosterone and dihydrotestosterone (DHT) are important not only in normal prostate growth but in the pathogenesis of prostate cancer. Similarly, the insulin-like growth factor-1 (IGF-1) signaling pathway regulates both cellular proliferation and apoptosis. Therefore, genes involved in the biosynthesis, activation, metabolism and degradation of androgens and the stimulation of mitogenic and antiapoptotic activities of prostate epithelial cells represent important candidates for affecting the development and progression of prostate cancer.

METHODS

Using resources from the Flint Men's Health Study, a population-based case control study of African-American men aged 40-79, we evaluated the associations between selected single-nucleotide polymorphisms (SNPs) in the CYP17, CYP3A4, CYP19A1, SDR5A2, IGF1, and IGFBP3 genes and prostate cancer diagnosis in 473 men (131 prostate cancer cases and 342 disease-free controls).

RESULTS

We found a significant association between prostate cancer and selected CYP17 SNP genotypes, with the heterozygous genotype conferring decreased risk. Suggestive evidence for association between IGF1 SNPs and prostate cancer were also found. No significant associations were observed between SNPs in the other genes and prostate cancer.

CONCLUSIONS

These findings suggest that variation in or around CYP17 and/or IGF1 may be associated with prostate cancer development in the African-American population. Additional studies are needed to determine whether these polymorphisms are indeed associated with prostate cancer risk in African Americans.

A three-dimensional model of CYP19 aromatase for structure-based drug design

Aromatase (CYP450(arom), CYP19) is an enzyme responsible for converting the aliphatic androgens androstenedione and testosterone to the aromatic estrogens estrone and estradiol, respectively. These endogenous hormones are a key factor in cancer tumor formation and proliferation through a cascade starting from estrogen binding to estrogen receptor. To interfere with the overproduction of estrogens especially in tumor tissue, it is possible to inhibit aromatase activity. This can be achieved using aromatase inhibitors. In order to design novel aromatase inhibitors, it is necessary to have an understanding of the active site of aromatase. As no crystal structure of the enzyme has yet been published, we built a homology model of aromatase using the first crystallized mammalian cytochrome enzyme, rabbit 21-progesterone hydroxylase 2C5, as a template structure. The initial model was validated with exhaustive molecular dynamics simulation with and without the natural substrate androstenedione. The resulting enzyme-substrate complex shows very good stability and only two of the residues are in disallowed regions in a Ramachandran plot.

Ensemble modeling of substrate binding to cytochromes P450: analysis of catalytic differences between CYP1A orthologs

A novel application of modeling and docking approaches involving ensembles of homology models is used to understand structural bases underlying subtle catalytic differences between related cytochromes P450 (CYPs). Mammalian CYP1A1s and fish CYP1As are orthologous enzymes with similar substrate preferences. With some substrates (3,3',4,4'-tetrachlorobiphenyl, TCB) oxidation rates differ by orders of magnitude, while others (e.g., benzo[a]pyrene; B[a]P) are oxidized at similar rates but with somewhat differing regiospecificity. These two environmental chemical substrates (TCB and B[a]P) as well as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) were docked to multiple models of rat, human, scup, and/or killifish CYP1As, based on multiple templates, retaining multiple poses from each model, giving ensembles of docked poses for each species. With TCB, more poses were observed closer to the heme in ensembles of rat or human CYP1A1 than of killifish CYP1A. Analysis of interacting residues suggested that differences in TCB pose distributions are due primarily to Leu387 and Val230 in killifish CYP1A. In silico mutations L387V and V230G enabled TCB to dock closer to the heme in killifish CYP1A. Mutating additional interacting residues (Ala127, Thr233, Asn317, and Tyr386) of killifish CYP1A to the corresponding residues of human CYP1A1 resulted in TCB pose distributions nearly identical with those of human CYP1A1. Docking of TCDD to sets of consensus models of killifish, rat, and human CYP1As showed species differences similar to those with TCB, but with further structural constraints possibly contributing to slower oxidation of TCDD. Docking B[a]P to sets of consensus models of the human and fish CYP1As yielded frequencies of substrate orientations correlating with known regiospecificities for metabolism of B[a]P by these enzymes. The results demonstrate the utility of this ensemble modeling method, which can account for uncertainty inherent in homology modeling and docking by producing statistical distributions of ligand positions.

Cytochrome P450 side-chain cleavage: insights gained from homology modeling

Cytochrome P450 side-chain cleavage (CYP11A1) catalyzes the conversion of cholesterol to pregnenolone, the first step in steroidogenesis. The absence of a solved crystal structure has complicated deductions pertaining to the structure/function relationships of this key enzyme. Although a number of techniques have been employed to identify domains and specific amino acid residues important for catalytic activity, these methods have been unsuccessful in predicting three-dimensional orientations in space and thus the mechanism by which they exert their kinetic effect. This review aims to demonstrate the significant contribution homology modelling, when employed as a tool in combination with other standard biochemical techniques, has made towards our understanding of CYP11A1.

Methods for Predicting Human Drug Metabolism

Drug metabolism information is a necessary component of drug discovery and development. The key issues in drug metabolism include identifying: the enzyme(s) involved, the site(s) of metabolism, the resulting metabolite(s), and the rate of metabolism. Methods for predicting human drug metabolism from in vitro and computational methodologies and determining relationships between the structure and metabolic activity of molecules are also critically important for understanding potential drug interactions and toxicity. There are numerous experimental and computational approaches that have been developed in order to predict human metabolism which have their own limitations. It is apparent that few of the computational tools for metabolism prediction alone provide the major integrated functions needed to assist in drug discovery. Similarly the different in vitro methods for human drug metabolism themselves have implicit limitations. The utilization of these methods for pharmaceutical and other applications as well as their integration is discussed as it is likely that hybrid methods will provide the most success.

Evidence for the functional role of residues in the B'-C loop of baboon cytochrome P450 side-chain cleavage (CYP11A1) obtained by site-directed mutagenesis, kinetic analysis and homology modelling

To gain further insight into the structure/function relationship of cytochrome P450 side-chain cleavage (CYP11A1), this enzyme was investigated in the Cape baboon (Papio ursinus). Four constructs were cloned and characterised in non-steroidogenic mammalian COS-1 cells. Wild type recombinant baboon CYP11A1 cDNA yielded a K(m) value of 1.6 microM for 25-hydroxycholesterol. The single amino acid substitutions, I98Q and I98K resulted in a 1.7- and 2.8-fold increases in K(m) values, respectively. Conversely, the introduction of the mutation, K103A, resulted in a 1.8-fold decrease in K(m). A homology model of CYP11A1, based on the crystal structures of CYP102 and CYP2C5, revealed that residues 98 and 103 lie within the B'-C loop and contribute to the spatial orientation and structural integrity of this domain. Based on these results we propose a topological model of the CYP11A1 active pocket, which is supported by substrate docking analysis and kinetic studies.

Molecular model of human CYP21 based on mammalian CYP2C5: structural features correlate with clinical severity of mutations causing congenital adrenal hyperplasia

Enhanced understanding of structure-function relationships of human 21-hydroxylase, CYP21, is required to better understand the molecular causes of congenital adrenal hyperplasia. To this end, a structural model of human CYP21 was calculated based on the crystal structure of rabbit CYP2C5. All but two known allelic variants of missense type, a total of 60 disease-causing mutations and six normal variants, were analyzed using this model. A structural explanation for the corresponding phenotype was found for all but two mutants for which available clinical data are also discrepant with in vitro enzyme activity. Calculations of protein stability of modeled mutants were found to correlate inversely with the corresponding clinical severity. Putative structurally important residues were identified to be involved in heme and substrate binding, redox partner interaction, and enzyme catalysis using docking calculations and analysis of structurally determined homologous cytochrome P450s (CYPs). Functional and structural consequences of seven novel mutations, V139E, C147R, R233G, T295N, L308F, R366C, and M473I, detected in Scandinavian patients with suspected congenital adrenal hyperplasia of different severity, were predicted using molecular modeling. Structural features deduced from the models are in good correlation with clinical severity of CYP21 mutants, which shows the applicability of a modeling approach in assessment of new CYP21 mutations.

Functional interaction of nitrogenous organic bases with cytochrome P450: A critical assessment and update of substrate features and predicted key active-site elements steering the access, binding, and orientation of amines

The widespread use of nitrogenous organic bases as environmental chemicals, food additives, and clinically important drugs necessitates precise knowledge about the molecular principles governing biotransformation of this category of substrates. In this regard, analysis of the topological background of complex formation between amines and P450s, acting as major catalysts in C- and N-oxidative attack, is of paramount importance. Thus, progress in collaborative investigations, combining physico-chemical techniques with chemical-modification as well as genetic engineering experiments, enables substantiation of hypothetical work resulting from the design of pharmacophores or homology modelling of P450s. Based on a general, CYP2D6-related construct, the majority of prospective amine-docking residues was found to cluster near the distal heme face in the six known SRSs, made up by the highly variant helices B', F and G as well as the N-terminal portion of helix C and certain beta-structures. Most of the contact sites examined show a frequency of conservation < 20%, hinting at the requirement of some degree of conformational versatility, while a limited number of amino acids exhibiting a higher level of conservation reside close to the heme core. Some key determinants may have a dual role in amine binding and/or maintenance of protein integrity. Importantly, a series of non-SRS elements are likely to be operative via long-range effects. While hydrophobic mechanisms appear to dominate orientation of the nitrogenous compounds toward the iron-oxene species, polar residues seem to foster binding events through H-bonding or salt-bridge formation. Careful uncovering of structure-function relationships in amine-enzyme association together with recently developed unsupervised machine learning approaches will be helpful in both tailoring of novel amine-type drugs and early elimination of potentially toxic or mutagenic candidates. Also, chimeragenesis might serve in the construction of more efficient P450s for activation of amine drugs and/or bioremediation.

Structural motif-based homology modeling of CYP27A1 and site-directed mutational analyses affecting vitamin D hydroxylation

Human CYP27A1 is a mitochondrial cytochrome P450, which is principally found in the liver and plays important roles in the biological activation of vitamin D(3) and in the biosynthesis of bile acids. We have applied a systematic analysis of hydrogen bonding patterns in 11 prokaryotic and mammalian CYP crystal structures to construct a homology-based model of CYP27A1. Docking of vitamin D(3) structures into the active site of this model identified potential substrate contact residues in the F-helix, the beta-3 sheet, and the beta-5 sheet. Site-directed mutagenesis and expression in COS-1 cells confirmed that these positions affect enzymatic activity, in some cases shifting metabolism of 1alpha-hydroxyvitamin D(3) to favor 25- or 27-hydroxylation. The results suggest that conserved hydrophobic residues in the beta-5 hairpin help define the shape of the substrate binding cavity and that this structure interacts with Phe-248 in the F-helix. Mutations directed toward the beta-3a strand suggested a possible heme-binding interaction centered on Asn-403 and a structural role for substrate contact residues Thr-402 and Ser-404.

Three-dimensional model of the human aromatase enzyme and density functional parameterization of the iron-containing protoporphyrin IX for a molecular dynamics study of heme-cysteinato cytochromes

Mammalian cytochromes P450 (CYP) are enzymes of great biological and pharmaco-toxicological relevance. Due to their membrane-bound nature, the structural characterization of these proteins is extremely difficult, and therefore computational techniques, such as comparative modeling, may help obtaining reliable structures of members of this family. An important feature of CYP is the presence of an iron-containing porphyrin group at the enzyme active site. This calls for quantum chemical calculations to derive charges and parameters suitable for classical force field-based investigations of this proteins family. In this report, we first carried out density functional theory (DFT) computations to derive suitable charges for the Fe2+-containing heme group of P450 enzymes. Then, by means of the homology modeling technique, and taking advantage of the recently published crystal structure of the human CYP2C9, we built a new model of the human aromatase (CYP19) enzyme. Furthermore, to study the thermal stability of the new model as well as to test the suitability of the new DFT-based heme parameters, molecular dynamics (MD) simulations were carried out on both CYP2C9 and CYP19. Finally, the last few ns of aromatase MD trajectories were investigated following the essential dynamics protocol that allowed the detection of some correlated motions among some protein domains.

The Use of Genetic Markers to Determine Risk for Prostate Cancer at Prostate Biopsy

PURPOSE

We examined a panel of 13 polymorphisms in 13 different genes to determine whether specific genotypes can help predict prostate cancer at the time of biopsy among men prescreened with prostate-specific antigen and digital rectal exam.

EXPERIMENTAL DESIGN

We examined 2,088 consecutive men who were referred for prostate biopsy from 1997 to 2003. Thirteen genes were examined, including TNF308, GSTT1, KLK2, endostatin, MCRA, MCRV, tyrosinase, MSR1, CHK2, RNasel, HOGG1-326, HOGG1-11657, and HRAS1. Odds ratio for detection of prostate cancer were adjusted for age, race, prostate-specific antigen, digital rectal exam, family history of prostate cancer, and urinary symptoms.

RESULTS

Of the 2,088 men, 996 (47.7%) had cancer detected. Four genes (TNF308, GSTT1, KLK2, and HOGG1-326) were significantly associated with prostate cancer. The adjusted odds ratios (OR) for prostate cancer for patients with the AA genotype of the TNF308 gene was 1.92 [95% confidence interval (95% CI), 1.0-1.5, P = 0.05], compared with those with the GG genotype, and for patients with the TT genotype of the KLK2 gene, the OR was 1.5 (95% confidence interval, 1.0-2.2, P = 0.04), compared with the CC genotype. The OR for patients with a homozygous deletion of the GSTT1 gene was 0.81 (95% CI, 0.6-1.0, P = 0.06) compared with those with the deletion, and the OR for patients with the GG genotype of the HOGG1-326 gene was 0.68 (95% CI, 0.5-1.0, P = 0.05) compared with the CC genotype. Patients who had all four alleles that were positively associated with prostate cancer had an OR of 9.33 (95% CI, 2.4-35.8, P = 0.0005) for prostate cancer compared with patients with alleles that were negatively associated with prostate cancer.

CONCLUSIONS

Of the 13 polymorphisms, two were found to be positively associated with prostate cancer (TNF308 and KLK2) and two were negatively associated with prostate cancer (GSTT1 and HOGG1-326). Future studies are required to confirm these results.

Structural and functional analysis of a novel mutation of CYP21B in a heterozygote carrier of 21-hydroxylase deficiency

Congenital adrenal hyperplasia (CAH) due to 21-hydroxylase deficiency is one of the most common autosomal recessive disorders and occurs in its non-classical form in up to 6% of hirsute women. We report on a young woman with the clinical diagnosis of non-classical CAH and a novel, heterozygous missense mutation CTG-->GTG in exon 8, codon 317, of the steroid 21-hydroxylase CYP21B and complete loss of pseudogenes. Protein sequences of closely related P450 cytochromes and a homology-based 3D model of CYP21B were used for further functional analyses. We found that the mutated residue is part of a large cluster of hydrophobic residues. This cluster has three important features: (1) it is located directly next to the binding pocket, in close vicinity of the heme-cofactor, (2) all amino acids of the cluster are directly connected to two important binding regions, and (3) the packing within the cluster is very dense. Due to the tight packing in the cluster and its direct connection to the binding pocket region, any changes induced by the mutation of residue 317 can be expected to lead to structural shifts within the binding pocket and can explain the clinically observed impairment of 21-hydroxylase activity. In conclusion, the novel mutation L317V of the steroid 21-hydroxylase gene is associated with reduced steroid 21-hydroxylase activity probably due to structural shifts within the binding pocket and a mild phenotype of steroid 21-hydroxylase deficiency. In addition, the results support previous findings in which heterozygous CYP21 mutations are associated with symptoms of hyperandrogenism in susceptible individuals.

Association of Cytochrome P450 Enzymes is a Determining Factor in their Catalytic Activity

Previously, our laboratory demonstrated that one cytochrome P450 isoenzyme can influence the catalytic properties of another P450 isoenzyme when combined in a reconstituted system. Moreover, our data and that of other investigators indicate that P450 interaction is required for catalytic activity even when one isoenzyme is present. The goal of the current study was to examine the possible mechanism of these interactions in more detail. Analyzing recently published X-ray data of microsomal P450 enzymes and protein docking studies, four types of dimer formations of P450 enzymes were examined in more detail. In case of two dimer types, the aggregating partner was shown to contribute to NADPH cytochrome P450 reductase (CPR) binding-a flavoprotein whose interaction with P450 is required for expressing P450 functional activity of the neighboring P450 moiety. Thus, it was shown that dimerization of P450 enzymes might result in an altered affinity towards the CPR. Two dimer types were shown to exist only in the presence of a substrate, while the other two types exist also without a substrate present. The molecular basis was established for the fact that the presence of a substrate and other P450 enzymes simultaneously determine the catalytic activity. Furthermore, a kinetic model was improved describing the catalytic activity of P450 enzymes as a function of CPR concentration based on equilibrium between different supramolecular organizations of P450 enzymes. This model was successfully applied in order to explain our experimental data and that of other investigators.

Quantitative structure-activity relationships (QSARs) within the cytochrome P450 system: QSARs describing substrate binding, inhibition and induction of P450s

Quantitative structure-activity relationships (QSARs) within substrates, inducers and inhibitors of cytochromes P450 involved in xenobiotic metabolism are reported, together with QSARs associated with induction, inhibition and metabolic rate. The importance of frontier orbitals and shape descriptors, such as planarity (estimated by the area/depth(2) parameter) and rectangularity (estimated by the length/width parameter) is discussed, particularly in the context of the COMPACT system which discriminates between several P450 families associated with the activation and detoxication of xenobiotics. The use of parameters, particularly those derived from homology modelling of mammalian (especially human) P450s that are involved in exogenous metabolism, in generating QSARs for P450 substrates is discussed in the context of explaining differences in the binding affinities of human P450 substrates which are pharmacologically active.

The planetary biology of cytochrome P450 aromatases

Background

Joining a model for the molecular evolution of a protein family to the paleontological and geological records (geobiology), and then to the chemical structures of substrates, products, and protein folds, is emerging as a broad strategy for generating hypotheses concerning function in a post-genomic world. This strategy expands systems biology to a planetary context, necessary for a notion of fitness to underlie (as it must) any discussion of function within a biomolecular system.

Results

Here, we report an example of such an expansion, where tools from planetary biology were used to analyze three genes from the pig Sus scrofa that encode cytochrome P450 aromatases–enzymes that convert androgens into estrogens. The evolutionary history of the vertebrate aromatase gene family was reconstructed. Transition redundant exchange silent substitution metrics were used to interpolate dates for the divergence of family members, the paleontological record was consulted to identify changes in physiology that correlated in time with the change in molecular behavior, and new aromatase sequences from peccary were obtained. Metrics that detect changing function in proteins were then applied, including K_A/K_S values and those that exploit structural biology. These identified specific amino acid replacements that were associated with changing substrate and product specificity during the time of presumed adaptive change. The combined analysis suggests that aromatase paralogs arose in pigs as a result of selection for Suoidea with larger litters than their ancestors, and permitted the Suoidea to survive the global climatic trauma that began in the Eocene.

Conclusions

This combination of bioinformatics analysis, molecular evolution, paleontology, cladistics, global climatology, structural biology, and organic chemistry serves as a paradigm in planetary biology. As the geological, paleontological, and genomic records improve, this approach should become widely useful to make systems biology statements about high-level function for biomolecular systems.

Comparative molecular field analysis and QSAR on substrates binding to cytochrome P450 2D6

In this study, we utilized comparative molecular field analysis (CoMFA) to gain a better understanding of the steric and electrostatic features of the cytochrome p450 2D6 (CYP2D6) active site. The training set consists of 24 substrates with reported K(M) values from liver microsomal CYP2D6 spanning an activity range of almost three log units. The low energy conformers were fit by root mean square (RMS) to minaprine at the site of metabolism and to the protonated nitrogen. In this manner, we constructed two CoMFA models, one model with a distance constraint and another without. The model with the distance parameter (non-cross-validated R(2)=0.99) was approximately equal to the CoMFA without a distance parameter (non-cross-validated R(2)=0.98). Validation of our CoMFA was accomplished by predicting the K(M) values of 15 diverse CYP2D6 substrates not in the original training set resulting in a predictive R(2)=0.62. Finally, we also pursued correlations of pK(a) and log P with CYP2D6 substrate K(M) in an effort to investigate other physicochemical properties.

Functional interaction of cytochrome P450 with its redox partners: a critical assessment and update of the topology of predicted contact regions

The problem of donor-acceptor recognition has been the most important and intriguing one in the area of P450 research. The present review outlines the topological background of electron-transfer complex formation, showing that the progress in collaborative investigations, combining physical techniques with chemical-modification and immunolocalization studies as well as site-directed mutagenesis experiments, has increasingly enabled the substantiation of hypothetical work resulting from homology modelling of P450s. Circumstantial analysis reveals the contact regions for redox proteins to cluster on the proximal face of P450s, constituting parts of the highly conserved, heme-binding core fold. However, more variable structural components located in the periphery of the hemoprotein molecules also participate in donor docking. The cross-reactivity of electron carriers, purified from pro- and eukaryotic sources, with a diversity of P450 species points at a possible evolutionary conservation of common anchoring domains. While electrostatic mechanisms appear to dominate orientation toward each other of the redox partners to generate pre-collisional encounter complexes, hydrophobic forces are likely to foster electron transfer events by through-bonding or pi-stacking interactions. Moreover, electron-tunneling pathways seem to be operative as well. The availability of new P450 crystal structures together with improved validation strategies will undoubtedly permit the production of increasingly satisfactory three-dimensional donor-acceptor models serving to better understand the molecular principles governing functional association of the redox proteins.

A compound heterozygote case of type II aldosterone synthase deficiency

An infant with failure to thrive, persistent hyponatremia and episodic vomiting and diarrhea was admitted to hospital at 9 months of age, and the diagnosis of type II aldosterone synthase deficiency was confirmed by plasma and urinary steroid determinations. The entire coding sequence of the aldosterone synthase gene (CYP11B2) was determined (both strands) in the affected infant, an unaffected sibling, and both parents. An exon 3 mutation (C554T, leading to amino acid T185I) was found in the father and both siblings, and an exon 9 mutation (A1492G, leading to T498A) was found in the affected infant and the mother. Expression of the mutant sequences in COS cells showed steroidogenic patterns typical of aldosterone synthase type II deficiency, including very low levels of aldosterone synthesis (< or =0.5% of wild-type enzyme) consistent with the low aldosterone levels in the patient's plasma. Both mutations in this compound heterozygote localize to the beta 3-sheet in the cytochrome P450 enzyme structure, as does the previously characterized R181W mutation. This region of the enzyme is not part of the putative structural core, but mutations to this region suggest that it is important for conferring the unique ability of aldosterone synthase to catalyze efficient oxygenation of the C(18) carbon of steroid substrates.

In vitro and in vivo models for the evaluation of potent inhibitors of male rat 17alpha-hydroxylase/C17,20-lyase

The C(17,20)-lyase is a key enzyme in the biosynthesis of androgens by both the testes and adrenals. A complete inhibition of this enzyme would provide an alternative means of androgen suppression for the treatment of prostatic cancers. In the present study, the inhibitory effects of new non-steroidal compounds were tested in vitro on rat C(17,20)-lyase versus abiraterone, a reference steroidal inhibitor. Their activities were also evaluated in vivo on plasma testosterone (T) and luteinizing hormone (LH) levels and on testes, adrenals, seminal vesicles (SV) and ventral prostate (VP) weights after 3 days of oral treatment to adult male rats (50mg/kg per day p.o.). Inhibition in the nanomolar range was obtained with TX 977, the lead racemate product in this series, and optimization is ongoing based on a slight dissociation observed between its two diastereoisomers, TX 1196-11 (S) and TX 1197-11 (R). These non-steroidal compounds (including YM 55208, a reference competitor) proved to be more active in vivo than abiraterone acetate in this model, but the observed impact on adrenal weight suggests that the specificity of lyase inhibition versus corticosteroid biosynthesis deserves further investigations with this new class of potentially useful agents for the treatment of androgen-dependent prostate cancer.

Membrane reconstitution of recombinant guinea pig cytochrome P45017alpha and the effects of site-directed mutagenesis on androgen formation

Although cytochrome P45017alpha catalyzes the formation of androgen from both pregnenolone and progesterone, the production of androstenedione from progesterone is a major pathway in the guinea pig, rat, mouse, and hamster. In contrast, human, bovine and sheep P45017alpha produce dehydroepiandrosterone from pregnenolone. Cytochrome P45017alphas from all of these animals have high homology in the amino acid sequence around the 340-370 region. To investigate the substrate preferences for androgen production, we replaced a few amino acids in the 340-370 region of guinea pig P45017alpha with those found in the other animals. The recombinant P45017alphas were expressed in E. coli DH5alpha, purified by column chromatography and incorporated into liposome membranes. The (His)(4) tag in the recombinant P45017alphas had little effect on the interaction with NADPH-P450 reductase in the membranes. The recombinant P45017alphas with a single-position mutation of F344I, H349R or M352L and with double-position mutations of F344I and H349R and triple-position mutations showed decreases in the production of 17alpha-hydroxypregnenolone, androstenedione and dehydroepiandrosterone. The activity for 17alpha-hydroxyprogeterone production was increased significantly by the F344I mutation. The addition of cytochrome b5 did not have much of an effect on the 17alpha-hydroxylation but had a significant effect on androgen production in both the nonmutated and mutated P45017alphas.

A comparative approach to structure-function studies of mammalian aromatases

To date, structure--function studies of aromatase cytochrome P450 (P450arom) have been advanced by point mutation analyses utilizing almost exclusively the human enzyme, in conjunction with computer-generated models of the three-dimensional form of the enzyme based on prokaryotic cytochromes P450. Recent studies have identified duplicated isozymes of porcine P450arom, the gonadal and placental forms of which appear to differ substantially in substrate utilization and inhibitor sensitivity. We present a comparative approach to define regions of P450arom responsible for specific functional characteristics using complimentary DNAs encoding the porcine isozymes. Constructs encoding the native and chimeric porcine and human P450arom enzymes were transiently expressed and activity was assessed using the tritiated water assay. Sensitivity to inhibition by the imidazole etomidate was investigated, and P450arom expression was assessed by immunoblot analysis. All constructs yielded active P450arom, suggesting that exchanging entire structural elements does not preclude catalytic function. The activity of the gonadal isozyme was shown to be inhibited by etomidate at concentrations 185 and 300-fold lower than those required to induce a similar inhibition of the placental and human enzymes, respectively. In contrast, there was only a two-fold difference in the sensitivity of the gonadal and placental isozymes to inhibition by CGS16949A. Analysis of chimeric constructs indicated that the sensitivity to etomidate was associated with residues in the B, B' and C helices of the gonadal P450arom encompassing only one of six putative substrate recognition sites. Additionally, sensitivity to etomidate was not correlated with enzyme activity among the chimeric enzymes. Therefore, it appears that residues of the porcine gonadal P450arom that are responsible for etomidate binding may be distinct from those involved in substrate recognition and metabolism. These data support the notion that a comparative approach employing the use of chimeric enzymes provides a useful tool in directing point mutational analysis to determine residues important in inhibitor and perhaps substrate recognition of P450 enzymes such as P450arom. These studies are currently in progress.

Interactions between redox partners in various cytochrome P450 systems: functional and structural aspects

The various types of redox partner interactions employed in cytochrome P450 systems are described. The similarities and differences between the redox components in the major categories of P450 systems present in bacteria, mitochondria and microsomes are discussed in the light of the accumulated evidence from X-ray crystallographic and NMR spectroscopic determinations. Molecular modeling of the interactions between the redox components in various P450 mono-oxygenase systems is proposed on the basis of structural and mutagenesis information, together with experimental findings based on chemical modification of key residues likely to be associated with complementary binding sites on certain typical P450 isoforms and their respective redox partners.

A polymorphism in the CYP17 gene is associated with prostate cancer risk

CYP17 encodes the enzyme cytochrome P-450c17 alpha, which mediates both 17 alpha-hydroxylase and 17,20-lyase in the steroid biosynthesis pathway. A polymorphism in the 5; promoter region of the CYP17 gene has been described. Steroid hormones, especially androgens, are believed to play a key role in the etiology of prostate cancer. Therefore, polymorphisms in genes involved in the androgen metabolism may affect the risk of prostate cancer. We conducted a case-control study of 63 patients with untreated histologically proven prostate cancer and 126 age-matched control men with benign prostatic hyperplasia (BPH) to determine whether a polymorphism in the CYP17 gene is associated with prostate cancer risk. This polymorphism was investigated by PCR/RFLP using DNA from lymphocytes. The transition (T-->C) in the risk allele (A2) creates a new recognition site for the restriction enzyme MspAI, which permits designation of the wildtype (A1) and the risk allele (A2). The prevalence of the A2/A2 genotype was significantly higher (P = 0.03) in the cancer group (23.8%) than in the BPH control group (9.5%). We found an increased risk in men carrying 2 A2 alleles (OR = 2.80, 95%CI = 1.02-77.76). For carrier with at least 1 A2 allele, the OR was 0.90 (95%CI = 0.43-1.89). After stratification by median age (66 years) at time of diagnosis, a marked increased risk was found in carriers of the A2/A2 genotype older than 66 years (OR = 8.93, 95%CI = 1.78-49.19, P = 0.01). Although the sample size is rather small and the controls are BPH patients, our results suggest that the CYP17A2/A2 genotype may be a biomarker for prostate cancer risk, especially for older men.

A three-dimensional model of lanosterol 14alpha-demethylase of Candida albicans and its interaction with azole antifungals

The three-dimensional structure of lanosterol 14alpha-demethylase (P450(14DM), CYP51) of Candida albicans was modeled on the basis of crystallographic coordinates of four prokaryotic P450s: P450BM3, P450cam, P450terp, and P450eryF. The P450(14DM) sequence was aligned to those of known proteins using a knowledge-based alignment method. The main chain coordinates of the core regions were transferred directly from the corresponding coordinates of P450BM3. The side chain conformations of the core regions were determined by the conformations of the equivalent residues with the highest homologous scores in four crystal structures. The model was then refined using molecular mechanics and molecular dynamics. The reliability of the resulting model was assessed by Ramachandran plots, Profile-3D, hydropathy plot analysis, and by analyzing the consistency of the model with the experimental data. The structurally and functionally important residues such as the heme binding residues, the residues interacting with redox-partner protein and/or involved in electron transfer, the residues lining substrate access channel, and the substrate binding residues were identified from the model. These residues are candidates for further site-directed mutagenesis and site-specific antipeptide antibody binding experiments. The active analogue approach was employed to search the pharmacophoric conformations for 14 azole antifungals. The resulting bioactive conformations were docked into the active site of lanosterol 14alpha-demethylase of Candida albicans. All 14 azole antifungals are shown to have a similar docking mode in the active site. The halogenated phenyl group of azole inhibitors is deep in the same hydrophobic binding cleft as the 17-alkyl chain of substrate. The pi-pi stacking interaction might exist between halogenated phenyl ring of inhibitors and the aromatic ring of residue Y132. The long side chains of some inhibitors such as itraconazole and ketoconazole surpass the active site and interact with the residues in the substrate access channel. To compare with mammalian enzymes, structurally selective residues of the active site of fungal lanosterol 14alpha-demethylase are distributed in the C terminus of F helix, beta6-1 sheet and beta6-2 sheet.

Congenital adrenal hyperplasia due to 21-hydroxylase deficiency

More than 90% of cases of congenital adrenal hyperplasia (CAH, the inherited inability to synthesize cortisol) are caused by 21-hydroxylase deficiency. Females with severe, classic 21-hydroxylase deficiency are exposed to excess androgens prenatally and are born with virilized external genitalia. Most patients cannot synthesize sufficient aldosterone to maintain sodium balance and may develop potentially fatal "salt wasting" crises if not treated. The disease is caused by mutations in the CYP21 gene encoding the steroid 21-hydroxylase enzyme. More than 90% of these mutations result from intergenic recombinations between CYP21 and the closely linked CYP21P pseudogene. Approximately 20% are gene deletions due to unequal crossing over during meiosis, whereas the remainder are gene conversions--transfers to CYP21 of deleterious mutations normally present in CYP21P. The degree to which each mutation compromises enzymatic activity is strongly correlated with the clinical severity of the disease in patients carrying it. Prenatal diagnosis by direct mutation detection permits prenatal treatment of affected females to minimize genital virilization. Neonatal screening by hormonal methods identifies affected children before salt wasting crises develop, reducing mortality from this condition. Glucocorticoid and mineralocorticoid replacement are the mainstays of treatment, but more rational dosing and additional therapies are being developed.

Mechanism-based inactivation of cytochrome P450 2B1 by 7-ethynylcoumarin: verification of apo-P450 adduction by electrospray ion trap mass spectrometry

7-Ethynylcoumarin was synthesized as a potential mechanism-based inhibitor, and it was found to be an effective inactivator of 7-ethoxy-4-(trifluoromethyl)coumarin (7EFC) O-deethylation catalyzed by purified, reconstituted P450 2B1. In contrast, 7-ethynylcoumarin demonstrated minimal inactivation of P450 2A6-mediated 7-hydroxycoumarin formation. The inactivation of P450 2B1 demonstrated pseudo-first-order kinetics and was NADPH- and inhibitor-dependent. The maximal rate constant for the inactivation of 2B1 was 0.39 min(-)(1) at 30 degrees C, and thus, the time required to inactivate 50% of the P450 2B1 that was present (t(1/2)) was 1.8 min. The estimated concentration which led to half-maximal inactivation (K(I)) was 25 microM. No protection from inactivation was seen in the presence of nucleophiles (glutathione and sodium cyanide), an iron chelator (deferroxamine), or superoxide dismutase and catalase. Addition of the substrate (7EFC) protected P450 2B1 from inactivation, in a concentration-dependent manner. The partition ratio for P450 2B1 was 25; i.e., the number of metabolic events was 25-fold higher than the number of inactivating events. Incubations of 7-ethynylcoumarin with P450 2B1 for 10 min resulted in an 80% loss in enzymatic activity, while 90% of the ability to form a reduced-CO complex remained. This activity loss was not recovered following dialysis, indicative of irreversible inactivation. Covalent attachment of the entire inhibitor and oxygen to apo-P450 2B1, in a 1:1 ratio, was shown via electrospray ion trap mass spectrometry. This method also verified the absence of modification to the heme or the cytochrome P450 reductase. Taken together, the characterization of the inhibition seen with P450 2B1 and 7-ethynylcoumarin was consistent with all of the criteria required to distinguish a mechanism-based inactivator. In addition, electrospray ion trap mass spectrometry has the potential to be applied to protein adducts above and beyond those associated with the mechanism-based inactivation of cytochrome P450s.

Molecular modelling of lanosterol 14 alpha-demethylase (CYP51) from Saccharomyces cerevisiae via homology with CYP102, a unique bacterial cytochrome P450 isoform: quantitative structure-activity relationships (QSARs) within two related series of antifungal azole derivatives

The construction of a three-dimensional molecular model of the fungal form of cytochrome P450 (CYP51) from Saccharomyces cerevisiae, based on homology with the haemoprotein domain of CYP102 from Bacillus megaterium (a unique bacterial P450 of known crystal structure) is described. It is found that the endogenous substrate, lanosterol, can readily occupy the putative active site of the CYP51 model such that the known mono-oxygenation reaction, leading to C14-demethylation of lanosterol, is the preferred route of metabolism for this particular substrate. Key amino acid contacts within the CYP51 active site appear to orientate lanosterol for oxidative attack at the C14-methyl group, and the position of the substrate relative to the haem moiety is consistent with the phenyl-iron complexation studies reported by Tuck et al. [J. Biol. Chem., 267, 13175-13179 (1992)]. Typical azole inhibitors, such as ketoconazole, are able to fit the putative active site of CYP51 by a combination of haem ligation, hydrogen bonding, pi-pi stacking and hydrophobic interactions within the enzyme's haem environment. The mode of action of azole antifungals, as described by the modelling studies, is supported by quantitative structure-activity relationship (QSAR) analyses on two groups of structurally related fungal inhibitors. Moreover, the results of molecular electrostatic isopotential (EIP) energy calculations are compatible with the proposed mode of binding between azole antifungal agents and the putative active site of CYP51, although membrane interactions may also have a role in the antifungal activity of azole derivatives.

Evolution of the cytochrome P450 superfamily: sequence alignments and pharmacogenetics

The evolution of the cytochrome P450 (CYP) superfamily is described, with particular reference to major events in the development of biological forms during geological time. It is noted that the currently accepted timescale for the elaboration of the P450 phylogenetic tree exhibits close parallels with the evolution of terrestrial biota. Indeed, the present human P450 complement of xenobiotic-metabolizing enzymes may have originated from coevolutionary 'warfare' between plants and animals during the Devonian period about 400 million years ago. A number of key correspondences between the evolution of P450 system and the course of biological development over time, point to a mechanistic molecular biology of evolution which is consistent with a steady increase in atmospheric oxygenation beginning over 2000 million years ago, whereas dietary changes during more recent geological time may provide one possible explanation for certain species differences in metabolism. Alignment between P450 protein sequences within the same family or subfamily, together with across-family comparisons, aid the rationalization of drug metabolism specificities for different P450 isoforms, and can assist in an understanding of genetic polymorphisms in P450-mediated oxidations at the molecular level. Moreover, the variation in P450 regulatory mechanisms and inducibilities between different mammalian species are likely to have important implications for current procedures of chemical safety evaluation, which rely on pure genetic strains of laboratory bred rodents for the testing of compounds destined for human exposure.