Orientation of caffeine within the active site of human cytochrome P450 1A2 based on NMR longitudinal (T1) relaxation measurements

Hyperpolarized 15N caffeine, a potential probe of liver function and perfusion

PURPOSE

To demonstrate hyperpolarization of 15N-caffeine and report exploratory findings as a potential probe of liver function and perfusion.

METHODS

An amorphous formulation of [1,3-15N2]caffeine was developed for hyperpolarization via dissolution dynamic nuclear polarization. Polarizer hardware was augmented to support monitoring of solid-state 15N MR signals during the buildup of hyperpolarization. Liquid state hyperpolarized 15N MR signals were obtained in a preclinical 3T magnet by interfacing an external spectrometer console with home-built RF surface coils. 15N signal decay constants were estimated in H2O and in vivo in liver and brain regions of rats at 3 T. Decays were also measured at 9.4 T to assess the effect of B0, and in the presence of albumin to assess the impact of protein binding.

RESULTS

Polarization levels of 3.5% and aqueous T1 relaxation times of nearly 200 s were attained for both N1 and N3 positions at 3 T. Shorter apparent decay constants were observed in vivo, ranging from 25 s to 43 s, with modest extensions possible by exploiting competitive binding of iophenoxate with plasma albumin. Downstream products of caffeine could not be detected on in vivo 15N-MR spectra of the liver region, even with metabolic stimulation byβ $$ \beta $$ -naphthoflavone treatment. Considering the high perfusion rate of brain, persistence of caffeine signal in this region is consistent with potential value as a perfusion imaging agent.

CONCLUSION

These results establish the feasibility of hyperpolarization of hyperpolarized 15N-caffeine, but further work is necessary to establish the role of this new agent to probe liver metabolism and perfusion.

Active Binding Modes of Caffeine with Cytochrome P450 1A2 Determine Its Metabolite Profiles

Caffeine is a very common kind of nervous stimulant, and it is primarily metabolized by Cytochrome P450 1A2 (CYP1A2) in the human body. Over the years, determining the interactions between caffeine and CYP1A2 has been a tough issue. The active binding modes and the catalytic regioselectivity of the metabolism between CYP1A2 and caffeine remain unclear. Here, to investigate the interactions between CYP1A2 and caffeine, we constructed the all-sequence CYP1A2-caffeine-membrane system using a multiple template approach. According to our simulation results, four active binding modes between CYP1A2 and caffeine that correspond to the four metabolic sites of caffeine are determined. What is more, a pre-reaction state for the CYP1A2-catalyzed reaction at caffeine's N3 site is identified. A more preponderant active binding mode might be the reason why the N3 site of caffeine becomes the primary metabolic site. Our findings could enhance our knowledge of the interactions between CYP1A2 and caffeine and help us better understand the regioselectivity of the metabolism between CYP1A2 and caffeine.

Single Mutations in Cytochrome P450 Oxidoreductase Can Alter the Specificity of Human Cytochrome P450 1A2-Mediated Caffeine Metabolism

A unique cytochrome P450 (CYP) oxidoreductase (CPR) sustains activities of human microsomal CYPs. Its function requires toggling between a closed conformation enabling electron transfers from NADPH to FAD and then FMN cofactors and open conformations forming complexes and transferring electrons to CYPs. We previously demonstrated that distinct features of the hinge region linking the FAD and FMN domain (FD) modulate conformer poses and their interactions with CYPs. Specific FD residues contribute in a CYP isoform-dependent manner to the recognition and electron transfer mechanisms that are additionally modulated by the structure of CYP-bound substrate. To obtain insights into the underlying mechanisms, we analyzed how hinge region and FD mutations influence CYP1A2-mediated caffeine metabolism. Activities, metabolite profiles, regiospecificity and coupling efficiencies were evaluated in regard to the structural features and molecular dynamics of complexes bearing alternate substrate poses at the CYP active site. Studies reveal that FD variants not only modulate CYP activities but surprisingly the regiospecificity of reactions. Computational approaches evidenced that the considered mutations are generally in close contact with residues at the FD-CYP interface, exhibiting induced fits during complexation and modified dynamics depending on caffeine presence and orientation. It was concluded that dynamic coupling between FD mutations, the complex interface and CYP active site exist consistently with the observed regiospecific alterations.

A Modified Arrhenius Approach to Thermodynamically Study Regioselectivity in Cytochrome P450-Catalyzed Substrate Conversion

The regio- (and stereo-)selectivity and specific activity of cytochrome P450s are determined by the accessibility of potential sites of metabolism (SOMs) of the bound substrate relative to the heme, and the activation barrier of the regioselective oxidation reaction(s). The accessibility of potential SOMs depends on the relative binding free energy (ΔΔGbind ) of the catalytically active substrate-binding poses, and the probability of the substrate to adopt a transition-state geometry. An established experimental method to measure activation energies of enzymatic reactions is the analysis of reaction rate constants at different temperatures and the construction of Arrhenius plots. This is a challenge for multistep P450-catalyzed processes that involve redox partners. We introduce a modified Arrhenius approach to overcome the limitations in studying P450 selectivity, which can be applied in multiproduct enzyme catalysis. Our approach gives combined information on relative activation energies, ΔΔGbind values, and collision entropies, yielding direct insight into the basis of selectivity in substrate conversion.

Binding Modes and Metabolism of Caffeine

A correct estimate of ligand binding modes and a ratio of their occupancies is crucial for calculations of binding free energies. The newly developed method BLUES combines molecular dynamics with nonequilibrium candidate Monte Carlo. Nonequilibrium candidate Monte Carlo generates a plethora of possible binding modes and molecular dynamics enables the system to relax. We used BLUES to investigate binding modes of caffeine in the active site of its metabolizing enzyme Cytochrome P450 1A2 with the aim of elucidating metabolite-formation profiles at different concentrations. Because the activation energies of all sites of metabolism do not show a clear preference for one metabolite over the others, the orientations in the active site must play a key role. In simulations with caffeine located in a spacious pocket above the I-helix, it points N3 and N1 to the heme iron, whereas in simulations where caffeine is in close proximity to the heme N7 and C8 are preferably oriented toward the heme iron. We propose a mechanism where at low caffeine concentrations caffeine binds to the upper part of the active site, leading to formation of the main metabolite paraxanthine. On the other hand, at high concentrations two molecules are located in the active site, forcing one molecule into close proximity to the heme and yielding metabolites theophylline and trimethyluretic acid. Our results offer an explanation of previously published experimental results.

Increased Phenacetin Oxidation upon the L382V Substitution in Cytochrome P450 1A2 is Associated with Altered Substrate Binding Orientation

Leucine382 of cytochrome P450 1A2 (CYP1A2) plays an important role in binding and O-dealkylation of phenacetin, with the L382V mutation increasing substrate oxidation (Huang and Szklarz, 2010, Drug Metab. Dispos. 38:1039–1045). This was attributed to altered substrate binding orientation, but no direct experimental evidence had been available. Therefore, in the current studies, we employed nuclear magnetic resonance (NMR) longitudinal (T₁) relaxation measurements to investigate phenacetin binding orientations within the active site of CYP1A2 wild type (WT) and mutants. Paramagnetic relaxation time (T_1P) for each proton of phenacetin was calculated from the T₁ value obtained from the enzymes in ferric and ferrous-CO state in the presence of phenacetin, and used to model the orientation of phenacetin in the active site. All aromatic protons of phenacetin were nearly equidistant from the heme iron (6.34–8.03 Å). In contrast, the distance between the proton of the –OCH₂– group, which is abstracted during phenacetin oxidation, and the heme iron, was much shorter in the L382V (5.93 Å) and L382V/N312L (5.96 Å) mutants compared to the N312L mutant (7.84 Å) and the wild type enzyme (6.55 Å), consistent with modeling results. These studies provide direct evidence for the molecular mechanism underlying increased oxidation of phenacetin upon the L382V mutation.

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.

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β.

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.

NMR-derived models of amidopyrine and its metabolites in complexes with rabbit cytochrome P450 2B4 reveal a structural mechanism of sequential N-dealkylation

To understand the molecular basis of sequential N-dealkylation by cytochrome P450 2B enzymes, we studied the binding of amidopyrine (AP) as well as the metabolites of this reaction, desmethylamidopyrine (DMAP) and aminoantipyrine (AAP), using the X-ray crystal structure of rabbit P450 2B4 and two nuclear magnetic resonance (NMR) techniques: saturation transfer difference (STD) spectroscopy and longitudinal (T(1)) relaxation NMR. Results of STD NMR of AP and its metabolites bound to P450 2B4 were similar, suggesting that they occupy similar niches within the enzyme's active site. The model-dependent relaxation rates (R(M)) determined from T(1) relaxation NMR of AP and DMAP suggest that the N-linked methyl is closest to the heme. To determine the orientation(s) of AP and its metabolites within the P450 2B4 active site, we used distances calculated from the relaxation rates to constrain the metabolites to the X-ray crystal structure of P450 2B4. Simulated annealing of the complex revealed that the metabolites do indeed occupy similar hydrophobic pockets within the active site, while the N-linked methyls are free to rotate between two binding modes. From these bound structures, a model of N-demethylation in which the N-linked methyl functional groups rotate between catalytic and noncatalytic positions was developed. This study is the first to provide a structural model of a drug and its metabolites complexed to a cytochrome P450 based on NMR and to provide a structural mechanism for how a drug can undergo sequential oxidations without unbinding. The rotation of the amide functional group might represent a common structural mechanism for N-dealkylation reactions for other drugs such as the local anesthetic lidocaine.

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.

Functional analysis of phenylalanine residues in the active site of cytochrome P450 2C9

The two published crystal structures of cytochrome P450 2C9, complexed with ( S)-warfarin or flurbiprofen, implicate a cluster of three active site phenylalanine residues (F100, F114, F476) in ligand binding. However, these three residues appear to interact differently with these two ligands based on the static crystal structures. To elucidate the importance of CYP2C9's active site phenylalanines on substrate binding, orientation, and catalytic turnover, a series of leucine and tryptophan mutants were constructed and their interactions with ( S)-warfarin and ( S)-flurbiprofen examined. The F100-->L mutation had minor effects on substrate binding and metabolism of each substrate. In contrast, the F114L and F476L mutants exhibited substantially reduced ( S)-warfarin metabolism and altered hydroxy metabolite profiles but only modestly decreased nonsteroidal antiinflammatory drug (NSAID) turnover while maintaining product regioselectivity. The F114-->W and F476-->W mutations also had opposing effects on ( S)-warfarin versus NSAID turnover. Notably, the F476W mutant increased the efficiency of ( S)-warfarin metabolism 5-fold, yet decreased the efficiency of ( S)-flurbiprofen turnover 20-fold. (1)H NMR T 1 relaxation studies suggested a slightly closer positioning of ( S)-warfarin to the heme in the F476W mutant relative to the wild-type enzyme, and stoichiometry studies indicated enhanced coupling of reducing equivalents to product formation for ( S)-warfarin, again in contrast to effects observed with ( S)-flurbiprofen. These data demonstrate that F114 and F476, but not F100, influence ( S)-warfarin's catalytic orientation. Differential interactions of F476 mutants with the two substrates suggest that their catalytically productive binding modes are not superimposable.

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.

Substrate proton to heme distances in CYP2C9 allelic variants and alterations by the heterotropic activator, dapsone

CYP2C9 polymorphisms result in reduced enzyme catalytic activity and greater activation by effector molecules as compared to wild-type protein, with the mechanism(s) for these changes in activity not fully elucidated. Through T(1) NMR and spectral binding analyses, mechanism(s) for these differences in behavior of the variant proteins (CYP2C9.2, CYP2C9.3, and CYP2C9.5) as compared to CYP2C9.1 were assessed. Neither altered binding affinity nor substrate (flurbiprofen) proton to heme-iron distances differed substantially among the four enzymes. Co-incubation with dapsone resulted in reduced substrate proton to heme-iron distances for all enzymes, providing at least a partial mechanism for the activation of CYP2C9 variants by dapsone. In summary, neither altered binding affinity nor substrate orientation appear to be major factors in the reduced catalytic activity noted in the CYP2C9 variants, but dapsone co-incubation caused similar changes in substrate proton to heme-iron distances suggesting at least partial common mechanisms in the activation of the CYP2C9 forms.

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.

Ensemble modeling of substrate binding to cytochromes P450: analysis of catalytic differences between CYP1A orthologs

A novel application of modeling and docking approaches involving ensembles of homology models is used to understand structural bases underlying subtle catalytic differences between related cytochromes P450 (CYPs). Mammalian CYP1A1s and fish CYP1As are orthologous enzymes with similar substrate preferences. With some substrates (3,3',4,4'-tetrachlorobiphenyl, TCB) oxidation rates differ by orders of magnitude, while others (e.g., benzo[a]pyrene; B[a]P) are oxidized at similar rates but with somewhat differing regiospecificity. These two environmental chemical substrates (TCB and B[a]P) as well as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) were docked to multiple models of rat, human, scup, and/or killifish CYP1As, based on multiple templates, retaining multiple poses from each model, giving ensembles of docked poses for each species. With TCB, more poses were observed closer to the heme in ensembles of rat or human CYP1A1 than of killifish CYP1A. Analysis of interacting residues suggested that differences in TCB pose distributions are due primarily to Leu387 and Val230 in killifish CYP1A. In silico mutations L387V and V230G enabled TCB to dock closer to the heme in killifish CYP1A. Mutating additional interacting residues (Ala127, Thr233, Asn317, and Tyr386) of killifish CYP1A to the corresponding residues of human CYP1A1 resulted in TCB pose distributions nearly identical with those of human CYP1A1. Docking of TCDD to sets of consensus models of killifish, rat, and human CYP1As showed species differences similar to those with TCB, but with further structural constraints possibly contributing to slower oxidation of TCDD. Docking B[a]P to sets of consensus models of the human and fish CYP1As yielded frequencies of substrate orientations correlating with known regiospecificities for metabolism of B[a]P by these enzymes. The results demonstrate the utility of this ensemble modeling method, which can account for uncertainty inherent in homology modeling and docking by producing statistical distributions of ligand positions.

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

Recent studies have indicated that CYP3A4 exhibits non-Michaelis-Menten kinetics for numerous substrates. Both homo- and heterotropic activation have been reported, and kinetic models have suggested multiple substrates within the active site. We provide some of the first physicochemical data supporting the hypothesis of allosteric substrate binding within the CYP3A4 active site. Midazolam (MDZ) is metabolized by CYP3A4 to two hydroxylated metabolites, 1'- and 4-hydroxymidazolam. Incubations using purified CYP3A4 and MDZ showed that both alpha-naphthoflavone (alpha-NF) and testosterone affect the ratio of formation rates of 1'- and 4-hydroxymidazolam. Similar to previous reports, alpha-NF was found to promote formation of 1'-hydroxymidazolam, while testosterone stimulated formation of 4-hydroxymidazolam. NMR was used to measure the closest approach of individual MDZ protons to the paramagnetic heme iron of CYP3A4 using paramagnetic T(1) relaxation measurements. Solutions of 0.2 microM CYP3A4 with 500 microM MDZ resulted in calculated distances between 7.4 and 8.3 A for all monitored MDZ protons. The distances were statistically equivalent for all protons except C3-H and were consistent with the rotation within the active site or sliding parallel to the heme plane. When 50 microM alpha-NF was added, proton-heme iron distances ranged from 7.3 to 10.0 A. Consistent with kinetics of activation, the 1' position was situated closest to the heme, while the fluorophenyl 5-H proton was the furthest. Proton-heme iron distances for MDZ with CYP3A4 and 50 microM testosterone ranged from 7.7 to 9.0 A, with the flourophenyl 5-H proton furthest from the heme iron and the C4-H closest to the heme, also consistent with kinetic observations. When titrated with CYP3A4 in the presence of MDZ, testosterone and alpha-NF resonances themselves exhibited significant broadening and enhanced relaxation rates, indicating that these effector molecules were also bound within the CYP3A4 active site near the paramagnetic heme iron. These results suggest that the effector exerts its cooperative effects on MDZ metabolism through simultaneous binding of MDZ and effector near the CYP3A4 heme.

Metabolic Regio- and Stereoselectivity of Cytochrome P450 2D6 towards 3,4-Methylenedioxy-N-alkylamphetamines: in Silico Predictions and Experimental Validation

A series of 3,4-methylenedioxy-N-alkylamphetamines (MDAAs) were automatically docked and subjected to molecular dynamics (MD) simulations in a cytochrome P450 2D6 (CYP2D6) protein model. The predicted substrate binding orientations, sites of oxidation, and relative reactivities were compared to the experimental data of wild-type and Phe120Ala mutant CYP2D6. Automated docking results were not sufficient to accurately rationalize experimental binding orientations of 3,4-methylenedioxy-N-methylamphetamine (MDMA) in the two enzymes as measured with spin lattice relaxation NMR. Nevertheless, the docking results could be used as starting structures for MD simulations. Predicted binding orientations of MDMA and sites of oxidation of the MDAAs derived from MD simulations matched well with the experimental data. It appeared the experimental results were best described in MD simulations considering the nitrogen atoms of the MDAAs in neutral form. Differences in regioselectivity and stereoselectivity in the oxidative metabolism of the MDAAs by the Phe120Ala mutant CYP2D6 were correctly predicted, and the effects of the Phe120Ala mutation could be rationalized as well.

CYP2C9 genotype-dependent effects on in vitro drug-drug interactions: switching of benzbromarone effect from inhibition to activation in the CYP2C9.3 variant

The CYP2C9.3 variant exhibits marked decreases in substrate turnover compared with the wild-type enzyme, but little is known regarding the effect this variant form may have on the occurrence of drug-drug interactions. To examine this possibility, the effect of the potent CYP2C9 inhibitor, benzbromarone, was studied with regard to CYP2C9.1- and CYP2C9.3-mediated flurbiprofen metabolism to evaluate whether the variant enzyme exhibits differential inhibition kinetics. Although benzbromarone inhibited CYP2C9.1 activity as expected, CYP2C9.3-mediated flurbiprofen 4'-hydroxylation was activated in the presence of benzbromarone. T1 relaxation studies revealed little change in distances of flurbiprofen protons from the heme iron of either CYP2C9.1 or CYP2C9.3 in the presence of benzbromarone compared with flurbiprofen alone. Spectral binding studies were also performed to investigate whether benzbromarone affected substrate binding, with the addition of benzbromarone having little effect on flurbiprofen-binding affinity in both CYP2C9.1 and CYP2C9.3. Docking studies with the 2C9.1 structure crystallized with a closed active site identified multiple but overlapping subsites with sufficient space for benzbromarone binding in the enzyme when flurbiprofen was positioned closest to the heme. If the closed conformation of 2C9.3 is structurally similar to 2C9.1, as expected for the conservative I359L mutation, then the dynamics of benzbromarone binding may account for the switching of drug interaction effects. In conclusion, the I359L amino acid substitution found in CYP2C9.3 not only reduces metabolism compared with CYP2C9.1 but can also dramatically alter inhibitor effects, suggesting that differential degrees of drug inhibition interactions may occur in individuals with this variant form of CYP2C9.

Effector-Mediated Alteration of Substrate Orientation in Cytochrome P450 2C9†

Cytochrome P450 2C9 (CYP2C9)-mediated flurbiprofen 4'-hydroxylation is activated by the presence of dapsone resulting in reduction of the K(m) for flurbiprofen hydroxylation and an increase in V(m). Previous spectral binding studies have demonstrated that the binding of flurbiprofen with CYP2C9 is increased (decrease in K(S)) by the presence of dapsone. We hypothesized that the two compounds are simultaneously in the active site with the presence of dapsone causing flurbiprofen to be oriented more closely to the heme. T(1) relaxation rates determined by NMR were used to estimate the distances of protons on these compounds from the paramagnetic heme-iron center. Samples contained 0.014 microM CYP2C9 and 145 microM flurbiprofen in the presence and absence of 100 microM dapsone. Estimated distances of various flurbiprofen protons from the heme ranged from 4.2 to 4.5 A in the absence of dapsone and from 3.2 to 3.8 A in the presence of dapsone. The 4' proton of flurbiprofen, the site of metabolism, showed one of the greatest differences in distance from the heme in the presence of dapsone, 3.50 A, as compared to the absence of dapsone, 4.41 A. Dapsone protons were less affected, being 4.40 A from the heme in the absence of flurbiprofen and 4.00-4.01 A from the heme in the presence of flurbiprofen. Molecular modeling studies were also performed to corroborate the relative orientations of flurbiprofen and dapsone in the active site of CYP2C9. Shift of the 4' proton of flurbiprofen closer to the heme iron of CYP2C9 in the presence of dapsone may play a role in activation.

Six Novel NonsynonymousCYP1A2Gene Polymorphisms: Catalytic Activities of the Naturally Occurring Variant Enzymes

Six novel nonsynonymous nucleotide alterations were found in the cytochrome P450 1A2 gene in a Japanese population, which resulted in the following amino acid substitutions: T83M, E168Q, F186L, S212C, G299A, and T438I. These individuals were heterozygous for the amino acid substitutions. The potential functional alterations caused by the amino acid substitutions were characterized by a cDNA-mediated expression system using Chinese hamster V79 cells. Among the six CYP1A2 variants, F186L showed the most profound and statistically significant reduction in O-deethylation of phenacetin and 7-ethoxyresorufin. Kinetic analyses performed for the ethoxyresorufin O-deethylation revealed that the Vmax of the F186L variant was approximately 5% of that of the CYP1A2 wild type, despite a 5-fold lower Km value of the variant, the consequence of which was reduced enzymatic activity toward the substrate. Thus, for the first time, phenylalanine at residue 186 is suggested to be a critical amino acid for catalytic activity.