Crystal Structure of Human Cytochrome P450 2D6

Marmoset cytochrome P450 2D8 in livers and small intestines metabolizes typical human P450 2D6 substrates, metoprolol, bufuralol and dextromethorphan

1. Although the New World non-human primate, the common marmoset (Callithrix jacchus), is a potentially useful animal model, comprehensive understanding of drug metabolizing enzymes is insufficient. 2. A cDNA encoding a novel cytochrome P450 (P450) 2D8 was identified in marmosets. The amino acid sequence deduced from P450 2D8 cDNA showed a high sequence identity (83-86%) with other primate P450 2Ds. Phylogenetic analysis showed that marmoset P450 2D8 was closely clustered with human P450 2D6, unlike P450 2Ds of miniature pig, dog, rabbit, guinea pig, mouse or rat. 3. Marmoset P450 2D8 mRNA was predominantly expressed in the liver and small intestine among the tissues types analyzed, whereas marmoset P450 2D6 mRNA was expressed predominantly in the liver where P450 2D protein was detected by immunoblotting. 4. By metabolic assays using marmoset P450 2D8 protein heterologously expressed in Escherichia coli, although P450 2D8 exhibits lower catalytic efficiency compared to marmoset and human P450 2D6 enzymes, P450 2D8 mediated O-demethylations of metoprolol and dextromethorphan and bufuralol 1'-hydroxylation. 5. These results suggest that marmoset P450 2D8 (also expressed in the extrahepatic tissues) has potential roles in drug metabolism in a similar manner to those of human and marmoset P450 2D6.

Comparison of intrinsic dynamics of cytochrome p450 proteins using normal mode analysis

Cytochrome P450 enzymes are hemeproteins that catalyze the monooxygenation of a wide-range of structurally diverse substrates of endogenous and exogenous origin. These heme monooxygenases receive electrons from NADH/NADPH via electron transfer proteins. The cytochrome P450 enzymes, which constitute a diverse superfamily of more than 8,700 proteins, share a common tertiary fold but < 25% sequence identity. Based on their electron transfer protein partner, cytochrome P450 proteins are classified into six broad classes. Traditional methods of pro are based on the canonical paradigm that attributes proteins' function to their three-dimensional structure, which is determined by their primary structure that is the amino acid sequence. It is increasingly recognized that protein dynamics play an important role in molecular recognition and catalytic activity. As the mobility of a protein is an intrinsic property that is encrypted in its primary structure, we examined if different classes of cytochrome P450 enzymes display any unique patterns of intrinsic mobility. Normal mode analysis was performed to characterize the intrinsic dynamics of five classes of cytochrome P450 proteins. The present study revealed that cytochrome P450 enzymes share a strong dynamic similarity (root mean squared inner product > 55% and Bhattacharyya coefficient > 80%), despite the low sequence identity (< 25%) and sequence similarity (< 50%) across the cytochrome P450 superfamily. Noticeable differences in Cα atom fluctuations of structural elements responsible for substrate binding were noticed. These differences in residue fluctuations might be crucial for substrate selectivity in these enzymes.

Prediction of cytochrome P450 mediated metabolism

Cytochrome P450 enzymes (CYPs) form one of the most important enzyme families involved in the metabolism of xenobiotics. CYPs comprise many isoforms, which catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be formed. However, it is often hard to rationalize what metabolites these enzymes generate. In recent years, many different in silico approaches have been developed to predict binding or regioselective product formation for the different CYP isoforms. These comprise ligand-based methods that are trained on experimental CYP data and structure-based methods that consider how the substrate is oriented in the active site or/and how reactive the part of the substrate that is accessible to the heme group is. We will review key aspects for various approaches that are available to predict binding and site of metabolism (SOM), what outcome can be expected from the predictions, and how they could potentially be improved.

Molecular evolution of the CYP2D subfamily in primates: purifying selection on substrate recognition sites without the frequent or long-tract gene conversion

The human cytochrome P450 (CYP) 2D6 gene is a member of the CYP2D gene subfamily, along with the CYP2D7P and CYP2D8P pseudogenes. Although the CYP2D6 enzyme has been studied extensively because of its clinical importance, the evolution of the CYP2D subfamily has not yet been fully understood. Therefore, the goal of this study was to reveal the evolutionary process of the human drug metabolic system. Here, we investigate molecular evolution of the CYP2D subfamily in primates by comparing 14 CYP2D sequences from humans to New World monkey genomes. Window analysis and statistical tests revealed that entire genomic sequences of paralogous genes were extensively homogenized by gene conversion during molecular evolution of CYP2D genes in primates. A neighbor-joining tree based on genomic sequences at the nonsubstrate recognition sites showed that CYP2D6 and CYP2D8 genes were clustered together due to gene conversion. In contrast, a phylogenetic tree using amino acid sequences at substrate recognition sites did not cluster the CYP2D6 and CYP2D8 genes, suggesting that the functional constraint on substrate specificity is one of the causes for purifying selection at the substrate recognition sites. Our results suggest that the CYP2D gene subfamily in primates has evolved to maintain the regioselectivity for a substrate hydroxylation activity between individual enzymes, even though extensive gene conversion has occurred across CYP2D coding sequences.

Overview of nanoscale NEXAFS performed with soft X-ray microscopes

Today, in material science nanoscale structures are becoming more and more important. Not only for the further miniaturization of semiconductor devices like carbon nanotube based transistors, but also for newly developed efficient energy storage devices, gas sensors or catalytic systems nanoscale and functionalized materials have to be analysed. Therefore, analytical tools like near-edge X-ray absorption fine structure (NEXAFS) spectroscopy has to be applied on single nanostructures. Scanning transmission X-ray microscopes (STXM) as well as full-field transmission X-ray microscopes (TXM) allow the required spatial resolution to study individual nanostructures. In the soft X-ray energy range only STXM was used so far for NEXAFS studies. Due to its unique setup, the TXM operated by the Helmholtz-Zentrum Berlin (HZB) at the electron storage ring BESSY II is the first one in the soft X-ray range which can be used for NEXAFS spectroscopy studies which will be shown in this review. Here we will give an overview of the different microscopes used for NEXAFS studies and describe their advantages and disadvantages for different samples.

Contributions of ionic interactions and protein dynamics to cytochrome P450 2D6 (CYP2D6) substrate and inhibitor binding

P450 2D6 contributes significantly to the metabolism of >15% of the 200 most marketed drugs. Open and closed crystal structures of P450 2D6 thioridazine complexes were obtained using different crystallization conditions. The protonated piperidine moiety of thioridazine forms a charge-stabilized hydrogen bond with Asp-301 in the active sites of both complexes. The more open conformation exhibits a second molecule of thioridazine bound in an expanded substrate access channel antechamber with its piperidine moiety forming a charge-stabilized hydrogen bond with Glu-222. Incubation of the crystalline open thioridazine complex with alternative ligands, prinomastat, quinidine, quinine, or ajmalicine, displaced both thioridazines. Quinine and ajmalicine formed charge-stabilized hydrogen bonds with Glu-216, whereas the protonated nitrogen of quinidine is equidistant from Asp-301 and Glu-216 with protonated nitrogen H-bonded to a water molecule in the access channel. Prinomastat is not ionized. Adaptations of active site side-chain rotamers and polypeptide conformations were evident between the complexes, with the binding of ajmalicine eliciting a closure of the open structure reflecting in part the inward movement of Glu-216 to form a hydrogen bond with ajmalicine as well as sparse lattice restraints that would hinder adaptations. These results indicate that P450 2D6 exhibits sufficient elasticity within the crystal lattice to allow the passage of compounds between the active site and bulk solvent and to adopt a more closed form that adapts for binding alternative ligands with different degrees of closure. These crystals provide a means to characterize substrate and inhibitor binding to the enzyme after replacement of thioridazine with alternative compounds.

Conversions of tricyclic antidepressants and antipsychotics with selected P450s from Sorangium cellulosum So ce56

Human cytochromes P450 (P450s) play a major role in the biotransformation of drugs. The generated metabolites are important for pharmaceutical, medical, and biotechnological applications and can be used for derivatization or toxicological studies. The availability of human drug metabolites is restricted and alternative ways of production are requested. For this, microbial P450s turned out to be a useful tool for the conversion of drugs and related derivatives. Here, we used 10 P450s from the myxobacterium Sorangium cellulosum So ce56, which have been cloned, expressed, and purified. The P450s were investigated concerning the conversion of the antidepressant drugs amitriptyline, clomipramine, imipramine, and promethazine; the antipsychotic drugs carbamazepine, chlorpromazine, and thioridazine, as well as their precursors, iminodibenzyl and phenothiazine. Amitriptyline, chlorpromazine, clomipramine, imipramine, and thioridazine are efficiently converted during the in vitro reaction and were chosen to upscale the production by an Escherichia coli-based whole-cell bioconversion system. Two different approaches, a whole-cell system using M9CA medium and a system using resting cells in buffer, were used for the production of sufficient amounts of metabolites for NMR analysis. Amitriptyline, clomipramine, and imipramine are converted to the corresponding 10-hydroxylated products, whereas the conversion of chlorpromazine and thioridazine leads to a sulfoxidation in position 5. It is shown for the first time that myxobacterial P450s are efficient to produce known human drug metabolites in a milligram scale, revealing their ability to synthesize pharmaceutically important compounds.

Molecular dynamics of CYP2D6 polymorphisms in the absence and presence of a mechanism-based inactivator reveals changes in local flexibility and dominant substrate access channels

Cytochrome P450 enzymes (CYPs) represent an important enzyme superfamily involved in metabolism of many endogenous and exogenous small molecules. CYP2D6 is responsible for ∼ 15% of CYP-mediated drug metabolism and exhibits large phenotypic diversity within CYPs with over 100 different allelic variants. Many of these variants lead to functional changes in enzyme activity and substrate selectivity. Herein, a molecular dynamics comparative analysis of four different variants of CYP2D6 was performed. The comparative analysis included simulations with and without SCH 66712, a ligand that is also a mechanism-based inactivator, in order to investigate the possible structural basis of CYP2D6 inactivation. Analysis of protein stability highlighted significantly altered flexibility in both proximal and distal residues from the variant residues. In the absence of SCH 66712, *34, *17-2, and *17-3 displayed more flexibility than *1, and *53 displayed more rigidity. SCH 66712 binding reversed flexibility in *17-2 and *17-3, through *53 remained largely rigid. Throughout simulations with docked SCH 66712, ligand orientation within the heme-binding pocket was consistent with previously identified sites of metabolism and measured binding energies. Subsequent tunnel analysis of substrate access, egress, and solvent channels displayed varied bottle-neck radii. Taken together, our results indicate that SCH 66712 should inactivate these allelic variants, although varied flexibility and substrate binding-pocket accessibility may alter its interaction abilities.

Mechanism-based inactivation of human cytochrome P450 3A4 by two piperazine-containing compounds

Human cytochrome P450 3A4 (CYP3A4) is responsible for the metabolism of more than half of pharmaceutic drugs, and inactivation of CYP3A4 can lead to adverse drug-drug interactions. The substituted imidazole compounds 5-fluoro-2-[4-[(2-phenyl-1H-imidazol-5-yl)methyl]-1-piperazinyl]pyrimidine (SCH 66712) and 1-[(2-ethyl-4-methyl-1H-imidazol-5-yl)methyl]-4-[4-(trifluoromethyl)-2-pyridinyl]piperazine (EMTPP) have been previously identified as mechanism-based inactivators (MBI) of CYP2D6. The present study shows that both SCH 66712 and EMTPP are also MBIs of CYP3A4. Inhibition of CYP3A4 by SCH 66712 and EMTPP was determined to be concentration, time, and NADPH dependent. In addition, inactivation of CYP3A4 by SCH 66712 was shown to be unaffected by the presence of electrophile scavengers. SCH 66712 displays type I binding to CYP3A4 with a spectral binding constant (Ks) of 42.9 ± 2.9 µM. The partition ratios for SCH 66712 and EMTPP were 11 and 94, respectively. Whole protein mass spectrum analysis revealed 1:1 binding stoichiometry of SCH 66712 and EMTPP to CYP3A4 and a mass increase consistent with adduction by the inactivators without addition of oxygen. Heme adduction was not apparent. Multiple mono-oxygenation products with each inactivator were observed; no other products were apparent. These are the first MBIs to be shown to be potent inactivators of both CYP2D6 and CYP3A4.

Computational approaches elucidate the allosteric mechanism of human aromatase inhibition: a novel possible route to Small-molecule regulation of CYP450s activities?

Human aromatase (HA) is a P450 cytochrome (CYP) with an essential role in estrogen biosynthesis. Since more than 70% of breast cancers are positive for estrogenic receptor (ER), the reduction of estrogen physiological concentrations through HA inhibition is one of most important therapeutic strategies against this cancer type. Recently, experimental evidence showed that selected taxmoxifen metabolites, which are typically used as estrogen receptor modulators (SERMs), inhibit HA through an allosteric mechanism. In this work, we present a computational protocol to (i) characterize the structural framework and (ii) define the atomistic details of the determinants for the noncompetitive inhibition mechanism. Our calculations identify two putative binding sites able to efficiently bind all tamoxifen metabolites. Analysis of long-scale molecular dynamics simulations reveal that endoxifen, the most effective noncompetitive inhibitor, induces significant enzyme rigidity by binding in one of the possible peripheral sites. The consequence of this binding event is the suppression of one of the functional enzymatic collective motions associated with breathing of the substrate access channel. Moreover, an internal dynamics-based alignment of HA with six other human cytochromes shows that this collective motion is common to other members of the CYP450 protein family. On this basis, our findings may thus be of help for the development of new (pan)inhibitors for the therapeutic treatment of cancer, targeting and modulating the activity of HA and of estrogen receptor, and may also stimulate the development of new drug design strategies for chemoprevention and chemoprotection via allosteric inhibition of CYP450 proteins.

Including ligand-induced protein flexibility into protein tunnel prediction

In proteins with buried active sites, understanding how ligands migrate through the tunnels that connect the exterior of the protein to the active site can shed light on substrate specificity and enzyme function. A growing body of evidence highlights the importance of protein flexibility in the binding site on ligand binding; however, the influence of protein flexibility throughout the body of the protein during ligand entry and egress is much less characterized. We have developed a novel tunnel prediction and evaluation method named IterTunnel, which includes the influence of ligand-induced protein flexibility, guarantees ligand egress, and provides detailed free energy information as the ligand proceeds along the egress route. IterTunnel combines geometric tunnel prediction with steered molecular dynamics in an iterative process to identify tunnels that open as a result of ligand migration and calculates the potential of mean force of ligand egress through a given tunnel. Applying this new method to cytochrome P450 2B6, we demonstrate the influence of protein flexibility on the shape and accessibility of tunnels. More importantly, we demonstrate that the ligand itself, while traversing through a tunnel, can reshape tunnels due to its interaction with the protein. This process results in the exposure of new tunnels and the closure of preexisting tunnels as the ligand migrates from the active site.

A computational approach predicting CYP450 metabolism and estrogenic activity of an endocrine disrupting compound (PCB-30)

Endocrine disrupting chemicals influence growth and development through interactions with the hormone system, often through binding to hormone receptors such as the estrogen receptor. Computational methods can predict endocrine disrupting chemical activity of unmodified compounds, but approaches predicting activity following metabolism are lacking. The present study uses a well-known environmental contaminant, PCB-30 (2,4,6-trichlorobiphenyl), as a prototype endocrine disrupting chemical and integrates predictive (computational) and experimental methods to determine its metabolic transformation by cytochrome P450 3A4 (CYP3A4) and cytochrome P450 2D6 (CYP2D6) into estrogenic byproducts. Computational predictions suggest that hydroxylation of PCB-30 occurs at the 3- or 4-phenol positions and leads to metabolites that bind more strongly than the parent molecule to the human estrogen receptor alpha (hER-α). Gas chromatography-mass spectrometry experiments confirmed that the primary metabolite for CYP3A4 and CYP2D6 is 4-hydroxy-PCB-30, and the secondary metabolite is 3-hydroxy-PCB-30. Cell-based bioassays (bioluminescent yeast expressing hER-α) confirmed that hydroxylated metabolites are more estrogenic than PCB-30. These experimental results support the applied model's ability to predict the metabolic and estrogenic fate of PCB-30, which could be used to identify other endocrine disrupting chemicals involved in similar pathways.

Exposition and reactivity optimization to predict sites of metabolism in chemicals

Chemical modifications of drugs induced by phase I biotransformations significantly affect their pharmacokinetic properties. Because the metabolites produced can themselves have a pharmacological effect and an intrinsic toxicity, medicinal chemists need to accurately predict the sites of metabolism (SoM) of drugs as early as possible. However, site of metabolism prediction is rarely accompanied by a prediction of the relative abundance of the various metabolites. Such a prediction would be a great help in the study of drug– drug interactions and in the process of reducing the toxicity of potential drug candidates. The aim of this paper is to present recent developments in the prediction of xenobiotic metabolism and to use concrete examples to explain the computational mechanism employed.

Entropic and enthalpic contributions to stereospecific ligand binding from enhanced sampling methods

The stereoselective binding of R- and S-propranolol to the metabolic enzyme cytochrome P450 2D6 and its mutant F483A was studied using various computational approaches. Previously reported free-energy differences from Hamiltonian replica exchange simulations, combined with thermodynamic integration, are compared to the one-step perturbation approach, combined with local-elevation enhanced sampling, and an excellent agreement between methods was obtained. Further, the free-energy differences are interpreted in terms of enthalpic and entropic contributions where it is shown that exactly compensating contributions obscure a molecular interpretation of differences in the affinity while various reduced terms allow a more detailed analysis, which agree with heuristic observations on the interactions.

The endogenous substrates of brain CYP2D

CYP2D6, one of the major cytochrome P450 isoforms present in the human brain, is associated with the incidence and prevalence of central nervous system (CNS) diseases. Human CYP2D6 and rat CYP2D are involved in the metabolism of various neurotransmitters and neurosteroids. Brain CYP2D can be regulated by endogenous steroids, including sex hormones. The alteration of CYP2D-mediated metabolism induced by endogenous steroids may cause changes in sensitivity to environmental and industrial toxins and carcinogens as well as physiological and pathophysiological processes controlled by biologically active compounds. This review summarizes the current knowledge regarding the distribution, endogenous substrates, and regulation of brain CYP2D.

CYP 2D6 binding affinity predictions using multiple ligand and protein conformations

Because of the large flexibility and malleability of Cytochrome P450 enzymes (CYPs), in silico prediction of CYP binding affinities to drugs and other xenobiotic compounds is a true challenge. In the current work, we use an iterative linear interaction energy (LIE) approach to compute CYP binding affinities from molecular dynamics (MD) simulation. In order to improve sampling of conformational space, we combine results from simulations starting with different relevant protein-ligand geometries. For calculated binding free energies of a set of thiourea compounds binding to the flexible CYP 2D6 isoform, improved correlation with experiment was obtained by combining results of MD runs starting from distinct protein conformations and ligand-binding orientations. This accuracy was obtained from relatively short MD simulations, which makes our approach computationally attractive for automated calculations of ligand-binding affinities to flexible proteins such as CYPs.

Dynamics and flexibility of human aromatase probed by FTIR and time resolved fluorescence spectroscopy

Human aromatase (CYP19A1) is a steroidogenic cytochrome P450 converting androgens into estrogens. No ligand-free crystal structure of the enzyme is available to date. The crystal structure in complex with the substrate androstenedione and the steroidal inhibitor exemestane shows a very compact conformation of the enzyme, leaving unanswered questions on the conformational changes that must occur to allow access of the ligand to the active site. As H/D exchange kinetics followed by FTIR spectroscopy can provide information on the conformational changes in proteins where solvent accessibility is affected, here the amide I region was used to measure the exchange rates of the different elements of the secondary structure for aromatase in the ligand-free form and in the presence of the substrate androstenedione and the inhibitor anastrozole. Biphasic exponential functions were found to fit the H/D exchange data collected as a function of time. Two exchange rates were assigned to two populations of protons present in different flexible regions of the protein. The addition of the substrate androstenedione and the inhibitor anastrozole lowers the H/D exchange rates of the α-helices of the enzyme when compared to the ligand-free form. Furthermore, the presence of the inhibitor anastrozole lowers exchange rate constant (k₁) for β-sheets from 0.22±0.06 min⁻¹ for the inhibitor-bound enzyme to 0.12±0.02 min⁻¹ for the free protein. Dynamics effects localised in helix F were studied by time resolved fluorescence. The data demonstrate that the fluorescence lifetime component associated to Trp224 emission undergoes a shift toward longer lifetimes (from ≈5.0 to ≈5.5 ns) when the substrate or the inhibitor are present, suggesting slower dynamics in the presence of ligands. Together the results are consistent with different degrees of flexibility of the access channel and therefore different conformations adopted by the enzyme in the free, substrate- and inhibitor-bound forms.

Investigation of fluorinated and bifunctionalized 3-phenylchroman-4-one (isoflavanone) aromatase inhibitors

Fluorinated isoflavanones and bifunctionalized isoflavanones were synthesized through a one-step gold(I)-catalyzed annulation reaction. These compounds were evaluated for their in vitro inhibitory activities against aromatase in a fluorescence-based enzymatic assay. Selected compounds were tested for their anti-proliferative effects on human breast cancer cell line MCF-7. Compounds 6-methoxy-3-(pyridin-3-yl)chroman-4-one (3c) and 6-fluoro-3-(pyridin-3-yl)chroman-4-one (3e) were identified as the most potent aromatase inhibitors with IC₅₀ values of 2.5 μM and 0.8 μM. Therefore, these compounds have great potential for the development of pharmaceutical agents against breast cancer.

Use of density functional theory in drug metabolism studies

INTRODUCTION

The cytochrome P450 enzymes (CYPs) metabolize many drug compounds. They catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be generated. Density functional theory (DFT) has, over the past decade, been shown to be a powerful tool to rationalize and predict the possible metabolites generated by the CYPs as well as other drug-metabolizing enzymes.

AREAS COVERED

We review applications of DFT on reactions performed by the CYPs and other drug-metabolizing enzymes able to perform oxidation reactions, with an emphasis on predicting which metabolites are produced. We also cover calculations of binding energies for complexes in which the ligands interact directly with the heme iron atom.

EXPERT OPINION

DFT is a useful tool for prediction of the site of metabolism. The use of small models of the enzymes work surprisingly well for most CYP isoforms. This is probably due to the fact that the binding of the substrates is not the major determinant. When binding of the substrate plays a significant role, the well-known issue of determining the free energy of binding is the challenge. How approaches taking the protein environment into account, like docking, MD and QM/MM, can be used are discussed.

Identifying cytochrome p450 functional networks and their allosteric regulatory elements

Cytochrome P450 (CYP) enzymes play key roles in drug metabolism and adverse drug-drug interactions. Despite tremendous efforts in the past decades, essential questions regarding the function and activity of CYPs remain unanswered. Here, we used a combination of sequence-based co-evolutionary analysis and structure-based anisotropic thermal diffusion (ATD) molecular dynamics simulations to detect allosteric networks of amino acid residues and characterize their biological and molecular functions. We investigated four CYP subfamilies (CYP1A, CYP2D, CYP2C, and CYP3A) that are involved in 90% of all metabolic drug transformations and identified four amino acid interaction networks associated with specific CYP functionalities, i.e., membrane binding, heme binding, catalytic activity, and dimerization. Interestingly, we did not detect any co-evolved substrate-binding network, suggesting that substrate recognition is specific for each subfamily. Analysis of the membrane binding networks revealed that different CYP proteins adopt different membrane-bound orientations, consistent with the differing substrate preference for each isoform. The catalytic networks were associated with conservation of catalytic function among CYP isoforms, whereas the dimerization network was specific to different CYP isoforms. We further applied low-temperature ATD simulations to verify proposed allosteric sites associated with the heme-binding network and their role in regulating metabolic fate. Our approach allowed for a broad characterization of CYP properties, such as membrane interactions, catalytic mechanisms, dimerization, and linking these to groups of residues that can serve as allosteric regulators. The presented combined co-evolutionary analysis and ATD simulation approach is also generally applicable to other biological systems where allostery plays a role.

Bioactivation of Aromatic Amines by Human CYP2W1, An Orphan Cytochrome P450 Enzyme

The human genome contains approximately 13 orphan cytochrome P450 (P450, CYP) genes, of which the apparent function or substrate has not been identified. However, they seem to possess their own biological relevance in some tissues or developmental stages. Here, we characterized the heterologously expressed CYP2W1, an orphan P450 enzyme. The recombinant CYP2W1 protein containing a 6 × (His) -tag at Nterminus has been expressed in Escherichia coli and purified. Expression level of CYP2W1 holoenzyme was around 500 nmol P450 holoenzyme per liter culture medium. The reduced CO difference spectrum of CYP2W1 showed a maximum absorption at 449 nm. CYP2W1 indicated the significant induction to bioactivate Trp-P-1, MeIQ, and IQ in E. coli DJ701 tester strain. However, the bioactivation of B[α]P, and NNK by CYP2W1 was relatively low. The model structure of CYP2W1 suggested the characteristic P450 folds with the lengths and orientations of the individual secondary elements. The F-G loop is situated on the distal side of heme to accommodate the flexibility of active site of CYP2W1. These studies can provide useful information for the finding of its biological roles and structure-function relationships of an orphan CYP2W1 enzyme.

Prediction of amino acid residues participated in substrate recognition by cytochrome P450 subfamilies with broad substrate specificity

Cytochromes P450 comprise a large superfamily and several of their isoforms play a crucial role in metabolism of xenobiotics, including drugs. Although these enzymes demonstrate broad and cross-substrate specificity, different cytochrome P450 subfamilies exhibit certain selectivity for some types of substrates. Analysis of amino acid residues of the active sites of six cytochrome subfamilies (CYP1А, CYP2А, CYP2С, CYP2D, CYP2E and CYP3А) enables to define subfamily-specific patterns that consist of four residues. These residues are located on the periphery of the active sites of these cytochromes. We suggest that they can form a primary binding site at the entrance to the active site, defining cytochrome substrate recognition.

Genetic variations of human CYP2D6 in the Chinese Han population

Aim: The purpose of this study was to determine the genetic polymorphisms of the CYP2D6 gene and to elucidate the allele distribution pattern in the Chinese Han population. Materials & methods: We used PCR and bidirectional sequencing methods to analyze all nine exons of the CYP2D6 gene in 2129 unrelated, healthy Chinese Han subjects from two geographical locations in China: the northern and southern regions. Results: In total, 165 mutated sites were detected in 2129 participants, of which 67 sites were reported for the first time. Among these novel mutation sites, 22 were nonsynonymous and 12 were named as novel alleles (*87–*93, *94A, *94B and *95–*98) by the Human CYP Allele Nomenclature Committee. In addition, 29 previously reported alleles and 84 genotypes were also detected in 1954 volunteers. Functional prediction of novel variants revealed that eight variants might have a deleterious effect on CYP2D6. Linkage disequilibrium analysis and tagSNP selection were performed separately. By using these methods, distinct differences were found between the two regions. Conclusion: This study provides the most comprehensive data concerning CYP2D6 polymorphisms in the Chinese Han population to date and increases the number of known alleles; these findings may greatly contribute to the development of personalized medicine for the Chinese Han population.

Original submitted 1 January 2013; Revision submitted 14 August 2013

Quantum mechanical modeling: a tool for the understanding of enzyme reactions

Most enzyme reactions involve formation and cleavage of covalent bonds, while electrostatic effects, as well as dynamics of the active site and surrounding protein regions, may also be crucial. Accordingly, special computational methods are needed to provide an adequate description, which combine quantum mechanics for the reactive region with molecular mechanics and molecular dynamics describing the environment and dynamic effects, respectively. In this review we intend to give an overview to non-specialists on various enzyme models as well as established computational methods and describe applications to some specific cases. For the treatment of various enzyme mechanisms, special approaches are often needed to obtain results, which adequately refer to experimental data. As a result of the spectacular progress in the last two decades, most enzyme reactions can be quite precisely treated by various computational methods.

Metoclopramide is metabolized by CYP2D6 and is a reversible inhibitor, but not inactivator, of CYP2D6

1. Metoclopramide is a widely used clinical drug in a variety of medical settings with rare acute dystonic events reported. The aim of this study was to assess a previous report of inactivation of CYP2D6 by metoclopramide, to determine the contribution of various CYPs to metoclopramide metabolism, and to identify the mono-oxygenated products of metoclopramide metabolism. 2. Metoclopramide interacted with CYP2D6 with Type I binding and a Ks value of 9.56 ± 1.09 µM. CYP2D6 was the major metabolizer of metoclopramide and the two major products were N-deethylation of the diethyl amine and N-hydroxylation on the phenyl ring amine. CYPs 1A2, 2C9, 2C19, and 3A4 also metabolized metoclopramide. 3. While reversible inhibition of CYP2D6 was noted, CYP2D6 inactivation by metoclopramide was not observed under conditions of varying concentration or varying time using Supersomes(TM) or pooled human liver microsomes. 4. The major metabolites of metoclopramide were N-hydroxylation and N-deethylation formed most efficiently by CYP2D6 but also formed by all CYPs examined. Also, while metoclopramide is metabolized primarily by CYP2D6, it is not a mechanism-based inactivator of CYP2D6 in vitro.

Molecular dynamics simulations give insight into the conformational change, complex formation, and electron transfer pathway for cytochrome P450 reductase

Cytochrome P450 reductase (CYPOR) undergoes a large conformational change to allow for an electron transfer to a redox partner to take place. After an internal electron transfer over its cofactors, it opens up to facilitate the interaction and electron transfer with a cytochrome P450. The open conformation appears difficult to crystallize. Therefore, a model of a human CYPOR in the open conformation was constructed to be able to investigate the stability and conformational change of this protein by means of molecular dynamics simulations. Since the role of the protein is to provide electrons to a redox partner, the interactions with cytochrome P450 2D6 (2D6) were investigated and a possible complex structure is suggested. Additionally, electron pathway calculations with a newly written program were performed to investigate which amino acids relay the electrons from the FMN cofactor of CYPOR to the HEME of 2D6. Several possible interacting amino acids in the complex, as well as a possible electron transfer pathway were identified and open the way for further investigation by site directed mutagenesis studies.

Functional influence of human CYP2D6 allelic variations: P34S, E418K, S486T, and R296C

CYP2D6 is responsible for the oxidative metabolism of 20-25 % of clinical drugs and its genetic polymorphisms can significantly influence the drug metabolism. In this study, we analyzed the functional activities of four nonsynonymous single nucleotide polymorphisms from CYP2D6*52 allele, which were recently found, and one found frequently in CYP2D6 alleles. Recombinant variant enzymes of E418K, S486T, and R296C were successfully expressed in Escherichia coli and purified. However, a CYP holoenzyme spectrum of P34S variant was not detected in E. coli whole cell level. Structural analysis indicated that P34S mutation seemed to perturb a highly conserved proline-rich N-terminus of CYP2D6. Steady state kinetic analyses showed the significant reductions of enzymatic activities in E418K and R296C variants. In the case of bufuralol 1'-hydroxylation, a novel mutant, E418K, showed 32 % decrease in catalytic efficiency (k cat/K m) mainly due to the decrease of k cat value. R296C showed much greater reduction in the catalytic efficiency (9 % of wild-type) due to both of a decrease of k cat value and an increase of K m value. In the case of dextromethorphan O-demethylation, E418K showed both of a decrease of k cat value and an increase K m value to result in ~43 % reduction of catalytic efficiency. A highly decreased catalytic efficiency (~6 % of wild-type) in the mutant of R296C also was observed mainly due to the dramatic change of k cat value of dextromethorphan O-demethylation. These results suggested that individuals carrying these allelic variants are likely to have the altered metabolic abilities of many clinical drugs therefore, these polymorphisms of CYP2D6 should be much concerned for reliable drug treatment.

Prediction of cytochrome P450 xenobiotic metabolism: tethered docking and reactivity derived from ligand molecular orbital analysis

Metabolism of xenobiotic and endogenous compounds is frequently complex, not completely elucidated, and therefore often ambiguous. The prediction of sites of metabolism (SoM) can be particularly helpful as a first step toward the identification of metabolites, a process especially relevant to drug discovery. This paper describes a reactivity approach for predicting SoM whereby reactivity is derived directly from the ground state ligand molecular orbital analysis, calculated using Density Functional Theory, using a novel implementation of the average local ionization energy. Thus each potential SoM is sampled in the context of the whole ligand, in contrast to other popular approaches where activation energies are calculated for a predefined database of molecular fragments and assigned to matching moieties in a query ligand. In addition, one of the first descriptions of molecular dynamics of cytochrome P450 (CYP) isoforms 3A4, 2D6, and 2C9 in their Compound I state is reported, and, from the representative protein structures obtained, an analysis and evaluation of various docking approaches using GOLD is performed. In particular, a covalent docking approach is described coupled with the modeling of important electrostatic interactions between CYP and ligand using spherical constraints. Combining the docking and reactivity results, obtained using standard functionality from common docking and quantum chemical applications, enables a SoM to be identified in the top 2 predictions for 75%, 80%, and 78% of the data sets for 3A4, 2D6, and 2C9, respectively, results that are accessible and competitive with other recently published prediction tools.

In vitro approaches to investigate cytochrome P450 activities: update on current status and their applicability

INTRODUCTION

Cytochromes P450 (CYPs) play a central role in the Phase I metabolism of drugs and other xenobiotics. It is estimated that CYPs can metabolize up to two-thirds of drugs present in humans. Over the past two decades, there have been numerous advances in in vitro methodologies to characterize drug metabolism and interaction involving CYPs.

AREAS COVERED

This review focuses on the use of in vitro methodologies to examine CYPs' role in drug metabolism and interaction. There is an emphasis on their current development, applicability, advantages and limitations as well as the use of in silico approaches in complementing and supporting in vitro data. The article also highlights the challenges in extrapolating in vitro data to in vivo situations.

EXPERT OPINION

Advances in in vitro methodologies have been made such that data can be used for in vivo prediction with comfortable degree of confidence. Improved assay designs and analytical techniques have permitted development of miniaturized assay format and automated system with improved sensitivity and throughput capacity. High-quality experimental designs and scientifically rigorous assessment/validation protocols remain crucial in developing reliable and robust in vitro models. With continued progress made in the field, in vitro methodologies will continually be employed in evaluating CYP activities in pharmaceutical industries and laboratories.

Modelling species selectivity in rat and human cytochrome P450 2D enzymes

Updated models of the Rat Cytochrome P450 2D enzymes are produced based on the recent x-ray structures of the Human P450 2D6 enzyme both with and without a ligand bound. The differences in species selectivity between the epimers quinine and quinidine are rationalised using these models and the results are discussed with regard to previous studies. A close approach to the heme is not observed in this study. The x-ray structure of the enzyme with a ligand bound is shown to be a better model for explaining the observed experimental binding of quinine and quinidine. Hence models with larger closed binding sites are recommended for comparative docking studies. This is consistent with molecular recognition in Cytochrome P450 enzymes being the result of a number of non-specific interactions in a large binding site.

Structural diversity of eukaryotic membrane cytochrome p450s

X-ray crystal structures are available for 29 eukaryotic microsomal, chloroplast, or mitochondrial cytochrome P450s, including two non-monooxygenase P450s. These structures provide a basis for understanding structure-function relations that underlie their distinct catalytic activities. Moreover, structural plasticity has been characterized for individual P450s that aids in understanding substrate binding in P450s that mediate drug clearance.

Cytochrome P450 enzymes in drug metabolism: regulation of gene expression, enzyme activities, and impact of genetic variation

Cytochromes P450 (CYP) are a major source of variability in drug pharmacokinetics and response. Of 57 putatively functional human CYPs only about a dozen enzymes, belonging to the CYP1, 2, and 3 families, are responsible for the biotransformation of most foreign substances including 70-80% of all drugs in clinical use. The highest expressed forms in liver are CYPs 3A4, 2C9, 2C8, 2E1, and 1A2, while 2A6, 2D6, 2B6, 2C19 and 3A5 are less abundant and CYPs 2J2, 1A1, and 1B1 are mainly expressed extrahepatically. Expression of each CYP is influenced by a unique combination of mechanisms and factors including genetic polymorphisms, induction by xenobiotics, regulation by cytokines, hormones and during disease states, as well as sex, age, and others. Multiallelic genetic polymorphisms, which strongly depend on ethnicity, play a major role for the function of CYPs 2D6, 2C19, 2C9, 2B6, 3A5 and 2A6, and lead to distinct pharmacogenetic phenotypes termed as poor, intermediate, extensive, and ultrarapid metabolizers. For these CYPs, the evidence for clinical significance regarding adverse drug reactions (ADRs), drug efficacy and dose requirement is rapidly growing. Polymorphisms in CYPs 1A1, 1A2, 2C8, 2E1, 2J2, and 3A4 are generally less predictive, but new data on CYP3A4 show that predictive variants exist and that additional variants in regulatory genes or in NADPH:cytochrome P450 oxidoreductase (POR) can have an influence. Here we review the recent progress on drug metabolism activity profiles, interindividual variability and regulation of expression, and the functional and clinical impact of genetic variation in drug metabolizing P450s.

The Effect of CYP2B6, CYP2D6, and CYP3A4 Alleles on Methadone Binding: A Molecular Docking Study

Current methadone maintenance therapy (MMT) is yet to ensure 100% successful treatment as the optimum dosage has yet to be determined. Overdose leads to death while lower dose causes the opioid withdrawal effect. Single-nucleotide polymorphisms (SNP) in cytochrome P450s (CYPs), the methadone metabolizers, have been showen to be the main factor for the interindividual variability of methadone clinical effects. In this study, we investigated the effect of SNPs in three major methadone metabolizers (CYP2B6, CYP2D6, and CYP3A4) on methadone binding affinity. Results showed thatCYP2B6*11,CYP2B6*12,CYP2B6*18, andCYP3A4*12have significantly higher binding affinity toR-methadone compared to wild type.S-methadone has higher binding affinity inCYP3A4*3,CYP3A4*11, andCYP3A4*12compared to wild type.R-methadone was shown to be the active form of methadone; thus individuals with CYP alleles that binds better toR-methadone will have higher methadone metabolism rate. Therefore, a higher dosage of methadone is necessary to obtain the opiate effect compared to a normal individual and vice versa. These results provide an initial prediction on methadone metabolism rate for individuals with mutant type CYP which enables prescription of optimum methadone dosage for individuals with CYP alleles.

Prediction of Cytochrome P450 Profiles of Environmental Chemicals with QSAR Models Built from Drug-like Molecules

The human cytochrome P450 (CYP) enzyme family is involved in the biotransformation of many xenobiotics. As part of the U.S. Tox21 Phase I effort, we profiled the CYP activity of approximately three thousand compounds, primarily those of environmental concern, against human CYP1A2, CYP2C19, CYP2C9, CYP2D6, and CYP3A4 isoforms in a quantitative high throughput screening (qHTS) format. In order to evaluate the extent to which computational models built from a drug-like library screened in these five CYP assays under the same conditions can accurately predict the outcome of an environmental compound library, five support vector machines (SVM) models built from over 17,000 drug-like compounds were challenged to predict the CYP activities of the Tox21 compound collection. Although a large fraction of the test compounds fall outside of the applicability domain (AD) of the models, as measured by k-nearest neighbor (k-NN) similarities, the predictions were largely accurate for CYP1A2, CYP2C9, and CYP3A4 ioszymes with area under the receiver operator characteristic curves (AUC-ROC) ranging between 0.82 and 0.84. The lower predictive power of the CYP2C19 model (AUC-ROC = 0.76) is caused by experimental errors and that of the CYP2D6 model (AUC-ROC = 0.76) can be rescued by rebalancing the training data. Our results demonstrate that decomposing molecules into atom types enhanced the coverage of the AD and that computational models built from drug-like molecules can be used to predict the ability of non-drug like compounds to interact with these CYPs.

Rationalization of stereospecific binding of propranolol to cytochrome P450 2D6 by free energy calculations

Abstract

Cytochrome P450 2D6 is a major drug-metabolising enzyme with a wide substrate range. A single-point mutation introduced in this enzyme induces stereoselective binding of R and S-propranolol whereas the wild type has no preference. The system has previously been studied both experimentally and computationally (de Graaf et al. in Eur Biophys J 36:589–599, 2007a). The in silico study reported hysteresis and significant deviations from closure of thermodynamic cycles, probably because of lack of sampling. Here, we focus on the effect of prolonged simulation time and enhanced sampling methods, such as Hamiltonian replica exchange, to reduce these problems and to improve the precision of free energy calculations. Finally we rationalize the results at a molecular level and compare data with experimental findings and previously estimated free energies.

Graphical Abstract

Propranolol (PPD) binds stereospecifically to cytochrome P450 2D6, if the mutation F483A is introduced, by reducing binding affinity to R-PPD. This observation has previously been studied by molecular dynamics, by use of several different thermodynamic cycles. Previous computational results suffered from internal inconsistencies, probably because of insufficient sampling. Here we use prolonged simulation time and Hamiltonian replica exchange to overcome the inconsistencies, and give a possible rationalization of the phenomenon at a molecular level.

In silico and intuitive predictions of CYP46A1 inhibition by marketed drugs with subsequent enzyme crystallization in complex with fluvoxamine

Cytochrome P450 46A1 (cholesterol 24-hydroxylase) is an important brain enzyme that may be inhibited by structurally distinct pharmaceutical agents both in vitro and in vivo. To identify additional inhibitors of CYP46A1 among U.S. Food and Drug Administration-approved therapeutic agents, we used in silico and intuitive predictions and evaluated some of the predicted binders in the enzyme and spectral binding assays. We tested a total of 298 marketed drugs for the inhibition of CYP46A1-mediated cholesterol hydroxylation in vitro and found that 13 of them reduce CYP46A1 activity by >50%. Of these 13 inhibitors, 7 elicited a spectral response in CYP46A1 with apparent spectral K(d) values in a low micromolar range. One of the identified tight binders, the widely used antidepressant fluvoxamine, was cocrystallized with CYP46A1. The structure of this complex was determined at a 2.5 Å resolution and revealed the details of drug binding to the CYP46A1 active site. The NH(2)-containing arm of the Y-shaped fluvoxamine coordinates the CYP46A1 heme iron, whereas the methoxy-containing arm points away from the heme group and has multiple hydrophobic interactions with aliphatic amino acid residues. The CF(3)-phenyl ring faces the entrance to the substrate access channel and has contacts with the aromatic side chains. The crystal structure suggests that only certain drug conformers can enter the P450 substrate access channel and reach the active site. Once inside the active site, the conformer probably further adjusts its configuration and elicits the movement of the protein side chains.