Substrate binding to cytochromes P450

Detection of a Kinetically Competent Compound-I Intermediate in the Vancomycin Biosynthetic Enzyme OxyB

Cytochrome P450 enzymes are abundantly encoded in microbial genomes. Their reactions have two general outcomes, one involving oxygen insertion via a canonical "oxygen rebound" mechanism and a second that diverts from this pathway and leads to a wide array of products, notably intramolecular oxidative cross-links. The antibiotic of-last-resort, vancomycin, contains three such cross-links, which are crucial for biological activity and are installed by the P450 enzymes OxyB, OxyA, and OxyC. The mechanisms of these enzymes have remained elusive in part because of the difficulty in spectroscopically capturing transient intermediates. Using stopped-flow UV/visible absorption and rapid freeze-quench electron paramagnetic resonance spectroscopies, we show that OxyB generates the highly reactive compound-I intermediate, which can react with a model vancomycin peptide substrate in a kinetically competent fashion to generate product. Our results have implications for the mechanism of OxyB and are in line with the notion that oxygen rebound and oxidative cross-links share early steps in their catalytic cycles.

Binding of steroid substrates reveals the key to the productive transition of the cytochrome P450 OleP

OleP is a bacterial cytochrome P450 involved in oleandomycin biosynthesis as it catalyzes regioselective epoxidation on macrolide intermediates. OleP has recently been reported to convert lithocholic acid (LCA) into murideoxycholic acid through a highly regioselective reaction and to unspecifically hydroxylate testosterone (TES). Since LCA and TES mainly differ by the substituent group at the C17, here we used X-ray crystallography, equilibrium binding assays, and molecular dynamics simulations to investigate the molecular basis of the diverse reactivity observed with the two steroids. We found that the differences in the structure of TES and LCA affect the capability of these molecules to directly form hydrogen bonds with N-terminal residues of OleP internal helix I. The establishment of these contacts, by promoting the bending of helix I, fosters an efficient trigger of the open-to-closed structural transition that occurs upon substrate binding to OleP and contributes to the selectivity of the subsequent monooxygenation reaction.

Engineering Electron Transfer Pathway of Cytochrome P450s

Cytochrome P450s (P450s), a superfamily of heme-containing enzymes, existed in animals, plants, and microorganisms. P450s can catalyze various regional and stereoselective oxidation reactions, which are widely used in natural product biosynthesis, drug metabolism, and biotechnology. In a typical catalytic cycle, P450s use redox proteins or domains to mediate electron transfer from NAD(P)H to heme iron. Therefore, the main factors determining the catalytic efficiency of P450s include not only the P450s themselves but also their redox-partners and electron transfer pathways. In this review, the electron transfer pathway engineering strategies of the P450s catalytic system are reviewed from four aspects: cofactor regeneration, selection of redox-partners, P450s and redox-partner engineering, and electrochemically or photochemically driven electron transfer.

Synthesis, in vitro and in silico anticancer evaluation of novel pyridin-2-yl estra-1,3,5(10)-triene derivatives

Aim: The aim of this study was the synthesis of steroid compounds with heterocyclic rings and good anticancer properties. Materials & methods: The synthesis, in silico and in vitro anticancer testing of novel pyridin-2-yl estra-1,3,5(10)-triene derivatives was performed. Results: All synthesized compounds have shown promising results for, antiproliferative activity, relative binding affinities for the ligand binding domains of estrogen receptors α, β and androgen receptor, aromatase binding potential, and inhibition of AKR1C3 enzyme. Conclusion: 3-Benzyloxy (17E)-pycolinilidene derivative 9 showed the best antitumor potential against MDA-MB-231 cell line, an activity that can be explained by its moderate inhibition of AKR1C3. Molecular docking simulation indicates that it binds to AKR1C3 in a very similar orientation and geometry as steroidal inhibitor EM1404.

Impact of Heterotropic Allosteric Modulation on the Time-Dependent Inhibition of Cytochrome P450 3A4

The current study was designed to investigate the influence of allosteric effectors on the metabolism of the prototypical cytochrome P450 (CYP) 3A4 substrate midazolam (MDZ), and on the determination in vitro time-dependent inhibition (TDI) of CYP3A4 using human liver microsomes (HLM). As the concentration of midazolam increased to 250 µM in HLMs, homotropic cooperativity resulted in a decrease in the 1'-hydroxymidazolam to 4-hydroxymidazolam ratio to a maximum of 1.1. The presence of varying concentrations of testosterone, progesterone (PGS), or carbamazepine (CBZ) in HLMs with MDZ could recapitulate the effect of homotropic cooperativity such that the formation rates of the 1'hydroxymidazolam and 4-hydroxymidazolam were equal even at low concentrations of MDZ. The presence of PGS (10 or 100 µM) and CBZ (100 or 1000 µM) in in vitro TDI determination of four known CYP3A4 time-dependent inactivators (clarithromycin, troleandomycin, mibefradil, raloxifene) simultaneously decreased potency and inactivation rate constant, resulting in fold changes in inactivation efficiency on average of 1.6-fold and 13-fold for the low and high concentrations of allosteric modulator tested, respectively. The formation of a metabolic-intermediate complex (MIC) for clarithromycin and troleandomycin decreased in the presence of the allosteric modulators in a concentration-dependent manner, reaching a new steady state formation that could not be overcome with increased incubation time. Maximum reduction of the MIC formed by clarithromycin was up to ∼91%, while troleandomycin MIC decreased up to ∼31%. These findings suggest that the absence of endogenous allosteric modulators may contribute to the poor translation of HLM-based drug-drug interaction predictions. SIGNIFICANCE STATEMENT: The reported overprediction of in vitro human liver microsome time-dependent inhibition of CYP3A4 and observed drug interactions in vivo remains an issue in drug development. We provide characterization of allosteric modulators on the CYP3A4 metabolism of the prototypical substrate midazolam, demonstrating the ability of the modulators to recapitulate the homotropic cooperativity of midazolam. Furthermore, we demonstrate that allosteric heterotropic cooperativity of CYP3A4 can impact the time-dependent inhibition kinetics of known mechanisms-based inhibitors, providing a potential mechanism to explain the overprediction.

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.

Pyridine-containing substrate analogs are restricted from accessing the human cytochrome P450 8B1 active site by tryptophan 281

The human oxysterol 12α-hydroxylase cytochrome P450 8B1 (CYP8B1) is a validated drug target for both type 2 diabetes and nonalcoholic fatty liver disease, but effective selective inhibitors are not yet available. Herein, steroidal substrate-mimicking compounds with a pyridine ring appended to the C12 site of metabolism were designed as inhibitors, synthesized, and evaluated in terms of their functional and structural interactions with CYP8B1. While the pyridine nitrogen was intended to coordinate the CYP8B1 active site heme iron, none of these compounds elicited shifts in the CYP8B1 Soret absorbance consistent with this type of interaction. However, when CYP8B1 was cocrystallized with the pyridine-containing compound with the 3-keto-Δ4 steroid backbone most similar to the endogenous substrate, it was apparent that this ligand was bound in a channel leading to the active site, instead of near the heme iron. Inspection of this structure suggested that tryptophan 281 directly above the heme might restrict active site binding of potential inhibitors with this design. This hypothesis was supported when a CYP8B1 W281F mutation did allow all three compounds to coordinate the heme iron as designed. These results indicated that the design of next-generation CYP8B1 inhibitors should be compatible with the low-ceiling tryptophan immediately above the heme iron.

Effect of Salts on the Conformational Dynamics of the Cytochrome P450 OleP

Cytochrome P450 OleP catalytic activity is strongly influenced by its structural dynamic conformational behavior. Here, we combine equilibrium-binding experiments with all-atom molecular dynamics simulations to clarify how different environments affect OleP conformational equilibrium between the open and the closed-catalytic competent-forms. Our data clearly show that at high-ionic strength conditions, the closed form is favored, and, very interestingly, different mechanisms, depending on the chemistry of the cations, can be used to rationalize such an effect.

Pseudomonas aeruginosa cytochrome P450 CYP168A1 is a fatty acid hydroxylase that metabolizes arachidonic acid to the vasodilator 19-HETE

Pseudomonas aeruginosa is a Gram-negative opportunistic human pathogen that is highly prevalent in individuals with cystic fibrosis (CF). A major problem in treating CF patients infected with P. aeruginosa is the development of antibiotic resistance. Therefore, the identification of novel P. aeruginosa antibiotic drug targets is of the utmost urgency. The genome of P. aeruginosa contains four putative cytochrome P450 enzymes (CYPs) of unknown function that have never before been characterized. Analogous to some of the CYPs from Mycobacterium tuberculosis, these P. aeruginosa CYPs may be important for growth and colonization of CF patients’ lungs. In this study, we cloned, expressed, and characterized CYP168A1 from P. aeruginosa and identified it as a subterminal fatty acid hydroxylase. Spectral binding data and computational modeling of substrates and inhibitors suggest that CYP168A1 has a large, expansive active site and preferentially binds long chain fatty acids and large hydrophobic inhibitors. Furthermore, metabolic experiments confirm that the enzyme is capable of hydroxylating arachidonic acid, an important inflammatory signaling molecule present in abundance in the CF lung, to 19-hydroxyeicosatetraenoic acid (19-HETE; K_m = 41 μM, V_max = 220 pmol/min/nmol P450), a potent vasodilator, which may play a role in the pathogen’s ability to colonize the lung. Additionally, we found that the in vitro metabolism of arachidonic acid is subject to substrate inhibition and is also inhibited by the presence of the antifungal agent ketoconazole. This study identifies a new metabolic pathway in this important human pathogen that may be of utility in treating P. aeruginosa infections.

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.

Multifunctional P450 Monooxygenase CftA Diversifies the Clifednamide Pool through Tandem C-H Bond Activations

Polycyclic tetramate macrolactams (PTMs) are a class of structurally complex hybrid polyketide-nonribosomal peptide (PK-NRP) natural products produced by diverse bacteria. Several PTMs display pharmaceutically interesting bioactivities, and the early stages of PTM biosynthesis involving polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) enzymology are well studied. However, the timing and mechanisms of post PKS-NRPS oxidations by P450 monooxygenases encoded in PTM biosynthetic gene clusters (BGCs) remain poorly characterized. Here we demonstrate that CftA, encoded in clifednamide-type PTM BGCs, is a multifunctional P450 monooxygenase capable of converting the C29-C30 ethyl side chain of ikarugamycin to either a C29-C30 methyl ketone or a C29-C30 hydroxymethyl ketone through C-H bond activation, resulting in the formation of clifednamide A or clifednamide C, respectively. We also report the complete structure of clifednamide C solved via multidimensional NMR (COSY, HSQC, HMBC, NOESY, and TOCSY) using material purified from an engineered Streptomyces strain optimized for production. Finally, the in vitro reconstitution of recombinant CftA catalytic activity revealed the oxidation cascade for sequential conversion of ikarugamycin to clifednamide A and clifednamide C. Our findings confirm prior genetics-based predictions on the origins of clifednamide complexity via P450s encoded in PTM BGCs and place CftA into a growing group of multifunctional P450s that tailor PTM natural products through late-stage regioselective C-H bond activation.

Point Mutations at a Key Site Alter the Cytochrome P450 OleP Structural Dynamics

Substrate binding to the cytochrome P450 OleP is coupled to a large open-to-closed transition that remodels the active site, minimizing its exposure to the external solvent. When the aglycone substrate binds, a small empty cavity is formed between the I and G helices, the BC loop, and the substrate itself, where solvent molecules accumulate mediating substrate-enzyme interactions. Herein, we analyzed the role of this cavity in substrate binding to OleP by producing three mutants (E89Y, G92W, and S240Y) to decrease its volume. The crystal structures of the OleP mutants in the closed state bound to the aglycone 6DEB showed that G92W and S240Y occupied the cavity, providing additional contact points with the substrate. Conversely, mutation E89Y induces a flipped-out conformation of this amino acid side chain, that points towards the bulk, increasing the empty volume. Equilibrium titrations and molecular dynamic simulations indicate that the presence of a bulky residue within the cavity impacts the binding properties of the enzyme, perturbing the conformational space explored by the complexes. Our data highlight the relevance of this region in OleP substrate binding and suggest that it represents a key substrate-protein contact site to consider in the perspective of redirecting its activity towards alternative compounds.

Molecular Cloning and Identification of NADPH Cytochrome P450 Reductase from Panax ginseng

Ginseng (Panax ginseng C.A. Mey.) is a precious Chinese traditional medicine, for which ginsenosides are the most important medicinal ingredients. Cytochrome P450 enzymes (CYP450) and their primary redox molecular companion NADPH cytochrome P450 reductase (CPR) play a key role in ginsenoside biosynthesis pathway. However, systematic studies of CPR genes in ginseng have not been reported. Numerous studies on ginsenoside synthesis biology still use Arabidopsis CPR (AtCPR1) as a reductase. In this study, we isolated two CPR genes (PgCPR1, PgCPR2) from ginseng adventitious roots. Phylogenetic tree analysis showed that both PgCPR1 and PgCPR2 are grouped in classⅡ of dicotyledonous CPR. Enzyme experiments showed that recombinant proteins PgCPR1, PgCPR2 and AtCPR1 can reduce cytochrome c and ferricyanide with NADPH as the electron donor, and PgCPR1 had the highest enzymatic activities. Quantitative real-time PCR analysis showed that PgCPR1 and PgCPR2 transcripts were detected in all examined tissues of Panax ginseng and both showed higher expression in stem and main root. Expression levels of the PgCPR1 and PgCPR2s were both induced after a methyl jasmonate (MeJA) treatment and its pattern matched with ginsenoside accumulation. The present investigation suggested PgCPR1 and PgCPR2 are associated with the biosynthesis of ginsenoside. This report will assist in future CPR family studies and ultimately improving ginsenoside production through transgenic engineering and synthetic biology.

19F-NMR reveals substrate specificity of CYP121A1 in Mycobacterium tuberculosis

Cytochromes P450 are versatile enzymes that function in endobiotic and xenobiotic metabolism and undergo meaningful structural changes that relate to their function. However, the way in which conformational changes inform the specific recognition of the substrate is often unknown. Here, we demonstrate the utility of fluorine (¹⁹F)-NMR spectroscopy to monitor structural changes in CYP121A1, an essential enzyme from Mycobacterium tuberculosis. CYP121A1 forms functional dimers that catalyze the phenol-coupling reaction of the dipeptide dicyclotyrosine. The thiol-reactive compound 3-bromo-1,1,1-trifluoroacetone was used to label an S171C mutation of the enzyme FG loop, which is located adjacent to the homodimer interface. Substrate titrations and inhibitor-bound ¹⁹F-NMR spectra indicate that ligand binding reduces conformational heterogeneity at the FG loop in both the dimer and in an engineered monomer of CYP121A1. However, only the dimer was found to promote a substrate-bound conformation that was preexisting in the substrate-free spectra, thus confirming a role for the dimer interface in dicyclotyrosine recognition. Moreover, ¹⁹F-NMR spectra in the presence of substrate analogs indicate the hydrogen-bonding feature of the dipeptide aromatic side chain as a dicyclotyrosine specificity criterion. This study demonstrates the utility of ¹⁹F-NMR as applied to a multimeric cytochrome P450, while also revealing mechanistic insights for an essential M. tuberculosis enzyme.

Biocatalytic synthesis of non-vicinal aliphatic diols

Biocatalysts are receiving increased attention in the field of selective oxyfunctionalization of C-H bonds, with cytochrome P450 monooxygenases (CYP450s), and the related peroxygenases, leading the field. Here we report on the substrate promiscuity of CYP505A30, previously characterized as a fatty acid hydroxylase. In addition to its regioselective oxyfunctionalization of saturated fatty acids (ω-1 - ω-3 hydroxylation), primary fatty alcohols are also accepted with similar regioselectivities. Moreover, alkanes such as n-octane and n-decane are also readily accepted, allowing for the production of non-vicinal diols through sequential oxygenation.

Dissecting the Cytochrome P450 OleP Substrate Specificity: Evidence for a Preferential Substrate

The cytochrome P450 OleP catalyzes the epoxidation of aliphatic carbons on both the aglycone 8.8a-deoxyoleandolide (DEO) and the monoglycosylated L-olivosyl-8.8a-deoxyoleandolide (L-O-DEO) intermediates of oleandomycin biosynthesis. We investigated the substrate versatility of the enzyme. X-ray and equilibrium binding data show that the aglycone DEO loosely fits the OleP active site, triggering the closure that prepares it for catalysis only on a minor population of enzyme. The open-to-closed state transition allows solvent molecules to accumulate in a cavity that forms upon closure, mediating protein–substrate interactions. In silico docking of the monoglycosylated L-O-DEO in the closed OleP–DEO structure shows that the L-olivosyl moiety can be hosted in the same cavity, replacing solvent molecules and directly contacting structural elements involved in the transition. X-ray structures of aglycone-bound OleP in the presence of L-rhamnose confirm the cavity as a potential site for sugar binding. All considered, we propose L-O-DEO as the optimal substrate of OleP, the L-olivosyl moiety possibly representing the molecular wedge that triggers a more efficient structural response upon substrate binding, favoring and stabilizing the enzyme closure before catalysis. OleP substrate versatility is supported by structural solvent molecules that compensate for the absence of a glycosyl unit when the aglycone is bound.

Misdiagnosis of CTX due to propofol: The interference of total intravenous propofol anaesthesia with bile acid profiling

BACKGROUND

Cerebrotendinous xanthomatosis (CTX) is a rare genetic disorder, characterised by chronic diarrhoea, xanthomas, cataracts, and neurological deterioration. CTX is caused by CYP27A1 deficiency, which leads to abnormal cholesterol and bile acid metabolism. Urinary bile acid profiling (increased m/z 627: glucuronide-5β-cholestane-pentol) serves as diagnostic screening for CTX. However, this led to a false positive CTX diagnosis in two patients, who had received total intravenous anaesthesia (TIVA) with propofol.

METHODS

To determine the influence of propofol on bile acid profiling, 10 urinary samples and 2 blood samples were collected after TIVA with propofol Fresenius 7 to 10 mg/kg/h from 12 subjects undergoing scoliosis correction. Urinary bile acids were analysed using flow injection negative electrospray mass spectrometry. Propofol binding to recombinant CYP27A1, the effects of propofol on recombinant CYP27A1 activity, and CYP27A1 expression in liver organoids were investigated using spectral binding, enzyme activity assays, and qPCR, respectively. Accurate masses were determined with high-resolution mass spectrometry.

RESULTS

Abnormal urinary profiles were identified in all subjects after TIVA, with a trend correlating propofol dose per kilogramme and m/z 627 peak intensity. Propofol only induced a weak CYP27A1 response in the spectral binding assay, minimally affected CYP27A1 activity and did not affect CYP27A1 expression. The accurate mass of m/z 627 induced by propofol differed >10 PPM from m/z 627 observed in CTX.

CONCLUSIONS

TIVA with propofol invariably led to a urinary profile misleadingly suggestive of CTX, but not through CYP27A1 inhibition. To avoid further misdiagnoses, propofol administration should be considered when interpreting urinary bile acid profiles.

Dual mechanisms suppress meloxicam bioactivation relative to sudoxicam

Thiazoles are biologically active aromatic heterocyclic rings occurring frequently in natural products and drugs. These molecules undergo typically harmless elimination; however, a hepatotoxic response can occur due to multistep bioactivation of the thiazole to generate a reactive thioamide. A basis for those differences in outcomes remains unknown. A textbook example is the high hepatotoxicity observed for sudoxicam in contrast to the relative safe use and marketability of meloxicam, which differs in structure from sudoxicam by the addition of a single methyl group. Both drugs undergo bioactivation, but meloxicam exhibits an additional detoxification pathway due to hydroxylation of the methyl group. We hypothesized that thiazole bioactivation efficiency is similar between sudoxicam and meloxicam due to the methyl group being a weak electron donator, and thus, the relevance of bioactivation depends on the competing detoxification pathway. For a rapid analysis, we modeled epoxidation of sudoxicam derivatives to investigate the impact of substituents on thiazole bioactivation. As expected, electron donating groups increased the likelihood for epoxidation with a minimal effect for the methyl group, but model predictions did not extrapolate well among all types of substituents. Through analytical methods, we measured steady-state kinetics for metabolic bioactivation of sudoxicam and meloxicam by human liver microsomes. Sudoxicam bioactivation was 6-fold more efficient than that for meloxicam, yet meloxicam showed a 6-fold higher efficiency of detoxification than bioactivation. Overall, sudoxicam bioactivation was 15-fold more likely than meloxicam considering all metabolic clearance pathways. Kinetic differences likely arise from different enzymes catalyzing respective metabolic pathways based on phenotyping studies. Rather than simply providing an alternative detoxification pathway, the meloxicam methyl group suppressed the bioactivation reaction. These findings indicate the impact of thiazole substituents on bioactivation is more complex than previously thought and likely contributes to the unpredictability of their toxic potential.

Comprehensive kinetic and modeling analyses revealed CYP2C9 and 3A4 determine terbinafine metabolic clearance and bioactivation

Terbinafine N-dealkylation pathways result in formation of 6,6-dimethyl-2-hepten-4-ynal (TBF-A), a reactive allylic aldehyde, that may initiate idiosyncratic drug-induced liver toxicity. Previously, we reported on the importance of CYP2C19 and 3A4 as major contributors to TBF-A formation. In this study, we expanded on those efforts to assess individual contributions of CYP1A2, 2B6, 2C8, 2C9, and 2D6 in terbinafine metabolism. The combined knowledge gained from these studies allowed us to scale the relative roles of the P450 isozymes in hepatic clearance of terbinafine including pathways leading to TBF-A, and hence, provide a foundation for assessing their significance in terbinafine-induced hepatotoxicity. We used in vitro terbinafine reactions with recombinant P450s to measure kinetics for multiple metabolic pathways and calculated contributions of all individual P450 isozymes to in vivo hepatic clearance for the average human adult. The findings confirmed that CYP3A4 was a major contributor (at least 30% total metabolism) to all three of the possible N-dealkylation pathways; however, CYP2C9, and not CYP2C19, played a critical role in terbinafine metabolism and even exceeded CYP3A4 contributions for terbinafine N-demethylation. A combination of their metabolic capacities accounted for at least 80% of the conversion of terbinafine to TBF-A, while CYP1A2, 2B6, 2C8, and 2D6 made minor contributions. Computational approaches provide a more rapid, less resource-intensive strategy for assessing metabolism, and thus, we additionally predicted terbinafine metabolism using deep neural network models for individual P450 isozymes. Cytochrome P450 isozyme models accurately predicted the likelihood for terbinafine N-demethylation, but overestimated the likelihood for a minor N-denaphthylation pathway. Moreover, the models were not able to differentiate the varying roles of the individual P450 isozymes for specific reactions with this particular drug. Taken together, the significance of CYP2C9 and 3A4 and to a lesser extent, CYP2C19, in terbinafine metabolism is consistent with reported drug interactions. This finding suggests that variations in individual P450 contributions due to other factors like polymorphisms may similarly contribute to terbinafine-related adverse health outcomes. Nevertheless, the impact of their metabolic capacities on formation of reactive TBF-A and consequent idiosyncratic hepatotoxicity will be mitigated by competing detoxification pathways, TBF-A decay, and TBF-A adduction to glutathione that remain understudied.

Heterologous coexpression of the benzoate-para-hydroxylase CYP53B1 with different cytochrome P450 reductases in various yeasts

Cytochrome P450 monooxygenases (P450) are enzymes with high potential as biocatalysts for industrial applications. Their large-scale applications are, however, limited by instability and requirement for coproteins and/or expensive cofactors. These problems are largely overcome when whole cells are used as biocatalysts. We previously screened various yeast species heterologously expressing self-sufficient P450s for their potential as whole-cell biocatalysts. Most P450s are, however, not self-sufficient and consist of two or three protein component systems. Therefore, in the present study, we screened different yeast species for coexpression of P450 and P450-reductase (CPR) partners, using CYP53B1 from Rhodotorula minuta as an exemplary P450. The abilities of three different coexpressed CPR partners to support P450 activity were investigated, two from basidiomycetous origin and one from an ascomycete. The various P450-CPR combinations were cloned into strains of Saccharomyces cerevisiae, Kluyveromyces marxianus, Hansenula polymorpha, Yarrowia lipolytica and Arxula adeninivorans, using a broad-range yeast expression vector. The results obtained supported the previous finding that recombinant A. adeninivorans strains perform excellently as whole-cell biocatalysts. This study also demonstrated for the first time the P450 reductase activity of the CPRs from R. minuta and U. maydis. A very interesting observation was the variation in the supportive activity provided by the different reductase partners tested and demonstrated better P450 activity enhancement by a heterologous CPR compared to its natural partner CPR. This study highlights reductase selection as a critical variable for consideration in the pursuit of optimal P450-based catalytic systems. The usefulness of A. adeninivorans as both a host for recombinant P450s and whole-cell biocatalyst was emphasized, supporting earlier findings.

Binding Modes and Metabolism of Caffeine

A correct estimate of ligand binding modes and a ratio of their occupancies is crucial for calculations of binding free energies. The newly developed method BLUES combines molecular dynamics with nonequilibrium candidate Monte Carlo. Nonequilibrium candidate Monte Carlo generates a plethora of possible binding modes and molecular dynamics enables the system to relax. We used BLUES to investigate binding modes of caffeine in the active site of its metabolizing enzyme Cytochrome P450 1A2 with the aim of elucidating metabolite-formation profiles at different concentrations. Because the activation energies of all sites of metabolism do not show a clear preference for one metabolite over the others, the orientations in the active site must play a key role. In simulations with caffeine located in a spacious pocket above the I-helix, it points N3 and N1 to the heme iron, whereas in simulations where caffeine is in close proximity to the heme N7 and C8 are preferably oriented toward the heme iron. We propose a mechanism where at low caffeine concentrations caffeine binds to the upper part of the active site, leading to formation of the main metabolite paraxanthine. On the other hand, at high concentrations two molecules are located in the active site, forcing one molecule into close proximity to the heme and yielding metabolites theophylline and trimethyluretic acid. Our results offer an explanation of previously published experimental results.

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.

Four Major Channels Detected in the Cytochrome P450 3A4: A Step toward Understanding Its Multispecificity

We computed the network of channels of the 3A4 isoform of the cytochrome P450 (CYP) on the basis of 16 crystal structures extracted from the Protein Data Bank (PDB). The calculations were performed with version 2 of the CCCPP software that we developed for this research project. We identified the minimal cost paths (MCPs) output by CCCPP as probable ways to access to the buried active site. The algorithm of calculation of the MCPs is presented in this paper, with its original method of visualization of the channels. We found that these MCPs constitute four major channels in CYP3A4. Among the many channels proposed by Cojocaru et al. in 2007, we found that only four of them open in 3A4. We provide a refined description of these channels together with associated quantitative data.

Rational and semi-rational engineering of cytochrome P450s for biotechnological applications

The cytochrome P450 enzymes are ubiquitous heme-thiolate proteins performing regioselective and stereoselective oxygenation reactions in cellular metabolism. Due to their broad substrate scope and catalytic versatility, P450 enzymes are also attractive candidates for many industrial and biopharmaceutical applications. For particular uses, enzyme properties of P450s can be further optimized through directed evolution, rational, and semi-rational engineering approaches, all of which introduce mutations within the P450 structures. In this review, we describe the recent applications of these P450 engineering approaches and highlight the key regions and residues that have been identified using such approaches. These “hotspots” lie within critical functional areas of the P450 structure, including the active site, the substrate access channel, and the redox partner interaction interface.

Substrate-induced conformational change in cytochrome P450 OleP

The regulation of cytochrome P450 activity is often achieved by structural transitions induced by substrate binding. We describe the conformational transition experienced upon binding by the P450 OleP, an epoxygenase involved in oleandomycin biosynthesis. OleP bound to the substrate analog 6DEB crystallized in 2 forms: one with an ensemble of open and closed conformations in the asymmetric unit and another with only the closed conformation. Characterization of OleP-6DEB binding kinetics, also using the P450 inhibitor clotrimazole, unveiled a complex binding mechanism that involves slow conformational rearrangement with the accumulation of a spectroscopically detectable intermediate where 6DEB is bound to open OleP. Data reported herein provide structural snapshots of key precatalytic steps in the OleP reaction and explain how structural rearrangements induced by substrate binding regulate activity.-Parisi, G., Montemiglio, L. C., Giuffrè, A., Macone, A., Scaglione, A., Cerutti, G., Exertier, C., Savino, C., Vallone, B. Substrate-induced conformational change in cytochrome P450 OleP.

Molecular Determinants of Substrate Affinity and Enzyme Activity of a Cytochrome P450BM3 Variant

Cytochrome P450_BM3 catalyzes the hydroxylation and/or epoxidation of fatty acids, fatty amides, and alcohols. Protein engineering has produced P450_BM3 variants capable of accepting drug molecules normally metabolized by human P450 enzymes. The enhanced substrate promiscuity has been attributed to the greater flexibility of the lid of the substrate channel. However, it is not well understood how structurally different and highly polar drug molecules can stably bind in the active site nor how the activity and coupling efficiency of the enzyme may be affected by the lack of enzyme-substrate complementarity. To address these important aspects of non-native small molecule binding, this study investigated the binding of drug molecules with different size, charge, polar surface area, and human P450 affinity on the promiscuous R47L/F87V/L188Q/E267V/F81I pentuple mutant of P450_BM3. Binding free energy data and energy decomposition analysis showed that pentuple mutant P450_BM3 stably binds (i.e., negative ΔG_b°) a broad range of substrate and inhibitor types because dispersion interactions with active site residues overcome unfavorable repulsive and electrostatic effects. Molecular dynamics simulations revealed that 1) acidic substrates tend to disrupt the heme propionate A-K69 salt bridge, which may reduce heme oxidizing ability, and 2) the lack of complementarity leads to high substrate mobility and water density in the active site, which may lead to uncoupling. These factors must be considered in future developments of P450_BM3 as a biocatalyst in the large-scale production of drug metabolites.

Catalytic modulation of human cytochromes P450 17A1 and P450 11B2 by phospholipid

Unlike most of the drug-metabolizing cytochrome P450s, microsomal P450 17A1 and mitochondrial P450 11B2 catalyze sequential multi-step reactions in steroid biosynthesis. The membrane phospholipid composition might be one parameter that modulates the efficiency and processivity of specific pathways. Here we systematically examined the effects of physiologically relevant phospholipids on the catalysis of purified P450 17A1, P450 11B2, and P450 11B1 in reconstituted assay systems. Both dioleoylphosphatidylcholine (DOPC, 18:1) and dilauroylphosphatidylcholine (DLPC, 12:0) were found to be very efficient in reconstituting 17-hydroxylase and 1720-lyase reactions of P450 17A1. Phosphatidylethanolamine (PE) specifically enhanced 1720-lyase activity up to 2.4-fold in the presence of phosphatidylcholine. On the other hand, P450 11B2-catalyzed production of aldosterone from 11-deoxycorticosterone was very low and from 18-hydroxycorticosterone nil, implying low processivity. DOPC or cardiolipin, which is exclusively located in the inner mitochondrial membrane, maximized aldosterone yield. In sharp contrast, reconstitution of homologous P450 11B1 with DOPC significantly decreased corticosterone formation without affecting the synthesis of 18-hydroxycorticosterone. The intrinsic fluorescence of P450 17A1 and 11B2 increased in the presence of DOPC, DLPC and PE. Acrylamide quenching studies showed that PE decreased solvent accessibility for tryptophan in P450 17A1, as did 20:4 PC or 18:2 PC for P450 11B2. A moderately positive correlation between the proportion of high-spin substrate-bound species and catalytic activity was only observed in the presence of phosphatidylcholines with low-temperature phase transition. These results demonstrate the potential for phospholipids to regulate the activity of steroidogenic P450 activities and thereby steroid hormone biosynthetic pathways.

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.

Pathways of Metabolite-Related Damage to a Synthetic p53 Gene Exon 7 Oligonucleotide Using Magnetic Enzyme Bioreactor Beads and LC-MS/MS Sequencing

Reactive metabolites of environmental chemicals and drugs can cause site specific damage to the p53 tumor suppressor gene in a major pathway for genotoxicity. We report here a high-throughput, cell-free, 96-well plate magnetic bead-enzyme system interfaced with LC-MS/MS sequencing for bioactivating test chemicals and identifying resulting adduction sites on genes. Bioactivated aflatoxin B1 was reacted with a 32 bp exon 7 fragment of the p53 gene using eight microsomal cytochrome (cyt) P450 enzymes from different organs coated on magnetic beads. All cyt P450s converted aflatoxin B1 to aflatoxin B1-8,9-epoxide that adducts guanine (G) in codon 249, with subsequent depurination to give abasic sites and then strand breaks. This is the first demonstration in a cell-free medium that the aflatoxin B1 metabolite selectively causes abasic site formation and strand breaks at codon 249 of the p53 probe, corresponding to the chemical pathway and mutations of p53 in human liver cells and tumors. Molecular modeling supports the view that binding of aflatoxin B1-8,9-epoxide to G in codon 249 precedes the S_N2 adduction reaction. Among a range of metabolic enzymes characteristic of different organs, human liver microsomes and cyt P450 3A5 supersomes showed the highest bioactivation rate for p53 exon 7 damage. This method of identifying metabolite-related gene damage sites may facilitate predictions of organ specific cancers for test chemicals via correlations with mutation sites.

[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.

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.

Drugs and Scaffold That Inhibit Cytochrome P450 27A1 In Vitro and In Vivo

Cytochrome P450 27A1 (CYP27A1) is a ubiquitous enzyme that hydroxylates cholesterol and other sterols. Complete CYP27A1 deficiency owing to genetic mutations is detrimental to human health, whereas 50% of activity retention is not and does not affect the whole body cholesterol levels. CYP27A1 is considered a potential therapeutic target in breast cancer and age-related neurodegenerative diseases; however, CYP27A1 inhibition should be ≤50%. Herein, 131 pharmaceuticals were tested for their effect on CYP27A1-mediated cholesterol 27-hydroxylation by in vitro enzyme assay. Of them, 14 drugs inhibited CYP27A1 by ≥75% and were evaluated for in vitro binding to the enzyme active site and for inhibition constants. All drugs except one (dasatinib) elicited a spectral response in CYP27A1 and had K_i values for cholesterol 27-hydroxylation either in the submicromolar (clevidipine, delavirdine, etravirine, felodipine, nicardipine, nilotinib, and sorafenib) or low micromolar range (abiratone, candesartan, celecoxib, dasatinib, nilvadipine, nimodipine, and regorafenib). Clevidipine, felodipine, nicardipine, nilvadipine, and nimodipine have the same 1,4-dihydropyridine scaffold and are indicated for hypertension. We used two of these antihypertensives (felodipine and nilvadipine) for administration to mice at a 1-mg/kg of body weight dose, daily, for 7 days. Mouse 27-hydroxycholesterol levels in the plasma, brain, and liver were reduced, whereas tissue levels of total cholesterol were unchanged. Structure-activity relationships within the 1,4-dihydropyridine scaffold were investigated, and features important for CY27A1 inhibition were identified. We confirmed our previous finding that CYP27A1 is a druggable enzyme and found additional drugs as well as the scaffold with potential for partial CYP27A1 inhibition in humans.

Polymerization of hexamethylene diisocyanate in solution and a 260.23 m/z [M+H]+ ion in exposed human cells

Hexamethylene diisocyanate (HDI) is an important industrial chemical that can cause asthma, however pathogenic mechanisms remain unclear. Upon entry into the respiratory tract, HDI's N=C=O groups may undergo nucleophilic addition (conjugate) to host molecules (e.g. proteins), or instead react with water (hydrolyze), releasing CO₂ and leaving a primary amine in place of the original N=C=O. We hypothesized that (primary amine groups present on) hydrolyzed or partially hydrolyzed HDI may compete with proteins and water as a reaction target for HDI in solution, resulting in polymers that could be identified and characterized using LC-MS and LC-MS/MS. Analysis of the reaction products formed when HDI was mixed with a pH buffered, isotonic, protein containing solution identified multiple [M+H]⁺ ions with m/z's and collision-induced dissociation (CID) fragmentation patterns consistent with those expected for dimers (259.25/285.23 m/z), and trimers (401.36/427.35 m/z) of partially hydrolyzed HDI (e.g. ureas/oligoureas). Human peripheral blood mononuclear cells (PBMCs) and monocyte-like U937, but not airway epithelial NCI-H292 cell lines cultured with these HDI ureas contained a novel 260.23 m/z [M+H]⁺ ion. LC-MS/MS analysis of the 260.23 m/z [M+H]⁺ ion suggest the formula C₁₃H₂₉N₃O₂ and a structure containing partially hydrolyzed HDI, however definitive characterization will require further orthogonal analyses.

Cross-linking of dicyclotyrosine by the cytochrome P450 enzyme CYP121 from Mycobacterium tuberculosis proceeds through a catalytic shunt pathway

CYP121, the cytochrome P450 enzyme in Mycobacterium tuberculosis that catalyzes a single intramolecular C-C cross-linking reaction in the biosynthesis of mycocyclosin, is crucial for the viability of this pathogen. This C-C coupling reaction represents an expansion of the activities carried out by P450 enzymes distinct from oxygen insertion. Although the traditional mechanism for P450 enzymes has been well studied, it is unclear whether CYP121 follows the general P450 mechanism or uses a different catalytic strategy for generating an iron-bound oxidant. To gain mechanistic insight into the CYP121-catalyzed reaction, we tested the peroxide shunt pathway by using rapid kinetic techniques to monitor the enzyme activity with its substrate dicyclotyrosine (cYY) and observed the formation of the cross-linked product mycocyclosin by LC-MS. In stopped-flow experiments, we observed that cYY binding to CYP121 proceeds in a two-step process, and EPR spectroscopy indicates that the binding induces active site reorganization and uniformity. Using rapid freeze-quenching EPR, we observed the formation of a high-spin intermediate upon the addition of peracetic acid to the enzyme-substrate complex. This intermediate exhibits a high-spin (S = 5/2) signal with g values of 2.00, 5.77, and 6.87. Likewise, iodosylbenzene could also produce mycocyclosin, implicating compound I as the initial oxidizing species. Moreover, we also demonstrated that CYP121 performs a standard peroxidase type of reaction by observing substrate-based radicals. On the basis of these results, we propose plausible free radical-based mechanisms for the C-C bond coupling reaction.

Cell foundry with high product specificity and catalytic activity for 21-deoxycortisol biotransformation

Background

21-deoxycortisol (21-DF) is the key intermediate to manufacture pharmaceutical glucocorticoids. Recently, a Japan patent has realized 21-DF production via biotransformation of 17-hydroxyprogesterone (17-OHP) by purified steroid 11β-hydroxylase CYP11B1. Due to the less costs on enzyme isolation, purification and stabilization as well as cofactors supply, whole-cell should be preferentially employed as the biocatalyst over purified enzymes. No reports as so far have demonstrated a whole-cell system to produce 21-DF. Therefore, this study aimed to establish a whole-cell biocatalyst to achieve 21-DF transformation with high catalytic activity and product specificity.

Results

In this study, Escherichia coli MG1655(DE3), which exhibited the highest substrate transportation rate among other tested chassises, was employed as the host cell to construct our biocatalyst by co-expressing heterologous CYP11B1 together with bovine adrenodoxin and adrenodoxin reductase. Through screening CYP11B1s (with mutagenesis at N-terminus) from nine sources, Homo sapiens CYP11B1 mutant (G25R/G46R/L52 M) achieved the highest 21-DF transformation rate at 10.6 mg/L/h. Furthermore, an optimal substrate concentration of 2.4 g/L and a corresponding transformation rate of 16.2 mg/L/h were obtained by screening substrate concentrations. To be noted, based on structural analysis of the enzyme-substrate complex, two types of site-directed mutations were designed to adjust the relative position between the catalytic active site heme and the substrate. Accordingly, 1.96-fold enhancement on 21-DF transformation rate (to 47.9 mg/L/h) and 2.78-fold improvement on product/by-product ratio (from 0.36 to 1.36) were achieved by the combined mutagenesis of F381A/L382S/I488L. Eventually, after 38-h biotransformation in shake-flask, the production of 21-DF reached to 1.42 g/L with a yield of 52.7%, which is the highest 21-DF production as known.

Conclusions

Heterologous CYP11B1 was manipulated to construct E. coli biocatalyst converting 17-OHP to 21-DF. Through the strategies in terms of (1) screening enzymes (with N-terminal mutagenesis) sources, (2) optimizing substrate concentration, and most importantly (3) rational design novel mutants aided by structural analysis, the 21-DF transformation rate was stepwise improved by 19.5-fold along with 4.67-fold increase on the product/byproduct ratio. Eventually, the highest 21-DF reported production was achieved in shake-flask after 38-h biotransformation. This study highlighted above described methods to obtain a high efficient and specific biocatalyst for the desired biotransformation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-017-0720-y) contains supplementary material, which is available to authorized users.

Ligand characterization of CYP4B1 isoforms modified for high-level expression in Escherichia coli and HepG2 cells

Human CYP4B1, a cytochrome P450 monooxygenase predominantly expressed in the lung, inefficiently metabolizes classical CYP4B1 substrates, such as the naturally occurring furan pro-toxin 4-ipomeanol (4-IPO). Highly active animal forms of the enzyme convert 4-IPO to reactive alkylating metabolite(s) that bind(s) to cellular macromolecules. By substitution of 13 amino acids, we restored the enzymatic activity of human CYP4B1 toward 4-IPO and this modified cDNA is potentially valuable as a suicide gene for adoptive T-cell therapies. In order to find novel pro-toxins, we tested numerous furan analogs in in vitro cell culture cytotoxicity assays by expressing the wild-type rabbit and variants of human CYP4B1 in human liver-derived HepG2 cells. To evaluate the CYP4B1 substrate specificities and furan analog catalysis, we optimized the N-terminal sequence of the CYP4B1 variants by modification/truncation and established their heterologous expression in Escherichia coli (yielding 70 and 800 nmol·l-1 of recombinant human and rabbit enzyme, respectively). Finally, spectral binding affinities and oxidative metabolism of the furan analogs by the purified recombinant CYP4B1 variants were analyzed: the naturally occurring perilla ketone was found to be the tightest binder to CYP4B1, but also the analog that was most extensively metabolized by oxidative processes to numerous non-reactive reaction products.

Nanodiscs in Membrane Biochemistry and Biophysics

Membrane proteins play a most important part in metabolism, signaling, cell motility, transport, development, and many other biochemical and biophysical processes which constitute fundamentals of life on the molecular level. Detailed understanding of these processes is necessary for the progress of life sciences and biomedical applications. Nanodiscs provide a new and powerful tool for a broad spectrum of biochemical and biophysical studies of membrane proteins and are commonly acknowledged as an optimal membrane mimetic system that provides control over size, composition, and specific functional modifications on the nanometer scale. In this review we attempted to combine a comprehensive list of various applications of nanodisc technology with systematic analysis of the most attractive features of this system and advantages provided by nanodiscs for structural and mechanistic studies of membrane proteins.

Conformational Mobility in Cytochrome P450 3A4 Explored by Pressure-Perturbation EPR Spectroscopy

We used high hydrostatic pressure as a tool for exploring the conformational landscape of human cytochrome P450 3A4 (CYP3A4) by electron paramagnetic resonance and fluorescence spectroscopy. Site-directed incorporation of a luminescence resonance energy transfer donor-acceptor pair allowed us to identify a pressure-dependent equilibrium between two states of the enzyme, where an increase in pressure increased the spatial separation between the two distantly located fluorophores. This transition is characterized by volume change (ΔV°) and P1/2 values of -36.8 ± 5.0 mL/mol and 1.45 ± 0.33 kbar, respectively, which corresponds to a Keq° of 0.13 ± 0.06, so that only 15% of the enzyme adopts the pressure-promoted conformation at ambient pressure. This pressure-promoted displacement of the equilibrium is eliminated by the addition of testosterone, an allosteric activator. Using site-directed spin labeling, we demonstrated that the pressure- and testosterone-sensitive transition is also revealed by pressure-induced changes in the electron paramagnetic resonance spectra of a nitroxide side chain placed at position 85 or 409 of the enzyme. Furthermore, we observed a pressure-induced displacement of the emission maxima of a solvatochromic fluorophore (7-diethylamino-3-((((2-maleimidyl)ethyl)amino)carbonyl) coumarin) placed at the same positions, which suggests a relocation to a more polar environment. Taken together, the results reveal an effector-dependent conformational equilibrium between open and closed states of CYP3A4 that involves a pronounced change at the interface between the region of α-helices A/A' and the meander loop of the enzyme, where residues 85 and 409 are located. Our study demonstrates the high potential of pressure-perturbation strategies for studying protein conformational landscapes.

Mechanism-Based Inhibition of CYP3A4 by Podophyllotoxin: Aging of an Intermediate Is Important for in Vitro/in Vivo Correlations

An in vitro observation of time-dependent inhibition (TDI) of metabolic enzymes often results in removing a potential drug from the drug pipeline. However, the accepted method for predicting TDIs of the important drug metabolizing cytochrome P450 enzymes often overestimates the drug interaction potential. Better models that take into account the complexities of the cytochrome P450 enzyme system will lead to better predictions. Herein we report the use of our previously described models for complex kinetics of podophyllotoxin. Spectral characterization of the kinetics indicates that an intermediate MI complex is formed, which slowly progresses to an essentially irreversible MI complex. The intermediate MI complex can release free enzyme during the time course of a typical 30 min TDI experiment. This slow rate of MI complex conversion results in an overprediction of the kinact value if this process is not included in the analysis of the activity versus time profile. In vitro kinetic experiments in rat liver microsomes predicted a lack of drug interaction between podophyllotoxin and midazolam. In vivo rat pharmacokinetic studies confirmed this lack of drug interaction.

Rational Development of Novel Activity Probes for the Analysis of Human Cytochromes P450

The identification and quantification of functional cytochromes P450 (CYPs) in biological samples is proving important for robust analyses of drug efficacy and metabolic disposition. In this study, a novel CYP activity-based probe was rationally designed and synthesised, demonstrating selective binding of CYP isoforms. The dependence of probe binding upon the presence of NADPH permits the selective detection of functionally active CYP. This allows the detection and analysis of these enzymes using biochemical and proteomic methodologies and approaches.

The Role of Protein-Protein and Protein-Membrane Interactions on P450 Function

This symposium summary, sponsored by the ASPET, was held at Experimental Biology 2015 on March 29, 2015, in Boston, Massachusetts. The symposium focused on: 1) the interactions of cytochrome P450s (P450s) with their redox partners; and 2) the role of the lipid membrane in their orientation and stabilization. Two presentations discussed the interactions of P450s with NADPH-P450 reductase (CPR) and cytochrome b5. First, solution nuclear magnetic resonance was used to compare the protein interactions that facilitated either the hydroxylase or lyase activities of CYP17A1. The lyase interaction was stimulated by the presence of b5 and 17α-hydroxypregnenolone, whereas the hydroxylase reaction was predominant in the absence of b5. The role of b5 was also shown in vivo by selective hepatic knockout of b5 from mice expressing CYP3A4 and CYP2D6; the lack of b5 caused a decrease in the clearance of several substrates. The role of the membrane on P450 orientation was examined using computational methods, showing that the proximal region of the P450 molecule faced the aqueous phase. The distal region, containing the substrate-access channel, was associated with the membrane. The interaction of NADPH-P450 reductase (CPR) with the membrane was also described, showing the ability of CPR to "helicopter" above the membrane. Finally, the endoplasmic reticulum (ER) was shown to be heterogeneous, having ordered membrane regions containing cholesterol and more disordered regions. Interestingly, two closely related P450s, CYP1A1 and CYP1A2, resided in different regions of the ER. The structural characteristics of their localization were examined. These studies emphasize the importance of P450 protein organization to their function.

Membrane Interactions, Ligand-Dependent Dynamics, and Stability of Cytochrome P4503A4 in Lipid Nanodiscs

Membrane-bound cytochrome P4503A4 (CYP3A4) is the major source of enzymatic drug metabolism. Although several structural models of CYP3A4 in various ligand complexes are available, none includes a lipid bilayer. Details of the effects of the membrane on protein dynamics and solvation, and access channels for ligands, remain uncertain. H/D exchange mass spectrometry (H/DXMS) with ligand free CYP3A4 containing a deletion of residues 3-12, compared to that of the full length wild type, in lipid nanodiscs afforded 91% sequence coverage. Deuterium exchange was fast in the F- and G-helices, HI loop, and C-terminal loop. In contrast, there is very low exchange in the F'- and G'-helices. The results are consistent with the overall membrane orientation of CYP3A4 suggested by published MD simulations and spectroscopic results, and the solvent accessibility of the F/G loop suggests that it is not deeply membrane-embedded. Addition of ketoconazole results in only modest, but global, changes in solvent accessibility. Interestingly, with ketoconazole bound some peptides become less solvent accessible or dynamic, including the F- and G-helices, but several peptides demonstrate modestly increased accessibility. Differential scanning calorimetry (DSC) of CYP3A4-nanodiscs suggests membrane-induced stabilization compared to that of aggregated CYP3A4 in buffer, and this stabilization is enhanced upon addition of the ligand ketoconazole. This ligand-induced stabilization is accompanied by a very large increase in ΔH for CYP3A4 denaturation in nanodiscs, possibly due to increased CYP3A4-membrane interactions. Together, the results suggest a distinct orientation of CYP3A4 on the lipid membrane, and they highlight likely solvent access channels, which are consistent with several MD simulations.

Functional analysis and crystallographic structure of clotrimazole bound OleP, a cytochrome P450 epoxidase from Streptomyces antibioticus involved in oleandomycin biosynthesis

BACKGROUND

OleP is a cyt P450 from Streptomyces antibioticus carrying out epoxigenation of the antibiotic oleandomycin during its biosynthesis. The timing of its reaction has not been fully clarified, doubts remain regarding its substrate and catalytic mechanism.

METHODS

The crystal structure of OleP in complex with clotrimazole, an inhibitor of P450s used in therapy, was solved and the complex formation dynamics was characterized by equilibrium and kinetic binding studies and compared to ketoconazole, another azole differing for the N1-substituent.

RESULTS

Clotrimazole coordinates the heme and occupies the active site. Most of the residues interacting with clotrimazole are conserved and involved in substrate binding in MycG, the P450 epoxigenase with the highest homology with OleP. Kinetic characterization of inhibitor binding revealed OleP to follow a simple bimolecular reaction, without detectable intermediates.

CONCLUSIONS

Clotrimazole-bound OleP adopts an open form, held by a π-π stacking chain that fastens helices F and G and the FG loop. Affinity is affected by the interactions of the N1 substituent within the active site, given the one order of magnitude difference of the off-rate constants between clotrimazole and ketoconazole. Based on structural similarities with MycG, we propose a binding mode for both oleandomycin intermediates, that are the candidate substrates of OleP.

GENERAL SIGNIFICANCE

Among P450 epoxigenases OleP is the only one that introduces an epoxide on a non-activated C–C bond. The data here presented are necessary to understand the rare chemistry carried out by OleP, to engineer it and to design more selective and potent P450-targeted drugs.