Homology modelling of human cytochromes P450 involved in xenobiotic metabolism and rationalization of substrate selectivity

Multidisciplinary Insights into the Structure-Function Relationship of the CYP2B6 Active Site

Cytochrome P450 2B6 (CYP2B6) is a highly polymorphic human enzyme involved in the metabolism of many clinically relevant drugs, environmental toxins, and endogenous molecules with disparate structures. Over the last 20-plus years, in silico and in vitro studies of CYP2B6 using various ligands have provided foundational information regarding the substrate specificity and structure-function relationship of this enzyme. Approaches such as homology modeling, X-ray crystallography, molecular docking, and kinetic activity assays coupled with CYP2B6 mutagenesis have done much to characterize this originally neglected monooxygenase. However, a complete understanding of the structural details that make new chemical entities substrates of CYP2B6 is still lacking. Surprisingly little in vitro data has been obtained about the structure-function relationship of amino acids identified to be in the CYP2B6 active site. Since much attention has already been devoted to elucidating the function of CYP2B6 allelic variants, here we review the salient findings of in silico and in vitro studies of the CYP2B6 structure-function relationship with a deliberate focus on the active site. In addition to summarizing these complementary approaches to studying structure-function relationships, we note gaps/challenges in existing data such as the need for more CYP2B6 crystal structures, molecular docking results with various ligands, and data coupling CYP2B6 active site mutagenesis with kinetic parameter measurement under standard expression conditions. Harnessing in silico and in vitro techniques in tandem to understand the CYP2B6 structure-function relationship will likely offer further insights into CYP2B6-mediated metabolism. SIGNIFICANCE STATEMENT: The apparent importance of cytochrome P450 2B6 (CYP2B6) in the metabolism of various xenobiotics and endogenous molecules has grown since its discovery with many in silico and in vitro studies offering a partial description of its structure-function relationship. Determining the structure-function relationship of CYP2B6 is difficult but may be aided by thorough biochemical investigations of the CYP2B6 active site that provide a more complete pharmacological understanding of this important enzyme.

The Development and Validation of a Novel "Dual Cocktail" Probe for Cytochrome P450s and Transporter Functions to Evaluate Pharmacokinetic Drug-Drug and Herb-Drug Interactions

This study was designed to develop and validate a 10 probe drug cocktail named "Dual Cocktail", composed of caffeine (Cyp1a2 in rat and CYP1A2 in human, 1 mg/kg), diclofenac (Cyp2c11 in rat and CYP2C9 in human, 2 mg/kg), omeprazole (Cyp2c11 in rat and CYP2C19 in human, 2 mg/kg), dextromethorphan (Cyp2d2 in rat and CYP2D6 in human, 10 mg/kg), nifedipine (Cyp3a1 in rat and CYP3A4 in human, 0.5 mg/kg), metformin (Oct1/2 in rat and OCT1/2 in human, 0.5 mg/kg), furosemide (Oat1/3 in rat and OAT1/3 in human, 0.1 mg/kg), valsartan (Oatp2 in rat and OATP1B1/1B3 in human, 0.2 mg/kg), digoxin (P-gp in rat and human, 2 mg/kg), and methotrexate (Mrp2 in rat and MRP2 in human, 0.5 mg/kg), for the evaluation of pharmacokinetic drug-drug and herb-drug interactions through the modulation of a representative panel of CYP enzymes or transporters in rats. To ensure no interaction among the ten probe substrates, we developed a 2-step evaluation protocol. In the first step, the pharmacokinetic properties of five individual CYP probe substrates and five individual transporter substrates were compared with the pharmacokinetics of five CYP cocktail or five transporters cocktails in two groups of randomly assigned rats. Next, a pharmacokinetic comparison was conducted between the CYP or transporter cocktail group and the dual cocktail group, respectively. None of the ten comparison groups was found to be statistically significant, indicating the CYP and transporter substrate sets or dual cocktail set could be concomitantly administered in rats. The "Dual Cocktail" was further validated by assessing the metabolism of nifedipine and omeprazole, which was significantly reduced by a single oral dose of ketoconazole (10 mg/kg); however, no changes were observed in the pharmacokinetic parameters of other probe substrates. Additionally, multiple oral doses of rifampin (20 mg/kg) reduced the plasma concentrations of nifedipine and digoxin, although not any of the other substrates. In conclusion, the dual cocktail can be used to characterize potential pharmacokinetic drug-drug interactions by simultaneously monitoring the activity of multiple CYP isoforms and transporters.

Structure, function, regulation and polymorphism and the clinical significance of human cytochrome P450 1A2

Human CYP1A2 is one of the major CYPs in human liver and metabolizes a number of clinical drugs (e.g., clozapine, tacrine, tizanidine, and theophylline; n > 110), a number of procarcinogens (e.g., benzo[a]pyrene and aromatic amines), and several important endogenous compounds (e.g., steroids). CYP1A2 is subject to reversible and/or irreversible inhibition by a number of drugs, natural substances, and other compounds. The CYP1A gene cluster has been mapped on to chromosome 15q24.1, with close link between CYP1A1 and 1A2 sharing a common 5'-flanking region. The human CYP1A2 gene spans almost 7.8 kb comprising seven exons and six introns and codes a 515-residue protein with a molecular mass of 58,294 Da. The recently resolved CYP1A2 structure has a relatively compact, planar active site cavity that is highly adapted for the size and shape of its substrates. The architecture of the active site of 1A2 is characterized by multiple residues on helices F and I that constitutes two parallel substrate binding platforms on either side of the cavity. A large interindividual variability in the expression and activity of CYP1A2 has been observed, which is largely caused by genetic, epigenetic and environmental factors (e.g., smoking). CYP1A2 is primarily regulated by the aromatic hydrocarbon receptor (AhR) and CYP1A2 is induced through AhR-mediated transactivation following ligand binding and nuclear translocation. Induction or inhibition of CYP1A2 may provide partial explanation for some clinical drug interactions. To date, more than 15 variant alleles and a series of subvariants of the CYP1A2 gene have been identified and some of them have been associated with altered drug clearance and response and disease susceptibility. Further studies are warranted to explore the clinical and toxicological significance of altered CYP1A2 expression and activity caused by genetic, epigenetic, and environmental factors.

New insights into the structural characteristics and functional relevance of the human cytochrome P450 2D6 enzyme

To date, the crystal structures of at least 12 human CYPs (1A2, 2A6, 2A13, 2C8, 2C9, 2D6, 2E1, 2R1, 3A4, 7A1, 8A1, and 46A1) have been determined. CYP2D6 accounts for only a small percentage of all hepatic CYPs (< 2%), but it metabolizes approximately 25% of clinically used drugs with significant polymorphisms. CYP2D6 also metabolizes procarcinogens and neurotoxins, such as 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine, 1,2,3,4-tetrahydroquinoline, and indolealkylamines. Moreover, the enzyme utilizes hydroxytryptamines and neurosteroids as endogenous substrates. Typical CYP2D6 substrates are usually lipophilic bases with an aromatic ring and a nitrogen atom, which can be protonated at physiological pH. Substrate binding is generally followed by oxidation (5-7 A) from the proposed nitrogen-Asp301 interaction. A number of homology models have been constructed to explore the structural features of CYP2D6, while antibody studies also provide useful structural information. Site-directed mutagenesis studies have demonstrated that Glu216, Asp301, Phe120, Phe481, and Phe483 play important roles in determining the binding of ligands to CYP2D6. The structure of human CYP2D6 has been recently determined and shows the characteristic CYP fold observed for other members of the CYP superfamily. The lengths and orientations of the individual secondary structural elements in the CYP2D6 structure are similar to those seen in other human CYP2 members, such as CYP2C9 and 2C8. The 2D6 structure has a well-defined active-site cavity located above the heme group with a volume of approximately 540 A(3), which is larger than equivalent cavities in CYP2A6 (260 A(3)), 1A2 (375 A(3)), and 2E1 (190 A(3)), but smaller than those in CYP3A4 (1385 A(3)) and 2C8 (1438 A(3)). Further studies are required to delineate the molecular mechanisms involved in CYP2D6 ligand interactions and their implications for drug development and clinical practice.

Insights into drug metabolism by cytochromes P450 from modelling studies of CYP2D6-drug interactions

The cytochromes P450 (CYPs) comprise a vast superfamily of enzymes found in virtually all life forms. In mammals, xenobiotic metabolizing CYPs provide crucial protection from the effects of exposure to a wide variety of chemicals, including environmental toxins and therapeutic drugs. Ideally, the information on the possible metabolism by CYPs required during drug development would be obtained from crystal structures of all the CYPs of interest. For some years only crystal structures of distantly related bacterial CYPs were available and homology modelling techniques were used to bridge the gap and produce structural models of human CYPs, and thereby obtain useful functional information. A significant step forward in the reliability of these models came seven years ago with the first crystal structure of a mammalian CYP, rabbit CYP2C5, followed by the structures of six human enzymes, CYP1A2, CYP2A6, CYP2C8, CYP2C9, CYP2D6 and CYP3A4, and a second rabbit enzyme, CYP2B4. In this review we describe as a case study the evolution of a CYP2D6 model, leading to the validation of the model as an in silico tool for predicting binding and metabolism. This work has led directly to the successful design of CYP2D6 mutants with novel activity-including creating a testosterone hydroxylase, converting quinidine from inhibitor to substrate, creating a diclofenac hydroxylase and creating a dextromethorphan O-demethylase. Our modelling-derived hypothesis-driven integrated interdisciplinary studies have given key insight into the molecular determinants of CYP2D6 and other important drug metabolizing enzymes.

The human cytochrome P450 sub-family: Transcriptional regulation, inter-individual variation and interaction networks

The Cytochrome P450 super-family is a fundamental requirement for the viability of most life, with Cytochrome P450 proteins having been identified in organisms ranging from bacteria to man. These enzymes may be subdivided into those that metabolise purely endogenous chemicals, and those that are involved in xenobiotic metabolism. Of the latter group it can be argued that CYP3A sub-family members rank as the most important; their high expression in the liver and wide substrate specificity mean that they are clinically important in the metabolism of many therapeutic drugs, and alteration in their activity is central to many clinically-relevant drug-drug interactions. In this review I will examine the human CYP3A enzymes, discussing their genome structure, common allelic variants and, in greatest detail, their transcriptional regulation. Through examination of these characteristics we will see both striking similarities and differences between the four human CYP3A enzymes, which may have important impacts on inter-individual response to chemical exposure. Finally, the role of nuclear receptors in regulating CYP3A gene expression, and indeed that of many other proteins involved in drug metabolism, will be examined: Such an examination will show the need to utilize a systems biology approach to understand fully how the human body responds to chemical exposure.

Inter-individual variation of several cytochrome P450 2D6 splice variants in human liver

To examine the possibility that inter-individual differences in splicing partially explain the observed differences in CYP2D6 activity, we amplified its full-length cDNA in 96 human liver RNA samples and discovered five splice variants: intron 5 retention, intron 6 retention, intron 5 and intron 6 double retention, exon 3 skipping, and partial intron 1 retention. All of the CYP2D6 splice variants we identified are probably nonfunctional transcripts. Substantial inter-individual variation in the proportions of the CYP2D6 transcript represented by splice variants, measured by real-time PCR, suggests that the presence of these splice variants contributes to the population variation in CYP2D6 activity. Relatively high levels of intron 6 retention were not correlated with the newly discovered single nucleotide polymorphism 2988G > A in intron 6 (CYP2D6*41) but did correlate with the more common CYP2D6*34 allele. Our study prompts further investigations to explore the effect of these splice variants on drug metabolism.

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.

A reaction network model for CYP2E1-mediated metabolism of toxicant mixtures

In this paper, we describe a modeling approach to predict the interlinked pathways and kinetics resulting from CYP2E1-mediated metabolism of both pure species and chemical mixtures. This approach is based on the concept of chemical reaction networks, an idea that has formed the basis for simulation tools that have shown good predictive capabilities in the petroleum industry, but also an idea that has heretofore seen minimal application in the biomedical research arena. Although the initial target for developing this reaction network approach was cytochrome P450 2E1 (CYP2E1) and its over 200 substrates, this technology has been used for other families of CYP enzymes and their substrates in our laboratory. Utilizing this approach, we have produced a modular 'predictive metabolomics' simulation framework comprising interdependent software components that perform such tasks as testing of substrate binding feasibility, performing virtual chemistry, formulating reaction-rate equations, computing reaction kinetics and predicting time-dependent species concentrations. As an illustrative example, we outline the application of this framework to the prediction of the reaction networks resulting from the Phase I metabolism of two compounds of important toxicological interest. The potential of this modeling technology is immense in providing a computer simulation platform for complex-chemical mixtures and complex-biological systems. It is possible that this technology will play an important role in formulating a 'Virtual Human'.

Molecular modelling of CYP2A enzymes: application to metabolism of the tobacco-specific nitrosamine 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK)

1. Tobacco-specific nitrosamine 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) is a lung carcinogen in a variety of animal models and a putative human lung carcinogen. Its tumorigenic potential is unmasked via cytochrome P450 (CYP)-mediated hydroxylation of the carbon atoms adjacent to the nitroso moiety (i.e. alpha-hydroxylation). Therefore, elucidation of enzyme-substrate interactions that facilitate alpha-hydroxylation is important to gain insight into the tumorigenic mechanism of NNK and to develop potent inhibitors of this detrimental reaction. 2. Molecular models of CYP2A enzymes from mice, rats and humans that are catalysts of NNK bioactivation were constructed and used, in conjunction with docking experiments, to identify active-site residues that make important substrate contacts. 3. Docking studies revealed that hydrophobic residues at positions 117, 209, 365 and 481, among others, play critical roles in orienting NNK in the active site to effect alpha-hydroxylation. These molecular models were then used to rationalize the stereo- and regioselectivity, as well as the efficiency, of CYP2A-mediated NNK metabolism.

Molecular modeling of cytochrome P450 1A1: enzyme-substrate interactions and substrate binding affinities

Human cytochrome P450 1A1, which is present in lungs, plays an important role in the metabolic activation of chemical carcinogens, and in particular, is thought to be linked to lung cancer. The mechanism of carcinogenesis is related to the enzyme's ability to oxidize highly toxic compounds, such as polycyclic aromatic hydrocarbons (PAHs), to their carcinogenic derivatives. In order to better understand P450 1A1 function, a homology model of this enzyme has been constructed. The model has been based on the structure of P450 2C5, the first mammalian P450 to be crystallized. The coordinates of the model have been calculated using a consensus strategy, and the resulting structure has been evaluated with the ProStat and Profiles-3D programs. P450 1A1 substrates, such as benzo[a]pyrene, ethoxyresorufin and methoxyresorufin, were then docked into the active site of the model, and key amino acid residues able to interact with the substrate, have been identified. The analysis of enzyme-substrate interactions indicated that hydrophobic interactions are mainly responsible for binding of these substrates in the active site. Moreover, the non-bond enzyme-substrate interaction energy for ethoxyresorufin was lower than that for methoxyresorufin, which is consistent with higher activity of 1A1 towards the former substrate. Key residue Val-382 may play an important role in these interactions. Additionally, we performed binding free energy calculations for the three substrates. The obtained values were similar to those observed experimentally, which suggests that this approach might be useful for prediction of binding constants.

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.

Arginine to lysine 108 substitution in recombinant CYP1A2 abolishes methoxyresorufin metabolism in lymphoblastoid cells

1. Cytochrome P4501A2 (CYP1A2) activates a large number of procarcinogens to carcinogens. Phytochemicals such as flavones can inhibit CYP1A2 activity competitively, and hydroxylated derivatives of flavone (galangin) may be potent, selective inhibitors of CYP1A2 activity relative to CYP1A1 activity. Molecular modelling of the CYP1A2 interaction with hydroxylated derivatives of flavone suggests that a number of hydrophobic residues of the substrate-binding domain engage in hydrogen bonding with such inhibitors. 2. We have tested this model using site-directed mutagenesis of these residues in expression plasmids transfected into the human B-lymphoblastoid cell line, AHH-1 TK+/-. 3. Consistent with the molecular model's predicted placement in the active site, amino acid substitutions at the predicted residues abolished CYP1A2 enzymatic activity. 4. Transfected cell lines contained equal amounts of immunoreactive CYP1A2. 5. Our results support the molecular model's prediction of the critical amino acid residues present in the hydrophobic active site, residues that can hydrogen bond with CYP1A2 inhibitors and modify substrate binding and/or turnover.

Cloning, sequencing, and characterization of CYP1A1 cDNA from leaping mullet (Liza Saliens) liver and implications for the potential functions of its conserved amino acids

A 2,037 bp CYP1A1 cDNA (GenBank AF072899) was cloned through screening of a lambdaZipLox cDNA library constructed from the liver of a leaping mullet (Liza saliens) fish captured from Izmir Bay on the Aegean coast of Turkey using rainbow trout CYP1A1 cDNA as a probe. This clone has a 130 bp 5'-flanking region, a 1,563 bp open reading frame (ORF) encoding a 521-amino acid protein (58,972 Da), and a 344 bp 3'-untranslated region without a poly (A) tail. Alignment of the deduced amino acids of CYP1A1 cDNAs showed 58% and 69-96% identities with human and 12 other fish species, respectively. Southern blot analysis suggested that this CYP1A1 cDNA was from a single-copy gene. Based on the comparison with CYP1A1 genes reported for fish and mammals, the leaping mullet CYP1A1 gene is probably split into 7 exons. The intron insertion sites were predicted. Alignment of the CYP1A1 cDNA encoded amino acids from 13 fish and 7 mammalian species disclosed differences in highly conserved amino acids between aquatic and land vertebrates. The possible associated secondary structure; conserved motifs and substrate-binding sites were discussed. The phylogenetic relationships of CYP1A1s among 13 fish species were analyzed by a distance method.

Comparative modelling of cytochromes P450

The superfamily of enzymes known as the cytochromes P450 (P450s) comprises a wide-ranging class of proteins with diverse functions. They are known, amongst other things, to be involved in the hormonal regulation of metabolism and reproduction, as well as having a major clinical significance through their association with diseases such as cancer, diabetes and hepatitis. Knowledge of the three-dimensional (3D) structure of a protein gives insight into its function. The 3D structures of P450s are therefore of considerable scientific interest. A number of high-resolution structures of P450s have been determined by X-ray crystallography and studies of these structures have provided valuable insights into the mechanism of these enzymes. Only one of these structures is mammalian and as yet there is no structural information on human P450s in the public domain. Until such a structure is solved it is necessary to employ alternative methods to gain structural insight into how human P450s perform their biological function. Here we report on the use of comparative modelling to predict the structure of human P450s based on knowledge of their amino acid sequences plus the 3D structures of other (not human) P450s. As an illustrative example of these techniques we have modelled the structure of P450 2C5 using five bacterial P450 structures as templates. We examine the importance of selecting suitable templates, obtaining a good amino acid sequence alignment, and evaluating the models generated. To improve the quality of the models an iterative cycle of sequence alignment, model building, and model evaluation is employed. The result is a model with excellent stereochemistry, good amino acid side chain environment properties, and a Calpha trace similar to the crystal structure.

Molecular modeling of human cytochrome P450-substrate interactions

The results of homology modeling of 10 human cytochrome P450 (CYP) enzymes involved in the Phase 1 metabolism of drugs and other foreign compounds are reported. The models have been constructed from the CYP102 hemoprotein domain template for which the substrate-bound crystallographic coordinates are available. Selective substrates of individual human P450s: CYP1A2, CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP2E1, CYP3A4, and CYP4A11 are all shown to fit within the corresponding enzymes' active sites in such a manner that is consistent with reported experimental data for both known pathways of substrate metabolism and from the results of site-directed mutagenesis, either in those particular human P450 enzymes concerned or for ones within the same subfamily. The self-consistency of these homology models indicates that they may have potential utility for the pre-screening of novel drug structures.

Factors influencing rates and clearance in P450-mediated reactions: QSARs for substrates of the xenobiotic-metabolizing hepatic microsomal P450s

Various contributory factors associated with the kinetics of cytochrome P450-mediated catalytic activity and the metabolic clearance of drug substrates are discussed and evaluated, based on literature data and physicochemical parameters. Quantitative relationships between molecular structure and biological activity for several series of P450 substrates are presented which point to certain commonalities in P450-catalyzed reactions. In particular, it appears that frontier orbital energies are especially important for the estimation of reaction rates and clearance for many P450 substrates, although occasionally these have to be combined with other descriptors, such as compound lipophilicity (in the form of logP or logD(74)).

Modelling human cytochromes P450 for evaluating drug metabolism: an update

Cytochrome P450 (CYP) enzymes represent the major catalysts for the Phase 1 metabolism of drugs and other xenobiotics in Mammalia, including Homo sapiens. There is considerable current interest in evaluating and, consequently, predicting the metabolic fate of new chemical entities (NCEs) via modelling molecular interactions with P450 constructs, such that sites of metabolism, particular CYP involvement and binding affinities, can be estimated. This paper focuses on the principles for homology modelling of typical enzyme-substrate interactions within the putative active sites of major P450s associated with drug metabolism in man. It also represents an update on previously published work in this journal /1/.

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.

Structural characteristics of human P450s involved in drug metabolism: QSARs and lipophilicity profiles

The factors with human P450 substrate selectivity are reviewed, together with extensive tabulation of quantitative relationships between structure and activity for many species of P450 substrates, inducers and inhibitors. In addition, the physicochemical characteristics (logP and pK(a) values) of heptic microsomal P450 substrates are collated, such that comparisons can be made on the grounds of compound lipophilicities.