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

Interaction of CYP3A4 with caffeine: First insights into multiple substrate binding

Human cytochrome P450 3A4 (CYP3A4) is a major drug-metabolizing enzyme that shows extreme substrate promiscuity. Moreover, its large and malleable active site can simultaneously accommodate several substrate molecules of the same or different nature, which may lead to cooperative binding and allosteric behavior. Due to difficulty of crystallization of CYP3A4-substrate complexes, it remains unknown how multiple substrates can arrange in the active site. We determined crystal structures of CYP3A4 bound to three and six molecules of caffeine, a psychoactive alkaloid serving as a substrate and modulator of CYP3A4. In the ternary complex, one caffeine binds to the active site suitably for C8-hydroxylation, most preferable for CYP3A4. In the senary complex, three caffeine molecules stack parallel to the heme with the proximal ligand poised for 3-N-demethylation. However, the caffeine stack forms extensive hydrophobic interactions that could preclude product dissociation and multiple turnovers. In both complexes, caffeine is also bound in the substrate channel and on the outer surface known as a peripheral site. At all sites, aromatic stacking with the caffeine ring(s) is likely a dominant interaction, while direct and water-mediated polar contacts provide additional stabilization for the substrate-bound complexes. Protein-ligand interactions via the active site R212, intrachannel T224, and peripheral F219 were experimentally confirmed, and the latter two residues were identified as important for caffeine association. Collectively, the structural, spectral, and mutagenesis data provide valuable insights on the ligand binding mechanism and help better understand how purine-based pharmaceuticals and other aromatic compounds could interact with CYP3A4 and mediate drug-drug interactions.

Structural Dynamics of Cytochrome P450 3A4 in the Presence of Substrates and Cytochrome P450 Reductase

Cytochrome P450 3A4 (CYP3A4) is the most important drug-metabolizing enzyme in humans and has been associated with harmful drug interactions. The activity of CYP3A4 is known to be modulated by several compounds and by the electron transfer partner, cytochrome P450 reductase (CPR). The underlying mechanism of these effects, however, is poorly understood. We have used hydrogen-deuterium exchange mass spectrometry to investigate the impact of binding of CPR and of three different substrates (7-benzyloxy-4-trifluoromethyl-coumarin, testosterone, and progesterone) on the conformational dynamics of CYP3A4. Here, we report that interaction of CYP3A4 with substrates or with the oxidized or reduced forms of CPR leads to a global rigidification of the CYP3A4 structure. This was evident from the suppression of deuterium exchange in several regions of CYP3A4, including regions known to be involved in protein-protein interactions (helix C) and substrate binding and specificity (helices B' and E, and loop K/β1). Furthermore, the bimodal isotopic distributions observed for some CYP3A4-derived peptides were drastically impacted upon binding to CPR and/or substrates, suggesting the existence of stable CYP3A4 conformational populations that are perturbed by ligand/CPR binding. The results have implications for understanding the mechanisms of ligand binding, allostery, and catalysis in CYP enzymes.

Using Small Molecules to Enhance P450 OleT Enzyme Activity in Situ

Cytochrome P450 OleT is a fatty acid decarboxylase that catalyzes the production of olefins with biofuel and synthetic applications. However, the relatively sluggish catalytic efficiency of the enzyme limits its applications. Here, we report the application of a novel class of benzene containing small molecules to improve the OleT activity. The UV-Vis spectroscopy study and molecular docking results confirmed the high proximity of the small molecules to the heme group of OleT. Up to 6-fold increase of product yield has been achieved in the small molecule-modulated enzymatic reactions. Our work thus sheds the light to the application of small molecules to increase the OleT catalytic efficiency, which could be potentially used for future olefin productions.

Deciphering Key Interactions of Ligands with CYP3A4-Template* system

Cytochrome P450 (CYP)-mediated metabolisms are often associated with biological and toxicological events of chemicals. A major hepatic enzyme, CYP3A4, showed clear distinctions on their catalyses even among ligands having resemble structures. To better understand mechanisms of their distinct catalyses, possible associations of ligand interactions at specific parts of CYP3A4 residues were investigated using CYP3A4-Template system developed (DMPK 2019 and 2020). A placement was available selectively for CYP3A4-mediated R-thalidomide 5-oxidation on Template, but not for the 5’-oxidation and the S-isomer oxidations. Similar placements were generated for pomalidomide (4-amino-thalidomide), but not for a poor ligand, lenalidomide (3-deoxy-pomalidomide). The latter ligand took placements lacking IJK-Interaction or sticking the 4-amino part beyond the facial-side wall on Template. A placement was available for the tert-butyl oxidation of terfenadine, but not for an analog, ebastine. Their interactions with upper-Cavity-2 residue were expected to differ at their sites of oxygen substituents. Some phenolic antioxidants behave distinctly toward biological oxidations in vitro and in vivo. Butylated hydroxytoluene is oxidized to the peroxy-derivative in vitro, but solely to the oxidized metabolites at the benzyl and tert-butyl methyl positions in vivo. Involvement of CYP3A4 were suggested for all the three reactions from the placements on Template. Tocopherols were also applied on Template for the oxidations for chroman and side-chain terminals. The primary placement was suggested to undergo the futile-recycling through formation of the peroxide intermediate subsequently to lead the substantial lack of the CYP3A4-mediated oxidation. These data suggest the effectiveness of CYP3A4-Template assessment to understand the causal basis of poor oxidations and also to verify the in vivo contribution of CYP3A4-mediated peroxidative reactions.

Human Family 1-4 cytochrome P450 enzymes involved in the metabolic activation of xenobiotic and physiological chemicals: an update

This is an overview of the metabolic activation of drugs, natural products, physiological compounds, and general chemicals by the catalytic activity of cytochrome P450 enzymes belonging to Families 1-4. The data were collected from > 5152 references. The total number of data entries of reactions catalyzed by P450s Families 1-4 was 7696 of which 1121 (~ 15%) were defined as bioactivation reactions of different degrees. The data were divided into groups of General Chemicals, Drugs, Natural Products, and Physiological Compounds, presented in tabular form. The metabolism and bioactivation of selected examples of each group are discussed. In most of the cases, the metabolites are directly toxic chemicals reacting with cell macromolecules, but in some cases the metabolites formed are not direct toxicants but participate as substrates in succeeding metabolic reactions (e.g., conjugation reactions), the products of which are final toxicants. We identified a high level of activation for three groups of compounds (General Chemicals, Drugs, and Natural Products) yielding activated metabolites and the generally low participation of Physiological Compounds in bioactivation reactions. In the group of General Chemicals, P450 enzymes 1A1, 1A2, and 1B1 dominate in the formation of activated metabolites. Drugs are mostly activated by the enzyme P450 3A4, and Natural Products by P450s 1A2, 2E1, and 3A4. Physiological Compounds showed no clearly dominant enzyme, but the highest numbers of activations are attributed to P450 1A, 1B1, and 3A enzymes. The results thus show, perhaps not surprisingly, that Physiological Compounds are infrequent substrates in bioactivation reactions catalyzed by P450 enzyme Families 1-4, with the exception of estrogens and arachidonic acid. The results thus provide information on the enzymes that activate specific groups of chemicals to toxic metabolites.

Conformational selection is present in ligand binding to cytochrome P450 19A1 lipoprotein nanodiscs

Cytochromes P450 (CYPs) display remarkable plasticity in their ability to bind substrates and catalyze a broad array of chemical reactions. Herein we evaluate binding of androstenedione, testosterone, and 7-hydroxyflavone to CYP19A1, also known as aromatase, in phospholipid nanodiscs by stopped-flow UV-vis spectroscopy. Exponential fitting of the kinetic traces supports the possibility of a multi-step binding mechanism. Subsequent global fitting of the data to the solutions of the coupled differential equations describing the fundamental mechanisms of induced fit and conformational selection, consistently support presence of the latter. To our knowledge, this is the first discrimination of conformational selection from induced fit for a mono-disperse CYP in a native-like membrane environment. In addition, 7-hydroxyflavone binds to CYP19A1 nanodiscs with comparable affinity to the substrates and induces an unusual spectral response likely attributable to hydrogen bonding to, rather than displacement of the heme-coordinated water molecule.

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.

Membrane-embedded substrate recognition by cytochrome P450 3A4

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

Spectroscopic studies of the cytochrome P450 reaction mechanisms

The cytochrome P450 monooxygenases (P450s) are thiolate heme proteins that can, often under physiological conditions, catalyze many distinct oxidative transformations on a wide variety of molecules, including relatively simple alkanes or fatty acids, as well as more complex compounds such as steroids and exogenous pollutants. They perform such impressive chemistry utilizing a sophisticated catalytic cycle that involves a series of consecutive chemical transformations of heme prosthetic group. Each of these steps provides a unique spectral signature that reflects changes in oxidation or spin states, deformation of the porphyrin ring or alteration of dioxygen moieties. For a long time, the focus of cytochrome P450 research was to understand the underlying reaction mechanism of each enzymatic step, with the biggest challenge being identification and characterization of the powerful oxidizing intermediates. Spectroscopic methods, such as electronic absorption (UV-Vis), electron paramagnetic resonance (EPR), nuclear magnetic resonance (NMR), electron nuclear double resonance (ENDOR), Mössbauer, X-ray absorption (XAS), and resonance Raman (rR), have been useful tools in providing multifaceted and detailed mechanistic insights into the biophysics and biochemistry of these fascinating enzymes. The combination of spectroscopic techniques with novel approaches, such as cryoreduction and Nanodisc technology, allowed for generation, trapping and characterizing long sought transient intermediates, a task that has been difficult to achieve using other methods. Results obtained from the UV-Vis, rR and EPR spectroscopies are the main focus of this review, while the remaining spectroscopic techniques are briefly summarized. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Drug-Drug Interactions between Atorvastatin and Dronedarone Mediated by Monomeric CYP3A4

Heterotropic interactions between atorvastatin (ARVS) and dronedarone (DND) have been deciphered using global analysis of the results of binding and turnover experiments for pure drugs and their mixtures. The in vivo presence of atorvastatin lactone (ARVL) was explicitly taken into account by using pure ARVL in analogous experiments. Both ARVL and ARVS inhibit DND binding and metabolism, while a significantly higher affinity of CYP3A4 for ARVL makes the latter the main modulator of activity (effector) in this system. Molecular dynamics simulations reveal significantly different modes of interactions of DND and ARVL with the substrate binding pocket and with a peripheral allosteric site. Interactions of both substrates with residues F213 and F219 at the allosteric site play a critical role in the communication of conformational changes induced by effector binding to productive binding of the substrate at the catalytic site.

Mathematical model of LsrR-binding and derepression in Escherichia coli K12

Quorum sensing (QS) enables bacterial communication and collective behavior in response to self-secreted signaling molecules. Unlocking its genetic regulation will provide insight towards understanding its influence on pathogenesis, formation of biofilms, and many other phenotypes. There are few datasets available that link QS-mediated gene expression to its regulatory components and even fewer mathematical models that incorporate known mechanistic detail. By integrating these data with annotated sequence information, mathematical inferences can be pieced together that shed light on regulatory structure. A first principles model, developed here for the E. coli QS system, builds on known mechanistic detail and is used to develop a working model of LuxS-regulated (Lsr) activity. That is, our model is meant to discriminate among hypothetical mechanisms governing lsr transcriptional regulation. Our simulations are in qualitative agreement with experimentally observed data. Importantly, our results point to the importance of transcriptional regulator, LsrR, cycling on genetic control. We also found several experimental observations in E. coli and homologous systems that were not explained by current mechanistic understanding. For example, by comparing simulations with reports of the integrating host factor in Aggrigatibacter actinomycetemcomitans, we conclude that additional transcriptional components are likely involved. An iterative process of simulation and experiment, therefore, is needed to inform new experiments and incorporate new model detail, the benefit of which will more rapidly validate mechanistic understanding.

Association of prenatal exposure to acetaminophen and coffee with childhood asthma

PURPOSE

Some studies have suggested that maternal acetaminophen use during pregnancy is associated with asthma in the offspring, and coffee consumption may modify the toxicity of acetaminophen. We aim to examine whether pregnancy maternal acetaminophen use increases the risk for offspring asthma, and whether such a potential association could be modified by maternal coffee consumption.

METHODS

We included 63,652 live-born singletons enrolled in the Danish National Birth Cohort. Maternal acetaminophen use and coffee consumption during pregnancy were assessed prospectively via the enrolment questionnaire and three computer-assisted telephone interviews. Asthma cases were identified by using the Danish National Patient Register and the Danish National Prescription Registry. We estimated the hazard ratios (HRs) for asthma according to prenatal acetaminophen and coffee exposure using Cox proportional hazards regression model.

RESULTS

After adjusting for potential confounders, acetaminophen use during pregnancy was associated with an increased risk of offspring asthma (HR = 1.16, 95% confidence interval (CI): 1.11-1.22). Coffee drinking during pregnancy was associated with a slightly decreased risk (HR = 0.94, 95%CI: 0.90-0.99). But there was no strong evidence of effect measure modification of acetaminophen use on offspring asthma by coffee consumption.

CONCLUSIONS

Acetaminophen use during pregnancy was associated with a modest increased risk for offspring asthma, which was not modified by coffee consumption.

Computational biotransformation profile of paracetamol catalyzed by cytochrome P450

The P450-catalyzed biotransformation of the analgesic drug paracetamol (PAR) is a long-debated topic, involving different mechanistic hypotheses as well as experimental evidence for the metabolites N-acetyl-p-benzoquinone imine (NAPQI), p-benzoquinone, acetamide, and 3-hydroxy-PAR. During the catalytic cycle of P450, a high-valent iron(IV)-oxo species known as Compound I (Cpd I) is formed as the ultimate oxidant, featuring two energetically close-lying ground states in the doublet (low-spin) and quartet (high-spin) spin states, respectively. In order to clarify the catalytic mechanism, a computational chemistry analysis has been undertaken for both the high- and low-spin routes, employing density functional theory (DFT) including PCM (polarized continuum-solvation model) that yields an approximate simulation of the bulk polarization exerted through the protein. The results demonstrate that hydrogen abstraction transfer (HAT) by the P450 oxidant Cpd I (FeO) is kinetically strongly preferred over the alternative pathways of an oxygen addition reaction (OAR) or two consecutive single-electron transfers (SET). Moreover, only the respective high-spin route yields N-acetyl-p-semiquinone imine (NAPSQI) as an intermediate that is converted to the electrophile N-acetyl-p-benzoquinone imine (NAPQI). By contrast, 3-hydroxy-PAR, acetamide, and p-benzoquinone as electrophilic and redox-active agent are formed predominantly in the low-spin state through reactions that do not involve NAPSQI. Thus, all experimentally observed PAR metabolites are in accord with an initial HAT from the phenolic oxygen, and NAPSQI should indeed be the precursor of NAPQI, both of which are generated only via the high-spin pathway.

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.

A cooperative mechanism of clotrimazoles in P450 revealed by the dissociation picture of clotrimazole from P450

The dissociation processes of clotrimazole (CLT) in several models are comparatively investigated by molecular dynamics simulations to explore the cooperative mechanism of clotrimazoles in P450. Our results suggest that when P450 only accommodates the active CLT (CLT1), CLT1 continually diffuses away from heme, and the partial BC loop (residues 73-88) and the extended FG loop (residues 173-186) first close and then open. When the enzyme binds to two CLT molecules, CLT1 basically keeps close to heme, and the partial BC loop and the extended FG loop move close to each other. Clearly, the effector CLT (CLT2) plays a cooperative role in the inhibition of CLT1 on P450. CLT2 restrains the dissociation of CLT1 first through direct π-π stacking interactions and then through the rearranged binding site induced by CLT2. The presence of CLT1 can help to stabilize the protein structure around CLT2 by interacting with M86, Q173, and M174.

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

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

Understanding a substrate's product regioselectivity in a family of enzymes: a case study of acetaminophen binding in cytochrome P450s

Product regioselectivity as influenced by molecular recognition is a key aspect of enzyme catalysis. We applied large-scale two-dimensional (2D) umbrella sampling (USP) simulations to characterize acetaminophen (APAP) binding in the active sites of the family of Cytochrome P450 (CYP) enzymes as a case study to show the different regioselectivity exhibited by a single substrate in comparative enzymes. Our results successfully explain the experimentally observed product regioselectivity for all five human CYPs included in this study, demonstrating that binding events play an important role in determining regioselectivity. In CYP2C9 and CYP3A4, weak interactions in an overall large active site cavity result in a fairly small binding free energy difference between APAP reactive binding states, consistent with experimental results that show little preference for resulting metabolites. In contrast, in CYP1A2 and CYP2E1, APAP is strongly restrained by a compact binding pocket, leading to a preferred binding conformation. The calculated binding equilibrium of APAP within the compact active site of CYP2A6 is able to predict the experimentally documented product ratios and is also applied to explain APAP regioselectivity in CYP1A2 and CYP2C9. APAP regioselectivity seems to be related to the selectivity for one binding conformation over another binding conformation as dictated by the size and shape of the active site. Additionally, unlike docking and molecular dynamics (MD), our free energy calculations successfully reproduced a unique APAP pose in CYP3A4 that had been reported experimentally, suggesting this approach is well suited to find the realistic binding pose and the lowest-energy starting structure for studying the chemical reaction step in the future.

π-π Stacking mediated drug-drug interactions in human CYP2E1

Because of having many low molecular mass substrates, CYP2E1 is of particular interests to the pharmaceutical industry. Many evidences showed that this enzyme can adopt multiple substrates to significantly reduce the oxidation rate of the substrates. The detailed mechanism for this observation is still unclear. In the current study, we employed GPU-accelerated molecular dynamics simulations to study the multiple-binding mode of human CYP2E1, with an aim of offering a mechanistic explanation for the unexplained multiple-substrate binding. Our results showed that Thr303 and Phe478 were key factors for the substrate recognition and multiple-substrate binding. The former can form a significant hydrogen bond to recognize and position the substrate in the productive binding orientation in the active site. The latter acted as a mediator for the substrate communications via π-π stacking interactions. In the multiple-binding mode, the aforementioned π-π stacking interactions formed by the aromatic rings of both substrates and Phe478 drove the first substrate far away from the catalytic center, orienting in an additional binding position and going against the substrate metabolism. All these findings could give atomic insights into the detailed mechanism for the multiple-substrate binding in human CYP2E1, providing useful information for the drug metabolism mechanism and personalized use of clinical drugs.

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.

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.

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.

Interactions of cytochrome P450s with their ligands

Cytochrome P450s (CYPs) are heme-containing monooxygenases that contribute to an enormous range of enzymatic function including biosynthetic and detoxification roles. This review summarizes recent studies concerning interactions of CYPs with ligands including substrates, inhibitors, and diatomic heme-ligating molecules. These studies highlight the complexity in the relationship between the heme spin state and active site occupancy, the roles of water in directing protein-ligand and ligand-heme interactions, and the details of interactions between heme and gaseous diatomic CYP ligands. Both kinetic and thermodynamic aspects of ligand binding are considered.

Spectroscopic features of cytochrome P450 reaction intermediates

Cytochromes P450 constitute a broad class of heme monooxygenase enzymes with more than 11,500 isozymes which have been identified in organisms from all biological kingdoms [1]. These enzymes are responsible for catalyzing dozens chemical oxidative transformations such as hydroxylation, epoxidation, N-demethylation, etc., with very broad range of substrates [2,3]. Historically these enzymes received their name from 'pigment 450' due to the unusual position of the Soret band in UV-vis absorption spectra of the reduced CO-saturated state [4,5]. Despite detailed biochemical characterization of many isozymes, as well as later discoveries of other 'P450-like heme enzymes' such as nitric oxide synthase and chloroperoxidase, the phenomenological term 'cytochrome P450' is still commonly used as indicating an essential spectroscopic feature of the functionally active protein which is now known to be due to the presence of a thiolate ligand to the heme iron [6]. Heme proteins with an imidazole ligand such as myoglobin and hemoglobin as well as an inactive form of P450 are characterized by Soret maxima at 420nm [7]. This historical perspective highlights the importance of spectroscopic methods for biochemical studies in general, and especially for heme enzymes, where the presence of the heme iron and porphyrin macrocycle provides rich variety of specific spectroscopic markers available for monitoring chemical transformations and transitions between active intermediates of catalytic cycle.

A risk-benefit assessment of paracetamol (acetaminophen) combined with caffeine

OBJECTIVE

To determine the risk: benefit of paracetamol combined with caffeine in the short-term management of acute pain conditions.

DESIGN

Database searches were conducted to identify double-blind trials comparing paracetamol/caffeine with paracetamol alone (benefit analysis) and any data pertaining to hepatotoxicity of paracetamol when combined with caffeine (risk analysis).

INTERVENTIONS

Paracetamol/caffeine (1,000 mg/130 mg) vs paracetamol (1,000 mg) alone.

OUTCOME MEASURES

Assessment of benefit has been derived by meta-analysis. Information on the pain condition and number of patients studied, dosing regimen, study design and analgesic outcome measures (total pain relief scores) was extracted and dichotomous outcomes were obtained by calculating the number of patients in each treatment group who achieved at least 50% of the maximum total pain relief score. Assessment of risk has been made by appraisal of the literature.

RESULTS

Eight studies from four papers provided sufficient quantitative data for satisfactory meta-analysis. The relative benefit (of achieving at least 50% pain relief) of paracetamol/caffeine vs paracetamol alone was 1.12 (95% Confidence Interval 1.05-1.19) across a number of acute pain states (dysmenorrhoea, headache, post-partum pain, and dental pain). Review of the effects of the combination of paracetamol and caffeine on the liver revealed no compelling data to suggest a clinically meaningful increase in hepatotoxicity with use of paracetamol/caffeine combinations.

CONCLUSIONS

Paracetamol/caffeine (1,000 mg/130 mg) is effective and safe for use in acute management of pain. The hepatotoxicity of overdoses of paracetamol results from its oxidative metabolism, caffeine does not produce any increase in oxidative metabolism of therapeutic concentrations of paracetamol.

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.

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

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

Alterations of chemotherapeutic pharmacokinetic profiles by drug-drug interactions

BACKGROUND

Drug interactions in oncology are common place and largely ignored as we tolerate high thresholds of 'toxic' drug responses in these patients. However, in the era of 'targeted' or seemingly 'less toxic' therapy, these interactions are more commonly flagged and contribute significantly towards poor 'quality of life' and medical fatalities.

OBJECTIVE

This review and opinion article focuses on alteration of chemotherapeutic pharmacokinetic profiles by drug interactions in the setting of polypharmacy. The assumption is that the drugs, with changes in their pharmacokinetics, will contribute towards changes in their pharmacodynamics.

METHODS

The examples cited for such drug-drug interactions are culled from published literature with an emphasis on those interactions that have been well characterized at the molecular level.

RESULTS

Although very few drug interaction studies have been performed on approved oncology based drugs, it is clear that drugs whose pharmacokinetics profiles are closely related to their pharmacodynamics will indeed result in clinically important drug interactions. Some newer mechanisms are described that involve interactions at the level of gene transcription, whereby, drug metabolism is significantly altered. However, for any given drug interaction, there does not seem to be a comprehensive model describing interactions.

CONCLUSIONS

Mechanisms based drug interactions are plentiful in oncology; however, there is an absolute lack of a comprehensive model that would predict drug-drug interactions.

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.