The structure of human cytochrome P450 2C9 complexed with flurbiprofen at 2.0-A resolution

Multilevel superposition for deciphering the conformational variability of protein ensembles

The dynamics and variability of protein conformations are directly linked to their functions. Many comparative studies of X-ray protein structures have been conducted to elucidate the relevant conformational changes, dynamics and heterogeneity. The rapid increase in the number of experimentally determined structures has made comparison an effective tool for investigating protein structures. For example, it is now possible to compare structural ensembles formed by enzyme species, variants or the type of ligands bound to them. In this study, the author developed a multilevel model for estimating two covariance matrices that represent inter- and intra-ensemble variability in the Cartesian coordinate space. Principal component analysis using the two estimated covariance matrices identified the inter-/intra-enzyme variabilities, which seemed to be important for the enzyme functions, with the illustrative examples of cytochrome P450 family 2 enzymes and class A $\beta$-lactamases. In P450, in which each enzyme has its own active site of a distinct size, an active-site motion shared universally between the enzymes was captured as the first principal mode of the intra-enzyme covariance matrix. In this case, the method was useful for understanding the conformational variability after adjusting for the differences between enzyme sizes. The developed method is advantageous in small ensemble-size problems and hence promising for use in comparative studies on experimentally determined structures where ensemble sizes are smaller than those generated, for example, by molecular dynamics simulations.

The molecular basis of dapsone activation of CYP2C9-catalyzed nonsteroidal anti-inflammatory drug oxidation

Positive heterotropic cooperativity, or "activation," results in an instantaneous increase in enzyme activity in the absence of an increase in protein expression. Thus, cytochrome P450 (CYP) enzyme activation presents as a potential drug-drug interaction mechanism. It has been demonstrated previously that dapsone activates the CYP2C9-catalyzed oxidation of a number of nonsteroidal anti-inflammatory drugs in vitro. Here, we conducted molecular dynamics simulations (MDS) together with enzyme kinetic investigations and site-directed mutagenesis to elucidate the molecular basis of the activation of CYP2C9-catalyzed S-flurbiprofen 4'-hydroxylation and S-naproxen O-demethylation by dapsone. Supplementation of incubations of recombinant CYP2C9 with dapsone increased the catalytic efficiency of flurbiprofen and naproxen oxidation by 2.3- and 16.5-fold, respectively. MDS demonstrated that activation arises predominantly from aromatic interactions between the substrate, dapsone, and the phenyl rings of Phe114 and Phe476 within a common binding domain of the CYP2C9 active site, rather than involvement of a distinct effector site. Mutagenesis of Phe114 and Phe476 abrogated flurbiprofen and naproxen oxidation, and MDS and kinetic studies with the CYP2C9 mutants further identified a pivotal role of Phe476 in dapsone activation. MDS additionally showed that aromatic stacking interactions between two molecules of naproxen are necessary for binding in a catalytically favorable orientation. In contrast to flurbiprofen and naproxen, dapsone did not activate the 4'-hydroxylation of diclofenac, suggesting that the CYP2C9 active site favors cooperative binding of nonsteroidal anti-inflammatory drugs with a planar or near-planar geometry. More generally, the work confirms the utility of MDS for investigating ligand binding in CYP enzymes.

Metabolism of testosterone and progesterone by cytochrome P450 2C19 allelic variants

CYP2C19 is a member of the human microsomal cytochrome P450 (CYP). Significant variation in CYP2C19 levels and activity can be attributed to polymorphisms in this gene. Wildtype CYP2C19 and 13 mutants (CYP2C19.1B, CYP2C19.5A, CYP2C19.5B, CYP2C19.6, CYP2C19.8, CYP2C19.9, CYP2C19.10, CYP2C19.11, CYP2C19.13, CYP2C19.16, CYP2C19.19, CYP2C19.23, CYP2C19.30, and CYP2C19.33) were coexpressed with NADPH-cytochrome P450 reductase in Escherichia coli. Hydroxylase activity toward testosterone and progesterone was also examined. Ten CYP2C19 variants showed Soret peaks (450 nm) typical of P450 in the reduced CO-difference spectra. CYP2C19.11 and CYP2C19.23 showed higher testosterone 11α, 16α-/17- and progesterone 6β-,21-,16α-/17α-hydroxylase activities than CYP2C19.1B. CYP2C19.6, CYP2C19.16, CYP2C19.19, and CYP2C19.30 showed lower activity than CYP2C19.1B. CYP2C19.9, CYP2C19.10. CYP2C19.13, and CYP2C19.33 showed different hydroxylation activities than CYP2C19.1B. These results indicated that CYP2C19 variants have very different substrate specificities for testosterone and progesterone.

Human Cytochrome P450 1, 2, 3 Families as Pharmacogenes with Emphases on Their Antimalarial and Antituberculosis Drugs and Prevalent African Alleles

Precision medicine gives individuals tailored medical treatment, with the genotype determining the therapeutic strategy, the appropriate dosage, and the likelihood of benefit or toxicity. Cytochrome P450 (CYP) enzyme families 1, 2, and 3 play a pivotal role in eliminating most drugs. Factors that affect CYP function and expression have a major impact on treatment outcomes. Therefore, polymorphisms of these enzymes result in alleles with diverse enzymatic activity and drug metabolism phenotypes. Africa has the highest CYP genetic diversity and also the highest burden of malaria and tuberculosis, and this review presents current general information on CYP enzymes together with variation data concerning antimalarial and antituberculosis drugs, while focusing on the first three CYP families. Afrocentric alleles such as CYP2A6*17, CYP2A6*23, CYP2A6*25, CYP2A6*28, CYP2B6*6, CYP2B6*18, CYP2C8*2, CYP2C9*5, CYP2C9*8, CYP2C9*9, CYP2C19*9, CYP2C19*13, CYP2C19*15, CYP2D6*2, CYP2D6*17, CYP2D6*29, and CYP3A4*15 are implicated in diverse metabolic phenotypes of different antimalarials such as artesunate, mefloquine, quinine, primaquine, and chloroquine. Moreover, CYP3A4, CYP1A1, CYP2C8, CYP2C18, CYP2C19, CYP2J2, and CYP1B1 are implicated in the metabolism of some second-line antituberculosis drugs such as bedaquiline and linezolid. Drug-drug interactions, induction/inhibition, and enzyme polymorphisms that influence the metabolism of antituberculosis, antimalarial, and other drugs, are explored. Moreover, a mapping of Afrocentric missense mutations to CYP structures and a documentation of their known effects provided structural insights, as understanding the mechanism of action of these enzymes and how the different alleles influence enzyme function is invaluable to the advancement of precision medicine.

Na,K-ATPase: A murzyme facilitating thermodynamic equilibriums at the membrane-interface

The redox metabolic paradigm of murburn concept advocates that diffusible reactive species (DRS, particularly oxygen-centric radicals) are mainstays of physiology, and not mere pathological manifestations. The murburn purview of cellular function also integrates the essential principles of bioenergetics, thermogenesis, homeostasis, electrophysiology, and coherence. In this context, any enzyme that generates/modulates/utilizes/sustains DRS functionality is called a murzyme. We have demonstrated that several water-soluble (peroxidases, lactate dehydrogenase, hemogoblin, etc.) and membrane-embedded (Complexes I-V in mitochondria, Photosystems I/II in chloroplasts, rhodopsin/transducin in rod cells, etc.) proteins serve as murzymes. The membrane protein of Na,K-ATPase (NKA, also known as sodium-potassium pump) is the focus of this article, owing to its centrality in neuro-cardio-musculo electrophysiology. Herein, via a series of critical queries starting from the geometric/spatio-temporal considerations of diffusion/mass transfer of solutes in cells to an update on structural/distributional features of NKA in diverse cellular systems, and from various mechanistic aspects of ion-transport (thermodynamics, osmoregulation, evolutionary dictates, etc.) to assays/explanations of inhibitory principles like cardiotonic steroids (CTS), we first highlight some unresolved problems in the field. Thereafter, we propose and apply a minimalist murburn model of trans-membrane ion-differentiation by NKA to address the physiological inhibitory effects of trans-dermal peptide, lithium ion, volatile anesthetics, confirmed interfacial DRS + proton modulators like nitrophenolics and unsaturated fatty acid, and the diverse classes of molecules like CTS, arginine, oximes, etc. These explanations find a pan-systemic connectivity with the inhibitions/uncouplings of other membrane proteins in cells.

Identification and drug metabolic characterization of four new CYP2C9 variants CYP2C9*72-*75 in the Chinese Han population

Cytochrome 2C9 (CYP2C9), one of the most important drug metabolic enzymes in the human hepatic P450 superfamily, is required for the metabolism of 15% of clinical drugs. Similar to other CYP2C family members, CYP2C9 gene has a high genetic polymorphism which can cause significant racial and inter-individual differences in drug metabolic activity. To better understand the genetic distribution pattern of CYP2C9 in the Chinese Han population, 931 individuals were recruited and used for the genotyping in this study. As a result, seven synonymous and 14 non-synonymous variations were identified, of which 4 missense variants were designated as new alleles CYP2C9*72, *73, *74 and *75, resulting in the amino acid substitutions of A149V, R150C, Q214H and N418T, respectively. When expressed in insect cell microsomes, all four variants exhibited comparable protein expression levels to that of the wild-type CYP2C9 enzyme. However, drug metabolic activity analysis revealed that these variants exhibited significantly decreased catalytic activities toward three CYP2C9 specific probe drugs, as compared with that of the wild-type enzyme. These data indicate that the amino acid substitution in newly designated variants can cause reduced function of the enzyme and its clinical significance still needs further investigation in the future.

The protein conformational basis of isoflavone biosynthesis

Isoflavonoids play important roles in plant defense and also exhibit a range of mammalian health-promoting activities. Their biosynthesis is initiated by two enzymes with unusual catalytic activities; 2-hydroxyisoflavanone synthase (2-HIS), a membrane-bound cytochrome P450 catalyzing a coupled aryl-ring migration and hydroxylation, and 2-hydroxyisoflavanone dehydratase (2-HID), a member of a large carboxylesterase family that paradoxically catalyzes dehydration of 2-hydroxyisoflavanones to isoflavone. Here we report the crystal structures of 2-HIS from Medicago truncatula and 2-HID from Pueraria lobata. The 2-HIS structure reveals a unique cytochrome P450 conformation and heme and substrate binding mode that facilitate the coupled aryl-ring migration and hydroxylation reactions. The 2-HID structure reveals the active site architecture and putative catalytic residues for the dual dehydratase and carboxylesterase activities. Mutagenesis studies revealed key residues involved in substrate binding and specificity. Understanding the structural basis of isoflavone biosynthesis will facilitate the engineering of new bioactive isoflavonoids.

The structure and function of two isoflavone biosynthetic enzymes are reported revealing a novel cytochrome P450 conformation and identification of key residues for dual dehydratase and carboxylesterase activities

Designing Artificial Metalloenzymes by Tuning of the Environment beyond the Primary Coordination Sphere

Metalloenzymes catalyze a variety of reactions using a limited number of natural amino acids and metallocofactors. Therefore, the environment beyond the primary coordination sphere must play an important role in both conferring and tuning their phenomenal catalytic properties, enabling active sites with otherwise similar primary coordination environments to perform a diverse array of biological functions. However, since the interactions beyond the primary coordination sphere are numerous and weak, it has been difficult to pinpoint structural features responsible for the tuning of activities of native enzymes. Designing artificial metalloenzymes (ArMs) offers an excellent basis to elucidate the roles of these interactions and to further develop practical biological catalysts. In this review, we highlight how the secondary coordination spheres of ArMs influence metal binding and catalysis, with particular focus on the use of native protein scaffolds as templates for the design of ArMs by either rational design aided by computational modeling, directed evolution, or a combination of both approaches. In describing successes in designing heme, nonheme Fe, and Cu metalloenzymes, heteronuclear metalloenzymes containing heme, and those ArMs containing other metal centers (including those with non-native metal ions and metallocofactors), we have summarized insights gained on how careful controls of the interactions in the secondary coordination sphere, including hydrophobic and hydrogen bonding interactions, allow the generation and tuning of these respective systems to approach, rival, and, in a few cases, exceed those of native enzymes. We have also provided an outlook on the remaining challenges in the field and future directions that will allow for a deeper understanding of the secondary coordination sphere a deeper understanding of the secondary coordintion sphere to be gained, and in turn to guide the design of a broader and more efficient variety of ArMs.

TEMPO-mediated late stage photochemical hydroxylation of biaryl sulfonium salts

The late stage photochemical hydroxylation of biaryl sulfonium salts was enabled with a TEMPO derivative as a simple oxygen source, in metal free conditions. The scope and mechanism of this exceptionally simple synthetic method, which constructs important arylated phenols from aromatic C-H bonds, are herein discussed.

Machine learning-driven identification of drugs inhibiting cytochrome P450 2C9

Cytochrome P450 2C9 (CYP2C9) is a major drug-metabolizing enzyme that represents 20% of the hepatic CYPs and is responsible for the metabolism of 15% of drugs. A general concern in drug discovery is to avoid the inhibition of CYP leading to toxic drug accumulation and adverse drug-drug interactions. However, the prediction of CYP inhibition remains challenging due to its complexity. We developed an original machine learning approach for the prediction of drug-like molecules inhibiting CYP2C9. We created new predictive models by integrating CYP2C9 protein structure and dynamics knowledge, an original selection of physicochemical properties of CYP2C9 inhibitors, and machine learning modeling. We tested the machine learning models on publicly available data and demonstrated that our models successfully predicted CYP2C9 inhibitors with an accuracy, sensitivity and specificity of approximately 80%. We experimentally validated the developed approach and provided the first identification of the drugs vatalanib, piriqualone, ticagrelor and cloperidone as strong inhibitors of CYP2C9 with IC values <18 μM and sertindole, asapiprant, duvelisib and dasatinib as moderate inhibitors with IC50 values between 40 and 85 μM. Vatalanib was identified as the strongest inhibitor with an IC50 value of 0.067 μM. Metabolism assays allowed the characterization of specific metabolites of abemaciclib, cloperidone, vatalanib and tarafenacin produced by CYP2C9. The obtained results demonstrate that such a strategy could improve the prediction of drug-drug interactions in clinical practice and could be utilized to prioritize drug candidates in drug discovery pipelines.

Massively parallel characterization of CYP2C9 variant enzyme activity and abundance

CYP2C9 encodes a cytochrome P450 enzyme responsible for metabolizing up to 15% of small molecule drugs, and CYP2C9 variants can alter the safety and efficacy of these therapeutics. In particular, the anti-coagulant warfarin is prescribed to over 15 million people annually and polymorphisms in CYP2C9 can affect individual drug response and lead to an increased risk of hemorrhage. We developed click-seq, a pooled yeast-based activity assay, to test thousands of variants. Using click-seq, we measured the activity of 6,142 missense variants in yeast. We also measured the steady-state cellular abundance of 6,370 missense variants in a human cell line by using variant abundance by massively parallel sequencing (VAMP-seq). These data revealed that almost two-thirds of CYP2C9 variants showed decreased activity and that protein abundance accounted for half of the variation in CYP2C9 function. We also measured activity scores for 319 previously unannotated human variants, many of which may have clinical relevance.

Descriptors of Cytochrome Inhibitors and Useful Machine Learning Based Methods for the Design of Safer Drugs

Roughly 2.8% of annual hospitalizations are a result of adverse drug interactions in the United States, representing more than 245,000 hospitalizations. Drug-drug interactions commonly arise from major cytochrome P450 (CYP) inhibition. Various approaches are routinely employed in order to reduce the incidence of adverse interactions, such as altering drug dosing schemes and/or minimizing the number of drugs prescribed; however, often, a reduction in the number of medications cannot be achieved without impacting therapeutic outcomes. Nearly 80% of drugs fail in development due to pharmacokinetic issues, outlining the importance of examining cytochrome interactions during preclinical drug design. In this review, we examined the physiochemical and structural properties of small molecule inhibitors of CYPs 3A4, 2D6, 2C19, 2C9, and 1A2. Although CYP inhibitors tend to have distinct physiochemical properties and structural features, these descriptors alone are insufficient to predict major cytochrome inhibition probability and affinity. Machine learning based in silico approaches may be employed as a more robust and accurate way of predicting CYP inhibition. These various approaches are highlighted in the review.

Assessment of CYP2C9 Structural Models for Site of Metabolism Prediction

Structure-based prediction of a compound's potential sites of metabolism (SOMs) mediated by cytochromes P450 (CYPs) is highly advantageous in the early stage of drug discovery. However, the accuracy of the SOMs prediction can be influenced by several factors. CYP2C9 is one of the major drug-metabolizing enzymes in humans and is responsible for the metabolism of ∼13 % of clinically used drugs. In this study, we systematically evaluated the effects of protein crystal structure models, scoring functions, heme forms, conserved active-site water molecules, and protein flexibility on SOMs prediction of CYP2C9 substrates. Our results demonstrated that, on average, ChemScore and GlideScore outperformed four other scoring functions: Vina, GoldScore, ChemPLP, and ASP. The performance of the crystal structure models with pentacoordinated heme was generally superior to that of the hexacoordinated iron-oxo heme (referred to as Compound I) models. Inclusion of the conserved active-site water molecule improved the prediction accuracy of GlideScore, but reduced the accuracy of ChemScore. In addition, the effect of the conserved water on SOMs prediction was found to be dependent on the receptor model and the substrate. We further found that one of snapshots from molecular dynamics simulations on the apo form can improve the prediction accuracy when compared to the crystal structural model.

Differential inhibition of human CYP2C8 and molecular docking interactions elicited by sorafenib and its major N-oxide metabolite

The tyrosine kinase inhibitor sorafenib (SOR) is being used increasingly in combination with other anticancer agents like paclitaxel, but this increases the potential for drug toxicity. SOR inhibits several human CYPs, including CYP2C8, which is a major enzyme in the elimination of oncology drugs like paclitaxel and imatinib. It has been reported that CYP2C8 inhibition by SOR in human liver microsomes is potentiated by NADPH-dependent biotransformation. This implicates a SOR metabolite in enhanced inhibition, although the identity of that metabolite is presently unclear. The present study evaluated the capacity of the major N-oxide metabolite of SOR (SNO) to inhibit CYP2C8-dependent paclitaxel 6α-hydroxylation. The IC₅₀ of SNO against CYP2C8 activity was found to be 3.7-fold lower than that for the parent drug (14 μM versus 51 μM). In molecular docking studies, both SOR and SNO interacted with active site residues in CYP2C8, but four additional major hydrogen and halogen bonding interactions were identified between SNO and amino acids in the B-B' loop region and helixes F' and I that comprise the catalytic region of the enzyme. In contrast, the binding of both SOR and SNO to active site residues in the closely related human CYP2C9 enzyme was similar, as were the IC₅₀s determined against CYP2C9-mediated losartan oxidation. These findings suggest that the active metabolite SNO could impair the elimination of coadministered drugs that are substrates for CYP2C8, and mediate toxic adverse events, perhaps in those individuals in whom SNO is formed extensively.

Crystal Structures of Drug-Metabolizing CYPs

The complex enzyme kinetics displayed by drug-metabolizing cytochrome P450 enzymes (CYPs) (see Chapter 9 ) can, in part, be explained by an examination of their crystallographic protein structures. Fortunately, despite low sequence similarity between different families of drug-metabolizing CYPs, there exists a high degree of structural homology within the superfamily. This similarity in the protein fold allows for a direct comparison of the structural features of CYPs that contribute toward differences in substrate binding, heterotropic and homotropic cooperativity, and genetic variability in drug metabolism. In this chapter, we first provide an overview of the nomenclature and the role of structural features that are common in all CYPs. We then apply these definitions to understand the different substrate specificities and functions in the CYP3A, CYP2C, and CYP2D families of enzymes.

Engineered human CYP2C9 and its main polymorphic variants for bioelectrochemical measurements of catalytic response

Polymorphism is an important aspect in drug metabolism responsible for different individual response to drug dosage, often leading to adverse drug reactions. Here human CYP2C9 as well as its polymorphic variants CYP2C9*2 and CYP2C9*3 present in approximately 35% of the Caucasian population have been engineered by linking their gene to the one of D. vulgaris flavodoxin (FLD) that acts as regulator of the electron flow from the electrode surface to the haem. The redox properties of the immobilised proteins were investigated by cyclic voltammetry and electrocatalysis was measured in presence of the largely used anticoagulant drug S-warfarin, marker substrate for CYP2C9. Immobilisation of the CYP2C9-FLD, CYP2C9*2-FLD and CYP2C9*3-FLD on DDAB modified glassy carbon electrodes showed well defined redox couples on the oxygen-free cyclic voltammograms and mid-point potentials of all enzymes were calculated. Electrocatalysis in presence of substrate and quantification of the product formed showed lower catalytic activities for the CYP2C9*3-FLD (2.73 ± 1.07 min^-1) and CYP2C9*2-FLD (12.42 ± 2.17 min^-1) compared to the wild type CYP2C9-FLD (18.23 ± 1.29 min^-1). These differences in activity among the CYP2C9 variants are in line with the reported literature data, and this set the basis for the use of the bio-electrode for the measurement of the different catalytic responses towards drugs very relevant in therapy.

CYP2J2 Molecular Recognition: A New Axis for Therapeutic Design

Cytochrome P450 (CYP) epoxygenases are a special subset of heme-containing CYP enzymes capable of performing the epoxidation of polyunsaturated fatty acids (PUFA) and the metabolism of xenobiotics. This dual functionality positions epoxygenases along a metabolic crossroad. Therefore, structure-function studies are critical for understanding their role in bioactive oxy-lipid synthesis, drug-PUFA interactions, and for designing therapeutics that directly target the epoxygenases. To better exploit CYP epoxygenases as therapeutic targets, there is a need for improved understanding of epoxygenase structure-function. Of the characterized epoxygenases, human CYP2J2 stands out as a potential target because of its role in cardiovascular physiology. In this review, the early research on the discovery and activity of epoxygenases is contextualized to more recent advances in CYP epoxygenase enzymology with respect to PUFA and drug metabolism. Additionally, this review employs CYP2J2 epoxygenase as a model system to highlight both the seminal works and recent advances in epoxygenase enzymology. Herein we cover CYP2J2's interactions with PUFAs and xenobiotics, its tissue-specific physiological roles in diseased states, and its structural features that enable epoxygenase function. Additionally, the enumeration of research on CYP2J2 identifies the future needs for the molecular characterization of CYP2J2 to enable a new axis of therapeutic design.

Inhibitory Effect of AB-PINACA, Indazole Carboxamide Synthetic Cannabinoid, on Human Major Drug-Metabolizing Enzymes and Transporters

Indazole carboxamide synthetic cannabinoid, AB-PINACA, has been placed into Schedule I of the Controlled Substances Act by the US Drug Enforcement Administration since 2015. Despite the possibility of AB-PINACA exposure in drug abusers, the interactions between AB-PINACA and drug-metabolizing enzymes and transporters that play crucial roles in the pharmacokinetics and efficacy of various substrate drugs have not been investigated. This study was performed to investigate the inhibitory effects of AB-PINACA on eight clinically important human major cytochrome P450s (CYPs) and six uridine 5′-diphospho-glucuronosyltransferases (UGT) in human liver microsomes and the activities of six solute carrier transporters and two efflux transporters in transporter-overexpressing cells. AB-PINACA reversibly inhibited the metabolic activities of CYP2C8 (K_i, 16.9 µM), CYP2C9 (K_i, 6.7 µM), and CYP2C19 (K_i, 16.1 µM) and the transport activity of OAT3 (K_i, 8.3 µM). It exhibited time-dependent inhibition on CYP3A4 (K_i, 17.6 µM; k_inact, 0.04047 min⁻¹). Other metabolizing enzymes and transporters such as CYP1A2, CYP2A6, CYP2B6, CYP2D6, UGT1A1, UGT1A3, UGT1A4, UGT1A6, UGT1A9, UGT2B7, OAT1, OATP1B1, OATP1B3, OCT1, OCT2, P-glycoprotein, and BCRP, exhibited only weak interactions with AB-PINACA. These data suggest that AB-PINACA can cause drug-drug interactions with CYP3A4 substrates but that the significance of drug interactions between AB-PINACA and CYP2C8, CYP2C9, CYP2C19, or OAT3 substrates should be interpreted carefully.

Arginine-259 of UGT2B7 Confers UDP-Sugar Selectivity

Enzymes of the human UDP-glycosyltransferase (UGT) superfamily typically catalyze the covalent addition of the sugar moiety from a UDP-sugar cofactor to relatively low-molecular weight lipophilic compounds. Although UDP-glucuronic acid (UDP-GlcUA) is most commonly employed as the cofactor by UGT1 and UGT2 family enzymes, UGT2B7 and several other enzymes can use both UDP-GlcUA and UDP-glucose (UDP-Glc), leading to the formation of glucuronide and glucoside conjugates. An investigation of UGT2B7-catalyzed morphine glycosidation indicated that glucuronidation is the principal route of metabolism because the binding affinity of UDP-GlcUA is higher than that of UDP-Glc. Currently, it is unclear which residues in the UGT2B7 cofactor binding domain are responsible for the preferential binding of UDP-GlcUA. Here, molecular dynamics (MD) simulations were performed together with site-directed mutagenesis and enzyme kinetic studies to identify residues within the UGT2B7 binding site responsible for the selective cofactor binding. MD simulations demonstrated that Arg259, which is located within the N-terminal domain, specifically interacts with UDP-GlcUA, whereby the side chain of Arg259 H-bonds and forms a salt bridge with the carboxylate group of glucuronic acid. Consistent with the MD simulations, substitution of Arg259 with Leu resulted in the loss of morphine, 4-methylumbelliferone, and zidovudine glucuronidation activity, but morphine glucosidation was preserved. SIGNIFICANCE STATEMENT: Despite the importance of uridine diphosphate glycosyltransferase (UGT) enzymes in drug and chemical metabolism, cofactor binding interactions are incompletely understood, as is the molecular basis for preferential glucuronidation by UGT1 and UGT2 family enzymes. The study demonstrated that long timescale molecular dynamics (MD) simulations with a UGT2B7 homology model can be used to identify critical binding interactions of a UGT protein with UDP-sugar cofactors. Further, the data provide a basis for the application of MD simulations to the elucidation of UGT-aglycone interactions.

In-Silico Modeling in Drug Metabolism and Interaction: Current Strategies of Lead Discovery

BACKGROUND

Drug metabolism is a complex mechanism of human body systems to detoxify foreign particles, chemicals, and drugs through bio alterations. It involves many biochemical reactions carried out by invivo enzyme systems present in the liver, kidney, intestine, lungs, and plasma. After drug administration, it crosses several biological membranes to reach into the target site for binding and produces the therapeutic response. After that, it may undergo detoxification and excretion to get rid of the biological systems. Most of the drugs and its metabolites are excreted through kidney via urination. Some drugs and their metabolites enter into intestinal mucosa and excrete through feces. Few of the drugs enter into hepatic circulation where they go into the intestinal tract. The drug leaves the liver via the bile duct and is excreted through feces. Therefore, the study of total methodology of drug biotransformation and interactions with various targets is costly.

METHODS

To minimize time and cost, in-silico algorithms have been utilized for lead-like drug discovery. Insilico modeling is the process where a computer model with a suitable algorithm is developed to perform a controlled experiment. It involves the combination of both in-vivo and in-vitro experimentation with virtual trials, eliminating the non-significant variables from a large number of variable parameters. Whereas, the major challenge for the experimenter is the selection and validation of the preferred model, as well as precise simulation in real physiological status.

RESULTS

The present review discussed the application of in-silico models to predict absorption, distribution, metabolism, and excretion (ADME) properties of drug molecules and also access the net rate of metabolism of a compound.

CONCLUSION

It helps with the identification of enzyme isoforms; which are likely to metabolize a compound, as well as the concentration dependence of metabolism and the identification of expected metabolites. In terms of drug-drug interactions (DDIs), models have been described for the inhibition of metabolism of one compound by another, and for the compound-dependent induction of drug-metabolizing enzymes.

Heterologous Expression and Functional Characterization of Novel CYP2C9 Variants Identified in the Alaska Native People

CYP2C9 is a major form of human liver cytochrome P450 that is responsible for the oxidative metabolism of several widely used low-therapeutic index drugs, including (S)-warfarin and phenytoin. In a cohort of Alaska Native people, ultrarare or novel CYP2C9 protein variants, M1L (rs114071557), N218I (rs780801862), and P279T (rs182132442, CYP2C9*29), are expressed with higher frequencies than the well characterized CYP2C9*2 and CYP2C9*3 alleles. We report here on their relative expression in lentivirus-infected HepG2 cells and the functional characterization of purified reconstituted enzyme variants expressed in Escherichia coli toward (S)-warfarin, phenytoin, flurbiprofen, and (S)-naproxen. In the infected HepG2 cells, robust mRNA and protein expression were obtained for wild-type, N218I, and P279T variants, but as expected, the M1L variant protein was not translated in this liver-derived cell line. His-tagged wild-type protein and the N218I and P279T variants, but not M1L, expressed well in E. coli and were highly purified after affinity chromatography. Upon reconstitution with cytochrome P450 oxidoreductase and cytochrome b5, the N218I and P279T protein variants metabolized (S)-warfarin, phenytoin, flurbiprofen, and (S)-naproxen to the expected monohydroxylated or O-demethylated metabolites. Steady-state kinetic analyses revealed that the relative catalytic efficiency ratios of (S)-warfarin metabolism by the P279T and N218I variants were 87% and 24%, respectively, of wild-type CYP2C9 protein. A similar rank ordering was observed for metabolism of phenytoin, flurbiprofen, and (S)-naproxen. We conclude that carriers of the variant N218I and, especially, the M1L alleles would be at risk of exacerbated therapeutic effects from drugs that rely on CYP2C9 for their metabolic clearance. SIGNIFICANCE STATEMENT: Novel gene variants of CYP2C9-M1L, and N218I, along with P279T (CYP2C9*29)-are expressed in Alaska Native people at relatively high frequencies. In vitro characterization of their functional effects revealed that each variant confers reduced catalytic efficiency toward several substrates, including the low-therapeutic index drugs (S)-warfarin and phenytoin. These data provide the first functional information for new, common CYP2C9 variants in this understudied population. The data may help guide dose adjustments in allele carriers, thus mitigating potential healthcare disparities.

The Effect of Force-Field Parameters on Cytochrome P450-Membrane Interactions: Structure and Dynamics

The simulation of membrane proteins requires compatible protein and lipid force fields that reproduce the properties of both the protein and the lipid bilayer. Cytochrome P450 enzymes are bitopic membrane proteins with a transmembrane helical anchor and a large cytosolic globular domain that dips into the membrane. As such, they are representative and challenging examples of membrane proteins for simulations, displaying features of both peripheral and integral membrane proteins. We performed molecular dynamics simulations of three cytochrome P450 isoforms (2C9, 2C19 and 1A1) in a 2-oleoyl-1-palmitoyl-sn-glycerol-3-phosphocholine bilayer using two AMBER force field combinations: GAFF-LIPID with ff99SB for the protein, and LIPID14 with ff14SB for the protein. Comparison of the structural and dynamic properties of the proteins, the lipids and the protein-membrane interactions shows differing sensitivity of the cytochrome P450 isoforms to the choice of force field, with generally better agreement with experiment for the LIPID14 + ff14SB combination.

Exploring the biotransformation of N-(2-hydroxyphenyl)-2-propylpentanamide (an aryl valproic acid derivative) by CYP2C11, using in silico predictions and in vitro studies

OBJECTIVES

N-(2-hydroxyphenyl)-2-propylpentanamide (HO-AAVPA), a derivative of valproic acid (VPA), has been proposed as a potential anticancer agent due to its improved antiproliferative effects in some cancer cell lines. Although there is evidence that VPA is metabolized by cytochrome P450 2C11 rat isoform, HO-AAVPA CYP-mediated metabolism has not yet been fully explored. Therefore, in this work, the biotransformation of HO-AAVPA by CYP2C11 was investigated.

METHODS

Kinetic parameters and spectral interaction between HO-AAVPA and CYP were evaluated using rat liver microsomes. The participation of CYP2C11 in metabolism of HO-AAVPA was confirmed by cimetidine (CIM) inhibition assay. Docking and molecular dynamics simulations coupled to MMGBSA methods were used in theoretical study.

KEY FINDINGS

HO-AAVPA is metabolized by CYP enzymes (K_M  = 38.94 µm), yielding a hydroxylated metabolite according to its HPLC retention time (5.4 min) and MS analysis (252.2 m/z). In addition, CIM inhibition in rat liver microsomes (K_i  = 59.23 µm) confirmed that CYP2C11 is mainly involved in HO-AAVPA metabolism. Furthermore, HO-AAVPA interacts with CYP2C11 as a type I ligand. HO-AAVPA is stabilized at the CYP2C11 ligand recognition site through a map of interactions similar to other typical CYP2C11 substrates.

CONCLUSION

Therefore, rat liver CYP2C11 isoform is able to metabolize HO-AAVPA.

Murburn concept: a paradigm shift in cellular metabolism and physiology

Two decades of evidence-based exploratory pursuits in heme-flavin enzymology led to the formulation of a new biological electron/moiety transfer paradigm, called murburn concept. Murburn is a novel literary abstraction from "mured burning" or "mild unrestricted burning". This concept was invoked to explain the longstanding conundrum of maverick physiological dose responses and also applied to remodel the prevailing understanding of drug metabolism and cellular respiration. A conglomeration of simple ideas grounded in the known principles of thermodynamics and reaction chemistry, murburn concept invokes catalytic/functional roles for diffusible reactive species or radicals. Hitherto, diffusible reactive species were primarily seen as toxic agents of chaos, non-conducible to the maintenance of life-order. Since the murburn paradigm offers a distinctly different perspective for several biological phenomena, researchers holding conventional views of cellular metabolism pose a direct conflict of interests to the advancement of murburn concept. Murburn schemes are poised to integrate numerous metabolic motifs with holistic physiological outcomes; redefining pursuits in biology and medicine. To advance this agenda, I present a brief account of murburn concept and point out how redundant ideas are still advocated in some prestigious journals.

Induced fit for cytochrome P450 3A4 based on molecular dynamics

The present study aims at numerically describing to what extent substrate - enzyme complexes in solution may change over time as a natural process of conformational changes for a liganded enzyme in comparison to those movements which occur independently from substrate interaction, i.e. without a ligand. To this end, we selected structurally known pairs of liganded / unliganded CYP450 3A4 enzymes with different geometries hinting at induced fit events. We carried out molecular dynamics simulations (MD) comparing the trajectories in a “cross-over” protocol: (i) we added the ligand to the unliganded crystal form which should adopt geometries similar to the known geometry of the liganded crystal structure during MD, and – conversely – (ii) we removed the bound ligand form the known liganded complex to test if a geometry similar to the known unliganded (apo-) form can be adopted during MD. To compare continues changes we measured root means square deviations and frequencies. Results for case (i) hint at larger conformational changes required for accepting the substrate during its approach to final position – in contrast to case (ii) when mobility is fairly reduced by ligand binding (strain energy). In conclusion, a larger conformational sampling prior to ligand binding and the freezing-in (rigidity) of conformations for bound ligands can be interpreted as two conditions linked to induced-fit.

Molecular Docking as a Promising Predictive Model for Silver Nanoparticle-Mediated Inhibition of Cytochrome P450 Enzymes

Cytochrome P450 (CYP) enzymes are responsible for oxidative metabolisms of a large number of xenobiotics. In this study, we investigated interactions of silver nanoparticles (AgNPs) and silver ions (Ag⁺) with six CYP isoforms, namely, CYP1A2, CYP2C9, CYP2C19, CYP2D6, CYP2E1, and CYP3A4, within CYP-specific inhibitor-binding pockets by molecular docking and quantum mechanical (QM) calculations. The docking results revealed that the Ag₃ cluster, not Ag⁺, interacted with key amino acids of CYP2C9, CYP2C19, and CYP2D6 within a distance of about 3 Å. Moreover, the QM analysis confirmed that the amino acid residues of these CYP enzymes strongly interacted with the Ag₃ cluster, providing more insight into the mechanism of the potential inhibition of CYP enzyme activities. Interestingly, these results are consistent with previous in vitro data indicating that AgNPs inhibited activities of CYP2C and CYP2D in rat liver microsomes. It is suggested that the Ag₃ cluster is a minimal unit of AgNPs for in silico modeling. In summary, we demonstrated that molecular docking, together with QM analysis, is a promising tool to predict AgNP-mediated CYP inhibition. These methods are useful for deeper understanding of reaction mechanisms and could be used for other nanomaterials.

Differing Membrane Interactions of Two Highly Similar Drug-Metabolizing Cytochrome P450 Isoforms: CYP 2C9 and CYP 2C19

The human cytochrome P450 (CYP) 2C9 and 2C19 enzymes are two highly similar isoforms with key roles in drug metabolism. They are anchored to the endoplasmic reticulum membrane by their N-terminal transmembrane helix and interactions of their cytoplasmic globular domain with the membrane. However, their crystal structures were determined after N-terminal truncation and mutating residues in the globular domain that contact the membrane. Therefore, the CYP-membrane interactions are not structurally well-characterized and their dynamics and the influence of membrane interactions on CYP function are not well understood. We describe herein the modeling and simulation of CYP 2C9 and CYP 2C19 in a phospholipid bilayer. The simulations revealed that, despite high sequence conservation, the small sequence and structural differences between the two isoforms altered the interactions and orientations of the CYPs in the membrane bilayer. We identified residues (including K72, P73, and I99 in CYP 2C9 and E72, R73, and H99 in CYP 2C19) at the protein-membrane interface that contribute not only to the differing orientations adopted by the two isoforms in the membrane, but also to their differing substrate specificities by affecting the substrate access tunnels. Our findings provide a mechanistic interpretation of experimentally observed effects of mutagenesis on substrate selectivity.

Exploring the Chemical Space of Cytochrome P450 Inhibitors Using Integrated Physicochemical Parameters, Drug Efficiency Metrics and Decision Tree Models

The cytochrome P450s (CYPs) play a central role in the metabolism of various endogenous and exogenous compounds including drugs. CYPs are vulnerable to inhibition and induction which can lead to adverse drug reactions. Therefore, insights into the underlying mechanism of CYP450 inhibition and the estimation of overall CYP inhibitor properties might serve as valuable tools during the early phases of drug discovery. Herein, we present a large data set of inhibitors against five major metabolic CYPs (CYP1A2, CYP2C9, CYP2C19, CYP2D6 and CYP3A4) for the evaluation of important physicochemical properties and ligand efficiency metrics to define property trends across various activity levels (active, efficient and inactive). Decision tree models for CYP inhibition were developed with an accuracy >90% for both the training set and 10-folds cross validation. Overall, molecular weight (MW), hydrogen bond acceptors/donors (HBA/HBD) and lipophilicity (clogP/logPo/w) represent important physicochemical descriptors for CYP450 inhibitors. However, highly efficient CYP inhibitors show mean MW, HBA, HBD and logP values between 294.18–482.40,5.0–8.2,1–7.29 and 1.68–2.57, respectively. Our results might help in optimization of toxicological profiles associated with new chemical entities (NCEs), through a better understanding of inhibitor properties leading to CYP-mediated interactions.

Oxidation of Flavone, 5-Hydroxyflavone, and 5,7-Dihydroxyflavone to Mono-, Di-, and Tri-Hydroxyflavones by Human Cytochrome P450 Enzymes

Biologically active plant flavonoids, including 5,7-dihydroxyflavone (57diOHF, chrysin), 4',5,7-trihydroxyflavone (4'57triOHF, apigenin), and 5,6,7-trihydroxyflavone (567triOHF, baicalein), have important pharmacological and toxicological significance, e.g., antiallergic, anti-inflammatory, antioxidative, antimicrobial, and antitumorgenic properties. In order to better understand the metabolism of these flavonoids in humans, we examined the oxidation of flavone, 5-hydroxyflavone (5OHF), and 57diOHF to various products by human cytochrome P450 (P450 or CYP) and liver microsomal enzymes. Individual human P450s and liver microsomes oxidized flavone to 6-hydroxyflavone, small amounts of 5OHF, and 11 other monohydroxylated products at different rates and also produced several dihydroxylated products (including 57diOHF and 7,8-dihydroxyflavone) from flavone. We also found that 5OHF was oxidized by several P450 enzymes and human liver microsomes to 57diOHF and further to 567triOHF, but the turnover rates in these reactions were low. Interestingly, both CYP1B1.1 and 1B1.3 converted 57diOHF to 567triOHF at turnover rates (on the basis of P450 contents) of >3.0 min^-1, and CYP1A1 and 1A2 produced 567triOHF at rates of 0.51 and 0.72 min^-1, respectively. CYP2A13 and 2A6 catalyzed the oxidation of 57diOHF to 4'57triOHF at rates of 0.7 and 0.1 min^-1, respectively. Our present results show that different P450s have individual roles in oxidizing these phytochemical flavonoids and that these reactions may cause changes in their biological and toxicological properties in mammals.

Investigation on drug-binding in heme pocket of CYP2C19 with UV-visible and resonance Raman spectroscopies

Cytochrome P450 (CYP) is a class of heme-containing enzymes which mainly catalyze a monooxygenation reaction of various chemicals, and hence CYP plays a key role in the drug metabolism. Although CYP2C19 isoform is a minor hepatic CYP, it metabolizes clinically important drugs such as omeprazole and S‑mephenytoin. In this work, the interaction of purified CYP2C19 WT (CYP2C19) with seven drugs (phenytoin, S‑mephenytoin, omeprazole, lansoprazole, cimetidine, propranolol, and warfarin) was investigated using spectroscopic methods. The binding of each drug and the induced structural change in the heme distal environment were evaluated. Ferric form of CYP2C19 was revealed to contain a six-coordinate low-spin heme with a water molecule as a sixth ligand in a distal site, and the addition of each drug caused varied minor fraction of five-coordinate heme. It was suggested that the ligated water molecule was partly moved away from the heme distal environment and that the degree of water removal was dependent on the type of drugs. The effect on the coordination was varied with the studied drugs with wide variation in the dissociation constants from 2.6 μM for lansoprazole to 5400 μM for warfarin. Phenytoin and S‑mephenytoin showed that binding to CYP2C19 occurred in a stepwise manner and that the coordination of a water molecule was facilitated in the second binding step. In the ferrous CO-bound state, ν(FeCO) stretching mode was clearly observed at 471 cm^-1 in the absence of drugs. The Raman line was greatly up-shifted by omeprazole (487 cm^-1) and lansoprazole (477 cm^-1) but was minimally affected by propranolol, phenytoin, and S‑mephenytoin. These results indicate that slight chemical modification of a drug greatly affects the heme distal environments upon binding.

Pharmacogenomic Next-Generation DNA Sequencing: Lessons from the Identification and Functional Characterization of Variants of Unknown Significance in CYP2C9 and CYP2C19

CYP2C9 and CYP2C19 are highly polymorphic pharmacogenes; however, clinically actionable genetic variability in drug metabolism due to these genes has been limited to a few common alleles. The identification and functional characterization of less-common open reading frame sequence variation might help to individualize therapy with drugs that are substrates for the enzymes encoded by these genes. The present study identified seven uncharacterized variants each in CYP2C9 and CYP2C19 using next-generation sequence data for 1013 subjects, and functionally characterized the encoded proteins. Constructs were created and transiently expressed in COS-1 cells for the assay of protein concentration and enzyme activities using fluorometric substrates and liquid chromatography- tandem mass spectrometry with tolbutamide (CYP2C9) and (S)-mephenytoin (CYP2C19) as prototypic substrates. The results were compared with the SIFT, Polyphen, and Provean functional prediction software programs. Cytochrome P450 oxidoreductase (CPR) activities were also determined. Positive correlations were observed between protein content and fluorometric enzyme activity for variants of CYP2C9 (P < 0.05) and CYP2C19 (P < 0.0005). However, CYP2C9 709G>C and CYP2C19 65A>G activities were much lower than predicted based on protein content. Substrate intrinsic clearance values for CYP2C9 218C>T, 343A>C, and CYP2C19 337G>A, 518C>T, 556C>T, and 557G>A were less than 25% of wild-type allozymes. CPR activity levels were similar for all variants. In summary, sequencing of CYP2C9 and CYP2C19 in 1013 subjects identified low-frequency variants that had not previously been functionally characterized. In silico predictions were not always consistent with functional assay results. These observations emphasize the need for high-throughput methods for pharmacogene variant mutagenesis and functional characterization.

Rational and semi-rational engineering of cytochrome P450s for biotechnological applications

The cytochrome P450 enzymes are ubiquitous heme-thiolate proteins performing regioselective and stereoselective oxygenation reactions in cellular metabolism. Due to their broad substrate scope and catalytic versatility, P450 enzymes are also attractive candidates for many industrial and biopharmaceutical applications. For particular uses, enzyme properties of P450s can be further optimized through directed evolution, rational, and semi-rational engineering approaches, all of which introduce mutations within the P450 structures. In this review, we describe the recent applications of these P450 engineering approaches and highlight the key regions and residues that have been identified using such approaches. These “hotspots” lie within critical functional areas of the P450 structure, including the active site, the substrate access channel, and the redox partner interaction interface.

The Cytochrome P450 Slow Metabolizers CYP2C9*2 and CYP2C9*3 Directly Regulate Tumorigenesis via Reduced Epoxyeicosatrienoic Acid Production

Increased expression of cytochrome P450 CYP2C9, together with elevated levels of its products epoxyeicosatrienoic acids (EET), is associated with aggressiveness in cancer. Cytochrome P450 variants CYP2C9*2 and CYP2C9*3 encode proteins with reduced enzymatic activity, and individuals carrying these variants metabolize drugs more slowly than individuals with wild-type CYP2C9*1, potentially affecting their response to drugs and altering their risk of disease. Although genetic differences in CYP2C9-dependent oxidation of arachidonic acid (AA) have been reported, the roles of CYP2C9*2 and CYP2C9*3 in EET biosynthesis and their relevance to disease are unknown. Here, we report that CYP2C9*2 and CYP2C9*3 metabolize AA less efficiently than CYP2C9*1 and that they play a role in the progression of non-small cell lung cancer (NSCLC) via impaired EET biosynthesis. When injected into mice, NSCLC cells expressing CYP2C9*2 and CYP2C9*3 produced lower levels of EETs and developed fewer, smaller, and less vascularized tumors than cells expressing CYP2C9*1. Moreover, endothelial cells expressing these two variants proliferated and migrated less than cells expressing CYP2C*1. Purified CYP2C9*2 and CYP2C9*3 exhibited attenuated catalytic efficiency in producing EETs, primarily due to impaired reduction of these two variants by NADPH-P450 reductase. Loss-of-function SNPs within CYP2C9*2 and CYP2C9*3 were associated with improved survival in female cases of NSCLC. Thus, decreased EET biosynthesis represents a novel mechanism whereby CYPC29*2 and CYP2C9*3 exert a direct protective role in NSCLC development.Significance: These findings report single nucleotide polymorphisms in the human CYP2C9 genes, CYP2C9*2 and CYP2C9*3, exert a direct protective role in tumorigenesis by impairing EET biosynthesis. Cancer Res; 78(17); 4865-77. ©2018 AACR.

Molecular Determinants of Substrate Affinity and Enzyme Activity of a Cytochrome P450BM3 Variant

Cytochrome P450_BM3 catalyzes the hydroxylation and/or epoxidation of fatty acids, fatty amides, and alcohols. Protein engineering has produced P450_BM3 variants capable of accepting drug molecules normally metabolized by human P450 enzymes. The enhanced substrate promiscuity has been attributed to the greater flexibility of the lid of the substrate channel. However, it is not well understood how structurally different and highly polar drug molecules can stably bind in the active site nor how the activity and coupling efficiency of the enzyme may be affected by the lack of enzyme-substrate complementarity. To address these important aspects of non-native small molecule binding, this study investigated the binding of drug molecules with different size, charge, polar surface area, and human P450 affinity on the promiscuous R47L/F87V/L188Q/E267V/F81I pentuple mutant of P450_BM3. Binding free energy data and energy decomposition analysis showed that pentuple mutant P450_BM3 stably binds (i.e., negative ΔG_b°) a broad range of substrate and inhibitor types because dispersion interactions with active site residues overcome unfavorable repulsive and electrostatic effects. Molecular dynamics simulations revealed that 1) acidic substrates tend to disrupt the heme propionate A-K69 salt bridge, which may reduce heme oxidizing ability, and 2) the lack of complementarity leads to high substrate mobility and water density in the active site, which may lead to uncoupling. These factors must be considered in future developments of P450_BM3 as a biocatalyst in the large-scale production of drug metabolites.

Cytochrome P450: Polymorphisms and Roles in Cancer, Diabetes and Atherosclerosis

Cytochromes P450s (CYPs) constitute a superfamily of enzymes that catalyze the metabolism of drugs and other substances. Endogenous substrates of CYPs include eicosanoids, estradiol, arachidonic acids, cholesterol, vitamin D and neurotransmitters. Exogenous substrates of CYPs include the polycyclic aromatic hydrocarbons and about 80% of currently used drugs. Some isoforms can activate procarcinogens to ultimate carcinogens. Genetic polymorphisms of CYPs may affect the enzyme catalytic activity and have been reported among different populations to be associated with various diseases and adverse drug reactions. With regard of drug metabolism, phenotypes for CYP polymorphism range from ultrarapid to poor metabolizers. In this review, we discuss some of the most clinically important CYPs isoforms (CYP2D6, CYP2A6, CYP2C19, CYP2C9, CYP1B1 and CYP1A2) with respect to gene polymorphisms and drug metabolism. Moreover, we review the role of CYPs in renal, lung, breast and prostate cancers and also discuss their significance for atherosclerosis and type 2 diabetes mellitus.

Effects of polymorphic variation on the thermostability of heterogenous populations of CYP3A4 and CYP2C9 enzymes in solution

The effect of non-synonymous single nucleotide polymorphisms (SNPs) on cytochrome P450 (CYP450) drug metabolism is currently poorly understood due to the large number of polymorphisms, the diversity of potential substrates and the complexity of CYP450 function. Previously we carried out in silico studies to explore the effect of SNPs on CYP450 function, using in silico calculations to predict the effect of mutations on protein stability. Here we have determined the effect of eight CYP3A4 and seven CYP2C9 SNPs on the thermostability of proteins in solution to test these predictions. Thermostability assays revealed distinct CYP450 sub-populations with only 65–70% of wild-type CYP3A4 and CYP2C9 susceptible to rapid heat-induced P450 to P420 conversion. CYP3A4 mutations G56D, P218R, S222P, I223R, L373F and M445T and CYP2C9 mutations V76M, I359L and I359T were destabilising, increasing the proportion of protein sensitive to the rapid heat-induced P450 to P420 conversion and/or reducing the half-life of this conversion. CYP2C9 Q214L was the only stabilising mutation. These results corresponded well with the in silico protein stability calculations, confirming the value of these predictions and together suggest that the changes in thermostability result from destabilisation/stabilisation of the protein fold, changes in the haem-binding environment or effects on oligomer formation/conformation.

Can the inhibition of cytochrome P450 in aquatic invertebrates due to azole fungicides be estimated with in silico and in vitro models and extrapolated between species?

Azole fungicides, designed to halt fungal growth by specific inhibition of fungal cytochrome P450 (CYP51), inhibit cytochrome P450s involved in the metabolism of xenobiotics in several non-target organisms thus raising environmental concern. The present study investigates the degree by which inhibition strengths of azoles toward cytochrome P450 in rat liver, the insect Chironomus riparius larvae and the snail Lymnaea stagnalis can be extrapolated from estimated in silico affinities. Azoles' affinities toward human cytochrome P450 isoforms involved in xenobiotic metabolism (CYP3A4, CYP2C9 and CYP2D6) as well as fungal CYP51 were estimated with a ligand-protein docking model based on the ChemScore scoring function. Estimated affinities toward the selected enzymatic structures correlated strongly with measured inhibition strengths in rat liver (ChemScore vs. logIC₅₀ among cytochrome P450 isoforms: -0.662 < r < -0.891, n = 17 azoles), while weaker correlations were found for C. riparius larvae (-0.167 < r < -0.733, n = 9) and L. stagnalis (-0.084 < r < -0.648, n = 8). Inhibition strengths toward C. riparius and rat liver activities were found to be highly correlated to each other (r: 0.857) while no significant relationship was found between either of the species and L. stagnalis. The inhibition of cytochrome P450 due to azole fungicides could be estimated in vitro and to a lesser extent in silico for C. riparius but not for L. stagnalis, possibly due to different enzymatic susceptibility toward azole inhibition among the species.

Insights into molecular mechanisms of drug metabolism dysfunction of human CYP2C9*30

Cytochrome P450 2C9 (CYP2C9) metabolizes about 15% of clinically administrated drugs. The allelic variant CYP2C9*30 (A477T) is associated to diminished response to the antihypertensive effects of the prodrug losartan and affected metabolism of other drugs. Here, we investigated molecular mechanisms involved in the functional consequences of this amino-acid substitution. Molecular dynamics (MD) simulations performed for the active species of the enzyme (heme in the Compound I state), in the apo or substrate-bound state, and binding energy analyses gave insights into altered protein structure and dynamics involved in the defective drug metabolism of human CYP2C9.30. Our data revealed an increased rigidity of the key Substrate Recognition Sites SRS1 and SRS5 and shifting of the β turn 4 of SRS6 toward the helix F in CYP2C9.30. Channel and binding substrate dynamics analyses showed altered substrate channel access and active site accommodation. These conformational and dynamic changes are believed to be involved in the governing mechanism of the reduced catalytic activity. An ensemble of representative conformations of the WT and A477T mutant properly accommodating drug substrates were identified, those structures can be used for prediction of new CYP2C9 and CYP2C9.30 substrates and drug-drug interactions.

Cytochrome P450 2C9-natural antiarthritic interactions: Evaluation of inhibition magnitude and prediction from in vitro data

Many dietary supplements are promoted to patients with osteoarthritis (OA) including the three naturally derived compounds, glucosamine, chondroitin and diacerein. Despite their wide spread use, research on interaction of these antiarthritic compounds with human hepatic cytochrome P450 (CYP) enzymes is limited. This study aimed to examine the modulatory effects of these compounds on CYP2C9, a major CYP isoform, using in vitro biochemical assay and in silico models. Utilizing valsartan hydroxylase assay as probe, all forms of glucosamine and chondroitin exhibited IC₅₀ values beyond 1000 μM, indicating very weak potential in inhibiting CYP2C9. In silico docking postulated no interaction with CYP2C9 for chondroitin and weak bonding for glucosamine. On the other hand, diacerein exhibited mixed-type inhibition with IC₅₀ value of 32.23 μM and K_i value of 30.80 μM, indicating moderately weak inhibition. Diacerein's main metabolite, rhein, demonstrated the same mode of inhibition as diacerein but stronger potency, with IC₅₀ of 6.08 μM and K_i of 1.16 μM. The docking of both compounds acquired lower CDOCKER interaction energy values, with interactions dominated by hydrogen and hydrophobic bondings. The ranking with respect to inhibition potency for the investigated compounds was generally the same in both in vitro enzyme assay and in silico modeling with order of potency being diacerein/rhein > various glucosamine/chondroitin forms. In vitro-in vivo extrapolation of inhibition kinetics (using 1 + [I]/K_i ratio) demonstrated negligible potential of diacerein to cause interaction in vivo, whereas rhein was predicted to cause in vivo interaction, suggesting potential interaction risk with the CYP2C9 drug substrates.

Identification of karanjin isolated from the Indian beech tree as a potent CYP1 enzyme inhibitor with cellular efficacy via screening of a natural product repository

CYP1A1 is thought to mediate carcinogenesis in oral, lung and epithelial cancers. In order to identify a CYP1A1 inhibitor from an edible plant, 394 natural products in the IIIM's natural product repository were screened, at 10 μM concentration, using CYP1A1-Sacchrosomes™ (i.e. microsomal enzyme isolated from recombinant baker's yeast). Twenty-seven natural products were identified that inhibited 40-97% of CYP1A1's 7-ethoxyresorufin-O-deethylase activity. The IC50 values of the 'hits', belonging to different chemical scaffolds, were determined. Their selectivity was studied against a panel of 8 CYP-Sacchrosomes™. In order to assess cellular efficacy, the 'hits' were screened for their capability to inhibit CYP enzymes expressed within live recombinant human embryonic kidney (HEK293) cells from plasmids encoding specific CYP genes (1A2, 1B1, 2C9, 2C19, 2D6, 3A4). Isopimpinellin (IN-475; IC50, 20 nM) and karanjin (IN-195; IC50, 30 nM) showed the most potent inhibition of CYP1A1 in human cells. Isopimpinellin is found in celery, parsnip, fruits and in the rind and pulp of limes whereas different parts of the Indian beech tree, which contain karanjin, have been used in traditional medicine. Both isopimpinellin and karanjin negate the cellular toxicity of CYP1A1-mediated benzo[a]pyrene. Molecular docking and molecular dynamic simulations with CYP isoforms rationalize the observed trends in the potency and selectivity of isopimpinellin and karanjin.

Molecular dynamics investigations of regioselectivity of anionic/aromatic substrates by a family of enzymes: a case study of diclofenac binding in CYP2C isoforms

The CYP2C subfamily is of particular importance in the metabolism of drugs, food toxins, and procarcinogens. Like other P450 subfamilies, 2C enzymes share a high sequence identity, but significantly contribute in different ways to hepatic capacity to metabolize drugs. They often metabolize the same substrate to more than one product with different catalytic sites. Because it is challenging to characterize experimentally, much still remains unknown about the reason for why the substrate regioselectivity of these closely related subfamily members is different. Here, we have investigated the structural features of CYP2C8, CYP2C9, and CYP2C19 bound with their shared substrate diclofenac to elucidate the underlying molecular mechanism for the substrate regioselectivity of CYP2C subfamily enzymes. The obtained results demonstrate how a sequence divergence for the active site residues causes heterogeneous variations in the secondary structures and in major tunnel selections, and further affects the shape and chemical properties of the substrate-binding site. Structural analysis and free energy calculations showed that the most important determinants of regioselectivity among the CYP2C isoforms are the geometrical features of the active sites, as well as the hydrogen bonds and the hydrophobic interactions, mainly presenting as the various locations of Arg108 and substitutions of Phe205 for Ile205 in CYP2C8. The MM-GB/SA calculations combined with PMF results accord well with the experimental KM values, bridging the gap between the theory and the experimentally observed results of binding affinity differences. The present study provides important insights into the structure-function relationships of CYP2C subfamily enzymes, the knowledge of ligand binding characteristics and key residue contributions could guide future experimental and computational work on the synthesis of drugs with better pharmacokinetic properties so that CYP interactions could be avoided.