Anion-Dependent Stimulation of CYP3A4 Monooxygenase

Unraveling the Ligand-Binding Sites of CYP3A4 by Molecular Dynamics Simulations with Solvent Probes

Cytochrome P450 3A4 (CYP3A4) is one of the most important drug-metabolizing enzymes in the human body and is well known for its complicated, atypical kinetic characteristics. The existence of multiple ligand-binding sites in CYP3A4 has been widely recognized as being capable of interfering with the active pocket through allosteric effects. The identification of ligand-binding sites other than the canonical active site above the heme is especially important for understanding the atypical kinetic characteristics of CYP3A4 and the intriguing association between the ligand and the receptor. In this study, we first employed mixed-solvent molecular dynamics (MixMD) simulations coupled with the online computational predictive tools to explore potential ligand-binding sites in CYP3A4. The MixMD approach demonstrates better performance in dealing with the receptor flexibility compared with other computational tools. From the sites identified by MixMD, we then picked out multiple sites for further exploration using ensemble docking and conventional molecular dynamics (cMD) simulations. Our results indicate that three extra sites are suitable for ligand binding in CYP3A4, including one experimentally confirmed site and two novel sites.

Interaction of CYP3A4 with caffeine: First insights into multiple substrate binding

Human cytochrome P450 3A4 (CYP3A4) is a major drug-metabolizing enzyme that shows extreme substrate promiscuity. Moreover, its large and malleable active site can simultaneously accommodate several substrate molecules of the same or different nature, which may lead to cooperative binding and allosteric behavior. Due to difficulty of crystallization of CYP3A4-substrate complexes, it remains unknown how multiple substrates can arrange in the active site. We determined crystal structures of CYP3A4 bound to three and six molecules of caffeine, a psychoactive alkaloid serving as a substrate and modulator of CYP3A4. In the ternary complex, one caffeine binds to the active site suitably for C8-hydroxylation, most preferable for CYP3A4. In the senary complex, three caffeine molecules stack parallel to the heme with the proximal ligand poised for 3-N-demethylation. However, the caffeine stack forms extensive hydrophobic interactions that could preclude product dissociation and multiple turnovers. In both complexes, caffeine is also bound in the substrate channel and on the outer surface known as a peripheral site. At all sites, aromatic stacking with the caffeine ring(s) is likely a dominant interaction, while direct and water-mediated polar contacts provide additional stabilization for the substrate-bound complexes. Protein-ligand interactions via the active site R212, intrachannel T224, and peripheral F219 were experimentally confirmed, and the latter two residues were identified as important for caffeine association. Collectively, the structural, spectral, and mutagenesis data provide valuable insights on the ligand binding mechanism and help better understand how purine-based pharmaceuticals and other aromatic compounds could interact with CYP3A4 and mediate drug-drug interactions.

On the inhibition of cytochrome P450 3A4 by structurally diversified flavonoids

Cytochrome P450 3A4 (CYP3A4) is the most versatile enzyme involved in drug metabolism. The time-dependent inhibition of CYP3A4 by acacetin, apigenin, chrysin, and pinocembrin was experimentally detected, but not entirely elaborated so far. Thus, a two-level QM/MM (Quantum Mechanics/Molecular Mechanics) model is developed to yield insights into the receptor-flavonoid recognition at the molecular scale. Active site residues and the flavonoid are modelled using SCC-DFTB-D (QM level), while the rest of the complex is treated using AMBER force field (MM level). QM/MM binding free energies are well correlated with experimental data, indicating the largest inhibitory effect of chrysin on enzyme activity at a submicromolar concentration. Consequently, quercetin (QUE) and flavopiridol (FLP) are observed as representative examples of structurally different flavonoids. The inhibition parameters for QUE and FLP are evaluated using the well-calibrated QM/MM strategy, thereby aiding to quantitatively conceive the functional behavior of the whole family of flavonoids. A kinetic threshold for further assessment of the drug-drug interactions underlying the time-dependent inhibition of CYP3A4 by flavonoids is explored.Communicated by Ramaswamy H. Sarma.

Molecular Insights into the Heterotropic Allosteric Mechanism in Cytochrome P450 3A4-Mediated Midazolam Metabolism

Cytochrome P450 3A4 (CYP3A4) is the main P450 enzyme for drug metabolism and drug-drug interactions (DDIs), as it is involved in the metabolic process of approximately 50% of drugs. A detailed mechanistic elucidation of DDIs mediated by CYP3A4 is commonly believed to be critical for drug optimization and rational use. Here, two typical probes, midazolam (MDZ, substrate) and testosterone (TST, allosteric effector), are used to investigate the molecular mechanism of CYP3A4-mediated heterotropic allosteric interactions, through conventional molecular dynamics (cMD) and well-tempered metadynamics (WT-MTD) simulations. Distance monitoring shows that TST can stably bind in two potential peripheral sites (Site 1 and Site 2) of CYP3A4. The binding of TST at these two sites can induce conformational changes in CYP3A4 flexible loops on the basis of conformational analysis, thereby promoting the transition of the MDZ binding mode and affecting the ratio of MDZ metabolites. According to the results of the residue interaction network, multiple allosteric communication pathways are identified that can provide vivid and applicable insights into the heterotropic allostery of TST on MDZ metabolism. Comparing the regulatory effects and the communication pathways, the allosteric effect caused by TST binding in Site 2 seems to be more pronounced than in Site 1. Our findings could provide a deeper understanding of CYP3A4-mediated heterotropic allostery at the atomic level and would be helpful for rational drug use as well as the design of new allosteric modulators.

Atypical kinetics of cytochrome P450 enzymes in pharmacology and toxicology

Atypical kinetics are observed in metabolic reactions catalyzed by cytochrome P450 enzymes (P450). Yet, this phenomenon is regarded as experimental artifacts in some instances despite increasing evidence challenging the assumptions of typical Michaelis-Menten kinetics. As P450 play a major role in the metabolism of a wide range of substrates including drugs and endogenous compounds, it becomes critical to consider the impact of atypical kinetics on the accuracy of estimated kinetic and inhibitory parameters which could affect extrapolation of pharmacological and toxicological implications. The first half of this book chapter will focus on atypical non-Michaelis-Menten kinetics (e.g. substrate inhibition, biphasic and sigmoidal kinetics) as well as proposed underlying mechanisms supported by recent insights in mechanistic enzymology. In particular, substrate inhibition kinetics in P450 as well as concurrent drug inhibition of P450 in the presence of substrate inhibition will be further discussed. Moreover, mounting evidence has revealed that despite the high degree of sequence homology between CYP3A isoforms (i.e. CYP3A4 and CYP3A5), they have the propensities to exhibit vastly different susceptibilities and potencies of mechanism-based inactivation (MBI) with a common drug inhibitor. These experimental observations pertaining to the presence of these atypical isoform- and probe substrate-specific complexities in CYP3A isoforms by several clinically-relevant drugs will therefore be expounded and elaborated upon in the second half of this book chapter.

Transgenic Expression of Haemonchus contortus Cytochrome P450 Hco-cyp-13A11 Decreases Susceptibility to Particular but Not All Macrocyclic Lactones in the Model Organism Caenorhabditis elegans

The number of reported macrocyclic lactones (ML) resistance cases across all livestock hosts is steadily increasing. Different studies in the parasitic nematode Haemonchus contortus assume the participation of cytochrome P450s (Cyps) enzymes in ML resistance. Still, functional data about their individual contribution to resistance or substrate specificity is missing. Via microinjection, transgenic Caenorhabditis elegans expressing HCON_00141052 (transgene-Hco-cyp-13A11) from extrachromosomal arrays were generated. After 24 h of exposure to different concentrations of ivermectin (IVM), ivermectin aglycone (IVMa), selamectin (SEL), doramectin (DRM), eprinomectin (EPR), and moxidectin (MOX), motility assays were performed to determine the impact of the H. contortus Cyp to the susceptibility of the worms against each ML. While transgene-Hco-cyp-13A11 significantly decreased susceptibility to IVM (four-fold), IVMa (2-fold), and SEL (3-fold), a slight effect for DRM and no effect for MOX, and EPR was observed. This substrate specificity of Hco-cyp-13A11 could not be explained by molecular modeling and docking studies. Hco-Cyp-13A11 molecular models were obtained for alleles from isolates with different resistance statuses. Although 14 amino acid polymorphisms were detected, none was resistance specific. In conclusion, Hco-cyp-13A11 decreased IVM, IVMa, and SEL susceptibility to a different extent, but its potential impact on ML resistance is not driven by polymorphisms.

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.

Midazolam as a Probe for Heterotropic Drug-Drug Interactions Mediated by CYP3A4

Human cytochrome P450 CYP3A4 is involved in the processing of more than 35% of current pharmaceuticals and therefore is responsible for multiple drug-drug interactions (DDI). In order to develop a method for the detection and prediction of the possible involvement of new drug candidates in CYP3A4-mediated DDI, we evaluated the application of midazolam (MDZ) as a probe substrate. MDZ is hydroxylated by CYP3A4 in two positions: 1-hydroxy MDZ formed at lower substrate concentrations, and up to 35% of 4-hydroxy MDZ at high concentrations. The ratio of the formation rates of these two products (the site of metabolism ratio, SOM) was used as a measure of allosteric heterotropic interactions caused by effector molecules using CYP3A4 incorporated in lipid nanodiscs. The extent of the changes in the SOM in the presence of effectors is determined by chemical structure and is concentration-dependent. MD simulations of CYP3A4 in the lipid bilayer suggest that experimental results can be explained by the movement of the F-F' loop and concomitant changes in the shape and volume of the substrate-binding pocket. As a result of PGS binding at the allosteric site, several residues directly contacting MDZ move away from the substrate molecule, enabling the repositioning of the latter for minor product formation.

A Comprehensive Analysis of Human CYP3A4 Crystal Structures as a Potential Tool for Molecular Docking-Based Site of Metabolism and Enzyme Inhibition Studies

The notable ability of human liver cytochrome P450 3A4 (CYP3A4) to metabolize diverse xenobiotics encourages researchers to explore in-depth the mechanism of enzyme action. Numerous CYP3A4 protein crystal structures have been deposited in protein data bank (PDB) and are majorly used in molecular docking analysis. The quality of the molecular docking results depends on the three-dimensional CYP3A4 protein crystal structures from the PDB. Present review endeavors to provide a brief outline of some technical parameters of CYP3A4 PDB entries as valuable information for molecular docking research. PDB entries between 22 April 2004 and 2 June 2021 were compiled and the active sites were thoroughly observed. The present review identified 76 deposited PDB entries and described basic information that includes CYP3A4 from human genetic, Escherichia coli (E. coli) use for protein expression, crystal structure obtained from X-ray diffraction method, taxonomy ID 9606, Uniprot ID P08684, ligand–protein structure description, co-crystal ligand, protein site deposit and resolution ranges between 1.7[Formula: see text]Å and 2.95[Formula: see text]Å. The observation of protein–ligand interactions showed the various residues on the active site depending on the ligand. The residues Ala305, Ser119, Ala370, Phe304, Phe108, Phe213 and Phe215 have been found to frequently interact with ligands from CYP3A4 PDB. Literature surveys of 17 co-crystal ligands reveal multiple mechanisms that include competitive inhibition, noncompetitive inhibition, mixed-mode inhibition, mechanism-based inhibition, substrate with metabolite, inducer, or combination modes of action. This overview may help researchers choose a trustworthy CYP3A4 protein structure from the PDB database to apply the protein in molecular docking analysis for drug discovery.

Stepwise binding of inhibitors to human cytochrome P450 17A1 and rapid kinetics of inhibition of androgen biosynthesis

Cytochrome P450 (P450) 17A1 catalyzes the 17α-hydroxylation of progesterone and pregnenolone as well as the subsequent lyase cleavage of both products to generate androgens. However, the selective inhibition of the lyase reactions, particularly with 17α-hydroxy pregnenolone, remains a challenge for the treatment of prostate cancer. Here, we considered the mechanisms of inhibition of drugs that have been developed to inhibit P450 17A1, including ketoconazole, seviteronel, orteronel, and abiraterone, the only approved inhibitor used for prostate cancer therapy, as well as clotrimazole, known to inhibit P450 17A1. All five compounds bound to P450 17A1 in a multistep process, as observed spectrally, over a period of 10 to 30 s. However, no lags were observed for the onset of inhibition in rapid-quench experiments with any of these five compounds. Furthermore, the addition of substrate to inhibitor–P450 17A1 complexes led to an immediate formation of product, without a lag that could be attributed to conformational changes. Although abiraterone has been previously described as showing slow-onset inhibition (t_1/2 = 30 min), we observed rapid and strong inhibition. These results are in contrast to inhibitors of P450 3A4, an enzyme with a larger active site in which complete inhibition is not observed with ketoconazole and clotrimazole until the changes are completed. Overall, our results indicate that both P450 17A1 reactions—17α-hydroxylation and lyase activity—are inhibited by the initial binding of any of these inhibitors, even though subsequent conformational changes occur.

Structural Dynamics of Cytochrome P450 3A4 in the Presence of Substrates and Cytochrome P450 Reductase

Cytochrome P450 3A4 (CYP3A4) is the most important drug-metabolizing enzyme in humans and has been associated with harmful drug interactions. The activity of CYP3A4 is known to be modulated by several compounds and by the electron transfer partner, cytochrome P450 reductase (CPR). The underlying mechanism of these effects, however, is poorly understood. We have used hydrogen-deuterium exchange mass spectrometry to investigate the impact of binding of CPR and of three different substrates (7-benzyloxy-4-trifluoromethyl-coumarin, testosterone, and progesterone) on the conformational dynamics of CYP3A4. Here, we report that interaction of CYP3A4 with substrates or with the oxidized or reduced forms of CPR leads to a global rigidification of the CYP3A4 structure. This was evident from the suppression of deuterium exchange in several regions of CYP3A4, including regions known to be involved in protein-protein interactions (helix C) and substrate binding and specificity (helices B' and E, and loop K/β1). Furthermore, the bimodal isotopic distributions observed for some CYP3A4-derived peptides were drastically impacted upon binding to CPR and/or substrates, suggesting the existence of stable CYP3A4 conformational populations that are perturbed by ligand/CPR binding. The results have implications for understanding the mechanisms of ligand binding, allostery, and catalysis in CYP enzymes.

Combining Small-Molecule Bioconjugation and Hydrogen-Deuterium Exchange Mass Spectrometry (HDX-MS) to Expose Allostery: the Case of Human Cytochrome P450 3A4

We report a novel approach to study allostery which combines the use of carefully selected bioconjugates and hydrogen-deuterium exchange mass spectrometry (HDX-MS). This strategy avoids issues related to weak substrate binding and ligand relocalization. The utility of our method is demonstrated using human cytochrome P450 3A4 (CYP3A4), the most important drug-metabolizing enzyme. Allosteric activation and inhibition of CYP3A4 by pharmaceuticals is an important mechanism of drug interactions. We performed HDX-MS analysis on several CYP3A4-effector bioconjugates, some of which mimic the allosteric effect of positive effectors, while others show activity enhancement even though the label does not occupy the allosteric pocket (agonistic) or do not show activation while still blocking the allosteric site (antagonistic). This allowed us to better define the position of the allosteric site, the protein structural dynamics associated with allosteric activation, and the presence of coexisting conformers.

Midazolam as a Probe for Drug-Drug Interactions Mediated by CYP3A4: Homotropic Allosteric Mechanism of Site-Specific Hydroxylation

We developed an efficient and sensitive probe for drug-drug interactions mediated by human CYP3A4 by using midazolam (MDZ) as a probe substrate. Using global analysis of four parameters over several experimental data sets, we demonstrate that the first MDZ molecule (MDZ1) binds with high affinity at the productive site near the heme iron and gives only hydroxylation at the 1 position (1OH). The second midazolam molecule (MDZ2) binds at an allosteric site at the membrane surface and perturbs the position and mobility of MDZ1 such that the minor hydroxylation product at the 4 position (4OH) is formed in a 1:2 ratio (35%). No increase in catalytic rate is observed after the second MDZ binding. Hence, the site of the 1OH:4OH metabolism ratio is a sensitive probe for drugs, such as progesterone, that bind with high affinity to the allosteric site and serve as effectors. We observe similar changes in the MDZ 1OH:4OH ratio in the presence of progesterone (PGS), suggesting a direct communication between the active and allosteric sites. Mutations introduced into the F-F' loop indicate that residues F213 and D214 are directly involved in allosteric interactions leading to MDZ homotropic cooperativity, and these same residues, together with L211, are involved in heterotropic allosteric interactions in which PGS is the effector and MDZ the substrate. Molecular dynamics simulations provide a mechanistic picture of the origin of this cooperativity. These results show that the midazolam can be used as a sensitive probe for drug-drug interactions in human P450 CYP3A4.

Homotropic Cooperativity of Midazolam Metabolism by Cytochrome P450 3A4: Insight from Computational Studies

Human cytochrome P450 3A4 (CYP3A4) is responsible for the metabolism of ∼50% clinically used drugs. Midazolam (MDZ) is a commonly used sedative drug and serves as a marker substrate for the CYP3A4 activity assessment. MDZ is metabolized by CYP3A4 to two hydroxylation products, 1′-OH-MDZ and 4-OH-MDZ. It has been reported that the ratio of 1′-OH-MDZ and 4-OH-MDZ is dependent on the MDZ concentration, which reflects the homotropic cooperative behavior in MDZ metabolism by CYP3A4. Here, we used quantum chemistry (QC), molecular docking, conventional molecular dynamics (cMD), and Gaussian accelerated molecular dynamics (GaMD) approaches to investigate the mechanism of the interactions between CYP3A4 and MDZ. QC calculations suggest that C1′ is less reactive for hydroxylation than C4, which is a pro-chirality carbon. However, the 4-OH-MDZ product is likely to be racemic due to the chirality inversion in the rebound step. The MD simulation results indicate that MDZ at the peripheral allosteric site is not stable and the binding modes of the MDZ molecules at the productive site are in line with the experimental observations.

CYP2J2 Molecular Recognition: A New Axis for Therapeutic Design

Cytochrome P450 (CYP) epoxygenases are a special subset of heme-containing CYP enzymes capable of performing the epoxidation of polyunsaturated fatty acids (PUFA) and the metabolism of xenobiotics. This dual functionality positions epoxygenases along a metabolic crossroad. Therefore, structure-function studies are critical for understanding their role in bioactive oxy-lipid synthesis, drug-PUFA interactions, and for designing therapeutics that directly target the epoxygenases. To better exploit CYP epoxygenases as therapeutic targets, there is a need for improved understanding of epoxygenase structure-function. Of the characterized epoxygenases, human CYP2J2 stands out as a potential target because of its role in cardiovascular physiology. In this review, the early research on the discovery and activity of epoxygenases is contextualized to more recent advances in CYP epoxygenase enzymology with respect to PUFA and drug metabolism. Additionally, this review employs CYP2J2 epoxygenase as a model system to highlight both the seminal works and recent advances in epoxygenase enzymology. Herein we cover CYP2J2's interactions with PUFAs and xenobiotics, its tissue-specific physiological roles in diseased states, and its structural features that enable epoxygenase function. Additionally, the enumeration of research on CYP2J2 identifies the future needs for the molecular characterization of CYP2J2 to enable a new axis of therapeutic design.

In vitro evaluation of the metabolic enzymes and drug interaction potential of triapine

PURPOSE

To investigate the metabolic pathways of triapine in primary cultures of human hepatocytes and human hepatic subcellular fractions; to investigate interactions of triapine with tenofovir and emtricitabine; and to evaluate triapine as a perpetrator of drug interactions. The results will better inform future clinical studies of triapine, a radiation sensitizer currently being studied in a phase III study.

METHODS

Triapine was incubated with human hepatocytes and subcellular fractions in the presence of a number of inhibitors of drug metabolizing enzymes. Triapine depletion was monitored by LC-MS/MS. Tenofovir and emtricitabine were co-incubated with triapine in primary cultures of human hepatocytes. Triapine was incubated with a CYP probe cocktail and human liver microsomes, followed by LC-MS/MS monitoring of CYP specific metabolite formation.

RESULTS

Triapine was not metabolized by FMO, AO/XO, MAO-A/B, or NAT-1/2, but was metabolized by CYP450s. CYP1A2 accounted for most of the depletion of triapine. Tenofovir and emtricitabine did not alter triapine depletion. Triapine reduced CYP1A2 activity and increased CYP2C19 activity.

CONCLUSION

CYP1A2 metabolism is the major metabolic pathway for triapine. Triapine may be evaluated in cancer patients in the setting of HIV with emtricitabine or tenofovir treatment. Confirmatory clinical trials may further define the in vivo triapine metabolic fate and quantify any drug-drug interactions.

The crucial role of metabolic regulation in differential hepatotoxicity induced by furanoids in Dioscorea bulbifera

Diterpenoid lactones (DLs), a group of furan-containing compounds found in Dioscorea bulbifera L. (DB), have been reported to be associated with hepatotoxicity. Different hepatotoxicities of these DLs have been observed in vitro, but reasonable explanations for the differential hepatotoxicity have not been provided. Herein, the present study aimed to confirm the potential factors that contribute to varied hepatotoxicity of four representative DLs (diosbulbins A, B, C, F). In vitro toxic effects were evaluated in various cell models and the interactions between DLs and CYP3A4 at the atomic level were simulated by molecular docking. Results showed that DLs exhibited varied cytotoxicities, and that CYP3A4 played a modulatory role in this process. Moreover, structural variation may cause different affinities between DLs and CYP3A4, which was positively correlated with the observation of cytotoxicity. In addition, analysis of the glutathione (GSH) conjugates indicated that reactive intermediates were formed by metabolic oxidation that occurred on the furan moiety of DLs, whereas, GSH consumption analysis reflected the consistency between the reactive metabolites and the hepatotoxicity. Collectively, our findings illustrated that the metabolic regulation played a crucial role in generating the varied hepatotoxicity of DLs.

Conformational selection is present in ligand binding to cytochrome P450 19A1 lipoprotein nanodiscs

Cytochromes P450 (CYPs) display remarkable plasticity in their ability to bind substrates and catalyze a broad array of chemical reactions. Herein we evaluate binding of androstenedione, testosterone, and 7-hydroxyflavone to CYP19A1, also known as aromatase, in phospholipid nanodiscs by stopped-flow UV-vis spectroscopy. Exponential fitting of the kinetic traces supports the possibility of a multi-step binding mechanism. Subsequent global fitting of the data to the solutions of the coupled differential equations describing the fundamental mechanisms of induced fit and conformational selection, consistently support presence of the latter. To our knowledge, this is the first discrimination of conformational selection from induced fit for a mono-disperse CYP in a native-like membrane environment. In addition, 7-hydroxyflavone binds to CYP19A1 nanodiscs with comparable affinity to the substrates and induces an unusual spectral response likely attributable to hydrogen bonding to, rather than displacement of the heme-coordinated water molecule.

Functional impact of cytochrome P450 3A (CYP3A) missense variants in cattle

Cytochrome P450 3A is the most important CYP subfamily in humans, and CYP3A4/CYP3A5 genetic variants contribute to inter-individual variability in drug metabolism. However, no information is available for bovine CYP3A (bCYP3A). Here we described bCYP3A missense single nucleotide variants (SNVs) and evaluated their functional effects. CYP3A28, CYP3A38 and CYP3A48 missense SNVs were identified in 300 bulls of Piedmontese breed through targeted sequencing. Wild-type and mutant bCYP3A cDNAs were cloned and expressed in V79 cells. CYP3A-dependent oxidative metabolism of testosterone (TST) and nifedipine (NIF) was assessed by LC-MS/MS. Finally, SNVs functional impact on TST hydroxylation was measured ex vivo in liver microsomes from individually genotyped animals. Thirteen missense SNVs were identified and validated. Five variants showed differences in CYP3A catalytic activity: three CYP3A28 SNVs reduced TST 6β-hydroxylation; one CYP3A38 variant increased TST 16β-hydroxylation, while a CYP3A48 SNV showed enhanced NIF oxidation. Individuals homozygous for rs384467435 SNV showed a reduced TST 6β-hydroxylation. Molecular modelling showed that most of SNVs were distal to CYP3A active site, suggesting indirect effects on the catalytic activity. Collectively, these findings demonstrate the importance of pharmacogenetics studies in veterinary species and suggest bCYP3A genotype variation might affect the fate of xenobiotics in food-producing species such as cattle.

Structural Perspectives of the CYP3A Family and Their Small Molecule Modulators in Drug Metabolism

Cytochrome P450 enzymes function to catalyze a wide range of reactions, many of which are critically important for drug response. Members of the human cytochrome P450 3A (CYP3A) family are particularly important in drug clearance, and they collectively metabolize more than half of all currently prescribed medications. The ability of these enzymes to bind a large and structurally diverse set of compounds increases the chances of their modulating or facilitating drug metabolism in unfavorable ways. Emerging evidence suggests that individual enzymes in the CYP3A family play discrete and important roles in catalysis and disease progression. Here we review the similarities and differences among CYP3A enzymes with regard to substrate recognition, metabolism, modulation by small molecules, and biological consequence, highlighting some of those with clinical significance. We also present structural perspectives to further characterize the basis of these comparisons.

Conformational Response of N-Terminally Truncated Cytochrome P450 3A4 to Ligand Binding in Solution

Human cytochrome P450 3A4 (CYP3A4) is a membrane-associated monooxygenase that is responsible for metabolizing >50% of the pharmaceuticals in the current market, so studying its chemical mechanism and structural changes upon ligand binding will help provide deeper insights into drug metabolism and further drug development. The best-characterized cytochrome P450 is a bacterial form, P450cam, which undergoes significant conformational changes upon binding substrate and its redox partner, putidaredoxin. In contrast, most crystal structures of CYP3A4 with or without ligands have shown few changes, although allosteric effects and multiple-substrate binding in solution are well-documented. In this study, we use double electron-electron resonance (DEER) to measure distances between spatially separated spin-labels on CYP3A4 and molecular dynamics to interpret the DEER data. These methods were applied to a soluble N-terminally truncated CYP3A4 form, and the results show that there are few changes in the average structure upon binding ketoconazole, ritonavir, or midazolam. However, binding of midazolam, but not ketoconazole or ritonavir, resulted in a significant change in the motion and/or disorder in the F/G helix region near the substrate binding pocket. These results suggest that soluble CYP3A4 behaves in a unique way in response to inhibitor and substrate binding.

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.

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.

Steroid bioconjugation to a CYP3A4 allosteric site and its effect on substrate binding and coupling efficiency

Human cytochrome P450 3A4 (CYP3A4) is an important drug metabolizing enzyme involved in a number of drug-drug and food-drug interactions. As such, much effort has been devoted into investigating its mechanism of interaction with ligands. CYP3A4 has one of the highest levels of substrate promiscuity for an enzyme, and can even bind multiple ligands simultaneously. The location and orientation of these ligands depend on the chemical structure and stoichiometry, and are generally poorly understood. In the case of the steroid testosterone, up to three copies of the molecule can associate with the enzyme at once, likely two in the active site and one at a postulated allosteric site. Recently, we demonstrated that steroid bioconjugation at the allosteric site results in an increase in activity of CYP3A4 toward testosterone and 7-benzyloxy-4-trifluoromethylcoumarin oxidation. Here, using the established bioconjugation methodology, we show how steroid bioconjugation at the allosteric site affects the heme spin state, the binding affinity (K_S) of CYP3A4 for testosterone, as well as the enzyme coupling efficiency.

Membrane-embedded substrate recognition by cytochrome P450 3A4

Cytochrome P450 3A4 (CYP3A4) is the dominant xenobiotic-metabolizing enzyme in the liver and intestine and is involved in the disposition of more than 50% of drugs. Because of its ability to bind multiple substrates, its reaction kinetics are complex, and its association with the microsomal membrane confounds our understanding of how this enzyme recognizes and recruits diverse substrates. Testosterone (TST) hydroxylation is the prototypical CYP3A4 reaction, displaying positive homotropic cooperativity with three binding sites. Here, exploiting the capability of accelerated molecular dynamics (aMD) to sample events in the millisecond regime, I performed >25-μs aMD simulations in the presence of three TST molecules. These simulations identified high-occupancy surface-binding sites as well as a pathway for TST ingress into the CYP3A4 active site originating in the membrane. Adaptive biasing force analysis of the latter pathway revealed a metastable intermediate that could constitute a third binding site at high TST concentrations. Prompted by the observation that interactions between TST and the G'-helix mobilize the ligand into the active site, a free-energy analysis of TST distribution in the membrane was conducted and revealed that the depth of the G'-helix is optimal for extracting TST. In summary, these simulations confirm separate, but adjacent substrate-binding sites within the enzyme and the existence of an auxiliary TST-binding site. The broader impact of these simulations is that they support a mechanism in which cytochromes P450 directly recruit membrane-solubilized substrates.

Inhibition of Human CYP3A4 by Rationally Designed Ritonavir-Like Compounds: Impact and Interplay of the Side Group Functionalities

Structure-function relationships of nine rationally designed ritonavir-like compounds were investigated to better understand the ligand binding and inhibitory mechanism in human drug-metabolizing cytochrome P450 3A4 (CYP3A4). The analogs had a similar backbone and pyridine and tert-butyloxycarbonyl (Boc) as the heme-ligating and terminal groups, respectively. N-Isopropyl, N-cyclopentyl, or N-phenyl were the R₁-side group substituents alone (compounds 5a-c) or in combination with phenyl or indole at the R₂ position (8a-c and 8d-f subseries, respectively). Our experimental and structural data indicate that (i) for all analogs, a decrease in the dissociation constant (K_s) coincides with a decrease in IC₅₀, but no relation with other derived parameters is observed; (ii) an increase in the R₁ volume, hydrophobicity, and aromaticity markedly lowers K_s and IC₅₀, whereas the addition of aromatic R₂ has a more pronounced positive effect on the inhibitory potency than the binding strength; (iii) the ligands' association mode is strongly influenced by the mutually dependent R₁-R₂ interplay, but the R₁-mediated interactions are dominant and define the overall conformation in the active site; (iv) formation of a strong H-bond with Ser119 is a prerequisite for potent CYP3A4 inhibition; and (v) the strongest inhibitor in the series, the R₁-phenyl/R₂-indole containing 8f (K_s and IC₅₀ of 0.08 and 0.43 μM, respectively), is still less potent than ritonavir, even under conditions that prevent the mechanism based inactivation of CYP3A4. Crystallographic data were essential for better understanding and interpretation of the experimental results, and suggested how the inhibitor design could be further optimized.

The Cytochrome P450 3A4 has three Major Conformations: New Clues to Drug Recognition by this Promiscuous Enzyme

We computed the channels of the 3A4 isoform of the cytochrome P450 3A4 (CYP) on the basis of 24 crystal structures extracted from the Protein Data Bank (PDB). We identified three major conformations (denoted C, O1 and O2) using an enhanced version of the CCCPP software that we developed for the present work, while only two conformations (C and O² ) are considered in the literature. We established the flowchart of definition of these three conformations in function of the structural and physicochemical parameters of the ligand. The channels are characterized with qualitative and quantitative parameters, and not only with their surrounding secondary structures as it is usually done in the literature.

Use of bioconjugation with cytochrome P450 enzymes

Bioconjugation, defined as chemical modification of biomolecules, is widely employed in biological and biophysical studies. It can expand functional diversity and enable applications ranging from biocatalysis, biosensing and even therapy. This review summarizes how chemical modifications of cytochrome P450 enzymes (P450s or CYPs) have contributed to improving our understanding of these enzymes. Genetic modifications of P450s have also proven very useful but are not covered in this review. Bioconjugation has served to gain structural information and investigate the mechanism of P450s via photoaffinity labeling, mechanism-based inhibition (MBI) and fluorescence studies. P450 surface acetylation and protein cross-linking have contributed to the investigation of protein complexes formation involving P450 and its redox partner or other P450 enzymes. Finally, covalent immobilization on polymer surfaces or electrodes has benefited the areas of biocatalysis and biosensor design. 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.

High-Level Production and Properties of the Cysteine-Depleted Cytochrome P450 3A4

Human drug-metabolizing cytochrome P450 3A4 (CYP3A4) is a dynamic enzyme with a large and highly malleable active site that can fit structurally diverse compounds. Despite extensive investigations, structure-function relationships and conformational dynamics in CYP3A4 are not fully understood. This study was undertaken to engineer a well-expressed and functionally active cysteine-depleted CYP3A4 that can be used in biochemical and biophysical studies. cDNA codon optimization and screening mutagenesis were utilized to boost the level of bacterial expression of CYP3A4 and identify the least harmful substitutions for all six non-heme-ligating cysteines. The C58A/C64M/C98A/C239T/C377A/C468S (Cys-less) mutant was found to be expressed as highly as the optimized wild-type (opt-WT) CYP3A4. The high-resolution X-ray structures of opt-WT and Cys-less CYP3A4 revealed that gene optimization leads to a different folding in the Phe108 and Phe189 regions and promotes binding of the active site glycerol that interlocks Ser119 and Arg212, critical for ligand association, and the hydrophobic cluster adjacent to Phe108. Crowding and decreased flexibility of the active site, as well as structural alterations observed at the C64M, C239T, and C468S mutational sites, might be responsible for the distinct ligand binding behavior of opt-WT and Cys-less CYP3A4. Nonetheless, the Cys-less mutant could be used for structure-function investigations because it orients bromoergocryptine and ritonavir (a high-affinity substrate and a high-potency inhibitor, respectively) like the WT and has a higher activity toward 7-benzyloxy(4-trifluoromethyl)coumarin.

Allosteric Activation of Cytochrome P450 3A4 via Progesterone Bioconjugation

Human cytochrome P450 3A4 (CYP3A4) is responsible for the metabolism of the majority of drugs. As such, it is implicated in many adverse drug-drug and food-drug interactions, and is of significant interest to the pharmaceutical industry. This enzyme is known to simultaneously bind multiple ligands and display atypical enzyme kinetics, suggestive of allostery and cooperativity. As well, evidence of a postulated peripheral allosteric binding site has provoked debate around its significance and location. We report the use of bioconjugation to study the significance of substrate binding at the proposed allosteric site and its effect on CYP3A4 activity. CYP3A4 mutants were created and covalently modified with various small molecules including progesterone. The labeled mutants displayed enhanced kinetic stability and improved activity in testosterone and 7-benzyloxy-(4-trifluoromethyl)coumarin oxidation assays. Our work applies a new strategy to study cytochrome P450 allostery and supports the hypothesis that substrate binding at the postulated allosteric site of CYP3A4 may induce functional cooperativity.

7,8-benzoflavone binding to human cytochrome P450 3A4 reveals complex fluorescence quenching, suggesting binding at multiple protein sites

Human cytochrome P450 (P450) 3A4 is involved in the metabolism of one-half of marketed drugs and shows cooperative interactions with some substrates and other ligands. The interaction between P450 3A4 and the known allosteric effector 7,8-benzoflavone (α-naphthoflavone, αNF) was characterized using steady-state fluorescence spectroscopy. The binding interaction of P450 3A4 and αNF effectively quenched the fluorescence of both the enzyme and ligand. The Hill Equation and Stern-Volmer fluorescence quenching models were used to evaluate binding of ligand to enzyme. P450 3A4 fluorescence was quenched by titration with αNF; at the relatively higher [αNF]/[P450 3A4] ratios in this experiment, two weaker quenching interactions were revealed (K_d 1.8-2.5 and 6.5 μM). A range is given for the stronger interaction since αNF quenching of P450 3A4 fluorescence changed the protein spectral profile: quenching of 315 nm emission was slightly more efficient (K_d 1.8 μM) than the quenching of protein fluorescence at 335 and 355 nm (K_d 2.5 and 2.1 μM, respectively). In the reverse titration, αNF fluorescence was quenched by P450 3A4; at the lower [αNF]/[P450 3A4] ratios here, two strong quenching interactions were revealed (K_d 0.048 and 1.0 μM). Thus, four binding interactions of αNF to P450 3A4 are suggested by this study, one of which may be newly recognized and which could affect studies of drug oxidations by this important enzyme.

Structural basis for regiospecific midazolam oxidation by human cytochrome P450 3A4

Human cytochrome P450 3A4 (CYP3A4) is a major hepatic and intestinal enzyme that oxidizes more than 60% of administered therapeutics. Knowledge of how CYP3A4 adjusts and reshapes the active site to regioselectively oxidize chemically diverse compounds is critical for better understanding structure-function relations in this important enzyme, improving the outcomes for drug metabolism predictions, and developing pharmaceuticals that have a decreased ability to undergo metabolism and cause detrimental drug-drug interactions. However, there is very limited structural information on CYP3A4-substrate interactions available to date. Despite the vast variety of drugs undergoing metabolism, only the sedative midazolam (MDZ) serves as a marker substrate for the in vivo activity assessment because it is preferentially and regioselectively oxidized by CYP3A4. We solved the 2.7 Å crystal structure of the CYP3A4-MDZ complex, where the drug is well defined and oriented suitably for hydroxylation of the C1 atom, the major site of metabolism. This binding mode requires H-bonding to Ser119 and a dramatic conformational switch in the F-G fragment, which transmits to the adjacent D, E, H, and I helices, resulting in a collapse of the active site cavity and MDZ immobilization. In addition to providing insights on the substrate-triggered active site reshaping (an induced fit), the crystal structure explains the accumulated experimental results, identifies possible effector binding sites, and suggests why MDZ is predominantly metabolized by the CYP3A enzyme subfamily.

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.

Recent Structural Insights into Cytochrome P450 Function

Cytochrome P450 (P450) enzymes are important in the metabolism of drugs, steroids, fat-soluble vitamins, carcinogens, pesticides, and many other types of chemicals. Their catalytic activities are important issues in areas such as drug-drug interactions and endocrine function. During the past 30 years, structures of P450s have been very helpful in understanding function, particularly the mammalian P450 structures available in the past 15 years. We review recent activity in this area, focusing on the past 2 years (2014-2015). Structural work with microbial P450s includes studies related to the biosynthesis of natural products and the use of parasitic and fungal P450 structures as targets for drug discovery. Studies on mammalian P450s include the utilization of information about 'drug-metabolizing' P450s to improve drug development and also to understand the molecular bases of endocrine dysfunction.

Cytochrome P450 Organization and Function Are Modulated by Endoplasmic Reticulum Phospholipid Heterogeneity

Cytochrome P450s (P450s) comprise a superfamily of proteins that catalyze numerous monooxygenase reactions in animals, plants, and bacteria. In eukaryotic organisms, these proteins not only carry out reactions necessary for the metabolism of endogenous compounds, but they are also important in the oxidation of exogenous drugs and other foreign compounds. Eukaryotic P450 system proteins generally reside in membranes, primarily the endoplasmic reticulum or the mitochondrial membrane. These membranes provide a scaffold for the P450 system proteins that facilitate interactions with their redox partners as well as other P450s. This review focuses on the ability of specific lipid components to influence P450 activities, as well as the role of the membrane in P450 function. These studies have shown that P450s and NADPH-cytochrome P450 reductase appear to selectively associate with specific phospholipids and that these lipid-protein interactions influence P450 activities. Finally, because of the heterogeneous nature of the endoplasmic reticulum as well as other biologic membranes, the phospholipids are not arranged randomly but associate to generate lipid microdomains. Together, these characteristics can affect P450 function by 1) altering the conformation of the proteins, 2) influencing the P450 interactions with their redox partners, and 3) affecting the localization of the proteins into specific membrane microdomains.

The use of isomeric testosterone dimers to explore allosteric effects in substrate binding to cytochrome P450 CYP3A4

Cytochrome P450 CYP3A4 is the main drug-metabolizing enzyme in the human liver, being responsible for oxidation of 50% of all pharmaceuticals metabolized by human P450 enzymes. Possessing a large substrate binding pocket, it can simultaneously bind several substrate molecules and often exhibits a complex pattern of drug-drug interactions. In order to better understand structural and functional aspects of binding of multiple substrate molecules to CYP3A4 we used resonance Raman and UV-VIS spectroscopy to document the effects of binding of synthetic testosterone dimers of different configurations, cis-TST2 and trans-TST2. We directly demonstrate that the binding of two steroid molecules, which can assume multiple possible configurations inside the substrate binding pocket of monomeric CYP3A4, can lead to active site structural changes that affect functional properties. Using resonance Raman spectroscopy, we have documented perturbations in the ferric and Fe-CO states by these substrates, and compared these results with effects caused by binding of monomeric TST. While the binding of trans-TST2 yields results similar to those obtained with monomeric TST, the binding of cis-TST2 is much tighter and results in significantly more pronounced conformational changes of the porphyrin side chains and Fe-CO unit. In addition, binding of an additional monomeric TST molecule in the remote allosteric site significantly improves binding affinity and the overall spin shift for CYP3A4 with trans-TST2 dimer bound inside the substrate binding pocket. This result provides the first direct evidence for an allosteric effect of the peripheral binding site at the protein-membrane interface on the functional properties of CYP3A4.