Allosteric activation of cytochrome P450 3A4 by α-naphthoflavone: branch point regulation revealed by isotope dilution analysis

Cytochrome P450 3A4 (CYP3A4) is the dominant xenobiotic metabolizing CYP. Despite great interest in CYP enzymology, two in vitro aspects of CYP3A4 catalysis are still not well understood, namely, sequential metabolism and allosteric activation. We have therefore investigated such a system in which both phenomena are present. Here we report that the sequential metabolism of Nile Red (NR) is accelerated by the heterotropic allosteric effector α-naphthoflavone (ANF). ANF increases the rates of formation for NR metabolites M1 and M2 and also perturbs the metabolite ratio in favor of M2. Thus, ANF has as an allosteric effect on a kinetic branch point. Co-incubating deuterium-labeled NR and unlabeled M1, we show that ANF increases k(cat)/k(off) ~1.8-fold in favor of the k(cat) of M2 production. Steady-state metabolic experiments are analyzed using a kinetic model in which the enzyme and substrates are not in rapid equilibrium, and this distinction allows for the estimation of rates of catalysis for the formation of both the primary (M1) and secondary (M2) products, as well as the partitioning of enzyme between these states. These results are compared with those of earlier spectroscopic investigations of NR and ANF cooperativity, and a mechanism of ANF heteroactivation is presented that involves effects on substrate off rate and coupling efficiency.     

Influence of major structural features of tocopherols and tocotrienols on their omega-oxidation by tocopherol-omega-hydroxylase

Human cytochrome P450 4F2 (CYP4F2) catalyzes the initial omega-hydroxylation reaction in the metabolism of tocopherols and tocotrienols to carboxychromanols and is, to date, the only enzyme shown to metabolize vitamin E. The objective of this study was to characterize this activity, particularly the influence of key features of tocochromanol substrate structure. The influence of the number and positions of methyl groups on the chromanol ring, and of stereochemistry and saturation of the side chain, were explored using HepG2 cultures and microsomal reaction systems. Human liver microsomes and microsomes selectively expressing recombinant human CYP4F2 exhibited substrate activity patterns similar to those of HepG2 cells. Although activity was strongly associated with substrate accumulation by cells or microsomes, substantial differences in specific activities between substrates remained under conditions of similar microsomal membrane substrate concentration. Methylation at C5 of the chromanol ring was associated with markedly low activity. Tocotrienols exhibited much higher Vmax values than their tocopherol counterparts. Side chain stereochemistry had no effect on omega-hydroxylation of alpha-tocopherol (alpha-TOH) by any system. Kinetic analysis of microsomal CYP4F2 activity revealed Michaelis-Menten kinetics for alpha-TOH but allosteric cooperativity for other vitamers, especially tocotrienols. Additionally, alpha-TOH was a positive effector of omega-hydroxylation of other vitamers. These results indicate that CYP4F2-mediated tocopherol-omega-hydroxylation is a central feature underlying the different biological half-lives, and therefore biopotencies, of the tocopherols and tocotrienols.     

Allosteric phenomena in cytochrome P450-catalyzed monooxygenations.

Allosteric regulation of monooxygenase activity is shown to occur with diverse cytochrome P450 isoforms and is characterized by kinetic patterns deviating from the Michaelis-Menten model. Homotropic and heterotropic phenomena are encountered in both substrate activation and productive coupling of the electron donors NADPH-cytochrome P450 reductase and cytochrome b5, and the lipid environment of the system also appears to play a role as an effector. Circumstantial analysis reveals the components of the electron transfer chain to be mutually beneficial in interactions with each other depending on the substrate used and type of cytochrome P450 operative. It is noteworthy that association of diatomic gaseous ligands may be amenable to allosteric regulation as well. Thus, dioxygen binding to cytochrome P450 displays nonhyperbolic kinetic profiles in the presence of certain substrates; the latter, together with redox proteins such as cytochrome b5, can exert efficient control of the abortive breakdown of the oxyferrous intermediates formed. Similarly, substrates may modulate the structural features of the access channel for solutes such as carbon monoxide in specific cytochrome P450 isozymes to either facilitate or impair ligand diffusion to the heme iron. The in vivo importance of allosteric regulation of enzyme activity is discussed in detail.     

α-氯丙酸脱卤酶发酵动力学模型

对α-氯丙酸脱卤酶发酵动力学进行了研究。基于Logistic方程和Luedeking-Piret方程,得到了描述Pseudomonas W20菌发酵过程菌体生长、α-氯丙酸脱卤酶生成及基质消耗的动力学数学模型和模型参数,对试验数据与模型进行了验证比较,模型计算值与试验结果拟合良好,平均相对误差大部分小于10%;对脱卤酶反应动力学进行了研究,结果表明脱卤酶的脱卤反应基本符合米氏方程,并求得最大反应速率V_(max)=1.11×10~(-5)mol/(g·min),表观米氏常数K_m=3.72×10~(-3)mol/L。     

Inhibition of CYP2B4 by the mechanism-based inhibitor 2-ethynylnaphthalene: inhibitory potential of 2EN is dependent on the size of the substrate

2-Ethynylnaphthalene (2EN) is a mechanism-based inhibitor of CYP2B4 with two components to the inhibition, (1) enzyme inactivation, which requires covalent binding of the 2EN metabolite, and (2) reversible inhibition by 2EN itself. Both inhibitory components were examined using several different CYP2B4 substrates. Preincubation of CYP2B4 with 2EN led to a time-dependent inactivation of each of the CYP2B4-dependent activities examined; however, the ability of 2EN to reversibly inhibit CYP2B4 depended on the substrate employed, which is inconsistent with classical inhibition patterns. The degree 2EN's reversible inhibition was shown not to correlate with the substrate affinity for the active site, but with parameters related to the molecular size of the substrate. The results are consistent with 2EN and the smaller substrates simultaneously fitting in the CYP2B4 active site, leading to very little inhibition. Larger substrates exhibited greater degrees of inhibition because of their inability to co-bind with inhibitor in the active site.     

Measurement of biodegradability parameters for single unsubstituted and methylated polycyclic aromatic hydrocarbons in liquid bacterial suspensions

Substrate depletion experiments were conducted to characterize aerobic biodegradation of 20 single polycyclic aromatic hydrocarbons (PAHs) by induced Sphingomonas paucimobilis strain EPA505 in liquid suspensions. PAHs consisted of low molecular weight, unsubstituted, and methyl-substituted homologs. A material balance equation containing the Andrews kinetic model, an extension of the Monod model accounting for substrate inhibition, was numerically fitted to batch depletion data to estimate extant kinetic parameters including the maximal specific uptake rates, q(max), the affinity coefficients, K(S), and the substrate inhibition coefficients, K(I). Strain EPA505 degraded all PAHs tested. Applied kinetic models adequately simulated experimental data. A cell proliferation assay involving reduction of the tetrazolium dye WST-1 was used to evaluate the ability of strain EPA505 to utilize individual PAHs as sole energy and carbon sources. Of the 22 PAHs tested, 9 supported bacterial growth. Evaluation of the biokinetic data showed that q(max) correlated highly with transmembrane flux as theoretically estimated by a diffusion model, pointing to transmembrane transport as a potential rate-determining process. The biodegradability data generated in this study is essential for the development of quantitative structure-activity relationships (QSARs) for biodegradability and for modeling biodegradation of simple PAH mixtures.     

Characterizing the Effect of Cytochrome P450 (CYP) 2C8, CYP2C9, and CYP2D6 Genetic Polymorphisms on Stereoselective N‐demethylation of Fluoxetine

Fluoxetine (FLX) is one of the most widely prescribed selective serotonin reuptake inhibitors. Although FLX is used as racemate in the clinic, the clinical pharmacokinetics of FLX and its N‐demethylation metabolite norfluoxetine (NFLX) show obvious cytochrome P450 (CYP) polymorphism dependency and exhibit marked stereoselectivity. However, the kinetic profiles of CYP variants to FLX remain unclear. In the present study, some variants of human CYP2C8, CYP2C9, and CYP2D6 were first expressed in insect cells, and their catalytic roles with respect to FLX enantiomers were then investigated. CYP2C8.4 and CYP2C9.10 showed significantly lower activity and CYP2C8.3 showed significantly higher activity toward both R‐ and S‐FLX compared with the wildtype, while CYP2C9.3, CYP2C9.13, and CYP2C9.16 showed significantly lower activity only toward R‐FLX. Five CYP2C9 variants and CYP2D6.1 exhibited significantly stereoselective kinetic profiles prior to R‐FLX, and CYP2C8.3 showed a slight stereoselectivity. Interestingly, obvious substrate inhibition was observed in the CYP2C9 wildtype and its three variants only in the case of R‐FLX. Together, these findings suggest that CYP2C9 and CYP2D6 polymorphism may play an important role in the clearance of FLX and also in the stereoselective kinetic profiles of FLX enantiomers. Chirality 26:166–173, 2014. © 2014 Wiley Periodicals, Inc.     

Effect of trans-resveratrol on 7-benzyloxy-4-trifluoromethylcoumarin O-dealkylation catalyzed by human recombinant CYP3A4 and CYP3A5

Red wine concentrate has been reported to inhibit the catalytic activity of human recombinant cytochrome P450 (CYP) 3A4. Wine contains many polyphenolic compounds, including trans-resveratrol, which is also available commercially as a nutraceutical product. In the present study, we examined the in vitro effect of trans-resveratrol on human CYP3A catalytic activity by employing recombinant CYP3A4 and CYP3A5 as model enzymes and 7-benzyloxy-4-trifluoromethylcoumarin (BFC) as a CYP3A substrate. Trans-resveratrol inhibited BFC O-dealkylation catalyzed by CYP3A4 and CYP3A5 in a concentration-dependent manner. In each case, the inhibition was noncompetitive, as determined by Lineweaver-Burk and Dixon plots of the enzyme kinetic data. The apparent Ki values (mean +/- SEM) for the inhibition by trans-resveratrol of BFC O-dealkylation catalyzed by CYP3A4 and CYP3A5 were 10.2+/-1.1 microM and 14.7+/-0.3 microM, respectively. Preincubation of trans-resveratrol with NADPH and CYP3A4 or CYP3A5 for 10 or 15 min prior to initiation of substrate oxidation did not enhance the inhibitory effect, suggesting that this compound was not a mechanism-based inactivator of CYP3A4 or CYP3A5 when BFC was used as the substrate. Overall, our study provides the first demonstration that trans-resveratrol inhibits, in vitro, a substrate oxidation reaction catalyzed by human recombinant CYP3A4 and CYP3A5.     

Alternative kinetic laws to describe the turnover of the microbial biomass

Many models that describe the turnover of the microbial biomass in soil use either first order kinetics where the rate of turnover is directly proportional to the microbial mass, or a variant of the Michaelis-Menten law that describes enzyme kinetics. To account for the different rates of microbial turnover observed at different times after the addition of substrate, some authors have suggested the existence of more than one pool of biomass. Each pool obeys the same kinetic law but with a different rate. In other experiments a disproportionately large increase in the turnover of native organisms has been observed relative to the amount of fresh substrate added. A change in the kinetic law describing the turnover of organisms can account for these observations and yet retain the simplicity of a single pool of micro-organisms. However where multiple pools of organisms are justified a mixed kinetic law with both first and second order terms may be more appropriate; in other words one pool of micro-organisms but two rate constants. The advantage of retaining a single pool of microbial biomass is that models may more readily be constructed in relation to the routine measurements of total microbial mass.     

Multistep phosphorylation systems: tunable components of biological signaling circuits

Multisite phosphorylation of proteins is a powerful signal processing mechanism that plays crucial roles in cell division and differentiation as well as in disease. We recently demonstrated a novel phenomenon in cell cycle regulation by showing that cyclin-dependent kinase–dependent multisite phosphorylation of a crucial substrate is performed sequentially in the N-to-C terminal direction along the disordered protein. The process is controlled by key parameters, including the distance between phosphorylation sites, the distribution of serines and threonines in sites, and the position of docking motifs. According to our model, linear patterns of phosphorylation along disordered protein segments determine the signal-response function of a multisite phosphorylation switch. Here we discuss the general advantages and engineering principles of multisite phosphorylation networks as processors of kinase signals. We also address the idea of using the mechanistic logic of linear multisite phosphorylation networks to design circuits for synthetic biology applications.     

Substrate inhibition during bio-filtration of TCE using diazotrophic bacterial community

The kinetics of biodegradation of TCE in the biofilter packed with wood charcoal and inoculated with diazotrophic bacterial community had been investigated. Use of Michaelis-Menten type model showed that substrate inhibition was present in the system. The kinetic model proposed by Edwards (1970) was used to calculate kinetic parameters-maximum elimination capacity (EC(max)), substrate constant (K(s)), and inhibition constant (K(I)). The model fitted well with the experimental data and the EC(max) was found to be in the range of 10.8-6.1 g/m(3) h. The K(s) values depended upon substrate concentration and ranged from 0.024 to 0.043 g/m(3) indicating the high affinity of diazotrophs for TCE. The K(I) values were low and nearly constant (0.011-0.015 g/m(3)) indicating a moderate substrate inhibition.     

Multi-random binding molecular interactions: a general kinetic model

A kinetic model is suggested to account for the interactions of several ligands with a target whose molecule possesses several independent equivalent receptor sites for each ligand (multiligand multisite model). To analyse the problem, we shall derive solutions for three elementary situations: (a) interactions of a ligand with a mono-receptor site target molecule (monosite model); (b) interactions of several ligands with a target whose molecule possesses one receptor site for each ligand involved (multiligand model); (c) interactions of a ligand with a target whose molecule possesses several receptor sites of the same kind for this ligand (multisite model). Throughout this study, every ligand molecule is assumed to offer one binding site to the target. The main implications of the corresponding analytical solutions are discussed from a molecular point of view. The results cover a great many well-known aspects of the molecular interactions in various fields such as enzymology, endocrinology, radio-immunology and saturation analysis. As suggested by the inhibition patterns obtained, this model may therefore provide a new point of view to interpret the relevant phenomena. Furthermore, a kinetic approach to the generalized mass action law can be deduced from this model, and experimental conditions in which the isotopic dilution law applies are examined.     

Progress Curves Analysis as an Alternative for Exploration of Activation-inhibition Phenomena in Cholinesterases

The kinetic behaviour of insect acetylcholinesterases deviates from the Michaelis-Menten pattern. These deviations are known as activation or inhibition at various substrate concentrations and can be more or less observable depending on mutations around the active site of the enzyme. Most kinetic studies on these enzymes still rely on initial rate measurements. It is demonstrated here that according to this method one of the deviations can be overlooked. We attempt to point out that in such cases a detailed step-by-step progress curves analysis is successful. The study is focused on two different methods of analysing progress curves: (i) the first one is based on an integrated initial rate equation which can sufficiently fit truncated progress curves under corresponding conditions; and (ii) the other one precludes the algebraic formulae, but uses numerical integration for searching a non analytical solution of ordinary differential equations describing a kinetic model. All methods are tested on three different acetylcholinesterase mutants from Drosophila melanogaster. The results indicate that kinetic parameters for the E107K mutant with highly expressive activation and inhibition can be well evaluated applying any analysis method. It is quite different for E107W and E107Y mutants where latent activation is present, but discovered only using one or the other progress curves analysis methods.     

Thymidine kinase isoenzymes in human acute monocytic leukemia

Summary Two thymidine kinase isoenzymes, TK 3 and TK 4, from mononuclear leucocytes from a patient with acute monocytic leukemia, were purified and characterized in regard to the molecular weights and kinetic properties.The molecular weights of TK 3 and TK 4 were 60 000 and 45 000, respectively. In the presence of 2 mM ATP, the molecular weight of TK 3 increased to 200 000, whereas the molecular weight of TK 4 was unchanged.Studies of the kinetic properties showed clear differences between TK 3 and TK 4. With thymidine as substrate, TK 3 showed biphasic kinetics with a K_m of 22 µM, and TK 4 showed Michaelis-Menten kinetics with a K_m of 0.33 µM With ATP as substrate, TK 3 showed Michaelis-Menten kinetics with a K_m of 100 µM, and TK 4 showed biphasic kinetics with a Km of 3.5 µM. With dTTP as inhibitor, TK 3 showed cooperative inhibition kinetics, and TK 4 showed non-cooperative competitive inhibition kinetics. The dTTP concentration at 50% inhibition was 75 µM for TK 3 but 380 µM for TK 4.Comparison of the molecular weights and the kinetic properties of TK 3 and TK 4 with the corresponding data previously obtained for TK 1 and TK 2 from normal human lymphocytes indicate the existence of four thymidine kinase isoenzymes in human leucocytes.     

Kinetics of multilayer immobilized enzyme-filter reactors: Behavior of urease-filter reactors in different buffers

A sequential model which describes the kinetic behavior of multilayer immobilized-enzyme filter reactors for systems obeying the Michaelis-Menten scheme and various types of inhibition by product and by substrate is proposed. The model was varified experimentaly using a single enzyme-filter-substrate system-urease bound to nylon-filter disks acting on urea-in different buffers to obtain the various inhibition patterns.     

Role of substrate inhibition kinetics in enzymatic chemical oscillations.

Two chemical kinetic models are investigated using standard nonlinear dynamics techniques to determine the conditions under which substrate inhibition kinetics can lead to oscillations. The first model is a classical substrate inhibition scheme based on Michaelis-Menten kinetics and involves a single substrate. Only when this reaction takes place in a flow reactor (i.e., both substrate and product are taken to follow reversible flow terms) are oscillations observed; however, the range of parameter values over which such oscillations occur is so narrow it is experimentally unobservable. A second model based on a general mechanism applied to the kinetics of many pH-dependent enzymes is also studied. This second model includes both substrate inhibition kinetics as well as autocatalysis through the activation of the enzyme by hydrogen ion. We find that it is the autocatalysis that is always responsible for oscillatory behavior in this scheme. The substrate inhibition terms affect the steady-state behavior but do not lead to oscillations unless product inhibition or multiple substrates are present; this is a general conclusion we can draw from our studies of both the classical substrate inhibition scheme and the pH-dependent enzyme mechanism. Finally, an analysis of the nullclines for these two models allows us to prove that the nullcline slopes must have a negative value for oscillatory behavior to exist; this proof can explain our results. From our analysis, we conclude with a brief discussion of other enzymes that might be expected to produce oscillatory behavior based on a pH-dependent substrate inhibition mechanism.     

Numerical modelling of label-structured cell population growth using CFSE distribution data

Background

Drug-drug interactions resulting from the inhibition of an enzymatic process can have serious implications for clinical drug therapy. Quantification of the drugs internal exposure increase upon administration with an inhibitor requires understanding to avoid the drug reaching toxic thresholds. In this study, we aim to predict the effect of the CYP3A4 inhibitors, itraconazole (ITZ) and its primary metabolite, hydroxyitraconazole (OH-ITZ) on the pharmacokinetics of the anesthetic, midazolam (MDZ) and its metabolites, 1' hydroxymidazolam (1OH-MDZ) and 1' hydroxymidazolam glucuronide (1OH-MDZ-Glu) using mechanistic whole body physiologically-based pharmacokinetic simulation models. The model is build on MDZ, 1OH-MDZ and 1OH-MDZ-Glu plasma concentration time data experimentally determined in 19 CYP3A5 genotyped adult male individuals, who received MDZ intravenously in a basal state. The model is then used to predict MDZ, 1OH-MDZ and 1OH-MDZ-Glu concentrations in an CYP3A-inhibited state following ITZ administration.

Results

For the basal state model, three linked WB-PBPK models (MDZ, 1OH-MDZ, 1OH-MDZ-Glu) for each individual were elimination optimized that resulted in MDZ and metabolite plasma concentration time curves that matched individual observed clinical data. In vivo K_m and V_max optimized values for MDZ hydroxylation were similar to literature based in vitro measures. With the addition of the ITZ/OH-ITZ model to each individual coupled MDZ + metabolite model, the plasma concentration time curves were predicted to greatly increase the exposure of MDZ as well as to both increase exposure and significantly alter the plasma concentration time curves of the MDZ metabolites in comparison to the basal state curves. As compared to the observed clinical data, the inhibited state curves were generally well described although the simulated concentrations tended to exceed the experimental data between approximately 6 to 12 hours following MDZ administration. This deviations appeared to be greater in the CYP3A5 *1/*1 and CYP3A5 *1/*3 group than in the CYP3A5 *3/*3 group and was potentially the result of assuming that ITZ/OH-ITZ inhibits both CYP3A4 and CYP3A5, whereas in vitro inhibition is due to CYP3A4.

Conclusion

This study represents the first attempt to dynamically simulate metabolic enzymatic drug-drug interactions via coupled WB-PBPK models. The workflow described herein, basal state optimization followed by inhibition prediction, is novel and will provide a basis for the development of other inhibitor models that can be used to guide, interpret, and potentially replace clinical drug-drug interaction trials.     

Incorporating qualitative knowledge in enzyme kinetic models using fuzzy logic

Modeling of metabolic pathway dynamics requires detailed kinetic equations at the enzyme level. In particular, the kinetic equations must account for metabolite effectors that contribute significantly to the pathway regulation in vivo. Unfortunately, most kinetic rate laws available in the literature do not consider all the effectors simultaneously, and much kinetic information exists in a qualitative or semiquantitative form. In this article, we present a strategy to incorporate such information into the kinetic equation. This strategy uses fuzzy logic‐based factors to modify algebraic rate laws that account for partial kinetic characteristics. The parameters introduced by the fuzzy factors are then optimized by use of a hybrid of simplex and genetic algorithms. The resulting model provides a flexible form that can simulate various kinetic behaviors. Such kinetic models are suitable for pathway modeling without complete enzyme mechanisms. Three enzymes in Escherichia coli central metabolism are used as examples: phosphoenolpyruvate carboxylase; phosphoenolpyruvate carboxykinase; and pyruvate kinase I. Results show that, with fuzzy logic‐augmented models, the kinetic data can be much better described. In particular, complex behavior, such as allosteric inhibition, can be captured using fuzzy rules. The resulting models, even though they do not provide additional physical meaning in enzyme mechanisms, allow the model to incorporate semiquantitative information in metabolic pathway models. © 1999 John Wiley & Sons, Inc. Biotechnol Bioeng 62: 722–729, 1999.