Cooperative Binding of Midazolam with Testosterone and α-Naphthoflavone within the CYP3A4 Active Site: A NMR T1 Paramagnetic Relaxation Study

Assessing the role of residue Phe108 of cytochrome P450 3A4 in allosteric effects of midazolam metabolism

Cytochrome P450 3A4 (CYP3A4) is involved in the metabolism of more drugs in clinical use than any other xenobiotic-metabolizing enzyme. CYP3A4-mediated drug metabolism is usually allosterically modulated by substrate concentration (homotropic allostery) and other drugs (heterotropic allostery), exhibiting unusual kinetic profiles and regiospecific metabolism. Recent studies suggest that residue Phe108 (F108) of CYP3A4 may have an important role in drug metabolism. In this work, residue mutations were coupled with well-tempered metadynamics simulations to assess the importance of F108 in the allosteric effects of midazolam metabolism. Comparing the simulation results of the wild-type and mutation systems, we identify that the π-π interaction and steric effect between the F108 side chain and midazolam is favorable for the stable binding of substrate in the active site. F108 also plays an important role in the transition of substrate binding mode, which mainly induces the transition of substrate binding mode by forming π-π interactions with multiple aromatic rings of the substrate. Moreover, the side chain of F108 is closely related to the radius and depth of the 2a and 2f channels, and F108 may further regulate drug metabolism by affecting the pathway, orientation, or time of substrate entry into the CYP3A4 active site or product egress from the active site. Altogether, we suggest that F108 affects drug metabolism and regulatory mechanisms by affecting substrate binding stability, binding mode transition, and channel characteristics of CYP3A4. Our findings could promote the understanding of complicated allosteric mechanisms in CYP3A4-mediated drug metabolism, and the knowledge could be used for drug development and disease treatment.

High-throughput assay to simultaneously evaluate activation of CYP3A and the direct and time-dependent inhibition of CYP3A, CYP2C9, and CYP2D6 using liquid chromatography-tandem mass spectrometry

In the early stages of drug discovery, adequate evaluation of the potential drug-drug interactions (DDIs) of drug candidates is important. Several CYP3A activators are known to lead to underestimation of DDIs. These compounds affect midazolam 1'-hydroxylation but not midazolam 4-hydroxylation.We used both metabolic reactions of midazolam to evaluate the activation and inhibition of CYP3A activators simultaneously. For our CYP inhibition assay using cocktail probe substrates, simultaneous liquid chromatography-tandem mass spectrometry monitoring of 1'-hydroxymidazolam and 4-hydroxymidazolam for CYP3A was established in addition to monitoring of 4-hydroxydiclofenac and 1'-hydroxybufuralol for CYP2C9 and CYP2D6.The results of our cocktail inhibition assay were well correlated with those of a single inhibition assay, as were the estimated inhibition parameters for typical CYP3A inhibitors. In our assay, a proprietary compound that activated midazolam 1'-hydroxylation and tended to inhibit 4-hydroxylation was evaluated along with known CYP3A activators. All compounds were well characterised by comparison of the results of midazolam 1'- and 4-hydroxylation.In conclusion, our CYP cocktail inhibition assay can detect CYP3A activation and assess the direct and time-dependent inhibition potentials for CYP3A, CYP2C9, and CYP2D6. This method is expected to be very efficient in the early stages of drug discovery.

The molecular basis of dapsone activation of CYP2C9-catalyzed nonsteroidal anti-inflammatory drug oxidation

Positive heterotropic cooperativity, or "activation," results in an instantaneous increase in enzyme activity in the absence of an increase in protein expression. Thus, cytochrome P450 (CYP) enzyme activation presents as a potential drug-drug interaction mechanism. It has been demonstrated previously that dapsone activates the CYP2C9-catalyzed oxidation of a number of nonsteroidal anti-inflammatory drugs in vitro. Here, we conducted molecular dynamics simulations (MDS) together with enzyme kinetic investigations and site-directed mutagenesis to elucidate the molecular basis of the activation of CYP2C9-catalyzed S-flurbiprofen 4'-hydroxylation and S-naproxen O-demethylation by dapsone. Supplementation of incubations of recombinant CYP2C9 with dapsone increased the catalytic efficiency of flurbiprofen and naproxen oxidation by 2.3- and 16.5-fold, respectively. MDS demonstrated that activation arises predominantly from aromatic interactions between the substrate, dapsone, and the phenyl rings of Phe114 and Phe476 within a common binding domain of the CYP2C9 active site, rather than involvement of a distinct effector site. Mutagenesis of Phe114 and Phe476 abrogated flurbiprofen and naproxen oxidation, and MDS and kinetic studies with the CYP2C9 mutants further identified a pivotal role of Phe476 in dapsone activation. MDS additionally showed that aromatic stacking interactions between two molecules of naproxen are necessary for binding in a catalytically favorable orientation. In contrast to flurbiprofen and naproxen, dapsone did not activate the 4'-hydroxylation of diclofenac, suggesting that the CYP2C9 active site favors cooperative binding of nonsteroidal anti-inflammatory drugs with a planar or near-planar geometry. More generally, the work confirms the utility of MDS for investigating ligand binding in CYP enzymes.

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

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

Molecular Insights into the Heterotropic Allosteric Mechanism in Cytochrome P450 3A4-Mediated Midazolam Metabolism

Cytochrome P450 3A4 (CYP3A4) is the main P450 enzyme for drug metabolism and drug-drug interactions (DDIs), as it is involved in the metabolic process of approximately 50% of drugs. A detailed mechanistic elucidation of DDIs mediated by CYP3A4 is commonly believed to be critical for drug optimization and rational use. Here, two typical probes, midazolam (MDZ, substrate) and testosterone (TST, allosteric effector), are used to investigate the molecular mechanism of CYP3A4-mediated heterotropic allosteric interactions, through conventional molecular dynamics (cMD) and well-tempered metadynamics (WT-MTD) simulations. Distance monitoring shows that TST can stably bind in two potential peripheral sites (Site 1 and Site 2) of CYP3A4. The binding of TST at these two sites can induce conformational changes in CYP3A4 flexible loops on the basis of conformational analysis, thereby promoting the transition of the MDZ binding mode and affecting the ratio of MDZ metabolites. According to the results of the residue interaction network, multiple allosteric communication pathways are identified that can provide vivid and applicable insights into the heterotropic allostery of TST on MDZ metabolism. Comparing the regulatory effects and the communication pathways, the allosteric effect caused by TST binding in Site 2 seems to be more pronounced than in Site 1. Our findings could provide a deeper understanding of CYP3A4-mediated heterotropic allostery at the atomic level and would be helpful for rational drug use as well as the design of new allosteric modulators.

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.

Investigation of the molecular and mechanistic basis for the regioselective metabolism of midazolam by cytochrome P450 3A4

Cytochrome P450 3A4 (CYP3A4) is the most important P450 enzyme for drug metabolism and drug-drug interaction, due to it being responsible for the biotransformation of approximately 50% of clinically used drugs. Advance knowledge of the molecular and mechanistic basis of CYP3A4 regioselective metabolism is beneficial for understanding the production of metabolites, and may allow personalized metabolic pathways or designing pathway-specific therapeutics. In this work, we focus on investigating the ligand-receptor interactions, substrate conformational transition, and key factors regulating the specificity of metabolic pathways using midazolam (MDZ) as a probe. Here, three types of substrate-binding conformations related to the diversity of MDZ metabolites are identified. The results also suggest that an allosteric site for MDZ is located near the F'-helix, A-anchor, and C-terminal loop of CYP3A4. The presence of an effector in the allosteric site can accelerate the conformational transition of the substrate via modulating a "sandwich" structure, and may affect the proportion of metabolites at high substrate concentration. We hope that the results can improve the understanding of the CYP3A4 structure and function, and provide a new perspective for drug development.

Complex Cytochrome P450 Kinetics Due to Multisubstrate Binding and Sequential Metabolism. Part 2. Modeling of Experimental Data

Three CYP3A4 substrates, midazolam, ticlopidine, and diazepam, display non-Michaelis-Menten kinetics, form multiple primary metabolites, and are sequentially metabolized to secondary metabolites. We generated saturation curves for these compounds and analyzed the resulting datasets using a number of single-substrate and multisubstrate binding models. These models were parameterized using rate equations and numerical solutions of the ordinary differential equations. Multisubstrate binding models provided results superior to single-substrate models, and simultaneous modeling of multiple metabolites provided better results than fitting the individual datasets independently. Although midazolam datasets could be represented using standard two-substrate models, more complex models that include explicit enzyme-product complexes were needed to model the datasets for ticlopidine and diazepam. In vivo clearance predictions improved markedly with the use of in vitro parameters from the complex models versus the Michaelis-Menten equation. The results highlight the need to use sufficiently complex kinetic schemes instead of the Michaelis-Menten equation to generate accurate kinetic parameters. SIGNIFICANCE STATEMENT: The metabolism of midazolam, ticlopidine, and diazepam by CYP3A4 results in multiple metabolites and sequential metabolism. This study evaluates the use of rate equations and numerical methods to characterize the in vitro enzyme kinetics. Use of complex cytochrome P450 kinetic models is necessary to obtain accurate parameter estimates for predicting in vivo disposition.

Enzyme Kinetics of Oxidative Metabolism-Cytochromes P450

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

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

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

Cytochrome P450 3A Enzymes Are Key Contributors for Hepatic Metabolism of Bufotalin, a Natural Constitute in Chinese Medicine Chansu

Bufotalin (BFT), one of the naturally occurring bufodienolides, has multiple pharmacological and toxicological effects including antitumor activity and cardiotoxicity. This study aimed to character the metabolic pathway(s) of BFT and to identify the key drug metabolizing enzyme(s) responsible for hepatic metabolism of BFT in human, as well as to explore the related molecular mechanism of enzymatic selectivity. The major metabolite of BFT in human liver microsomes (HLMs) was fully identified as 5β-hydroxylbufotalin by LC-MS/MS and NMR techniques. Reaction phenotyping and chemical inhibition assays showed that CYP3A4 and CYP3A5 were key enzymes responsible for BFT 5β-hydroxylation. Kinetic analyses demonstrated that BFT 5β-hydroxylation in both HLMs and human CYP3A4 followed the biphasic kinetics, while BFT 5β-hydroxylation in CYP3A5 followed substrate inhibition kinetics. Furthermore, molecular docking simulations showed that BFT could bind on two different ligand-binding sites on both CYP3A4 and CYP3A5, which partially explained the different kinetic behaviors of BFT in CYP3A4 and CYP3A5. These findings are very helpful for elucidating the phase I metabolism of BFT in human and for deeper understanding the key interactions between CYP3A enzymes and bufadienolides, as well as for the development of bufadienolide-type drugs with improved pharmacokinetic and safety profiles.

[Methods for determining of cytochrome P450 isozymes functional activity]

The review is dedicated to modern methods and technologies for determining of cytochrome P450 isozymes functional activity, such as absorbance and fluorescent spectroscopy, electron paramagnetic resonance (EPR), nuclear magnetic resonance (NMR), Raman, Mossbauer, and X-ray spectroscopy, surface plasmon resonance (SPR), atomic force microscopy (AFM). Methods of molecular genetic analysis were reviewed from personalized medicine point of view. The use of chromate-mass-spectrometric methods for cytochrome P450-dependent catalytic reactions' products was discussed. The review covers modern electrochemical systems based on cytochrome P450 isozymes for their catalytic activity analysis, their use in practice and further development perspectives for experimental pharmacology, biotechnology and translational medicine.

Evaluation of modafinil as a perpetrator of metabolic drug-drug interactions using a model informed cocktail reaction phenotyping trial protocol

AIM

To evaluate the capacity for modafinil to be a perpetrator of metabolic drug-drug interactions by altering cytochrome P450 activity following a single dose and dosing to steady state.

METHODS

A single centre, open label, single sequence cocktail drug interaction trial. On days 0, 2 and 8 participants were administered an oral drug cocktail comprising 100 mg caffeine, 30 mg dextromethorphan, 25 mg losartan, 1 mg midazolam and 20 mg enteric-coated omeprazole. Timed blood samples were collected prior to and for up to 6 h post cocktail dosing. Between days 2 and 8 participants orally self-administered 200 mg modafinil each morning.

RESULTS

Following a single 200 mg dose of modafinil mean (± 95% CI) AUC ratios for caffeine, dextromethorphan, losartan, midazolam and omeprazole were 0.95 (± 0.08), 1.01 (± 0.35), 0.97 (± 0.10), 0.98 (± 0.10) and 1.36 (± 0.06), respectively. Following dosing of modafinil to steady state (200 mg for 7 days), AUC ratios for caffeine, dextromethorphan, losartan, midazolam and omeprazole were 0.90 (± 0.16), 0.79 (± 0.09), 0.98 (± 0.11), 0.66 (± 0.12) and 1.90 (± 0.53), respectively.

CONCLUSIONS

These data support consideration of the risk of clinically relevant metabolic drug-drug interactions perpetrated by modafinil when this drug is co-administered with drugs that are primarily cleared by CYP2C19 (single modafinil dose or steady state modafinil dosing) or CYP3A4 (steady state modafinil dosing only) catalysed metabolic pathways.

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

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

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.

Enzyme kinetics of oxidative metabolism: cytochromes P450

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

Paramagnetic nuclear magnetic resonance relaxation and molecular mechanics studies of the chloroperoxidase-indole complex: insights into the mechanism of chloroperoxidase-catalyzed regioselective oxidation of indole

To unravel the mechanism of chloroperoxidase (CPO)-catalyzed regioselective oxidation of indole, we studied the structure of the CPO-indole complex using nuclear magnetic resonance (NMR) relaxation measurements and computational techniques. The dissociation constant (KD) of the CPO-indole complex was calculated to be approximately 21 mM. The distances (r) between protons of indole and the heme iron calculated via NMR relaxation measurements and molecular docking revealed that the pyrrole ring of indole is oriented toward the heme with its 2-H pointing directly at the heme iron. Both KD and r values are independent of pH in the range of 3.0-6.5. The stability and structure of the CPO-indole complex are also independent of the concentration of chloride or iodide ion. Molecular docking suggests the formation of a hydrogen bond between the NH group of indole and the carboxyl O of Glu 183 in the binding of indole to CPO. Simulated annealing of the CPO-indole complex using r values from NMR experiments as distance restraints reveals that the van der Waals interactions were much stronger than the Coulomb interactions in the binding of indole to CPO, indicating that the association of indole with CPO is primarily governed by hydrophobic rather than electrostatic interactions. This work provides the first experimental and theoretical evidence of the long-sought mechanism that leads to the "unexpected" regioselectivity of the CPO-catalyzed oxidation of indole. The structure of the CPO-indole complex will serve as a lighthouse in guiding the design of CPO mutants with tailor-made activities for biotechnological applications.

Jahn-Teller dynamics in a series of high-symmetry Co(II) chelates determine paramagnetic relaxation enhancements

NMR paramagnetic relaxation enhancements (PREs) of a series of structurally characterized, trigonal bis-trispyrazolylborate (Tp) chelates of high-spin Co(II), spanning 100-850 MHz in field, are reported. Prior knowledge of the metal-nucleus distances allows numerical extraction of position-dependent electron spin relaxation rates (τ(c)(-1)) from direct measurement of the individual PREs of the four symmetry distinct protons in Co(Tp)(2), using available closed-form expressions. The data for this electronically complex system where spin-orbit coupling defines the ground state electronic structure are analyzed in terms of the Solomon-Bloembergen-Morgan (SBM) relations, as well as available zero-field splitting limit theories. A simple angular correction is shown to be sufficient to reconcile the individual τ(c)(T) data for the four classes of protons. The data identify a previously unrecognized dynamic Jahn-Teller effect in these historically important complexes, with a barrier of ~230 cm(-1), pointing to a level of dynamics in trispyrazolylborate chemistry that has not been described before, and further show that it is the Jahn-Teller that is responsible for the PREs in fluid solution. A field-dependent component is also identified for the two protons nearest g(//), which is suggested to arise due to Zeeman mixing of excited state character into the ground level.

Preferred binding orientations of phenacetin in CYP1A1 and CYP1A2 are associated with isoform-selective metabolism

Human cytochromes P450 1A1 and 1A2 play important roles in drug metabolism and chemical carcinogenesis. Although these two enzymes share high sequence identity, they display different substrate specificities and inhibitor susceptibilities. In the present studies, we investigated the structural basis for these differences with phenacetin as a probe using a number of complementary approaches, such as enzyme kinetics, stoichiometric assays, NMR, and molecular modeling. Kinetic and stoichiometric analyses revealed that substrate specificity (k(cat)/K(m)) of CYP1A2 was approximately 18-fold greater than that of CYP1A1, as expected. Moreover, despite higher H₂O₂ production, the coupling efficiency of reducing equivalents to acetaminophen formation in CYP1A2 was tighter than that in CYP1A1. CYP1A1, in contrast to CYP1A2, displayed much higher uncoupling, producing more water. The subsequent NMR longitudinal (T₁) relaxation studies with the substrate phenacetin and its product acetaminophen showed that both compounds displayed similar binding orientations within the active site of CYP1A1 and CYP1A2. However, the distance between the OCH₂ protons of the ethoxy group (site of phenacetin O-deethylation) and the heme iron was 1.5 Å shorter in CYP1A2 than in CYP1A1. The NMR findings are thus consistent with our kinetic and stoichiometric results, providing a likely molecular basis for more efficient metabolism of phenacetin by CYP1A2.

Plasticity of CYP2B enzymes: structural and solution biophysical methods

In the past three years, major advances in understanding cytochrome P450 2B (CYP2B) structure-function relationships have been made through determination of multiple ligand-bound and one ligand-free X-ray crystal structure of CYP2B4 and one ligand-bound X-ray crystal structure of CYP2B6. These structures have provided insight into the features that provide the high degree of plasticity of the enzymes. A combination of a phenylalanine cluster that allows for concerted movement of helices F through G and a conserved set of electrostatic interactions involving Arg(262) facilitates movement of this region to accommodate binding of ligands of various sizes without perturbing most of the P450 fold. Integrating solution based techniques such as NMR or deuterium exchange mass spectrometry (DXMS) with computational methods including molecular docking has provided further insight into enzyme behavior upon ligand binding. In addition, extended molecular dynamics simulations have provided a link between an open and a closed conformation of ligand-free CYP2B4 found in crystal structures. Other studies revealed the utility of rational engineering in improving stability of P450s to facilitate structural studies. The solution and computational results combined with the X-ray crystal structures yield a comprehensive picture of how these enzymes adopt different conformations to bind various ligands.

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.

Active site substitution A82W improves the regioselectivity of steroid hydroxylation by cytochrome P450 BM3 mutants as rationalized by spin relaxation nuclear magnetic resonance studies

Cytochrome P450 BM3 from Bacillus megaterium is a monooxygenase with great potential for biotechnological applications. In this paper, we present engineered drug-metabolizing P450 BM3 mutants as a novel tool for regioselective hydroxylation of steroids at position 16β. In particular, we show that by replacing alanine at position 82 with a tryptophan in P450 BM3 mutants M01 and M11, the selectivity toward 16β-hydroxylation for both testosterone and norethisterone was strongly increased. The A82W mutation led to a ≤42-fold increase in V(max) for 16β-hydroxylation of these steroids. Moreover, this mutation improves the coupling efficiency of the enzyme, which might be explained by a more efficient exclusion of water from the active site. The substrate affinity for testosterone increased at least 9-fold in M11 with tryptophan at position 82. A change in the orientation of testosterone in the M11 A82W mutant as compared to the orientation in M11 was observed by T(1) paramagnetic relaxation nuclear magnetic resonance. Testosterone is oriented in M11 with both the A- and D-ring protons closest to the heme iron. Substituting alanine at position 82 with tryptophan results in increased A-ring proton-iron distances, consistent with the relative decrease in the level of A-ring hydroxylation at position 2β.

Binding free energies of inhibitors to iron porphyrin complex as a model for Cytochrome P450

The binding free energies of the inhibitor-heme model complexes are calculated using the density functional methods and the implicit solvation models in water, where the 16 structurally diverse compounds with a spectrum of IC(50) values from 0.05 (clotrimazole) to 1000 (piroxicam) μM are chosen as inhibitors for Cytochrome P450 3A4 (CYP3A4). CYP3A4 is the most predominant constituent of the human hepatic CYP enzymes that play a role in metabolizing structurally diverse xenobiotics. The observed free energy change for each inhibitory binding, ΔG inh0, is obtained from its IC(50) value. The total binding free energy (ΔG b0) of each inhibitor-heme model complex is calculated by the sum of its relative free energy (ΔG(0) ) in the gas phase and solvation free energy to the water-heme model complex. The UB3LYP/LanL2DZ level of theory provides the correct relative stabilities of the high- and low-spin states for the penta- and hexa-coordinated ferric complexes, respectively. The optimized distances of the inhibitor nitrogen (or water oxygen) and the methyl mercaptide S to the ferric iron of the inhibitor-heme model complexes at the same level of theory are consistent with the values of the corresponding X-ray structures, except for the econazole complex. The correlation coefficient r(2) values of 0.91 and 0.75 are obtained from the ΔG b0-ΔG inh0 and ΔG(0) -ΔG inh0 plots, respectively, at the UM06/LanL2DZ:CPCM_UB3LYP/LanL2DZ//UB3LYP/LanL2DZ level of theory in water. This indicates that the total binding free energies calculated for the inhibitor-heme model complexes can be a good descriptor in interpreting the inhibitor binding to CYP3A4 and the relative free energies in the gas phase are mainly responsible for the total binding free energies in water, although the desolvation can be a factor to affect the binding affinity of the inhibitors to CYP3A4. From the theozyme analysis of the X-ray structures for ketoconazole- and metyrapone-CYP3A4 complexes, the interaction free energy of the neighboring residues with each inhibitor in the active site is calculated to be about -3 kcal mol(-1) in water, whose the interaction energy and the desolvation free energy change are about -5 and 2 kcal mol(-1) , respectively.

The structural basis for homotropic and heterotropic cooperativity of midazolam metabolism by human cytochrome P450 3A4

Human cytochrome P450 3A4 (CYP3A4) metabolizes a significant portion of clinically relevant drugs and often exhibits complex steady-state kinetics that can involve homotropic and heterotropic cooperativity between bound ligands. In previous studies, the hydroxylation of the sedative midazolam (MDZ) exhibited homotropic cooperativity via a decrease in the ratio of 1'-OH-MDZ to 4-OH-MDZ at higher drug concentrations. In this study, MDZ exhibited heterotropic cooperativity with the antiepileptic drug carbamazepine (CBZ) with characteristic decreases in the 1'-OH-MDZ to 4-OH-MDZ ratios. To unravel the structural basis of MDZ cooperativity, we probed MDZ and CBZ bound to CYP3A4 using longitudinal T(1) nuclear magnetic resonance (NMR) relaxation and molecular docking with AutoDock 4.2. The distances calculated from longitudinal T(1) NMR relaxation were used during simulated annealing to constrain the molecules to the substrate-free X-ray crystal structure of CYP3A4. These simulations revealed that either two MDZ molecules or an MDZ molecule and a CBZ molecule assume a stacked configuration within the CYP3A4 active site. In either case, the proton at position 4 of the MDZ molecule was closer to the heme than the protons of the 1'-CH(3) group. In contrast, molecular docking of a single molecule of MDZ revealed that the molecule was preferentially oriented with the 1'-CH(3) position closer to the heme than position 4. This study provides the first detailed molecular analysis of heterotropic and homotropic cooperativity of a human cytochrome P450 from an NMR-based model. Cooperativity of ligand binding through direct interaction between stacked molecules may represent a common motif for homotropic and heterotropic cooperativity.

Evidence of CYP3A allosterism in vivo: analysis of interaction between fluconazole and midazolam

The allosteric effect of fluconazole (effector) on the formation of 1’-hydroxymidazolam (1’-OH-MDZ) and 4-hydroxymidazolam (4-OH-MDZ) from the CYP3A4/5 substrate, midazolam (MDZ), was examined in healthy volunteers. Following pre-treatment of fluconazole, AUC_4-OH/AUC_MDZ increased 35–62%, while AUC_1’-OH/AUC_MDZ decreased 5–37%; AUC_1’-OH/AUC_4-OH ratio decreased 46–58% by fluconazole and had no association with CYP3A5 genotype. 1’-OH-MDZ formation in vitro was more susceptible than 4-OH-MDZ formation to inhibition by fluconazole. Fluconazole decreased the intrinsic formation clearance ratio of 1’-OH-MDZ/4-OH-MDZ to an extent that was quantitatively comparable to in vivo observations. The elimination clearance of midazolam metabolites appeared unaffected by fluconazole. This study demonstrated that fluconazole alters midazolam product formation both in vivo and in vitro in a manner consistent with an allosteric interaction. The 1'-OH-MDZ/4-OH-MDZ ratio may serve as a biomarker of such interactions between midazolam, CYP3A4/5 and other putative effectors.

Structural features of cytochromes P450 and ligands that affect drug metabolism as revealed by X-ray crystallography and NMR

Cytochromes P450 (P450s) play a major role in the clearance of drugs, toxins, and environmental pollutants. Additionally, metabolism by P450s can result in toxic or carcinogenic products. The metabolism of pharmaceuticals by P450s is a major concern during the design of new drug candidates. Determining the interactions between P450s and compounds of very diverse structures is complicated by the variability in P450-ligand interactions. Understanding the protein structural elements and the chemical attributes of ligands that dictate their orientation in the P450 active site will aid in the development of effective and safe therapeutic agents. The goal of this review is to describe P450-ligand interactions from two perspectives. The first is the various structural elements that microsomal P450s have at their disposal to assume the different conformations observed in X-ray crystal structures. The second is P450-ligand dynamics analyzed by NMR relaxation studies.

Structure and function of cytochromes P450 2B: from mechanism-based inactivators to X-ray crystal structures and back

This article reviews work from the author dating back to 1978 and focuses on the structural basis of cytochrome P450 (P450) function using available contemporary techniques. Early studies used mechanism-based inactivators that bound to the protein moiety of hepatic P450s to try to localize the active site. Subsequent studies used cDNA cloning, heterologous expression, site-directed mutagenesis, and homology modeling based on multiple bacterial P450 X-ray crystal structures to predict the active sites of CYP2B enzymes with considerable accuracy. Breakthroughs in engineering and expression of mammalian P450s enabled us to determine our first X-ray crystal structure of ligand-free rabbit CYP2B4. To date, we have solved 11 CYP2B4 and three human CYP2B6 structures, which represent four significantly different conformations. The plasticity of CYP2B4 has been confirmed by deuterium exchange mass spectrometry and is substantiated by molecular dynamics simulations. In addition to major movement of secondary structure elements, more subtle reorientation of active site side chains, especially Phe206, Phe297, and Glu301, contributes to the ability of CYP2B enzymes to bind various ligands. Isothermal titration calorimetry has proven to be a useful tool for studying the thermodynamics of ligand binding to CYP2B4 and CYP2B6, and NMR has enabled study of ligand binding orientation in solution as an adjunct to X-ray crystallography. A major challenge remains to harness the power of the various approaches to facilitate prediction of CYP2B specificity and inhibition.

Conformational plasticity and structure/function relationships in cytochromes P450

The cytochrome P450s are a superfamily of enzymes that are found in all kingdoms of living organisms, and typically catalyze the oxidative addition of atomic oxygen to an unactivated C-C or C-H bond. Over 8000 nonredundant sequences of putative and confirmed P450 enzymes have been identified, but three-dimensional structures have been determined for only a small fraction of these. While all P450 enzymes for which structures have been determined share a common global fold, the flexibility and modularity of structure around the active site account for the ability of P450 enzymes to accommodate a vast number of structurally dissimilar substrates and support a wide range of selective oxidations. In this review, known P450 structures are compared, and some structural criteria for prediction of substrate selectivity and reaction type are suggested. The importance of dynamic processes such as redox-dependent and effector-induced conformational changes in determining catalytic competence and regio- and stereoselectivity is discussed, and noncrystallographic methods for characterizing P450 structures and dynamics, in particular, mass spectrometry and nuclear magnetic resonance spectroscopy are reviewed.

Structures of cytochrome P450 2B4 complexed with the antiplatelet drugs ticlopidine and clopidogrel

Prior X-ray crystal structures of rabbit cytochrome P450 2B4 (2B4) in complexes with various imidazoles have demonstrated markedly different enzyme conformations depending on the size of the inhibitor occupying the active site. In this study, structures of 2B4 were determined with the antiplatelet drugs clopidogrel and ticlopidine, which were expected to have greater freedom of movement in the binding pocket. Ticlopidine could be modeled into the electron density maps in two distinct orientations, both of which are consistent with metabolic data gathered with other mammalian P450 enzymes. Results of ligand docking and heme-induced NMR relaxation of drug protons showed that ticlopidine was preferentially oriented with the chlorophenyl group closest to the heme. Because of its stereocenter, clopidogrel was easier to fit in the electron density and exhibited a single orientation, which points the chlorophenyl ring toward the heme. The C(α) traces of both complexes aligned very well with each other and revealed a compact, closed structure that resembles the conformation observed in two previously determined 2B4 structures with the small molecule inhibitors 4-(4-chlorophenyl)imidazole and 1-(4-chlorophenyl)imidazole. The 2B4 active site is able to accommodate small ligands by moving only a small number of side chains, suggesting that ligand reorientation is energetically favored over protein conformational changes for binding of these similarly sized molecules. Adjusting both protein conformation and ligand orientation in the active site gives 2B4 the flexibility to bind to the widest range of molecules, while also being energetically favorable.

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.

The effect of chlorpyrifos-oxon and other xenobiotics on the human cytochrome P450-dependent metabolism of naphthalene and deet

Chlorpyrifos-oxon (CPO), a metabolite of chlorpyrifos, is a potent inhibitor of acetylcholinesterase and, although the neurotoxicological impact of this organophosphorus compound has been broadly studied both in vitro and in vivo, there are few studies of metabolic interactions of CPO with other xenobiotics. CPO significantly activated the production of 1-naphthol (5-fold), 2-naphthol (10-fold), trans-1,2-dihydro-1,2-naphthalenediol (1.5-fold), and 1,4-naphthoquinone from naphthalene by human liver microsomes (HLM). It was further demonstrated that the production of naphthalene metabolites by CYP2C8, 2C9*(1), 2C19, 2D6*(1), 3A4, 3A5, and 3A7 was activated by CPO, while the production of naphthalene metabolites by CYP1A1, 1A2, 1B1, and 2B6 was inhibited by CPO. CPO inhibited CYP1A2 production of naphthalene metabolites, while activating their production by CYP3A4. Similarly, CPO inhibited the production of N,N-diethyl-m-hydroxymethylbenzamide (BALC) from DEET by human liver microsomes, but activated the production of N-ethyl-m-toluamide (ET) from this substrate. CYP2B6, the most efficient isoform for BALC production, was inhibited by CPO, while CYP3A4, the most efficient isoform for ET production, was activated by CPO. CPO inhibited CYP2B6 production of both BALC and ET from DEET, but activated CYP3A4 production of ET, while inhibiting CYP3A4 BALC production. CPO appears to facilitate the binding of naphthalene to CYP3A4. This metabolic activation is independent of cytochrome b5, suggesting that activation of CYP3A4 by CPO is associated with a conformational change of the isoform rather than facilitating electron transfer.

High quality binding modes in docking ligands to proteins

Multiple near-optimal conformations of protein-ligand complexes provide a better chance for accurate representation of biomolecular interactions, compared with a single structure. We present ISE-dock--a docking program which is based on the iterative stochastic elimination (ISE) algorithm. ISE eliminates values that consistently lead to the worst results, thus optimizing the search for docking poses. It constructs large sets of such poses with no additional computational cost compared with single poses. ISE-dock is validated using 81 protein-ligand complexes from the PDB and its performance was compared with those of Glide, GOLD, and AutoDock. ISE-dock has a better chance than the other three to find more than 60% top single poses under RMSD = 2.0 A and more than 80% under RMSD = 3.0 A from experimental. ISE alone produced at least one 3.0 A or better solutions among the top 20 poses in the entire test set. In 98% of the examined molecules, ISE produced solutions that are closer than 2.0 A from experimental. Paired t-tests (PTT) were used throughout to assess the significance of comparisons between the performances of the different programs. ISE-dock provides more than 100-fold docking solutions in a similar time frame as LGA in AutoDock. We demonstrate the usefulness of the large near optimal populations of ligand poses by showing a correlation between the docking results and experiments that support multiple binding modes in p38 MAP kinase (Pargellis et al., Nat Struct Biol 2002;9:268-272] and in Human Transthyretin (Hamilton, Benson, Cell Mol Life Sci 2001;58:1491-1521).

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.

Structure-activity relationship and substrate-dependent phenomena in effects of ginsenosides on activities of drug-metabolizing P450 enzymes

Ginseng, a traditional herbal medicine, may interact with several co-administered drugs in clinical settings, and ginsenosides, the major active components of ginseng, may be responsible for these ginseng-drug interactions (GDIs). Results from previous studies on ginsenosides' effects on human drug-metabolizing P450 enzymes are inconsistent and confusing. Herein, we first evaluated the inhibitory effects of fifteen ginsenosides and sapogenins on human CYP1A2, CYP2C9, CYP2C19, CYP2D6 and CYP3A4 enzymes by using commercially available fluorescent probes. The structure-activity relationship of their effects on the P450s was also explored and a pharmacophore model was established for CYP3A4. Moreover, substrate-dependent phenomena were found in ginsenosides' effects on CYP3A4 when another fluorescent probe was used, and were further confirmed in tests with conventional drug probes and human liver microsomes. These substrate-dependent effects of the ginsenosides may provide an explanation for the inconsistent results obtained in previous GDI reports.

Substrate binding to cytochromes P450

P450s have attracted tremendous attention owing to not only their involvement in the metabolism of drug molecules and endogenous substrates but also the unusual nature of the reaction they catalyze, namely, the oxidation of unactivated C-H bonds. The binding of substrates to P450s, which is usually viewed as the first step in the catalytic cycle, has been studied extensively via a variety of biochemical and biophysical approaches. These studies were directed towards answering different questions related to P450s, including mechanism of oxidation, substrate properties, unusual substrate oxidation kinetics, function, and active-site features. Some of the substrate binding studies extending over a period of more than 40 years of dedicated work have been summarized in this review and categorized by the techniques employed in the binding studies.

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.

Effect of glutathione on homo- and heterotropic cooperativity in cytochrome P450 3A4

Glutathione (GSH) exerted a profound effect on the oxidation of 7-benzyloxy-4-(trifluoromethyl)coumarin (BFC) and 7-benzyloxyquinoline (BQ) by human liver microsomes as well as by CYP3A4-containing insect cell microsomes (Baculosomes). The cooperativity in O-debenzylation of both substrates is eliminated in the presence of 1-4mM GSH. Addition of GSH also increased the amplitude of the 1-PB induced spin shift with purified CYP3A4 and abolished the cooperativity of 1-PB or BFC binding. Changes in fluorescence of 6-bromoacetyl-2-dimethylaminonaphthalene attached to the cysteine-depleted mutant CYP3A4(C58,C64) suggest a GSH-induced conformational changes in proximity of alpha-helix A. Importantly, the K(S) value for formation of the GSH complex and the concentrations in which GSH decreases CYP3A4 cooperativity are consistent with the physiological concentrations of GSH in hepatocytes. Therefore, the allosteric effect of GSH on CYP3A4 may play an important role in regulation of microsomal monooxygenase activity in vivo.

QM/MM study of the active species of the human cytochrome P450 3A4, and the influence thereof of the multiple substrate binding

Cytochrome P450 3A4 is involved in the metabolism of 50% of all swallowed drugs. The enzyme functions by means of a high-valent iron-oxo species, called compound I (Cpd I), which is formed after entrance of the substrate to the active site. We explored the features of Cpd I using hybrid quantum mechanical/molecular mechanical calculations on various models that are either substrate-free or containing one and two molecules of diazepam as a substrate. Mössbauer parameters of Cpd I were computed. Our major finding shows that without the substrate, Cpd I tends to elongate its Fe-S bond, localize the radical on the sulfur, and form hydrogen bonds with A305 and T309, which may hypothetically lead to Cpd I consumption by H-abstraction. However, the positioning of diazepam close to Cpd I, as enforced by the effector molecule, was found to strengthen the NH...S interactions of the conserved I443 and G444 residues with the proximal cysteinate ligand. These interactions are known to stabilize the Fe-S bond, and as such, the presence of the substrate leads to a shorter Fe-S bond and it prevents the localization of the radical on the sulfur. This diazepam-Cpd I stabilization was manifested in the 1W0E conformer. The effector substrate did not influence Cpd I directly but rather by positioning the active substrate close to Cpd I, thus displacing the hydrogen bonds with A305 and T309, and thereby giving preference to substrate oxidation. It is hypothesized that these effects on Cpd I, promoted by the restrained substrate, may be behind the special metabolic behavior observed in cases of multiple substrate binding (also called cooperative binding). This restraint constitutes a mechanism whereby substrates stabilize Cpd I sufficiently long to affect monooxygenation by P450s at the expense of Cpd I destruction by the protein residues.

Cooperative binding of acetaminophen and caffeine within the P450 3A4 active site

Acetaminophen (N-acetyl-p-aminophenol, APAP) is a commonly used analgesic/antipyretic. When oxidized by P450, a toxic APAP metabolite is generated. Human P450 3A4 was expressed in Escherichia coli , purified, and reconstituted using artificial liposomes. Oxidation of APAP by P450 3A4, as detected by the formation of its glutathione adduct, was found to exhibit negative homotropic cooperativity with a Hill coefficient of 0.7. In the presence of caffeine, the observed kinetics were close to classical Michaelis-Menten kinetics with a Hill coefficient approaching 1. In order to probe for a potential repositioning of APAP within the P450 3A4 pocket in the presence of caffeine, NMR T1 paramagnetic relaxation techniques were used to calculate distances from the P450 3A4 heme iron to protons of APAP alone and in the presence of caffeine. Both APAP and caffeine were found to bind at the active site in proximity to the heme iron. When APAP was incubated with P450 3A4, the acetamido group of APAP was found to be closest to the heme iron consistent with the amide group of APAP weakly associating with the heme iron. The addition of caffeine disrupted the ability of APAP to coordinate with the heme iron of P450 3A4 and enhanced the rate of oxidation to its toxic metabolite.

In Vitro Prediction and in Vivo Verification of Enantioselective Human Tofisopam Metabolite Profiles

In vitro studies were conducted to elucidate the metabolic profiles of and the enzymes responsible for the metabolism of (R)- and (S)-tofisopam (1-(3,4-dimethoxyphenyl)-5-ethyl-7,8-dimethoxy-4-methyl-5H-2,3-benzodiazepine). Large differences were observed between the two enantiomers. The major metabolite in incubations of 500 ng/ml (approximately 1.3 microM) (R)-tofisopam in human liver microsomes corresponded to demethylation of the methoxy group at the 4-position of the phenyl ring (M3). Incubating (R)-tofisopam with recombinant cytochrome P450 (P450) or with human liver microsomes and isoform-selective P450 chemical inhibitors indicated that M3 was primarily catalyzed by CYP2C9. Similar incubations with S-tofisopam indicated that the primary metabolite was due to demethylation of the methoxy group at the 7-position of the benzodiazepine ring (M1), and this reaction was catalyzed primarily by CYP3A4. The primary metabolites of both enantiomers were further demethylated to form a common didemethylated metabolite (M5) where the methoxy groups at positions 4 and 7 are demethylated. Analysis of plasma and urine samples from human clinical trials confirmed the in vitro observations. Subjects orally treated with 200 mg b.i.d. (R)-tofisopam had a 2-h M1/M3 plasma ratio of 1:29 and a ratio of 1:123 in urine, whereas patients orally administered (S)-tofisopam at 150 mg/kg t.i.d. had opposite M1 to M3 ratios of 8:1 in plasma and 6:1 in urine.

Ligand binding to cytochrome P450 3A4 in phospholipid bilayer nanodiscs: the effect of model membranes

The membrane-bound protein cytochrome P450 3A4 (CYP3A4) is a major drug-metabolizing enzyme. Most studies of ligand binding by CYP3A4 are currently carried out in solution, in the absence of a model membrane. Therefore, there is little information concerning the membrane effects on CYP3A4 ligand binding behavior. Phospholipid bilayer Nanodiscs are a novel model membrane system derived from high density lipoprotein particles, whose stability, monodispersity, and consistency are ensured by their self-assembly. We explore the energetics of four ligands (6-(p-toluidino)-2-naphthalenesulfonic acid (TNS), alpha-naphthoflavone (ANF), miconazole, and bromocriptine) binding to CYP3A4 incorporated into Nanodiscs. Ligand binding to Nanodiscs was monitored by a combination of environment-sensitive ligand fluorescence and ligand-induced shifts in the fluorescence of tryptophan residues present in the scaffold proteins of Nanodiscs; binding to the CYP3A4 active site was monitored by ligand-induced shifts in the heme Soret band absorbance. The dissociation constants for binding to the active site in CYP3A4-Nanodiscs were 4.0 microm for TNS, 5.8 microm for ANF, 0.45 microm for miconazole, and 0.45 microm for bromocriptine. These values are for CYP3A4 incorporated into a lipid bilayer and are therefore presumably more biologically relevant that those measured using CYP3A4 in solution. In some cases, affinity measurements using CYP3A4 in Nanodiscs differ significantly from solution values. We also studied the equilibrium between ligand binding to CYP3A4 and to the membrane. TNS showed no marked preference for either environment; ANF preferentially bound to the membrane, and miconazole and bromocriptine preferentially bound to the CYP3A4 active site.

Structural dynamics of the cooperative binding of organic molecules in the human cytochrome P450 3A4

Cytochrome P450 3A4 (CYP3A4) is a key enzyme responsible for the metabolism of 50% of all orally administered drugs which exhibit an intriguing kinetic behavior typified by a sigmoidal dependence of the reaction velocity on the substrate concentration. There is evidence for the binding of two substrates in the active site of the enzyme, but the mechanism of this cooperative binding is unclear. Diazepam is such a drug that undergoes metabolism by CYP3A4 with sigmoidal dependence. Metabolism is initiated by hydrogen atom abstraction from the drug. To understand the factors that determine the cooperative binding and the juxtaposition of the C-H bond undergoing abstraction, we carried out molecular dynamics simulations for two enzymatic conformers and examined the differences between the substrate-free and the bound enzymes, with one and two diazepam molecules. Our results indicate that the effector substrate interacts both with the active substrate and with the enzyme, and that this interaction results in side chain reorientation with relatively minor long-range effects. In accord with experiment, we find that F304, in the interface between the active and effector binding sites, is a key residue in the mechanism of cooperative binding. The addition of the effector substrate stabilizes F304 and its environment, especially F213, and induces a favorable orientation of the active substrate, leading to a short distance between the targeted hydrogen for abstraction and the active species of the enzyme. In addition, in one conformer of the enzyme, residue R212 may strongly interact with F304 and counteract the effector's impact on the enzyme.