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.     

Essential requirements for substrate binding affinity and selectivity toward human CYP2 family enzymes

A detailed analysis of substrate selectivity within the cytochrome P450 2 (CYP2) family is reported. From a consideration of specific interactions between drug substrates for human CYP2 family enzymes and the putative active sites of CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, and CYP2E1, it is likely that the number and disposition of hydrogen bond donor/acceptors and aromatic rings within the various P450 substrate molecules determines their enzyme selectivity and binding affinity, together with directing their preferred routes of metabolism by the CYP2 enzymes concerned. Although many aliphatic residues are present in most P450 active sites, it would appear that their main contribution centers around hydrophobic interactions and desolvation processes accompanying substrate binding. Molecular modeling studies based on the recent CYP2C5 crystal structure appear to show close agreement with site-directed mutagenesis experiments and with information on substrate metabolism and selectivity within the CYP2 family.     

Structure conservation in cytochromes P450

The recent availability of crystal structures for several diverse cytochromes P450 (CYPs) offers the possibility to perform an up-to-date comparative analysis to identify the degree of structure conservation among this superfamily of enzymes specially relevant for their involvement in drug metabolism and toxicity. A set of 9 CYPs sharing between 10% and 27% sequence identity was selected, including 7 class I (CYP 101, 107, 108, 119, 121, 51, and 55) and two class II (CYP 102, and 2C5) structures. After obtaining a multiprotein structure superimposition, a structure-based sequence alignment was derived. Mapping the level of three-dimensional structural conservation onto the sequence alignment revealed that over 28% of the alignment positions have the Calpha carbons of their residues within a root-mean-square deviation (RMSD) of 2 A. This degree of structure conservation is found to be generally preserved, even when the structure undergoes dramatic conformational changes. Performing the analysis on 4 members of the CYP2 family (CYP 2B4, 2C5, 2C8, and 2C9), the percentage of alignment positions within 2 A RMSD amounted to 73%, increasing to over 85% when only structures in a closed conformation are considered. The present findings suggest that it should be plausible to derive models of overall good quality for the major CYP2 metabolizing forms (CYP 2A6, 2C19, 2D6, and 2E1), whereas high levels of uncertainty are still likely to be expected in models for the remaining 2 major P450 metabolizing forms (CYP 1A2 and 3A4), with the corresponding implications for their potential applicability in drug design activities.     

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

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

Engineering human cytochrome P450 enzymes into catalytically self-sufficient chimeras using molecular Lego

The membrane-bound human cytochrome P450s have essential roles in the metabolism of endogenous compounds and drugs. Presented here are the results on the construction and characterization of three fusion proteins containing the N-terminally modified human cytochrome P450s CYP2C9, CY2C19 and CYP3A4 fused to the soluble NADPH-dependent oxidoreductase domain of CYP102A1 from Bacillus megaterium. The constructs, CYP2C9/BMR, CYP2C19/BMR and CYP3A4/BMR are well expressed in Escherichia coli as holo proteins. The chimeras can be purified in the absence of detergent and the purified enzymes are both active and correctly folded in the absence of detergent, as demonstrated by circular dichroism and functional studies. Additionally, in comparison with the parent P450 enzyme, these chimeras have greatly improved solubility properties. The chimeras are catalytically self-sufficient and present turnover rates similar to those reported for the native enzymes in reconstituted systems, unlike previously reported mammalian cytochrome P450 fusion proteins. Furthermore the specific activities of these chimeras are not dependent on the enzyme concentration present in the reaction buffer and they do not require the addition of accessory proteins, detergents or phospholipids to be fully active. The solubility, catalytic self-sufficiency and wild-type like activities of these chimeras would greatly simplify the studies of cytochrome P450 mediated drug metabolism in solution.     

细胞色素P450基因多态性与药物代谢研究进展

细胞色素P450（cytochrome P450,CYP450）在人体药物代谢过程中起着非常重要的作用并参与代谢80%以上的临床药物。由于CYP450在不同种族和不同人群中存在基因多态性,从而造成药物反应的个体差异,一度成为药物基因组学研究的热点。通过查阅国外相关文献,综述了近年来关于CYP1A2、CYP2C9、CYP2C19、CYP2D6和CYP3A4五种主要的药物代谢酶的基因多态性和药物代谢的研究进展,为临床指导个体化用药、避免药物不良反应和新药研发提供科学参考依据。     

Effect of acetylcholinesterase oxime-type reactivators K-48 and HI-6 on human liver microsomal cytochromes P450 in vitro

Substances K-48 and HI-6, oxime-type acetylcholinesterase (AChE) reactivators, were tested for their potential to inhibit the activities of human liver microsomal cytochromes P450 (CYP). The compounds were shown to bind to microsomal cytochromes P450 with spectral binding constants of 0.25 ± 0.05 μM (K-48) and 0.54 ± 0.15 μM (HI-6). To find which cytochrome P450 from the human liver microsomal fraction interacts with these compounds, an inhibition of enzyme activities specific for nine individual CYP enzymes (CYP1A2, CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP2E1, and CYP3A4) was studied. The results have shown no prominent inhibition of individual CYP activities with both compounds except the CYP2E1 activity and the HI-6 reactivator. However, the inhibition of this activity was less than 50% which makes the possible drug interactions highly unlikely. Hence, the interaction of K-48 and HI-6 oxime-type AChE reactivators with human liver microsomal CYP enzymes does not seem to be clinically significant and both compounds could be taken in this respect as antidotal drugs with low risk of drug interactions.     

Prediction of binding modes for ligands in the cytochromes P450 and other heme-containing proteins

The cytochromes P450 (P450s) are a family of heme-containing monooxygenase enzymes involved in a variety of functions, including the metabolism of endogenous and exogenous substances in the human body. During lead optimization, and in drug development, many potential drug candidates are rejected because of the affinity they display for drug-metabolising P450s. Recently, crystal structures of human enzymes involved in drug metabolism have been determined, significantly augmenting the prospect of using structure-based design to modulate the binding and metabolizing properties of compounds against P450 proteins. An important step in the application of structure-based metabolic optimization is the accurate prediction of docking modes in heme binding proteins. In this paper we assess the performance of the docking program GOLD at predicting the binding mode of 45 heme-containing complexes. We achieved success rates of 64% and 57% for Chemscore and Goldscore respectively; these success rates are significantly lower than the value of 79% observed with both scoring functions for the full GOLD validation set. Re-parameterization of metal-acceptor interactions and lipophilicity of planar nitrogen atoms in the scoring functions resulted in a significant increase in the percentage of successful dockings against the heme binding proteins (Chemscore 73%, Goldscore 65%). The modified scoring functions will be useful in docking applications on P450 enzymes and other heme binding proteins.     

CYP3A4 and pregnane X receptor humanized mice

Marked species differences exist in P450 expression and activities. In order to produce mouse models that can be used to more accurately predict human drug and carcinogen metabolism, P450- and xenobiotic receptor humanized mice are being prepared using bacterial artificial chromosomes (BAC) and P1 phage artificial chromosomes (PAC) genomic clones. In some cases, transgenic mice carrying the human genes are bred with null-mice to produce fully humanized mice. Mice expressing human CYP1A1, CYP1A2, CYP2E1, CYP2D6, CYP3A4, and CYP3A7 were generated and characterized. Studies with the CYP3A4-humanized (hCYP3A4) mouse line revealed new information on the physiological function of this P450 and its role in drug metabolism in vivo. With this mouse line, CYP3A4, under certain circumstances, was found to alter the serum levels of estrogen resulting in deficient lactation and low pup survival as a result of underdeveloped mammary glands. This hCYP3A4 mouse established the importance of intestinal CYP3A4 in the pharmacokinetics of orally administered drugs. The hCYP3A4 mice were also used to establish the mechanisms of potential gender differences in CYP3A4 expression (adult female > adult male) that could account for human gender differences in drug metabolism and response. The pregnane X receptor (PXR) is also involved in induction of drug metabolism through its target genes including CYP3A4. Since species differences exist in ligand specificity between human and mice, a PXR-humanized mouse (hPXR) was produced that responds to human PXR activators such as rifampicin but does not respond to the rodent activator pregnenalone 16alpha-carbonitrile.     

Quantitative structure-activity relationships (QSARs) for inhibitors and substrates of CYP2B enzymes: importance of compound lipophilicity in explanation of potency differences

The results of quantitative structure-activity relationship (QSAR) studies on inhibitors and substrates of cytochrome P450 2B (CYP2B) subfamily enzymes are reported. It was found that lipophilicity (in the form of log P) is the most important property for explaining the variations in inhibitory activity, and there are similarities between QSARs for both substrates and inhibitors for CYP2B6 (human), and also between those of other CYP2B enzymes, such as CYP2B1 (rat) and CYP2B4 (rabbit). Both linear and quadratic lipophilicity relationships are evidenced in human and other mammalian species, and the particular type of expression found is probably due to the nature of the compounds under investigation, as it is usually the homologous series which tend to show quadratic relationships in log P. The findings from QSAR studies can be rationalized by molecular modelling of the active site interactions with both P450 crystal structures and homology models of CYP2B subfamily enzymes.     

SOMEViz: A Web Service for Site of Metabolism Estimating and Visualizing

Phase I metabolism is an important consideration in drug discovery because it profoundly affects the toxicity and activity profile of a drug candidate. In these metabolic processes, CYP450 family is responsible for the majority of biotransformation events. However, it is still an important challenge to predict sites of metabolism (SOM) of a new chemical entity due to the complex reaction mechanism and variety in CYP450 enzymes. SOMEViz is an online service designed for predicting and visualizing human cytochromes P450 (CYP450)-mediated sites of metabolism (SOM) of a molecule, on the basis of a previously reported model [1]. The service provides an access for predicting sites of metabolism of molecules with reasonable accuracy, and predicted results are shown in a user-friendly as well as interactive way, which may help chemists explore metabolism properties of chemicals in the early stage of drug discovery. The web-based GUI of SOMEViz offers user a straightforward way to manage and visualize the sites of metabolism (SOM) prediction results. The service and examples are available free of charge at http://www.dddc.ac.cn/some.     

Cytochromes P450 and experimental models of drug metabolism

For the development of new drugs, evaluation of drug-drug interactions with already known compounds, as well as for better understanding of metabolism pathways of various toxicants and pollutants, we studied the drug metabolism mediated by cytochromes P450. The experimental approach is based on animal drug-metabolising systems. From the ethical as well as rational reasons, the selection of an appropriate system is crucial. Here, it is necessary to decide on the basis of expected CYP system involved. For CYP1A-mediated pathways, all the commonly used experimental models are appropriate except probably the dog. On the contrary, the dog seems to be suitable for modelling of processes depending on the CYP2D. With CYP2C, which is possibly the most large and complicated subfamily, the systems based on monkey ( Maccacus rhesus ) may be a good representative. The CYP3A seems to be well modelled by pig or minipig CYP3A29. Detailed studies on activities with individual isolated CYP forms are needed to understand in full all aspects of inter-species differences and variations.     

Wild-type CYP102A1 as a biocatalyst: turnover of drugs usually metabolised by human liver enzymes

This work provides functional data showing that the bacterial CYP102A1 recognises compounds metabolised by human CYP3A4, CYP2E1 and CYP1A2 and is able to catalyse different reactions. Wild-type cytochrome CYP102A1 from Bacillus megaterium is a catalytically self-sufficient enzyme, containing an NADPH-dependent reductase and a P450 haem domain fused in a single polypeptidie chain. An NADPH-dependent method (Tsotsou et al. in Biosens. Bioelectron. 17:119–131, 2002) together with spectroscopic assays were applied to investigate the catalytic activity of CYP102A1 towards 19 xenobiotics, including 17 commercial drugs. These molecules were chosen to represent typical substrates of the five main families of drug-metabolising human cytochromes P450. Liquid chromatography–mass spectrometry analysis showed that CYP102A1 catalyses the hydroxylation of chlorzoxazone, aniline and p-nitrophenol, as well as the N-dealkylation of propranolol and the dehydrogenation of nifedipine. These drugs are typical substrates of human CYP2E1 and CYP3A4. The K _M values calculated for these compounds were in the millimolar range: 1.21 ± 0.07 mM for chlorzoxazone, 2.52 ± 0.08 mM for aniline, 0.81 ± 0.04 mM for propranolol. The values of v _max for chlorzoxazone and propranolol were 46.0 ± 9.0 and 7.6 ± 3.4 nmol min⁻¹ nmol⁻¹, respectively. These values are higher then those measured for the human enzymes. The v _max value for aniline was 9.4 ± 1.3 nmol min⁻¹ nmol⁻¹, comparable to that calculated for human cytochromes P450. The functional data were found to be in line with the sequence alignments, showing that the identity percentage of CYP102A1 with CYP3A4 and CYP2E1 is higher than that found for CYP1A2, CYP2C9 and CYP2D6 families.     

Homology modeling and substrate binding study of human CYP2C9 enzyme

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

Key Mutations Alter the Cytochrome P450 BM3 Conformational Landscape and Remove Inherent Substrate Bias

Cytochrome P450 monooxygenases (P450s) have enormous potential in the production of oxychemicals, due to their unparalleled regio- and stereoselectivity. The Bacillus megaterium P450 BM3 enzyme is a key model system, with several mutants (many distant from the active site) reported to alter substrate selectivity. It has the highest reported monooxygenase activity of the P450 enzymes, and this catalytic efficiency has inspired protein engineering to enable its exploitation for biotechnologically relevant oxidations with structurally diverse substrates. However, a structural rationale is lacking to explain how these mutations have such effects in the absence of direct change to the active site architecture. Here, we provide the first crystal structures of BM3 mutants in complex with a human drug substrate, the proton pump inhibitor omeprazole. Supported by solution data, these structures reveal how mutation alters the conformational landscape and decreases the free energy barrier for transition to the substrate-bound state. Our data point to the importance of such “gatekeeper” mutations in enabling major changes in substrate recognition. We further demonstrate that these mutants catalyze the same 5-hydroxylation reaction as performed by human CYP2C19, the major human omeprazole-metabolizing P450 enzyme.     

Exploring CYP1A1 as anticancer target: homology modeling and in silico inhibitor design

Microsomal cytochrome P450 family 1 enzymes has great importance in the bioactivation of mutagens. P450 catalyzed reactions involve a wide range of substrates, and this versatality is reflected in a structural diversity, evident in the active sites of available P450 structures. This structure offers a template to study further structure-function relationships of alternative substrates and other cytochrome P450 family 1 members. In this paper, we document a homology model of CYP P450 1A1 from Homo sapiens, developed on the basis of template crystal structure of human microsomal P450 1a2 in complex with inhibitor (PDB Id: 2HI4). Homology modeling is performed at the programs, both in the commercial and public realms. We tried to explore CYP1A1 as a potential target for anticancer chemotherapy. To gain an insight into the binding of ligands with enzyme, protein-ligand complex was developed by including information about the known ligand as spatial restraints and optimizing the mutual interactions between the ligand and the binding site. Active site characterization and the study for involvement of specific aminoacids in binding with ligand, facilitates structure based inhibitor design. This study should prove useful in the design and development of potential novel anticancer agents. Figure The refined structure of homology model of CYP1A1     

Some distinctions between flavin-containing and cytochrome P450 monooxygenases

This minireview summarizes information concerning the differences and similarities of the human flavin-containing- (FMO, E.C. 1.14.13.8) and the cytochrome P450-monooxygenases (CYP, E.C. 1.14.14.1). Human FMO oxygenates soft nucleophiles. CYP mainly catalyzes C-H abstraction but also oxidizes nitrogen- and sulfur-containing compounds. Both FMO and CYP generally convert lipophilic compounds into more hydrophilic metabolites. The mechanism by which each monooxygenase operates is quite distinct. Sometimes, CYP or FMO bioactivate chemicals to reactive metabolites but to date, drug toxicity thus far observed in the clinic is mainly the result of CYP-dependent oxidation. Both FMO and CYP possess genetic variability that may contribute to inter-individual variability observed for drug metabolism. In contrast to CYP, FMO is not induced or readily inhibited and potential adverse drug-drug interactions are minimized for drugs prominently metabolized by FMO. These properties may provide advantages in drug design, and by incorporating FMO detoxication pathways into drug candidates, more drug-like materials may emerge.     

Cytochrome P450 expression in human keratinocytes: an aryl hydrocarbon receptor perspective

The goal of this review is to stress the importance of the cytochrome P450 (CYP) superfamily that is expressed in human skin in the hope that it may stimulate further study in an intriguing topic that currently suffers from a relative dearth of information. Like the cells that line the respiratory and GI tracts [X. Ding, L.S. Kaminsky, Human extrahepatic cytochromes P450: function in xenobiotic metabolism and tissue-selective chemical toxicity in the respiratory and gastrointestinal tracts, Annu. Rev. Pharmacol. Toxicol. 43 (2003) 149-173] those present in human skin express a variety of CYPs that play important roles in xenobiotic, drug and steroid metabolism. In addition, a few CYPs, with potentially novel roles in metabolism and keratinocyte function, have recently been discovered that appear to be expressed in a keratinocyte-specific manner [L. Du, S.M. Hoffman, D.S. Keeney, Epidermal CYP2 family cytochromes P450, Toxicol. Appl. Pharmacol. 195 (2004) 278-287]. However, in preparing this review, it soon became apparent that in contrast to the progress made in understanding these events in the liver, relatively little is known in the human skin. Thus, while a number of tantalizing stories are beginning to emerge, they are far from complete. In this review, a brief synopsis of the structure of skin and methods of culturing keratinocytes will be presented. This will be followed by an overview of the various CYPs and their putative regulators that have been currently identified to be expressed in human keratinocytes. Then, a more detailed analysis of CYP regulation that involves the aryl hydrocarbon receptor (AHR) signaling pathway will be offered in the hope that it may serve as a paradigm for other CYP regulatory studies in the skin. Finally, several clinical implications that may arise due to altered regulation of CYPs will be considered.     

Identification of human liver cytochromes P450 by using MALDI-TOF mass spectrometry

Proteomic approaches have been used for detection and identification of cytochromes P450 forms from highly purified membrane preparations of human liver. These included the protein separation by 2D-and/or 1D-electrophoresis and molecular scanning of a SDS-PAGE gel fragment in a range 45–66 kDa (this area corresponds molecular weights of cytochromes P450). The analysis of protein content was statistically evaluated by means of an original 1D-ZOOMER software package which allowed to carry out the processing of mass spectra mixture instead of individual mass spectra used by standard techniques. In the range 45–66 kDa we identified 13 microsomal membrane proteins including such cytochrome P450 forms as CYPs 1A2, 1B1, 2A6, 2E1, 2C8, 2C9, 2C10, 2D6, 3A4, 4A11, 4F2. Study of enzymatic activities of human liver microsomal cytochrome P450 isoforms CYP 1A, 2B, 3A, and 2E revealed the decrease in the rates of O-dealkylation and N-demethylation catalyzed by CYP 450 1A1/1A2 and 3A4 under pathological conditions, whereas 7-benzyloxyresorufin-O-debenzylase activity (which characterizes the total activity of CYP 2B and CYP 2C), the activities of CYP 2E1 (methanol oxidation), 7-pentoxyresorufin-O-dealkylation (CYP 2B), 7-ethoxy-and 7-methoxycoumarin-O-dealkylases (CYP 2B1) remained basically unchanged.