NON-MICHAELIS-MENTEN KINETICS IN CYTOCHROME P450-CATALYZED REACTIONS

Mechanistic studies on substrate inhibition and substrate activation of ribonuclease A: experimental and in silico investigation

Ribonuclease A (RNase A) is an endonuclease that plays a significant role in antimicrobial activity by the cleavage and hydrolysis of ssRNA. In this study, the interactions between RNase A and cCMP have been investigated, through molecular docking and a 200 ns molecular dynamics simulation. The enzyme kinetic properties were analyzed using saturation curve, Eadie-Hofstee, and Hill plots. The docking results indicate that the cCMP-RNase A complexes are stabilized through hydrophobic interaction, hydrogen bonding, and π-π stacking interaction. The most binding site is observed in the catalytic cleft, especially at residue His12 and His119. Enzyme-ligand docking study indicates that cCMP binds to residues located in the internal cavity of the catalytic site and shows three phases of binding modes. The analysis of MD simulations shows that RMSD, Rg, and RMSF reach equilibrium. The ΔGbinding of the cCMP-RNase A was negative (-619.673 kJ/mol), The comparison between the results pre and post MD simulation showed that ΔGbinding after MD simulation causes to more stable the structure than before simulation. Experimental methods such as saturation, Eadie-Hofstee, and Hill plots confirm theoretical data and show three phases of binding modes as well. Two significant events are demonstrated in the interaction between RNase A and cCMP: substrate activation and substrate inhibition are observed in concentrations below and above 0.8 mM, respectively, for cCMP. Choosing the appropriate concentration of cCMP is very important in investigation of ribonuclease A's catalytic behavour, especially for exploration of the effects of some drugs on biological behaviours related to ribonuclease A.Communicated by Ramaswamy H. Sarma.

Interspecies comparison of metabolism of two novel prototype PFAS

As a result of proposed global restrictions and regulations on current-use per-and polyfluoroalkyl substances (PFAS), research on possible alternatives is highly required. In this study, phase I in vitro metabolism of two novel prototype PFAS in human and rat was investigated. These prototype chemicals are intended to be safer-by-design and expected to mineralize completely, and thus be less persistent in the environment compared to the PFAS available on the market. Following incubation with rat liver S9 (RL-S9) fractions, two main metabolites per initial substance were produced, namely an alcohol and a short-chain carboxylic acid. While with human liver S9 (HL-S9) fractions, only the short-chain carboxylic acid was detected. Beyond these major metabolites, two and five additional metabolites were identified at very low levels by non-targeted screening for the ether- and thioether-linked prototype chemicals, respectively. Overall, complete mineralization during the in vitro hepatic metabolism of these novel PFAS by HL-S9 and RL-S9 fractions was not observed. The reaction kinetics of the surfactants was determined by using the metabolite formation, rather than the substrate depletion approach. With rat liver enzymes, the formation rates of primary metabolite alcohols were at least two orders of magnitude higher than those of secondary metabolite carboxylic acids. When incubating with human liver enzymes, the formation rates of single metabolite carboxylic acids, were similar or smaller than those experienced in rat. It also indicates that the overall metabolic rate and clearance of surfactants are significantly higher in rat liver than in human liver. The maximum formation rate of the thioether congener exceeded 10-fold that of the ether in humans but were similar in rats. Overall, the results suggest that metabolism of the prototype chemicals followed a similar trend to those reported in studies of fluorotelomer alcohols.

Inhibition Effect of Okanin Toward Human Cytochrome P450 3A4 and 2D6 with Multi-spectroscopic Studies and Molecular Docking

Okanin, a major flavonoid of a popular herb tea, Coreopsis tinctoria Nutt., showed strong inhibition on CYP3A4 and CYP2D6. The strong interaction between okanin and CYPs were determined by enzyme kinetics, multispectral technique and molecular docking. The inhibition type of two enzymes, CYP3A4 and CYP2D6, by okanin are mixed and non-competitive inhibition type, respectively. The IC50 values and the binding constant of okanin to CYP3A4 can be deduced that the interaction was stronger than that of CYP2D6. The Conformations of CYP3A4 and CYP2D6 were changed by okanin. The evidence from fluorescence measurement along with molecular docking verified that these two CYPs were bound with okanin by hydrogen bonds and hydrophobic forces. Our investigation suggested that okanin may lead to interactions between herb and drug by inhibiting CYP3A4 and CYP2D6 activities, thus its consumption should be taken with caution.

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.

Nanodisc-embedded cytochrome P450 P3A4 binds diverse ligands by distributing conformational dynamics to its flexible elements

Cytochrome P450 3A4 (CYP3A4) metabolizes a wide range of drugs and toxins. Interactions of CYP3A4 with ligands are difficult to predict due to promiscuity and conformational flexibility. To better understand CYP3A4 conformational responses to ligands we use hydrogen deuterium exchange mass spectrometry (HDX-MS) to investigate the effect of ligands on nanodisc-embedded CYP3A4. For a subset of CYP3A4-ligand complexes, differences in the low-frequency modes derived by principal component analyses of molecular dynamics trajectories mirrored the HDX-MS results. The effects of ligands are distributed to flexible elements of CYP3A4 between stretches of secondary structure. The largest effects occur in the F- and G-helices, where most ligands increase the flexibility of the F-helix and connecting loops and decrease the flexibility of the C-term of the G-helix. Most ligands affect the E-F-G, CD and HI regions of the protein. Ligand-dependent differences are observed in the A"-A' loop, BC region, E-helix, K-β1 region, proximal loop, and C-term loop. Correlated HDX responses were observed in the CD region and the C-term of the G-helix that were most pronounced for Type II ligands. Collectively, the HDX and molecular dynamics results suggest that CYP3A4 accommodates diverse binding partners by propagating local backbone fluctuations from the binding site onto the flexible regions of the enzyme via long-range interactions that are differentially modulated by ligands. In contrast to the paradigm wherein ligands decrease protein dynamics at their binding site, a wide range of ligands modestly increase CYP3A4 dynamics throughout the protein including effects remote from the active site.

Phytate-Induced Dose-Response Auto-Activation of Enzyme in Commercial Recombinant Phytase From Escherichia coli

Background

Microbial phytase is one of the most widely used enzymes in food industries like cattle, poultry, and aquaculture food. Therefore, understanding the kinetic properties of the enzyme is very important to evaluate and predict its behavior in the digestive system of livestock. Working on phytase is one of the most challenging experiments because of some problems, including free inorganic phosphate (FIP) impurity in phytate (substrate) and interference reaction of the reagent with both phosphates (product and phytate impurity).

Objective

In the present study, FIP impurity of phytate was removed, and then it was shown that the substrate (phytate) has a dual role in enzyme kinetics: substrate and activator.

Material and Methods

phytate impurity was decreased by two-step recrystallization prior to the enzyme assay. The impurity removal was estimated by the ISO30024:2009 method and confirmed by Fourier-transform infrared (FTIR) spectroscopy. The kinetic behavior of phytase activity was evaluated using the purified phytate as substrate by non-Michaelis-Menten analysis, including Eadie-Hofstee, Clearance, and Hill plots. The possibility of an allosteric site on phytase was assessed by molecular docking.

Results

The results showed a 97.2% decrease in FIP due to recrystallization. The phytase saturation curve had a sigmoidal appearance, and Lineweaver-Burk plot with a negative y-intercept indicated the positive homotropic effect of the substrate on the enzyme activity. A right-side concavity of Eadie-Hofstee plot confirmed it. Hill coefficient was calculated to be 2.26. Molecular docking also showed that Escherichia coli phytase molecule has another binding site for phytate very close to the active site, called "allosteric site".

Conclusions

The observations strongly propose the existence of an intrinsic molecular mechanism in Escherichia coli phytase molecules to be promoted for more activity by its substrate, phytate (positive homotropic allosteric effect). In silico analysis showed that phytate binding to the allosteric site caused new substrate-mediated inter-domain interactions, which seems to lead to a more active conformation of phytase. Our results provide a strong basis for animal feed development strategies, especially in the case of poultry food and supplements, regarding a short food passage time in their gastrointestinal tract and variable concentration of phytate along with it. Additionally, the results strengthen our understanding of phytase auto-activation as well as allosteric regulation of monomeric proteins in general.

Crystal Structure of CYP3A4 Complexed with Fluorol Identifies the Substrate Access Channel as a High-Affinity Ligand Binding Site

Cytochrome P450 3A4 (CYP3A4) is a major human drug-metabolizing enzyme, notoriously known for its extreme substrate promiscuity, allosteric behavior, and implications in drug–drug interactions. Despite extensive investigations, the mechanism of ligand binding to CYP3A4 is not fully understood. We determined the crystal structure of CYP3A4 complexed with fluorol, a small fluorescent dye that can undergo hydroxylation. In the structure, fluorol associates to the substrate channel, well suited for the binding of planar polyaromatic molecules bearing polar groups, through which stabilizing H-bonds with the polar channel residues, such as Thr224 and Arg372, can be established. Mutagenesis, spectral, kinetic, and functional data confirmed the involvement but not strict requirement of Thr224 for the association of fluorol. Collectively, our data identify the substrate channel as a high-affinity ligand binding site and support the notion that hydrophobic ligands first dock to the nearby peripheral surface, before migrating to the channel and, subsequently, into the active site.

Novel Unspecific Peroxygenase from Truncatella angustata Catalyzes the Synthesis of Bioactive Lipid Mediators

Lipid mediators, such as epoxidized or hydroxylated eicosanoids (EETs, HETEs) of arachidonic acid (AA), are important signaling molecules and play diverse roles at different physiological and pathophysiological levels. The EETs and HETEs formed by the cytochrome P450 enzymes are still not fully explored, but show interesting anti-inflammatory properties, which make them attractive as potential therapeutic target or even as therapeutic agents. Conventional methods of chemical synthesis require several steps and complex separation techniques and lead only to low yields. Using the newly discovered unspecific peroxygenase TanUPO from the ascomycetous fungus Truncatella angustata, 90% regioselective conversion of AA to 14,15-EET could be achieved. Selective conversion of AA to 18-HETE, 19-HETE as well as to 11,12-EET and 14,15-EET was also demonstrated with known peroxygenases, i.e., AaeUPO, CraUPO, MroUPO, MweUPO and CglUPO. The metabolites were confirmed by HPLC-ELSD, MS1 and MS2 spectrometry as well as by comparing their analytical data with authentic standards. Protein structure simulations of TanUPO provided insights into its substrate access channel and give an explanation for the selective oxyfunctionalization of AA. The present study expands the scope of UPOs as they can now be used for selective syntheses of AA metabolites that serve as reference material for diagnostics, for structure-function elucidation as well as for therapeutic and pharmacological purposes.

Structural Modeling of Drosophila melanogaster Gut Cytochrome P450s and Docking Comparison of Fruit Fly Gut and Human Cytochrome P450s

Drosophila melanogaster is a prominent model organism in developmental biology research and in studies related to pathophysiological conditions like cancer and Alzheimer's disease. The fruit fly gut contains several cytochrome P450s (CYP450s) which have central roles in Drosophila development and in the normal physiology of the gut. Since the crystal structures of these proteins have not been deciphered yet, we modeled the structure of 29 different D. melanogaster gut CYP450s using Prime (Schrödinger). The sequences of chosen D. melanogaster gut CYP450s were compared with that of their human counterparts. The common gut (and liver) microsomal CYP450s in humans were chosen for structural comparison to find the homology and identity % of D. melanogaster CYPs with that of their human counterparts. The modeled structures were validated using PROCHECK and the best fit models were used for docking several known human pharmacological agents/drugs to the modeled D. melanogaster gut CYP450s. Based on the binding affinities (ΔG values) of the selected drug molecules with the modeled fly gut CYPs, the plausible differences in metabolism of the prominent drugs in humans and fly were projected. The gut is involved in absorption of oral drugs/pharmacological agents and hence, upregulation of intestinal CYP450 and their reactions with endobiotics and xenobiotics is envisaged. The insights gleaned from this work can validate D. melanogaster as a model organism for studying intestinal drug metabolism, particularly in the context of a) toxicology of pharmacological agents to the gut cells and b) how gut P450 metabolites/products can influence gut homeostasis. This work can help establish a platform for further in vitro investigations on how intestinal CYP450 metabolism can influence gut health. The data from this work can be used for further in silico studies and this work can serve as a platform for future in vitro investigations on intestinal CYP450-mediated metabolism of endo- and xeno-biotics in D. melanogaster.

On the reproducibility of enzyme reactions and kinetic modelling

Enzyme reactions are highly dependent on reaction conditions. To ensure reproducibility of enzyme reaction parameters, experiments need to be carefully designed and kinetic modeling meticulously executed. Furthermore, to enable quality control of enzyme reaction parameters, the experimental conditions, the modeling process as well as the raw data need to be reported comprehensively. By taking these steps, enzyme reaction parameters can be open and FAIR (findable, accessible, interoperable, re-usable) as well as repeatable, replicable and reproducible. This review discusses these requirements and provides a practical guide to designing initial rate experiments for the determination of enzyme reaction parameters and gives an open, FAIR and re-editable example of the kinetic modeling of an enzyme reaction. Both the guide and example are scripted with Python in Jupyter Notebooks and are publicly available (https://fairdomhub.org/investigations/483/snapshots/1). Finally, the prerequisites of automated data analysis and machine learning algorithms are briefly discussed to provide further motivation for the comprehensive, open and FAIR reporting of enzyme reaction parameters.

Complex Cytochrome P450 Kinetics Due to Multisubstrate Binding and Sequential Metabolism. Part 1. Theoretical Considerations

Complexities in P450-mediated metabolism kinetics include multisubstrate binding, multiple-product formation, and sequential metabolism. Saturation curves and intrinsic clearances were simulated for single-substrate and multisubstrate models using derived velocity equations and numerical solutions of ordinary differential equations (ODEs). Multisubstrate models focused on sigmoidal kinetics because of their dramatic impact on clearance predictions. These models were combined with multiple-product formation and sequential metabolism, and simulations were performed with random error. Use of single-substrate models to characterize multisubstrate data can result in inaccurate kinetic parameters and poor clearance predictions. Comparing results for use of standard velocity equations with ODEs clearly shows that ODEs are more versatile and provide better parameter estimates. It would be difficult to derive concentration-velocity relationships for complex models, but these relationships can be easily modeled using numerical methods and ODEs. SIGNIFICANCE STATEMENT: The impact of multisubstrate binding, multiple-product formation, and sequential metabolism on the P450 kinetics was investigated. Numerical methods are capable of characterizing complicated P450 kinetics.

Characterization of fluorescent probe substrates to develop an efficient high-throughput assay for neonatal hepatic CYP3A7 inhibition screening

CYP3A7 is a member of the cytochrome P450 (CYP) 3A enzyme sub-family that is expressed in the fetus and neonate. In addition to its role metabolizing retinoic acid and the endogenous steroid dehydroepiandrosterone sulfate (DHEA-S), it also has a critical function in drug metabolism and disposition during the first few weeks of life. Despite this, it is generally ignored in the preclinical testing of new drug candidates. This increases the risk for drug-drug interactions (DDI) and toxicities occurring in the neonate. Therefore, screening drug candidates for CYP3A7 inhibition is essential to identify chemical entities with potential toxicity risks for neonates. Currently, there is no efficient high-throughput screening (HTS) assay to assess CYP3A7 inhibition. Here, we report our testing of various fluorescent probes to assess CYP3A7 activity in a high-throughput manner. We determined that the fluorescent compound dibenzylfluorescein (DBF) is superior to other compounds in meeting the criteria considered for an efficient HTS assay. Furthermore, a preliminary screen of an HIV/HCV antiviral drug mini-library demonstrated the utility of DBF in a HTS assay system. We anticipate that this tool will be of great benefit in screening drugs that may be used in the neonatal population in the future.

Enzyme Kinetics of Oxidative Metabolism-Cytochromes P450

The cytochrome P450 enzymes (CYPs) are the most important enzymes in the oxidative metabolism of hydrophobic drugs and other foreign compounds (xenobiotics). The versatility of these enzymes results in some unusual kinetic properties, stemming from the simultaneous interaction of multiple substrates with the CYP active site. Often, the CYPs display kinetics that deviate from standard hyperbolic saturation or inhibition kinetics. Non-Michaelis-Menten or "atypical" saturation kinetics include sigmoidal, biphasic, and substrate inhibition kinetics (see Chapter 2 ). Interactions between substrates include competitive inhibition, noncompetitive inhibition, mixed inhibition, partial inhibition, activation, and activation followed by inhibition (see Chapters 4 and 6 ). Models and equations that can result in these kinetic profiles will be presented and discussed.

Numerical Methods for Modeling Enzyme Kinetics

Differential equations are used to describe time-dependent changes in enzyme kinetics and pharmacokinetics. Analytical and numerical methods can be used to solve differential equations. This chapter describes the use of numerical methods in solving differential equations and its applications in characterizing the complexities observed in enzyme kinetics. A discussion is included on the use of numerical methods to overcome limitations of explicit equations in the analysis of metabolism kinetics, reversible inhibition kinetics, and inactivation kinetics. The chapter describes the advantages of using numerical methods when Michaelis-Menten assumptions do not hold.

Cooperative Kinetics of the Glucan Phosphatase Starch Excess4

Glucan phosphatases are members of a functionally diverse family of dual-specificity phosphatase (DSP) enzymes. The plant glucan phosphatase Starch Excess4 (SEX4) binds and dephosphorylates glucans, contributing to processive starch degradation in the chloroplast at night. Little is known about the complex kinetics of SEX4 when acting on its complex physiologically relevant glucan substrate. Therefore, we explored the kinetics of SEX4 against both insoluble starch and soluble amylopectin glucan substrates. SEX4 displays robust activity and a unique sigmoidal kinetic response to amylopectin, characterized by a Hill coefficient of 2.77 ± 0.63, a signature feature of cooperativity. We investigated the basis for this positive kinetic cooperativity and determined that the SEX4 carbohydrate-binding module (CBM) dramatically influences the binding cooperativity and substrate transformation rates. These findings provide insights into a previously unknown but important regulatory role for SEX4 in reversible starch phosphorylation and further advances our understanding of atypical kinetic mechanisms.

Homotropic Cooperativity of Midazolam Metabolism by Cytochrome P450 3A4: Insight from Computational Studies

Human cytochrome P450 3A4 (CYP3A4) is responsible for the metabolism of ∼50% clinically used drugs. Midazolam (MDZ) is a commonly used sedative drug and serves as a marker substrate for the CYP3A4 activity assessment. MDZ is metabolized by CYP3A4 to two hydroxylation products, 1′-OH-MDZ and 4-OH-MDZ. It has been reported that the ratio of 1′-OH-MDZ and 4-OH-MDZ is dependent on the MDZ concentration, which reflects the homotropic cooperative behavior in MDZ metabolism by CYP3A4. Here, we used quantum chemistry (QC), molecular docking, conventional molecular dynamics (cMD), and Gaussian accelerated molecular dynamics (GaMD) approaches to investigate the mechanism of the interactions between CYP3A4 and MDZ. QC calculations suggest that C1′ is less reactive for hydroxylation than C4, which is a pro-chirality carbon. However, the 4-OH-MDZ product is likely to be racemic due to the chirality inversion in the rebound step. The MD simulation results indicate that MDZ at the peripheral allosteric site is not stable and the binding modes of the MDZ molecules at the productive site are in line with the experimental observations.

Anti-inflammatory dopamine- and serotonin-based endocannabinoid epoxides reciprocally regulate cannabinoid receptors and the TRPV1 channel

The endocannabinoid system is a promising target to mitigate pain as the endocannabinoids are endogenous ligands of the pain-mediating receptors—cannabinoid receptors 1 and 2 (CB1 and CB2) and TRPV1. Herein, we report on a class of lipids formed by the epoxidation of N-arachidonoyl-dopamine (NADA) and N-arachidonoyl-serotonin (NA5HT) by epoxygenases. EpoNADA and epoNA5HT are dual-functional rheostat modulators of the endocannabinoid-TRPV1 axis. EpoNADA and epoNA5HT are stronger modulators of TRPV1 than either NADA or NA5HT, and epoNA5HT displays a significantly stronger inhibition on TRPV1-mediated responses in primary afferent neurons. Moreover, epoNA5HT is a full CB1 agonist. These epoxides reduce the pro-inflammatory biomarkers IL-6, IL-1β, TNF-α and nitrous oxide and raise anti-inflammatory IL-10 cytokine in activated microglial cells. The epoxides are spontaneously generated by activated microglia cells and their formation is potentiated in the presence of anandamide. Detailed kinetics and molecular dynamics simulation studies provide evidence for this potentiation using the epoxygenase human CYP2J2. Taken together, inflammation leads to an increase in the metabolism of NADA, NA5HT and other eCBs by epoxygenases to form the corresponding epoxides. The epoxide metabolites are bioactive lipids that are potent, multi-faceted molecules, capable of influencing the activity of CB1, CB2 and TRPV1 receptors.

Endocannabinoids are ligands of cannabinoid receptors and a promising target for pain management. Here, the authors report a class of lipids formed by the epoxidation of N-arachidonoyl dopamine and N-arachidonoyl serotonin by cytochrome P450 epoxygenases, which reciprocally regulate canabinoid receptors and display anti-inflammatory activity.

Explicit Treatment of Non-Michaelis-Menten and Atypical Kinetics in Early Drug Discovery*

Biological systems are highly regulated. They are also highly resistant to sudden perturbations enabling them to maintain the dynamic equilibrium essential to sustain life. This robustness is conferred by regulatory mechanisms that influence the activity of enzymes/proteins within their cellular context to adapt to changing environmental conditions. However, the initial rules governing the study of enzyme kinetics were mostly tested and implemented for cytosolic enzyme systems that were easy to isolate and/or recombinantly express. Moreover, these enzymes lacked complex regulatory modalities. Now, with academic labs and pharmaceutical companies turning their attention to more-complex systems (for instance, multiprotein complexes, oligomeric assemblies, membrane proteins and post-translationally modified proteins), the initial axioms defined by Michaelis-Menten (MM) kinetics are rendered inadequate, and the development of a new kind of kinetic analysis to study these systems is required. This review strives to present an overview of enzyme kinetic mechanisms that are atypical and, oftentimes, do not conform to the classical MM kinetics. Further, it presents initial ideas on the design and analysis of experiments in early drug-discovery for such systems, to enable effective screening and characterisation of small-molecule inhibitors with desirable physiological outcomes.

Meloxicam methyl group determines enzyme specificity for thiazole bioactivation compared to sudoxicam

Meloxicam is a thiazole-containing NSAID that was approved for marketing with favorable clinical outcomes despite being structurally similar to the hepatotoxic sudoxicam. Introduction of a single methyl group on the thiazole results in an overall lower toxic risk, yet the group's impact on P450 isozyme bioactivation is unclear. Through analytical methods, we used inhibitor phenotyping and recombinant P450s to identify contributing P450s, and then measured steady-state kinetics for bioactivation of sudoxicam and meloxicam by the recombinant P450s to determine relative efficiencies. Experiments showed that CYP2C8, 2C19, and 3A4 catalyze sudoxicam bioactivation, and CYP1A2 catalyzes meloxicam bioactivation, indicating that the methyl group not only impacts enzyme affinity for the drugs, but also alters which isozymes catalyze the metabolic pathways. Scaling of relative P450 efficiencies based on average liver concentration revealed that CYP2C8 dominates the sudoxicam bioactivation pathway and CYP2C9 dominates meloxicam detoxification. Dominant P450s were applied for an informatics assessment of electronic health records to identify potential correlations between meloxicam drug-drug interactions and drug-induced liver injury. Overall, our findings provide a cautionary tale on assumed impacts of even simple structural modifications on drug bioactivation while also revealing specific targets for clinical investigations of predictive factors that determine meloxicam-induced idiosyncratic liver injury.

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.

Mechanisms of Herb-Drug Interactions Involving Cinnamon and CYP2A6: Focus on Time-Dependent Inhibition by Cinnamaldehyde and 2-Methoxycinnamaldehyde

Information is scarce regarding pharmacokinetic-based herb-drug interactions (HDI) with trans-cinnamaldehyde (CA) and 2-methoxycinnamaldehyde (MCA), components of cinnamon. Given the presence of cinnamon in food and herbal treatments for various diseases, HDIs involving the CYP2A6 substrates nicotine and letrozole with MCA (KS = 1.58 µM; Hill slope = 1.16) and CA were investigated. The time-dependent inhibition (TDI) by MCA and CA of CYP2A6-mediated nicotine metabolism is a complex process involving multiple mechanisms. Molecular dynamic simulations showed that CYP2A6's active site accommodates two dynamic ligands. The preferred binding orientations for MCA and CA were consistent with the observed metabolism: epoxidation, O-demethylation, and aromatic hydroxylation of MCA and cinnamic acid formation from CA. The percent remaining activity plots for TDI by MCA and CA were curved, and they were analyzed with a numerical method using models of varying complexity. The best-fit models support multiple inactivator binding, inhibitor depletion, and partial inactivation. Deconvoluted mass spectra indicated that MCA and CA modified CYP2A6 apoprotein with mass additions of 156.79 (142.54-171.04) and 132.67 (123.37-141.98), respectively, and it was unaffected by glutathione. Heme degradation was observed in the presence of MCA (48.5% ± 13.4% loss; detected by liquid chromatography-tandem mass spectrometry). In the absence of clinical data, HDI predictions were made for nicotine and letrozole using inhibition parameters from the best-fit TDI models and parameters scaled from rats. Predicted area under the concentration-time curve fold changes were 4.29 (CA-nicotine), 4.92 (CA-letrozole), 4.35 (MCA-nicotine), and 5.00 (MCA-letrozole). These findings suggest that extensive exposure to cinnamon (corresponding to ≈ 275 mg CA) would lead to noteworthy interactions. SIGNIFICANCE STATEMENT: Human exposure to cinnamon is common because of its presence in food and cinnamon-based herbal treatments. Little is known about the risk for cinnamaldehyde and methoxycinnamaldehyde, two components of cinnamon, to interact with drugs that are eliminated by CYP2A6-mediated metabolism. The interactions with CYP2A6 are complex, involving multiple-ligand binding, time-dependent inhibition of nicotine metabolism, heme degradation, and apoprotein modification. An herb-drug interaction prediction suggests that extensive exposure to cinnamon would lead to noteworthy interactions with nicotine.

Conformational selection is present in ligand binding to cytochrome P450 19A1 lipoprotein nanodiscs

Cytochromes P450 (CYPs) display remarkable plasticity in their ability to bind substrates and catalyze a broad array of chemical reactions. Herein we evaluate binding of androstenedione, testosterone, and 7-hydroxyflavone to CYP19A1, also known as aromatase, in phospholipid nanodiscs by stopped-flow UV-vis spectroscopy. Exponential fitting of the kinetic traces supports the possibility of a multi-step binding mechanism. Subsequent global fitting of the data to the solutions of the coupled differential equations describing the fundamental mechanisms of induced fit and conformational selection, consistently support presence of the latter. To our knowledge, this is the first discrimination of conformational selection from induced fit for a mono-disperse CYP in a native-like membrane environment. In addition, 7-hydroxyflavone binds to CYP19A1 nanodiscs with comparable affinity to the substrates and induces an unusual spectral response likely attributable to hydrogen bonding to, rather than displacement of the heme-coordinated water molecule.

Time-dependent enzyme inactivation: Numerical analyses of in vitro data and prediction of drug-drug interactions

Cytochrome P450 (CYP) enzyme kinetics often do not conform to Michaelis-Menten assumptions, and time-dependent inactivation (TDI) of CYPs displays complexities such as multiple substrate binding, partial inactivation, quasi-irreversible inactivation, and sequential metabolism. Additionally, in vitro experimental issues such as lipid partitioning, enzyme concentrations, and inactivator depletion can further complicate the parameterization of in vitro TDI. The traditional replot method used to analyze in vitro TDI datasets is unable to handle complexities in CYP kinetics, and numerical approaches using ordinary differential equations of the kinetic schemes offer several advantages. Improvement in the parameterization of CYP in vitro kinetics has the potential to improve prediction of clinical drug-drug interactions (DDIs). This manuscript discusses various complexities in TDI kinetics of CYPs, and numerical approaches to model these complexities. The extrapolation of CYP in vitro TDI parameters to predict in vivo DDIs with static and dynamic modeling is discussed, along with a discussion on current gaps in knowledge and future directions to improve the prediction of DDI with in vitro data for CYP catalyzed drug metabolism.

Functional impact of cytochrome P450 3A (CYP3A) missense variants in cattle

Cytochrome P450 3A is the most important CYP subfamily in humans, and CYP3A4/CYP3A5 genetic variants contribute to inter-individual variability in drug metabolism. However, no information is available for bovine CYP3A (bCYP3A). Here we described bCYP3A missense single nucleotide variants (SNVs) and evaluated their functional effects. CYP3A28, CYP3A38 and CYP3A48 missense SNVs were identified in 300 bulls of Piedmontese breed through targeted sequencing. Wild-type and mutant bCYP3A cDNAs were cloned and expressed in V79 cells. CYP3A-dependent oxidative metabolism of testosterone (TST) and nifedipine (NIF) was assessed by LC-MS/MS. Finally, SNVs functional impact on TST hydroxylation was measured ex vivo in liver microsomes from individually genotyped animals. Thirteen missense SNVs were identified and validated. Five variants showed differences in CYP3A catalytic activity: three CYP3A28 SNVs reduced TST 6β-hydroxylation; one CYP3A38 variant increased TST 16β-hydroxylation, while a CYP3A48 SNV showed enhanced NIF oxidation. Individuals homozygous for rs384467435 SNV showed a reduced TST 6β-hydroxylation. Molecular modelling showed that most of SNVs were distal to CYP3A active site, suggesting indirect effects on the catalytic activity. Collectively, these findings demonstrate the importance of pharmacogenetics studies in veterinary species and suggest bCYP3A genotype variation might affect the fate of xenobiotics in food-producing species such as cattle.

Evaluation of aromatic hydrocarbon decomposition catalyzed by the dioxygenase system and substitution of ferredoxin and ferredoxin reductase

In this study, the catalytic activity and kinetic characteristics of the aromatic hydrocarbon dioxygenase system and the possibility of substituting its ferredoxin and ferredoxin reductase components were evaluated. The genes encoding toluene dioxygenase and toluene dihydrodiol dehydrogenase were cloned from Pseudomonas putida F1, and the corresponding enzymes were overexpressed and purified to homogeneity. Oxidative hydroxylation of toluene to cis-toluene dihydrodiol was catalyzed by toluene dioxygenase, and its subsequent dehydrogenation to 3-methylcatechol was catalyzed by toluene dihydrodiol dehydrogenase. The specific activity of the dioxygenase was 2.82 U/mg-protein, which is highly remarkable compared with the values obtained in previous researches conducted with crude extracts or insoluble forms of enzymes. Kinetic parameters, as characterized by the Hill equation, were v_max = 497.2 μM/min, K_M = 542.4 μM, and n_H = 2.2, suggesting that toluene dioxygenase has at least three cooperative binding sites for toluene. In addition, the use of alternative ferredoxins and reductases was examined. Ferredoxin cloned from CYP153 could transfer electrons to the iron sulfur protein component of toluene dioxygenase. The ferredoxin could be reduced by ferredoxin, rubredoxin, and putidaredoxin reductases of CYP153, alkane-1 monooxygenase, and camphor 5-monooxygenase, respectively. The results provide useful information regarding the effective enzymatic biotreatment of hazardous aromatic hydrocarbon contaminants.

Perspectives from the Innovation and Quality Consortium Induction Working Group on Factors Impacting Clinical Drug-Drug Interactions Resulting from Induction: Focus on Cytochrome 3A Substrates

A recent publication from the Innovation and Quality Consortium Induction Working Group collated a large clinical data set with the goal of evaluating the accuracy of drug-drug interaction (DDI) prediction from in vitro data. Somewhat surprisingly, comparison across studies of the mean- or median-reported area under the curve ratio showed appreciable variability in the magnitude of outcome. This commentary explores the possible drivers of this range of outcomes observed in clinical induction studies. While recommendations on clinical study design are not being proposed, some key observations were informative during the aggregate analysis of clinical data. Although DDI data are often presented using median data, individual data would enable evaluation of how differences in study design, baseline expression, and the number of subjects contribute. Since variability in perpetrator pharmacokinetics (PK) could impact the overall DDI interpretation, should this be routinely captured? Maximal induction was typically observed after 5-7 days of dosing. Thus, when the half-life of the inducer is less than 30 hours, are there benefits to a more standardized study design? A large proportion of CYP3A4 inducers were also CYP3A4 inhibitors and/or inactivators based on in vitro data. In these cases, using CYP3A selective substrates has limitations. More intensive monitoring of changes in area under the curve over time is warranted. With selective CYP3A substrates, the net effect was often inhibition, whereas less selective substrates could discern induction through mechanisms not susceptible to inhibition. The latter included oral contraceptives, which raise concerns of reduced efficacy following induction. Alternative approaches for modeling induction, such as applying biomarkers and physiologically based pharmacokinetic modeling (PBPK), are also considered. SIGNIFICANCE STATEMENT: The goal of this commentary is to stimulate discussion on whether there are opportunities to optimize clinical drug-drug interaction study design. The overall aim is to reduce, understand and contextualize the variability observed in the magnitude of induction across reported clinical studies. A large clinical CYP3A induction dataset was collected and further analyzed to identify trends and gaps. Reporting individual victim PK data, characterizing perpetrator PK and including additional PK assessments for mixed-mechanism perpetrators may provide insights into how these factors impact differences observed in clinical outcomes. The potential utility of biomarkers and PBPK modeling are discussed in considering future directions.

Impact of Lipid Partitioning on the Design, Analysis, and Interpretation of Microsomal Time-Dependent Inactivation

Nonspecific drug partitioning into microsomal membranes must be considered for in vitro-in vivo correlations. This work evaluated the effect of including lipid partitioning in the analysis of complex TDI kinetics with numerical methods. The covariance between lipid partitioning and multiple inhibitor binding was evaluated. Simulations were performed to test the impact of lipid partitioning on the interpretation of TDI kinetics, and experimental TDI datasets for paroxetine (PAR) and itraconazole (ITZ) were modeled. For most kinetic schemes, modeling lipid partitioning results in statistically better fits. For MM-IL simulations (KI,u = 0.1 µM, kinact = 0.1 minute-1), concurrent modeling of lipid partitioning for an fumic range (0.01, 0.1, and 0.5) resulted in better fits compared with post hoc correction (AICc: -526 vs. -496, -579 vs. -499, and -636 vs. -579, respectively). Similar results were obtained with EII-IL. Lipid partitioning may be misinterpreted as double binding, leading to incorrect parameter estimates. For the MM-IL datasets, when fumic = 0.02, MM-IL, and EII model fits were indistinguishable (δAICc = 3). For less partitioned datasets (fumic = 0.1 or 0.5), the inclusion of partitioning resulted in better models. The inclusion of lipid partitioning can lead to markedly different estimates of KI,u and kinact A reasonable alternate experimental design is nondilution TDI assays, with post hoc fumic incorporation. The best fit models for PAR (MIC-M-IL) and ITZ (MIC-EII-M-IL and MIC-EII-M-Seq-IL) were consistent with their reported mechanism and kinetics. Overall, experimental fumic values should be concurrently incorporated into TDI models with complex kinetics, when dilution protocols are used.

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.

Interaction of Human Drug-Metabolizing CYP3A4 with Small Inhibitory Molecules

Binding of small inhibitory compounds to human cytochrome P450 3A4 (CYP3A4) could interfere with drug metabolism and lead to drug-drug interactions, the underlying mechanism of which is not fully understood due to insufficient structural information. This study investigated the interaction of recombinant CYP3A4 with a nonspecific inhibitor metyrapone, antifungal drug fluconazole, and protease inhibitor phenylmethanesulfonyl fluoride (PMSF). Metyrapone and fluconazole are classic type II ligands that inhibit CYP3A4 with medium strength by ligating to the heme iron, whereas PMSF, lacking the heme-ligating moiety, acts as a weak type I ligand and inhibitor of CYP3A4. High-resolution crystal structures revealed that the orientation of metyrapone is similar but not identical to that in the previously reported 1W0G model, whereas the flexible fluconazole adapts a conformer markedly different from that observed in the target CYP51 enzymes, which could explain its high potential for cross-reactivity. Besides hydrophobic and aromatic interactions with the heme and active site residues, both drugs establish water-mediated contacts that stabilize the inhibitory complexes. PMSF also binds near the catalytic center, with the phenyl group parallel to the heme. However, it does not displace the water ligand and is held in place via strong H-bonds formed by the sulfofluoride moiety with Ser119 and Arg212. Collectively, our data suggest that PMSF might have multiple binding sites and likely occupies the high-affinity site in the crystal structure. Moreover, its hydrolysis product, phenylmethanesulfonic acid, can also access and be retained in the CYP3A4 active site. Therefore, to avoid experimental artifacts, PMSF should be excluded from purification and assay solutions.

Characterization of kinetics of human cytochrome P450s involved in bioactivation of flucloxacillin: inhibition of CYP3A-catalysed hydroxylation by sulfaphenazole

Background and Purpose

The aim of this study was to characterize the human cytochrome P450s (CYPs) involved in oxidative bioactivation of flucloxacillin to 5‐hydroxymethyl flucloxacillin, a metabolite with high cytotoxicity towards biliary epithelial cells.

Experimental Approach

The CYPs involved in hydroxylation of flucloxacillin were characterized using recombinant human CYPs, pooled liver microsomes in the presence of CYP‐specific inhibitors and by correlation analysis using a panel of liver microsomes from 16 donors.

Key Results

Recombinant CYPs showing the highest specific activity were CYP3A4, CYP3A7 and to lower extent CYP2C9 and CTP2C8. Michaelis–Menten enzyme kinetics were determined for pooled human liver microsomes, recombinant CYP3A4, CYP3A7 and CYP2C9. Surprisingly, sulfaphenazole appeared to be a potent inhibitor of 5′‐hydroxylation of flucloxacillin by both recombinant CYP3A4 and CYP3A7.

Conclusions and Implications

The combined results show that the 5′‐hydroxylation of flucloxacillin is primarily catalysed by CYP3A4, CYP3A7 and CYP2C9. The large variability of the hepatic expression of these enzymes could affect the formation of 5′‐hydroxymethyl flucloxacillin, which may determine the differences in susceptibility to flucloxacillin‐induced liver injury. Additionally, the strong inhibition in CYP3A‐catalysed flucloxacillin metabolism by sulfaphenazole suggests that unanticipated drug–drug interactions could occur with coadministered drugs.

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.

Complete Substrate Inhibition of Cytochrome P450 2C8 by AZD9496, an Oral Selective Estrogen Receptor Degrader

AZD9496 ((E)-3-(3,5-difluoro-4-((1R,3R)-2-(2-fluoro-2-methylpropyl)-3-methyl-2,3,4,9-tetrahydro-1H-pyrido[3,4-b]indol-1-yl)phenyl)acrylic acid) is an oral selective estrogen receptor degrader currently in clinical development for treatment of estrogen receptor-positive breast cancer. In a first-in-human phase 1 study, AZD9496 exhibited dose nonlinear pharmacokinetics, the mechanistic basis of which was investigated in this study. The metabolism kinetics of AZD9496 were studied using human liver microsomes (HLMs), recombinant cytochrome P450s (rP450s), and hepatocytes. In addition, modeling approaches were used to gain further mechanistic insights. CYP2C8 was predominantly responsible for biotransformation of AZD9496 to its two main metabolites whose rate of formation with increasing AZD9496 concentrations exhibited complete substrate inhibition in HLM, rCYP2C8, and hepatocytes. Total inhibition by AZD9496 of amodiaquine N-deethylation, a specific probe of CYP2C8 activity, confirmed the completeness of this inhibition. The commonly used substrate inhibition model analogous to uncompetitive inhibition fit poorly to the data. However, using the same model but without constraints on the number of molecules occupying the inhibitory binding site (i.e., nS₁ES) provided a significantly better fit (F test, P< 0.005). With the improved model, up to three AZD9496 molecules were predicted to bind the inhibitory site of CYP2C8. In contrast to previous studies showing substrate inhibition of P450s to be partial, our results demonstrate complete substrate inhibition of CYP2C8 via binding of more than one molecule of AZD9496 to the inhibitory site. As CYP2C8 appears to be the sole isoform catalyzing formation of the main metabolites, the substrate inhibition might explain the observed dose nonlinearity in the clinic at higher doses.

Improved Predictions of Drug-Drug Interactions Mediated by Time-Dependent Inhibition of CYP3A

Time-dependent inactivation (TDI) of cytochrome P450s (CYPs) is a leading cause of clinical drug-drug interactions (DDIs). Current methods tend to overpredict DDIs. In this study, a numerical approach was used to model complex CYP3A TDI in human-liver microsomes. The inhibitors evaluated included troleandomycin (TAO), erythromycin (ERY), verapamil (VER), and diltiazem (DTZ) along with the primary metabolites N-demethyl erythromycin (NDE), norverapamil (NV), and N-desmethyl diltiazem (NDD). The complexities incorporated into the models included multiple-binding kinetics, quasi-irreversible inactivation, sequential metabolism, inhibitor depletion, and membrane partitioning. The resulting inactivation parameters were incorporated into static in vitro-in vivo correlation (IVIVC) models to predict clinical DDIs. For 77 clinically observed DDIs, with a hepatic-CYP3A-synthesis-rate constant of 0.000 146 min-1, the average fold difference between the observed and predicted DDIs was 3.17 for the standard replot method and 1.45 for the numerical method. Similar results were obtained using a synthesis-rate constant of 0.000 32 min-1. These results suggest that numerical methods can successfully model complex in vitro TDI kinetics and that the resulting DDI predictions are more accurate than those obtained with the standard replot approach.

CW EPR parameters reveal cytochrome P450 ligand binding modes

Cytochrome P450 (CYP) monoxygenses utilize heme cofactors to catalyze oxidation reactions. They play a critical role in metabolism of many classes of drugs, are an attractive target for drug development, and mediate several prominent drug interactions. Many substrates and inhibitors alter the spin state of the ferric heme by displacing the heme's axial water ligand in the resting enzyme to yield a five-coordinate iron complex, or they replace the axial water to yield a nitrogen-ligated six-coordinate iron complex, which are traditionally assigned by UV-vis spectroscopy. However, crystal structures and recent pulsed electron paramagnetic resonance (EPR) studies find a few cases where molecules hydrogen bond to the axial water. The water-bridged drug-H₂O-heme has UV-vis spectra similar to nitrogen-ligated, six-coordinate complexes, but are closer to "reverse type I" complexes described in older liteature. Here, pulsed and continuous wave (CW) EPR demonstrate that water-bridged complexes are remarkably common among a range of nitrogenous drugs or drug fragments that bind to CYP3A4 or CYP2C9. Principal component analysis reveals a distinct clustering of CW EPR spectral parameters for water-bridged complexes. CW EPR reveals heterogeneous mixtures of ligated states, including multiple directly-coordinated complexes and water-bridged complexes. These results suggest that water-bridged complexes are under-represented in CYP structural databases and can have energies similar to other ligation modes. The data indicates that water-bridged binding modes can be identified and distinguished from directly-coordinated binding by CW EPR.

Mutational Analysis of a Conserved Glutamate Reveals Unique Mechanistic and Structural Features of the Phosphatase PRL-3

Phosphatase of regenerating liver (PRL)-3 (PTP4A3) has gained much attention in cancer research due to its involvement in tumor promoting and metastatic processes. It belongs to the protein tyrosine phosphatase (PTP) superfamily and is thought to follow the catalytic mechanism shared by this family, which aside from the conserved active-site amino acids includes a conserved glutamic acid residue that is usually required for the integrity of the active site in PTPs. We noted that in structures of PRL-3, PRL-1, and PTEN these residues do not clearly align and therefore we sought to investigate if the glutamic acid residue fulfills its usual function in these proteins. Although this residue was essential for PTEN's catalytic activity, it was nonessential for PRL-1 and PRL-3. Surprisingly, the mutation E50R increased PRL-3 activity against all tested in vitro substrates and also enhanced PRL-3-promoted cell adhesion and migration. We show that the introduction of Arg50 leads to an enhancement of substrate turnover for both PRL-3 and, to a lesser extent, PRL-1, and that the stronger gain in activity correlates with a higher structural flexibility of PRL-3, likely allowing for conformational adaptation during catalysis. Thus, in contrast to its crucial functions in other PTPs, this conserved glutamic acid can be replaced in PRL-3 without impairing the structural integrity. The variant with enhanced activity might serve as a tool to study PRL-3 in the future.

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.

Interaction and Transport of Methamphetamine and its Primary Metabolites by Organic Cation and Multidrug and Toxin Extrusion Transporters

Methamphetamine is one of the most abused illicit drugs with roughly 1.2 million users in the United States alone. A large portion of methamphetamine and its metabolites is eliminated by the kidney with renal clearance larger than glomerular filtration clearance. Yet the mechanism of active renal secretion is poorly understood. The goals of this study were to characterize the interaction of methamphetamine and its major metabolites with organic cation transporters (OCTs) and multidrug and toxin extrusion (MATE) transporters and to identify the major transporters involved in the disposition of methamphetamine and its major metabolites, amphetamine and para-hydroxymethamphetamine (p-OHMA). We used cell lines stably expressing relevant transporters to show that methamphetamine and its metabolites inhibit human OCTs 1-3 (hOCT1-3) and hMATE1/2-K with the greatest potencies against hOCT1 and hOCT2. Methamphetamine and amphetamine are substrates of hOCT2, hMATE1, and hMATE2-K, but not hOCT1 and hOCT3. p-OHMA is transported by hOCT1-3 and hMATE1, but not hMATE2-K. In contrast, organic anion transporters 1 and 3 do not interact with or transport these compounds. Methamphetamine and its metabolites exhibited complex interactions with hOCT1 and hOCT2, suggesting the existence of multiple binding sites. Our studies suggest the involvement of the renal OCT2/MATE pathway in tubular secretion of methamphetamine and its major metabolites and the potential of drug-drug interactions with substrates or inhibitors of the OCTs. This information may be considered when prescribing medications to suspected or known abusers of methamphetamine to mitigate the risk of increased toxicity or reduced therapeutic efficacy.

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.

Structural basis for regiospecific midazolam oxidation by human cytochrome P450 3A4

Human cytochrome P450 3A4 (CYP3A4) is a major hepatic and intestinal enzyme that oxidizes more than 60% of administered therapeutics. Knowledge of how CYP3A4 adjusts and reshapes the active site to regioselectively oxidize chemically diverse compounds is critical for better understanding structure-function relations in this important enzyme, improving the outcomes for drug metabolism predictions, and developing pharmaceuticals that have a decreased ability to undergo metabolism and cause detrimental drug-drug interactions. However, there is very limited structural information on CYP3A4-substrate interactions available to date. Despite the vast variety of drugs undergoing metabolism, only the sedative midazolam (MDZ) serves as a marker substrate for the in vivo activity assessment because it is preferentially and regioselectively oxidized by CYP3A4. We solved the 2.7 Å crystal structure of the CYP3A4-MDZ complex, where the drug is well defined and oriented suitably for hydroxylation of the C1 atom, the major site of metabolism. This binding mode requires H-bonding to Ser119 and a dramatic conformational switch in the F-G fragment, which transmits to the adjacent D, E, H, and I helices, resulting in a collapse of the active site cavity and MDZ immobilization. In addition to providing insights on the substrate-triggered active site reshaping (an induced fit), the crystal structure explains the accumulated experimental results, identifies possible effector binding sites, and suggests why MDZ is predominantly metabolized by the CYP3A enzyme subfamily.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Regorafenib induced severe toxic hepatitis: characterization and discussion

BACKGROUND

Regorafenib is the first small-molecule multikinase inhibitor which showed survival benefits in pretreated metastatic colorectal cancer (mCRC) patients. Besides classical adverse events of this drug class, hepatotoxicity has been described as a frequent side effect.

MATERIAL AND METHODS

Patients with refractory mCRC treated with regorafenib in our institution were reviewed. Severe treatment-related liver toxicity was investigated. Clinical history, liver histology and genetic assessment (sequence analysis) of cytochrome P3A4 (CYP3A4) and uridine diphosphate-glucuronosyltransferase 1A9 (UGT1A9) involved in regorafenib metabolization were here reported for patients with severe hepatotoxicity.

RESULTS

Among the 93 reviewed patients, 3 presented severe and icteric toxic hepatitis which was fatal for 1 patient. Histopathological liver lesions were different depending on the onset of hepatotoxicity (acute or subacute): acinar zone 3 necrosis in case of acute symptoms, and portal tract inflammation with porto-central bridging and fibrosis in the delayed presentation. None of the patients had CYP3A4 gene mutations. Similar polymorphisms in UGT1A9 gene promoter region (UGT1A9 variant -118T_9>10 [rs3832043]) were found in both patients who presented acute hepatitis. Moreover, it appears retrospectively that both of them already experienced significant toxicity under irinotecan-based chemotherapy.

CONCLUSION

This is the first report of severe hepatotoxicity with available liver histology and genetic assessment of enzymes involved in regorafenib metabolization. This report also reminds the importance of close liver tests monitoring during regorafenib treatment.