Modelling atypical CYP3A4 kinetics: principles and pragmatism

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.

Atypical kinetics of cytochrome P450 enzymes in pharmacology and toxicology

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

INTERSYSTEM EXTRAPOLATION FACTORS (ISEF) ARE SUBSTRATE-DEPENDENT FOR CYP3A4: IMPACT ON CYTOCHROME P450 REACTION PHENOTYPING

The use of intersystem extrapolation factors (ISEF) is required for the quantitative scaling of drug metabolism data generated in individually expressed cytochrome P450 enzymes when estimating fractional contribution to metabolism by P450 enzymes in vivo (f_m,_CYP). For successful prediction of f_m, ISEF values must be universal across all substrates for any individual enzyme. In this study, ISEF values were generated for ten CYP3A4 selective substrates using a common source of recombinant heterologously expressed CYP3A4 and a pool of human liver microsomes. The resulting ISEF values for CYP3A4 were substrate-dependent and ranged 8-fold, with the highest value generated from intrinsic clearance of midazolam depletion (0.36) and the lowest from quinidine depletion (0.044). Application of these ISEF values for estimation of the fractional contribution of CYP3A4 and CYP2C19 to omeprazole clearance yielded values that ranged from 0.21-0.63 and 0.37-0.79, respectively, as compared to back-extrapolated in vivo f_m values of 0.27 (CYP3A4) and 0.85 (CYP2C19) from clinical pharmacokinetic data. For risperidone, estimated f_m values for CYP3A4 and CYP2D6 ranged from 0.87-0.98 and 0.02-0.13, respectively, as compared to in vivo values of 0.36 (CYP3A4) and 0.63-0.88 (CYP2D6), showing that the importance of CYP3A4 was over-estimated and the importance of CYP2D6 under-estimated. Overall, these findings suggest that ISEF values for CYP3A4 can vary with the marker substrate used to derive them, thereby reducing the effectiveness of the approach of using metabolism data from rCYP3A4 with ISEF values for the prediction of f_m values in vivo. Significance Statement Intersystem extrapolation factors (ISEF) are utilized for assigning fractional contributions of individual enzymes to drug clearance (f_m) from drug metabolism data generated in recombinant P450s. The present data shows that ISEF values for cytochrome P4503A4 vary with the substrate. This can lead to variable and erroneous prediction of f_m.

Fundamentals of Enzyme Kinetics: Michaelis-Menten and Non-Michaelis-Type (Atypical) Enzyme Kinetics

This chapter will provide a general introduction to the kinetics of enzyme-catalyzed reactions, including a general discussion of catalysts, reaction rates, and binding constants. This section will be followed by a discussion of various types of enzyme kinetics observed in drug metabolism reactions. A large number of enzymatic reactions can be adequately described by Michaelis-Menten kinetics. The Michaelis-Menten equation represents a rectangular hyperbola, with a y-asymptote at the V_max value. However, in other cases, more complex kinetic models are required to explain the observed data. Atypical kinetic profiles are believed to arise from the simultaneous binding of multiple molecules within the active site of the enzyme (Tracy and Hummel, Drug Metab Rev 36:231-242, 2004). Several cytochromes P450 (CYPs) have large active sites that enable binding of multiple molecules (Yano et al., J Biol Chem 279:38091-38094, 2004; Wester et al., J Biol Chem 279:35630-35637, 2004). Thus, atypical kinetics are not uncommon in in vitro drug metabolism studies.

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.

Differential Reversible and Irreversible Interactions between Benzbromarone and Human Cytochrome P450s 3A4 and 3A5

Mounting evidence has revealed that despite the high degree of sequence homology between cytochrome P450 3A isoforms (i.e., CYP3A4 and CYP3A5), they have the propensities to exhibit vastly different irreversible and reversible interactions with a single substrate. We have previously established that benzbromarone (BBR), a potent uricosuric agent used in the management of gout, irreversibly inhibits CYP3A4 via mechanism-based inactivation (MBI). However, it remains unelucidated if CYP3A5-its highly homologous counterpart-is susceptible to inactivation by BBR. Using three structurally distinct probe substrates, we consistently demonstrated that MBI was not elicited in CYP3A5 by BBR. Our in silico covalent docking models and molecular dynamics simulations suggested that disparities in the susceptibilities toward MBI could be attributed to the specific effects of BBR covalent adducts on the F-F' loop. Serendipitously, we also discovered that BBR reversibly activated CYP3A5-mediated rivaroxaban hydroxylation wherein apparent V _max increased and K _m decreased with increasing BBR concentration. Fitting data to the two-site model yielded interaction factors α and β of 0.44 and 5.88, respectively, thereby confirming heterotropic activation of CYP3A5 by BBR. Furthermore, heteroactivation was suppressed by the CYP3A inhibitor ketoconazole in a concentration-dependent manner and decreased with increasing preincubation time, implying that activation was incited via binding of parent BBR molecule within the enzymatic active site. Finally, noncovalent docking revealed that CYP3A5 can more favorably accommodate both BBR and rivaroxaban in concert as compared with CYP3A4, which further substantiated our experimental observations. SIGNIFICANCE STATEMENT: Although it has been previously demonstrated that benzbromarone (BBR) inactivates CYP3A4, it remains uninterrogated whether it also elicits mechanism-based inactivation in CYP3A5, which shares ∼85% sequence similarity with CYP3A4. This study reported that BBR exhibited differential irreversible and reversible interactions with both CYP3A isoforms and further unraveled the molecular determinants underpinning their diverging interactions. These data offer important insight into differential kinetic behavior of CYP3A4 and CYP3A5, which potentially contributes to interindividual variabilities in drug disposition.

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.

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

Background and Purpose

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

Experimental Approach

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

Key Results

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

Conclusions and Implications

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

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

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

Metabolism of megestrol acetate in vitro and the role of oxidative metabolites

1. There is limited knowledge regarding the metabolism of megestrol acetate (MA), as it was approved by FDA in 1971, prior to the availability of modern tools for identifying specific drug-metabolizing enzymes. We determined the cytochrome P450s (P450s) and UDP-glucuronosyltransferases (UGTs) that metabolize MA, identified oxidative metabolites and determined pharmacologic activity at the progesterone, androgen and glucocorticoid receptors (PR, AR and GR, respectively). 2. Oxidative metabolites were produced using human liver microsomes (HLMs), and isolated for mass spectral (MS) and nuclear magnetic resonance (NMR) analyses. We screened recombinant P450s using MA at 62 μM (HLM K_m for metabolite 1; M1) and 28 μM (HLM K_m for metabolite 2; M2). UGT isoforms were simultaneously incubated with UDPGA, nicotinamide adenine dinucleotide phosphate (NADPH), CYP3A4 and MA. Metabolites were evaluated for pharmacologic activity on the PR, AR and GR. CYP3A4 and CYP3A5 are responsible for oxidative metabolism of 62 μM MA. 3. At 28 μM substrate concentration, CYP3A4 was the only contributing enzyme. Mass spectral and NMR data suggest metabolism of MA to two alcohols. After oxidation, MA is converted into two secondary glucuronides by UGT2B17 among other UGTs. MA, M1 and M2 had significant pharmacologic activity on the PR while only MA showed activity on the AR and GR.

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.

Physical Studies of P450-P450 Interactions: Predicting Quaternary Structures of P450 Complexes in Membranes from Their X-ray Crystal Structures

Cytochrome P450 enzymes, which catalyze oxygenation reactions of both exogenous and endogenous chemicals, are membrane bound proteins that require interaction with their redox partners in order to function. Those responsible for drug and foreign compound metabolism are localized primarily in the endoplasmic reticulum of liver, lung, intestine, and other tissues. More recently, the potential for P450 enzymes to exist as supramolecular complexes has been shown by the demonstration of both homomeric and heteromeric complexes. The P450 units in these complexes are heterogeneous with respect to their distribution and function, and the interaction of different P450s can influence P450-specific metabolism. The goal of this review is to examine the evidence supporting the existence of physical complexes among P450 enzymes. Additionally, the review examines the crystal lattices of different P450 enzymes derived from X-ray diffraction data to make assumptions regarding possible quaternary structures in membranes and in turn, to predict how the quaternary structures could influence metabolism and explain the functional effects of specific P450-P450 interactions.

In vitro oxidative decarboxylation of free fatty acids to terminal alkenes by two new P450 peroxygenases

BACKGROUND

P450 fatty acid decarboxylases represented by the unusual CYP152 peroxygenase family member OleT_JE have been receiving great attention recently since these P450 enzymes are able to catalyze the simple and direct production of 1-alkenes for potential applications in biofuels and biomaterials. To gain more mechanistic insights, broader substrate spectra, and improved decarboxylative activities, it is demanded to discover and investigate more P450 fatty acid decarboxylases.

RESULTS

Here, we describe for the first time the expression, purification, and in vitro biochemical characterization of two new CYP152 peroxygenases, CYP-Aa162 and CYP-Sm46Δ29, that are capable of decarboxylating straight-chain saturated fatty acids. Both enzymes were found to catalyze the decarboxylation and hydroxylation of a broad range of free fatty acids (C₁₀-C₂₀) with overlapping substrate specificity, yet distinct chemoselectivity. CYP-Sm46Δ29 works primarily as a fatty (lauric) acid decarboxylase (66.1 ± 3.9% 1-undecene production) while CYP-Aa162 more as a fatty (lauric) acid hydroxylase (72.2 ± 0.9% hydroxy lauric acid production). Notably, the optical spectroscopic analysis of functional CYP-Sm46Δ29 revealed no characteristic P450 band, suggesting a unique heme coordination environment. Active-site mutagenesis analysis showed that substitution with the proposed key decarboxylation-modulating residues, His85 and Ile170, enhanced the decarboxylation activity of CYP-Aa162 and P450_BSβ, emphasizing the importance of these residues in directing the decarboxylation pathway. Furthermore, the steady-state kinetic analysis of CYP-Aa162 and CYP-Sm46Δ29 revealed both cooperative and substrate inhibition behaviors which are substrate carbon chain length dependent.

CONCLUSIONS

Our data identify CYP-Sm46Δ29 as an efficient OleT_JE-like fatty acid decarboxylase. Oxidative decarboxylation chemoselectivity of the CYP152 decarboxylases is largely dependent upon the carbon chain length of fatty acid substrates and their precise positioning in the enzyme active site. Finally, the kinetic mode analysis of the enzymes could provide important guidance for future process design.

Key to Opening Kidney for In Vitro-In Vivo Extrapolation Entrance in Health and Disease: Part I: In Vitro Systems and Physiological Data

The programme for the 2015 AAPS Annual Meeting and Exhibition (Orlando, FL; 25-29 October 2015) included a sunrise session presenting an overview of the state-of-the-art tools for in vitro-in vivo extrapolation (IVIVE) and mechanistic prediction of renal drug disposition. These concepts are based on approaches developed for prediction of hepatic clearance, with consideration of scaling factors physiologically relevant to kidney and the unique and complex structural organisation of this organ. Physiologically relevant kidney models require a number of parameters for mechanistic description of processes, supported by quantitative information on renal physiology (system parameters) and in vitro/in silico drug-related data. This review expands upon the themes raised during the session and highlights the importance of high quality in vitro drug data generated in appropriate experimental setup and robust system-related information for successful IVIVE of renal drug disposition. The different in vitro systems available for studying renal drug metabolism and transport are summarised and recent developments involving state-of-the-art technologies highlighted. Current gaps and uncertainties associated with system parameters related to human kidney for the development of physiologically based pharmacokinetic (PBPK) model and quantitative prediction of renal drug disposition, excretion, and/or metabolism are identified.

Impact of CYP2C8*3 polymorphism on in vitro metabolism of imatinib to N-desmethyl imatinib

1. Imatinib is metabolized to N-desmethyl imatinib by CYPs 3A4 and 2C8. The effect of CYP2C8*3 genotype on N-desmethyl imatinib formation was unknown. 2. We examined imatinib N-demethylation in human liver microsomes (HLMs) genotyped for CYP2C8*3, in CYP2C8*3/*3 pooled HLMs and in recombinant CYP2C8 and CYP3A4 enzymes. Effects of CYP-selective inhibitors on N-demethylation were also determined. 3. A single-enzyme Michaelis-Menten model with autoinhibition best fitted CYP2C8*1/*1 HLM (n = 5) and recombinant CYP2C8 kinetic data (median ± SD Ki = 139 ± 61 µM and 149 µM, respectively). Recombinant CYP3A4 showed two-site enzyme kinetics with no autoinhibition. Three of four CYP2C8*1/*3 HLMs showed single-enzyme kinetics with no autoinhibition. Binding affinity was higher in CYP2C8*1/*3 than CYP2C8*1/*1 HLM (median ± SD Km = 6 ± 2 versus 11 ± 2 µM, P=0.04). CYP2C8*3/*3 (pooled HLM) also showed high binding affinity (Km = 4 µM) and single-enzyme weak autoinhibition (Ki = 449 µM) kinetics. CYP2C8 inhibitors reduced HLM N-demethylation by 47-75%, compared to 0-30% for CYP3A4 inhibitors. 4. In conclusion, CYP2C8*3 is a gain-of-function polymorphism for imatinib N-demethylation, which appears to be mainly mediated by CYP2C8 and not CYP3A4 in vitro in HLM.

In vitro metabolism of the alkaloid piplartine by rat liver microsomes

Because piplartine (PPT) has demonstrated biological activities, such as cytotoxic, anxiolytic, antidepressant, antifungal and antiplatelet activities, this molecule is a relevant drug candidate. The metabolic fate of drug candidates is an essential requirement in assessing their safety and efficacy. Based on this requirement, the biotransformation of PPT by cytochrome P450 enzymes (CYP) was investigated for the first time. To determine the in vitro enzymatic kinetic parameters, an HPLC method was developed and validated to quantify PPT. All samples were separated on a reversed-phase C18 column using a mobile phase of acetonitrile:water (40:60, v/v). The method exhibited a linear range of 2.4-157.7 μmol/L, with the following calibration curve: y=0.0934 (±0.0010)x+0.0027, r=0.9975. The lower limit of quantitation was verified to be 2.4 μmol/L, with an RSD below 7%. The precision and accuracy were assessed for both within-day and between-day determinations; neither relative standard (RSD%) deviations nor relative errors (RER) exceeded a value of 15%. The mean absolute recovery was 85%, with an RSD value below 6%. The enzymatic kinetic parameters revealed a sigmoidal profile, with V(max)=4.7±0.3 μmol/mg mL⁻¹/min, h=2.5±0.4, S₅₀=44.7±0.3 μmol/L and CL(max)=0.054 μL/min/mg protein, indicating cooperativity behavior. Employing a mammalian model, PPT metabolism yielded two unreported monohydroxylated products (m/z 334). The identification and structural elucidation of the metabolites were performed by comparing their mass spectra with those spectra of the parent drug. For the first time, the in vitro metabolism studies employing microsomes were demonstrated to be a suitable tool for data regarding enzymatic kinetics and for the metabolites formed in the PPT mammalian metabolism.

Different enzyme kinetic models

As described in Chapter 2 , a large number of enzymatic reactions can be adequately described by Michaelis-Menten kinetics. The Michaelis-Menten equation represents a rectangular hyperbola, with a y-asymptote at the V max value. In many cases, more complex kinetic models are required to explain the observed data. Atypical kinetic profiles are believed to arise from the simultaneous binding of multiple molecules within the active site of the enzyme (Tracy and Hummel, Drug Metab Rev 36:231-242, 2004). Several cytochromes P450 have large active sites that enable binding of multiple molecules (Wester et al. J Biol Chem 279:35630-35637, 2004; Yano et al. J Biol Chem 279:38091-38094, 2004). Thus, atypical kinetics are not uncommon in in vitro drug metabolism studies. This chapter covers enzyme kinetic reactions in which a single enzyme has multiple binding sites for substrates and/or inhibitors as well as reactions catalyzed by multiple enzymes.

Pivotal role of P450-P450 interactions in CYP3A4 allostery: the case of α-naphthoflavone

We investigated the relationship between oligomerization of CYP3A4 (cytochrome P450 3A4) and its response to ANF (α-naphthoflavone), a prototypical heterotropic activator. The addition of ANF resulted in over a 2-fold increase in the rate of CYP3A4-dependent debenzylation of 7-BFC [7-benzyloxy-4-(trifluoromethyl)coumarin] in HLM (human liver microsomes), but failed to produce activation in BD Supersomes or Baculosomes containing recombinant CYP3A4 and NADPH-CPR (cytochrome P450 reductase). However, incorporation of purified CYP3A4 into Supersomes containing only recombinant CPR reproduced the behaviour observed with HLM. The activation in this system was dependent on the surface density of the enzyme. Although no activation was detectable at an L/P (lipid/P450) ratio ≥750, it reached 225% at an L/P ratio of 140. To explore the relationship between this effect and CYP3A4 oligomerization, we probed P450-P450 interactions with a new technique that employs LRET (luminescence resonance energy transfer). The amplitude of LRET in mixed oligomers of the haem protein labelled with donor and acceptor fluorophores exhibited a sigmoidal dependence on the surface density of CYP3A4 in Supersomes™. The addition of ANF eliminated this sigmoidal character and increased the degree of oligomerization at low enzyme concentrations. Therefore the mechanisms of CYP3A4 allostery with ANF involve effector-dependent modulation of P450-P450 interactions.

Electrochemical detection of anti-breast-cancer agents in human serum by cytochrome P450-coated carbon nanotubes

We report on the electrochemical detection of anti-cancer drugs in human serum with sensitivity values in the range of 8–925 nA/μM. Multi-walled carbon nanotubes were functionalized with three different cytochrome P450 isoforms (CYP1A2, CYP2B6, and CYP3A4). A model used to effectively describe the cytochrome P450 deposition onto carbon nanotubes was confirmed by Monte Carlo simulations. Voltammetric measurements were performed in phosphate buffer saline (PBS) as well as in human serum, giving well-defined current responses upon addition of increasing concentrations of anti-cancer drugs. The results assert the capability to measure concentration of drugs in the pharmacological ranges in human serum. Another important result is the possibility to detect pairs of drugs present in the same sample, which is highly required in case of therapies with high side-effects risk and in anti-cancer pharmacological treatments based on mixtures of different drugs. Our technology holds potentials for inexpensive multi-panel drug-monitoring in personalized therapy.

Modeling chemical interaction profiles: II. Molecular docking, spectral data-activity relationship, and structure-activity relationship models for potent and weak inhibitors of cytochrome P450 CYP3A4 isozyme

Polypharmacy increasingly has become a topic of public health concern, particularly as the U.S. population ages. Drug labels often contain insufficient information to enable the clinician to safely use multiple drugs. Because many of the drugs are bio-transformed by cytochrome P450 (CYP) enzymes, inhibition of CYP activity has long been associated with potentially adverse health effects. In an attempt to reduce the uncertainty pertaining to CYP-mediated drug-drug/chemical interactions, an interagency collaborative group developed a consensus approach to prioritizing information concerning CYP inhibition. The consensus involved computational molecular docking, spectral data-activity relationship (SDAR), and structure-activity relationship (SAR) models that addressed the clinical potency of CYP inhibition. The models were built upon chemicals that were categorized as either potent or weak inhibitors of the CYP3A4 isozyme. The categorization was carried out using information from clinical trials because currently available in vitro high-throughput screening data were not fully representative of the in vivo potency of inhibition. During categorization it was found that compounds, which break the Lipinski rule of five by molecular weight, were about twice more likely to be inhibitors of CYP3A4 compared to those, which obey the rule. Similarly, among inhibitors that break the rule, potent inhibitors were 2–3 times more frequent. The molecular docking classification relied on logistic regression, by which the docking scores from different docking algorithms, CYP3A4 three-dimensional structures, and binding sites on them were combined in a unified probabilistic model. The SDAR models employed a multiple linear regression approach applied to binned 1D ¹³C-NMR and 1D ¹⁵N-NMR spectral descriptors. Structure-based and physical-chemical descriptors were used as the basis for developing SAR models by the decision forest method. Thirty-three potent inhibitors and 88 weak inhibitors of CYP3A4 were used to train the models. Using these models, a synthetic majority rules consensus classifier was implemented, while the confidence of estimation was assigned following the percent agreement strategy. The classifier was applied to a testing set of 120 inhibitors not included in the development of the models. Five compounds of the test set, including known strong inhibitors dalfopristin and tioconazole, were classified as probable potent inhibitors of CYP3A4. Other known strong inhibitors, such as lopinavir, oltipraz, quercetin, raloxifene, and troglitazone, were among 18 compounds classified as plausible potent inhibitors of CYP3A4. The consensus estimation of inhibition potency is expected to aid in the nomination of pharmaceuticals, dietary supplements, environmental pollutants, and occupational and other chemicals for in-depth evaluation of the CYP3A4 inhibitory activity. It may serve also as an estimate of chemical interactions via CYP3A4 metabolic pharmacokinetic pathways occurring through polypharmacy and nutritional and environmental exposures to chemical mixtures.

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.

Substrate inhibition kinetics in drug metabolism reactions

Inhibition of enzyme activity at high substrate concentrations, so-called "substrate inhibition," is commonly observed and has been recognized in drug metabolism reactions since the last decade. Although the importance of such "atypical" kinetics in vivo remains poorly understood, a substrate with substrate inhibition kinetics has been shown to unconventionally alter the metabolism of other substrates. In recent years, it is becoming increasingly evident that the mechanisms for substrate inhibition are highly complex, which are possibly contributed by multiple (at least two) binding sites within the enzyme protein, the formation of a ternary dead-end enzyme complex, and/or the ligand-induced changes in enzyme conformation. This review primarily discusses the mechanisms for substrate inhibition displayed by the important drug-metabolizing enzymes, such as cytochrome p450s, UDP-glucuronyltransferases, and sulfotransferases. Kinetic modeling of substrate inhibition in the absence or presence of a modifier is another central issue in this review because of its importance in the determination of kinetic parameters and in vitro/in vivo predictions.

Application of homology modeling to generate CYP1A1 mutants with enhanced activation of the cancer chemotherapeutic prodrug dacarbazine

The chemotherapeutic prodrug dacarbazine (DTIC) has limited efficacy in human malignancies and exhibits numerous adverse effects that arise from systemic exposure to the cytotoxic metabolite. DTIC is activated by CYP1A1 and CYP1A2 catalyzed N-demethylation. However, structural features of these enzymes that confer DTIC N-demethylation have not been characterized. A validated homology model of CYP1A1 was employed to elucidate structure-activity relationships and to engineer CYP1A1 enzymes with altered DTIC activation. In silico docking demonstrated that DTIC orientates proximally to Ser122, Phe123, Asp313, Ala317, Ile386, Tyr259, and Leu496 of human CYP1A1. The site of metabolism is positioned 5.6 Å from the heme iron at an angle of 105.3°. Binding in the active site is stabilized by H-bonding between Tyr259 and the N(2) position of the imidazole ring. Twenty-seven CYP1A1 mutants were generated and expressed in Escherichia coli in yields ranging from 9 to 225 pmol P450/mg. DTIC N-demethylation by the E161K, E256K, and I458V mutants exhibited Michaelis-Menten kinetics, with decreases in K(m) (183-249 μM) that doubled the catalytic efficiency (p < 0.05) relative to wild-type CYP1A1 (K(m), 408 ± 43 μM; V(max), 28 ± 4 pmol · min(-1) · pmol of P450(-1)). The generation of enzymes with catalytically enhanced DTIC activation highlights the potential use of mutant CYP1A1 proteins in P450-based gene-directed enzyme prodrug therapy for the treatment of metastatic malignant melanoma.

Mechanism of cytochrome P450-3A inhibition by ketoconazole

OBJECTIVES

Ketoconazole is extensively used as an index inhibitor of cytochrome P450-3A (CYP3A) activity in vitro and in vivo, but the mechanism of ketoconazole inhibition of CYP3A still is not clearly established.

METHODS

Inhibition of metabolite formation by ketoconazole (seven concentrations from 0.01 to 1.0 µm) was studied in human liver microsomes (n = 4) at six to seven substrate concentrations for triazolam, midazolam, and testosterone, and at two substrate concentrations for nifedipine.

KEY FINDINGS

Analysis of multiple data points per liver sample based on a mixed competitive-noncompetitive model yielded mean inhibition constant K(i) values in the range of 0.011 to 0.045 µm. Ketoconazole IC50 increased at higher substrate concentrations, thereby excluding pure noncompetitive inhibition. For triazolam, testosterone, and midazolam α-hydroxylation, mean values of α (indicating the 'mix' of competitive and noncompetitive inhibition) ranged from 2.1 to 6.3. However, inhibition of midazolam 4-hydroxylation was consistent with a competitive process. Determination of K(i) and α based on the relation between 50% inhibitory concentration values and substrate concentration yielded similar values. Pre-incubation of ketoconazole with microsomes before addition of substrate did not enhance inhibition, whereas inhibition by troleandomycin was significantly enhanced by pre-incubation.

CONCLUSIONS

Ketoconazole inhibition of triazolam α- and 4-hydroxylation, midazolam α-hydroxylation, testosterone 6β-hydroxylation, and nifedipine oxidation appeared to be a mixed competitive-noncompetitive process, with the noncompetitive component being dominant but not exclusive. Quantitative estimates of K(i) were in the low nanomolar range for all four substrates.

First-pass metabolism via UDP-glucuronosyltransferase: a barrier to oral bioavailability of phenolics

Glucuronidation mediated by UDP-glucuronosyltransferases (UGTs) is a significant metabolic pathway that facilitates efficient elimination of numerous endobiotics and xenobiotics, including phenolics. UGT genetic deficiency and polymorphisms or inhibition of glucuronidation by concomitant use of drugs are associated with inherited physiological disorders or drug-induced toxicities. Moreover, extensive glucuronidation can be a barrier to oral bioavailability as the first-pass glucuronidation (or premature clearance by UGTs) of orally administered agents usually results in the poor oral bioavailability and lack of efficacies. This review focused on the first-pass glucuronidation of phenolics including natural polyphenols and pharmaceuticals. The complexity of UGT-mediated metabolism of phenolics is highlighted with species-, gender-, organ- and isoform-dependent specificity, as well as functional compensation between UGT1A and 2B subfamily. In addition, recent advances are discussed with respect to the mechanisms of enzymatic actions, including the important properties such as binding pocket size and phosphorylation requirements.

Correlation between bilirubin glucuronidation and estradiol-3-gluronidation in the presence of model UDP-glucuronosyltransferase 1A1 substrates/inhibitors

Inhibition of UDP-glucuronosyltransferase (UGT) 1A1-catalyzed bilirubin glucuronidation by drug compounds may potentially be of clinical concern. However, in drug discovery and development settings, bilirubin is less than an ideal in vitro probe for assessing the potential of a chemical entity to inhibit bilirubin glucuronidation. In part, this is due to the propensity of bilirubin to photodegrade and to the instability of its metabolites. To this end, the utility of estradiol-3-glucuronidation as a surrogate in vitro predictor for interactions with bilirubin was evaluated. The glucuronidation kinetics of bilirubin and estradiol were carefully characterized with recombinant UGT1A1 expressed in human embryonic kidney 293 cells. Consistent with previous reports, estradiol-3-glucuronidation displayed sigmoidal kinetics, whereas bilirubin glucuronidation exhibited typical hyperbolic kinetics. The two compounds also mutually inhibited the metabolism of the other. Sixteen UGT1A1 substrates/inhibitors were evaluated as effectors of each reaction. Fourteen compounds inhibited both bilirubin and estradiol glucuronidation. However, two compounds (ethinylestradiol and daidzein) exhibited mixed effects (concentration-dependent activation and inhibition) on estradiol-3-glucuronidation, whereas bilirubin glucuronidation was inhibited by both compounds. In addition, 7-ethyl-10-hydroxycamptothecin, a substrate of UGT1A1 (reported K(m) = 24 μM) seemed to be a weak inhibitor of bilirubin glucuronidation (IC(50) = 356.4 μM) but a partial inhibitor of estradiol-3-glucuronidation. The IC(50) values of the inhibitors against estradiol-3-glucuronidation were strongly correlated with IC(50) values against bilirubin glucuronidation, resulting in an R(2) value of 0.9604 (activator excluded) or 0.8287 (activator included). Thus, estradiol-3-glucuronidation can serve as a good surrogate for predicting inhibition of bilirubin glucuronidation with the caveat that occasionally compounds may demonstrate activation of estradiol-3-glucuronidation.

Inhibition of aromatase activity in MCF-7aro human breast cancer cells by the natural androgens testosterone and androstenedione

: The human breast contains all the enzymes responsible for local bioformation of estradiol (E

: The cells were incubated with physiological concentrations of [

: The MCF-7aro cells have a very high aromatase activity because conversion of [

: In breast cancer cells, the natural androgens testosterone and androstenedione, have the capacity to control bioformation of estradiol by blocking aromatase activity. The data can provide important information on the control mechanism of estrogen intratumoral levels and open new possibilities in breast cancer treatment.

Prediction of pharmacokinetic drug-drug interaction caused by changes in cytochrome P450 activity using in vivo information

The aim of the present paper was to present an overview of the current status of the methods used to predict the magnitude of pharmacokinetic drug-drug interactions (DDIs) which are caused by apparent changes in cytochrome P450 (CYP) activity with an emphasis on a method using in vivo information. In addition, more than a hundred representative CYP substrates, inhibitor and inducer drugs involved in significant pharmacokinetic DDIs were selected from the literature and are listed. Although the magnitude of DDIs has been conventionally predicted based on in vitro experiments, their predictability is restricted occasionally due to several difficulties, including a precise determination of the unbound inhibitor concentrations at the enzyme site and a reliable in vitro measurement of the inhibition constant (K(i)). Alternatively, a simple method has been recently proposed for the prediction of the magnitude of DDIs based on information fully available from in vivo clinical studies. The new in vivo-based method would be applicable to the adjustment of dose regimens in actual pharmacotherapy situations although it requires a prior clinical study for the prediction. In this review, theoretical and quantitative relationships between the in vivo- and the in vitro-based prediction methods are considered. One of the interesting outcomes of the consideration is that the K(i)-normalized dose (dose/in vitro K(i)) of larger than approximately 20L (2-200L, when variability is considered) may be a pragmatic index which predicts significant in vivo DDIs. In the last part of the article, the relevance of the inclusion of the in vivo-based method into the process of new drug development is discussed for good prediction of in vivo DDIs.

Cooperative properties of cytochromes P450

Cytochromes P450 form a large and important class of heme monooxygenases with a broad spectrum of substrates and corresponding functions, from steroid hormone biosynthesis to the metabolism of xenobiotics. Despite decades of study, the molecular mechanisms responsible for the complex non-Michaelis behavior observed with many members of this superfamily during metabolism, often termed 'cooperativity', remain to be fully elucidated. Although there is evidence that oligomerization may play an important role in defining the observed cooperativity, some monomeric cytochromes P450, particularly those involved in xenobiotic metabolism, also display this behavior due to their ability to simultaneously bind several substrate molecules. As a result, formation of distinct enzyme-substrate complexes with different stoichiometry and functional properties can give rise to homotropic and heterotropic cooperative behavior. This review aims to summarize the current understanding of cooperativity in cytochromes P450, with a focus on the nature of cooperative effects in monomeric enzymes.

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

According to the initial hypothesis on the mechanisms of cooperativity in drug-metabolizing cytochromes P450, a loose fit of a single substrate molecule in the P450 active site results in a requirement for the binding of multiple ligand molecules for efficient catalysis. Although simultaneous occupancy of the active site by multiple ligands is now well established, there is increasing evidence that the mechanistic basis of cooperativity also involves an important ligand-induced conformational transition. Moreover, recent studies demonstrate that the conformational heterogeneity of the enzyme is stabilized by ligand-dependent interactions of several P450 molecules. Application of the concept of an oligomeric allosteric enzyme to microsomal cytochromes P450 in combination with a general paradigm of multiple ligand occupancy of the active site provides an excellent explanation for complex manifestations of the atypical kinetic behavior of the enzyme.

Curcuminoid biosynthesis by two type III polyketide synthases in the herb Curcuma longa

Curcuminoids found in the rhizome of turmeric, Curcuma longa, possess various biological activities. Despite much attention regarding the biosynthesis of curcuminoids because of their pharmaceutically important properties and biosynthetically intriguing structures, no enzyme systems have been elucidated. Here we propose a pathway for curcuminoid biosynthesis in the herb C. longa, which includes two novel type III polyketide synthases. One of the type III polyketide synthases, named diketide-CoA synthase (DCS), catalyzed the formation of feruloyldiketide-CoA by condensing feruloyl-CoA and malonyl-CoA. The other, named curcumin synthase (CURS), catalyzed the in vitro formation of curcuminoids from cinnamoyldiketide-N-acetylcysteamine (a mimic of the CoA ester) and feruloyl-CoA. Co-incubation of DCS and CURS in the presence of feruloyl-CoA and malonyl-CoA yielded curcumin at high efficiency, although CURS itself possessed low activity for the synthesis of curcumin from feruloyl-CoA and malonyl-CoA. These findings thus revealed the curcumin biosynthetic route in turmeric, in which DCS synthesizes feruloyldiketide-CoA, and CURS then converts the diketide-CoA esters into a curcuminoid scaffold.

Potential role of intestinal first-pass metabolism in the prediction of drug-drug interactions

BACKGROUND

The contribution of intestine to the magnitude of drug-drug interactions (DDI) may be significant, considering high levels of inhibitors in the gut lumen achieved during absorption and the abundance of metabolic enzymes in the mature enterocytes. Intestinal inhibition is incorporated in the DDI prediction models as the ratio of the intestinal wall availability in the presence and absence of the inhibitor (F(G)(') and F(G), respectively).

OBJECTIVE

This review will focus on the ability of the current approaches to estimate the extent of intestinal DDI accurately, addressing predominantly the most abundant intestinal P450 enzyme, CYP3A4.

METHODS

Considering the sensitivity of the DDI prediction models to the accuracy of the F(G) estimates, the current study focuses on 3 different in vitro and in vivo approaches to assess this parameter.

RESULTS/CONCLUSION

The advantages and limitations of each of F(G) methods are outlined. Accurate assessment of this parameter is essential for the prediction of human drug clearance and drug-drug interaction potential.

Enzyme kinetics of GTI-2040, a phosphorothioate oligonucleotide targeting ribonucleotide reductase

Enzyme kinetics of GTI-2040 (5'-GGC TAA ATC GCT CCA CCA AG-3'), a phosphorothioate ribonucleotide reductase antisense, were investigated for the first time in 3' exonuclease solution and human liver microsomes (HLMs), using the ion-pair high-performance liquid chromatogram method for quantification of the parent drug and two major 3'N-1 and 3'N-2 metabolites. Enzyme kinetics of GTI-2040 in 3'-exonuclease solution were found to be well characterized by the Michaelis-Menten model, using the sum of formation rates of 3'N-1 and 3'N-2 (approximately total metabolism) because of sequential metabolism. In HLMs, a biphasic binding was observed for GTI-2040 with high- and low-affinity constants (K(d)s) of 0.03 and 3.8 microM, respectively. Enzyme kinetics of GTI-2040 in HLMs were found to deviate from Michaelis-Menten kinetics when the total GTI-2040 substrate was used. However, after correction for the unbound fractions, the formation rate of total metabolites could be described by Michaelis-Menten kinetics. Using the free substrate fraction, the K(m) and V(max) of GTI-2040 were determined to be 6.33 +/- 3.2 microM and 16.5 +/- 8.4 nmol/mg/h, respectively. Using these values, in vitro hepatic intrinsic clearance (CL(int)) in HLM was estimated to be 2.61 +/- 0.56 ml/h. The CL(int) was then used to predict GTI-2040's in vivo intrinsic clearance in humans by a microsomal protein scaling factor, which gave a mean value of 182.7 l/h, representing 24.1% of the observed in vivo mean scaled hepatic intrinsic clearance of 758.7 l/h in patients with acute myeloid leukemia. We concluded that the saturable nonspecific binding of GTI-2040 in HLMs complicated the interpretation of its enzyme kinetics, and scaled intrinsic clearance from HLMs only partially predicted the in vivo intrinsic clearance.

Multiple substrate binding by cytochrome P450 3A4: estimation of the number of bound substrate molecules

Cytochrome P450 3A4, a major drug-metabolizing enzyme in man, is well known to show non-Michaelis-Menten steady-state kinetics for a number of substrates, indicating that more than one substrate can bind to the enzyme simultaneously, but it has proved difficult to obtain reliable estimates of exactly how many substrate molecules can bind. We have used a simple method involving studies of the effect of large inhibitors on the Hill coefficient to provide improved estimates of substrate stoichiometry from simple steady-state kinetics. Using a panel of eight inhibitors, we show that at least four molecules of the widely used CYP3A4 substrate 7-benzyloxyquinoline can bind simultaneously to the enzyme. Computational docking studies show that this is consistent with the recently reported crystal structures of the enzyme. In the case of midazolam, which shows simple Michaelis-Menten kinetics, the inhibitor effects demonstrate that two molecules must bind simultaneously, consistent with earlier evidence, whereas for diltiazem, the experiments provide no evidence for the binding of more than one molecule. The consequences of this "inhibitor-induced cooperativity" for the prediction of pharmacokinetics and drug-drug interactions are discussed.

Inhibition and induction of human cytochrome P450 enzymes: current status

Variability of drug metabolism, especially that of the most important phase I enzymes or cytochrome P450 (CYP) enzymes, is an important complicating factor in many areas of pharmacology and toxicology, in drug development, preclinical toxicity studies, clinical trials, drug therapy, environmental exposures and risk assessment. These frequently enormous consequences in mind, predictive and pre-emptying measures have been a top priority in both pharmacology and toxicology. This means the development of predictive in vitro approaches. The sound prediction is always based on the firm background of basic research on the phenomena of inhibition and induction and their underlying mechanisms; consequently the description of these aspects is the purpose of this review. We cover both inhibition and induction of CYP enzymes, always keeping in mind the basic mechanisms on which to build predictive and preventive in vitro approaches. Just because validation is an essential part of any in vitro-in vivo extrapolation scenario, we cover also necessary in vivo research and findings in order to provide a proper view to justify in vitro approaches and observations.

In vitro evaluation of reversible and irreversible cytochrome P450 inhibition: current status on methodologies and their utility for predicting drug-drug interactions

It is widely accepted that today's practice of polypharmacy inevitably increases the incidence of drug-drug interactions (DDIs). Serious DDI is a major liability for any new chemical entity (NCE) entering the pharmaceutical market. As such, pharmaceutical companies employ various strategies to avoid problematic compounds for clinical development. A key cause for DDIs is the inhibition of cytochrome P450 enzymes (CYPs) that are responsible for metabolic clearance of many drugs. Screening for inhibition potency of CYPs by NCEs has therefore become a routine practice during the drug discovery stage. However, in order to make proper use of DDI data, an understanding of the strengths and weaknesses of the various experimental systems in current use is required. An illustrated review of experimental practices is presented with discussion of likely future developments. The combination of high quality in vitro data generation and the application of in vivo CYP inhibition modelling approaches should allow more informed decisions to be made in the search for drug molecules with acceptable DDI characteristics.

Modeling the catalysis of anti-cocaine catalytic antibody: competing reaction pathways and free energy barriers

The competing reaction pathways and the corresponding free energy barriers for cocaine hydrolysis catalyzed by an anti-cocaine catalytic antibody, mAb15A10, were studied by using a novel computational strategy based on the binding free energy calculations on the antibody binding with cocaine and transition states. The calculated binding free energies were used to evaluate the free energy barrier shift from the cocaine hydrolysis in water to the antibody-catalyzed cocaine hydrolysis for each reaction pathway. The free energy barriers for the antibody-catalyzed cocaine hydrolysis were predicted to be the corresponding free energy barriers for the cocaine hydrolysis in water plus the calculated free energy barrier shifts. The calculated free energy barrier shift of -6.87 kcal/mol from the dominant reaction pathway of the cocaine benzoyl ester hydrolysis in water to the dominant reaction pathway of the antibody-catalyzed cocaine hydrolysis is in good agreement with the experimentally derived free energy barrier shift of -5.93 kcal/mol. The calculated mutation-caused shifts of the free energy barrier are also reasonably close to the available experimental activity data. The good agreement suggests that the protocol for calculating the free energy barrier shift from the cocaine hydrolysis in water to the antibody-catalyzed cocaine hydrolysis may be used in future rational design of possible high-activity mutants of the antibody as anti-cocaine therapeutics. The general strategy of the free energy barrier shift calculation may also be valuable in studying a variety of chemical reactions catalyzed by other antibodies or proteins through noncovalent bonding interactions with the substrates.

Kinetics of electron transfer in the complex of cytochrome P450 3A4 with the flavin domain of cytochrome P450BM-3 as evidence of functional heterogeneity of the heme protein

We used a rapid scanning stop-flow technique to study the kinetics of reduction of cytochrome P450 3A4 (CYP3A4) by the flavin domain of cytochrome P450-BM3 (BMR), which was shown to form a stoichiometric complex (K(D)=0.48 microM) with CYP3A4. In the absence of substrates only about 50% of CYP3A4 was able to accept electrons from BMR. Whereas the high-spin fraction was completely reducible, the reducibility of the low-spin fraction did not exceed 42%. Among four substrates tested (testosterone, 1-pyrenebutanol, bromocriptine, or alpha-naphthoflavone (ANF)) only ANF is capable of increasing the reducibility of the low-spin fraction to 75%. Our results demonstrate that the pool of CYP3A4 is heterogeneous, and not all P450 is competent for electron transfer in the complex with reductase. The increase in the reducibility of the enzyme in the presence of ANF may represent an important element of the mechanism of action of this activator.