Allosteric activation of cytochrome P450 3A4 by α-naphthoflavone: branch point regulation revealed by isotope dilution analysis

Unravelling drug-drug interactions in pigs: Induction of hepatic cytochrome P450 1A (CYP1A) metabolism after the in-feed medication with the anthelmintic fenbendazole

The anthelmintic fenbendazole (FBZ) undergoes hepatic S‑oxygenation by monooxygenases belonging to the cytochrome P450 (CYP) and flavin-monooxygenase (FMO) families. The in-feed medication with FBZ induced CYP1A-dependent metabolism in pig liver. This fact may alter the metabolism of the anthelmintic itself, and of CYP1A substrates like aflatoxin B1 (AFB1). This work evaluated the effect of the in-feed administration of FBZ on CYP1A-dependent metabolism, on its own pattern of hepatic S‑oxygenation, and on the metabolism of AFB1. Landrace piglets remained untreated (n = 5) or received a pre-mix of FBZ (n = 6) in feed for 9 days. Pigs were slaughtered for preparation of liver microsomes used for: CYP content determination; monitoring the CYP1A-dependent enzyme activities, 7-ethoxyresorufin O-deethylase (EROD) and 7-methoxyresorufin O-demethylase (MROD); measurement of FBZ (50 μM) S‑oxygenation, and AFB1 (16 nM) disappearance from the incubation medium. In microsomes of FBZ-treated animals, EROD and MROD increased 19-fold (p = 0.002) and 14-fold (p = 0.003), respectively. An enhanced (3-fold, p = 0.004) participation of the CYP pathway in FBZ S‑oxygenation was observed in the liver of piglets treated with the anthelmintic (210 ± 69 pmol/min.nmol CYP) compared to untreated animals (68 ± 34 pmol/min.nmol CYP). AFB1 metabolism was 93% higher (p = 0.009) in the liver of FBZ-treated compared to untreated pigs. Positive and significant (p < 0.05) correlations were observed between CYP1A-dependent enzyme activities and FBZ or AFB1 metabolism. The sustained administration of FBZ caused an auto-induction of the CYP1A-dependent S‑oxygenation of this anthelmintic. The CYP1A induction triggered by the anthelmintic could amplify the production of AFB1 metabolites in pig liver, including the hepatotoxic AFB1-derived epoxide.+.

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.

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.

Mechanisms of allosteric and mixed mode aromatase inhibitors

Aromatase (CYP19) catalyzes the last biosynthetic step of estrogens in mammals and is a primary drug target for hormone-related breast cancer. However, treatment with aromatase inhibitors is often associated with adverse effects and drug resistance. In this study, we used virtual screening targeting a predicted cytochrome P450 reductase binding site on aromatase to discover four novel non-steroidal aromatase inhibitors. The inhibitors have potencies comparable to the noncompetitive tamoxifen metabolite, endoxifen. Our two most potent inhibitors, AR11 and AR13, exhibit both mixed-type and competitive-type inhibition. The cytochrome P450 reductase-CYP19 coupling interface likely acts as a transient binding site. Our modeling shows that our inhibitors bind better at different sites near the catalytic site. Our results predict the location of multiple ligand binding sites on aromatase. The combination of modeling and experimental results supports the important role of the reductase binding interface as a low affinity, promiscuous ligand binding site. Our new inhibitors may be useful as alternative chemical scaffolds that may show different adverse effects profiles than current clinically used aromatase inhibitors.

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.

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.

Dynamics and Location of the Allosteric Midazolam Site in Cytochrome P4503A4 in Lipid Nanodiscs

Promiscuous and allosteric drug interactions with cytochrome P450 3A4 (CYP3A4) are ubiquitous but incompletely understood at the molecular level. A classic allosteric CYP3A4 drug interaction includes the benzodiazepine midazolam (MDZ). MDZ exhibits homotropic and heterotropic allostery when metabolized to 1'-hydroxy and 4-hydroxy metabolites in varying ratios. The combination of hydrogen-deuterium exchange mass spectrometry (HDX-MS) and Gaussian accelerated molecular dynamics (GaMD) simulations of CYP3A4 in lipid nanodiscs and in a lipid bilayer, respectively, reveals MDZ-dependent changes in dynamics in a membrane environment. The F-, G-, and intervening helices, as well as the loop preceding the β1-sheets, display the largest observed changes in HDX. The GaMD suggests a potential allosteric binding site for MDZ in the F'- and G'-regions, which undergo significant increases in HDX at near-saturating MDZ concentrations. The HDX-MS and GaMD results confirm that changes in dynamics are most significant near the developing consensus allosteric site, and these changes are distinct from those observed previously with the nonallosteric inhibitor ketoconazole. The results suggest that the allosteric MDZ remains mobile in its binding site at the Phe-cluster. The results further suggest that this binding site remains dynamic or changes the depth of insertion in the membrane.

Steroid bioconjugation to a CYP3A4 allosteric site and its effect on substrate binding and coupling efficiency

Human cytochrome P450 3A4 (CYP3A4) is an important drug metabolizing enzyme involved in a number of drug-drug and food-drug interactions. As such, much effort has been devoted into investigating its mechanism of interaction with ligands. CYP3A4 has one of the highest levels of substrate promiscuity for an enzyme, and can even bind multiple ligands simultaneously. The location and orientation of these ligands depend on the chemical structure and stoichiometry, and are generally poorly understood. In the case of the steroid testosterone, up to three copies of the molecule can associate with the enzyme at once, likely two in the active site and one at a postulated allosteric site. Recently, we demonstrated that steroid bioconjugation at the allosteric site results in an increase in activity of CYP3A4 toward testosterone and 7-benzyloxy-4-trifluoromethylcoumarin oxidation. Here, using the established bioconjugation methodology, we show how steroid bioconjugation at the allosteric site affects the heme spin state, the binding affinity (K_S) of CYP3A4 for testosterone, as well as the enzyme coupling efficiency.

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.

CYP26C1 Is a Hydroxylase of Multiple Active Retinoids and Interacts with Cellular Retinoic Acid Binding Proteins

The clearance of retinoic acid (RA) and its metabolites is believed to be regulated by the CYP26 enzymes, but the specific roles of CYP26A1, CYP26B1, and CYP26C1 in clearing active vitamin A metabolites have not been defined. The goal of this study was to establish the substrate specificity of CYP26C1, and determine whether CYP26C1 interacts with cellular retinoic acid binding proteins (CRABPs). CYP26C1 was found to effectively metabolize all-trans retinoic acid (atRA), 9-cis-retinoic acid (9-cis-RA), 13-cis-retinoic acid, and 4-oxo-atRA with the highest intrinsic clearance toward 9-cis-RA. In comparison with CYP26A1 and CYP26B1, CYP26C1 resulted in a different metabolite profile for retinoids, suggesting differences in the active-site structure of CYP26C1 compared with other CYP26s. Homology modeling of CYP26C1 suggested that this is attributable to the distinct binding orientation of retinoids within the CYP26C1 active site. In comparison with other CYP26 family members, CYP26C1 was up to 10-fold more efficient in clearing 4-oxo-atRA (intrinsic clearance 153 μl/min/pmol) than CYP26A1 and CYP26B1, suggesting that CYP26C1 may be important in clearing this active retinoid. In support of this, CRABPs delivered 4-oxo-atRA and atRA for metabolism by CYP26C1. Despite the tight binding of 4-oxo-atRA and atRA with CRABPs, the apparent Michaelis-Menten constant in biological matrix (K_m) value of these substrates with CYP26C1 was not increased when the substrates were bound with CRABPs, in contrast to what is predicted by free drug hypothesis. Together these findings suggest that CYP26C1 is a 4-oxo-atRA hydroxylase and may be important in regulating the concentrations of this active retinoid in human tissues.

Membrane-embedded substrate recognition by cytochrome P450 3A4

Cytochrome P450 3A4 (CYP3A4) is the dominant xenobiotic-metabolizing enzyme in the liver and intestine and is involved in the disposition of more than 50% of drugs. Because of its ability to bind multiple substrates, its reaction kinetics are complex, and its association with the microsomal membrane confounds our understanding of how this enzyme recognizes and recruits diverse substrates. Testosterone (TST) hydroxylation is the prototypical CYP3A4 reaction, displaying positive homotropic cooperativity with three binding sites. Here, exploiting the capability of accelerated molecular dynamics (aMD) to sample events in the millisecond regime, I performed >25-μs aMD simulations in the presence of three TST molecules. These simulations identified high-occupancy surface-binding sites as well as a pathway for TST ingress into the CYP3A4 active site originating in the membrane. Adaptive biasing force analysis of the latter pathway revealed a metastable intermediate that could constitute a third binding site at high TST concentrations. Prompted by the observation that interactions between TST and the G'-helix mobilize the ligand into the active site, a free-energy analysis of TST distribution in the membrane was conducted and revealed that the depth of the G'-helix is optimal for extracting TST. In summary, these simulations confirm separate, but adjacent substrate-binding sites within the enzyme and the existence of an auxiliary TST-binding site. The broader impact of these simulations is that they support a mechanism in which cytochromes P450 directly recruit membrane-solubilized substrates.

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.

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.

7,8-benzoflavone binding to human cytochrome P450 3A4 reveals complex fluorescence quenching, suggesting binding at multiple protein sites

Human cytochrome P450 (P450) 3A4 is involved in the metabolism of one-half of marketed drugs and shows cooperative interactions with some substrates and other ligands. The interaction between P450 3A4 and the known allosteric effector 7,8-benzoflavone (α-naphthoflavone, αNF) was characterized using steady-state fluorescence spectroscopy. The binding interaction of P450 3A4 and αNF effectively quenched the fluorescence of both the enzyme and ligand. The Hill Equation and Stern-Volmer fluorescence quenching models were used to evaluate binding of ligand to enzyme. P450 3A4 fluorescence was quenched by titration with αNF; at the relatively higher [αNF]/[P450 3A4] ratios in this experiment, two weaker quenching interactions were revealed (K_d 1.8-2.5 and 6.5 μM). A range is given for the stronger interaction since αNF quenching of P450 3A4 fluorescence changed the protein spectral profile: quenching of 315 nm emission was slightly more efficient (K_d 1.8 μM) than the quenching of protein fluorescence at 335 and 355 nm (K_d 2.5 and 2.1 μM, respectively). In the reverse titration, αNF fluorescence was quenched by P450 3A4; at the lower [αNF]/[P450 3A4] ratios here, two strong quenching interactions were revealed (K_d 0.048 and 1.0 μM). Thus, four binding interactions of αNF to P450 3A4 are suggested by this study, one of which may be newly recognized and which could affect studies of drug oxidations by this important enzyme.

Direct protein-protein interactions and substrate channeling between cellular retinoic acid binding proteins and CYP26B1

Cellular retinoic acid binding proteins (CRABPs) bind all-trans-retinoic acid (atRA) tightly. This study aimed to determine whether atRA is channeled directly to cytochrome P450 (CYP) CYP26B1 by CRABPs, and whether CRABPs interact directly with CYP26B1. atRA bound to CRABPs (holo-CRABP) was efficiently metabolized by CYP26B1. Isotope dilution experiments showed that delivery of atRA to CYP26B1 in solution was similar with or without CRABP. Holo-CRABPs had higher affinity for CYP26B1 than free atRA, but both apo-CRABPs inhibited the formation of 4-OH-RA by CYP26B1. Similar protein-protein interactions between soluble binding proteins and CYPs may be important for other lipophilic CYP substrates.

Membrane Interactions, Ligand-Dependent Dynamics, and Stability of Cytochrome P4503A4 in Lipid Nanodiscs

Membrane-bound cytochrome P4503A4 (CYP3A4) is the major source of enzymatic drug metabolism. Although several structural models of CYP3A4 in various ligand complexes are available, none includes a lipid bilayer. Details of the effects of the membrane on protein dynamics and solvation, and access channels for ligands, remain uncertain. H/D exchange mass spectrometry (H/DXMS) with ligand free CYP3A4 containing a deletion of residues 3-12, compared to that of the full length wild type, in lipid nanodiscs afforded 91% sequence coverage. Deuterium exchange was fast in the F- and G-helices, HI loop, and C-terminal loop. In contrast, there is very low exchange in the F'- and G'-helices. The results are consistent with the overall membrane orientation of CYP3A4 suggested by published MD simulations and spectroscopic results, and the solvent accessibility of the F/G loop suggests that it is not deeply membrane-embedded. Addition of ketoconazole results in only modest, but global, changes in solvent accessibility. Interestingly, with ketoconazole bound some peptides become less solvent accessible or dynamic, including the F- and G-helices, but several peptides demonstrate modestly increased accessibility. Differential scanning calorimetry (DSC) of CYP3A4-nanodiscs suggests membrane-induced stabilization compared to that of aggregated CYP3A4 in buffer, and this stabilization is enhanced upon addition of the ligand ketoconazole. This ligand-induced stabilization is accompanied by a very large increase in ΔH for CYP3A4 denaturation in nanodiscs, possibly due to increased CYP3A4-membrane interactions. Together, the results suggest a distinct orientation of CYP3A4 on the lipid membrane, and they highlight likely solvent access channels, which are consistent with several MD simulations.

Individualization of treatments with drugs metabolized by CES1: combining genetics and metabolomics

CES1 is involved in the hydrolysis of ester group-containing xenobiotic and endobiotic compounds including several essential and commonly used drugs. The individual variation in the efficacy and tolerability of many drugs metabolized by CES1 is considerable. Hence, there is a large interest in individualizing the treatment with these drugs. The present review addresses the issue of individualized treatment with drugs metabolized by CES1. It describes the composition of the gene encoding CES1, reports variants of this gene with focus upon those with a potential effect on drug metabolism and provides an overview of the protein structure of this enzyme bringing notice to mechanisms involved in the regulation of enzyme activity. Subsequently, the review highlights drugs metabolized by CES1 and argues that individual differences in the pharmacokinetics of these drugs play an important role in determining drug response and tolerability suggesting prospects for individualized drug therapies. Our review also discusses endogenous substrates of CES1 and assesses the potential of using metabolomic profiling of blood to identify proxies for the hepatic activity of CES1 that predict the rate of drug metabolism. Finally, the combination of genetics and metabolomics to obtain an accurate prediction of the individual response to CES1-dependent drugs is discussed.

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

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

Drug modulation of water-heme interactions in low-spin P450 complexes of CYP2C9d and CYP125A1

Azoles and pyridines are commonly incorporated into small molecule inhibitor scaffolds that target cytochromes P450 (CYPs) as a strategy to increase drug binding affinity, impart isoform-dependent selectivity, and improve metabolic stability. Optical absorbance spectra of the CYP-inhibitor complex are widely used to infer whether these inhibitors are ligated directly to the heme iron as catalytically inert, low-spin (type II) complexes. Here, we show that the low-spin complex between a drug-metabolizing CYP2C9 variant and 4-(3-phenylpropyl)-1H-1,2,3-triazole (PPT) retains an axial water ligand despite exhibiting elements of "classic" type II optical behavior. Hydrogens of the axial water ligand are observed by pulsed electron paramagnetic resonance (EPR) spectroscopy for both inhibitor-free and inhibitor-bound species and show that inhibitor binding does not displace the axial water. A (15)N label incorporated into PPT is 0.444 nm from the heme iron, showing that PPT is also in the active site. The reverse type I inhibitor, LP10, of CYP125A1 from Mycobacterium tuberculosis, known from X-ray crystal structures to form a low-spin water-bridged complex, is found by EPR and by visible and near-infrared magnetic circular dichroism spectroscopy to retain the axial water ligand in the complex in solution.

Current Approaches for Investigating and Predicting Cytochrome P450 3A4-Ligand Interactions

Cytochrome P450 3A4 (CYP3A4) is the major and most important drug-metabolizing enzyme in humans that oxidizes and clears over a half of all administered pharmaceuticals. This is possible because CYP3A4 is promiscuous with respect to substrate binding and has the ability to catalyze diverse oxidative chemistries in addition to traditional hydroxylation reactions. Furthermore, CYP3A4 binds and oxidizes a number of substrates in a cooperative manner and can be both induced and inactivated by drugs. In vivo, CYP3A4 inhibition could lead to undesired drug-drug interactions and drug toxicity, a major reason for late-stage clinical failures and withdrawal of marketed pharmaceuticals. Owing to its central role in drug metabolism, many aspects of CYP3A4 catalysis have been extensively studied by various techniques. Here, we give an overview of experimental and theoretical methods currently used for investigation and prediction of CYP3A4-ligand interactions, a defining factor in drug metabolism, with an emphasis on the problems addressed and conclusions derived from the studies.

Strength of axial water ligation in substrate-free cytochrome P450s is isoform dependent

The heme-containing cytochrome P450s exhibit isoform-dependent ferric spin equilibria in the resting state and differential substrate-dependent spin equilibria. The basis for these differences is not well understood. Here, magnetic circular dichroism (MCD) reveals significant differences in the resting low spin ligand field of CYPs 3A4, 2E1, 2C9, 125A1, and 51B1, which indicates differences in the strength of axial water ligation to the heme. The near-infrared bands that specifically correspond to charge-transfer porphyrin-to-metal transitions span a range of energies of nearly 2 kcal/mol. In addition, the experimentally determined MCD bands are not entirely in agreement with the expected MCD energies calculated from electron paramagnetic resonance parameters, thus emphasizing the need for the experimental data. MCD marker bands of the high spin heme between 500 and 680 nm were also measured and suggest only a narrow range of energies for this ensemble of high spin Cys(S(-)) → Fe(3+) transitions among these isoforms. The differences in axial ligand energies between CYP isoforms of the low spin states likely contribute to the energetics of substrate-dependent spin state perturbation. However, these ligand field energies do not correlate with the fraction of high spin vs low spin in the resting state enzyme, suggestive of differences in water access to the heme or isoform-dependent differences in the substrate-free high spin states as well.

Allosteric modulation of substrate motion in cytochrome P450 3A4-mediated xylene oxidation

Many cytochrome P450 enzymes (CYPs) exhibit allosteric behavior reflecting a complex ligand-binding process involving numerous factors: conformational selection, protein-protein interactions, substrate/effector/protein structure, and multiple-ligand binding. The interplay of CYP plasticity and rigidity contributes to substrate/product selectivity and to allosterism. Detailed evidence describing how protein motion modulates product selectivity is incomplete as are descriptions of effector-induced modulation of substrate dynamics. Our intent was to discover details of allosteric behavior and CYP3A4 flexibility and rigidity by investigating substrate motion using low-molecular weight ligands. Steady state kinetics and product ratios were measured for oxidation of m-xylene-(2)H3 and p-xylene; intramolecular isotope effects were measured for m-xylene-(2)H3 oxidation as a function of m-xylene-(2)H3 and p-xylene concentration. Biphasic kinetic plots indicated homotropic cooperative behavior with xylene isomers. Selectivity for aromatic hydroxylation over benzylic hydroxylation of m-xylene-(2)H3 supports a model in which the region near the CYP3A4 active oxidizing species limits substrate dynamics. p-Xylene impedes the motion of m-xylene-(2)H3 substrates that have access to the active oxidizing species: (kH/kD)obs values for m-xylene-(2)H3 decreased with p-xylene concentration. m-Xylene-(2)H3 and p-xylene do not have simultaneous access to the active oxidizing species: deuterium-labeled and unlabeled p-xylene exhibited similar effects on the (kH/kD)obs values for m-xylene-(2)H3 oxidation. p-Xylene and m-xylene-(2)H3 bind at different sites: m-xylene-(2)H3 oxidation rates and product selectivity were consistent across the p-xylene concentration range. Overall, this study indicates that the intramolecular isotope effect experimental design provides a unique opportunity to investigate allosteric mechanisms as it provides information about substrate motion when the enzyme is primed to oxidize substrates.

Identifying cytochrome p450 functional networks and their allosteric regulatory elements

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

Cytochrome P450-dependent catabolism of vitamin K: ω-hydroxylation catalyzed by human CYP4F2 and CYP4F11

Vitamin K plays an essential role in many biological processes including blood clotting, maintenance of bone health, and inhibition of arterial calcification. A menaquinone form of vitamin K, MK4, is increasingly recognized for its key roles in mitochondrial electron transport, as a ligand for the nuclear receptor SXR, which controls the expression of genes involved in transport and metabolism of endo- and xenobiotics, and as a pharmacotherapeutic in the treatment of osteoporosis. Although cytochrome P450 (CYP) 4F2 activity is recognized as an important determinant of phylloquinone (K1) metabolism, the enzymes involved in menaquinone catabolism have not been studied previously. CYP4F2 and CYP4F11 were expressed and purified and found to be equally efficient as in vitro catalysts of MK4 ω-hydroxylation. CYP4F2, but not CYP4F11, catalyzed sequential metabolism of MK4 to the ω-acid without apparent release of the intermediate aldehyde. The ω-alcohol could also be metabolized to the acid by microsomal NAD(+)-dependent alcohol and aldehyde dehydrogenases. LC-MS/MS analysis of trypsinized human liver microsomes (using a surrogate peptide approach) revealed the mean concentrations of CYP4F2 and CYP4F11 to be 14.3 and 8.4 pmol/mg protein, respectively. Microsomal MK4 ω-hydroxylation activities correlated with the CYP4F2 V433M genotype but not the CYP4F11 D446N genotype. Collectively, these data expand the lexicon of vitamin K ω-hydroxylases to include the 'orphan' P450 CYP4F11 and identify a common variant, CYP4F2 (rs2108622), as a major pharmacogenetic variable influencing MK4 catabolism.

Cooperative effects for CYP2E1 differ between styrene and its metabolites

Cooperative interactions are frequently observed in the metabolism of drugs and pollutants by cytochrome P450s; nevertheless, the molecular determinants for cooperativity remain elusive. Previously, we demonstrated that steady-state styrene metabolism by CYP2E1 exhibits positive cooperativity. We hypothesized that styrene metabolites have lower affinity than styrene toward CYP2E1 and limited ability to induce cooperative effects during metabolism. To test the hypothesis, we determined the potency and mechanism of inhibition for styrene and its metabolites toward oxidation of 4-nitrophenol using CYP2E1 Supersomes® and human liver microsomes. Styrene inhibited the reaction through a mixed cooperative mechanism with high affinity for the catalytic site (67 µM) and lower affinity for the cooperative site (1100 µM), while increasing substrate turnover at high concentrations. Styrene oxide and 4-vinylphenol possessed similar affinity for CYP2E1. Styrene oxide behaved cooperatively like styrene, but 4-vinylphenol decreased turnover at high concentrations. Styrene glycol was a very poor competitive inhibitor. Among all compounds, there was a positive correlation with binding and hydrophobicity. Taken together, these findings for CYP2E1 further validate contributions of cooperative mechanisms to metabolic processes, demonstrate the role of molecular structure on those mechanisms and underscore the potential for heterotropic cooperative effects between different compounds.

The biochemistry of acetaminophen hepatotoxicity and rescue: a mathematical model

BACKGROUND

Acetaminophen (N-acetyl-para-aminophenol) is the most widely used over-the-counter or prescription painkiller in the world. Acetaminophen is metabolized in the liver where a toxic byproduct is produced that can be removed by conjugation with glutathione. Acetaminophen overdoses, either accidental or intentional, are the leading cause of acute liver failure in the United States, accounting for 56,000 emergency room visits per year. The standard treatment for overdose is N-acetyl-cysteine (NAC), which is given to stimulate the production of glutathione.

METHODS

We have created a mathematical model for acetaminophen transport and metabolism including the following compartments: gut, plasma, liver, tissue, urine. In the liver compartment the metabolism of acetaminophen includes sulfation, glucoronidation, conjugation with glutathione, production of the toxic metabolite, and liver damage, taking biochemical parameters from the literature whenever possible. This model is then connected to a previously constructed model of glutathione metabolism.

RESULTS

We show that our model accurately reproduces published clinical and experimental data on the dose-dependent time course of acetaminophen in the plasma, the accumulation of acetaminophen and its metabolites in the urine, and the depletion of glutathione caused by conjugation with the toxic product. We use the model to study the extent of liver damage caused by overdoses or by chronic use of therapeutic doses, and the effects of polymorphisms in glucoronidation enzymes. We use the model to study the depletion of glutathione and the effect of the size and timing of N-acetyl-cysteine doses given as an antidote. Our model accurately predicts patient death or recovery depending on size of APAP overdose and time of treatment.

CONCLUSIONS

The mathematical model provides a new tool for studying the effects of various doses of acetaminophen on the liver metabolism of acetaminophen and glutathione. It can be used to study how the metabolism of acetaminophen depends on the expression level of liver enzymes. Finally, it can be used to predict patient metabolic and physiological responses to APAP doses and different NAC dosing strategies.

Optimization of the bacterial cytochrome P450 BM3 system for the production of human drug metabolites

Drug metabolism in human liver is a process involving many different enzymes. Among them, a number of cytochromes P450 isoforms catalyze the oxidation of most of the drugs commercially available. Each P450 isoform acts on more than one drug, and one drug may be oxidized by more than one enzyme. As a result, multiple products may be obtained from the same drug, and as the metabolites can be biologically active and may cause adverse drug reactions (ADRs), the metabolic profile of a new drug has to be known before this can be commercialized. Therefore, the metabolites of a certain drug must be identified, synthesized and tested for toxicity. Their synthesis must be in sufficient quantities to be used for metabolic tests. This review focuses on the progresses done in the field of the optimization of a bacterial self-sufficient and efficient cytochrome P450, P450 BM3 from Bacillus megaterium, used for the production of metabolites of human enzymes. The progress made in the improvement of its catalytic performance towards drugs, the substitution of the costly NADPH cofactor and its immobilization and scale-up of the process for industrial application are reported.

Cross-linking mass spectrometry and mutagenesis confirm the functional importance of surface interactions between CYP3A4 and holo/apo cytochrome b(5)

Cytochrome b(5) (cyt b(5)) is one of the key components in the microsomal cytochrome P450 monooxygenase system. Consensus has not been reached about the underlying mechanism of cyt b(5) modulation of CYP catalysis. Both cyt b(5) and apo b(5) are reported to stimulate the activity of several P450 isoforms. In this study, the surface interactions of both holo and apo b(5) with CYP3A4 were investigated and compared for the first time. Chemical cross-linking coupled with mass spectrometric analysis was used to identify the potential electrostatic interactions between the protein surfaces. Subsequently, the models of interaction of holo/apo b(5) with CYP3A4 were built using the identified interacting sites as constraints. Both cyt b(5) and apo b(5) were predicted to bind to the same groove on CYP3A4 with close contacts to the B-B' loop of CYP3A4, a substrate recognition site. Mutagenesis studies further confirmed that the interacting sites on CYP3A4 (Lys96, Lys127, and Lys421) are functionally important. Mutation of these residues reduced or abolished cyt b(5) binding affinity. The critical role of Arg446 on CYP3A4 in binding to cyt b(5) and/or cytochrome P450 reductase was also discovered. The results indicated that electrostatic interactions on the interface of the two proteins are functionally important. The results indicate that apo b(5) can dock with CYP3A4 in a manner analogous to that of holo b(5), so electron transfer from cyt b(5) is not required for its effects.

Understanding the mechanism of cytochrome P450 3A4: recent advances and remaining problems

Cytochromes P450 (CYPs) represent a diverse group of heme-thiolate proteins found in almost all organisms. CYPs share a common protein fold but differ in substrate selectivity and catalyze a wide variety of monooxygenation reactions via activation of molecular oxygen. Among 57 human P450s, the 3A4 isoform (CYP3A4) is the most abundant and the most important because it metabolizes the majority of administered drugs. A remarkable feature of CYP3A4 is its extreme promiscuity in substrate specificity and cooperative substrate binding, which often leads to undesirable drug-drug interactions and toxic side effects. Owing to its importance in drug development and therapy, CYP3A4 has been the most extensively studied mammalian P450. In this review we provide an overview on recent progress and remaining problems in the CYP3A4 research.

1,2,3-Triazole-heme interactions in cytochrome P450: functionally competent triazole-water-heme complexes

In comparison to imidazole (IMZ) and 1,2,4-triazole (1,2,4-TRZ), the isosteric 1,2,3-triazole (1,2,3-TRZ) is unrepresented among cytochrome P450 (CYP) inhibitors. This is surprising because 1,2,3-TRZs are easily obtained via "click" chemistry. To understand this underrepresentation of 1,2,3-TRZs among CYP inhibitors, thermodynamic and density functional theory computational studies were performed with unsubstituted IMZ, 1,2,4-TRZ, and 1,2,3-TRZ. The results indicate that the lower affinity of 1,2,3-TRZ for the heme iron includes a large unfavorable entropy term likely originating in solvent-1,2,3-TRZ interactions; the difference is not solely due to differences in the enthalpy of heme-ligand interactions. In addition, the 1,2,3-TRZ fragment was incorporated into a well-established CYP3A4 substrate and mechanism-based inactivator, 17-α-ethynylestradiol (17EE), via click chemistry. This derivative, 17-click, yielded optical spectra consistent with low-spin ferric heme iron (type II) in contrast to 17EE, which yields a high-spin complex (type I). Furthermore, the rate of CYP3A4-mediated metabolism of 17-click was comparable to that of 17EE, with a different regioselectivity. Surprisingly, continuous-wave electron paramagnetic resonance (EPR) and HYSCORE EPR spectroscopy indicate that 17-click does not displace water from the sixth axial ligand position of CYP3A4 as expected for a type II ligand. We propose a binding model in which 17-click pendant 1,2,3-TRZ hydrogen bonds with the sixth axial water ligand. The results demonstrate the potential for 1,2,3-TRZ to form metabolically labile water-bridged low-spin heme complexes, consistent with recent evidence that nitrogenous type II ligands of CYPs can be efficiently metabolized. The specific case of [CYP3A4·17-click] highlights the risk of interpreting CYP-ligand complex structure on the basis of optical spectra.

A novel type of allosteric regulation: functional cooperativity in monomeric proteins

Cooperative functional properties and allosteric regulation in cytochromes P450 play an important role in xenobiotic metabolism and define one of the main mechanisms of drug-drug interactions. Recent experimental results suggest that ability to bind simultaneously two or more small organic molecules can be the essential feature of cytochrome P450 fold, and often results in rich and complex pattern of allosteric behavior. Manifestations of non-Michaelis kinetics include homotropic and heterotropic activation and inhibition effects depending on the stoichiometric ratios of substrate and effector, changes in the regio- and stereospecificity of catalytic transformations, and often give rise to the clinically important drug-drug interactions. In addition, functional response of P450 systems is modulated by the presence of specific and non-specific effector molecules, metal ions, membrane incorporation, formation of homo- and hetero-oligomers, and interactions with the protein redox partners. In this article we briefly overview the main factors contributing to the allosteric effects in cytochromes P450 with the main focus on the sources of cooperative behavior in xenobiotic metabolizing monomeric heme enzymes with their conformational flexibility and extremely broad substrate specificity. The novel mechanism of functional cooperativity in P450 enzymes does not require substantial binding cooperativity, rather it implies the presence of one or more binding sites with higher affinity than the single catalytically active site in the vicinity of the heme iron.