Allosteric P450 mechanisms: multiple binding sites, multiple conformers or both?

Effect of Lipid Bilayer Anchoring on the Conformational Properties of the Cytochrome P450 2D6 Binding Site

Human cytochrome P450 (CYP) proteins metabolize 75% of small-molecule pharmaceuticals, which makes structure-based modeling of CYP metabolism and inhibition, bolstered by the current availability of X-ray crystal structures of CYP globular catalytic domains, an attractive prospect. Accounting for this broad metabolic capacity is a combination of the existence of multiple different CYP proteins and the capacity of a single CYP protein to metabolize multiple different small molecules. It is thought that structural plasticity and flexibility contribute to this latter property; therefore, incorporating diverse conformational states of a particular CYP is likely an important consideration in structure-based CYP metabolism and inhibition modeling. While all-atom explicit-solvent molecular dynamics simulations can be used to generate conformational ensembles under biologically relevant conditions, existing CYP crystal structures are of the globular domain only, whereas human CYPs contain N-terminal transmembrane and linker peptides that anchor the globular catalytic domain to the endoplasmic reticulum. To determine whether this can cause significant differences in the sampled binding site conformations, microsecond scale all-atom explicit-solvent molecular dynamics simulations of the CYP2D6 globular domain in an aqueous environment were compared with those of the full-length protein anchored in a POPC lipid bilayer. While bilayer-anchoring damped some structural fluctuations in the globular domain relative to the aqueous simulations, none of the affected residues included binding site pocket residues. Furthermore, clustering of molecular dynamics snapshots based on either pairwise binding site pocket RMSD or volume differences demonstrated a lack of separation of snapshots from the two simulation conditions into different clusters. These results suggest the substantially simpler and computationally cheaper aqueous simulation approach can be used to generate a relevant conformational ensemble of the CYP2D6 binding site for structure-based metabolism and inhibition modeling.

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.

16R-HETE and 16S-HETE alter human cytochrome P450 1B1 enzyme activity probably through an allosteric mechanism

Cytochrome P450 1B1 (CYP1B1) has been widely associated with the development of cardiac pathologies due to its ability to produce cardiotoxic metabolites like midchain hydroxyeicosatetraenoic acids (HETEs) from arachidonic acid (AA) through an allylic oxidation reaction. 16-HETE is a subterminal HETE that is also produced by CYP-mediated AA metabolism. 19-HETE is another subterminal HETE that was found to inhibit CYP1B1 activity, lower midchain HETEs, and have cardioprotective effects. However, the effect of 16-HETE enantiomers on CYP1B1 has not yet been investigated. We hypothesized that 16(R/S)-HETE could alter the activity of CYP1B1 and other CYP enzymes. Therefore, this study was carried out to investigate the modulatory effect of 16-HETE enantiomers on CYP1B1 enzyme activity, and to examine the mechanisms by which they exert these modulatory effects. To investigate whether these effects are specific to CYP1B1, we also investigated 16-HETE modulatory effects on CYP1A2. Our results showed that 16-HETE enantiomers significantly increased CYP1B1 activity in RL-14 cells, recombinant human CYP1B1, and human liver microsomes, as seen by the significant increase in 7-ethoxyresorufin deethylation rate. On the contrary, 16-HETE enantiomers significantly inhibited CYP1A2 catalytic activity mediated by the recombinant human CYP1A2 and human liver microsomes. 16R-HETE showed stronger effects than 16S-HETE. The sigmoidal binding mode of the enzyme kinetics data demonstrated that CYP1B1 activation and CYP1A2 inhibition occurred through allosteric regulation. In conclusion, our study provides the first evidence that 16R-HETE and 16S-HETE increase CYP1B1 catalytic activity through an allosteric mechanism.

Pleiotropy of Progesterone Receptor Membrane Component 1 in Modulation of Cytochrome P450 Activity

Progesterone receptor membrane component 1 (PGRMC1) is one of few proteins that have been recently described as direct modulators of the activity of human cytochrome P450 enzymes (CYP)s. These enzymes form a superfamily of membrane-bound hemoproteins that metabolize a wide variety of physiological, dietary, environmental, and pharmacological compounds. Modulation of CYP activity impacts the detoxification of xenobiotics as well as endogenous pathways such as steroid and fatty acid metabolism, thus playing a central role in homeostasis. This review is focused on nine main topics that include the most relevant aspects of past and current PGRMC1 research, focusing on its role in CYP-mediated drug metabolism. Firstly, a general overview of the main aspects of xenobiotic metabolism is presented (I), followed by an overview of the role of the CYP enzymatic complex (IIa), a section on human disorders associated with defects in CYP enzyme complex activity (IIb), and a brief account of cytochrome b5 (cyt b5)'s effect on CYP activity (IIc). Subsequently, we present a background overview of the history of the molecular characterization of PGRMC1 (III), regarding its structure, expression, and intracellular location (IIIa), and its heme-binding capability and dimerization (IIIb). The next section reflects the different effects PGRMC1 may have on CYP activity (IV), presenting a description of studies on the direct effects on CYP activity (IVa), and a summary of pathways in which PGRMC1's involvement may indirectly affect CYP activity (IVb). The last section of the review is focused on the current challenges of research on the effect of PGRMC1 on CYP activity (V), presenting some future perspectives of research in the field (VI).

Unraveling the Ligand-Binding Sites of CYP3A4 by Molecular Dynamics Simulations with Solvent Probes

Cytochrome P450 3A4 (CYP3A4) is one of the most important drug-metabolizing enzymes in the human body and is well known for its complicated, atypical kinetic characteristics. The existence of multiple ligand-binding sites in CYP3A4 has been widely recognized as being capable of interfering with the active pocket through allosteric effects. The identification of ligand-binding sites other than the canonical active site above the heme is especially important for understanding the atypical kinetic characteristics of CYP3A4 and the intriguing association between the ligand and the receptor. In this study, we first employed mixed-solvent molecular dynamics (MixMD) simulations coupled with the online computational predictive tools to explore potential ligand-binding sites in CYP3A4. The MixMD approach demonstrates better performance in dealing with the receptor flexibility compared with other computational tools. From the sites identified by MixMD, we then picked out multiple sites for further exploration using ensemble docking and conventional molecular dynamics (cMD) simulations. Our results indicate that three extra sites are suitable for ligand binding in CYP3A4, including one experimentally confirmed site and two novel sites.

3site Multisubstrate-Bound State of Cytochrome P450cam

We identified a multisubstrate-bound state, hereby referred as a 3site state, in cytochrome P450cam via integrating molecular dynamics simulation with nuclear magnetic resonance (NMR) pseudocontact shift measurements. The 3site state is a result of simultaneous binding of three camphor molecules in three locations around P450cam: (a) in a well-established "catalytic" site near heme, (b) in a kink-separated "waiting" site along channel-1, and (c) in a previously reported "allosteric" site at E, F, G, and H helical junctions. These three spatially distinct binding modes in the 3site state mutually communicate with each other via homotropic allostery and act cooperatively to render P450cam functional. The 3site state shows a significantly superior fit with NMR pseudo contact shift (PCS) data with a Q-score of 0.045 than previously known bound states and consists of D251 free of salt-bridges with K178 and R186, rendering the enzyme functionally primed. To date, none of the reported cocomplex of P450cam with its redox partner putidaredoxin (pdx) has been able to match solution NMR data and controversial pdx-induced opening of P450cam's channel-1 remains a matter of recurrent discourse. In this regard, inclusion of pdx to the 3site state is able to perfectly fit the NMR PCS measurement with a Q-score of 0.08 and disfavors the pdx-induced opening of channel-1, reconciling previously unexplained remarkably fast hydroxylation kinetics with a koff of 10.2 s-1. Together, our findings hint that previous experimental observations may have inadvertently captured the 3site state as an in vitro solution state, instead of the catalytic state alone, and provided a distinct departure from the conventional understanding of cytochrome P450.

Impact of Heterotropic Allosteric Modulation on the Time-Dependent Inhibition of Cytochrome P450 3A4

The current study was designed to investigate the influence of allosteric effectors on the metabolism of the prototypical cytochrome P450 (CYP) 3A4 substrate midazolam (MDZ), and on the determination in vitro time-dependent inhibition (TDI) of CYP3A4 using human liver microsomes (HLM). As the concentration of midazolam increased to 250 µM in HLMs, homotropic cooperativity resulted in a decrease in the 1'-hydroxymidazolam to 4-hydroxymidazolam ratio to a maximum of 1.1. The presence of varying concentrations of testosterone, progesterone (PGS), or carbamazepine (CBZ) in HLMs with MDZ could recapitulate the effect of homotropic cooperativity such that the formation rates of the 1'hydroxymidazolam and 4-hydroxymidazolam were equal even at low concentrations of MDZ. The presence of PGS (10 or 100 µM) and CBZ (100 or 1000 µM) in in vitro TDI determination of four known CYP3A4 time-dependent inactivators (clarithromycin, troleandomycin, mibefradil, raloxifene) simultaneously decreased potency and inactivation rate constant, resulting in fold changes in inactivation efficiency on average of 1.6-fold and 13-fold for the low and high concentrations of allosteric modulator tested, respectively. The formation of a metabolic-intermediate complex (MIC) for clarithromycin and troleandomycin decreased in the presence of the allosteric modulators in a concentration-dependent manner, reaching a new steady state formation that could not be overcome with increased incubation time. Maximum reduction of the MIC formed by clarithromycin was up to ∼91%, while troleandomycin MIC decreased up to ∼31%. These findings suggest that the absence of endogenous allosteric modulators may contribute to the poor translation of HLM-based drug-drug interaction predictions. SIGNIFICANCE STATEMENT: The reported overprediction of in vitro human liver microsome time-dependent inhibition of CYP3A4 and observed drug interactions in vivo remains an issue in drug development. We provide characterization of allosteric modulators on the CYP3A4 metabolism of the prototypical substrate midazolam, demonstrating the ability of the modulators to recapitulate the homotropic cooperativity of midazolam. Furthermore, we demonstrate that allosteric heterotropic cooperativity of CYP3A4 can impact the time-dependent inhibition kinetics of known mechanisms-based inhibitors, providing a potential mechanism to explain the overprediction.

Functional characterization of CYP1 enzymes: Complex formation, membrane localization and function

CYP1A1, CYP1A2, and CYP1B1 have a high degree of sequence similarity, similar substrate selectivities and induction characteristics. However, experiments suggest that there are significant differences in their quaternary structures and function. The goal of this study was to characterize the CYP1 proteins regarding their ability to form protein-protein complexes, lipid microdomain localization, and ultimately function. This was accomplished by examining (1) substrate metabolism of the CYP1s as a function of NADPH-cytochrome P450 reductase (POR) concentration, and (2) quaternary structure, using bioluminescence resonance energy transfer (BRET). Both CYP1As were able to form BRET-detectable homomeric complexes, which was not observed with CYP1B1. When activities were measured as a function of [POR], CYP1A1 and CYP1B1 showed a hyperbolic response, consistent with mass-action binding; however, CYP1A2 produced a sigmoidal response, suggesting that the homomeric complex affected its function. Differences were observed in their ability to form heteromeric complexes. Whereas CYP1B1 and CYP1A1 formed a complex, neither the CYP1A1/CYP1A2 nor the CYP1B1/CYP1A2 pair formed BRET-detectable complexes. These proteins also differed in their lipid microdomain localization, with CYP1A2 and CYP1B1 residing in ordered membranes, and CYP1A1 in the disordered lipid regions. Taken together, despite their sequence similarities, there are substantial differences in quaternary structures and microdomain localization that can influence enzymatic activities. As these proteins exist in the endoplasmic reticulum with other ER-resident proteins, the P450s need to be considered as part of multi-enzyme systems rather than simply monomeric proteins interacting with their redox partners.

Structural basis of bidirectional allostery across the heme in a cytochrome P450 enzyme

Cytochromes P450 (CYPs) are heme-containing enzymes that are present in all kingdoms of life and share a structurally homologous, globular protein fold. CYPs utilize structures distal to the heme to recognize and coordinate substrates, while the necessary interactions with redox partner proteins are mediated at the opposite, proximal surface. In the current study, we investigated the functional allostery across the heme for the bacterial enzyme CYP121A1, which utilizes a non-polar distal-to-distal dimer interface for specific binding of its dicyclotyrosine substrate. Fluorine-detected Nuclear Magnetic Resonance (19F-NMR) spectroscopy was combined with site-specific labeling of a distal surface residue (S171C of the FG-loop), one residue of the B-helix (N84C), and two proximal surface residues (T103C and T333C) with a thiol-reactive fluorine label. Adrenodoxin was used as a substitute redox protein and was found to promote a closed arrangement of the FG-loop, similar to the addition of substrate alone. Disruption of the protein-protein interface by mutagenesis of two CYP121 basic surface residues removed the allosteric effect. Moreover, 19F-NMR spectra of the proximal surface indicate that ligand-induced allostery modulates the environment at the C-helix but not the meander region of the enzyme. In light of the high degree of structural homology in this family of enzymes, we interpret the findings from this work to represent a conserved allosteric network in CYPs.

Cooperative Substrate Binding Controls Catalysis in Bacterial Cytochrome P450terp (CYP108A1)

Despite being one of the most well-studied aspects of cytochrome P450 chemistry, important questions remain regarding the nature and ubiquity of allosteric regulation of catalysis. The crystal structure of a bacterial P450, P450terp, in the presence of substrate reveals two binding sites, one above the heme in position for regioselective hydroxylation and another in the substrate access channel. Unlike many bacterial P450s, P450terp does not exhibit an open to closed conformational change when substrate binds; instead, P450terp uses the second substrate molecule to hold the first substrate molecule in position for catalysis. Spectral titrations clearly show that substrate binding to P450terp is cooperative with a Hill coefficient of 1.4 and is supported by isothermal titration calorimetry. The importance of the allosteric site was explored by a series of mutations that weaken the second site and that help hold the first substrate in position for proper catalysis. We further measured the coupling efficiency of both the wild-type (WT) enzyme and the mutant enzymes. While the WT enzyme exhibits 97% efficiency, each of the variants showed lower catalytic efficiency. Additionally, the variants show decreased spin shifts upon binding of substrate. These results are the first clear example of positive homotropic allostery in a class 1 bacterial P450 with its natural substrate. Combined with our recent results from P450cam showing complex substrate allostery and conformational dynamics, our present study with P450terp indicates that bacterial P450s may not be as simple as once thought and share complex substrate binding properties usually associated with only mammalian P450s.

Cardiomyocyte-specific CYP2J2 and its therapeutic implications

INTRODUCTION

Cytochrome P450s (CYPs) are a superfamily of monooxygenases with diverse biological roles. CYP2J2 is an isozyme highly expressed in the heart where it metabolizes endogenous substrates such as N-3/N-6 polyunsaturated fatty acids (PUFA) to produce lipid mediators involved in homeostasis and cardioprotective responses. Expanding our knowledge of the role CYP2J2 has within the heart is important for understanding its impact on cardiac health and disease.

AREAS COVERED

The objective of this review was to assess the state of knowledge regarding cardiac CYP2J2. A literature search was conducted using PubMed-MEDLINE (from 2022 and earlier) to evaluate relevant studies regarding CYP2J2 mediated cardioprotection, small molecule modulators, effects of CYP2J2 substrates toward biologically relevant effects and implications of CYP2J2 polymorphisms and sexual dimorphism in the heart.

EXPERT OPINION

Cardiac CYP2J2-mediated metabolism of endogenous and exogenous substrates have been shown to impact cardiac function. Identifying individual factors, like sex and age, that affect CYP2J2 require further elucidation to better understand CYP2J2's clinical relevance. Resolving the biological targets and activities of CYP2J2-derived PUFA metabolites will be necessary to safely target CYP2J2 and design novel analogues. Targeting CYP2J2 for therapeutic aims offers a potential novel approach to regulating cardiac homeostasis, drug metabolism and cardioprotection.

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.

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.

Harnessing P450 Enzyme for Biotechnology and Synthetic Biology

Cytochrome P450 enzymes (P450s, CYPs) catalyze the oxidative transformation of a wide range of organic substrates. Their functions are crucial to xenobiotic metabolism and steroid transformation in humans and other organisms. The enzymes are promising for synthetic biology applications but limited by several drawbacks including low turnover rates, poor stability, the dependance of expensive cofactors and redox partners, and the narrow substrate scope. To conquer these obstacles, emerging strategies including substrate engineering, usage of decoy and decoy-based small molecules auxiliaries, designing of artificial enzyme cascades and the incorporation of materials have been explored based on the unique properties of P450s. These strategies can be applied to a wide range of P450s and can be combined with protein engineering to improve the enzymatic activities. This minireview will focus on some recent developments of these strategies which have been used to leverage P450 catalysis. Remaining challenges and future opportunities will also be discussed.

Development of Robust Quantitative Structure-Activity Relationship Models for CYP2C9, CYP2D6, and CYP3A4 Catalysis and Inhibition

Cytochrome P450 enzymes are responsible for the metabolism of >75% of marketed drugs, making it essential to identify the contributions of individual cytochromes P450 to the total clearance of a new candidate drug. Overreliance on one cytochrome P450 for clearance levies a high risk of drug-drug interactions; and considering that several human cytochrome P450 enzymes are polymorphic, it can also lead to highly variable pharmacokinetics in the clinic. Thus, it would be advantageous to understand the likelihood of new chemical entities to interact with the major cytochrome P450 enzymes at an early stage in the drug discovery process. Typical screening assays using human liver microsomes do not provide sufficient information to distinguish the specific cytochromes P450 responsible for clearance. In this regard, we experimentally assessed the metabolic stability of ∼5000 compounds for the three most prominent xenobiotic metabolizing human cytochromes P450, i.e., CYP2C9, CYP2D6, and CYP3A4, and used the data sets to develop quantitative structure-activity relationship models for the prediction of high-clearance substrates for these enzymes. Screening library included the NCATS Pharmaceutical Collection, comprising clinically approved low-molecular-weight compounds, and an annotated library consisting of drug-like compounds. To identify inhibitors, the library was screened against a luminescence-based cytochrome P450 inhibition assay; and through crossreferencing hits from the two assays, we were able to distinguish substrates and inhibitors of these enzymes. The best substrate and inhibitor models (balanced accuracies ∼0.7), as well as the data used to develop these models, have been made publicly available (https://opendata.ncats.nih.gov/adme) to advance drug discovery across all research groups. SIGNIFICANCE STATEMENT: In drug discovery and development, drug candidates with indiscriminate cytochrome P450 metabolic profiles are considered advantageous, since they provide less risk of potential issues with cytochrome P450 polymorphisms and drug-drug interactions. This study developed robust substrate and inhibitor quantitative structure-activity relationship models for the three major xenobiotic metabolizing cytochromes P450, i.e., CYP2C9, CYP2D6, and CYP3A4. The use of these models early in drug discovery will enable project teams to strategize or pivot when necessary, thereby accelerating drug discovery research.

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.

Enantioselective inhibition of human CYP2C19 by the chiral pesticide ethofumesate: Prediction of pesticide-drug interactions in humans

Ethofumesate is a chiral herbicide that may display enantioselective behavior in humans. For this reason, the enantioselective potential of ethofumesate and its main metabolite ethofumesate-2-hydroxy to cause pesticide-drug interactions on cytochrome P450 forms (CYPs) has been evaluated by using human liver microsomes. Among the evaluated CYPs, CYP2C19 had its activity decreased by the ethofumesate racemic mixture (rac-ETO), (+)-ethofumesate ((+)-ETO), and (-)-ethofumesate ((-)-ETO). CYP2C19 inhibition was not time-dependent, but a strong inhibition potential was observed for rac-ETO (IC₅₀ = 5 ± 1 μmol L^-1), (+)-ETO (IC₅₀ = 1.6 ± 0.4 μmol L^-1), and (-)-ETO (IC₅₀ = 1.8 ± 0.4 μmol L^-1). The reversible inhibition mechanism was competitive, and the inhibition constant (K_i) values for rac-ETO (2.6 ± 0.4 μmol L^-1), (+)-ETO (1.5 ± 0.2 μmol L^-1), and (-)-ETO (0.7 ± 0.1 μmol L^-1) were comparable to the K_i values of strong CYP2C19 inhibitors. Inhibition of CYP2C19 by ethofumesate was enantioselective, being almost twice higher for (-)-ETO than for (+)-ETO, which indicates that this enantiomer may be a more potent inhibitor of this CYP form. For an in vitro-in vivo correlation, the Food and Drug Administration's (FDA) guideline on the assessment of drug-drug interactions used in the early stages of drug development was used. The FDA's R₁ values were estimated on the basis of the obtained ethofumesate K_i and distribution volume, metabolism, unbound plasma fraction, gastrointestinal and dermal absorption data available in the literature. The correlation revealed that ethofumesate probably inhibits CYP2C19 in vivo for both chronic (oral) and occupational (dermal) exposure scenarios.

Structural Basis for the Diminished Ligand Binding and Catalytic Ability of Human Fetal-Specific CYP3A7

Cytochrome P450 3A7 (CYP3A7) is a fetal/neonatal liver enzyme that participates in estriol synthesis, clearance of all-trans retinoic acid, and xenobiotic metabolism. Compared to the closely related major drug-metabolizing enzyme in adult liver, CYP3A4, the ligand binding and catalytic capacity of CYP3A7 are substantially reduced. To better understand the structural basis for these functional differences, the 2.15 Å crystal structure of CYP3A7 has been solved. Comparative analysis of CYP3A enzymes shows that decreased structural plasticity rather than the active site microenvironment defines the ligand binding ability of CYP3A7. In particular, a rotameric switch in the gatekeeping amino acid F304 triggers local and long-range rearrangements that transmit to the F-G fragment and alter its interactions with the I-E-D-helical core, resulting in a more rigid structure. Elongation of the β₃-β₄ strands, H-bond linkage in the substrate channel, and steric constraints in the C-terminal loop further increase the active site rigidity and limit conformational ensemble. Collectively, these structural distinctions lower protein plasticity and change the heme environment, which, in turn, could impede the spin-state transition essential for optimal reactivity and oxidation of substrates.

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.

Insights into oral bioavailability enhancement of therapeutic herbal constituents by cytochrome P450 3A inhibition

Herbal plants typically have complex compositions and diverse mechanisms. Among them, bioactive constituents with relatively high exposure in vivo are likely to exhibit therapeutic efficacy. On the other hand, their bioavailability may be influenced by the synergistic effects of different bioactive components. Cytochrome P450 3A (CYP3A) is one of the most abundant CYP enzymes, responsible for the metabolism of 50% of approved drugs. In recent years, many therapeutic herbal constituents have been identified as CYP3A substrates. It is more evident that CYP3A inhibition derived from the herbal formula plays a critical role in improving the oral bioavailability of therapeutic constituents. CYP3A inhibition may be the mechanism of the synergism of herbal formula. In this review, we explored the multiplicity of CYP3A, summarized herbal monomers with CYP3A inhibitory effects, and evaluated herb-mediated CYP3A inhibition, thereby providing new insights into the mechanisms of CYP3A inhibition-mediated oral herb bioavailability.

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.

Effects of alcohol-induced increase in CYP2E1 content in human liver microsomes on the activity and cooperativity of CYP3A4

We investigate the effect of the alcohol-induced increase in the content of CYP2E1 in human liver microsomes (HLM) on the function of CYP3A4. Membrane incorporation of the purified CYP2E1 into HLM considerably increases the rate of metabolism of 7-benzyloxyquinoline (BQ) and attenuates the homotropic cooperativity observed with this CYP3A4-specific substrate. It also eliminates the activating effect of α-naphthoflavone (ANF) seen in some HLM samples. To probe the physiological relevance of these effects, we compared three pooled preparations of HLM from normal donors (HLM-N) with a pooled preparation from ten heavy alcohol consumers (HLM-A). The composition of the P450 pool in all samples was characterized by the mass-spectrometric determination of 11 cytochrome P450 species. The fractional content of CYP2E1 in HLM-A was from 2.0 to 3.4 times higher than in HLM-N. In contrast, the content of CYP3A4 in HLM-A was the lowest among all samples. Despite that, HLM-A exhibited a much higher metabolism rate and a lower homotropic cooperativity with BQ, similar to CYP2E1-enriched HLM-N. To substantiate the involvement of interactions between CYP2E1 and CYP3A4 in these effects, we probed hetero-association of these proteins in CYP3A4-containing Supersomes™ with a technique employing CYP2E1 labeled with BODIPY-618 maleimide. These experiments evinced the interactions between the two enzymes and revealed an inhibitory effect of ANF on their association. Our results demonstrate that the functional properties of CYP3A4 are fundamentally dependent on the composition of the cytochrome P450 ensemble and suggest a possible impact of chronic alcohol exposure on the pharmacokinetics of drugs metabolized by CYP3A4.

Demystifying Chronic Kidney Disease of Unknown Etiology (CKDu): Computational Interaction Analysis of Pesticides and Metabolites with Vital Renal Enzymes

Chronic kidney disease of unknown etiology (CKDu) has been recognized as a global non-communicable health issue. There are many proposed risk factors for CKDu and the exact reason is yet to be discovered. Understanding the inhibition or manipulation of vital renal enzymes by pesticides can play a key role in understanding the link between CKDu and pesticides. Even though it is very important to take metabolites into account when investigating the relationship between CKDu and pesticides, there is a lack of insight regarding the effects of pesticide metabolites towards CKDu. In this study, a computational approach was used to study the effects of pesticide metabolites on CKDu. Further, interactions of selected pesticides and their metabolites with renal enzymes were studied using molecular docking and molecular dynamics simulation studies. It was evident that some pesticides and metabolites have affinity to bind at the active site or at regulatory sites of considered renal enzymes. Another important discovery was the potential of some metabolites to have higher binding interactions with considered renal enzymes compared to the parent pesticides. These findings raise the question of whether pesticide metabolites may be a main risk factor towards CKDu.

Nonadditivity in human microsomal drug metabolism revealed in a study with coumarin 152, a polyspecific cytochrome P450 substrate

We closely characterized 7-Dimethylamino-4-trifluromethylcoumarin (Coumarin 152, C152), a substrate metabolized by multiple P450 species, to establish a new fluorogenic probe for the studies of functional integration in the cytochrome P450 ensemble. Scanning fluorescence spectroscopy and LC/MS-MS were used to characterize the products of N-demethylation of C152 and optimize their fluorometric detection. The metabolism of C152 by the individual P450 species was characterized using the microsomes containing cDNA-expressed enzymes. C152 metabolism in human liver microsomes (HLM) was studied in a preparation with quantified content of eleven P450 species. C152 is metabolized by CYP2B6, CYP3A4, CYP3A5, CYP2C19, CYP1A2, CYP2C9, and CYP2C8 listed in the order of decreasing turnover. The affinities exhibited by CYP3A5, CYP2C9, and CYP2C8 were lower than those characteristic to the other enzymes. The presumption of additivity suggests the participation of CYP3A4, CYP2B6, and CYP2C19 to be 84, 8, and 0.2%, respectively. Contrary to this prediction, inhibitory analysis identified CYP2C19 as the principal C152-metabolizing enzyme. We thoroughly characterize C152 for the studies of drug metabolism in HLM and demonstrate the limitations of the proportional projection approach by providing an example, where the involvement of individual P450 species cannot be predicted from their content.

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.

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.

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.

Kinetic mechanism of time-dependent inhibition of CYP2D6 by 3,4-methylenedioxymethamphetamine (MDMA): Functional heterogeneity of the enzyme and the reversibility of its inactivation

We investigate the mechanism of time-dependent inhibition (TDI) of human cytochrome P450 2D6 (CYP2D6) by 3,4-methylenedioxymethamphetamine (MDMA, ecstasy), one of the most widespread recreational drugs of abuse. In an effort to unravel the kinetic mechanism of the formation of metabolic inhibitory complex (MIC) of CYP2D6 with MDMA-derived carbene we carried out a series of spectrophotometric studies paralleled with registration of the kinetics of time-dependent inhibition (TDI) in CYP2D6-incorporated proteoliposomes. The high amplitude of spectral signal in this system allowed us to characterize the spectral properties of the formed MIC in details and obtain an accurate spectral signature of MIC formation. This information was then used in the studies with CYP2D6-containing microsomes of insect cells (CYP2D6 Supersomes™). Our results demonstrate that in both systems the formation of the ferrous carbene-derived MIC is relatively slow, reversible and is not associated with the accumulation of the ferric carbene intermediate, as takes place in the case of CYP3A4 and podophylotoxin. Furthermore, the limited amplitude of MIC formation suggests that only a fraction (∼50%) of spectrally detectable CYP2D6 in both proteoliposomes and Supersomes participates in the formation of MIC and is therefore involved in the MDMA metabolism. This observation reveals yet another example of a cytochrome P450 that exhibits persistent functional heterogeneity of its population in microsomal membranes. Our study provides a solid methodological background for further mechanistic studies of MIC formation in human liver microsomes and demonstrates that the potency and physiological relevance of MDMA-dependent TDI of CYP2D6 may be overestimated.

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.

The functional effects of physical interactions involving cytochromes P450: putative mechanisms of action and the extent of these effects in biological membranes

Cytochromes P450 represent a family of enzymes, which are responsible for the oxidative metabolism of a wide variety of xenobiotics. Although the mammalian P450s require interactions with their redox partners in order to function, more recently, P450 system proteins have been shown to exist as multi-protein complexes that include the formation of P450•P450 complexes. Evidence has shown that the metabolism of some substrates by a given P450 can be influenced by the specific interaction of the enzyme with other forms of P450. Detailed kinetic analysis of these reactions in vitro has shown that the P450-P450 interactions can alter metabolism by changing the ability of a P450 to bind to its cognate redox partner, NADPH-cytochrome P450 reductase; by altering substrate binding to the affected P450; and/or by changing the rate of a catalytic step of the reaction cycle. This review summarizes the known examples of P450-P450 interactions that have been shown in vitro to influence metabolism and categorizes them according to the mechanism(s) causing the effects. P450-P450 interactions have the potential to cause major changes in the metabolism and elimination of drugs in vivo. This review summarizes the evidence that the P450-P450 interactions influence metabolism in biological membranes and discusses the studies, which will provide further insight into the extent of these effects 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.

Deciphering common recognition principles of nucleoside mono/di and tri-phosphates binding in diverse proteins via structural matching of their binding sites

Nucleoside triphosphate (NTP) ligands are of high biological importance and are essential for all life forms. A pre-requisite for them to participate in diverse biochemical processes is their recognition by diverse proteins. It is thus of great interest to understand the basis for such recognition in different proteins. Towards this, we have used a structural bioinformatics approach and analyze structures of 4677 NTP complexes available in Protein Data Bank (PDB). Binding sites were extracted and compared exhaustively using PocketMatch, a sensitive in-house site comparison algorithm, which resulted in grouping the entire dataset into 27 site-types. Each of these site-types represent a structural motif comprised of two or more residue conservations, derived using another in-house tool for superposing binding sites, PocketAlign. The 27 site-types could be grouped further into 9 super-types by considering partial similarities in the sites, which indicated that the individual site-types comprise different combinations of one or more site features. A scan across PDB using the 27 structural motifs determined the motifs to be specific to NTP binding sites, and a computational alanine mutagenesis indicated that residues identified to be highly conserved in the motifs are also most contributing to binding. Alternate orientations of the ligand in several site-types were observed and rationalized, indicating the possibility of some residues serving as anchors for NTP recognition. The presence of multiple site-types and the grouping of multiple folds into each site-type is strongly suggestive of convergent evolution. Knowledge of determinants obtained from this study will be useful for detecting function in unknown proteins. Proteins 2017; 85:1699-1712. © 2017 Wiley Periodicals, Inc.

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.

Conformational Mobility in Cytochrome P450 3A4 Explored by Pressure-Perturbation EPR Spectroscopy

We used high hydrostatic pressure as a tool for exploring the conformational landscape of human cytochrome P450 3A4 (CYP3A4) by electron paramagnetic resonance and fluorescence spectroscopy. Site-directed incorporation of a luminescence resonance energy transfer donor-acceptor pair allowed us to identify a pressure-dependent equilibrium between two states of the enzyme, where an increase in pressure increased the spatial separation between the two distantly located fluorophores. This transition is characterized by volume change (ΔV°) and P1/2 values of -36.8 ± 5.0 mL/mol and 1.45 ± 0.33 kbar, respectively, which corresponds to a Keq° of 0.13 ± 0.06, so that only 15% of the enzyme adopts the pressure-promoted conformation at ambient pressure. This pressure-promoted displacement of the equilibrium is eliminated by the addition of testosterone, an allosteric activator. Using site-directed spin labeling, we demonstrated that the pressure- and testosterone-sensitive transition is also revealed by pressure-induced changes in the electron paramagnetic resonance spectra of a nitroxide side chain placed at position 85 or 409 of the enzyme. Furthermore, we observed a pressure-induced displacement of the emission maxima of a solvatochromic fluorophore (7-diethylamino-3-((((2-maleimidyl)ethyl)amino)carbonyl) coumarin) placed at the same positions, which suggests a relocation to a more polar environment. Taken together, the results reveal an effector-dependent conformational equilibrium between open and closed states of CYP3A4 that involves a pronounced change at the interface between the region of α-helices A/A' and the meander loop of the enzyme, where residues 85 and 409 are located. Our study demonstrates the high potential of pressure-perturbation strategies for studying protein conformational landscapes.

Structural basis of steroid binding and oxidation by the cytochrome P450 CYP109E1 from Bacillus megaterium

Cytochrome P450 monooxygenases (P450s) are attractive enzymes for the pharmaceutical industry, in particular, for applications in steroidal drug synthesis. Here, we report a comprehensive functional and structural characterization of CYP109E1, a novel steroid‐converting cytochrome P450 enzyme identified from the genome of Bacillus megaterium DSM319. In vitro and whole‐cell in vivo turnover experiments, combined with binding assays, revealed that CYP109E1 is able to hydroxylate testosterone at position 16β. Related steroids with bulky substituents at carbon C17, like corticosterone, bind to the enzyme without being converted. High‐resolution X‐ray structures were solved of a steroid‐free form of CYP109E1 and of complexes with testosterone and corticosterone. The structural analysis revealed a highly dynamic active site at the distal side of the heme, which is wide open in the absence of steroids, can bind four ordered corticosterone molecules simultaneously, and undergoes substantial narrowing upon binding of single steroid molecules. In the crystal structures, the single bound steroids adopt unproductive binding modes coordinating the heme‐iron with their C3‐keto oxygen. Molecular dynamics (MD) simulations suggest that the steroids may also bind in ~180° reversed orientations with the C16 carbon and C17‐substituents pointing toward the heme, leading to productive binding of testosterone explaining the observed regio‐ and stereoselectivity. The X‐ray structures and MD simulations further identify several residues with important roles in steroid binding and conversion, which could be confirmed by site‐directed mutagenesis. Taken together, our results provide unique insights into the CYP109E1 activity, substrate specificity, and regio/stereoselectivity.

Database

The atomic coordinates and structure factors have been deposited in the Protein Data Bank with accession codes 5L90 (steroid‐free CYP109E1), 5L91 (CYP109E1‐COR4), 5L94 (CYP109E1‐TES), and 5L92 (CYP109E1‐COR).

Enzymes

Cytochrome P450 monooxygenase CYP109E1, EC 1.14.14.1, UniProt ID: D5DKI8, Adrenodoxin reductase EC 1.18.1.6.

Quantum Chemical and Docking Insights into Bioavailability Enhancement of Curcumin by Piperine in Pepper

We combine quantum chemical and molecular docking techniques to provide new insights into how piperine molecule in various forms of pepper enhances bioavailability of a number of drugs including curcumin in turmeric for which it increases its bioavailability by a 20-fold. We have carried out docking studies of quantum chemically optimized piperine structure binding to curcumin, CYP3A4 in cytochrome P450, p-Glycoprotein and UDP-glucuronosyltransferase (UGT), the enzyme responsible for glucuronosylation, which increases the solubility of curcumin. All of these studies establish that piperine binds to multiple sites on the enzymes and also intercalates with curcumin forming a hydrogen bonded complex with curcumin. The conjugated network of double bonds and the presence of multiple charge centers of piperine offer optimal binding sites for piperine to bind to enzymes such as UDP-GDH, UGT, and CYP3A4. Piperine competes for curcumin's intermolecular hydrogen bonding and its stacking propensity by hydrogen bonding with enolic proton of curcumin. This facilitates its metabolic transport, thereby increasing its bioavailability both through intercalation into curcumin layers through intermolecular hydrogen bonding, and by inhibiting enzymes that cause glucuronosylation of curcumin.

Coumarin Derivatives as Substrate Probes of Mammalian Cytochromes P450 2B4 and 2B6: Assessing the Importance of 7-Alkoxy Chain Length, Halogen Substitution, and Non-Active Site Mutations

Using a combined structural and biochemical approach, the functional importance of a recently described peripheral pocket bounded by the E-, F-, G-, and I-helices in CYP2B4 and 2B6 was probed. Three series of 4-substituted-7-alkoxycoumarin derivatives with -H, -CH3, or -CF3 at the 4 position of the coumarin core were used initially to monitor functional differences between CYP2B4 and 2B6. 7-Ethoxy-4-(trifluoromethyl)coumarin (7-EFC) displayed the highest catalytic efficiency among these substrates. Mutants were made to alter side-chain polarity (V/E194Q) or bulk (F/Y244W) to alter access to the peripheral pocket. Modest increases in catalytic efficiency of 7-EFC O-deethylation by the mutants were magnified considerably by chlorination or bromination of the substrate ethoxy chain. A structure of CYP2B6 Y244W in complex with (+)-α-pinene was solved at 2.2 Å and showed no CYMAL-5 in the peripheral pocket. A ligand free structure of CYP2B4 F244W was solved at 3.0 Å with CYMAL-5 in the peripheral pocket. In both instances, comparison of the respective wild-type and mutant CYP2B enzymes revealed that CYMAL-5 occupancy of the peripheral pocket had little effect on the topology of active site residue side-chains, despite the fact that the peripheral pocket and active site are located on opposite sides of the I-helix. Analysis of available CYP2B structures suggest that the effect of the amino acid substitutions within the peripheral pocket derive from altered interactions between the F and G helices.

Effect of Ethanol on the Metabolic Characteristics of HIV-1 Integrase Inhibitor Elvitegravir and Elvitegravir/Cobicistat with CYP3A: An Analysis Using a Newly Developed LC-MS/MS Method

Elvitegravir (EVG), an integrase inhibitor for the treatment HIV infection, is increasingly becoming the part of first-line antiretroviral therapy (ART) regimen. EVG is mainly metabolized through cytochrome P450 (CYP) 3A4. Previously, we have shown that ethanol alters ART-CYP3A4 interactions with protease inhibitors thereby altering their metabolisms. However, as EVG is a fairly new class of drug, its kinetic characteristics and the effect of ethanol on EVG-CYPP3A4 interaction is poorly understood. In this study, we characterized EVG and cobicistat (COBI)-boosted EVG metabolism in human microsomes followed by ethanol-EVG, ethanol-COBI-EVG interaction with CYP3A. First, we developed and validated a simple, sensitive, and robust liquid chromatography–tandem mass spectrometry (LC-MS/MS) method for the quantification of EVG in the human liver microsomes. The lower limit of quantification for the drug was at 0.003 μM (1.34ng/ml). Extraction yield, matrix effects, drug stability, and calibration curves for the proposed method were validated according to the FDA guidelines. Time dependent kinetics data showed that 20mM ethanol decreases the apparent half-life of EVG degradation by ~50% compared to EVG alone. Our substrate kinetic results revealed that ethanol mildly decreases the catalytic efficiency for EVG metabolism. Inhibition studies demonstrated that EVG inhibits CYP3A4, and 20 mM ethanol causes a decrease in the IC₅₀ of EVG. However, in the presence of COBI we were unable to determine these parameters effectively because COBI, being a strong inhibitor of CYP3A4, blocked the EVG/ethanol-CYP3A4 interactions. Docking studies predicted a shift of EVG or COBI binding to the active site of CYP3A4 in the presence of ethanol. Taken together, these results suggest that ethanol interacts with microsomal CYP3A and alters EVG-CYP3A4 interaction thereby altering EVG metabolism and inhibition of CYP3A4 by EVG. This finding has clinical significance because alcohol use is highly prevalent in HIV population, and there are no separate guidelines for these patients while they are on ART medication.

The use of isomeric testosterone dimers to explore allosteric effects in substrate binding to cytochrome P450 CYP3A4

Cytochrome P450 CYP3A4 is the main drug-metabolizing enzyme in the human liver, being responsible for oxidation of 50% of all pharmaceuticals metabolized by human P450 enzymes. Possessing a large substrate binding pocket, it can simultaneously bind several substrate molecules and often exhibits a complex pattern of drug-drug interactions. In order to better understand structural and functional aspects of binding of multiple substrate molecules to CYP3A4 we used resonance Raman and UV-VIS spectroscopy to document the effects of binding of synthetic testosterone dimers of different configurations, cis-TST2 and trans-TST2. We directly demonstrate that the binding of two steroid molecules, which can assume multiple possible configurations inside the substrate binding pocket of monomeric CYP3A4, can lead to active site structural changes that affect functional properties. Using resonance Raman spectroscopy, we have documented perturbations in the ferric and Fe-CO states by these substrates, and compared these results with effects caused by binding of monomeric TST. While the binding of trans-TST2 yields results similar to those obtained with monomeric TST, the binding of cis-TST2 is much tighter and results in significantly more pronounced conformational changes of the porphyrin side chains and Fe-CO unit. In addition, binding of an additional monomeric TST molecule in the remote allosteric site significantly improves binding affinity and the overall spin shift for CYP3A4 with trans-TST2 dimer bound inside the substrate binding pocket. This result provides the first direct evidence for an allosteric effect of the peripheral binding site at the protein-membrane interface on the functional properties of CYP3A4.

[Molecular organization of the microsomal oxidative system: a new connotation for an old term]

The central role that cytochromes P450 play in the metabolism of drugs and other xenobiotics makes these enzymes a major subject for studies of drug disposition, adverse drug effects and drug-drug interactions. Although there has been tremendous success in delineating P450 mechanisms, the concept of the drug-metabolizing ensemble as a functionally integrated system remains undeveloped. However, eukaryotic cells typically possess a multitude of different P450 enzymes that are co-localized in the membrane of endoplasmic reticulum (ER) and interact with each other with the formation of dynamic heteromeric complexes (mixed oligomers). Appreciation of the importance of developing an integral, systems approach to the ensemble of cytochromes P450 as an integral system inspired growing interest of researchers to the molecular organization of microsomal monooxygenase, which remained in the focus of research of academician Archakov for over 40 years. Fundamental studies carried out under his guidance have an important impact on our current concepts in this area. Further exploration of the molecular organization of the system of microsomal monooxygenase as an integral multienzyme and multifunctional system will have an essential impact on our understanding of the key factors that determine the changes in human drug metabolism and other P450-related functions in development, aging, and disease, as well as under influence of drugs, food ingredients, and environmental contaminants.

Highly Efficient CYP167A1 (EpoK) dependent Epothilone B Formation and Production of 7-Ketone Epothilone D as a New Epothilone Derivative

Since their discovery in the soil bacterium Sorangium cellulosum, epothilones have emerged as a valuable substance class with promising anti-tumor activity. Because of their benefits in the treatment of cancer and neurodegenerative diseases, epothilones are targets for drug design and pharmaceutical research. The final step of their biosynthesis – a cytochrome P450 mediated epoxidation of epothilone C/D to A/B by CYP167A1 (EpoK) – needs significant improvement, in particular regarding the efficiency of its redox partners. Therefore, we have investigated the ability of various hetero- and homologous redox partners to transfer electrons to EpoK. Hereby, a new hybrid system was established with conversion rates eleven times higher and V_max of more than seven orders of magnitudes higher as compared with the previously described spinach redox chain. This hybrid system is the most efficient redox chain for EpoK described to date. Furthermore, P450s from So ce56 were identified which are able to convert epothilone D to 14-OH, 21-OH, 26-OH epothilone D and 7-ketone epothilone D. The latter one represents a novel epothilone derivative and is a suitable candidate for pharmacological tests. The results revealed myxobacterial P450s from S. cellulosum So ce56 as promising candidates for protein engineering for biotechnological production of epothilone derivatives.