Visible spectra of type II cytochrome P450-drug complexes: evidence that "incomplete" heme coordination is common

Emerging trends in small molecule inhibitors targeting aldosterone synthase: A new paradigm in cardiovascular disease treatment

Aldosterone synthase (CYP11B2) is the rate-limiting enzyme in aldosterone production. In recent years, CYP11B2 has become an appealing target for treating conditions associated with excess aldosterone, such as hypertension, heart failure, and cardiometabolic diseases. Several small-molecule inhibitors of CYP11B2 have demonstrated efficacy in both preclinical studies and clinical trials. Among them, the tetrahydroisoquinoline derivative Baxdrostat has entered clinical trial phases and demonstrated efficacy in treating patients with hypertension. However, the high homology (>93 %) between CYP11B2 and steroid-11β-hydroxylase (CYP11B1), which catalyzes cortisol production, implies that insufficient drug specificity can lead to severe side effects. Developing selective inhibitors for CYP11B2 remains a considerable challenge that requires ongoing attention. This review summarizes recent research progress on small-molecule inhibitors targeting CYP11B2, focusing on structure-activity relationships (SAR) and structural optimization. It discusses strategies for enhancing the specificity and inhibitory activity of inhibitors, while also exploring potential applications and future prospects for CYP11B2 inhibitors, providing a theoretical foundation for developing the new generation of CYP11B2-targeted medications.

On the Possible Effect of Phytic Acid (Myo-Inositol Hexaphosphoric Acid, IP6) on Cytochromes P450 and Systems of Xenobiotic Metabolism in Different Hepatic Models

As compounds of natural origin enter human body, it is necessary to investigate their possible interactions with the metabolism of drugs and xenobiotics in general, namely with the cytochrome P450 (CYP) system. Phytic acid (myo-inositol hexaphosphoric acid, IP6) is mainly present in plants but is also an endogenous compound present in mammalian cells and tissues. It has been shown to exhibit protective effect in many pathological conditions. For this paper, its interaction with CYPs was studied using human liver microsomes, primary human hepatocytes, the HepG2 cell line, and molecular docking. Docking experiments and absorption spectra demonstrated the weak ability of IP6 to interact in the heme active site of CYP1A. Molecular docking suggested that IP6 preferentially binds to the protein surface, whereas binding to the active site of CYP1A2 was found to be less probable. Subsequently, we investigated the ability of IP6 to modulate the metabolism of xenobiotics for both the mRNA expression and enzymatic activity of CYP1A enzymes. Our findings revealed that IP6 can slightly modulate the mRNA levels and enzyme activity of CYP1A. However, thanks to the relatively weak interactions of IP6 with CYPs, the chances of the mechanisms of clinically important drug-drug interactions involving IP6 are low.

The binding of nitrogen-donor ligands to the ferric and ferrous forms of cytochrome P450 enzymes

The cytochrome P450 superfamily of heme-thiolate monooxygenase enzymes can catalyse various oxidation reactions. The addition of a substrate or an inhibitor ligand induces changes in the absorption spectrum of these enzymes and UV-visible (UV-vis) absorbance spectroscopy is the most common and readily available technique used to interrogate their heme and active site environment. Nitrogen-containing ligands can inhibit the catalytic cycle of heme enzymes by interacting with the heme. Here we evaluate the binding of imidazole and pyridine-based ligands to the ferric and ferrous forms of a selection of bacterial cytochrome P450 enzymes using UV-visible absorbance spectroscopy. The majority of these ligands interact with the heme as one would expect for type II nitrogen directly coordinated to a ferric heme-thiolate species. However, the spectroscopic changes observed in the ligand-bound ferrous forms indicated differences in the heme environment across these P450 enzyme/ligand combinations. Multiple species were observed in the UV-vis spectra of the ferrous ligand-bound P450s. None of the enzymes gave rise to the isolation of a single species with a Soret band at ∼442-447 nm, indicative of a 6-coordinate ferrous thiolate species with a nitrogen-donor ligand. A ferrous species with Soret band at ∼427 nm coupled with an α-band of increased intensity was observed with the imidazole ligands. With some enzyme-ligand combinations reduction resulted in breaking of the iron‑nitrogen bond yielding a 5-coordinate high-spin ferrous species. In other instances, the ferrous form was readily oxidised back to the ferric form on addition of the ligand.

Structural insights into the effects of glycerol on ligand binding to cytochrome P450

New antitubercular drugs are vital due to the spread of resistant strains. Carbethoxyhexyl imidazole (CHImi) inhibits cytochrome P450 CYP124, which is a steroid-metabolizing enzyme that is important for the survival of Mycobacterium tuberculosis in macrophages. The available crystal structure of the CYP124-CHImi complex reveals two glycerol molecules in the active site. A 1.15 Å resolution crystal structure of the glycerol-free CYP124-CHimi complex reported here shows multiple conformations of CHImi and the CYP124 active site which were previously restricted by glycerol. Complementary molecular dynamics simulations show coherence of the ligand and enzyme conformations. Spectrophotometric titration confirmed the influence of glycerol on CHImi binding: the affinity decreases more than tenfold in glycerol-containing buffer. In addition, it also showed that glycerol has a similar effect on other azole and triazole CYP124 ligands. Together, these data show that glycerol may compromise structural-functional studies and impede rational drug-design campaigns.

Interaction of CYP3A4 with Rationally Designed Ritonavir Analogues: Impact of Steric Constraints Imposed on the Heme-Ligating Group and the End-Pyridine Attachment

Controlled inhibition of drug-metabolizing cytochrome P450 3A4 (CYP3A4) is utilized to boost bioavailability of anti-viral and immunosuppressant pharmaceuticals. We investigate structure–activity relationships (SARs) in analogues of ritonavir, a potent CYP3A4 inhibitor marketed as pharmacoenhancer, to determine structural elements required for potent inhibition and whether the inhibitory potency can be further improved via a rational structure-based design. This study investigated eight (series VI) inhibitors differing in head- and end-moieties and their respective linkers. SAR analysis revealed the multifactorial regulation of inhibitory strength, with steric constraints imposed on the tethered heme-ligating moiety being a key factor. Minimization of these constraints by changing the linkers’ length/flexibility and N-heteroatom position strengthened heme coordination and markedly improved binding and/or inhibitory strength. Impact of the end-pyridine attachment was not uniform due to influence of other determinants controlling the ligand-binding mode. This interplay between pharmacophoric determinants and the end-group enlargement can be used for further inhibitor optimization.

Cytochrome P450 168A1 from Pseudomonas aeruginosa is involved in the hydroxylation of biologically relevant fatty acids

The cytochrome P450 CYP168A1 from Pseudomonas aeruginosa was cloned and expressed in Escherichia coli followed by purification and characterization of function. CYP168A1 is a fatty acid hydroxylase that hydroxylates saturated fatty acids, including myristic (0.30 min-1), palmitic (1.61 min-1) and stearic acids (1.24 min-1), at both the ω-1- and ω-2-positions. However, CYP168A1 only hydroxylates unsaturated fatty acids, including palmitoleic (0.38 min-1), oleic (1.28 min-1) and linoleic acids (0.35 min-1), at the ω-1-position. CYP168A1 exhibited a catalytic preference for palmitic, oleic and stearic acids as substrates in keeping with the phosphatidylcholine-rich environment deep in the lung that is colonized by P. aeruginosa.

Spectroscopic observations of β-eudesmol binding to human cytochrome P450 isoforms 3A4 and 1A2, but not to isoforms 2C9, 2C19 and 2D6

β-Eudesmol is a sesquiterpenoid component o Atractylodes lancea with cytotoxic activity against cholangiocarcinoma. Its lipophilic nature makes β-eudesmol a likely substrate of human cytochrome P450 (P450) enzymes.Using ligand-binding difference spectroscopy, the affinities of this compound to recombinant CYP1A2, CYP2C9, CYP2C19, CYP2D6 and CYP3A4 were investigated in Escherichia coli membrane preparations.CYP3A4 showed a type I spectral change, with a binding constant Ks of 77 ± 23 (mean ± SD) μM at 0.5 μM P450 (Ks/[P450] ≈ 155). The reference substrate testosterone and the inhibitor fluconazole bound to the enzyme with apparent affinities of 86 ± 4 μM (type I) and 21 μM (type II), respectively. β-Eudesmol was bound to CYP3A4 in a non-cooperative manner (Hill coefficient n ≈ 0.8). CYP1A2 showed reverse type I difference spectra with either β-eudesmol or caffeine. The CYP1A2 affinity for β-eudesmol was higher (0.23 mM) than for caffeine (0.37 mM) but lower than for phenacetin (0.11 mM, type I). β-Eudesmol did not bind significantly to CYP2C9, CYP2C19, and CYP2D6.Confirmation of metabolic activity and studies on the involvement of other human P450 isoforms studies are required. Double-beam spectrometry is needed to validate Ks measurements made with a plate reader.

To Be, or Not to Be, an Inhibitor: A Comparison of Azole Interactions with and Oxidation by a Cytochrome P450 Enzyme

The cytochrome P450 (CYP) superfamily of heme monooxygenases is involved in a range of important chemical biotransformations across nature. Azole-containing molecules have been developed as drugs that bind to the heme center of these enzymes, inhibiting their function. The optical spectrum of CYP enzymes after the addition of these inhibitors is used to assess how the molecules bind. Here we use the bacterial CYP199A4 enzyme, from Rhodopseudomonas palustris HaA2, to compare how imidazolyl and triazolyl inhibitors bind to ferric and ferrous heme. 4-(Imidazol-1-yl)benzoic acid induced a red shift in the Soret wavelength (424 nm) in the ferric enzyme along with an increase and a decrease in the intensities of the δ and α bands, respectively. 4-(1H-1,2,4-Triazol-1-yl)benzoic acid binds to CYP199A4 with a 10-fold lower affinity and induces a smaller red shift in the Soret band. The crystal structures of CYP199A4 with these two inhibitors confirmed that these differences in the optical spectra were due to coordination of the imidazolyl ligand to the ferric Fe, but the triazolyl inhibitor interacts with, rather than displaces, the ferric aqua ligand. Additional water molecules were present in the active site of 4-(1H-1,2,4-triazol-1-yl)benzoic acid-bound CYP199A4. The space required to accommodate these additional water molecules in the active site necessitates changes in the position of the hydrophobic phenylalanine 298 residue. Upon reduction of the heme, the imidazole-based inhibitor Fe-N ligation was not retained. A 5-coordinate heme was also the predominant species in 4-(1H-1,2,4-triazol-1-yl)benzoic acid-bound ferrous CYP199A4, but there was an obvious shoulder at 447 nm indicative of some degree of Fe-N coordination. Rather than inhibit CYP199A4, 4-(imidazol-1-yl)benzoic acid was a substrate and was oxidized to generate a metabolite derived from ring opening of the imidazolyl ring: 4-[[2-(formylamino)acetyl]amino]benzoic acid.

Biophysical Techniques for Distinguishing Ligand Binding Modes in Cytochrome P450 Monooxygenases

The cytochrome P450 superfamily of heme monooxygenases catalyzes important chemical reactions across nature. The changes in the optical spectra of these enzymes, induced by the addition of substrates or inhibitors, are critical for assessing how these molecules bind to the P450, enhancing or inhibiting the catalytic cycle. Here we use the bacterial CYP199A4 enzyme (Uniprot entry Q2IUO2), from Rhodopseudomonas palustris HaA2, and a range of substituted benzoic acids to investigate different binding modes. 4-Methoxybenzoic acid elicits an archetypal type I spectral response due to a ≥95% switch from the low- to high-spin state with concomitant dissociation of the sixth aqua ligand. 4-(Pyridin-3-yl)- and 4-(pyridin-2-yl)benzoic acid induced different type II ultraviolet-visible (UV-vis) spectral responses in CYP199A4. The former induced a greater red shift in the Soret wavelength (424 nm vs 422 nm) along with a larger overall absorbance change and other differences in the α-, β-, and δ-bands. There were also variations in the ferrous UV-vis spectra of these two substrate-bound forms with a spectrum indicative of Fe-N bond formation with 4-(pyridin-3-yl)benzoic acid. The crystal structures of CYP199A4, with the pyridinyl compounds bound, revealed that while the nitrogen of 4-(pyridin-3-yl)benzoic acid is coordinated to the heme, with 4-(pyridin-2-yl)benzoic acid an aqua ligand remains. Continuous wave and pulse electron paramagnetic resonance data in frozen solution revealed that the substrates are bound in the active site in a form consistent with the crystal structures. The redox potential of each CYP199A4-substrate combination was measured, allowing correlation among binding modes, spectroscopic properties, and the observed biochemical activity.

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.

Structure-Activity Relationships of Rationally Designed Ritonavir Analogues: Impact of Side-Group Stereochemistry, Headgroup Spacing, and Backbone Composition on the Interaction with CYP3A4

In a continuing effort to identify structural attributes required for strong binding and potent inhibition of human drug-metabolizing CYP3A4, we designed ten ritonavir-like analogues differing in the side-group stereochemistry, backbone atomic composition, and headgroup spacing. All analogues had pyridine and tert-butyloxycarbonyl (Boc) as the heme-ligating head and tail groups, respectively, phenyl side groups, and either a methyl- or ethyl-pyridyl linker. Each linker subseries had S/ R, R/ S, R/ R, and S/S side-group conformers (4a-d and 4e-h, respectively), and one S/S stereoisomer with the backbone S-to-N-heteroatom substitution (6a and 6b). To elucidate structure-activity relationships, ligand-dependent changes in optical spectra, dissociation constant ( K_s), inhibitory potency (IC₅₀), thermostability, and heme ligation and reduction kinetics were analyzed. Comparison of the subseries and individual compounds showed that CYP3A4 only weakly discriminates between side-group configurations, associates more tightly with the pyridyl-ethyl-linker analogues, and strongly disfavors the N-containing backbone. K_s and IC₅₀ for the pyridyl-ethyl R/ R conformer, 4g, were the lowest and close to those for ritonavir: 0.04 and 0.31 μM versus 0.02 and 0.13 μM, respectively. Determination of the X-ray structures of the inhibitory complexes was critical for experimental data interpretation, especially for the uniquely oriented 4a and 4e. Based on structural analysis, we conclude that, for this series of analogues, the ligand-mediated interactions near the heme are dominant and define the binding mode and that fine-tuning of these interactions as well as the backbone spacing could further improve the affinity and inhibitory strength.

Critical Differences between Enzyme Sources in Sensitivity to Detect Time-Dependent Inactivation of CYP2C8

Clopidogrel acyl-β-d-glucuronide is a mechanism-based inhibitor of cytochrome P450 2C8 in human liver microsomes (HLMs). However, time-dependent inactivation (TDI) of CYP2C8 could not be detected in an earlier study in human recombinant CYP2C8 (Supersomes). Here, we investigate whether different enzyme sources exhibit differences in detection of CYP2C8 TDI under identical experimental conditions. Inactivation of CYP2C8 by amiodarone (100 μM), clopidogrel acyl-β-d-glucuronide (100 μM), gemfibrozil 1-O-β-glucuronide (100 μM), and phenelzine (100 μM) was investigated in HLMs and three recombinant human CYP2C8 preparations (Supersomes, Bactosomes, and EasyCYP Bactosomes) using amodiaquine N-deethylation as the marker reaction. Furthermore, the inactivation kinetics of CYP2C8 by clopidogrel glucuronide (5-250 μM) was determined in Supersomes and Bactosomes. Amiodarone caused weak TDI in all enzyme preparations tested, while the extent of inactivation by clopidogrel glucuronide, gemfibrozil glucuronide, and phenelzine varied markedly between preparations, and even different Supersome lots. Both glucuronides caused strong inactivation of CYP2C8 in HLMs, Bactosomes and in one Supersome lot (>50% inhibition), but significant inactivation could not be reliably detected in other Supersome lots or EasyCYP Bactosomes. In Bactosomes, the concentration producing half of kinact (KI) and maximal inactivation rate (kinact) of clopidogrel glucuronide (14 μM and 0.054 minute-1) were similar to those determined previously in HLMs. Phenelzine caused strong inactivation of CYP2C8 in one Supersome lot (91% inhibition) but not in HLMs or other recombinant CYP2C8 preparations. In conclusion, different enzyme sources and different lots of the same recombinant enzyme preparation are not equally sensitive to detect inactivation of CYP2C8, suggesting that recombinant CYPs should be avoided when identifying mechanism-based inhibitors.

Interaction of Human Drug-Metabolizing CYP3A4 with Small Inhibitory Molecules

Binding of small inhibitory compounds to human cytochrome P450 3A4 (CYP3A4) could interfere with drug metabolism and lead to drug-drug interactions, the underlying mechanism of which is not fully understood due to insufficient structural information. This study investigated the interaction of recombinant CYP3A4 with a nonspecific inhibitor metyrapone, antifungal drug fluconazole, and protease inhibitor phenylmethanesulfonyl fluoride (PMSF). Metyrapone and fluconazole are classic type II ligands that inhibit CYP3A4 with medium strength by ligating to the heme iron, whereas PMSF, lacking the heme-ligating moiety, acts as a weak type I ligand and inhibitor of CYP3A4. High-resolution crystal structures revealed that the orientation of metyrapone is similar but not identical to that in the previously reported 1W0G model, whereas the flexible fluconazole adapts a conformer markedly different from that observed in the target CYP51 enzymes, which could explain its high potential for cross-reactivity. Besides hydrophobic and aromatic interactions with the heme and active site residues, both drugs establish water-mediated contacts that stabilize the inhibitory complexes. PMSF also binds near the catalytic center, with the phenyl group parallel to the heme. However, it does not displace the water ligand and is held in place via strong H-bonds formed by the sulfofluoride moiety with Ser119 and Arg212. Collectively, our data suggest that PMSF might have multiple binding sites and likely occupies the high-affinity site in the crystal structure. Moreover, its hydrolysis product, phenylmethanesulfonic acid, can also access and be retained in the CYP3A4 active site. Therefore, to avoid experimental artifacts, PMSF should be excluded from purification and assay solutions.

Ruthenium-containing P450 inhibitors for dual enzyme inhibition and DNA damage

Cytochrome P450s are key players in drug metabolism, and overexpression in tumors is associated with significant resistance to many medicinal agents. Consequently, inhibition of P450s could serve as a strategy to restore drug efficacy. However, the widespread expression of P450s throughout the human body and the critical roles they play in various biosynthetic pathways motivates the development of P450 inhibitors capable of controlled local administration. Ruthenium complexes containing P450 inhibitors as ligands were synthesized in order to develop pro-drugs that can be triggered to release the inhibitors in a spatially and temporally controlled fashion. Upon light activation the compounds release ligands that directly bind and inhibit P450 enzymes, while the ruthenium center is able to directly damage DNA.

Increased Phenacetin Oxidation upon the L382V Substitution in Cytochrome P450 1A2 is Associated with Altered Substrate Binding Orientation

Leucine382 of cytochrome P450 1A2 (CYP1A2) plays an important role in binding and O-dealkylation of phenacetin, with the L382V mutation increasing substrate oxidation (Huang and Szklarz, 2010, Drug Metab. Dispos. 38:1039–1045). This was attributed to altered substrate binding orientation, but no direct experimental evidence had been available. Therefore, in the current studies, we employed nuclear magnetic resonance (NMR) longitudinal (T₁) relaxation measurements to investigate phenacetin binding orientations within the active site of CYP1A2 wild type (WT) and mutants. Paramagnetic relaxation time (T_1P) for each proton of phenacetin was calculated from the T₁ value obtained from the enzymes in ferric and ferrous-CO state in the presence of phenacetin, and used to model the orientation of phenacetin in the active site. All aromatic protons of phenacetin were nearly equidistant from the heme iron (6.34–8.03 Å). In contrast, the distance between the proton of the –OCH₂– group, which is abstracted during phenacetin oxidation, and the heme iron, was much shorter in the L382V (5.93 Å) and L382V/N312L (5.96 Å) mutants compared to the N312L mutant (7.84 Å) and the wild type enzyme (6.55 Å), consistent with modeling results. These studies provide direct evidence for the molecular mechanism underlying increased oxidation of phenacetin upon the L382V mutation.

CW EPR parameters reveal cytochrome P450 ligand binding modes

Cytochrome P450 (CYP) monoxygenses utilize heme cofactors to catalyze oxidation reactions. They play a critical role in metabolism of many classes of drugs, are an attractive target for drug development, and mediate several prominent drug interactions. Many substrates and inhibitors alter the spin state of the ferric heme by displacing the heme's axial water ligand in the resting enzyme to yield a five-coordinate iron complex, or they replace the axial water to yield a nitrogen-ligated six-coordinate iron complex, which are traditionally assigned by UV-vis spectroscopy. However, crystal structures and recent pulsed electron paramagnetic resonance (EPR) studies find a few cases where molecules hydrogen bond to the axial water. The water-bridged drug-H₂O-heme has UV-vis spectra similar to nitrogen-ligated, six-coordinate complexes, but are closer to "reverse type I" complexes described in older liteature. Here, pulsed and continuous wave (CW) EPR demonstrate that water-bridged complexes are remarkably common among a range of nitrogenous drugs or drug fragments that bind to CYP3A4 or CYP2C9. Principal component analysis reveals a distinct clustering of CW EPR spectral parameters for water-bridged complexes. CW EPR reveals heterogeneous mixtures of ligated states, including multiple directly-coordinated complexes and water-bridged complexes. These results suggest that water-bridged complexes are under-represented in CYP structural databases and can have energies similar to other ligation modes. The data indicates that water-bridged binding modes can be identified and distinguished from directly-coordinated binding by CW EPR.

Spectroscopic studies of the cytochrome P450 reaction mechanisms

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

Dissociation Constants of Cytochrome P450 2C9/Cytochrome P450 Reductase Complexes in a Lipid Bilayer Membrane Depend on NADPH: A Single-Protein Tracking Study

Cytochrome P450-reductase (CPR) is a versatile NADPH-dependent electron donor located in the cytoplasmic side of the endoplasmic reticulum. It is an electron transferase that is able to deliver electrons to a variety of membrane-bound oxidative partners, including the drug-metabolizing enzymes of the cytochrome P450s (P450). CPR is also stoichiometrically limited compared to its oxidative counterparts, and hypotheses have arisen about possible models that can overcome the stoichiometric imbalance, including quaternary organization of P450 and diffusion-limited models. Described here are results from a single-protein tracking study of fluorescently labeled CPR and cytochrome P450 2C9 (CYP2C9) molecules in which stochastic analysis was used to determine the dissociation constants of CPR/CYP2C9 complexes in a lipid bilayer membrane for the first time. Single-protein trajectories demonstrate the transient nature of these CPR-CYP2C9 interactions, and the measured K_d values are highly dependent on the redox state of CPR. It is shown that CPR^ox/CYP2C9 complexes have a much higher dissociation constant than CPR^2-/CYP2C9 or CPR^4-/CYP2C9 complexes, and a model is presented to account for these results. An Arrhenius analysis of diffusion constants was also carried out, demonstrating that the reduced forms of CPR and CYP2C9 interact differently with the biomimetic ER and may, in addition to protein conformational changes, contribute to the observed NADPH-dependent shift in K_d. Finally, it is also shown that the CPR^ox/CYP2C9 affinity depends on the nature of the ligand, being higher when a substrate is bound, compared to an inhibitor.

Interaction of Azole-Based Environmental Pollutants with the Coelomic Hemoglobin from Amphitrite ornata: A Molecular Basis for Toxicity

The toxicities of azole pollutants that have widespread agricultural and industrial uses are either poorly understood or unknown, particularly with respect to how infaunal organisms are impacted by this class of persistent organic pollutant. To identify a molecular basis by which azole compounds may have unforeseen toxicity on marine annelids, we examine here their impact on the multifunctional dehaloperoxidase (DHP) hemoglobin from the terebellid polychaete Amphitrite ornata. Ultraviolet-visible and resonance Raman spectroscopic studies showed an increase in the six-coordinate low-spin heme population in DHP isoenzyme B upon binding of imidazole, benzotriazole, and benzimidazole (K_d values of 52, 82, and 110 μM, respectively), suggestive of their direct binding to the heme-Fe. Accordingly, atomic-resolution X-ray crystal structures, supported by computational studies, of the DHP B complexes of benzotriazole (1.14 Å), benzimidazole (1.08 Å), imidazole (1.08 Å), and indazole (1.12 Å) revealed two ligand binding motifs, one with direct ligand binding to the heme-Fe, and another in which the ligand binds in the hydrophobic distal pocket without coordinating the heme-Fe. Taken together, the results demonstrate a new mechanism by which azole pollutants can potentially disrupt hemoglobin function, thereby improving our understanding of their impact on infaunal organisms in marine and aquatic environments.

Separation of enzymatic functions and variation of spin state of rice allene oxide synthase-1 by mutation of Phe-92 and Pro-430

Rice allene oxide synthase-1 mutants carrying F92L, P430A or F92L/P430A amino acid substitution mutations were constructed, recombinant mutant and wild type proteins were purified and their substrate preference, UV-vis spectra and heme iron spin state were characterized. The results show that the hydroperoxide lyase activities of F92L and F92L/P430A mutants prefer 13-hydroperoxy substrate to other hydroperoxydienoic acids or hydroperoxytrienoic acids. The Soret maximum was completely red-shifted in P430A and F92L/P430A mutants, but it was partially shifted in the F92L mutant. ESR spectral data showed that wild type, F92L and P430A mutants occupied high and low spin states, while the F92L/P430A mutant occupied only low spin state. The extent of the red shift of the Soret maximum increased as the population of low spin heme iron increased, suggesting that the spectral shift reflects the high to low transition of heme iron spin state in rice allene oxide synthase-1. Relative to wild type allene oxide synthase-1, the hydroperoxide lyase activities of F92L and F92L/P430A are less sensitive to inhibition by imidazole with (13S or 9S)-hydroperoxydienoic acid as substrate and more sensitive than wild type with (13S)-hydroperoxytrienoic acid as substrate. Our results suggest that hydroperoxydienoic acid is the preferred substrate for the hydroperoxide lyase activity and (13S)-hydroperoxytrienoic acid is the preferred substrate for allene oxide synthase activity of allene oxide synthase-1.

Rapid kinetic methods to dissect steroidogenic cytochrome P450 reaction mechanisms

All cytochrome P450 enzyme reactions involve a catalytic cycle with several discreet physical or chemical steps. This cycle ends with the formation of the reactive heme iron-oxygen complex, which oxygenates substrate. While the steps might be very similar for each P450 enzyme, the rates of each step varies tremendously for each enzyme and sometimes even for different reactions catalyzed by the same enzyme. For example, the rate-limiting step for most bacterial P450 enzymes, with turnover numbers over 1000s(-1), is the second electron transfer. In contrast, steroidogenic P450s from eukaryotes catalyze much slower reactions, with turnover numbers of ∼5-250min(-1); therefore, assumptions about kinetic properties for the mammalian P450 enzymes based on the bacterial enzymes are tenuous. In order to dissect the rates for individual steps, special techniques that isolate individual steps and/or single turnovers are required. This article will review the theoretical principles and practical considerations for several of these techniques, with illustrative published examples. The reader should gain an appreciation for the appropriate methods used to interrogate particular steps in the P450 reaction cycle.

Selective Targeting of Heme Protein in Cytochrome P450 and Nitric Oxide Synthase by Diphenyleneiodonium

Cytochrome P450 (CYP) enzymes mediate mixed-function oxidation reactions important in drug metabolism. The aromatic heterocyclic cation, diphenyleneiodonium (DPI), binds flavin in cytochrome P450 reductase and inhibits CYP-mediated activity. DPI also inhibits CYP by directly interacting with heme. Herein, we report that DPI effectively inhibits a number of CYP-related monooxygenase reactions including NADPH oxidase, a microsomal enzyme activity that generates hydrogen peroxide in the absence of metabolizing substrates. Inhibition of monooxygenase by DPI was time and concentration dependent with IC50's ranging from 0.06 to 1.9 μM. Higher (4.6-23.9 μM), but not lower (0.06-1.9 μM), concentrations of DPI inhibited electron flow via cytochrome P450 reductase, as measured by its ability to reduce cytochrome c and mediate quinone redox cycling. Similar results were observed with inducible nitric oxide synthase (iNOS), an enzyme containing a C-terminal reductase domain homologous to cytochrome P450 reductase that mediates reduction of cytochrome c, and an N-terminal heme-thiolate oxygenase domain mediating nitric oxide production. Significantly greater concentrations of DPI were required to inhibit cytochrome c reduction by iNOS (IC50 = 3.5 µM) than nitric oxide production (IC50 = 0.16 µM). Difference spectra of liver microsomes, recombinant CYPs, and iNOS demonstrated that DPI altered heme-carbon monoxide interactions. In the presence of NADPH, DPI treatment of microsomes and iNOS yielded a type II spectral shift. These data indicate that DPI interacts with both flavin and heme in CYPs and iNOS. Increased sensitivity for inhibition of CYP-mediated metabolism and nitric oxide production by iNOS indicates that DPI targets heme moieties within the enzymes.

The Effect of Genetic Polymorphism on the Inhibition of Azole Antifungal Agents Against CYP2C9-Mediated Metabolism

We investigated the effect of cytochrome P450 (CYP) 2C9 polymorphism on the inhibition of methylhydroxylation activity of tolbutamide, a typical CYP2C9 substrate, by triazole antifungal agents, fluconazole and voriconazole. Although the Michaelis constants (Km), maximal velocities (Vmax), and Vmax/Km values for CYP2C9.1 (wild type) and CYP2C9.2 (Arg144Cys) were similar and CYP2C9.3 (Ile359Leu) had a higher Km and a lower Vmax than CYP2C9.1 and CYP2C9.2, the inhibition constants of fluconazole and voriconazole against CYP2C9.2 were lower than that against CYP2C9.1 and CYP2C9.3. These results suggest that more careful administration of azole antifungals to patients with the CYP2C9*2 allele might be required because of the strong inhibitory effects.

Structural Insights into Binding of the Antifungal Drug Fluconazole to Saccharomyces cerevisiae Lanosterol 14α-Demethylase

Infections by fungal pathogens such as Candida albicans and Aspergillus fumigatus and their resistance to triazole drugs are major concerns. Fungal lanosterol 14α-demethylase belongs to the CYP51 class in the cytochrome P450 superfamily of enzymes. This monospanning bitopic membrane protein is involved in ergosterol biosynthesis and is the primary target of azole antifungal drugs, including fluconazole. The lack of high-resolution structural information for this drug target from fungal pathogens has been a limiting factor for the design of modified triazole drugs that will overcome resistance. Here we report the X-ray structure of full-length Saccharomyces cerevisiae lanosterol 14α-demethylase in complex with fluconazole at a resolution of 2.05 Å. This structure shows the key interactions involved in fluconazole binding and provides insight into resistance mechanisms by revealing a water-mediated hydrogen bonding network between the drug and tyrosine 140, a residue frequently found mutated to histidine or phenylalanine in resistant clinical isolates.

Nanoscale electron transport measurements of immobilized cytochrome P450 proteins

Gold nanopillars, functionalized with an organic self-assembled monolayer, can be used to measure the electrical conductance properties of immobilized proteins without aggregation. Measurements of the conductance of nanopillars with cytochrome P450 2C9 (CYP2C9) proteins using conducting probe atomic force microscopy demonstrate that a correlation exists between the energy barrier height between hopping sites and CYP2C9 metabolic activity. Measurements performed as a function of tip force indicate that, when subjected to a large force, the protein is more stable in the presence of a substrate. This agrees with the hypothesis that substrate entry into the active site helps to stabilize the enzyme. The relative distance between hopping sites also increases with increasing force, possibly because protein functional groups responsible for electron transport (ETp) depend on the structure of the protein. The inhibitor sulfaphenazole, in addition to the previously studied aniline, increased the barrier height for electron transfer and thereby makes CYP2C9 reduction more difficult and inhibits metabolism. This suggests that P450 Type II binders may decrease the ease of ETp processes in the enzyme, in addition to occupying the active site.

Drug modulation of water-heme interactions in low-spin P450 complexes of CYP2C9d and CYP125A1

Azoles and pyridines are commonly incorporated into small molecule inhibitor scaffolds that target cytochromes P450 (CYPs) as a strategy to increase drug binding affinity, impart isoform-dependent selectivity, and improve metabolic stability. Optical absorbance spectra of the CYP-inhibitor complex are widely used to infer whether these inhibitors are ligated directly to the heme iron as catalytically inert, low-spin (type II) complexes. Here, we show that the low-spin complex between a drug-metabolizing CYP2C9 variant and 4-(3-phenylpropyl)-1H-1,2,3-triazole (PPT) retains an axial water ligand despite exhibiting elements of "classic" type II optical behavior. Hydrogens of the axial water ligand are observed by pulsed electron paramagnetic resonance (EPR) spectroscopy for both inhibitor-free and inhibitor-bound species and show that inhibitor binding does not displace the axial water. A (15)N label incorporated into PPT is 0.444 nm from the heme iron, showing that PPT is also in the active site. The reverse type I inhibitor, LP10, of CYP125A1 from Mycobacterium tuberculosis, known from X-ray crystal structures to form a low-spin water-bridged complex, is found by EPR and by visible and near-infrared magnetic circular dichroism spectroscopy to retain the axial water ligand in the complex in solution.

Myoglobin-catalyzed intermolecular carbene N-H insertion with arylamine substrates

Engineered variants of the heme-containing protein myoglobin can efficiently catalyze the insertion of α-diazo esters into the N-H bond of arylamines, featuring a combination of high chemoselectivity, elevated turnover numbers, and broad substrate scope.

The clinical candidate VT-1161 is a highly potent inhibitor of Candida albicans CYP51 but fails to bind the human enzyme

The binding and cytochrome P45051 (CYP51) inhibition properties of a novel antifungal compound, VT-1161, against purified recombinant Candida albicans CYP51 (ERG11) and Homo sapiens CYP51 were compared with those of clotrimazole, fluconazole, itraconazole, and voriconazole. VT-1161 produced a type II binding spectrum with Candida albicans CYP51, characteristic of heme iron coordination. The binding affinity of VT-1161 for Candida albicans CYP51 was high (dissociation constant [Kd], ≤ 39 nM) and similar to that of the pharmaceutical azole antifungals (Kd, ≤ 50 nM). In stark contrast, VT-1161 at concentrations up to 86 μM did not perturb the spectrum of recombinant human CYP51, whereas all the pharmaceutical azoles bound to human CYP51. In reconstitution assays, VT-1161 inhibited Candida albicans CYP51 activity in a tight-binding fashion with a potency similar to that of the pharmaceutical azoles but failed to inhibit the human enzyme at the highest concentration tested (50 μM). In addition, VT-1161 (MIC = 0.002 μg ml(-1)) had a more pronounced fungal sterol disruption profile (increased levels of methylated sterols and decreased levels of ergosterol) than the known CYP51 inhibitor voriconazole (MIC = 0.004 μg ml(-1)). Furthermore, VT-1161 weakly inhibited human CYP2C9, CYP2C19, and CYP3A4, suggesting a low drug-drug interaction potential. In summary, VT-1161 potently inhibited Candida albicans CYP51 and culture growth but did not inhibit human CYP51, demonstrating a >2,000-fold selectivity. This degree of potency and selectivity strongly supports the potential utility of VT-1161 in the treatment of Candida infections.

Modelling cytochromes P450 binding modes to predict P450 inhibition, metabolic stability and isoform selectivity

The cytochromes P450 (P450) superfamily is a diverse group of enzymes involved in the metabolism of xenobiotics, whose orientations within the catalytic site can lead to different binding modes, namely productive, nonproductive, and inhibitory. This article collects the most recent approaches that individually study P450- ligand interactions, including a novel in silico technology, developed in the framework of the Human Cytochrome P450 Consortium initiative, that provides reliable in silico predictions of P450 inhibition, metabolic stability and isoform selectivity.

Dual positional substrate specificity of rice allene oxide synthase-1: insight into mechanism of inhibition by type II ligand imidazole

Phylogenetic and amino acid sequence analysis indicated that rice allene oxide synthase-1 (OsAOS1) is CYP74, and is clearly distinct from CYP74B, C and D subfamilies. Regio- and stereo-chemical analysis revealed the dual substrate specificity of OsAOS1 for (cis,trans)-configurational isomers of 13(S)- and 9(S)-hydroperoxyoctadecadienoic acid. GC-MS analysis showed that OsAOS1 converts 13(S)- and 9(S)-hydroperoxyoctadecadi(tri)enoic acid into their corresponding allene oxide. UV-Visible spectral analysis of native OsAOS1 revealed a Soret maximum at 393 nm, which shifted to 424 nm with several clean isobestic points upon binding of OsAOS1 to imidazole. The spectral shift induced by imidazole correlated with inhibition of OsAOS1 activity, implying that imidazole may coordinate to ferric heme iron, triggering a heme-iron transition from high spin state to low spin state. The implications and significance of a putative type II ligand-induced spin state transition in OsAOS1 are discussed. [BMB Reports 2013; 46(3):151-156]

Effect of CYP3A5 expression on the inhibition of CYP3A-catalyzed drug metabolism: impact on modeling CYP3A-mediated drug-drug interactions

The purpose of this study was to determine the impact of CYP3A5 expression on inhibitory potency (Ki or IC50 values) of CYP3A inhibitors, using recombinant CYP3A4 and CYP3A5 (rCYP3A4 and rCYP3A5) and CYP3A5 genotyped human liver microsomes (HLMs). IC50 ratios between rCYP3A4 and rCYP3A5 (rCYP3A5/rCYP3A4) of ketoconazole (KTZ) and itraconazole (ITZ) were 8.5 and 8.8 for midazolam (MDZ), 4.7 and 9.1 for testosterone (TST), 1.3 and 2.8 for terfenadine, and 0.6 and 1.7 for vincristine, respectively, suggesting substrate- and inhibitor-dependent selectivity of the two azoles. Due to the difference in the IC50 values for CYP3A4 and CYP3A5, nonconcordant expression of CYP3A4 and CYP3A5 protein can significantly affect the observed magnitude of CYP3A-mediated drug-drug interactions in humans. Indeed, the IC50 values of KTZ and ITZ for CYP3A-catalyzed MDZ and TST metabolism were significantly higher in HLMs with CYP3A5*1/*1 and CYP3A5*1/*3 genotypes than in HLMs with the CYP3A5*3/*3 genotype, showing CYP3A5 expression-dependent IC50 values. Moreover, when IC50 values of KTZ and ITZ for different HLMs were kinetically simulated based on CYP3A5 expression level and enzyme-specific IC50 values, a good correlation between the simulated and the experimentally measured IC50 values was observed. Further simulation analysis revealed that both the Ki ratio (for inhibitors) and Vmax/Km ratio (for substrates) between CYP3A4 and CYP3A5 were major factors for CYP3A5 expression-dependent IC50 values. In conclusion, the present study demonstrates that CYP3A5 genotype and expression level have a significant impact on inhibitory potency for CYP3A-catalyzed drug metabolism, but that the magnitude of its effect is inhibitor-substrate pair specific.

Carboetomidate: an analog of etomidate that interacts weakly with 11β-hydroxylase

BACKGROUND

Carboetomidate is a pyrrole etomidate analog that is 3 orders of magnitude less potent an inhibitor of in vitro cortisol synthesis than etomidate (an imidazole) and does not inhibit in vivo steroid production. Although carboetomidate's reduced functional effect on steroid synthesis is thought to reflect lower binding affinity to 11β-hydroxylase, differential binding to this enzyme has never been experimentally demonstrated. In the current study, we tested the hypothesis that carboetomidate and etomidate bind with differential affinity to 11β-hydroxylase by comparing their abilities to inhibit photoaffinity labeling of purified enzyme by a photoactivatable etomidate analog and to modify the enzyme's absorption spectrum in a way that is indicative of ligand binding. In addition, we made a preliminary exploration of the manner in which etomidate and carboetomidate might differentially interact with this site using spectroscopic methods as well as molecular modeling techniques to better understand the structural basis for their selectivity.

METHODS

The ability of azi-etomidate to inhibit cortisol synthesis was tested by assessing its ability to inhibit cortisol synthesis by H295R cells. The binding affinities of etomidate and carboetomidate to 11β-hydroxylase were compared by assessing their abilities to (1) inhibit photoincorporation of the photolabile etomidate analog [(3)H]azi-etomidate into the enzyme and (2) modify the absorption spectrum of the enzyme's heme group. In silico docking studies of etomidate, carboetomidate, and azi-etomidate binding to 11β-hydroxylase were performed using the computer software GOLD.

RESULTS

Similar to etomidate, azi-etomidate potently inhibits in vitro cortisol synthesis. Etomidate inhibited [(3)H]azi-etomidate photolabeling of 11β-hydroxylase in a concentration-dependent manner. At a concentration of 40 µM, etomidate reduced photoincorporation of [(3)H]azi-etomidate by 96% ± 1% whereas carboetomidate had no experimentally detectable effect. On addition of etomidate to 11β-hydroxylase, a type 2 difference spectrum was produced indicative of etomidate complexation with the enzyme's heme iron; carboetomidate had no effect whereas azi-etomidate produced a reverse type 1 spectrum. Computer modeling studies predicted that etomidate, carboetomidate, and azi-etomidate can fit into the heme-containing pocket that forms 11β-hydroxylase's active site and pose with their carbonyl oxygens interacting with the heme iron and their phenyl rings stacking with phenylalanine-80. However, additional unique poses were identified for etomidate and azi-etomidate that likely account for their higher affinities.

CONCLUSIONS

Carboetomidate's reduced ability to suppress in vitro and in vivo steroid synthesis as compared with etomidate reflects its lower binding affinity to 11β-hydroxylase and may be attributed to carboetomidate's inability to form a coordination bond with the heme iron located at the enzyme's active site.

Measurement of electron transfer through cytochrome P450 protein on nanopillars and the effect of bound substrates

Electron transfer in cytochrome P450 enzymes is a fundamental process for activity. It is difficult to measure electron transfer in these enzymes because under the conditions typically used they exist in a variety of states. Using nanotechnology-based techniques, gold conducting nanopillars were constructed in an indexed array. The P450 enzyme CYP2C9 was attached to each of these nanopillars, and conductivity measurements made using conducting probe atomic force microscopy under constant force conditions. The conductivity measurements were made on CYP2C9 alone and with bound substrates, a bound substrate-effector pair, and a bound inhibitor. Fitting of the data with the Poole-Frenkel model indicates a correlation between the barrier height for electron transfer and the ease of CYP2C9-mediated metabolism of the bound substrates, though the spin state of iron is not well correlated. The approach described here should have broad application to the measurement of electron transfer in P450 enzymes and other metalloenzymes.

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.

1,2,3-Triazole-heme interactions in cytochrome P450: functionally competent triazole-water-heme complexes

In comparison to imidazole (IMZ) and 1,2,4-triazole (1,2,4-TRZ), the isosteric 1,2,3-triazole (1,2,3-TRZ) is unrepresented among cytochrome P450 (CYP) inhibitors. This is surprising because 1,2,3-TRZs are easily obtained via "click" chemistry. To understand this underrepresentation of 1,2,3-TRZs among CYP inhibitors, thermodynamic and density functional theory computational studies were performed with unsubstituted IMZ, 1,2,4-TRZ, and 1,2,3-TRZ. The results indicate that the lower affinity of 1,2,3-TRZ for the heme iron includes a large unfavorable entropy term likely originating in solvent-1,2,3-TRZ interactions; the difference is not solely due to differences in the enthalpy of heme-ligand interactions. In addition, the 1,2,3-TRZ fragment was incorporated into a well-established CYP3A4 substrate and mechanism-based inactivator, 17-α-ethynylestradiol (17EE), via click chemistry. This derivative, 17-click, yielded optical spectra consistent with low-spin ferric heme iron (type II) in contrast to 17EE, which yields a high-spin complex (type I). Furthermore, the rate of CYP3A4-mediated metabolism of 17-click was comparable to that of 17EE, with a different regioselectivity. Surprisingly, continuous-wave electron paramagnetic resonance (EPR) and HYSCORE EPR spectroscopy indicate that 17-click does not displace water from the sixth axial ligand position of CYP3A4 as expected for a type II ligand. We propose a binding model in which 17-click pendant 1,2,3-TRZ hydrogen bonds with the sixth axial water ligand. The results demonstrate the potential for 1,2,3-TRZ to form metabolically labile water-bridged low-spin heme complexes, consistent with recent evidence that nitrogenous type II ligands of CYPs can be efficiently metabolized. The specific case of [CYP3A4·17-click] highlights the risk of interpreting CYP-ligand complex structure on the basis of optical spectra.

Cytochrome P450 3A4 inhibition by ketoconazole: tackling the problem of ligand cooperativity using molecular dynamics simulations and free-energy calculations

Cytochrome P450 3A4 (CYP3A4) metabolizes more than 50% of clinically used drugs and is often involved in adverse drug-drug interactions. It displays atypical binding and kinetic behavior toward a number of ligands characterized by a sigmoidal shape of the corresponding titration curves, which is indicative of a positive homotropic cooperativity. This requires a participation of at least two ligand molecules, whereby the binding of the first ligand molecule increases the affinity of CYP3A4 for the binding of the second ligand molecule. In the current study, a combination of molecular dynamics simulations and free-energy calculations was applied to elucidate the physicochemical origin of the observed positive homotropic cooperativity in ketoconazole binding to CYP3A4. The binding of the first ketoconazole molecule was established to increase the affinity for the binding of the second ketoconazole molecule by 5 kJ mol(-1), which explains and quantifies the experimentally observed cooperative behavior of CYP3A4. Shape complementarity through nonpolar van der Waals interactions was identified as the main driving force of this binding, which seems to be in line with the promiscuous nature of CYP3A4. Moreover, the calculated binding free energies were found to be in good agreement with the values predicted from a simple 2-ligand binding kinetic model as well as to successfully reproduce the experimental titration curve. This confirms the general applicability of rapid free-energy methods to study challenging biomolecular systems like cytochromes P450, which are characterized by a large flexibility and malleability of their active sites.

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.

NMR-derived models of amidopyrine and its metabolites in complexes with rabbit cytochrome P450 2B4 reveal a structural mechanism of sequential N-dealkylation

To understand the molecular basis of sequential N-dealkylation by cytochrome P450 2B enzymes, we studied the binding of amidopyrine (AP) as well as the metabolites of this reaction, desmethylamidopyrine (DMAP) and aminoantipyrine (AAP), using the X-ray crystal structure of rabbit P450 2B4 and two nuclear magnetic resonance (NMR) techniques: saturation transfer difference (STD) spectroscopy and longitudinal (T(1)) relaxation NMR. Results of STD NMR of AP and its metabolites bound to P450 2B4 were similar, suggesting that they occupy similar niches within the enzyme's active site. The model-dependent relaxation rates (R(M)) determined from T(1) relaxation NMR of AP and DMAP suggest that the N-linked methyl is closest to the heme. To determine the orientation(s) of AP and its metabolites within the P450 2B4 active site, we used distances calculated from the relaxation rates to constrain the metabolites to the X-ray crystal structure of P450 2B4. Simulated annealing of the complex revealed that the metabolites do indeed occupy similar hydrophobic pockets within the active site, while the N-linked methyls are free to rotate between two binding modes. From these bound structures, a model of N-demethylation in which the N-linked methyl functional groups rotate between catalytic and noncatalytic positions was developed. This study is the first to provide a structural model of a drug and its metabolites complexed to a cytochrome P450 based on NMR and to provide a structural mechanism for how a drug can undergo sequential oxidations without unbinding. The rotation of the amide functional group might represent a common structural mechanism for N-dealkylation reactions for other drugs such as the local anesthetic lidocaine.

Interactions of cytochrome P450s with their ligands

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

Spectroscopic features of cytochrome P450 reaction intermediates

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