Surface plasmon resonance analysis of antifungal azoles binding to CYP3A4 with kinetic resolution of multiple binding orientations

Co-based carbon material as CYP3A4-like nanozyme with both biocatalytic activity and inhibition behaviors

Up to now, the biocatalytic activity of nanozymes has been extensively studied, while little research focus on their inhibitory behaviors. Here, Co-based carbon material (Co-DMOF) containing abundant carboxylic acid groups was prepared, with defects introduced by COx escape during pyrolysis to achieve controllable activity. As a result, Co-DMOF exhibited biocatalytic activity similar to cytochrome P450 3A4 (CYP3A4) in the metabolism of 1,4-Dihydropyridine (1,4-DHP, a calcium channel blocker). Excitingly, studies on IC50 and drug-drug interaction (DDI) suggested that Co-DMOF had similar inhibitory behaviors to CYP3A4. Moreover, Co-DMOF displayed excellent stability even under high temperature (100 °C), organic solvents, and a wide range of pH (4-9). Additionally, it can be reused for at least 7 times with only slight loss of activity. Therefore, Co-DMOF has great potential to become a low-cost alternative to CYP3A4 for drug dosage guideline, drug metabolism and DDI. This work provides more possibilities for expanding the CYP3A4-like nanozyme library.

Multi-temperature experiments to ease analysis of heterogeneous binder solutions by surface plasmon resonance biosensing

Surface Plasmon Resonance (SPR) biosensing is a well-established tool for the investigation of binding kinetics between a soluble species and an immobilized (bio)molecule. While robust and accurate data analysis techniques are readily available for single species, methods to exploit data collected with a solution containing multiple interactants are scarce. In a previous study, our group proposed two data analysis algorithms for (1) the precise and reliable identification of the kinetic parameters of N interactants present at different ratios in N mixtures and (2) the estimation of the composition of a given mixture, assuming that the kinetic parameters and the total concentration of all interactants are known. Here, we extend the first algorithm by reducing the number of necessary mixtures. This is achieved by conducting experiments at different temperatures. Through the Van't Hoff and Eyring equations, identifying the kinetic and thermodynamic parameters of N binders becomes possible with M mixtures with M comprised between 2 and N and at least N/M temperatures. The second algorithm is improved by adding the total analyte concentration as a supplementary variable to be identified in an optimization routine. We validated our analysis framework experimentally with a system consisting of mixtures of low molecular weight drugs, each competing to bind to an immobilized protein. We believe that the analysis of mixtures and composition estimation could pave the way for SPR biosensing to become a bioprocess monitoring tool, on top of expanding its already substantial role in drug discovery and development.

Numerical Methods for Modeling Enzyme Kinetics

Differential equations are used to describe time-dependent changes in enzyme kinetics and pharmacokinetics. Analytical and numerical methods can be used to solve differential equations. This chapter describes the use of numerical methods in solving differential equations and its applications in characterizing the complexities observed in enzyme kinetics. A discussion is included on the use of numerical methods to overcome limitations of explicit equations in the analysis of metabolism kinetics, reversible inhibition kinetics, and inactivation kinetics. The chapter describes the advantages of using numerical methods when Michaelis-Menten assumptions do not hold.

On the Use of Surface Plasmon Resonance Biosensing to Understand IgG-FcγR Interactions

Surface plasmon resonance (SPR)-based optical biosensors offer real-time and label-free analysis of protein interactions, which has extensively contributed to the discovery and development of therapeutic monoclonal antibodies (mAbs). As the biopharmaceutical market for these biologics and their biosimilars is rapidly growing, the role of SPR biosensors in drug discovery and quality assessment is becoming increasingly prominent. One of the critical quality attributes of mAbs is the N-glycosylation of their Fc region. Other than providing stability to the antibody, the Fc N-glycosylation influences immunoglobulin G (IgG) interactions with the Fcγ receptors (FcγRs), modulating the immune response. Over the past two decades, several studies have relied on SPR-based assays to characterize the influence of N-glycosylation upon the IgG-FcγR interactions. While these studies have unveiled key information, many conclusions are still debated in the literature. These discrepancies can be, in part, attributed to the design of the reported SPR-based assays as well as the methodology applied to SPR data analysis. In fact, the SPR biosensor best practices have evolved over the years, and several biases have been pointed out in the development of experimental SPR protocols. In parallel, newly developed algorithms and data analysis methods now allow taking into consideration complex biomolecular kinetics. In this review, we detail the use of different SPR biosensing approaches for characterizing the IgG-FcγR interactions, highlighting their merit and inherent experimental complexity. Furthermore, we review the latest SPR-derived conclusions on the influence of the N-glycosylation upon the IgG-FcγR interactions and underline the differences and similarities across the literature. Finally, we explore new avenues taking advantage of novel computational analysis of SPR results as well as the latest strategies to control the glycoprofile of mAbs during production, which could lead to a better understanding and modelling of the IgG-FcγRs interactions.

Determination of the composition of heterogeneous binder solutions by surface plasmon resonance biosensing

Surface plasmon resonance-based biosensors have been extensively applied to the characterization of the binding kinetics between purified (bio)molecules, thanks to robust data analysis techniques. However, data analysis for solutions containing multiple interactants is still at its infancy. We here present two algorithms for (1) the reliable and accurate determination of the kinetic parameters of N interactants present at different ratios in N mixtures and (2) the estimation of the ratios of each interactant in a given mixture, assuming that their kinetic parameters are known. Both algorithms assume that the interactants compete to bind to an immobilized ligand in a 1:1 fashion and necessitate prior knowledge of the total concentration of all interactants combined. The effectiveness of these two algorithms was experimentally validated with a model system corresponding to mixtures of four small molecular weight drugs binding to an immobilized protein. This approach enables the in-depth characterization of mixtures using SPR, which may be of considerable interest for many drug discovery or development applications, notably for protein glycovariant analysis.

Comparison of the inhibitory effects of azole antifungals on cytochrome P450 3A4 genetic variants

Cytochrome P450 (CYP) 3A4 is one of the major drug-metabolizing enzymes. Genetic variants of CYP3A4 with altered activity are one of the factors responsible for interindividual differences in drug metabolism. Azole antifungals inhibit CYP3A4 to cause clinically significant drug-drug interactions. In the present quantitative study, we investigated the inhibitory effects of three azole antifungals (ketoconazole, voriconazole, and fluconazole) on testosterone metabolism by recombinant CYP3A4 genetic variants (CYP3A4.1 (WT), CYP3A4.2, CYP3A4.7, CYP3A4.16, and CYP3A4.18) and compared them with those previously reported for itraconazole. The inhibition constants (K_i) of ketoconazole, voriconazole, and fluconazole for rCYP3A4.1 were 3.6 nM, 3.2 μM, and 16.1 μM, respectively. The K_i values of these azoles for rCYP3A4.16 were 13.9-, 13.6-, and 6.2-fold higher than those for rCYP3A4.1, respectively, whereas the K_i value of itraconazole for rCYP3A4.16 was 0.54-fold of that for rCYP3A4.1. The other genetic variants had similar effects on the K_i values of the three azoles, whereas a very different pattern was seen for itraconazole. In conclusion, itraconazole has unique characteristics that are distinct from those shared by the other azole anti-fungal drugs ketoconazole, voriconazole, and fluconazole with regard to the influence of genetic variations on the inhibition of CYP3A4.

Kinetics of cytochrome P450 3A4 inhibition by heterocyclic drugs defines a general sequential multistep binding process

Cytochrome P450 (P450) 3A4 is the enzyme most involved in the metabolism of drugs and can also oxidize numerous steroids. This enzyme is also involved in one-half of pharmacokinetic drug-drug interactions, but details of the exact mechanisms of P450 3A4 inhibition are still unclear in many cases. Ketoconazole, clotrimazole, ritonavir, indinavir, and itraconazole are strong inhibitors; analysis of the kinetics of reversal of inhibition with the model substrate 7-benzoyl quinoline showed lag phases in several cases, consistent with multiple structures of P450 3A4 inhibitor complexes. Lags in the onset of inhibition were observed when inhibitors were added to P450 3A4 in 7-benzoyl quinoline O-debenzylation reactions, and similar patterns were observed for inhibition of testosterone 6β-hydroxylation by ritonavir and indinavir. Upon mixing with inhibitors, P450 3A4 showed rapid binding as judged by a spectral shift with at least partial high-spin iron character, followed by a slower conversion to a low-spin iron-nitrogen complex. The changes were best described by two intermediate complexes, one being a partial high-spin form and the second another intermediate, with half-lives of seconds. The kinetics could be modeled in a system involving initial loose binding of inhibitor, followed by a slow step leading to a tighter complex on a multisecond time scale. Although some more complex possibilities cannot be dismissed, these results describe a system in which conformationally distinct forms of P450 3A4 bind inhibitors rapidly and two distinct P450-inhibitor complexes exist en route to the final enzyme-inhibitor complex with full inhibitory activity.

Structure-guided approaches to targeting stress responses in human fungal pathogens

Fungi inhabit extraordinarily diverse ecological niches, including the human body. Invasive fungal infections have a devastating impact on human health worldwide, killing ∼1.5 million individuals annually. The majority of these deaths are attributable to species of Candida, Cryptococcus, and Aspergillus Treating fungal infections is challenging, in part due to the emergence of resistance to our limited arsenal of antifungal agents, necessitating the development of novel therapeutic options. Whereas conventional antifungal strategies target proteins or cellular components essential for fungal growth, an attractive alternative strategy involves targeting proteins that regulate fungal virulence or antifungal drug resistance, such as regulators of fungal stress responses. Stress response networks enable fungi to adapt, grow, and cause disease in humans and include regulators that are highly conserved across eukaryotes as well as those that are fungal-specific. This review highlights recent developments in elucidating crystal structures of fungal stress response regulators and emphasizes how this knowledge can guide the design of fungal-selective inhibitors. We focus on the progress that has been made with highly conserved regulators, including the molecular chaperone Hsp90, the protein phosphatase calcineurin, and the small GTPase Ras1, as well as with divergent stress response regulators, including the cell wall kinase Yck2 and trehalose synthases. Exploring structures of these important fungal stress regulators will accelerate the design of selective antifungals that can be deployed to combat life-threatening fungal diseases.

A Versatile Chemoenzymatic Synthesis for the Discovery of Potent Cryptophycin Analogs

The cryptophycins are a family of macrocyclic depsipeptide natural products that display exceptionally potent antiproliferative activity against drug-resistant cancers. Unique challenges facing the synthesis and derivatization of this complex group of molecules motivated us to investigate a chemoenzymatic synthesis designed to access new analogs for biological evaluation. The cryptophycin thioesterase (CrpTE) and the cryptophycin epoxidase (CrpE) are a versatile set of enzymes that catalyze macrocyclization and epoxidation of over 20 natural cryptophycin metabolites. Thus, we envisioned a drug development strategy involving their use as standalone biocatalysts for production of unnatural derivatives. Herein, we developed a scalable synthesis of 12 new unit A-B-C-D linear chain elongation intermediates containing heterocyclic aromatic groups as alternatives to the native unit A benzyl group. N-Acetyl cysteamine activated forms of each intermediate were assessed for conversion to macrocyclic products using wild type CrpTE, which demonstrated the exceptional flexibility of this enzyme. Semipreparative scale reactions were conducted for isolation and structural characterization of new cryptophycins. Each was then evaluated as a substrate for CrpE P450 and its ability to generate the epoxidized products from these substrates that possess altered electronics at the unit A styrenyl double bond position. Finally, biological evaluation of the new cryptophycins revealed a des-β-epoxy analog with low picomolar potency, previously limited to cryptophycins bearing epoxide functionality.

Modelling of substrate access and substrate binding to cephalosporin acylases

Semisynthetic cephalosporins are widely used antibiotics currently produced by different chemical steps under harsh conditions, which results in a considerable amount of toxic waste. Biocatalytic synthesis by the cephalosporin acylase from Pseudomonas sp. strain N176 is a promising alternative. Despite intensive engineering of the enzyme, the catalytic activity is still too low for a commercially viable process. To identify the bottlenecks which limit the success of protein engineering efforts, a series of MD simulations was performed to study for two acylase variants (WT, M6) the access of the substrate cephalosporin C from the bulk to the active site and the stability of the enzyme-substrate complex. In both variants, cephalosporin C was binding to a non-productive substrate binding site (E86α, S369β, S460β) at the entrance to the binding pocket, preventing substrate access. A second non-productive binding site (G372β, W376β, L457β) was identified within the binding pocket, which competes with the active site for substrate binding. Noteworthy, substrate binding to the protein surface followed a Langmuir model resulting in binding constants K = 7.4 and 9.2 mM for WT and M6, respectively, which were similar to the experimentally determined Michaelis constants KM = 11.0 and 8.1 mM, respectively.

Impact of Lipid Partitioning on the Design, Analysis, and Interpretation of Microsomal Time-Dependent Inactivation

Nonspecific drug partitioning into microsomal membranes must be considered for in vitro-in vivo correlations. This work evaluated the effect of including lipid partitioning in the analysis of complex TDI kinetics with numerical methods. The covariance between lipid partitioning and multiple inhibitor binding was evaluated. Simulations were performed to test the impact of lipid partitioning on the interpretation of TDI kinetics, and experimental TDI datasets for paroxetine (PAR) and itraconazole (ITZ) were modeled. For most kinetic schemes, modeling lipid partitioning results in statistically better fits. For MM-IL simulations (KI,u = 0.1 µM, kinact = 0.1 minute-1), concurrent modeling of lipid partitioning for an fumic range (0.01, 0.1, and 0.5) resulted in better fits compared with post hoc correction (AICc: -526 vs. -496, -579 vs. -499, and -636 vs. -579, respectively). Similar results were obtained with EII-IL. Lipid partitioning may be misinterpreted as double binding, leading to incorrect parameter estimates. For the MM-IL datasets, when fumic = 0.02, MM-IL, and EII model fits were indistinguishable (δAICc = 3). For less partitioned datasets (fumic = 0.1 or 0.5), the inclusion of partitioning resulted in better models. The inclusion of lipid partitioning can lead to markedly different estimates of KI,u and kinact A reasonable alternate experimental design is nondilution TDI assays, with post hoc fumic incorporation. The best fit models for PAR (MIC-M-IL) and ITZ (MIC-EII-M-IL and MIC-EII-M-Seq-IL) were consistent with their reported mechanism and kinetics. Overall, experimental fumic values should be concurrently incorporated into TDI models with complex kinetics, when dilution protocols are used.

Human cytochrome P450 enzymes bind drugs and other substrates mainly through conformational-selection modes

Cytochrome P450 (P450) enzymes are major catalysts involved in the oxidations of most drugs, steroids, carcinogens, fat-soluble vitamins, and natural products. The binding of substrates to some of the 57 human P450s and other mammalian P450s is more complex than a two-state system and has been proposed to involve mechanisms such as multiple ligand occupancy, induced-fit, and conformational-selection. Here, we used kinetic analysis of binding with multiple concentrations of substrates and computational modeling of these data to discern possible binding modes of several human P450s. We observed that P450 2D6 binds its ligand rolapitant in a mechanism involving conformational-selection. P450 4A11 bound the substrate lauric acid via conformational-selection, as did P450 2C8 with palmitic acid. Binding of the steroid progesterone to P450 21A2 was also best described by a conformational-selection model. Hexyl isonicotinate binding to P450 2E1 could be described by either a conformational-selection or an induced-fit model. Simulation of the binding of the ligands midazolam, bromocriptine, testosterone, and ketoconazole to P450 3A4 was consistent with an induced-fit or a conformational-selection model, but the concentration dependence of binding rates for varying both P450 3A4 and midazolam concentrations revealed discordance in the parameters, indicative of conformational-selection. Binding of the P450s 2C8, 2D6, 3A4, 4A11, and 21A2 was best described by conformational-selection, and P450 2E1 appeared to fit either mode. These findings highlight the complexity of human P450-substrate interactions and that conformational-selection is a dominant feature of many of these interactions.

Novel Tetrazole-Containing Analogues of Itraconazole as Potent Antiangiogenic Agents with Reduced Cytochrome P450 3A4 Inhibition

Itraconazole has been found to possess potent antiangiogenic activity, exhibiting promising antitumor activity in several human clinical studies. The wider use of itraconazole in the treatment of cancer, however, has been limited by its potent inhibition of the drug metabolizing enzyme cytochrome P450 3A4 (CYP3A4). In an effort to eliminate the CYP3A4 inhibition while retaining its antiangiogenic activity, we designed and synthesized a series of derivatives in which the 1,2,4-triazole ring is replaced with various azoles and nonazoles. Among these analogues, 15n with tetrazole in place of 1,2,4-triazole exhibited optimal inhibition of human umbilical vein endothelial cell proliferation with an IC₅₀ of 73 nM without a significant effect on CYP3A4 (EC₅₀ > 20 μM). Similar to itraconazole, 15n induced Niemann-Pick C phenotype (NPC phenotype) and blocked AMPK/mechanistic target of rapamycin signaling. These results suggest that 15n is a promising angiogenesis inhibitor that can be used in combination with most other known anticancer drugs.

[Methods for determining of cytochrome P450 isozymes functional activity]

The review is dedicated to modern methods and technologies for determining of cytochrome P450 isozymes functional activity, such as absorbance and fluorescent spectroscopy, electron paramagnetic resonance (EPR), nuclear magnetic resonance (NMR), Raman, Mossbauer, and X-ray spectroscopy, surface plasmon resonance (SPR), atomic force microscopy (AFM). Methods of molecular genetic analysis were reviewed from personalized medicine point of view. The use of chromate-mass-spectrometric methods for cytochrome P450-dependent catalytic reactions' products was discussed. The review covers modern electrochemical systems based on cytochrome P450 isozymes for their catalytic activity analysis, their use in practice and further development perspectives for experimental pharmacology, biotechnology and translational medicine.

Advances in the Understanding of Protein-Protein Interactions in Drug Metabolizing Enzymes through the Use of Biophysical Techniques

In recent years, a growing appreciation has developed for the importance of protein-protein interactions to modulate the function of drug metabolizing enzymes. Accompanied with this appreciation, new methods and technologies have been designed for analyzing protein-protein interactions both in vitro and in vivo. These technologies have been applied to several classes of drug metabolizing enzymes, including: cytochrome P450's (CYPs), monoamine oxidases (MAOs), UDP-glucuronosyltransferases (UGTs), glutathione S-transferases (GSTs), and sulfotransferases (SULTs). In this review, we offer a brief description and assessment of the impact of many of these technologies to the study of protein-protein interactions in drug disposition. The still expanding list of these techniques and assays has the potential to revolutionize our understanding of how these enzymes carry out their important functions in vivo.

Membrane Fluidity Modulates Thermal Stability and Ligand Binding of Cytochrome P4503A4 in Lipid Nanodiscs

Cytochrome P4503A4 (CYP3A4) is a peripheral membrane protein that plays a major role in enzymatic detoxification of many drugs and toxins. CYP3A4 has an integral membrane N-terminal helix and a localized patch comprised of the G' and F' helix regions that are embedded in the membrane, but the effects of membrane composition on CYP3A4 function are unknown. Here, circular dichroism and differential scanning calorimetry were used to compare the stability of CYP3A4 in lipid bilayer nanodiscs with varying ratios of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine to 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). These lipids differ in the acyl-chain length and their degree of unsaturation. The thermal denaturation of CYP3A4 in nanodiscs occurs in a temperature range distinct from that of the nanodisc denaturation so it can be monitored calorimetrically. Melting temperatures (T_m), heat capacities (ΔC_p), and calorimetric enthalpies (ΔH_cal) for denaturation of CYP3A4 each increased with an increasing fraction of DMPC, with a maximum at 50% DMPC, before decreasing at 75% DMPC. Addition of the inhibitor ketoconazole results in increased thermal stability, and larger ΔC_p and ΔH_cal values, with different sensitivities to lipid composition. Effects of lipid composition on ligand binding dynamics were also studied. Equilibrium binding affinities of both ketoconazole (KTZ) and testosterone (TST) were minimally affected by lipid composition. However, stopped-flow analyses indicate that the rates of KTZ binding reach a maximum in membranes containing 50% DMPC, whereas the rate of TST binding decreases continuously with an increasing DMPC concentration. These results indicate that CYP3A4 is highly sensitive to the acyl-chain composition of the lipids and fluidity of the membrane in which it is embedded.

Optimizing Multiple Analyte Injections in Surface Plasmon Resonance Biosensors with Analytes having Different Refractive Index Increments

Surface plasmon resonance-based biosensors have been successfully applied to the study of the interactions between macromolecules and small molecular weight compounds. In an effort to increase the throughput of these SPR-based experiments, we have already proposed to inject multiple compounds simultaneously over the same surface. When specifically applied to small molecular weight compounds, such a strategy would however require prior knowledge of the refractive index increment of each compound in order to correctly interpret the recorded signal. An additional experiment is typically required to obtain this information. In this manuscript, we show that through the introduction of an additional global parameter corresponding to the ratio of the saturating signals associated with each molecule, the kinetic parameters could be identified with similar confidence intervals without any other experimentation.

Mucin-like Region of Herpes Simplex Virus Type 1 Attachment Protein Glycoprotein C (gC) Modulates the Virus-Glycosaminoglycan Interaction

Glycoprotein C (gC) mediates the attachment of HSV-1 to susceptible host cells by interacting with glycosaminoglycans (GAGs) on the cell surface. gC contains a mucin-like region located near the GAG-binding site, which may affect the binding activity. Here, we address this issue by studying a HSV-1 mutant lacking the mucin-like domain in gC and the corresponding purified mutant protein (gCΔmuc) in cell culture and GAG-binding assays, respectively. The mutant virus exhibited two functional alterations as compared with native HSV-1 (i.e. decreased sensitivity to GAG-based inhibitors of virus attachment to cells and reduced release of viral particles from the surface of infected cells). Kinetic and equilibrium binding characteristics of purified gC were assessed using surface plasmon resonance-based sensing together with a surface platform consisting of end-on immobilized GAGs. Both native gC and gCΔmuc bound via the expected binding region to chondroitin sulfate and sulfated hyaluronan but not to the non-sulfated hyaluronan, confirming binding specificity. In contrast to native gC, gCΔmuc exhibited a decreased affinity for GAGs and a slower dissociation, indicating that once formed, the gCΔmuc-GAG complex is more stable. It was also found that a larger number of gCΔmuc bound to a single GAG chain, compared with native gC. Taken together, our data suggest that the mucin-like region of HSV-1 gC is involved in the modulation of the GAG-binding activity, a feature of importance both for unrestricted virus entry into the cells and release of newly produced viral particles from infected cells.

Current Approaches for Investigating and Predicting Cytochrome P450 3A4-Ligand Interactions

Cytochrome P450 3A4 (CYP3A4) is the major and most important drug-metabolizing enzyme in humans that oxidizes and clears over a half of all administered pharmaceuticals. This is possible because CYP3A4 is promiscuous with respect to substrate binding and has the ability to catalyze diverse oxidative chemistries in addition to traditional hydroxylation reactions. Furthermore, CYP3A4 binds and oxidizes a number of substrates in a cooperative manner and can be both induced and inactivated by drugs. In vivo, CYP3A4 inhibition could lead to undesired drug-drug interactions and drug toxicity, a major reason for late-stage clinical failures and withdrawal of marketed pharmaceuticals. Owing to its central role in drug metabolism, many aspects of CYP3A4 catalysis have been extensively studied by various techniques. Here, we give an overview of experimental and theoretical methods currently used for investigation and prediction of CYP3A4-ligand interactions, a defining factor in drug metabolism, with an emphasis on the problems addressed and conclusions derived from the studies.

Use of chemical auxiliaries to control p450 enzymes for predictable oxidations at unactivated C-h bonds of substrates

Cytochrome P450 enzymes (P450s) have the ability to oxidize unactivated C-H bonds of substrates with remarkable regio- and stereoselectivity. Comparable selectivity for chemical oxidizing agents is typically difficult to achieve. Hence, there is an interest in exploiting P450s as potential biocatalysts. Despite their impressive attributes, the current use of P450s as biocatalysts is limited. While bacterial P450 enzymes typically show higher activity, they tend to be highly selective for one or a few substrates. On the other hand, mammalian P450s, especially the drug-metabolizing enzymes, display astonishing substrate promiscuity. However, product prediction continues to be challenging. This review discusses the use of small molecules for controlling P450 substrate specificity and product selectivity. The focus will be on two approaches in the area: (1) the use of decoy molecules, and (2) the application of substrate engineering to control oxidation by the enzyme.

Probing the origins of aromatase inhibitory activity of disubstituted coumarins via QSAR and molecular docking

This study investigated the quantitative structure-activity relationship (QSAR) of imidazole derivatives of 4,7-disubstituted coumarins as inhibitors of aromatase, a potential therapeutic protein target for the treatment of breast cancer. Herein, a series of 3,7- and 4,7-disubstituted coumarin derivatives (1-34) with R1 and R2 substituents bearing aromatase inhibitory activity were modeled as a function of molecular and quantum chemical descriptors derived from low-energy conformer geometrically optimized at B3LYP/6-31G(d) level of theory. Insights on origins of aromatase inhibitory activity was afforded by the computed set of 7 descriptors comprising of F10[N-O], Inflammat-50, Psychotic-80, H-047, BELe1, B10[C-O] and MAXDP. Such significant descriptors were used for QSAR model construction and results indicated that model 4 afforded the best statistical performance. Good predictive performance were achieved as verified from the internal (comprising the training and the leave-one-out cross-validation (LOO-CV) sets) and external sets affording the following statistical parameters: R (2) Tr = 0.9576 and RMSETr = 0.0958 for the training set; Q (2) CV = 0.9239 and RMSECV = 0.1304 for the LOO-CV set as well as Q (2) Ext = 0.7268 and RMSEExt = 0.2927 for the external set. Significant descriptors showed correlation with functional substituents, particularly, R1 in governing high potency as aromatase inhibitor. Molecular docking calculations suggest that key residues interacting with the coumarins were predominantly lipophilic or non-polar while a few were polar and positively-charged. Findings illuminated herein serve as the impetus that can be used to rationally guide the design of new aromatase inhibitors.

Controlling substrate specificity and product regio- and stereo-selectivities of P450 enzymes without mutagenesis

P450 enzymes (P450s) are well known for their ability to oxidize unactivated CH bonds with high regio- and stereoselectivity. Hence, there is emerging interest in exploiting P450s as potential biocatalysts. Although bacterial P450s typically show higher activity than their mammalian counterparts, they tend to be more substrate selective. Most drug-metabolizing P450s on the other hand, display remarkable substrate promiscuity, yet product prediction remains challenging. Protein engineering is one established strategy to overcome these issues. A less explored, yet promising alternative involves substrate engineering. This review discusses the use of small molecules for controlling the substrate specificity and product selectivity of P450s. The focus is on two approaches, one taking advantage of non-covalent decoy molecules, and the other involving covalent substrate modifications.

A numerical method for analysis of in vitro time-dependent inhibition data. Part 1. Theoretical considerations

Inhibition of cytochromes P450 by time-dependent inhibitors (TDI) is a major cause of clinical drug-drug interactions. It is often difficult to predict in vivo drug interactions based on in vitro TDI data. In part 1 of these manuscripts, we describe a numerical method that can directly estimate TDI parameters for a number of kinetic schemes. Datasets were simulated for Michaelis-Menten (MM) and several atypical kinetic schemes. Ordinary differential equations were solved directly to parameterize kinetic constants. For MM kinetics, much better estimates of KI can be obtained with the numerical method, and even IC50 shift data can provide meaningful estimates of TDI kinetic parameters. The standard replot method can be modified to fit non-MM data, but normal experimental error precludes this approach. Non-MM kinetic schemes can be easily incorporated into the numerical method, and the numerical method consistently predicts the correct model at errors of 10% or less. Quasi-irreversible inactivation and partial inactivation can be modeled easily with the numerical method. The utility of the numerical method for the analyses of experimental TDI data is provided in our companion manuscript in this issue of Drug Metabolism and Disposition (Korzekwa et al., 2014b).

SPR-biosensor assay for analysis of small compounds interaction with human cytochrome P450 51A1 (CYP51A1)

The SPR assay for human cytochrome P450 51A1's (CYP51A1) ligand screening was developed. Assay has been validated with known azole inhibitors of cytochrome P450s. The studied azoles selectively interacted with human cytochrome P450 51A1, which showed the highest affinity towards ketoconazole. The efficiency of the SPR assay was showed with 19 steroid and triterpene compounds, which were not investigated as potential ligands of CYP51A1.

Understanding the mechanism of cytochrome P450 3A4: recent advances and remaining problems

Cytochromes P450 (CYPs) represent a diverse group of heme-thiolate proteins found in almost all organisms. CYPs share a common protein fold but differ in substrate selectivity and catalyze a wide variety of monooxygenation reactions via activation of molecular oxygen. Among 57 human P450s, the 3A4 isoform (CYP3A4) is the most abundant and the most important because it metabolizes the majority of administered drugs. A remarkable feature of CYP3A4 is its extreme promiscuity in substrate specificity and cooperative substrate binding, which often leads to undesirable drug-drug interactions and toxic side effects. Owing to its importance in drug development and therapy, CYP3A4 has been the most extensively studied mammalian P450. In this review we provide an overview on recent progress and remaining problems in the CYP3A4 research.

Binding free energies of inhibitors to iron porphyrin complex as a model for Cytochrome P450

The binding free energies of the inhibitor-heme model complexes are calculated using the density functional methods and the implicit solvation models in water, where the 16 structurally diverse compounds with a spectrum of IC(50) values from 0.05 (clotrimazole) to 1000 (piroxicam) μM are chosen as inhibitors for Cytochrome P450 3A4 (CYP3A4). CYP3A4 is the most predominant constituent of the human hepatic CYP enzymes that play a role in metabolizing structurally diverse xenobiotics. The observed free energy change for each inhibitory binding, ΔG inh0, is obtained from its IC(50) value. The total binding free energy (ΔG b0) of each inhibitor-heme model complex is calculated by the sum of its relative free energy (ΔG(0) ) in the gas phase and solvation free energy to the water-heme model complex. The UB3LYP/LanL2DZ level of theory provides the correct relative stabilities of the high- and low-spin states for the penta- and hexa-coordinated ferric complexes, respectively. The optimized distances of the inhibitor nitrogen (or water oxygen) and the methyl mercaptide S to the ferric iron of the inhibitor-heme model complexes at the same level of theory are consistent with the values of the corresponding X-ray structures, except for the econazole complex. The correlation coefficient r(2) values of 0.91 and 0.75 are obtained from the ΔG b0-ΔG inh0 and ΔG(0) -ΔG inh0 plots, respectively, at the UM06/LanL2DZ:CPCM_UB3LYP/LanL2DZ//UB3LYP/LanL2DZ level of theory in water. This indicates that the total binding free energies calculated for the inhibitor-heme model complexes can be a good descriptor in interpreting the inhibitor binding to CYP3A4 and the relative free energies in the gas phase are mainly responsible for the total binding free energies in water, although the desolvation can be a factor to affect the binding affinity of the inhibitors to CYP3A4. From the theozyme analysis of the X-ray structures for ketoconazole- and metyrapone-CYP3A4 complexes, the interaction free energy of the neighboring residues with each inhibitor in the active site is calculated to be about -3 kcal mol(-1) in water, whose the interaction energy and the desolvation free energy change are about -5 and 2 kcal mol(-1) , respectively.

Stereospecific metabolism of itraconazole by CYP3A4: dioxolane ring scission of azole antifungals

Itraconazole (ITZ) is a mixture of four cis-stereoisomers that inhibit CYP3A4 potently and coordinate CYP3A4 heme via the triazole nitrogen. However, (2R,4S,2'R)-ITZ and (2R,4S,2'S)-ITZ also undergo stereoselective sequential metabolism by CYP3A4 at a site distant from the triazole ring to 3'-OH-ITZ, keto-ITZ, and N-desalkyl-ITZ. This stereoselective metabolism demonstrates specific interactions of ITZ within the CYP3A4 active site. To further investigate this process, the binding and metabolism of the four trans-ITZ stereoisomers by CYP3A4 were characterized. All four trans-ITZ stereoisomers were tight binding inhibitors of CYP3A4-mediated midazolam hydroxylation (IC(50) 16-26 nM), and each gave a type II spectrum upon binding to CYP3A4. However, instead of formation of 3'-OH-ITZ, they were oxidized at the dioxolane ring, leading to ring scission and formation of two new metabolites of ITZ. These two metabolites were also formed from the four cis-ITZ stereoisomers, although not as efficiently. The catalytic rates of dioxolane ring scission were similar to the dissociation rates of ITZ stereoisomers from CYP3A4, suggesting that the heme iron is reduced while the triazole moiety coordinates to it and no dissociation of ITZ is necessary before catalysis. The triazole containing metabolite [1-(2,4-dichlorophenyl)-2-(1H-1,2,4-triazol-1-yl)ethanone] also inhibited CYP3A4 (IC(50) >15 μM) and showed type II binding with CYP3A4. The dioxolane ring scission appears to be clinically relevant because this metabolite was detected in urine samples from subjects that had been administered the mixture of cis-ITZ isomers. These data suggest that the dioxolane ring scission is a metabolic pathway for drugs that contain this moiety.

An enlarged, adaptable active site in CYP164 family P450 enzymes, the sole P450 in Mycobacterium leprae

CYP164 family P450 enzymes are found in only a subset of mycobacteria and include CYP164A1, which is the sole P450 found in Mycobacterium leprae, the causative agent of leprosy. This has previously led to interest in this enzyme as a potential drug target. Here we describe the first crystal structure of a CYP164 enzyme, CYP164A2 from Mycobacterium smegmatis. CYP164A2 has a distinctive, enlarged hydrophobic active site that extends above the porphyrin ring toward the access channels. Unusually, we find that CYP164A2 can simultaneously bind two econazole molecules in different regions of the enlarged active site and is accompanied by the rearrangement and ordering of the BC loop. The primary location is through a classic interaction of the azole group with the porphyrin iron. The second econazole molecule is bound to a unique site and is linked to a tetracoordinated metal ion complexed to one of the heme carboxylates and to the side chains of His 105 and His 364. All of these features are preserved in the closely homologous M. leprae CYP164A1. The computational docking of azole compounds to a homology model of CYP164A1 suggests that these compounds will form effective inhibitors and is supported by the correlation of parallel docking with experimental binding studies of CYP164A2. The binding of econazole to CYP164A2 occurs primarily through the high-spin "open" conformation of the enzyme (K(d) [dissociation constant] of 0.1 μM), with binding to the low-spin "closed" form being significantly hindered (K(d) of 338 μM). These studies support previous suggestions that azole derivatives may provide an effective strategy to improve the treatment of leprosy.

The effects of nitrogen-heme-iron coordination on substrate affinities for cytochrome P450 2E1

A descriptor based computational model was developed for cytochrome P450 2E1 (CYP2E1) based on inhibition constants determined for inhibition of chlorzoxazone, or 4-nitrophenol, metabolism. An empirical descriptor for type II binding was developed and tested for a series of CYP2E1 inhibitors. Inhibition constants where measured for 51 different compounds. A fast 2-dimensional predictive model was developed based on 40 compounds, and tested on 8 compounds of diverse structure. The trained model (n=40) had an r(2) value of 0.76 and an RMSE of 0.48. The correlation between the predicted and actual pK(i) values of the test set of compounds not included in the model gives an r(2) value of 0.78. The features that described binding include heme coordination (type II binding), molecular volume, octanol/water partition coefficient, solvent accessible surface area, and the sum of the atomic polarizabilities. The heme coordination parameter assigns an integer between 0 and 6 depending on structure, and is a new descriptor, based on simple quantum chemical calculations with correction for steric effects. The type II binding parameter was found to be important in obtaining a good correlation between predicted and experimental inhibition constants increasing the r(2) value from 0.38 to 0.77.

The kinetic mechanism for cytochrome P450 metabolism of type II binding compounds: evidence supporting direct reduction

The metabolic stability of a drug is an important property that should be optimized during drug design and development. Nitrogen incorporation is hypothesized to increase the stability by coordination of nitrogen to the heme iron of cytochrome P450, a binding mode that is referred to as type II binding. However, we noticed that the type II binding compound 1 has less metabolic stability at sub-saturating conditions than a closely related type I binding compound 3. Three kinetic models will be presented for type II binder metabolism; (1) Dead-end type II binding, (2) a rapid equilibrium between type I and II binding modes before reduction, and (3) a direct reduction of the type II coordinated heme. Data will be presented on reduction rates of iron, the off rates of substrate (using surface plasmon resonance) and the catalytic rate constants. These data argue against the dead-end, and rapid equilibrium models, leaving the direct reduction kinetic mechanism for metabolism of the type II binding compound 1.

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.

Anilinic N-oxides support cytochrome P450-mediated N-dealkylation through hydrogen-atom transfer

The mechanism of N-dealkylation mediated by cytochrome P450 (P450) has long been studied and argued as either a single electron transfer (SET) or a hydrogen atom transfer (HAT) from the amine to the oxidant of the P450, the reputed iron-oxene. In our study, tertiary anilinic N-oxides were used as oxygen surrogates to directly generate a P450-mediated oxidant that is capable of N-dealkylating the dimethylaniline derived from oxygen donation. These surrogates were employed to probe the generated reactive oxygen species and the subsequent mechanism of N-dealkylation to distinguish between the HAT and SET mechanisms. In addition to the expected N-demethylation of the product aniline, 2,3,4,5,6-pentafluoro-N,N-dimethylaniline N-oxide (PFDMAO) was found to be capable of N-dealkylating both N,N-dimethylaniline (DMA) and N-cyclopropyl-N-methylaniline (CPMA). Rate comparisons of the N-demethylation of DMA supported by PFDMAO show a 27-fold faster rate than when supported by N,N-dimethylaniline N-oxide (DMAO). Whereas intermolecular kinetic isotope effects were masked, intramolecular measurements showed values reflective of those seen previously in DMAO- and the native NADPH/O(2)-supported systems (2.33 and 2.8 for the N-demethylation of PFDMA and DMA from the PFDMAO system, respectively). PFDMAO-supported N-dealkylation of CPMA led to the ring-intact product N-cyclopropylaniline (CPA), similar to that seen with the native system. The formation of CPA argues against a SET mechanism in favor of a P450-like HAT mechanism. We suggest that the similarity of KIEs, in addition to the formation of the ring-intact CPA, argues for a similar mechanism of Compound I (Cpd I) formation followed by HAT for N-dealkylation by the native and N-oxide-supported systems and demonstrate the ability of the N-oxide-generated oxidant to act as an accurate mimic of the native P450 oxidant.

2-aminoimidazole amino acids as inhibitors of the binuclear manganese metalloenzyme human arginase I

Arginase, a key metalloenzyme of the urea cycle that converts L-arginine into L-ornithine and urea, is presently considered a pharmaceutical target for the management of diseases associated with aberrant l-arginine homeostasis, such as asthma, cardiovascular diseases, and erectile dysfunction. We now report the design, synthesis, and evaluation of a series of 2-aminoimidazole amino acid inhibitors in which the 2-aminoimidazole moiety serves as a guanidine mimetic. These compounds represent a new class of arginase inhibitors. The most potent inhibitor identified in this study, 2-(S)-amino-5-(2-aminoimidazol-1-yl)pentanoic acid (A1P, 10), binds to human arginase I with K(d) = 2 microM and significantly attenuates airways hyperresponsiveness in a murine model of allergic airways inflammation. These findings suggest that 2-aminoimidazole amino acids represent new leads for the development of arginase inhibitors with promising pharmacological profiles.

The effects of type II binding on metabolic stability and binding affinity in cytochrome P450 CYP3A4

One goal in drug design is to decrease clearance due to metabolism. It has been suggested that a compound's metabolic stability can be increased by incorporation of a sp(2) nitrogen into an aromatic ring. Nitrogen incorporation is hypothesized to increase metabolic stability by coordination of nitrogen to the heme-iron (termed type II binding). However, questions regarding binding affinity, metabolic stability, and how metabolism of type II binders occurs remain unanswered. Herein, we use pyridinyl quinoline-4-carboxamide analogs to answer these questions. We show that type II binding can have a profound influence on binding affinity for CYP3A4, and the difference in binding affinity can be as high as 1200-fold. We also find that type II binding compounds can be extensively metabolized, which is not consistent with the dead-end complex kinetic model assumed for type II binders. Two alternate kinetic mechanisms are presented to explain the results. The first involves a rapid equilibrium between the type II bound substrate and a metabolically oriented binding mode. The second involves direct reduction of the nitrogen-coordinated heme followed by oxygen binding.

Decreased susceptibility of the cytochrome P450 2B6 variant K262R to inhibition by several clinically important drugs

Cytochrome P450 (P450) 2B6 metabolizes a number of clinically relevant drugs and is one of the most highly polymorphic human P450 enzymes, with the Lys(262)-->Arg substitution being especially common in several genetic variants. Therefore, K262R (2B6*4) was created in the CYP2B6dH background (N-terminal-modified and C-terminal His-tagged) and expressed in Escherichia coli. The recombinant CYP2B6dH and K262R were purified and studied to investigate the effect of the Lys(262)-->Arg substitution with six of the most potent drug inhibitors of CYP2B6, namely, clopidogrel, clotrimazole, itraconazole, raloxifene, sertraline, and ticlopidine. K262R showed a >3-fold increase in the K(i) values with clopidogrel, itraconazole, and raloxifene and approximately 6-fold increase in K(i) with sertraline compared with CYP2B6dH. Likewise, K262R showed 2-, 4-, and >20-fold higher K(s) values than CYP2B6dH with clopidogrel, sertraline, and itraconazole, respectively. In contrast, when tested with several known type II inhibitors of CYP2B enzymes, K262R showed a 10-fold lower IC(50) with 4-(phenyl)pyridine and approximately 2-fold lower IC(50) with 4-(4-nitrobenzyl)pyridine or 1-(4-phenyl)benzylimidazole than CYP2B6dH. Subsequent analysis predicted possible in vivo drug-drug interactions between the CYP2B6 substrate efavirenz and drug inhibitors clopidogrel, clotrimazole, itraconazole, sertraline, and ticlopidine. Furthermore, Q172H/K262R (2B6*6), which is the most common genetic variant of CYP2B6 harboring K262R, was created in CYP2B6dH, expressed, purified, and characterized for inhibition. Q172H/K262R showed a >6-fold increase in K(i) with sertraline and clopidogrel compared with CYP2B6dH. The results suggest that individuals, especially homozygotes, with the 2B6*4 or 2B6*6 allele might be less susceptible to drug interactions resulting from P450 inhibition.

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.

Role of biotransformation studies in minimizing metabolism-related liabilities in drug discovery

Metabolism-related liabilities continue to be a major cause of attrition for drug candidates in clinical development. Such problems may arise from the bioactivation of the parent compound to a reactive metabolite capable of modifying biological materials covalently or engaging in redox-cycling reactions leading to the formation of other toxicants. Alternatively, they may result from the formation of a major metabolite with systemic exposure and adverse pharmacological activity. To avert such problems, biotransformation studies are becoming increasingly important in guiding the refinement of a lead series during drug discovery and in characterizing lead candidates prior to clinical evaluation. This article provides an overview of the methods that are used to uncover metabolism-related liabilities in a pre-clinical setting and offers suggestions for reducing such liabilities via the modification of structural features that are used commonly in drug-like molecules.

Cytochrome P450 expression and regulation in CYP3A4/CYP2D6 double transgenic humanized mice

Analysis of the developmental and sexual expression of cytochrome P450 drug-metabolizing enzymes is impeded by multiple and varied external factors that influence its regulation. In the present study, a CYP2D6/CYP3A4-double transgenic (Tg-CYP2D6/CYP3A4) mouse model was employed to investigate hepatic CYP2D6 and CYP3A4 ontogeny and sexual dimorphism. Both age and sex have considerable effects on hepatic CYP3A4 protein expression in 3- to 8-week-old transgenic mice, whereas neither factor alters CYP2D6 content. Constitutive CYP2D6 expression resulted in 2- to 3-fold higher dextromethorphan O-demethylase activity in Tg-CYP2D6/CYP3A4 mouse liver microsomes compared with wild-type mice. In contrast, expression of CYP3A4 in transgenic mouse livers did not increase dextromethorphan N-demethylase and midazolam 1'-hydroxylase activities. Pretreatment with pregnenolone 16alpha-carbonitrile (PCN) and 1,4-bis-2-(3, 5-dichloropyridyloxy)-benzene (TCPOBOP) elevated CYP3A4 expression in double transgenic mice. Interestingly, induction of hepatic CYP3A4 was greater in females than age- and treatment-matched males. Consequently, the increase in midazolam 1'-hydroxylase activity was markedly higher in 8-week-old female mice than in corresponding males (8-fold versus 6-fold for PCN treatment and 6-fold versus 5-fold for TCPOBOP). Furthermore, increases in testosterone 6beta-hydroxylase activity after CYP3A induction were relatively lower compared with those in midazolam 1'-hydroxylation for age-, sex-, and treatment-matched mice. The difference in CYP3A4 expression and induction between male and female mice suggests that women may be more susceptible to CYP3A4-mediated drug-drug interactions, and the extent of drug-drug interactions could be substrate dependent.

Mechanisms of cytochrome P450 substrate oxidation: MiniReview

Cytochrome P450 (P450) enzymes catalyze a variety of oxidation and some reduction reactions, collectively involving thousands of substrates. A general chemical mechanism can be used to rationalize most of the oxidations and involves a perfenyl intermediate (FeO3+) and odd-electron chemistry, i.e. abstraction of a hydrogen atom or electron followed by oxygen rebound and sometimes rearrangement. This general mechanism can explain carbon hydroxylation, heteroatom oxygenation and dealkylation, epoxidation, desaturation, heme destruction, and other reactions. Another approach to understanding catalysis involves analysis of the more general catalytic cycle, including substrate specificity, because complex patterns of cooperativity are observed with several P450s. Some of the complexity is due to slow conformational changes in the proteins that occur on the same timescale as other steps.

Ligand binding to cytochrome P450 3A4 in phospholipid bilayer nanodiscs: the effect of model membranes

The membrane-bound protein cytochrome P450 3A4 (CYP3A4) is a major drug-metabolizing enzyme. Most studies of ligand binding by CYP3A4 are currently carried out in solution, in the absence of a model membrane. Therefore, there is little information concerning the membrane effects on CYP3A4 ligand binding behavior. Phospholipid bilayer Nanodiscs are a novel model membrane system derived from high density lipoprotein particles, whose stability, monodispersity, and consistency are ensured by their self-assembly. We explore the energetics of four ligands (6-(p-toluidino)-2-naphthalenesulfonic acid (TNS), alpha-naphthoflavone (ANF), miconazole, and bromocriptine) binding to CYP3A4 incorporated into Nanodiscs. Ligand binding to Nanodiscs was monitored by a combination of environment-sensitive ligand fluorescence and ligand-induced shifts in the fluorescence of tryptophan residues present in the scaffold proteins of Nanodiscs; binding to the CYP3A4 active site was monitored by ligand-induced shifts in the heme Soret band absorbance. The dissociation constants for binding to the active site in CYP3A4-Nanodiscs were 4.0 microm for TNS, 5.8 microm for ANF, 0.45 microm for miconazole, and 0.45 microm for bromocriptine. These values are for CYP3A4 incorporated into a lipid bilayer and are therefore presumably more biologically relevant that those measured using CYP3A4 in solution. In some cases, affinity measurements using CYP3A4 in Nanodiscs differ significantly from solution values. We also studied the equilibrium between ligand binding to CYP3A4 and to the membrane. TNS showed no marked preference for either environment; ANF preferentially bound to the membrane, and miconazole and bromocriptine preferentially bound to the CYP3A4 active site.

Contribution of itraconazole metabolites to inhibition of CYP3A4 in vivo

Itraconazole (ITZ) is metabolized in vitro to three inhibitory metabolites: hydroxy-itraconazole (OH-ITZ), keto-itraconazole (keto-ITZ), and N-desalkyl-itraconazole (ND-ITZ). The goal of this study was to determine the contribution of these metabolites to drug-drug interactions caused by ITZ. Six healthy volunteers received 100 mg ITZ orally for 7 days, and pharmacokinetic analysis was conducted at days 1 and 7 of the study. The extent of CYP3A4 inhibition by ITZ and its metabolites was predicted using this data. ITZ, OH-ITZ, keto-ITZ, and ND-ITZ were detected in plasma samples of all volunteers. A 3.9-fold decrease in the hepatic intrinsic clearance of a CYP3A4 substrate was predicted using the average unbound steady-state concentrations (C(ss,ave,u)) and liver microsomal inhibition constants for ITZ, OH-ITZ, keto-ITZ, and ND-ITZ. Accounting for circulating metabolites of ITZ significantly improved the in vitro to in vivo extrapolation of CYP3A4 inhibition compared to a consideration of ITZ exposure alone.

Role of subunit interactions in P450 oligomers in the loss of homotropic cooperativity in the cytochrome P450 3A4 mutant L211F/D214E/F304W

The contribution of conformational heterogeneity to cooperativity in cytochrome P450 3A4 was investigated using the mutant L211F/D214E/F304W. Initial spectral studies revealed a loss of cooperativity of the 1-pyrenebutanol (1-PB) induced spin shift (S(50)=5.4 microM, n=1.0) but retained cooperativity of alpha-naphthoflavone binding. Continuous variation (Job's titration) experiments showed the existence of two pools of enzyme with different 1-PB binding characteristics. Monitoring of 1-PB binding by fluorescence resonance energy transfer from the substrate to the heme confirmed that the high-affinity site (K(D)=0.3 microM) is retained in at least some fraction of the enzyme, although cooperativity is masked. Removal of apoprotein on a second column increased the high-spin content and restored cooperativity of 1-PB binding and of progesterone and testosterone 6beta-hydroxylation. The loss of cooperativity in the mutant is, therefore, mediated by the interaction of holo- and apo-P450 in mixed oligomers.

Multiple Sequential Steps Involved in the Binding of Inhibitors to Cytochrome P450 3A4

Cytochrome P450 (P450) 3A4 is an extensively studied human enzyme involved in the metabolism of >50% of drugs. The mechanism of the observed homotropic and heterotropic cooperativity in P450 3A4-catalyzed oxidations is not well understood, and together with the cooperative behavior, a detailed understanding of interaction of drug inhibitors with P450 3A4 is important in predicting clinical drug-drug interactions. The interactions of P450 3A4 with several structurally diverse inhibitors were investigated using both kinetic and thermodynamic approaches to resolve the steps involved in binding of these ligands. The results of pre-steady-state absorbance and fluorescence experiments demonstrate that inhibitor binding is clearly a multistep process, even more complex than the binding of substrates. Based on spectrophotometric equilibrium binding titrations as well as isothermal titration calorimetry experiments, the stoichiometry of binding appears to be 1:1 in the concentration ranges studied. Using a sequential-mixing stopped-flow approach, we were also able to show that the observed multiphasic binding kinetics is the result of sequential events as opposed to the existence of multiple enzyme populations in dynamic equilibrium that interact with ligands at different rates. We propose a three-step minimal model for inhibitor binding, developed with kinetic simulations, consistent with our previously reported model for the binding of substrates, although it is possible that even more steps are involved.