Concurrent cooperativity and substrate inhibition in the epoxidation of carbamazepine by cytochrome P450 3A4 active site mutants inspired by molecular dynamics simulations

CYP2J2-mediated metabolism of arachidonic acid in heart: A review of its kinetics, inhibition and role in heart rhythm control

Cytochrome P450 2 J2 (CYP2J2) is primarily expressed extrahepatically and is the predominant epoxygenase in human cardiac tissues. This highlights its key role in the metabolism of endogenous substrates. Significant scientific interest lies in cardiac CYP2J2 metabolism of arachidonic acid (AA), an omega-6 polyunsaturated fatty acid, to regioisomeric bioactive epoxyeicosatrienoic acid (EET) metabolites that show cardioprotective effects including regulation of cardiac electrophysiology. From an in vitro perspective, the accurate characterization of the kinetics of CYP2J2 metabolism of AA including its inhibition and inactivation by drugs could be useful in facilitating in vitro-in vivo extrapolations to predict drug-AA interactions in drug discovery and development. In this review, background information on the structure, regulation and expression of CYP2J2 in human heart is presented alongside AA and EETs as its endogenous substrate and metabolites. The in vitro and in vivo implications of the kinetics of this endogenous metabolic pathway as well as its perturbation via inhibition and inactivation by drugs are elaborated. Additionally, the role of CYP2J2-mediated metabolism of AA to EETs in cardiac electrophysiology will be expounded.

[Metabolic Activities Catalyzed by Human Cytochrome P450 (CYP) 2D6 and CYP3A Subfamily Members and Effect of Various Compounds, Including Endogenous Steroid Hormones, on These Activities]

My research focused on the effects of various drugs on (1) dopamine formation from p-tyramine catalyzed by polymorphic cytochrome P450 (CYP or P450) 2D6 variants and (2) endogenous steroid hormone hydroxylation catalyzed by CYP3A subfamily members (CYP3A4, CYP3A5, CYP3A7). The activation (cooperativity) of metabolic reactions catalyzed by P450s was especially emphasized. The effects of various psychotropic agents on dopamine formation from p-tyramine, catalyzed by wild-type CYP2D6.1 and CYP2D6 variants, including CYP2D6.2 (Arg296Cys;Ser486Thr), CYP2D6.10 (Pro34Ser;Ser486Thr), and CYP2D6.39 (Ser486Thr) were compared. Michaelis (Km) and inhibition (Ki) constants of the psychotropic agents in the presence of CYP2D6.10 were higher than those observed in the presence of other CYP2D6 variants. Fluvoxamine, fluoxetine, milnacipran, and haloperidol activated CYP2D6-catalyzed dopamine formation [decreasing the Km and/or increasing the maximal velocity (kcat)], and this activation was CYP2D6 variant-dependent. Regarding the CYP3A subfamily, the effects of various compounds including endogenous steroid hormones on the 6β-hydroxylation of steroid hormones, such as testosterone, progesterone, and cortisol, were determined; it was found that testosterone, dehydroepiandrosterone, and/or α-naphthoflavone activated 6β-hydroxylation of cortisol and/or progesterone, but the effects varied in the presence of different CYP3A subfamily members. Further studies are required to confirm the mechanisms and therapeutic relevance of these activation phenomena.

300-fold higher neuro- and immunotoxicity from low-redox transformation of carbamazepine

Current challenges in (eco)toxicology are in understanding the transformation of (reactive) substances, and how transformation affects toxic modes of action. Empirical assessment of transformation products of, practically an infinite number of substances, via experimentation, is impossible. Predicting transformation products for (benchmarking) compounds from conditions, facilitates risk analyses. This study applied calculus to predict transformation products of an important environmental and medicinal/toxicological marker, carbamazepine. As radicals are ubiquitous in humans and the environment, we looked into radical-mediated transformations of carbamazepine as a benchmark. We calculated proportions of their speciation states as function of redox conditions, which we took as pH and O2 concentration, describing transformation via covalent and ionic interactions. Formation of ring-contracted products with neuro-immunological activity is thermodynamically favored under anaerobic conditions and at low pH. Experimentally observed product distributions and toxicities reflect that pattern. Our predictive method may support toxicity predictions for other substances and conditions 'similar' to the current case study via interpolation. This paves the way for a more coherent, effective and easier risk assessment of transformation products.

Atypical kinetics of cytochrome P450 enzymes in pharmacology and toxicology

Atypical kinetics are observed in metabolic reactions catalyzed by cytochrome P450 enzymes (P450). Yet, this phenomenon is regarded as experimental artifacts in some instances despite increasing evidence challenging the assumptions of typical Michaelis-Menten kinetics. As P450 play a major role in the metabolism of a wide range of substrates including drugs and endogenous compounds, it becomes critical to consider the impact of atypical kinetics on the accuracy of estimated kinetic and inhibitory parameters which could affect extrapolation of pharmacological and toxicological implications. The first half of this book chapter will focus on atypical non-Michaelis-Menten kinetics (e.g. substrate inhibition, biphasic and sigmoidal kinetics) as well as proposed underlying mechanisms supported by recent insights in mechanistic enzymology. In particular, substrate inhibition kinetics in P450 as well as concurrent drug inhibition of P450 in the presence of substrate inhibition will be further discussed. Moreover, mounting evidence has revealed that despite the high degree of sequence homology between CYP3A isoforms (i.e. CYP3A4 and CYP3A5), they have the propensities to exhibit vastly different susceptibilities and potencies of mechanism-based inactivation (MBI) with a common drug inhibitor. These experimental observations pertaining to the presence of these atypical isoform- and probe substrate-specific complexities in CYP3A isoforms by several clinically-relevant drugs will therefore be expounded and elaborated upon in the second half of this book chapter.

Midazolam as a Probe for Heterotropic Drug-Drug Interactions Mediated by CYP3A4

Human cytochrome P450 CYP3A4 is involved in the processing of more than 35% of current pharmaceuticals and therefore is responsible for multiple drug-drug interactions (DDI). In order to develop a method for the detection and prediction of the possible involvement of new drug candidates in CYP3A4-mediated DDI, we evaluated the application of midazolam (MDZ) as a probe substrate. MDZ is hydroxylated by CYP3A4 in two positions: 1-hydroxy MDZ formed at lower substrate concentrations, and up to 35% of 4-hydroxy MDZ at high concentrations. The ratio of the formation rates of these two products (the site of metabolism ratio, SOM) was used as a measure of allosteric heterotropic interactions caused by effector molecules using CYP3A4 incorporated in lipid nanodiscs. The extent of the changes in the SOM in the presence of effectors is determined by chemical structure and is concentration-dependent. MD simulations of CYP3A4 in the lipid bilayer suggest that experimental results can be explained by the movement of the F-F' loop and concomitant changes in the shape and volume of the substrate-binding pocket. As a result of PGS binding at the allosteric site, several residues directly contacting MDZ move away from the substrate molecule, enabling the repositioning of the latter for minor product formation.

Pseudomonas aeruginosa cytochrome P450 CYP168A1 is a fatty acid hydroxylase that metabolizes arachidonic acid to the vasodilator 19-HETE

Pseudomonas aeruginosa is a Gram-negative opportunistic human pathogen that is highly prevalent in individuals with cystic fibrosis (CF). A major problem in treating CF patients infected with P. aeruginosa is the development of antibiotic resistance. Therefore, the identification of novel P. aeruginosa antibiotic drug targets is of the utmost urgency. The genome of P. aeruginosa contains four putative cytochrome P450 enzymes (CYPs) of unknown function that have never before been characterized. Analogous to some of the CYPs from Mycobacterium tuberculosis, these P. aeruginosa CYPs may be important for growth and colonization of CF patients’ lungs. In this study, we cloned, expressed, and characterized CYP168A1 from P. aeruginosa and identified it as a subterminal fatty acid hydroxylase. Spectral binding data and computational modeling of substrates and inhibitors suggest that CYP168A1 has a large, expansive active site and preferentially binds long chain fatty acids and large hydrophobic inhibitors. Furthermore, metabolic experiments confirm that the enzyme is capable of hydroxylating arachidonic acid, an important inflammatory signaling molecule present in abundance in the CF lung, to 19-hydroxyeicosatetraenoic acid (19-HETE; K_m = 41 μM, V_max = 220 pmol/min/nmol P450), a potent vasodilator, which may play a role in the pathogen’s ability to colonize the lung. Additionally, we found that the in vitro metabolism of arachidonic acid is subject to substrate inhibition and is also inhibited by the presence of the antifungal agent ketoconazole. This study identifies a new metabolic pathway in this important human pathogen that may be of utility in treating P. aeruginosa infections.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Atypical kinetics of cytochrome P450 2J2: Epoxidation of arachidonic acid and reversible inhibition by xenobiotic inhibitors

Extrahepatic CYP2J2 metabolism of arachidonic acid (AA) to bioactive regioisomeric epoxyeicosatrienoic acids (EETs) is implicated in both physiological and pathological conditions. Here, we aimed to characterize atypical substrate inhibition kinetics of this endogenous metabolic pathway and its reversible inhibition by xenobiotic inhibitors when AA is used as the physiologically-relevant substrate vis-à-vis conventional probe substrate astemizole (AST). As compared to typical Michaelis-Menten kinetics observed for AST, complete substrate inhibition was observed for CYP2J2 metabolism of AA to 14,15-EET whereby velocity of the reaction declined significantly at concentrations of AA above 20-30 µM with an estimated substrate inhibition constant (K_s) of 31 µM. In silico sequential docking of two AA substrates to orthosteric (OBS) and adjacent secondary binding sites (SBS) within a 3-dimensional homology model of CYP2J2 revealed favorable and comparable binding poses of glide-scores -3.1 and -3.8 respectively. Molecular dynamics (MD) simulations ascertained CYP2J2 conformational stability with dual AA substrate binding as time-dependent root mean squared deviation (RMSD) of protein Cα atoms and ligand heavy atoms stabilized to a plateau in all but one trajectory (n=6). The distance between heme-iron and ω6 (C14, C15) double bond of AA in OBS also increased from 7.5 ± 1.4 Å to 8.5 ± 1.8 Å when CYP2J2 was simulated with only AA in OBS versus the presence of AA in both OBS and SBS (p<0.001), supporting the observed in vitro substrate inhibition phenomenon. Poor correlation was observed between inhibitory constants (K_i) determined for a panel of nine competitive and mixed mode xenobiotic inhibitors against CYP2J2 metabolism of AA as compared to AST, whereby 4 out of 9 drugs had a greater than 5-fold difference between K_i values. Nonlinear Eadie-Hofstee plots illustrated that complete substrate inhibition of CYP2J2 by AA was not attenuated even at high concentrations of xenobiotic inhibitors which further corroborates that CYP2J2 may accommodate three or more ligands simultaneously. In light of the atypical kinetics, our results highlight the importance of using physiologically-relevant substrates in in vitro enzymatic inhibition assays for the characterization of xenobiotic-endobiotic interactions which is applicable to other complex endogenous metabolic pathways beyond CYP2J2 metabolism of AA to EETs. The accurate determination of K_i would further facilitate the association of xenobiotic-endobiotic interactions to observed therapeutic or toxic outcomes.

Combining Small-Molecule Bioconjugation and Hydrogen-Deuterium Exchange Mass Spectrometry (HDX-MS) to Expose Allostery: the Case of Human Cytochrome P450 3A4

We report a novel approach to study allostery which combines the use of carefully selected bioconjugates and hydrogen-deuterium exchange mass spectrometry (HDX-MS). This strategy avoids issues related to weak substrate binding and ligand relocalization. The utility of our method is demonstrated using human cytochrome P450 3A4 (CYP3A4), the most important drug-metabolizing enzyme. Allosteric activation and inhibition of CYP3A4 by pharmaceuticals is an important mechanism of drug interactions. We performed HDX-MS analysis on several CYP3A4-effector bioconjugates, some of which mimic the allosteric effect of positive effectors, while others show activity enhancement even though the label does not occupy the allosteric pocket (agonistic) or do not show activation while still blocking the allosteric site (antagonistic). This allowed us to better define the position of the allosteric site, the protein structural dynamics associated with allosteric activation, and the presence of coexisting conformers.

Midazolam as a Probe for Drug-Drug Interactions Mediated by CYP3A4: Homotropic Allosteric Mechanism of Site-Specific Hydroxylation

We developed an efficient and sensitive probe for drug-drug interactions mediated by human CYP3A4 by using midazolam (MDZ) as a probe substrate. Using global analysis of four parameters over several experimental data sets, we demonstrate that the first MDZ molecule (MDZ1) binds with high affinity at the productive site near the heme iron and gives only hydroxylation at the 1 position (1OH). The second midazolam molecule (MDZ2) binds at an allosteric site at the membrane surface and perturbs the position and mobility of MDZ1 such that the minor hydroxylation product at the 4 position (4OH) is formed in a 1:2 ratio (35%). No increase in catalytic rate is observed after the second MDZ binding. Hence, the site of the 1OH:4OH metabolism ratio is a sensitive probe for drugs, such as progesterone, that bind with high affinity to the allosteric site and serve as effectors. We observe similar changes in the MDZ 1OH:4OH ratio in the presence of progesterone (PGS), suggesting a direct communication between the active and allosteric sites. Mutations introduced into the F-F' loop indicate that residues F213 and D214 are directly involved in allosteric interactions leading to MDZ homotropic cooperativity, and these same residues, together with L211, are involved in heterotropic allosteric interactions in which PGS is the effector and MDZ the substrate. Molecular dynamics simulations provide a mechanistic picture of the origin of this cooperativity. These results show that the midazolam can be used as a sensitive probe for drug-drug interactions in human P450 CYP3A4.

Comparison of the Stimulatory and Inhibitory Effects of Steroid Hormones and α-Naphthoflavone on Steroid Hormone Hydroxylation Catalyzed by Human Cytochrome P450 3A Subfamilies

The inhibitory and stimulatory effects of steroid hormones and related compounds on the hydroxylation activity at the 6β-position of two steroid hormones, progesterone and testosterone, by CYP3A4, polymorphically expressed CYP3A5, and fetal CYP3A7 were compared to clarify the catalytic properties of the predominant forms of the human CYP3A subfamily. Hydroxylation activities of progesterone and testosterone by CYP3A4, CYP3A5, and CYP3A7 were estimated using HPLC. The Michaelis constants (K_m) for progesterone 6β-hydroxylation by CYP3A5 were markedly decreased in the presence of dehydroepiandrosterone (DHEA) and α-naphthoflavone (ANF), whereas progesterone and DHEA competitively inhibited testosterone 6β-hydroxylation mediated by CYP3A4, and progesterone competitively inhibited CYP3A5-mediated activity, which was weaker than that for CYP3A4. ANF noncompetitively inhibited testosterone 6β-hydroxylation mediated by both CYP3A4 and CYP3A5. Progesterone and testosterone 6β-hydroxylation mediated by CYP3A7 was inhibited or unaffected by DHEA, pregnenolone, and ANF. These results suggested that DHEA and ANF stimulated progesterone 6β-hydroxylation by CYP3A5 but not by CYP3A4 and CYP3A7; however, progesterone, DHEA, and ANF inhibited testosterone 6β-hydroxylation mediated by all CYP3A subfamily members. The inhibitory/stimulatory pattern of steroid-steroid interactions is different among CYP3A subfamily members and CYP3A5 is the most sensitive in terms of activation among the CYP3A subfamily members investigated.

Kinetic Modeling of Steady-State Situations in Cytochrome P450 Enzyme Reactions

In the course of investigations of the kinetics of individual reactions of cytochrome P450 (P450) enzymes, a number of points about the complexity of P450 enzyme kinetics have become apparent. Several of these are of particular relevance to work with P450 enzymes in the course of drug development and lead optimization, particularly with regard to estimating in vitro kinetic parameters and dealing with enzyme inhibitors. Modern simulation modeling has been applied to situations involving issues of preincubation time with moderate strength and strong inhibitors, inhibition by tightly bound ligands that have been identified in P450 enzymes, extensive substrate depletion, P450 reactions with a rate-limiting step after product formation, and the consumption of an inhibitor during a reaction by either a P450 enzyme being monitored or another one in a mixture. The results all follow from first principles, and simulations reveal the extent of their significance in various settings. The order of addition of substrate and inhibitors can change the apparent outcome (inhibition constant, K _i), and the effect of the order is more pronounced with a stronger inhibitor. Substrate depletion alters parameters (Michaelis constant, K _m) and can generate apparently sigmoidal plots. A rate-limiting step after product formation lowers the apparent K _m and distorts K _i Consumption of an inhibitor during a reaction affects K _i and differs depending on which enzyme is involved. The results are relevant with P450 enzymes and other enzymes as well. SIGNIFICANCE STATEMENT: Kinetic simulations have been used to address several potential problems in enzyme kinetic analysis. Although the simulations done here are general for enzyme reactions, several problems addressed here are particularly relevant to cytochrome P450 reactions encountered in drug development work.

Efavirenz Metabolism: Influence of Polymorphic CYP2B6 Variants and Stereochemistry

Efavirenz (more specifically the S-enantiomer) is a cornerstone antiretroviral therapy for treatment of HIV infection. The major primary metabolite is S-8-hydroxyefavirenz, which does not have antiretroviral activity but is neurotoxic. Cytochrome P450 2B6 (CYP2B6) is the major enzyme catalyzing S-8-hydroxyefavirenz formation. CYP2B6 genetics and drug interactions are major determinants of clinical efavirenz disposition and dose adjustment. In addition, as a prototypic CYP2B6 substrate, S-efavirenz and analogs can inform on the structure, activity, catalytic mechanisms, and stereoselectivity of CYP2B6. Metabolism of R-efavirenz by CYP2B6 remains unexplored. This investigation assessed S-efavirenz metabolism by clinically relevant CYP2B6 genetic variants. This investigation also evaluated R-efavirenz hydroxylation by wild-type CYP2B6.1 and CYP2B6 variants. S-Efavirenz 8-hydroxylation by wild-type CYP2B6.1 and variants exhibited positive cooperativity and apparent cooperative substrate inhibition. On the basis of Cl_max values, relative activities for S-efavirenz 8-hydroxylation were in the order CYP2B6.4 > CYP2B6.1 ≈ CYP2B6.5 ≈ CYP2B6.17 > CYP2B6.6 ≈ CYP2B6.7 ≈ CYP2B6.9 ≈ CYP2B6.19 ≈ CYP2B6.26; CYP2B6.16 and CYP2B6.18 showed minimal activity. Rates of R-efavirenz metabolism were approximately 1/10 those of S-efavirenz for wild-type CYP2B6.1 and variants. On the basis of Cl_max values, there was 14-fold enantioselectivity (S > R-efavirenz) for wild-type CYP2B6.1, and 5- to 22-fold differences for other CYP2B6 variants. These results show that both CYP2B6 516G > T (CYP2B6*6 and CYP2B6*9) and 983T > C (CYP2B6*16 and CYP2B6*18) polymorphisms cause canonical diminishment or loss-of-function variants for S-efavirenz 8-hydroxylation, provide a mechanistic basis for known clinical pharmacogenetic differences in efavirenz disposition, and may predict additional clinically important variant alleles. Efavirenz is the most stereoselective CYP2B6 drug substrate yet identified and may be a useful probe for the CYP2B6 active site and catalytic mechanisms. SIGNIFICANCE STATEMENT: Clinical disposition of the antiretroviral S-efavirenz is affected by CYP2B6 polymorphisms. Expressed CYP2B6 with 516G>T (CYP2B6*6 and CYP2B6*9), and 983T>C (CYP2B6*16 and CYP2B6*18) polymorphisms had a diminishment or loss of function for efavirenz 8-hydroxylation. This provides a mechanistic basis for efavirenz clinical pharmacogenetics and may predict additional clinically important variant alleles. Efavirenz metabolism showed both cooperativity and cooperative substrate inhibition. With greater than 10-fold enantioselectivity (S- vs. R- metabolism), efavirenz is the most stereoselective CYP2B6 drug substrate yet identified. These findings may provide mechanistic insights.

Theoretical explanation for the pharmaceutical incompatibility through the cooperativity effect of the drug-drug intermolecular interactions in the phenobarbital∙∙∙paracetamol∙∙∙H2O complex

In order to reveal the essence of the pharmaceutical incompatibility, the cooperativity effects of the drug-drug intermolecular π∙∙∙π and H∙∙∙O H-bonding interactions involving hydration were evaluated in the phenobarbital∙∙∙paracetamol∙∙∙H₂O complex at the M06-2X/6-311++G** and MP2/6-311++G** levels. The thermodynamic cooperativity effects were also investigated by the statistical thermodynamic method. The results show that the π∙∙∙π stacking ternary complexes with the moderate anti-cooperativity effects are dominant in controling the aggregation process of phenobarbital, paracetamol, and H₂O, as is confirmed by the atoms-in-molecules (AIM) and reduced density gradient (RDG) analyses. Therefore, it can be inferred that the anti-cooperativity effect plays an important role in forming the pharmaceutical incompatibility, and thus a deduction on the formation process of the pharmaceutical incompatibility between phenobarbital and paracetamol, with the hydration effect, is given. Several valuable models that relate the features of molecular surface electrostatic potentials or their statistical parameters, such as the surface areas, average values ([Formula: see text]), variances ([Formula: see text], [Formula: see text] and [Formula: see text]), and product of [Formula: see text] and electrostatic balance parameter (ν) ([Formula: see text]ν), to the values of the cooperativity effects were predicted. The formation of the pharmaceutical incompatibility is a thermodynamic cooperativity process driven by the enthalpy change. Graphical abstract Anti-cooperativity effect plays an important role in forming the pharmaceutical incompatibility.

Allosteric Interactions in Human Cytochrome P450 CYP3A4: The Role of Phenylalanine 213

The role of Phe213 in the allosteric mechanism of human cytochrome P450 CYP3A4 was studied using a combination of progesterone (PGS) and carbamazepine (CBZ) as probe substrates. We expressed, purified, and incorporated into POPC Nanodiscs three mutants, F213A, F213S, and F213Y, and compared them with wild-type (WT) CYP3A4 by monitoring spectral titration, the rate of NADPH oxidation, and steady-state product turnover rates with pure substrates and substrate mixtures. All mutants demonstrated higher activity with CBZ, lower activity with PGS, and a reduced level of activation of CBZ epoxidation by PGS, which was most pronounced in the F213A mutant. Using all-atom molecular dynamics simulations, we compared the dynamics of WT CYP3A4 and the F213A mutant incorporated into the lipid bilayer and the effect of the presence of the PGS molecule at the allosteric peripheral site and evaluated the critical role of Phe213 in mediating the heterotropic allosteric interactions in CYP3A4.

Drug-Drug Interactions between Atorvastatin and Dronedarone Mediated by Monomeric CYP3A4

Heterotropic interactions between atorvastatin (ARVS) and dronedarone (DND) have been deciphered using global analysis of the results of binding and turnover experiments for pure drugs and their mixtures. The in vivo presence of atorvastatin lactone (ARVL) was explicitly taken into account by using pure ARVL in analogous experiments. Both ARVL and ARVS inhibit DND binding and metabolism, while a significantly higher affinity of CYP3A4 for ARVL makes the latter the main modulator of activity (effector) in this system. Molecular dynamics simulations reveal significantly different modes of interactions of DND and ARVL with the substrate binding pocket and with a peripheral allosteric site. Interactions of both substrates with residues F213 and F219 at the allosteric site play a critical role in the communication of conformational changes induced by effector binding to productive binding of the substrate at the catalytic site.

Digging Deeper into CYP3A Testosterone Metabolism: Kinetic, Regioselectivity, and Stereoselectivity Differences between CYP3A4/5 and CYP3A7

The metabolism of testosterone to 6β-hydroxytestosterone (6β-OH-T) is a commonly used assay to evaluate human CYP3A enzyme activities. However, previous reports have indicated that CYP3A7 also produces 2α-hydroxytestosterone (2α-OH-T) and that a 2α-OH-T/6β-OH-T ratio may be a unique endogenous biomarker of the activity of the enzyme. Until now, the full metabolite and kinetic profile for testosterone hydroxylation by CYP3A7 has not been fully examined. To this end, we performed a complete kinetic analysis of the 6β-OH-T, 2α-OH-T, and 2β-hydroxytestosterone metabolites for recombinant Supersome CYP3A4, CYP3A5, and CYP3A7 enzymes and monitored metabolism in fetal and adult human liver microsomes for comparison. In general, a decrease in the velocity of the reaction was observed between CYP3A4 and the two other enzymes, with CYP3A7 showing the lowest metabolic capacity. Interestingly, we found that the 2α-OH-T/6β-OH-T ratio varied with substrate concentration when testosterone was incubated with CYP3A7, suggesting that this ratio would likely not function well as a biomarker for CYP3A7 activity. In silico docking studies revealed at least two different binding modes for testosterone between CYP3A4 and CYP3A7. In CYP3A4, the most energetically favorable docking mode places testosterone in a position with the methyl groups directed toward the heme iron, which is more favorable for oxidation at C6β, whereas for CYP3A7 the testosterone methyl groups are positioned away from the heme, which is more favorable for an oxidation event at C2α In conclusion, our data indicate an alternative binding mode for testosterone in CYP3A7 that favors the 2α-hydroxylation, suggesting significant structural differences in its active site compared with CYP3A4/5.

Use of bioconjugation with cytochrome P450 enzymes

Bioconjugation, defined as chemical modification of biomolecules, is widely employed in biological and biophysical studies. It can expand functional diversity and enable applications ranging from biocatalysis, biosensing and even therapy. This review summarizes how chemical modifications of cytochrome P450 enzymes (P450s or CYPs) have contributed to improving our understanding of these enzymes. Genetic modifications of P450s have also proven very useful but are not covered in this review. Bioconjugation has served to gain structural information and investigate the mechanism of P450s via photoaffinity labeling, mechanism-based inhibition (MBI) and fluorescence studies. P450 surface acetylation and protein cross-linking have contributed to the investigation of protein complexes formation involving P450 and its redox partner or other P450 enzymes. Finally, covalent immobilization on polymer surfaces or electrodes has benefited the areas of biocatalysis and biosensor design. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Solution Conformations and Dynamics of Substrate-Bound Cytochrome P450 MycG

MycG is a P450 monooxygenase that catalyzes the sequential hydroxylation and epoxidation of mycinamicin IV (M-IV), the last two steps in the biosynthesis of mycinamicin II, a macrolide antibiotic isolated from Micromonospora griseorubida. The crystal structure of MycG with M-IV bound was previously determined but showed the bound substrate in an orientation that did not rationalize the observed regiochemistry of M-IV hydroxylation. Nuclear magnetic resonance paramagnetic relaxation enhancements provided evidence of an orientation of M-IV in the MycG active site more compatible with the observed chemistry, but substrate-induced changes in the enzyme structure were not characterized. We now describe the use of amide ¹H-¹⁵N residual dipolar couplings as experimental restraints in solvated "soft annealing" molecular dynamics simulations to generate solution structural ensembles of M-IV-bound MycG. Chemical shift perturbations, hydrogen-deuterium exchange, and ¹⁵N relaxation behavior provide insight into the dynamic and electronic perturbations in the MycG structure in response to M-IV binding. The solution and crystallographic structures are compared, and the possibility that the crystallographic orientation of bound M-IV represents an inhibitory mode is discussed.

Allosteric Activation of Cytochrome P450 3A4 via Progesterone Bioconjugation

Human cytochrome P450 3A4 (CYP3A4) is responsible for the metabolism of the majority of drugs. As such, it is implicated in many adverse drug-drug and food-drug interactions, and is of significant interest to the pharmaceutical industry. This enzyme is known to simultaneously bind multiple ligands and display atypical enzyme kinetics, suggestive of allostery and cooperativity. As well, evidence of a postulated peripheral allosteric binding site has provoked debate around its significance and location. We report the use of bioconjugation to study the significance of substrate binding at the proposed allosteric site and its effect on CYP3A4 activity. CYP3A4 mutants were created and covalently modified with various small molecules including progesterone. The labeled mutants displayed enhanced kinetic stability and improved activity in testosterone and 7-benzyloxy-(4-trifluoromethyl)coumarin oxidation assays. Our work applies a new strategy to study cytochrome P450 allostery and supports the hypothesis that substrate binding at the postulated allosteric site of CYP3A4 may induce functional cooperativity.

In Silico Identification of Bioremediation Potential: Carbamazepine and Other Recalcitrant Personal Care Products

Emerging contaminants are principally personal care products not readily removed by conventional wastewater treatment and, with an increasing reliance on water recycling, become disseminated in drinking water supplies. Carbamazepine, a widely used neuroactive pharmaceutical, increasingly escapes wastewater treatment and is found in potable water. In this study, a mechanism is proposed by which carbamazepine resists biodegradation, and a previously unknown microbial biodegradation was predicted computationally. The prediction identified biphenyl dioxygenase from Paraburkholderia xenovorans LB400 as the best candidate enzyme for metabolizing carbamazepine. The rate of degradation described here is 40 times greater than the best reported rates. The metabolites cis-10,11-dihydroxy-10,11-dihydrocarbamazepine and cis-2,3-dihydroxy-2,3-dihydrocarbamazepine were demonstrated with the native organism and a recombinant host. The metabolites are considered nonharmful and mitigate the generation of carcinogenic acridine products known to form when advanced oxidation methods are used in water treatment. Other recalcitrant personal care products were subjected to prediction by the Pathway Prediction System and tested experimentally with P. xenovorans LB400. It was shown to biodegrade structurally diverse compounds. Predictions indicated hydrolase or oxygenase enzymes catalyzed the initial reactions. This study highlights the potential for using the growing body of enzyme-structural and genomic information with computational methods to rapidly identify enzymes and microorganisms that biodegrade emerging contaminants.

Spectroscopic studies and molecular docking on the interaction of organotin antitumor compound bis[2,4-difluoro-N-(hydroxy-⟨κ⟩O)benzamidato-⟨κ⟩O]diphenyltin(IV) with human cytochrome P450 3A4 protease

A novel organotin DFDPT was synthesized and characterized by elemental analysis, IR, (1)H, (13)C, (119)Sn, NMR techniques,etc. In order to investigate profoundly the relationship between DFDPT with human CYP3A4 proteaset and anticancer molecular mechanism of DFDPT, the intercalative mode of binding of DFDPT with CYP3A4 under physiological conditions were comprehensively evaluated using steady state, synchronous, three-dimensional fluorescence spectroscopy,circular dichroism and molecular docking. Fluorescence emission data showed that CYP3A4 fluorescence affected by DFDPT was a static quenching procedure, which implied that DFDPT-CYP3A4 complex had been formed. Apparent binding constants Kb of CYP3A4 with compound at 298 and 310K were 2.51×10(7) and 3.09×10(5), respectively. The binding sites number n was 1.64 and 1.22, respectively. The thermodynamic parameters ΔH and ΔS of the DFDPT-CYP3A4 complex were negative, which indicated that their interaction was driven mainly by hydrogen bonding and van der Waals force. The binding of DFDPT-CYP3A4 was spontaneous process in which ΔG was negative. The synchronous results showed DFDPT induced conformational changes of CYP3A4 protein. Three-dimensional fluorescence and circular dichroism spectra results also revealed conformation of CYP3A4 protein had been possible changed in the presence of DFDPT. Molecular docking was used to study the interaction orientation between DFDPT and CYP3A4 protease. The results indicated that DFDPT interacted with a panel of amino acids in the active sites of CYP3A4 protein mainly through formation of hydrogen bond. Furthermore, the predicted binding mode of DFDPT into CYP3A4 appeared to adopt an orientation with interactions among Arg105, Ser119 and Thr309.

Mechanism of drug-drug interactions mediated by human cytochrome P450 CYP3A4 monomer

Using Nanodiscs, we quantitate the heterotropic interaction between two different drugs mediated by monomeric CYP3A4 incorporated into a nativelike membrane environment. The mechanism of this interaction is deciphered by global analysis of multiple-turnover experiments performed under identical conditions using the pure substrates progesterone (PGS) and carbamazepine (CBZ) and their mixtures. Activation of CBZ epoxidation and simultaneous inhibition of PGS hydroxylation are measured and quantitated through differences in their respective affinities for both a remote allosteric site and the productive catalytic site near the heme iron. Preferred binding of PGS at the allosteric site and a stronger preference for CBZ binding at the productive site give rise to a nontrivial drug-drug interaction. Molecular dynamics simulations indicate functionally important conformational changes caused by PGS binding at the allosteric site and by two CBZ molecules positioned inside the substrate binding pocket. Structural changes involving Phe-213, Phe-219, and Phe-241 are thought to be responsible for the observed synergetic effects and positive allosteric interactions between these two substrates. Such a mechanism is likely of general relevance to the mutual heterotropic effects caused by biologically active compounds that exhibit different patterns of interaction with the distinct allosteric and productive sites of CYP3A4, as well as other xenobiotic metabolizing cytochromes P450 that are also involved in drug-drug interactions. Importantly, this work demonstrates that a monomeric CYP3A4 can display the full spectrum of activation and cooperative effects that are observed in hepatic membranes.