Implications of the allosteric kinetics of cytochrome P450s

Explicit Treatment of Non-Michaelis-Menten and Atypical Kinetics in Early Drug Discovery*

Biological systems are highly regulated. They are also highly resistant to sudden perturbations enabling them to maintain the dynamic equilibrium essential to sustain life. This robustness is conferred by regulatory mechanisms that influence the activity of enzymes/proteins within their cellular context to adapt to changing environmental conditions. However, the initial rules governing the study of enzyme kinetics were mostly tested and implemented for cytosolic enzyme systems that were easy to isolate and/or recombinantly express. Moreover, these enzymes lacked complex regulatory modalities. Now, with academic labs and pharmaceutical companies turning their attention to more-complex systems (for instance, multiprotein complexes, oligomeric assemblies, membrane proteins and post-translationally modified proteins), the initial axioms defined by Michaelis-Menten (MM) kinetics are rendered inadequate, and the development of a new kind of kinetic analysis to study these systems is required. This review strives to present an overview of enzyme kinetic mechanisms that are atypical and, oftentimes, do not conform to the classical MM kinetics. Further, it presents initial ideas on the design and analysis of experiments in early drug-discovery for such systems, to enable effective screening and characterisation of small-molecule inhibitors with desirable physiological outcomes.

The involvement of human organic anion transporting polypeptides (OATPs) in drug-herb/food interactions

Organic anion transporting polypeptides (OATPs) are important transporter proteins that are expressed at the plasma membrane of cells, where they mediate the influx of endogenous and exogenous substances including hormones, natural compounds and many clinically important drugs. OATP1A2, OATP2B1, OATP1B1 and OATP1B3 are the most important OATP isoforms and influence the pharmacokinetic performance of drugs. These OATPs are highly expressed in the kidney, intestine and liver, where they determine the distribution of drugs to these tissues. Herbal medicines are increasingly popular for their potential health benefits. Humans are also exposed to many natural compounds in fruits, vegetables and other food sources. In consequence, the consumption of herbal medicines or food sources together with a range of important drugs can result in drug-herb/food interactions via competing specific OATPs. Such interactions may lead to adverse clinical outcomes and unexpected toxicities of drug therapies. This review summarises the drug-herb/food interactions of drugs and chemicals that are present in herbal medicines and/or food in relation to human OATPs. This information can contribute to improving clinical outcomes and avoiding unexpected toxicities of drug therapies in patients.

Simulated Microgravity Altered the Metabolism of Loureirin B and the Expression of Major Cytochrome P450 in Liver of Rats

Loureirin B (LB) is the marker compound of dragon blood (DB), which exhibits great potentials in protecting astronauts' health against radiation and simulated microgravity (SM). Pharmacokinetics of LB is reported to be significantly altered by SM. Here, we investigated key metabolic features of LB in rat liver microsome (RLM) and the effects of SM on LB metabolism as well as on expression of major hepatic cytochrome P450 (CYP450) isoforms. Ten metabolites were tentatively identified based on fragmentation pathways using LC-MS/MS method and elimination kinetics of LB followed a typical Michaelis-Menten equation (V max was 1.32 μg/min/mg and K m was 13.33 μg/mL). CYP1A2, CYP2C11, CYP2D1, and CYP3A2 were involved in the metabolism of LB and the relative strength was: CYP3A2 > CYP2C11 > CYP2D1 > CYP1A2. Comparative studies suggested that elimination of LB in RLM was remarkably increased by 3-day and 14-day SM, and the generation of identified metabolites was affected as well. Additionally, 3-day and 14-day SM showed obvious regulatory effects on the expression of major CYP450 isoforms, which might contribute to the increased elimination of LB. The data provided supports for the application of DB as a protective agent and the reasonable use of current medications metabolized by hepatic CYP450 in space missions.

Use of bioconjugation with cytochrome P450 enzymes

Bioconjugation, defined as chemical modification of biomolecules, is widely employed in biological and biophysical studies. It can expand functional diversity and enable applications ranging from biocatalysis, biosensing and even therapy. This review summarizes how chemical modifications of cytochrome P450 enzymes (P450s or CYPs) have contributed to improving our understanding of these enzymes. Genetic modifications of P450s have also proven very useful but are not covered in this review. Bioconjugation has served to gain structural information and investigate the mechanism of P450s via photoaffinity labeling, mechanism-based inhibition (MBI) and fluorescence studies. P450 surface acetylation and protein cross-linking have contributed to the investigation of protein complexes formation involving P450 and its redox partner or other P450 enzymes. Finally, covalent immobilization on polymer surfaces or electrodes has benefited the areas of biocatalysis and biosensor design. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Genome-Wide Transcription and Functional Analyses Reveal Heterogeneous Molecular Mechanisms Driving Pyrethroids Resistance in the Major Malaria Vector Anopheles funestus Across Africa

Pyrethroid resistance in malaria vector, An. funestus is increasingly reported across Africa, threatening the sustainability of pyrethroid-based control interventions, including long lasting insecticidal nets (LLINs). Managing this problem requires understanding of the molecular basis of the resistance from different regions of the continent, to establish whether it is being driven by a single or independent selective events. Here, using a genome-wide transcription profiling of pyrethroid resistant populations from southern (Malawi), East (Uganda), and West Africa (Benin), we investigated the molecular basis of resistance, revealing strong differences between the different African regions. The duplicated cytochrome P450 genes (CYP6P9a and CYP6P9b) which were highly overexpressed in southern Africa are not the most upregulated in other regions, where other genes are more overexpressed, including GSTe2 in West (Benin) and CYP9K1 in East (Uganda). The lack of directional selection on both CYP6P9a and CYP6P9b in Uganda in contrast to southern Africa further supports the limited role of these genes outside southern Africa. However, other genes such as the P450 CYP9J11 are commonly overexpressed in all countries across Africa. Here, CYP9J11 is functionally characterized and shown to confer resistance to pyrethroids and moderate cross-resistance to carbamates (bendiocarb). The consistent overexpression of GSTe2 in Benin is coupled with a role of allelic variation at this gene as GAL4-UAS transgenic expression in Drosophila flies showed that the resistant 119F allele is highly efficient in conferring both DDT and permethrin resistance than the L119. The heterogeneity in the molecular basis of resistance and cross-resistance to insecticides in An. funestus populations throughout sub-Saharan African should be taken into account in designing resistance management strategies.

Allosteric Activation of Cytochrome P450 3A4 via Progesterone Bioconjugation

Human cytochrome P450 3A4 (CYP3A4) is responsible for the metabolism of the majority of drugs. As such, it is implicated in many adverse drug-drug and food-drug interactions, and is of significant interest to the pharmaceutical industry. This enzyme is known to simultaneously bind multiple ligands and display atypical enzyme kinetics, suggestive of allostery and cooperativity. As well, evidence of a postulated peripheral allosteric binding site has provoked debate around its significance and location. We report the use of bioconjugation to study the significance of substrate binding at the proposed allosteric site and its effect on CYP3A4 activity. CYP3A4 mutants were created and covalently modified with various small molecules including progesterone. The labeled mutants displayed enhanced kinetic stability and improved activity in testosterone and 7-benzyloxy-(4-trifluoromethyl)coumarin oxidation assays. Our work applies a new strategy to study cytochrome P450 allostery and supports the hypothesis that substrate binding at the postulated allosteric site of CYP3A4 may induce functional cooperativity.

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

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

Predicting Inhibitory Drug—Drug Interactions and Evaluating Drug Interaction Reports Using Inhibition Constants

OBJECTIVE:

To review the use of inhibitory constants (Ki) determined from in vitro experiments in the prediction of the significance of inhibitory drug—drug interactions (DDIs).

DATA SOURCES:

Searches of MEDLINE (1966—August 2004) and manual review of journals, conference proceedings, reference textbooks, and Web sites were performed using the key search terms cytochrome P450, drug—drug interaction, inhibition constant, and Ki.

STUDY SELECTION AND DATA EXTRACTION:

All articles identified from the data sources were evaluated, and information deemed relevant was included for this review.

DATA SYNTHESIS:

The cytochrome P450 isoenzymes factor prominently in the explanation of numerous DDIs. Although the regulation of these enzymes by one drug can affect the pharmacokinetics of other drugs, the consequences may not necessarily be significant either in terms of pharmacokinetic or clinical outcomes. Yet, many DDI monographs originate as unconfirmed case reports that implicate the influence of one drug on the CYP-mediated metabolism of another, and these often uncorroborated mechanisms can eventually become regarded as dogma. One consequence of this process is the overprediction of potentially important DDIs. The pharmaceutical industry, Food and Drug Administration, and pharmaceutical scientists have developed a strategy for predicting the significance of inhibitory DDIs at the earliest possible stages of drug development based on a new chemical entity's Ki value, determined in vitro.

CONCLUSIONS:

We suggest that the use of Ki values of drugs purported to behave as CYP inhibitors be incorporated in the assessment of case reports that ascribe DDIs to inhibition of metabolism of one drug by another.

Metabolism of BYZX in human liver microsomes and cytosol: identification of the metabolites and metabolic pathways of BYZX

BYZX, [(E)-2-(4-((diethylamino)methyl)benzylidene)-5,6-dimethoxy-2,3-dihydroinden-one], belongs to a series of novel acetylcholinesterase inhibitors and has been synthesized as a new chemical entity for the treatment of Alzheimer’s disease symptoms. When incubated with human liver microsomes (HLMs), BYZX was rapidly transformed into its metabolites M1, M2, and M3. The chemical structures of these metabolites were identified using liquid chromatography tandem mass spectrometry and nuclear magnetic resonance, which indicated that M1 was an N-desethylated and C = C hydrogenation metabolite of BYZX. M2 and M3 were 2 precursor metabolites, which resulted from the hydrogenation and desethylation of BYZX, respectively. Further studies with chemical inhibitors and human recombinant cytochrome P450s (CYPs), and correlation studies were performed. The results indicated that the N-desethylation of BYZX and M2 was mediated by CYP3A4 and CYP2C8. The reduced form of β-nicotinamide adenine dinucleotide 2′-phosphate was involved in the hydrogenation of BYZX and M3, and this reaction occurred in the HLMs and in the human liver cytosol. The hydrogenation reaction was not inhibited by any chemical inhibitors of CYPs, but it was significantly inhibited by some substrates of α,β-ketoalkene C = C reductases and their inhibitors such as benzylideneacetone, dicoumarol, and indomethacin. Our results suggest that α,β-ketoalkene C = C reductases may play a role in the hydrogenation reaction, but this issue requires further clarification.

CYP2E1 metabolism of styrene involves allostery

We are the first to report allosterism during styrene oxidation by recombinant CYP2E1 and human liver microsomes. At low styrene concentrations, oxidation is inefficient because of weak binding to CYP2E1 (K(s) = 830 μM). A second styrene molecule then binds CYP2E1 with higher affinity (K(ss) = 110 μM) and significantly improves oxidation to achieve a k(cat) of 6.3 nmol · min(-1) · nmol CYP2E1(-1). The transition between these metabolic cycles coincides with reported styrene concentrations in blood from exposed workers; thus, this CYP2E1 mechanism may be relevant in vivo. Scaled modeling of the in vitro-positive allosteric mechanism for styrene metabolism to its in vivo clearance led to significant deviations from the traditional model based on Michaelis-Menten kinetics. Low styrene levels were notably much less toxic than generally assumed. We interrogated the allosteric mechanism using the CYP2E1-specific inhibitor and drug 4-methylpyrazole, which we have shown binds two CYP2E1 sites. From the current studies, styrene was a positive allosteric effector on 4-methylpyrazole binding, based on a 10-fold increase in 4-methylpyrazole binding affinity from K(i) 0.51 to K(si) 0.043 μM. The inhibitor was a negative allosteric effector on styrene oxidation, because k(cat) decreased 6-fold to 0.98 nmol · min(-1) · nmol CYP2E1(-1). Consequently, mixtures of styrene and other molecules can induce allosteric effects on binding and metabolism by CYP2E1 and thus mitigate the efficiency of their metabolism and corresponding effects on human health. Taken together, our elucidation of mechanisms for these allosteric reactions provides a powerful tool for further investigating the complexities of CYP2E1 metabolism of drugs and pollutants.

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.

Does compound I vary significantly between isoforms of cytochrome P450?

The cytochrome P450 (CYP) enzymes are important in many areas, including pharmaceutical development. Subtle changes in the electronic structure of the active species, Compound I, have been postulated previously to account partly for the experimentally observed differences in reactivity between isoforms. Current predictive models of CYP metabolism typically assume an identical Compound I in all isoforms. Here we present a method to calculate the electronic structure and to estimate the Fe-O bond enthalpy of Compound I, and apply it to several human and bacterial CYP isoforms. Conformational flexibility is accounted for by sampling large numbers of structures from molecular dynamics simulations, which are subsequently optimized with density functional theory (B3LYP) based quantum mechanics/molecular mechanics. The observed differences in Compound I between human isoforms are small: They are generally smaller than the spread of values obtained for the same isoform starting from different initial structures. Hence, it is unlikely that the variation in activity between human isoforms is due to differences in the electronic structure of Compound I. A larger difference in electronic structure is observed between the human isoforms and P450(cam) and may be explained by the slightly different hydrogen-bonding environment surrounding the cysteinyl sulfur. The presence of substrate in the active site of all isoforms studied appears to cause a slight decrease in the Fe-O bond enthalpy, apparently due to displacement of water out of the active site, suggesting that Compound I is less stable in the presence of substrate.

Quantifying and predicting the promiscuity and isoform specificity of small-molecule cytochrome P450 inhibitors

Drug promiscuity (i.e., inhibition of multiple enzymes by a single compound) is increasingly recognized as an important pharmacological consideration in the drug development process. However, systematic studies of functional or physicochemical characteristics that correlate with drug promiscuity are handicapped by the lack of a good way of quantifying promiscuity. In this article, we present a new entropy-based index of drug promiscuity. We apply this index to two high-throughput data sets describing inhibition of cytochrome P450 isoforms by small-molecule drugs and drug candidates, and we demonstrate how drug promiscuity or specificity can be quantified. For these drug-metabolizing enzymes, we find that there is essentially no correlation between a drug's potency and specificity. We also present an index to quantify the susceptibilities of different enzymes to inhibition by diverse substrates. Finally, we use partial least-squares regression to successfully predict isoform specificity and promiscuity of small molecules, using a set of fingerprint-based descriptors.

Single-molecule fluorescence spectroscopy using phospholipid bilayer nanodiscs

Nanodiscs are a new class of model membranes that are being used to solubilize and study a range of integral membrane proteins and membrane-associated proteins. Unlike other model membranes, the Nanodisc bilayer is bounded by a scaffold protein coat that confers enhanced stability and a narrow particle size distribution. The bilayer diameter can be precisely controlled by changing the diameter of the protein coat. All these properties make Nanodiscs excellent model membranes for single-molecule fluorescence applications. In this chapter, we describe our work using Nanodiscs to apply total internal reflection fluorescence microscopy (TIRFM), fluorescence correlation spectroscopy (FCS), and Förster resonance energy transfer (FRET) to study the integral membrane protein cytochrome P450 3A4 and the peripheral membrane-binding proteins islet amyloid polypeptide (IAPP) and alpha-synuclein, respectively. The monodisperse size distribution of Nanodiscs enhances control over the oligomeric state of the membrane protein of interest, and facilitates accurate solution-based measurements as well. Nanodiscs also comprise an excellent system to stably immobilize integral membrane proteins in a bilayer without covalent modification, enabling a range of surface-based experiments where accurate localization of the protein of interest is required.

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.

Quantitative structure–activity relationships (QSARs) in CYP3A4 inhibitors: The importance of lipophilic character and hydrogen bonding

The results of Quantitative Structure-Activity Relationship (QSAR) analyses on three series of CYP3A4 inhibitors are reported for enzyme inhibition expressed as Ki values. These include a small series of structurally related statins and two larger groupings of structurally diverse compounds, some of which display competitive inhibition of CYP3A4 whereas others act via heme iron ligation. In all cases, however, it is apparent that there are lipophilicity relationships associated with CYP3A4 inhibitory activity in the total of 46 compounds investigated. This is evidenced by linear correlations between inhibition of CYP3A4 and the octanol-water partition coefficient (P value) when expressed logarithmically (ie., log P). In the case of the statins, however, the distribution coefficient (D) at pH 7 is used due to the effect of compound ionization. Conversion of equilibrium constants (ie. Ki and P) to the corresponding free energy changes (deltaG values) facilitates exploration of the likely intermolecular forces of interaction between the inhibitors and the active site region of CYP3A4. In this respect, there appears to be good agreement between QSAR analyses and molecular modelling of the CYP3A4 enzyme itself, and both are consistent with the known mechanisms of inhibition displayed.

Analysis of the cryptophycin P450 epoxidase reveals substrate tolerance and cooperativity

Cryptophycins are potent anticancer agents isolated from Nostoc sp. ATCC 53789 and Nostoc sp. GSV 224. The most potent natural cryptophycin analogues retain a beta-epoxide at the C2'-C3' position of the molecule. A P450 epoxidase encoded by c rpE recently identified from the cryptophycin gene cluster was shown to install this key functional group into cryptophycin-4 (Cr-4) to produce cryptophycin-2 (Cr-2) in a regio- and stereospecific manner. Here we report a detailed characterization of the CrpE epoxidase using an engineered maltose binding protein (MBP)-CrpE fusion. The substrate tolerance of the CrpE polypeptide was investigated with a series of structurally related cryptophycin analogues generated by chemoenzymatic synthesis. The enzyme specifically installed a beta-epoxide between C2' and C3' of cyclic cryptophycin analogues. The kcat/Km values of the enzyme were determined to provide further insights into the P450 epoxidase catalytic efficiency affected by substrate structural variation. Finally, binding analysis revealed cooperativity of MBP-CrpE toward natural and unnatural desepoxy cryptophycin substrates.

Thermodynamics of ligand binding to P450 2B4 and P450eryF studied by isothermal titration calorimetry

Structural plasticity and cooperativity in ligand recognition are two key aspects of the catalytic diversity of cytochrome P450 enzymes. As more mammalian P450 crystal structures have emerged, computational modeling has become a major tool to predict drug metabolism and interactions. There is a need for real solution thermodynamic data to support modeling and crystallographic observations. Using isothermal titration calorimetry (ITC) we successfully evaluated the conformational flexibility of P450 2B4 in binding imidazole inhibitors of different size and chemistry and dissected the stoichiometry and energetics of ligand binding allostery in P450eryF. Thermodynamic signatures obtained by ITC nicely correlated with structural and modeling results. Thus, ITC is a powerful tool to study structure-function relationships in P450s.

Energetics of heterotropic cooperativity between alpha-naphthoflavone and testosterone binding to CYP3A4

Cytochrome P450 3A4 (CYP3A4) is involved in the metabolism of a majority of drugs. Heterotropic cooperativity of drug binding to CYP3A4 was examined with the flavanoid, alpha-naphthoflavone (ANF) and the steroid, testosterone (TST). UV-vis and EPR spectroscopy of CYP3A4 show that ANF binding to CYP3A4 occurs with apparent negative cooperativity and that there are at least two binding sites: (1) a relatively tight spin-state insensitive binding site (CYP.ANF) and (2) a relatively low affinity spin-state sensitive binding site (CYP.ANF.ANF). Since binding to the spin-state insensitive binding site is considerably tighter for ANF than TST, the spin-state insensitive binding site could be occupied by ANF, while titrating TST at the other site(s). The spin-state insensitive binding site of ANF appears to compete with the spin-state insensitive binding site of TST. The formation of the spin-state insensitive CYP.ANF complex is strongly temperature dependent, when compared to the formation of the CYP.TST complex, suggesting that the formation of the CYP3A4.ANF complex leads to long-range conformational changes within the protein. When the CYP.ANF complex is titrated with TST, the formation of CYP.ANF.TST is favored by 3:1 over the formation of CYP.TST.TST, suggesting that there is an allosteric interaction between ANF and TST. A model of heterotropic cooperativity of CYP3A4 is presented, where the spin-state insensitive binding of ANF occurs at the same peripheral binding site of CYP3A4 as TST.

Chemiluminescent picture of diphenyleneiodonium-inhibited NADPH oxidase: a bimodal process and its logistic-exponential model-based description

A chemiluminescence (CL) study of diphenyleneiodonium-inhibited NADPH oxidase was performed on a cellular system containing neutrophils stimulated by phorbol myristate acetate, indicating a complex bimodal structure of CL processes corresponding to different stages of the inhibition. The complex structure of these processes was described by a superposition of two logistic-exponential (LE) models, characterizing these processes as bimodal ones. To determine the mechanistic foundation of the LE model-described processes, a generalized form of the second-order dynamic system of CL reactions, the solution to which corresponds to the LE model, was constructed. The diphenyleneiodonium effects on neutrophil NADPH oxidase were separated from the total bimodal CL of the whole measurement system by the use of difference CL processes. These difference processes were also found to be bimodal; thus, inhibitor-induced reduction of CL could be described by a second-order dynamic system. The rate constants and initial concentrations in this dynamic system were determined by the least squares method applied to numerical solutions approximating the difference processes. Using interrelations between the parameters of the dynamic system, cooperative effects in the inhibitor reactions with NADPH oxidase were found and described quantitatively. Other evidences of cooperativity were obtained from integral characteristics of the CL reduction process, i.e. dose-response and progress curves, determined by numerical integration of the LE models constituting the superposition. On this basis, it was also possible to detect a specific binding of the inhibitor to the enzyme. Finally, putative reaction mechanisms suggested by the model obtained were considered and compared with those known at present.

Mechanism of interactions of alpha-naphthoflavone with cytochrome P450 3A4 explored with an engineered enzyme bearing a fluorescent probe

Design of a partially cysteine-depleted C98S/C239S/C377S/C468A cytochrome P450 3A4 mutant designated CYP3A4(C58,C64) allowed site-directed incorporation of thiol-reactive fluorescent probes into alpha-helix A. The site of modification was identified as Cys-64 with the help of CYP3A4(C58) and CYP3A4(C64), each bearing only one accessible cysteine. Changes in the fluorescence of CYP3A4(C58,C64) labeled with 6-(bromoacetyl)-2-(dimethylamino)naphthalene (BADAN), 7-(diethylamino)-3-(4'-maleimidylphenyl)-4-methylcoumarin (CPM), or monobromobimane (mBBr) were used to study the interactions with bromocriptine (BCT), 1-pyrenebutanol (1-PB), testosterone (TST), and alpha-naphthoflavone (ANF). Of these substrates only ANF has a specific effect, causing a considerable decrease in fluorescence intensity of BADAN and CPM and increasing the fluorescence of mBBr. This ANF-binding event in the case of the BADAN-modified enzyme is characterized by an S50 of 18.2 +/- 0.7, compared with the value of 2.2 +/- 0.3 for the ANF-induced spin transition, thus revealing an additional low-affinity binding site. Studies of the effect of TST, 1-PB, and BCT on the interactions of ANF monitored by changes in fluorescence of CYP3A4(C58,C64)-BADAN or by the ANF-induced spin transition revealed no competition by these substrates. Investigation of the kinetics of fluorescence increase upon H2O2-dependent heme depletion suggests that labeled CYP3A4(C58,C64) is represented by two conformers, one of which has the fluorescence of the BADAN and CPM labels completely quenched, presumably by photoinduced electron transfer from the neighboring Trp-72 and/or Tyr-68 residues. The binding of ANF to the newly discovered binding site appears to affect the interactions of the label with the above residue(s), thus modulating the fraction of the fluorescent conformer.

Dissecting the Thermodynamics and Cooperativity of Ligand Binding in Cytochrome P450eryF

Conformational flexibility and cooperativity in ligand recognition are two key aspects of the catalytic diversity of cytochrome P450 enzymes. In this study, we dissect the ligand binding stoichiometry and energetics of the soluble bacterial P450eryF by isothermal titration calorimetry (ITC) using three allosteric and two non-allosteric ligands of diverse chemistry. Complementary spectral binding studies and sequential, two-ligand docking simulations were performed to help assign the binding sites. Binding of 4-phenylpyridine (4-PP) or 4-(4-chlorophenyl)imidazole (4-CPI) showed 1:1 stoichiometry in ITC, consistent with the lack of cooperativity observed in spectral binding studies. The larger ligands 9-aminophenanthrene (9-AP), 1-pyrenebutanol (1-PB), and alpha-naphthoflavone (ANF) show cooperative spectral binding and yielded 2:1 stoichiometry. The associated thermodynamic parameters for the sites were calculated using a sequential binding mechanism. The binding constant (KD) for the first site was almost two times lower than that of the second site for all three compounds. Ligand binding at site 1 was entropically favored, whereas binding at site 2 was weak and entropically unfavorable. Simulations showed that two molecules of 9-AP, ANF or 1-PB can be adequately docked to two individual sub-sites within a large binding pocket. The absence of hydrophobic tethering and ligand stacking are consistent with the single low affinity binding site observed for 4-CPI and 4-PP. Competitive binding studies with P450eryF preloaded with either 1-PB or ANF showed a decrease in the affinities for 9-AP at both the sites, with large entropy-enthalpy compensation, indicating the ability of the binding pocket to accommodate two ligands of diverse chemistry and enable cooperativity.

Capillary electrophoretic investigation of the enantioselective metabolism of propafenone by human cytochrome P-450 SUPERSOMES: Evidence for atypical kinetics by CYP2D6 and CYP3A4

An enantioselective CE method was used to identify the ability of CYP450 enzymes and their stereoselectivity in catalyzing the transformation of propafenone (PPF) to 5-hydroxy-propafenone (5OH-PPF) and N-despropyl-propafenone (NOR-PPF). Using in vitro incubations with single CYP450 enzymes (SUPERSOMES), 5OH-PPF is shown to be selectively produced by CYP2D6 and N-dealkylation is demonstrated to be mediated by CYP2D6, CYP3A4, CYP1A2, and CYP1A1. For the elucidation of kinetic aspects of the metabolism with CYP2D6 and CYP3A4, incubations with individual PPF enantiomers and racemic PPF were investigated. With the exception of the dealkylation in presence of R-PPF only, which can be described by the Michaelis-Menten model, all CYP2D6-induced reactions were found to follow autoactivation kinetics. For CYP3A4, all NOR-PPF enantiomer formation rates as function of PPF enantiomer concentration were determined to follow substrate inhibition kinetics. The formation of NOR-PPF by the different enzymes is stereoselective and is reduced significantly when racemic PPF is incubated. Clearance values obtained for CYP3A4 dealkylation are stereoselective whereas those of CYP2D6 hydroxylation are not. This paper reports the first investigation of the PPF hydroxylation and dealkylation kinetics by the CYP2D6 enzyme and represents the first report in which enantioselective CE data provide the complete in vitro kinetics of metabolic steps of a drug.

Heterotropic and homotropic cooperativity by a drug-metabolising mutant of cytochrome P450 BM3

Recently, we described a triple mutant of the bacterial cytochrome P450 BM3 as the first mutant with affinity for drug-like compounds. In this paper, we show that this mutant, but not wild-type BM3, is able to metabolise testosterone and several drug-like molecules such as amodiaquine, dextromethorphan, acetaminophen, and 3,4-methylenedioxymethylamphetamine that are known substrates of human P450s. Interestingly, the metabolism of 3,4-methylenedioxymethylamphetamine and acetaminophen could be stimulated up to 70-fold by the addition of caffeine, a known activator of rat P450 3A2. With testosterone metabolism, homotropic cooperativity was observed. This shows that heterotropic and homotropic cooperativity, known to occur in the P450 3A family, can also take place in BM3. BM3 therefore can be used as a model system to study atypical kinetics in mammalian P450s. Second, this study shows that BM3 can be engineered to a drug-metabolising enzyme, making it a promising candidate to use as biocatalyst in drug discovery and synthesis.

Resolution of multiple substrate binding sites in cytochrome P450 3A4: the stoichiometry of the enzyme-substrate complexes probed by FRET and Job's titration

To explore the mechanism of homotropic cooperativity in human cytochrome P450 3A4 (CYP3A4) we studied the interactions of the enzyme with 1-pyrenebutanol (1-PB), 1-pyrenemethylamine (PMA), and bromocriptine by FRET from the substrate fluorophore to the heme, and by absorbance spectroscopy. These approaches combined with an innovative setup of titration-by-dilution and continuous variation (Job's titration) experiments allowed us to probe the relationship between substrate binding and the subsequent spin transition caused by 1-PB or bromocriptine or the type-II spectral changes caused by PMA. The 1-PB-induced spin shift in CYP3A4 reveals prominent homotropic cooperativity, which is characterized by a Hill coefficient of 1.8 +/- 0.3 (S50 = 8.0 +/- 1.1 microM). In contrast, the interactions of CYP3A4 with bromocriptine or PMA reveal no cooperativity, exhibiting KD values of 0.31 +/- 0.08 microM and 7.1 +/- 2.3 microM, respectively. The binding of all three substrates monitored by FRET in titration-by-dilution experiments at an enzyme:substrate ratio of 1 reveals a simple bimolecular interaction with KD values of 0.16 +/- 0.09, 4.8 +/- 1.4, and 0.18 +/- 0.09 microM for 1-PB, PMA, and bromocriptine, respectively. Correspondingly, Job's titration experiments showed that the 1-PB-induced spin shift reflects the formation of a complex of the enzyme with two substrate molecules, while bromocriptine and PMA exhibit 1:1 binding stoichiometry. Combining the results of Job's titrations with the value of KD obtained in our FRET experiments, we demonstrate that the interactions of CYP3A4 with 1-PB obey a sequential binding mechanism, where the spin transition is triggered by the binding of 1-PB to the low-affinity site, which becomes possible only upon saturation of the high-affinity site.

Matrix metalloproteinase inhibitor RECK expression in canine tumors

Matrix metalloproteinases (MMPs) selectively degrade the extracellular matrix, and they have been reported to play an important role in tumor invasion, metastasis and angiogenesis. These enzymes are closely related to tumor malignancy and patient survival time. Recently, reversion-inducing cysteine-rich protein with Kazal motifs (RECK) gene was identified as an endogenous membrane-anchored MMP inhibitor. The down-regulation of RECK has been implicated in tumor progression. In this study, the expression levels of the RECK messenger ribonucleic acid (mRNA) in various spontaneously developed canine tumors were investigated by using quantitative reverse transcriptase-polymerase chain reaction (RT-PCR), and the correlation between RECK and clinicopathological factors, as well as MMP-9 expression were analyzed. The median age of 36 dogs investigated in this study was 9 years old (range, 1-15 years old). Quantitative RT-PCR could detect low levels of expression of RECK mRNA in the tumor samples. The expression levels of RECK mRNA in some tumor tissue samples were significantly lower than those in normal tissue samples. No significant associations of RECK with clinicopathological factors were observed. Using the Mann-Whitney U test, the expression level of the MMP-9 mRNA was observed to be significantly correlated to RECK expression (p<0.05).

IN VITRO METABOLISM OF THE CALMODULIN ANTAGONIST DY-9760e (3-[2-[4-(3-CHLORO-2-METHYLPHENYL)-1-PIPERAZINYL]ETHYL]-5,6-DIMETHOXY-1-(4-IMIDAZOLYLMETHYL)-1H-INDAZOLE DIHYDROCHLORIDE 3.5 HYDRATE) BY HUMAN LIVER MICROSOMES: INVOLVEMENT OF CYTOCHROMES P450 IN ATYPICAL KINETICS AND POTENTIAL DRUG INTERACTIONS

Human cytochrome P450 (P450) isozyme(s) responsible for metabolism of the calmodulin antagonist 3-[2-[4-(3-chloro-2-methylphenyl)-1-piperazinyl]ethyl]-5,6-dimethoxy-1-(4-imidazolylmethyl)-1H-indazole dihydrochloride 3.5 hydrate (DY-9760e) and kinetic profiles for formation of its six primary metabolites [M3, M5, M6, M7, M8, and DY-9836 (3-[2-[4-(3-chloro-2-methylphenyl)piperazinyl]ethyl]-5,6-dimethoxyindazole)] were identified using human liver microsomes and recombinant P450 enzymes. In vitro experiments, including an immunoinhibition study, correlation analysis, and reactions with recombinant P450 enzymes, revealed that CYP3A4 is the primary P450 isozyme responsible for the formation of the DY-9760e metabolites, except for M5, which is metabolized by CYP2C9. Additionally, at clinically relevant concentrations, CYP2C8 and 2C19 make some contribution to the formation of M3 and M5, respectively. The formation rates of DY-9760e metabolites except for M8 by human liver microsomes are not consistent with a Michaelis-Menten kinetics model, but are better described by a substrate inhibition model. In contrast, the enzyme kinetics for all metabolites formed by recombinant CYP3A4 can be described by an autoactivation model or a mixed model of autoactivation and biphasic kinetics. Inhibition of human P450 enzymes by DY-9760e in human liver microsomes was also investigated. DY-9760e is a very potent competitive inhibitor of CYP2C8, 2C9 and 2D6 (Ki 0.25-1.7 microM), a mixed competitive and noncompetitive inhibitor of CYP2C19 (Ki 2.4 microM) and a moderate inhibitor of CYP1A2 and 3A4 (Ki 11.4-20.1 microM), suggesting a high possibility for human drug-drug interaction.

Fromin vivotoin vitro/in silicoADME: progress and challenges

High-throughput screening technologies in biological sciences of large libraries of compounds obtained via combinatorial or parallel chemistry approaches, as well as the application of design rules for drug-likeness, have resulted in more hits to be evaluated with respect to their ADME or drug metabolism and pharmacokinetic properties. The traditional in vivo methods using preclinical species, such as rat, dog or monkey, are no longer sufficient to cope with this demand. This editorial discusses the changes towards medium- to high-throughput in vitro and in silico ADME screening. In addition, much more attention is now put on early safety and risk assessment of promising lead series and potential clinical candidates.

Multivariate modeling of cytochrome P450 3A4 inhibition

In the early phases of current pharmaceutical research projects, huge numbers of compounds are tested on their biological activity with respect to a certain target by experimental or virtual screening campaigns. To reduce the attrition rate in later stages of a project, other relevant properties such as physicochemical and ADMET (absorption, distribution, metabolism, excretion, toxicity) properties should be assessed as early as possible in lead discovery and optimization. The present study describes the development of in silico models to predict the inhibition of human cytochrome P450 3A4 (CYP3A4) from calculated molecular descriptors. The models were trained and validated using a set of 967 structural diverse drug-like research compounds with an experimentally determined CYP3A4 inhibition potency (IC50 value) which was carefully split into a training and a test set. For classification models, the data sets were further subdivided into strong, medium, and weak inhibitors. Different descriptor sets were used to cover various aspects of molecular properties, including properties derived from the 2D structure, the interaction of the molecule with its environment, and properties derived from quantum-mechanical calculations. The descriptors were related to the CYP3A4 inhibition potency by multivariate data analysis methods such as partial least-squares projection to latent structures (PLS), PLS discriminant analysis (PLS-DA), and soft independent class modeling (SIMCA). The squared correlation between experimental and predicted IC50 values of the previously unseen test set compounds was Qext2=0.6 for the best PLS models, corresponding to a root mean squared error (RMSE) of RMSE=0.45 (logarithm of IC50). The best PLS-DA models were able to correctly classify more than 60% of the test set compounds, whereas almost no strong inhibitors were wrongly classified as weak inhibitors and vice versa. Furthermore, relevant molecular properties were identified which are closely related to the CYP3A4 inhibition potency of a compound. The results presented here are very encouraging since our models could, for instance, serve to flag problematic compounds or to guide further synthesis efforts.