Molecular modelling of human CYP2C subfamily enzymes CYP2C9 and CYP2C19: rationalization of substrate specificity and site-directed mutagenesis experiments in the CYP2C subfamily

Deep mutational scanning of CYP2C19 in human cells reveals a substrate specificity-abundance tradeoff

The cytochrome P450s enzyme family metabolizes ∼80% of small molecule drugs. Variants in cytochrome P450s can substantially alter drug metabolism, leading to improper dosing and severe adverse drug reactions. Due to low sequence conservation, predicting variant effects across cytochrome P450s is challenging. Even closely related cytochrome P450s like CYP2C9 and CYP2C19, which share 92% amino acid sequence identity, display distinct phenotypic properties. Using variant abundance by massively parallel sequencing, we measured the steady-state protein abundance of 7,660 single amino acid variants in CYP2C19 expressed in cultured human cells. Our findings confirmed critical positions and structural features essential for cytochrome P450 function, and revealed how variants at conserved positions influence abundance. We jointly analyzed 4,670 variants whose abundance was measured in both CYP2C19 and CYP2C9, finding that the homologs have different variant abundances in substrate recognition sites within the hydrophobic core. We also measured the abundance of all single and some multiple wild type amino acid exchanges between CYP2C19 and CYP2C9. While most exchanges had no effect, substitutions in substrate recognition site 4 reduced abundance in CYP2C19. Double and triple mutants showed distinct interactions, highlighting a region that points to differing thermodynamic properties between the 2 homologs. These positions are known contributors to substrate specificity, suggesting an evolutionary tradeoff between stability and enzymatic function. Finally, we analyzed 368 previously unannotated human variants, finding that 43% had decreased abundance. By comparing variant effects between these homologs, we uncovered regions underlying their functional differences, advancing our understanding of this versatile family of enzymes.

Functioning of drug-metabolizing microsomal cytochrome P450s: In silico probing of proteins suggests that the distal heme 'active site' pocket plays a relatively 'passive role' in some enzyme-substrate interactions

PURPOSE

The currently held mechanistic understanding of microsomal cytochrome P450s (CYPs) seeks that diverse drug molecules bind within the deep-seated distal heme pocket and subsequently react at the heme centre. To explain a bevy of experimental observations and meta-analyses, we indulge a hypothesis that involves a "diffusible radical mediated" mechanism. This new hypothesis posits that many substrates could also bind at alternate loci on/within the enzyme and be reacted without the pertinent moiety accessing a bonding proximity to the purported catalytic Fe-O enzyme intermediate.

METHODS

Through blind and heme-distal pocket centered dockings of various substrates and non-substrates (drug molecules of diverse sizes, classes, topographies etc.) of microsomal CYPs, we explored the possibility of access of substrates via the distal channels, its binding energies, docking orientations, distance of reactive moieties (or molecule per se) to/from the heme centre, etc. We investigated specific cases like- (a) large drug molecules as substrates, (b) classical marker drug substrates, (c) class of drugs as substrates (Sartans, Statins etc.), (d) substrate preferences between related and unrelated CYPs, (e) man-made site-directed mutants' and naturally occurring mutants' reactivity and metabolic disposition, (f) drug-drug interactions, (g) overall affinities of drug substrate versus oxidized product, (h) meta-analysis of in silico versus experimental binding constants and reaction/residence times etc.

RESULTS

It was found that heme-centered dockings of the substrate/modulator drug molecules with the available CYP crystal structures gave poor docking geometries and distances from Fe-heme centre. In conjunction with several other arguments, the findings discount the relevance of erstwhile hypothesis in many CYP systems. Consequently, the newly proposed hypothesis is deemed a viable alternate, as it satisfies Occam's razor.

CONCLUSIONS

The new proposal affords expanded scope for explaining the mechanism, kinetics and overall phenomenology of CYP mediated drug metabolism. It is now understood that the heme-iron and the hydrophobic distal pocket of CYPs serve primarily to stabilize the reactive intermediate (diffusible radical) and the surface or crypts of the apoprotein bind to the xenobiotic substrate (and in some cases, the heme distal pocket could also serve the latter function). Thus, CYPs enhance reaction rates and selectivity/specificity via a hitherto unrecognized modality.

Conformational adaptation of human cytochrome P450 2B6 and rabbit cytochrome P450 2B4 revealed upon binding multiple amlodipine molecules

Structures of human cytochrome P450 2B6 and rabbit cytochrome P450 2B4 in complex with two molecules of the calcium channel blocker amlodipine have been determined by X-ray crystallography. The presence of two drug molecules suggests clear substrate access channels in each P450. According to a previously established nomenclature, amlodipine molecules were trapped in access pathway 2f in P450 2B6 and in pathway 2a or 2f in P450 2B4. These pathways overlap for part of the length and then diverge as they extend toward the protein surface. A previously described solvent channel was also found in each enzyme. The results indicate that key residues located on the surface and at the entrance of the substrate access channels in each of these P450s may play a crucial role in guiding substrate entry. In addition, the region of P450 2B6 and 2B4 involving helices B', F, F', and G' and part of helix G is substantially more open in the amlodipine complexes than in the corresponding 4-(4-chlorophenyl)imidazole complexes. The increased active site volume observed results from the major retraction of helices F, F', and B' and the β4 sheet region located close to the binding cavity to accommodate amlodipine. These structures demonstrate novel insight into distinct conformational states not observed with previous P450 2B structures and provide clear evidence of the substrate access channels in two drug-metabolizing P450s. In addition, the structures exhibit the versatility that can be exploited via in silico studies with other P450 2B6 ligands as large as raloxifene and itraconazole.

Molecular dynamics simulation studies of betulinic acid with human serum albumin

Betulinic acid (BA) is a naturally occurring pentacyclictriterpenoid possessing anti-retroviral, anti-cancer, and anti-inflammatory properties. Here, we studied the interaction of BA with human serum albumin (HSA) by using molecular docking, and molecular dynamic simulation methods. Molecular docking studies revealed that BA can bind in the large hydrophobic cavity of drug binding site I of sub-domain IIA and IIB, mainly by the hydrophobic interactions and also by hydrogen bond interactions. In which several cyclohexyl groups of BA are interacting with Phe(206), Arg(209), Ala(210), Ala(213), Leu(327), Gly(328), Leu(331), Ala(350), and Lys(351), residues of sub-domain IIA and IIB of HSA by hydrophobic interactions. Also, hydrogen bond interactions were observed between the hydroxyl (OH) group of BA with Phe(206) and Glu(354) of HSA, with hydrogen bond distances of 0.24 nm,0.28 nm respectively. Further, specifically, the molecular dynamics study makes an important contribution in understanding the effect of the binding of BA on conformational changes of HSA and the stability of the protein-drug complex system in aqueous solution. The root mean square deviation values of atoms in the free HSA molecule were calculated from 3000 ps to 5000 ps trajectory and the results were obtained as 0.72 ± 0.036 nm and 0.81 ± 0.032 nm for free HSA and HSA-BA, respectively. The radius of gyration (Rg) values of both unliganded HSA and HSA-BA were stabilized at 2.59 ± 0.03 nm, 2.51 ± 0.01 nm, respectively. Thus, this work may play an important role in the design of new BA inspired drugs with desired HSA binding affinity.

An evaluation of the cytochrome P450 inhibition potential of selected pesticides in human hepatic microsomes

The goal of this work was to study the ability of 18 pesticides to inhibit selective model activities for all major xenobiotic-metabolizing enzymes, namely CYP1A1/2, 2A6, 2B6, 2C8, 2C9, 2C19, 2D6, 2E1 and 3A4. Generally organophosphorus insecticides were the most potent and extensive inhibitors, especially towards CYP1A1/2 (IC(50) values of chlorpyrifos, fenitrothion and profenofos approximately 3 micro M), CYP2B6 (IC(50) values of chlorpyrifos and fenitrothion 2.5 micro M), CYP2C8 (fenitrothion 4.3 micro M), CYP2C9 (fenitrothion and malathion 4.8 and 2.5 micro M, respectively), CYP2D6 (chlorpyrifos and phenthoate approximately 3 micro M) and CYP3A4 (chlorpyrifos, fenitrothion and phenthoate 3-4 micro M). Otherwise there were quite considerable differences in potency and extent of inhibition between different organophosphates. Pyrethroids were in general very weak or inactive. Deltamethrin and fenvalerate were potent inhibitors of CYP2D6 (IC(50) values of approximately 3 micro M) while lambda-cyhalothrin potently inhibited both CYP2D6 and CYP3A4-mediated activities (IC(50)'s about 3-4 micro M). Some pesticides caused relatively potent inhibitions sporadically (carbendazim, CYP2D6, IC(50) = 12 micro M; atrazine, CYP3A4, IC(50) = 2.8 micro M; glyphosate, CYP2C9, IC(50) = 3.7 micro M; hexaflumuron, IC(50) = 6.0 micro M). With the exceptions of alpha-cypermethrin, cypermethrin, isoproturon, carbaryl and abamectin, most pesticides inhibited relatively potently at least one CYP-selective activity, which may have relevance for potential interactions in occupational exposures and for further studies on the CYP-associated metabolism of respective pesticides.

Betulinic acid binding to human serum albumin: a study of protein conformation and binding affinity

Betulinic acid (BA) has anti cancer and anti-HIV activity and has been proved to be therapeutically effective against cancerous and HIV-infected cells. Human serum albumin (HSA) is the predominant protein in the blood. Most drugs that bind to HSA will be transported to other parts of the body. Using micro TOF-Q mass spectrometry, we have shown, for the first time that BA isolated from a plant (Tephrosia calophylla) binds to HSA. The binding constant of BA to HSA was calculated from fluorescence data and found to be K(BA)=1.685+/-0.01 x 10(6) M(-1), indicating a strong binding affinity. The secondary structure of the HSA-BA complex was determined by circular dichroism. The results indicate that the HSA in this complex is partially unfolded. Further, binding of BA at nanomolar concentrations of BA to free HSA was detected using micro TOF-Q mass spectrometry. The study revealed a mass increase from 65199 Da (free HSA) to 65643 Da (HSA+drug), where the additional mass of 444 Da was due to bound BA. Based on the results of this study, it is suggested that micro TOF-Q mass spectrometry is useful technique for drug binding studies.

Pharmacological evaluation of C-3 modified Betulinic acid derivatives with potent anticancer activity

In vitro and in vivo pharmacological screening of Betulinic acid (BA) and five dihydro-BA derivatives modified at C-3 position [4-nitrobenzyl-oximino (1), 2-4-difluoro-benzoyloxy (2), 2-4-difluoro-benzylidene-amino (3), benzoyl-hydrazono (4), and 4-fluorophenyl-hydrazono (5)], having potent in vitro anti-cancer activity was carried out using ADME, animal PK and tumor studies. We found that BA and the derivatives had poor aqueous solubility (<0.1 microg/ml), low to moderate permeability (log Pe<-5.0) and high plasma protein binding (>70%). Although BA and 5 were metabolized by human liver microsomes, derivatives 1, 2, 3 and 4 possessed good in vitro metabolic stability. Except 3 which inhibited CYP1A2 isoform by more than 50% none of the other compounds inhibited key cytochrome P450 enzyme isoforms (CYP1A2, CYP2C9, CYP2D6 and CYP3A4) at 10 microM. Based on in vitro results one derivative 1 was tested in rodent PK and tumor studies. We found that 1 exhibited favorable pharmacokinetic characteristics of a systemically administered drug and showed better in vivo anti-tumor efficacy as compared to BA in a human colon cancer xenograft model. Our results show that BA derivatives are potential anti-cancer compounds which need to be explored in detail.

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.

Cynomolgus Monkey Cytochrome P450 2C43: cDNA Cloning, Heterologous Expression, Purification and Characterization

The cDNA of cytochrome P450 (CYP) 2C43 was cloned from cynomolgus monkey liver by RT-PCR. The deduced amino acid sequence showed 93% and 91% identity to human CYP2C9 and CYP2C19, respectively. The cDNA was expressed in Escherichia coli and purified by a series of chromatography steps, yielding a specific content of 11.5 nmol P450/mg protein. The substrate specificity of the purified CYP2C43 was examined in a reconstitution system comprising NADPH-P450 reductase, lipid, cytochrome b(5) and CYP2C marker substrates. The purified CYP2C43 showed high activity for testosterone 17-oxidation and progesterone 21-hydroxylation, which were also observed for CYP2C19 but not CYP2C9. In addition, CYP2C43 showed activity for (S)-mephenytoin 4'-hydroxylation, a marker reaction for CYP2C19. With CYP2C9 marker substrates, CYP2C43 exhibited low activity for diclofenac 4'-hydroxylation and no activity for tolbutamide p-methylhydroxylation. Therefore, in terms of substrate specificity, our results indicate that CYP2C43 is similar to CYP2C19, rather than CYP2C9.

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.

Computational prediction of human drug metabolism

There is an urgent requirement within the pharmaceutical and biotechnology industries, regulatory authorities and academia to improve the success of molecules that are selected for clinical trials. Although absorption, distribution, metabolism, excretion and toxicity (ADME/Tox) properties are some of the many components that contribute to successful drug discovery and development, they represent factors for which we currently have in vitro and in vivo data that can be modelled computationally. Understanding the possible toxicity and the metabolic fate of xenobiotics in the human body is particularly important in early drug discovery. There is, therefore, a need for computational methodologies for uncovering the relationships between the structure and the biological activity of novel molecules. The convergence of numerous technologies, including high-throughput techniques, databases, ADME/Tox modelling and systems biology modelling, is leading to the foundation of systems-ADME/Tox. Results from experiments can be integrated with predictions to globally simulate and understand the likely complete effects of a molecule in humans. The development and early application of major components of MetaDrug (GeneGo, Inc.) software will be described, which includes rule-based metabolite prediction, quantitative structure-activity relationship models for major drug metabolising enzymes, and an extensive database of human protein-xenobiotic interactions. This represents a combined approach to predicting drug metabolism. MetaDrug can be readily used for visualising Phase I and II metabolic pathways, as well as interpreting high-throughput data derived from microarrays as networks of interacting objects. This will ultimately aid in hypothesis generation and the early triaging of molecules likely to have undesirable predicted properties or measured effects on key proteins and cellular functions.

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.

Assessment of arginine 97 and lysine 72 as determinants of substrate specificity in cytochrome P450 2C9 (CYP2C9)

CYP2C9 is distinguished by a preference for substrates bearing a negative charge at physiological pH. Previous studies have suggested that CYP2C9 residues R97 and K72 may play roles in determining preference for anionic substrates by interaction at the active site or in the access channel. The aim of the present study was to assess the role of these two residues in determining substrate selectivity. R97 and K72 were substituted with negative, uncharged polar and hydrophobic residues using a degenerate polymerase chain reaction-directed strategy. Wild-type and mutant enzymes were expressed in bicistronic format with human cytochrome P450 reductase in Escherichia coli. Mutation of R97 led to a loss of holoenzyme expression for R97A, R97V, R97L, R97T, and R97E mutants. Low levels of hemoprotein were detected for R97Q, R97K, R97I, and R97P mutants. Significant apoenzyme was observed, suggesting that heme insertion or protein stability was compromised in R97 mutants. These observations are consistent with a structural role for R97 in addition to any role in substrate binding. By contrast, all K72 mutants examined (K72E, K72Q, K72V, and K72L) could be expressed as hemoprotein at levels comparable to wild-type. Type I binding spectra were obtained with wild-type and K72 mutants using diclofenac and ibuprofen. Mutation of K72 had little or no effect on the interaction with these substrates, arguing against a critical role in determining substrate specificity. Thus, neither residue appears to play a role in determining substrate specificity, but a structural role for R97 can be proposed consistent with recently published crystallographic data for CYP2C9 and CYP2C5.

Studies of Binding Modes of (S)-Mephenytoin to Wild Types and Mutants of Cytochrome P450 2C19 and 2C9 Using Homology Modeling and Computational Docking

PURPOSE

This study investigated the structural features of CYP2C19 complexed with (S)-mephenytoin, using computational methods. In addition to wild-type CYP2C19 proteins (1A and 1B), which have selective 4-hydroxylase activities of (S)-mephenytoin, CYP2C19 mutants were also studied, together with a wild type and artificial mutants of CYP2C19.

METHODS

Three-dimensional structures of wild-type and mutant proteins of CYP2C19 and CYP2C9 were estimated from homology modeling using the crystal structure of rabbit CYP2C5 as a reference. The binding mode of (S)-mephenytoin to CYP2C19 was investigated using computational docking.

RESULTS

The results reproduced the specific bindings between (S)-mephenytoin and the wild types of CYP2C19. Our findings suggest that Asp293 of CYP2C19 plays an important role in the binding of (S)-mephenytoin, which was surrounded by Vall13 and Ala297, and points the phenyl ring at the heme iron. In addition the wild types of CYP2C19, the computational docking studies also accounted for the experimental activities of CYP2C19 mutants, and wild-type and mutant CYP2C19 proteins.

CONCLUSIONS

These results confirm that the predicted three-dimensional structure of the CYP2C19-(S)-mephenytoin complex is reasonable, and that this strategy is useful for investigating complex structures. Virtual screening for drug discovery can also be carried out using these methods.

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.

In vitro sulfoxidation of thioether compounds by human cytochrome P450 and flavin-containing monooxygenase isoforms with particular reference to the CYP2C subfamily

Cytochrome P450 (P450) and flavin-containing monooxygenase (FMO) enzymes are major catalysts involved in the metabolism of xenobiotics. The sulfoxidation of the thioether pesticides, phorate, disulfoton, sulprofos, and methiocarb, was investigated. Using pooled human liver microsomes (HLMs), thioether compounds displayed similar affinities; however, phorate and disulfoton displayed higher intrinsic clearance rates than either sulprofos or methiocarb. The sulfoxidation of thioethers by HLMs was found to be predominantly P450-driven (85-90%) compared with FMO (10-15%). Among 16 cDNA-expressed human P450 isoforms and 3 human FMO isoforms examined, the following isoforms and their polymorphisms had the highest rates for sulfoxidation, as follows: phorate, CYP1A2, 3A4, 2B6, 2C9*1, 2C18, 2C19, 2D6*1, and FMO1; disulfoton, CYP1A2, 3A4, 2B6, 2C9*1, 2C9*2, 2C18, 2C19, 2D6*1, and FMO1; sulprofos, CYP1A1, 1A2, 3A4, 2C9*1, 2C9*2, 2C9*3, 2C18, 2C19, 2D6*1, and FMO1; methiocarb, CYP1A1, 1A2, 3A4, 2B6, 2C9*1, 2C19, 2D6*1, and FMO1. Among these isoforms, members of the CYP2C subfamily often had the highest affinities and clearance rates. Moreover, sulfaphenazole, a CYP2C9 competitive inhibitor, inhibited disulfoton sulfoxidation by CYP2C9 (IC50 0.84 microM) as well as in HLMs. Ticlopidine, a CYP2C19 mechanism-based inhibitor, inhibited disulfoton sulfoxidation by CYP2C19 (IC50 after coincubation, 43.5 microM; IC50 after preincubation, 4.3 microM) and also in HLMs. Our results indicate that current models of the substrate binding site of the CYP2C subfamily would not effectively predict thioether pesticide metabolism. Thus, the substrate specificity of CYP2Cs is more extensive than is currently believed, and some reevaluation of structure-activity relationships may be required.

Active-site characteristics of CYP2C19 and CYP2C9 probed with hydantoin and barbiturate inhibitors

Three series of N-3 alkyl substituted phenytoin, nirvanol, and barbiturate derivatives were synthesized and their inhibitor potencies were tested against recombinant CYP2C19 and CYP2C9 to probe the interaction of these ligands with the active sites of these enzymes. All compounds were found to be competitive inhibitors of both enzymes, although the degree of inhibitory potency was generally much greater towards CYP2C19. Inhibitor stereochemistry did not markedly influence K(i) towards CYP2C9, and log P adequately predicted inhibitor potency for this enzyme. In contrast, stereochemistry was an important factor in determining inhibitor potency towards CYP2C19. (S)-(+)-N-3-Benzylnirvanol and (R)-(-)-N-3-benzylphenobarbital emerged as the most potent and selective CYP2C19 inhibitors, with K(i) values of < 250nM--at least two orders of magnitude greater inhibitor potency than towards CYP2C9. Both inhibitors were metabolized preferentially at their C-5 phenyl substituents, indicating that CYP2C19 prefers to orient the N-3 substituents away from the active oxygen species. These features were incorporated into expanded CoMFA models for CYP2C9, and a new, validated CoMFA model for CYP2C19.

Chemistry, biological activity, and chemotherapeutic potential of betulinic acid for the prevention and treatment of cancer and HIV infection

3beta-Hydroxy-lup-20(29)-en-28-oic acid (betulinic acid) is a pentacyclic lupane-type triterpene that is widely distributed throughout the plant kingdom. A variety of biological activities have been ascribed to betulinic acid including anti-inflammatory and in vitro antimalarial effects. However, betulinic acid is most highly regarded for its anti-HIV-1 activity and specific cytotoxicity against a variety of tumor cell lines. Interest in developing even more potent anti-HIV agents based on betulinic acid has led to the discovery of a host of highly active derivatives exhibiting greater potencies and better therapeutic indices than some current clinical anti-HIV agents. While its mechanism of action has not been fully determined, it has been shown that some betulinic acid analogs disrupt viral fusion to the cell in a post-binding step through interaction with the viral glycoprotein gp41 whereas others disrupt assembly and budding of the HIV-1 virus. With regard to its anticancer properties, betulinic acid was previously reported to exhibit selective cytotoxicity against several melanoma-derived cell lines. However, more recent work has demonstrated that betulinic acid is cytotoxic against other non-melanoma (neuroectodermal and malignant brain tumor) human tumor varieties. Betulinic acid appears to function by means of inducing apoptosis in cells irrespective of their p53 status. Because of its selective cytotoxicity against tumor cells and favorable therapeutic index, even at doses up to 500 mg/kg body weight, betulinic acid is a very promising new chemotherapeutic agent for the treatment of HIV infection and cancer.

Quantitative structure–activity relationships (QSARs) for substrates of human cytochromes P450 CYP2 family enzymes

The results of quantitative structure-activity relationship (QSAR) studies on substrates of human CYP2 family enzymes are reported, together with those of a small number of CYP2A6, CYP2C19 and CYP2D6 inhibitors. In general, there are good correlations (R = 0.90-0.99) between binding affinity (based on Km or KD values) and various parameters relating to active site interactions such as hydrogen bonding and pi-pi stacking. There is also evidence for the role of compound lipophilicity (as determined by either log P or log D7.4 values) in overall substrate binding affinity, and this could reflect the desolvation energy involved in substrate interaction within the enzyme active site. It is possible to estimate the substrate binding energy for a given P450 from a combination of energy terms relating to hydrogen bonding, pi-pi stacking, desolvation and loss in rotatable bond energy, which agree closely (R = 0.98) with experimental data based on either Km or KD values. Consequently, it is likely that active site interactions represent the major contributory factors to the overall binding affinities for human CYP2 family substrates and, therefore, their estimation is of potential importance for the development of new chemical entities (NCEs) as this can facilitate an assessment of likely metabolic clearance.

STUDIES ON THE METABOLISM OF THE NOVEL, SELECTIVE CYCLOOXYGENASE-2 INHIBITOR INDOMETHACIN PHENETHYLAMIDE IN RAT, MOUSE, AND HUMAN LIVER MICROSOMES: IDENTIFICATION OF ACTIVE METABOLITES

The metabolism of 2-[1-(4-chlorobenzoyl)-5-methoxy-2-methyl-1H-indol-3-yl]-N-phenethyl-acetamide (indomethacin phenethylamide, LM-4108), a highly selective cyclooxygenase-2 inhibitor, was studied in rat, mouse, and human liver microsomes. The primary site of oxidation in all species examined was on the methylene carbons of the phenethyl side chain to form the 1'- and 2'-hydroxy and 2'-oxo metabolites as determined by electrospray ionization liquid chromatography-tandem mass spectrometry. Half-lives for the disappearance of 10 microM LM-4108 in rat, human, and mouse liver microsomes (0.15 pmol P450/ml) were 11 min, 21 min, and 51 min, respectively. Indomethacin formation was not observed in incubations with rat, mouse, or human liver microsomes. Both the 2'-hydroxy-LM-4108 and 2'-oxo-LM-4108 metabolites were synthesized and found to be equipotent to the parent compound with regard to COX-2 inhibitory potency and selectivity [2'-hydroxy-LM-4108: IC(50)(COX-2) = 0.06 microM, IC(50)(COX-1) >66 microM; 2'-oxo-LM-4108: IC(50)(COX-2) = 0.05 microM, IC(50)(COX-1) >66 microM]. The formation of the metabolites was strongly inhibited by specific CYP3A4 inhibitors ketoconazole and troleandomycin but not by other isoform-selective inhibitors. These findings were confirmed by demonstrating that cloned, expressed CYP3A4 catalyzed side chain oxidation. O-Demethylation was a minor oxidative pathway in contrast to the metabolism of indomethacin and was catalyzed by CYP2D6. Upon intravenous administration of LM-4108 to Sprague-Dawley rats, oxidative metabolism on the phenethyl side chain constituted the rate-limiting steps in its clearance. The active metabolites, 2'-oxo- and 2'-hydroxy-LM-4108, as well as 1'-hydroxy-LM-4108, were all observed in rat plasma and thus may contribute to COX-2 inhibition in vivo. The glucuronides of 2'hydroxy-LM-4108 and O-desmethyl-2'-hydroxy-LM-4108 were also identified in rat bile.

Substrate selectivity of human cytochrome P450 2C9: importance of residues 476, 365, and 114 in recognition of diclofenac and sulfaphenazole and in mechanism-based inactivation by tienilic acid

A series of six site-directed mutants of CYP 2C9 were constructed with the aim to better define the amino acid residues that play a critical role in substrate selectivity of CYP 2C9, particularly in three distinctive properties of this enzyme: (i) its selective mechanism-based inactivation by tienilic acid (TA), (ii) its high affinity and hydroxylation regioselectivity toward diclofenac, and (iii) its high affinity for the competitive inhibitor sulfaphenazole (SPA). The S365A mutant exhibited kinetic characteristics for the 5-hydroxylation of TA very similar to those of CYP 2C9; however, this mutant did not undergo any detectable mechanism-based inactivation by TA, which indicates that the OH group of Ser 365 could be the nucleophile forming a covalent bond with an electrophilic metabolite of TA in TA-dependent inactivation of CYP 2C9. The F114I mutant was inactive toward the hydroxylation of diclofenac; moreover, detailed analyses of its interaction with a series of SPA derivatives by difference visible spectroscopy showed that the high affinity of SPA to CYP 2C9 (K(s)=0.4 microM) was completely lost when the phenyl substituent of Phe 114 was replaced with the alkyl group of Ile (K(s)=190+/-20 microM), or when the phenyl substituent of SPA was replaced with a cyclohexyl group (K(s)=120+/-30 microM). However, this cyclohexyl derivative of SPA interacted well with the F114I mutant (K(s)=1.6+/-0.5 microM). At the opposite end, the F94L and F110I mutants showed properties very similar to those of CYP 2C9 toward TA and diclofenac. Finally, the F476I mutant exhibited at least three main differences compared to CYP 2C9: (i) big changes in the k(cat) and K(m) values for TA and diclofenac hydroxylation, (ii) a 37-fold increase of the K(i) value found for the inhibition of CYP 2C9 by SPA, and (iii) a great change in the regioselectivity of diclofenac hydroxylation, the 5-hydroxylation of this substrate by CYP 2C9 F476I exhibiting a k(cat) of 28min(-1). These data indicate that Phe 114 plays an important role in recognition of aromatic substrates of CYP 2C9, presumably via Pi-stacking interactions. They also provide the first experimental evidence showing that Phe 476 plays a crucial role in substrate recognition and hydroxylation by CYP 2C9.

Clinical significance of the cytochrome P450 2C19 genetic polymorphism

Cytochrome P450 2C19 (CYP2C19) is the main (or partial) cause for large differences in the pharmacokinetics of a number of clinically important drugs. On the basis of their ability to metabolise (S)-mephenytoin or other CYP2C19 substrates, individuals can be classified as extensive metabolisers (EMs) or poor metabolisers (PMs). Eight variant alleles (CYP2C19*2 to CYP2C19*8) that predict PMs have been identified. The distribution of EM and PM genotypes and phenotypes shows wide interethnic differences. Nongenetic factors such as enzyme inhibition and induction, old age and liver cirrhosis can also modulate CYP2C19 activity. In EMs, approximately 80% of doses of the proton pump inhibitors (PPIs) omeprazole, lansoprazole and pantoprazole seem to be cleared by CYP2C19, whereas CYP3A is more important in PMs. Five-fold higher exposure to these drugs is observed in PMs than in EMs of CYP2C19, and further increases occur during inhibition of CYP3A-catalysed alternative metabolic pathways in PMs. As a result, PMs of CYP2C19 experience more effective acid suppression and better healing of duodenal and gastric ulcers during treatment with omeprazole and lansoprazole compared with EMs. The pharmacoeconomic value of CYP2C19 genotyping remains unclear. Our calculations suggest that genotyping for CYP2C19 could save approximately 5000 US dollars for every 100 Asians tested, but none for Caucasian patients. Nevertheless, genotyping for the common alleles of CYP2C19 before initiating PPIs for the treatment of reflux disease and H. pylori infection is a cost effective tool to determine appropriate duration of treatment and dosage regimens. Altered CYP2C19 activity does not seem to increase the risk for adverse drug reactions/interactions of PPIs. Phenytoin plasma concentrations and toxicity have been shown to increase in patients taking inhibitors of CYP2C19 or who have variant alleles and, because of its narrow therapeutic range, genotyping of CYP2C19 in addition to CYP2C9 may be needed to optimise the dosage of phenytoin. Increased risk of toxicity of tricyclic antidepressants is likely in patients whose CYP2C19 and/or CYP2D6 activities are diminished. CYP2C19 is a major enzyme in proguanil activation to cycloguanil, but there are no clinical data that suggest that PMs of CYP2C19 are at a greater risk for failure of malaria prophylaxis or treatment. Diazepam clearance is clearly diminished in PMs or when inhibitors of CYP2C19 are coprescribed, but the clinical consequences are generally minimal. Finally, many studies have attempted to identify relationships between CYP2C19 genotype and phenotype and susceptibility to xenobiotic-induced disease, but none of these are compelling.

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.

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.

Structure-activity relationship for human cytochrome P450 substrates and inhibitors

Criteria governing the avidity of substrate binding to human hepatic cytochromes P450 (CYP) associated with Phase 1 metabolism of drugs are described. The results of extensive quantitative structure-activity relationship (QSAR) analyses are reported for substrates of human P450s: CYPIA2, CYP2A6, CYP2B6, CYP2C9, CYP2C19, CYP2D6, CYP2E1, and CYP3A4, representing the enzymes exhibiting major involvement in the metabolism of drug substrates in Homo sapiens. In particular, it is shown that hydrogen bond properties in each class of enzyme-substrate complex are especially important factors in determining substrate binding affinity towards those human P450s which are involved in drug metabolism.

Quantitative structure-activity relationships (QSARs) within substrates of human cytochromes P450 involved in drug metabolism

The results of quantitative structure-activity relationship (QSAR) analyses are reported for structurally diverse series of chemicals which act as substrates or inhibitors for human hepatic microsomal cytochromes P450 (CYP). In particular, this study focuses on the major catalysts of drug metabolism in man, namely CYP1A2, CYP2B6, CYP2C9, CYP2C19, CYP2D6 and CYP3A4. It is found that good correlations (with correlation coefficients ranging from R = 0.94 to 0.99) with P450 binding affinity (Km and K(D)) or competitive inhibition (Ki) values are obtained in each case, especially when consideration of hydrogen bonding parameters are included in the QSAR analysis, together with the number of pi-pi stacking interactions.

Interaction of new sulfaphenazole derivatives with human liver cytochrome p450 2Cs: structural determinants required for selective recognition by CYP 2C9 and for inhibition of human CYP 2Cs

A series of new derivatives of sulfaphenazole (SPA), in which the NH(2) and phenyl substituents of SPA are replaced by various groups or in which the sulfonamide function of SPA is N-alkylated, were synthesized in order to further explore CYP 2C9 active site and to determine the structural factors explaining the selectivity of SPA for CYP 2C9 within the human P450 2C subfamily. Compounds in which the NH(2) group of SPA was replaced with R(1) = CH(3), Br, CH = CH(2), CH(2)CH = CH(2), and CH(2)CH(2)OH exhibited a high affinity for CYP 2C9, as shown by the dissociation constant of their CYP 2C9 complexes, K(s), which was determined by difference visible spectroscopy (K(s) between 0.1 and 0.4 microM) and their constant of CYP 2C9 inhibition (K(i) between 0.3 and 0.6 microM). This indicates that the CYP 2C9-iron(III)-NH(2)R bond previously described to exist in the CYP 2C9-SPA complex does not play a key role in the high affinity of SPA for CYP 2C9. Compounds in which the phenyl group of SPA was replaced with various aryl or alkyl R(2) substituents only exhibited a high affinity for CYP 2C9 if R(2) is a freely rotating and sufficiently electron-rich aryl substituent. Finally, compounds resulting from a N-alkylation of the SPA sulfonamide function (R(3) = CH(3), C(2)H(5), or C(3)H(7)) did not retain the selective inhibitory properties of SPA toward CYP 2C9. However, they are reasonably good inhibitors of CYP 2C8 and CYP 2C18 (IC(50) approximately 20 microM). These data allow one to better understand the structural factors that are important for selective binding in the CYP 2C9 active site. They also provide us with clues towards new selective inhibitors of CYP 2C8 and CYP 2C18.

High-throughput screening approaches for investigating drug metabolism and pharmacokinetics

1. High-throughput screening approaches have been adopted throughout the pharmaceutical industry to aid in the rapid discovery of new chemical entities. Because it is now well recognized that the selection of a robust candidate requires a balance of potency, safety and pharmacokinetics, the role of drug metabolism departments has widened from their traditional one of supporting drug development to include the screening of compounds during the discovery process. To put drug metabolism and pharmacokinetic (DMPK) studies in context, the evolving role of DMPK screening in the drug discovery strategy of pharmaceutical companies will be discussed and a generalized approach will be presented. 2. With the increasing numbers of compounds requiring screening, DMPK optimization methods have had to be adapted for high throughput. There have been many developments in this field over the past decade and this review will focus on the high-throughput DMPK screening methodologies used today and in the recent past. 3. In vitro and in silico (computer-based) methods have proven most amenable to high-throughput approaches and these will firm the bulk of the review, but some advances with in vivo methods will also be discussed. As there has been a vast increase in published material on the topic of high-throughput DMPK methodologies in the past 10 years, it would be impossible to cover every method in detail, so this review will concentrate on the key areas and refer the reader to other, more detailed reviews wherever possible. 4. Most high-throughput methods would not be possible without the enabling technologies of computing, automation, new sample preparation technologies, and highly sensitive and selective detection systems, and these will also be reviewed. 5. The advantages and disadvantages of the screening methods will be presented, in particular the issue of handling the false-positives and -negatives that arise. 6. In concluding the review, future developments in this field will be discussed along with key issues that will need to be addressed.

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/.

In-vitro metabolism of celecoxib, a cyclooxygenase-2 inhibitor, by allelic variant forms of human liver microsomal cytochrome P450 2C9: correlation with CYP2C9 genotype and in-vivo pharmacokinetics

In-vitro studies were conducted to assess the impact of CYP2C9 genotype on the metabolism (methyl hydroxylation) and pharmacokinetics of celecoxib, a novel cyclooxygenase-2 inhibitor and CYP2C9 substrate. When compared to cDNA-expressed wild-type CYP2C9 (CYP2C9*1), the Vmax/Km ratio for celecoxib methyl hydroxylation was reduced by 34% and 90% in the presence of recombinant CYP2C9*2 and CYP2C9*3, respectively. These data indicated that the amino acid substitution at position 359 (Ile to Leu) elicited a more pronounced effect on the metabolism of celecoxib than did a substitution at position 144 (Arg to Cys). The Vmax/Km ratio was also decreased in microsomes of livers genotyped CYP2C9*1/*2 (47% decrease, mean of two livers), or CYP2C9*1/*3 (59% decrease, one liver). In all cases, these changes were largely reflective of a decrease in Vmax, with a minimal change in Km. Based on simulations of the in-vitro data obtained with the recombinant CYP2C9 proteins, it was anticipated that the pharmacokinetics of celecoxib (as a much as a five-fold increase in plasma AUC) would be altered (versus CYP2C9*1/*1 subjects) in subjects genotyped heterozygous or homozygous for the CYP2C9*2 (Cys144) or CYP2C9*3 (Leu359) allele. In a subsequent clinical study, the AUC of celecoxib was increased (versus CYP2C9*1/*1 subjects) approximately 2.2-fold (range, 1.6-3-fold) in two CYP2C9*1/*3 subjects and one CYP2C9*3/*3 subject receiving a single oral dose (200 mg) of the drug. In contrast, there was no significant change in celecoxib AUC in two subjects genotyped CYP2C9*1/*2.

The CYP2C19 enzyme polymorphism

The genetic test is gradually replacing probe drugs as the primary tool for screening populations for the CYP2C19 polymorphism. A full appreciation for the clinical and toxicological relevance of this genetic variation is presently limited. Further research is needed in several areas. The development and use of 3-D models of the CYP2C19 enzyme to automate and increase the rate at which CYP2C19 substrates are identified could reap great benefits. Meanwhile, clinical research should begin to determine whether the CYP2C19 polymorphism affects therapeutic outcomes and toxicity of drugs in actual patient settings. Combining research efforts in molecular modeling, genetic testing, clinical and epidemiological research will be required if better appreciation of this genetic variation and its importance in the population at large is to emerge.

High-throughput screening in drug metabolism and pharmacokinetic support of drug discovery

The application of rapid methods currently used for screening discovery drug candidates for metabolism and pharmacokinetic characteristics is discussed. General considerations are given for screening in this context, including the criteria for good screens, the use of counterscreens, the proper sequencing of screens, ambiguity in the interpretation of results, strategies for false positives and negatives, and the special difficulties encountered in drug metabolism and pharmacokinetic screening. Detailed descriptions of the present status of screening are provided for absorption potential, blood-brain barrier penetration, inhibition and induction of cytochrome P450, pharmacokinetics, biotransformation, and computer modeling. Although none of the systems currently employed for drug metabolism and pharmacokinetic screening can be considered truly high-throughput, several of them are rapid enough to be a practical part of the screening paradigm for modern, fast-moving discovery programs.

On the recognition of mammalian microsomal cytochrome P450 substrates and their characteristics: towards the prediction of human p450 substrate specificity and metabolism

The characteristics of mammalian microsomal P450 xenobiotic substrates are described, particularly with reference to the major P450 isoforms associated with drug metabolism in humans. It is further reported that a relatively small number of molecular, electronic, and physico-chemical properties are required to discriminate between chemicals that exhibit specificity for human P450 isoforms: CYP1A2, CYP2A6, CYP2B6, CYP2C9, CYP2C19, CYP2D6, CYP2E1, and CYP3A4. Molecular templates of superimposed substrates are shown to be complementary with the putative active sites of the relevant enzymes, thus enabling a possible prediction of P450 specificity from structure. Factors contributing to metabolic clearance and binding affinity are also discussed, and methods for their calculation are described.

Arginines 97 and 108 in CYP2C9 are important determinants of the catalytic function

Human cytochrome P450 2C9 (CYP2C9) is one of the major drug metabolising enzymes which exhibits a broad substrate specificity. The B-C loop is located in the active-site but has been difficult to model, owing to its diverse and flexible structure. To elucidate the function of the B-C loop we used homology modelling based on the Cyp102 structure in combination with functional studies of mutants using diclofenac as a model substrate for CYP2C9. The study shows the importance of the conserved arginine in position 97 and the arginine in position 108 for the catalytic function. The R97A mutant had a 13-fold higher K(m) value while the V(max) was in the same order as the wild type. The R108 mutant had a 100-fold lower activity with diclofenac compared to the wild-type enzyme. The other six mutants (S95A, F100A, L102A, E104A, R105A, and N107A) had kinetic parameters similar to the CYP2C9 wild-type. Our homology model based on the CYP102 structure as template indicates that R97, L102, and R105 are directed into the active site, whereas R108 is not. The change in catalytic function when arginine 97 was replaced with alanine and the orientation of this amino acid in our homology model indicates its importance for substrate interaction.

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.

Cytochrome P450 substrate specificities, substrate structural templates and enzyme active site geometries

The structural characteristics of human cytochrome P450 substrates are outlined in the light of extensive studies on P450 substrate specificity. Templates of superimposed substrates for individual P450 isozymes are shown to fit the corresponding enzyme active sites, where contacts with specific amino acid residues appear to be involved in the interaction with each structural template. Procedures leading to the evaluation of likely P450 specificity, binding affinity and rate of metabolism are described in the context of key examples in which molecular modelling appears to rationalize experimentally observed findings.

Theoretical investigation of substrate specificity for cytochromes P450 IA2, P450 IID6 and P450 IIIA4

Three-dimensional models of the cytochromes P450 IA2, P450 IID6 and P450 IIIA4 were built by means of comparative modeling using the X-ray crystallographic structures of P450 CAM, P450 BM-3, P450 TERP and P450 ERYF as templates. The three cytochromes were analyzed both in their intrinsic structural features and in their interaction properties with fifty specific and non-specific substrates. Substrate/enzyme complexes were obtained by means of both automated rigid and flexible body docking. The comparative analysis of the three cytochromes and the selected substrates, in their free and bound forms, allowed for the building of semi-quantitative models of substrate specificity based on both molecular and intermolecular interaction descriptors. The results of this study provide new insights into the molecular determinants of substrate specificity for the three different eukaryotic P450 isozymes and constitute a useful tool for predicting the specificity of new compounds.

Molecular modelling of CYP1 family enzymes CYP1A1, CYP1A2, CYP1A6 and CYP1B1 based on sequence homology with CYP102

Molecular modelling of a number of CYP1 family enzymes from rat, plaice and human is described based on amino acid sequence homology with the haemoprotein domain of CYP102, a unique bacterial P450 of known structure. The interaction of various substrates and inhibitors within the putative active sites of rat CYP1A1, human CYP1A2, a fish CYP1 enzyme CYP1A6 (from plaice) and human CYP1B1, is shown to be consistent with P450-mediated oxidation in each example or, in the case of inhibitors, mechanism of inhibition. It is reported that relatively small changes between the enzymes' active site regions assist in the rationalization of CYP1 enzyme preferences for particular substrate types, and a template of superimposed CYP1A2 substrates is shown to fit the putative active site of the human CYP1A2 enzyme.

In vitro inhibition of human liver drug metabolizing enzymes by second generation antihistamines

Cetirizine, terfenadine, loratadine, astemizole and mizolastine were compared for their ability to inhibit marker activities for CYP1A2, CYP2C9, CYP2C19, CYP2D6, CYP3A4 and for some glucuronidation isoenzymes in human liver microsomes. The most pronounced effects were observed with terfenadine, astemizole and loratadine which inhibited CYP3A4-mediated testosterone 6beta-hydroxylation (IC50 of 23, 21 and 32 microM, respectively) and CYP2D6-mediated dextromethorphan O-demethylation (IC50 of 18, 36 and 15 microM, respectively). In addition, loratadine markedly inhibited the CYP2C19 marker activity, (S)-mephenytoin 4-hydroxylation (Ki of 0.17 microM). Furthermore, loratadine activated the CYP2C9-catalyzed tolbutamide hydroxylation (ca. 3-fold increase at 30 microM) and inhibited some glucuronidation enzymes. Mizolastine appeared to be a relatively weak and unspecific inhibitor of CYP2E1, CYP2C9, CYP2D6 and CYP3A4 (IC50Ss in the 100 micromolar range). Cetirizine demonstrated no effect on the investigated activities. A comparison of the inhibitory potencies of cetirizine, terfenadine, loratidine, astemizole and mizolastine with their corresponding plasma concentrations in humans suggests that these antihistamines are not likely to interfere with the metabolic clearance of coadministered drugs, with the exception of loratidine, which appears to inhibit CYP2C19 with sufficient potency to warrant additional investigation.

Use and value of metabolism databases

A key objective during drug discovery is to ensure selection of lead compounds for development that have the optimum pharmacodynamic profile without undesirable toxicity. Metabolism plays a key role in determining the biological activity of compounds in vivo, and is therefore an important factor to consider during the discovery phase. The ability to predict the biotransformation pathways for specific chemical structures provides the opportunity to implement structural modifications that might advantageously modify these processes. Potentially valuable tools include knowledge databases that contain searchable information on the known metabolites of existing compounds and databases that are designed to predict the metabolites.

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

Molecular modelling of human cytochrome P450 (CYP) isoforms is described, based on amino acid sequence homology with a unique bacterial P450 (CYP102) of known crystal structure. It is found that for the human hepatic P450s involved in the metabolism of xenobiotics, ie. CYPIA2, CYP 1A6, CYP2B6, CYP2C9, CYP2C 19, CYP2D6, CYP2E1 and CYP3A4, there is a satisfactory agreement between specific substrate characteristics and topographical features of the putative active sites, including complementarity with key amino acid residues in the P450 haem environments. A combination of homology model interactions with substrates and certain molecular properties of the compounds themselves provides a methodology for the evaluation of potential P450 selectivity in new chemical entities (NCEs).

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.

Molecular modelling and quantitative structure-activity relationship studies on the interaction of omeprazole with cytochrome P450 isozymes

Molecular modelling of the anti-ulcerative agent, omeprazole, with the putative active sites of cytochromes P4503A4 and P4502C19, enzymes which are the major catalysts of omeprazole metabolism in man, are reported. Interactive docking of omeprazole in both CYP3A4 and CYP2C19 gives rise to binding orientations which are consistent with both the known sites of metabolism reported for these isoforms and with evidence from site-directed mutagenesis experiments on CYP2C19, a P450 associated with genetic polymorphism in human drug metabolism. The potential P450 enzymic interactions, inhibition and induction of omeprazole are discussed in the light of molecular modelling and QSAR (quantitative structure-activity relationship) studies on related compounds.