Many promising drug candidates metabolized by aldehyde oxidase (AOX) fail during clinical trial owing to underestimation of their clearance. AOX is species-specific, which makes traditional allometric studies a poor choice for estimating human clearance. Other studies have suggested using half-life calculated by measuring substrate depletion to measure clearance. In this study, we proposed using numerical fitting to enzymatic pathways other than Michaelis-Menten (MM) to avoid missing the initial high turnover rate of product formation. Here, product formation over a 240-minute time course of six AOX substrates-, and zoniporide-have been provided to illustrate enzyme deactivation over time to help better understand why MM kinetics sometimes leads to underestimation of rate constants. Based on the data provided in this article, the total velocity for substrates becomes slower than the initial velocity by 3.1-, 6.5-, 2.9-, 32.2-, 2.7-, and 0.2-fold, respectively, in human expressed purified enzyme, whereas the K m remains constant. Also, our studies on the role of reactive oxygen species (ROS), such as superoxide and hydrogen peroxide, show that ROS did not significantly alter the change in enzyme activity over time. Providing a new electron acceptor, 5-nitroquinoline, did, however, alter the change in rate over time for mumerous compounds. The data also illustrate the difficulties in using substrate disappearance to estimate intrinsic clearance.

Functional cross-talk among human drug-metabolizing cytochromes P450 through their association is a topic of emerging importance. Here we studied the interactions of human CYP2D6, a major metabolizer of psychoactive drugs, with one of the most prevalent human P450 enzymes, ethanol-inducible CYP2E1. Detection of P450-P450 interactions was accomplished through luminescence resonance energy transfer (LRET) between labeled proteins incorporated into human liver microsomes and the microsomes of insect cells containing NADPH-cytochrome P450 reductase. The potential of CYP2D6 to form oligomers in the microsomal membrane is among the highest observed with human cytochromes P450 studied up to date. We also observed the formation of heteromeric complexes of CYP2D6 with CYP2E1 and CYP3A4, and found a significant modulation of these interactions by MDMA, a widespread drug of abuse metabolized by CYP2D6. Our results demonstrate an ample alteration of the catalytic properties of CYP2D6 and CYP2E1 caused by their association. In particular, we demonstrated that pre-incubation of microsomes containing co-incorporated CYP2D6 and CYP2E1 with CYP2D6-specific substrates resulted in considerable time-dependent activation of CYP2D6, which presumably occurs via a slow substrate-induced reorganization of CYP2E1-CYP2D6 heterooligomers. Furthermore, we demonstrated that the formation of heteromeric complexes between CYP2E1 and CYP2D6 affects the stoichiometry of futile cycling and substrate oxidation by CYP2D6 by means of decreasing the electron leakage through the peroxide-generating pathways. Our results further emphasize the role of P450-P450 interactions in regulatory cross-talk in human drug-metabolizing ensemble and suggest a role of interactions of CYP2E1 with CYP2D6 in pharmacologically important instances of alcohol-drug interactions.

We investigate the mechanism of time-dependent inhibition (TDI) of human cytochrome P450 2D6 (CYP2D6) by 3,4-methylenedioxymethamphetamine (MDMA, ecstasy), one of the most widespread recreational drugs of abuse. In an effort to unravel the kinetic mechanism of the formation of metabolic inhibitory complex (MIC) of CYP2D6 with MDMA-derived carbene we carried out a series of spectrophotometric studies paralleled with registration of the kinetics of time-dependent inhibition (TDI) in CYP2D6-incorporated proteoliposomes. The high amplitude of spectral signal in this system allowed us to characterize the spectral properties of the formed MIC in details and obtain an accurate spectral signature of MIC formation. This information was then used in the studies with CYP2D6-containing microsomes of insect cells (CYP2D6 Supersomes™). Our results demonstrate that in both systems the formation of the ferrous carbene-derived MIC is relatively slow, reversible and is not associated with the accumulation of the ferric carbene intermediate, as takes place in the case of CYP3A4 and podophylotoxin. Furthermore, the limited amplitude of MIC formation suggests that only a fraction (∼50%) of spectrally detectable CYP2D6 in both proteoliposomes and Supersomes participates in the formation of MIC and is therefore involved in the MDMA metabolism. This observation reveals yet another example of a cytochrome P450 that exhibits persistent functional heterogeneity of its population in microsomal membranes. Our study provides a solid methodological background for further mechanistic studies of MIC formation in human liver microsomes and demonstrates that the potency and physiological relevance of MDMA-dependent TDI of CYP2D6 may be overestimated.

Methylenedioxymethamphetamine (MDMA) is a known drug of abuse and schedule 1 narcotic under the Controlled Substances Act. Previous pharmacokinetic work on MDMA used classic linearization techniques to conclude irreversible mechanism-based inhibition of CYP2D6. The current work challenges this outcome by assessing the possibility of two alternative reversible kinetic inhibition mechanisms known as the quasi-irreversible (QI) model and equilibrium model (EM). In addition, progress curve experiments were used to investigate the residual metabolism of MDMA by liver microsomes and CYP2D6 baculosomes over incubation periods up to 30 minutes. These experiments revealed activity in a terminal linear phase at the fractional rates with respect to initial turnover of 0.0354 6 0.0089 in human liver microsomes and 0.0114 6 0.0025 in baculosomes. Numerical model fits to percentage of remaining activity (PRA) data were consistent with progress curve modeling results, wherein an irreversible inhibition pathway was found unnecessary for good fit scoring. Both QI and EM kinetic mechanisms fit the PRA data well, although in CYP2D6 baculosomes the inclusion of an irreversible inactivation pathway did not allow for convergence to a reasonable fit. The kinetic complexity accessible to numerical modeling has been used to determine that MDMA is not an irreversible inactivator of CYP2D6, and instead follows what can be generally referred to as slowly reversible inhibition. SIGNIFICANCE STATEMENTThe work herein describes the usage of computational models to delineate between irreversible and slowly reversible time-dependent inhibition. Such models are used in the paper to analyze MDMA and classify it as a reversible time-dependent inhibitor.

In this work, we present the 3D Kinetics approach as a step forward in the field of enzyme kinetics. Normally in enzyme kinetics, it is first assumed that the kinetics will conform to the assumption of Michaelis and Menten and an experiment is conducted at various concentrations around the concentration that gives half-maximal velocity. Often, these experiments could be compromised by having too much substrate, nonlinear reaction over time, time-dependent or substrate inhibition, or several other kinetic models. Herein, we present a general strategy that will decrease the number of experiments required to develop an accurate representation of the kinetics of an enzymatic reaction. We show that with a single experimental protocol, we can fit a number of the most common kinetic models associated with enzyme-catalyzed reactions. Through this experiment, we introduce the effect of time on saturation curves by modeling the reaction velocity over time and across a set of substrate concentrations. Michaelis-Menten (MM) kinetics and other analytical solutions used to solve more complex kinetic models were introduced to the field of enzymology over a hundred years ago and have only marginally changed over the years with each analytical model requiring a different set of experiments and concentration ranges. Although this approach was necessary at the time, the computational power today makes any such simplifying and limiting efforts unnecessary and avoidable. In this study, we use a single experimental protocol and fit a number of different models to the resulting data. We present four different case studies to compare and contrast the outcomes of 3D Kinetics with MM analysis for different kinetic scenarios such as enzymatic reactions with linear kinetics, biphasic kinetics, substrate inhibition, and time-dependent inhibition (TDI) to confirm the advantage of the 3D Kinetics method to the long-established MM analysis.