Peripheral ligand-binding site in cytochrome P450 3A4 located with fluorescence resonance energy transfer (FRET)

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.

Equilibrium landscape of ingress/egress channels and gating residues of the Cytochrome P450 3A4

The Cytochrome P450 (CYP) enzymes metabolize a variety of drugs, which may potentially lead to toxicity or reduced efficacy when drugs are co-administered. These drug-drug interactions are often manifested by CYP3A4, the most prevalent of all CYP isozymes. We carried out multiple MD simulations employing CAVER to quantify the channels, and Hidden Markov Models (HMM) to characterize the behavior of the gating residues. We discuss channel properties, bottleneck residues with respect to their likelihood to deem the respective channel ingress or egress, gating residues regarding their open or closed states, and channel location relative to the membrane. Channels do not display coordinated motion and randomly transition between different conformations. Gateway residues also behave in a random fashion. Our findings shed light on the equilibrium behavior of the gating residues and channels in the apo state.

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.

Ni2+ and Cu2+ complexes of N-(2,6-dichlorophenyl)-N-mesityl formamidine dithiocarbamate structural and functional properties as CYP3A4 potential substrates

Metal compounds continued to attract diverse applications due to their malleability in several capacities. In this study, we present our findings on the crystal structures and functional properties of Ni2+ and Cu2+ complexes of N'-(2,6-dichlorophenyl)-N-mesitylformamidine dithiocarbamate (L) comprising [Ni-(L)2] (1) and [Cu-(L)2] (2) with a four-coordinate metal center. We established the two complex structures through 1H and 13C nuclear magnetic resonance (NMR), elemental, and single-crystal X-ray analysis. The analyses showed that the two complexes are isomorphous, having P21/c as a space group and a unit-cell similarity index (π) of 0.002. The two complexes conform to a distorted square planar geometry around the metal centers. The calculated and experimental data, including bond lengths, angles, and NMR values, are similar. Hirshfeld surface analysis revealed the variational contribution of the different types of intermolecular contacts driven by the crystal lattice of the two solvated complexes. Our knowledge of the potential biological implication of these structures enabled us to probe the compounds as prospective CYP3A4 inhibitors. This approach mimics current trends in pharmaceutical design and biomedicine by incorporating potentially active molecules into various media to predict their biological efficacies. The simulations show appreciable binding of compounds 1 and 2 to CYP3A4 with average interaction energies of -97 and -87 kcal/mol, respectively. The protein attains at least five conformational states in the three studied models using a Gaussian Mixture Model-based clustering and free energy prediction. Electric field analysis shows the crucial residues to substrate binding at the active site, enabling CYP3A4 structure to function prediction. The predicted inhibition with these Ni2+ and Cu2+ complexes indicates that CYP3A4 overexpression in a diseased state like cancer would reduce, thereby increasing the chemotherapeutic compounds' shelf-lives for adsorption. This multidimensional study addresses various aspects of molecular metal electronics, including their application as substrate-mimicking inhibitors. The outcome would enable further research on bio-metal compounds of critical potential.

Human cytochrome P450 17A1 structures with metabolites of prostate cancer drug abiraterone reveal substrate-binding plasticity and a second binding site

Abiraterone acetate is a first-line therapy for castration-resistant prostate cancer. This prodrug is deacetylated in vivo to abiraterone, which is a potent and specific inhibitor of cytochrome P450 17A1 (CYP17A1). CYP17A1 performs two sequential steps that are required for the biosynthesis of androgens that drive prostate cancer proliferation, analogous to estrogens in breast cancer. Abiraterone can be further metabolized in vivo on the steroid A ring to multiple metabolites that also inhibit CYP17A1. Despite its design as an active-site-directed substrate analog, abiraterone and its metabolites demonstrate mixed competitive/noncompetitive inhibition. To understand their binding, we solved the X-ray structures of CYP17A1 with three primary abiraterone metabolites. Despite different conformations of the steroid A ring and substituents, all three bound in the CYP17A1 active site with the steroid core packed against the I helix and the A ring C3 keto or hydroxyl oxygen forming a hydrogen bond with N202 similar to abiraterone itself. The structure of CYP17A1 with 3-keto, 5α-abiraterone was solved to 2.0 Å, the highest resolution to date for a CYP17A1 complex. This structure had additional electron density near the F/G loop, which is likely a second molecule of the inhibitor and which may explain the noncompetitive inhibition. Mutation of the adjacent Asn52 to Tyr positions its side chain in this space, maintains enzyme activity, and prevents binding of the peripheral ligand. Collectively, our findings provide further insight into abiraterone metabolite binding and CYP17A1 function.

Crystal Structure of CYP3A4 Complexed with Fluorol Identifies the Substrate Access Channel as a High-Affinity Ligand Binding Site

Cytochrome P450 3A4 (CYP3A4) is a major human drug-metabolizing enzyme, notoriously known for its extreme substrate promiscuity, allosteric behavior, and implications in drug–drug interactions. Despite extensive investigations, the mechanism of ligand binding to CYP3A4 is not fully understood. We determined the crystal structure of CYP3A4 complexed with fluorol, a small fluorescent dye that can undergo hydroxylation. In the structure, fluorol associates to the substrate channel, well suited for the binding of planar polyaromatic molecules bearing polar groups, through which stabilizing H-bonds with the polar channel residues, such as Thr224 and Arg372, can be established. Mutagenesis, spectral, kinetic, and functional data confirmed the involvement but not strict requirement of Thr224 for the association of fluorol. Collectively, our data identify the substrate channel as a high-affinity ligand binding site and support the notion that hydrophobic ligands first dock to the nearby peripheral surface, before migrating to the channel and, subsequently, into the active site.

Vanillic acid attenuates cell proliferation, xenobiotic enzyme activity, and the status of pulmonary mitochondrial enzymes in lung carcinoma

The purpose of the study is to determine the anti-proliferative and mitochondrial status of benzo(a)pyrene-induced lung cancer in Swiss albino mice, as well as the modulatory effect of vanillic acid on it. B(a)P had altered levels of lysosomal enzymes, xenobiotic-metabolizing enzymes, cell proliferation, inflammation, and mitochondrial abnormalities, whereas treatment with VA treatment significantly reversed the aforementioned activities. According to the findings, VA greatly reduces lung carcinogenesis by restoring antioxidants and xenobiotic-enzyme levels, consequently proving to be an anti-proliferative and anti-inflammatory drug against lung cancer in mice. PRACTICAL APPLICATIONS: As we all know, lung cancer is on the rise all over the world. A recent study demonstrated that vanillic acid protects against B(a)P in experimental mice. According to the findings, VA considerably suppresses lung carcinogenesis by restoring lysosomal enzyme levels, xenobiotic-metabolizing enzyme levels, and mitochondrial activities, effectively functioning as an anti-proliferative and anti-inflammatory therapy against lung cancer. According to the most recent study, vanillic acid can be used as a defensive medicine in the treatment of lung cancer.

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.

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.

Identification of the contact region responsible for the formation of the homomeric CYP1A2•CYP1A2 complex

Previous studies showed that cytochrome P450 1A2 (CYP1A2) forms a homomeric complex that influences its metabolic characteristics. Specifically, CYP1A2 activity exhibits a sigmoidal response as a function of NADPH-cytochrome P450 reductase (POR) concentration and is consistent with an inhibitory CYP1A2•CYP1A2 complex that is disrupted by increasing [POR] (Reed et al. (2012) Biochem. J. 446, 489-497). The goal of this study was to identify the CYP1A2 contact regions involved in homomeric complex formation. Examination of X-ray structure of CYP1A2 implicated the proximal face in homomeric complex formation. Consequently, the involvement of residues L91-K106 (P1 region) located on the proximal face of CYP1A2 was investigated. This region was replaced with the homologous region of CYP2B4 (T81-S96) and the protein was expressed in HEK293T/17 cells. Complex formation and its disruption was observed using bioluminescence resonance energy transfer (BRET). The P1-CYP1A2 (CYP1A2 with the modified P1 region) exhibited a decreased BRET signal as compared with wild-type CYP1A2 (WT-CYP1A2). On further examination, P1-CYP1A2 was much less effective at disrupting the CYP1A2•CYP1A2 homomeric complex, when compared with WT-CYP1A2, thereby demonstrating impaired binding of P1-CYP1A2 to WT-CYP1A2 protein. In contrast, the P1 substitution did not affect its ability to form a heteromeric complex with CYP2B4. P1-CYP1A2 also showed decreased activity as compared with WT-CYP1A2, which was consistent with a decrease in the ability of P1-CYP1A2 to associate with WT-POR, again implicating the P1 region in POR binding. These results indicate that the contact region responsible for the CYP1A2•CYP1A2 homomeric complex resides in the proximal region of the protein.

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.

Effects of alcohol-induced increase in CYP2E1 content in human liver microsomes on the activity and cooperativity of CYP3A4

We investigate the effect of the alcohol-induced increase in the content of CYP2E1 in human liver microsomes (HLM) on the function of CYP3A4. Membrane incorporation of the purified CYP2E1 into HLM considerably increases the rate of metabolism of 7-benzyloxyquinoline (BQ) and attenuates the homotropic cooperativity observed with this CYP3A4-specific substrate. It also eliminates the activating effect of α-naphthoflavone (ANF) seen in some HLM samples. To probe the physiological relevance of these effects, we compared three pooled preparations of HLM from normal donors (HLM-N) with a pooled preparation from ten heavy alcohol consumers (HLM-A). The composition of the P450 pool in all samples was characterized by the mass-spectrometric determination of 11 cytochrome P450 species. The fractional content of CYP2E1 in HLM-A was from 2.0 to 3.4 times higher than in HLM-N. In contrast, the content of CYP3A4 in HLM-A was the lowest among all samples. Despite that, HLM-A exhibited a much higher metabolism rate and a lower homotropic cooperativity with BQ, similar to CYP2E1-enriched HLM-N. To substantiate the involvement of interactions between CYP2E1 and CYP3A4 in these effects, we probed hetero-association of these proteins in CYP3A4-containing Supersomes™ with a technique employing CYP2E1 labeled with BODIPY-618 maleimide. These experiments evinced the interactions between the two enzymes and revealed an inhibitory effect of ANF on their association. Our results demonstrate that the functional properties of CYP3A4 are fundamentally dependent on the composition of the cytochrome P450 ensemble and suggest a possible impact of chronic alcohol exposure on the pharmacokinetics of drugs metabolized by CYP3A4.

Dynamics and Location of the Allosteric Midazolam Site in Cytochrome P4503A4 in Lipid Nanodiscs

Promiscuous and allosteric drug interactions with cytochrome P450 3A4 (CYP3A4) are ubiquitous but incompletely understood at the molecular level. A classic allosteric CYP3A4 drug interaction includes the benzodiazepine midazolam (MDZ). MDZ exhibits homotropic and heterotropic allostery when metabolized to 1'-hydroxy and 4-hydroxy metabolites in varying ratios. The combination of hydrogen-deuterium exchange mass spectrometry (HDX-MS) and Gaussian accelerated molecular dynamics (GaMD) simulations of CYP3A4 in lipid nanodiscs and in a lipid bilayer, respectively, reveals MDZ-dependent changes in dynamics in a membrane environment. The F-, G-, and intervening helices, as well as the loop preceding the β1-sheets, display the largest observed changes in HDX. The GaMD suggests a potential allosteric binding site for MDZ in the F'- and G'-regions, which undergo significant increases in HDX at near-saturating MDZ concentrations. The HDX-MS and GaMD results confirm that changes in dynamics are most significant near the developing consensus allosteric site, and these changes are distinct from those observed previously with the nonallosteric inhibitor ketoconazole. The results suggest that the allosteric MDZ remains mobile in its binding site at the Phe-cluster. The results further suggest that this binding site remains dynamic or changes the depth of insertion in the membrane.

So many roads traveled: A career in science and administration

I have traveled many roads during my career. After spending my first 19 years in Los Angeles, I became somewhat of an academic nomad, studying and/or working in six universities in the United States and three in Sweden. In chronological order, I have a B.A. in Scandinavian languages and literature from UCLA, a Ph.D. in biochemistry from Uppsala University, and an M.S. in toxicology from the Karolinska Institute. I have been in schools of natural science, pharmacy, and medicine and have worked in multiple basic science departments and one clinical department. I have served as a research-track and tenured faculty member, department chair, associate dean, and dean. My research has spanned toxinology, biochemistry, toxicology, and pharmacology. Through all the moves, I have gained much and lost some. For the past 40 years, my interest has been cytochrome P450 structure-function and structure-activity relationships. My lab has focused on CYP2B enzymes using X-ray crystallography, site-directed mutagenesis, deuterium-exchange MS, isothermal titration calorimetry, and computational methods in conjunction with a variety of functional assays. This comprehensive approach has enabled detailed understanding of the structural basis of the remarkable substrate promiscuity of CYP2B enzymes. We also have investigated the mechanisms of CYP3A4 allostery using biophysical and advanced spectroscopic techniques, and discovered a pivotal role of P450-P450 interactions and of multiple-ligand binding. A major goal of this article is to provide lessons that may be useful to scientists in the early and middle stages of their careers and those more senior scientists contemplating an administrative move.

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.

Allosteric Interactions in Human Cytochrome P450 CYP3A4: The Role of Phenylalanine 213

The role of Phe213 in the allosteric mechanism of human cytochrome P450 CYP3A4 was studied using a combination of progesterone (PGS) and carbamazepine (CBZ) as probe substrates. We expressed, purified, and incorporated into POPC Nanodiscs three mutants, F213A, F213S, and F213Y, and compared them with wild-type (WT) CYP3A4 by monitoring spectral titration, the rate of NADPH oxidation, and steady-state product turnover rates with pure substrates and substrate mixtures. All mutants demonstrated higher activity with CBZ, lower activity with PGS, and a reduced level of activation of CBZ epoxidation by PGS, which was most pronounced in the F213A mutant. Using all-atom molecular dynamics simulations, we compared the dynamics of WT CYP3A4 and the F213A mutant incorporated into the lipid bilayer and the effect of the presence of the PGS molecule at the allosteric peripheral site and evaluated the critical role of Phe213 in mediating the heterotropic allosteric interactions in CYP3A4.

Cytochrome P450 3A Enzymes Are Key Contributors for Hepatic Metabolism of Bufotalin, a Natural Constitute in Chinese Medicine Chansu

Bufotalin (BFT), one of the naturally occurring bufodienolides, has multiple pharmacological and toxicological effects including antitumor activity and cardiotoxicity. This study aimed to character the metabolic pathway(s) of BFT and to identify the key drug metabolizing enzyme(s) responsible for hepatic metabolism of BFT in human, as well as to explore the related molecular mechanism of enzymatic selectivity. The major metabolite of BFT in human liver microsomes (HLMs) was fully identified as 5β-hydroxylbufotalin by LC-MS/MS and NMR techniques. Reaction phenotyping and chemical inhibition assays showed that CYP3A4 and CYP3A5 were key enzymes responsible for BFT 5β-hydroxylation. Kinetic analyses demonstrated that BFT 5β-hydroxylation in both HLMs and human CYP3A4 followed the biphasic kinetics, while BFT 5β-hydroxylation in CYP3A5 followed substrate inhibition kinetics. Furthermore, molecular docking simulations showed that BFT could bind on two different ligand-binding sites on both CYP3A4 and CYP3A5, which partially explained the different kinetic behaviors of BFT in CYP3A4 and CYP3A5. These findings are very helpful for elucidating the phase I metabolism of BFT in human and for deeper understanding the key interactions between CYP3A enzymes and bufadienolides, as well as for the development of bufadienolide-type drugs with improved pharmacokinetic and safety profiles.

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.

Binding Specificity Determines the Cytochrome P450 3A4 Mediated Enantioselective Metabolism of Metconazole

Cytochrome P450 3A4 (CYP3A4) is a promiscuous enzyme, mediating the biotransformations of ∼50% of clinically used drugs, many of which are chiral molecules. Probing the interactions between CYP3A4 and chiral chemicals is thus essential for the elucidation of molecular mechanisms of enantioselective metabolism. We developed a stepwise-restrained-molecular-dynamics (MD) method to model human CYP3A4 in a complex with cis-metconazole (MEZ) isomers and performed conventional MD simulations with a total simulation time of 2.2 μs to probe the molecular interactions. Our current study, which employs a combined experimental and theoretical approach, reports for the first time on the distinct conformational changes of CYP3A4 that are induced by the enantioselective binding of cis-MEZ enantiomers. CYP3A4 preferably metabolizes cis-RS MEZ over the cis-SR isomer, with the resultant enantiomer fraction for cis-MEZ increasing rapidly from 0.5 to 0.82. cis-RS MEZ adopts a more extended structure in the active pocket with its Cl atom exposed to the solvent, whereas cis-SR MEZ sits within the hydrophobic core of the active pocket. Free-energy-perturbation calculations indicate that unfavorable van der Waals interactions between the cis-MEZ isomers and the CYP3A4 binding pocket predominantly contribute to their binding-affinity differences. These results demonstrate that binding specificity determines the cytochrome P450 3A4 mediated enantioselective metabolism of cis-MEZ.

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.

Physical Studies of P450-P450 Interactions: Predicting Quaternary Structures of P450 Complexes in Membranes from Their X-ray Crystal Structures

Cytochrome P450 enzymes, which catalyze oxygenation reactions of both exogenous and endogenous chemicals, are membrane bound proteins that require interaction with their redox partners in order to function. Those responsible for drug and foreign compound metabolism are localized primarily in the endoplasmic reticulum of liver, lung, intestine, and other tissues. More recently, the potential for P450 enzymes to exist as supramolecular complexes has been shown by the demonstration of both homomeric and heteromeric complexes. The P450 units in these complexes are heterogeneous with respect to their distribution and function, and the interaction of different P450s can influence P450-specific metabolism. The goal of this review is to examine the evidence supporting the existence of physical complexes among P450 enzymes. Additionally, the review examines the crystal lattices of different P450 enzymes derived from X-ray diffraction data to make assumptions regarding possible quaternary structures in membranes and in turn, to predict how the quaternary structures could influence metabolism and explain the functional effects of specific P450-P450 interactions.

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.

A FRET-based method for monitoring septin polymerization and binding of septin-associated proteins

Much about septin function has been inferred from in vivo studies using mainly genetic methods, and much of what we know about septin organization has been obtained through examination of static structures in vitro primarily by electron microscopy. Deeper mechanistic insight requires real-time analysis of the dynamics of the assembly of septin-based structures and how other proteins associate with them. We describe here a Förster resonance energy transfer (FRET)-based approach for measuring in vitro the rate and extent of filament formation from septin complexes, binding of other proteins to septin structures, and the apparent affinities of these interactions. FRET is particularly well suited for interrogating protein-protein interactions, especially on a rapid timescale; the spectral change provides an unambiguous indication of whether two elements within the system under study are associating and serves as a molecular-level "ruler" because it is very sensitive to the separation between the donor and acceptor fluorophores over biologically relevant distances (≤10nm). The necessary procedures involve generation of appropriate cysteine-less and single cysteine-containing septin variants, expression and purification of the heterooctameric complexes containing them, efficient labeling of the purified complexes with desired fluorophores, fluorimetric measurement of FRET, and appropriate safeguards and controls in data acquisition and analysis. Our methods can be used to interrogate the effects of buffer conditions, small molecules, and septin-binding proteins on septin filament assembly or stability; determine the effect of alternative septin subunits, mutational alterations, or posttranslational modifications on assembly; and, delineate the location of septin-binding proteins.

Functional importance of a peripheral pocket in mammalian cytochrome P450 2B enzymes

The functional importance of a peripheral pocket found in previously published X-ray crystal structures of CYP2B4 and CYP2B6 was probed using a biophysical approach. Introduction of tryptophan within the pocket of CYP2B4 at F202 or I241 leads to marked impairment of 7-ethoxy-4-(trifluoromethyl)coumarin (7-EFC) or 7-benzyloxyresorufin O-dealkylation efficiency; a similar substitution at F195, near the surface access to the pocket, does not affect these activities. The analogous CYP2B6 F202W mutant is inactive in the 7-EFC O-dealkylation assay. The stoichiometry of 7-EFC deethylation suggested that the decreased activity of F202W and I241W in CYP2B4 and lack of activity of F202W in CYP2B6 coincided with a sharp increase in the flux of reducing equivalents through the oxidase shunt to produce excess water. The results indicate that the chemical identity of residues within this peripheral pocket, but not at the mouth of the pocket, is important in substrate turnover and redox coupling, likely through effects on active site topology.

A Combined Molecular Docking/Dynamics Approach to Probe the Binding Mode of Cancer Drugs with Cytochrome P450 3A4

Cytarabine, daunorubicin, doxorubicin and vincristine are clinically used for combinatorial therapies of cancers in different combinations. However, the knowledge about the interaction of these drugs with the metabolizing enzyme cytochrome P450 is limited. Therefore, we utilized computational methods to predict and assess the drug-binding modes. In this study, we performed docking, MD simulations and free energy landscape analysis to understand the drug-enzyme interactions, protein domain motions and the most populated free energy minimum conformations of the docked protein-drug complexes, respectively. The outcome of docking and MD simulations predicted the productive, as well as the non-productive binding modes of the selected drugs. Based on these interaction studies, we observed that S119, R212 and R372 are the major drug-binding residues in CYP3A4. The molecular mechanics Poisson–Boltzmann surface area analysis revealed the dominance of hydrophobic forces in the CYP3A4-drug association. Further analyses predicted the residues that may contain favorable drug-specific interactions. The probable binding modes of the cancer drugs from this study may extend the knowledge of the protein-drug interaction and pave the way to design analogs with reduced toxicity. In addition, they also provide valuable insights into the metabolism of the cancer drugs.

Anion-Dependent Stimulation of CYP3A4 Monooxygenase

We co-crystallized human cytochrome P450 3A4 (CYP3A4) with progesterone (PRG) under two different conditions, but the resulting complexes contained only one PRG molecule bound to the previously identified peripheral site. A novel feature in one of our structures is a citrate ion, originating from the crystallization solution. The citrate-binding site is located in an area where the N-terminus splits from the protein core and, thus, is suitable for the interaction with the anionic phospholipids of the microsomal membrane. We investigated how citrate affects the function of a soluble CYP3A4 monooxygenase system consisting of equimolar amounts of CYP3A4 and cytochrome P450 reductase (CPR). Citrate was found to affect the properties of both redox partners and stimulated their catalytic activities in a concentration-dependent manner via a complex mechanism. CYP3A4-substrate binding, reduction of CPR with NADPH, and interflavin and interprotein electron transfer were identified as citrate-sensitive steps. Comparative analysis of various negatively charged organic compounds indicated that, in addition to alterations caused by changes in ionic strength, anions modulate the properties of CYP3A4 and CPR through specific anion-protein interactions. Our results help to better understand previous observations and provide new mechanistic insights into CYP3A4 function.

Interactions between CYP3A4 and Dietary Polyphenols

The human cytochrome P450 enzymes (P450s) catalyze oxidative reactions of a broad spectrum of substrates and play a critical role in the metabolism of xenobiotics, such as drugs and dietary compounds. CYP3A4 is known to be the main enzyme involved in the metabolism of drugs and most other xenobiotics. Dietary compounds, of which polyphenolics are the most studied, have been shown to interact with CYP3A4 and alter its expression and activity. Traditionally, the liver was considered the prime site of CYP3A-mediated first-pass metabolic extraction, but in vitro and in vivo studies now suggest that the small intestine can be of equal or even greater importance for the metabolism of polyphenolics and drugs. Recent studies have pointed to the role of gut microbiota in the metabolic fate of polyphenolics in human, suggesting their involvement in the complex interactions between dietary polyphenols and CYP3A4. Last but not least, all the above suggests that coadministration of drugs and foods that are rich in polyphenols is expected to stimulate undesirable clinical consequences. This review focuses on interactions between dietary polyphenols and CYP3A4 as they relate to structural considerations, food-drug interactions, and potential negative consequences of interactions between CYP3A4 and polyphenols.

Mechanism of drug-drug interactions mediated by human cytochrome P450 CYP3A4 monomer

Using Nanodiscs, we quantitate the heterotropic interaction between two different drugs mediated by monomeric CYP3A4 incorporated into a nativelike membrane environment. The mechanism of this interaction is deciphered by global analysis of multiple-turnover experiments performed under identical conditions using the pure substrates progesterone (PGS) and carbamazepine (CBZ) and their mixtures. Activation of CBZ epoxidation and simultaneous inhibition of PGS hydroxylation are measured and quantitated through differences in their respective affinities for both a remote allosteric site and the productive catalytic site near the heme iron. Preferred binding of PGS at the allosteric site and a stronger preference for CBZ binding at the productive site give rise to a nontrivial drug-drug interaction. Molecular dynamics simulations indicate functionally important conformational changes caused by PGS binding at the allosteric site and by two CBZ molecules positioned inside the substrate binding pocket. Structural changes involving Phe-213, Phe-219, and Phe-241 are thought to be responsible for the observed synergetic effects and positive allosteric interactions between these two substrates. Such a mechanism is likely of general relevance to the mutual heterotropic effects caused by biologically active compounds that exhibit different patterns of interaction with the distinct allosteric and productive sites of CYP3A4, as well as other xenobiotic metabolizing cytochromes P450 that are also involved in drug-drug interactions. Importantly, this work demonstrates that a monomeric CYP3A4 can display the full spectrum of activation and cooperative effects that are observed in hepatic membranes.

Interactions among cytochromes P450 in microsomal membranes: oligomerization of cytochromes P450 3A4, 3A5, and 2E1 and its functional consequences

The body of evidence of physiologically relevant P450-P450 interactions in microsomal membranes continues to grow. Here we probe oligomerization of human CYP3A4, CYP3A5, and CYP2E1 in microsomal membranes. Using a technique based on luminescence resonance energy transfer, we demonstrate that all three proteins are subject to a concentration-dependent equilibrium between the monomeric and oligomeric states. We also observed the formation of mixed oligomers in CYP3A4/CYP3A5, CYP3A4/CYP2E1, and CYP3A5/CYP2E1 pairs and demonstrated that the association of either CYP3A4 or CYP3A5 with CYP2E1 causes activation of the latter enzyme. Earlier we hypothesized that the intersubunit interface in CYP3A4 oligomers is similar to that observed in the crystallographic dimers of some microsomal drug-metabolizing cytochromes P450 (Davydov, D. R., Davydova, N. Y., Sineva, E. V., Kufareva, I., and Halpert, J. R. (2013) Pivotal role of P450-P450 interactions in CYP3A4 allostery: the case of α-naphthoflavone. Biochem. J. 453, 219-230). Here we report the results of intermolecular cross-linking of CYP3A4 oligomers with thiol-reactive bifunctional reagents as well as the luminescence resonance energy transfer measurements of interprobe distances in the oligomers of labeled CYP3A4 single-cysteine mutants. The results provide compelling support for the physiological relevance of the dimer-specific peripheral ligand-binding site observed in certain CYP3A4 structures. According to our interpretation, these results reveal an important general mechanism that regulates the activity and substrate specificity of the cytochrome P450 ensemble through interactions between multiple P450 species. As a result of P450-P450 cross-talk, the catalytic properties of the cytochrome P450 ensemble cannot be predicted by simple summation of the properties of the individual P450 species.

Allosteric modulation of substrate motion in cytochrome P450 3A4-mediated xylene oxidation

Many cytochrome P450 enzymes (CYPs) exhibit allosteric behavior reflecting a complex ligand-binding process involving numerous factors: conformational selection, protein-protein interactions, substrate/effector/protein structure, and multiple-ligand binding. The interplay of CYP plasticity and rigidity contributes to substrate/product selectivity and to allosterism. Detailed evidence describing how protein motion modulates product selectivity is incomplete as are descriptions of effector-induced modulation of substrate dynamics. Our intent was to discover details of allosteric behavior and CYP3A4 flexibility and rigidity by investigating substrate motion using low-molecular weight ligands. Steady state kinetics and product ratios were measured for oxidation of m-xylene-(2)H3 and p-xylene; intramolecular isotope effects were measured for m-xylene-(2)H3 oxidation as a function of m-xylene-(2)H3 and p-xylene concentration. Biphasic kinetic plots indicated homotropic cooperative behavior with xylene isomers. Selectivity for aromatic hydroxylation over benzylic hydroxylation of m-xylene-(2)H3 supports a model in which the region near the CYP3A4 active oxidizing species limits substrate dynamics. p-Xylene impedes the motion of m-xylene-(2)H3 substrates that have access to the active oxidizing species: (kH/kD)obs values for m-xylene-(2)H3 decreased with p-xylene concentration. m-Xylene-(2)H3 and p-xylene do not have simultaneous access to the active oxidizing species: deuterium-labeled and unlabeled p-xylene exhibited similar effects on the (kH/kD)obs values for m-xylene-(2)H3 oxidation. p-Xylene and m-xylene-(2)H3 bind at different sites: m-xylene-(2)H3 oxidation rates and product selectivity were consistent across the p-xylene concentration range. Overall, this study indicates that the intramolecular isotope effect experimental design provides a unique opportunity to investigate allosteric mechanisms as it provides information about substrate motion when the enzyme is primed to oxidize substrates.

The role of cytochrome P450 2B6 and 2B4 substrate access channel residues predicted based on crystal structures of the amlodipine complexes

Recent X-ray crystal structures of human cytochrome P450 2B6 and rabbit cytochrome P450 2B4 in complex with amlodipine showed two bound ligand molecules, one in the active site and one in the substrate access channel. Based on the X-ray crystal structures, we investigated the interactions of P450 2B4 and 2B6 with amlodipine using absorbance spectroscopy, and determined the steady-state kinetics of 7-ethoxy-4-(trifluoromethyl)coumarin and 7-benzyloxyresorufin oxidation by some access channel mutants to evaluate the functional role of these residues in substrate turnover. The results of absorbance titrations are consistent with a simple mechanism with two parallel binding events that result in the formation of the enzyme complex with two molecules of amlodipine. Using this model we were able to resolve two separate ligand-binding events, which are characterized by two distinct KD values in each enzyme. The access channel mutants R73K in P450 2B6 and R73K, V216W, L219W, and F220W in P450 2B4 showed a significant decrease in kcat/KM with the both substrates. Overall, the results suggest that P450 2B4 and 2B6 form an enzyme complex with two molecules of amlodipine in solution, and R73, V216, L219 and F220 in P450 2B4 may play an important role in substrate metabolism.

A large-scale allosteric transition in cytochrome P450 3A4 revealed by luminescence resonance energy transfer (LRET)

Effector-induced allosteric transitions in cytochrome P450 3A4 (CYP3A4) were investigated by luminescence resonance energy transfer (LRET) between two SH-reactive probes attached to various pairs of distantly located cysteine residues, namely the double-cysteine mutants CYP3A4(C64/C468), CYP3A4(C377/C468) and CYP3A4(C64/C121). Successive equimolar labeling of these proteins with the phosphorescent probe erythrosine iodoacetamide (donor) and the near-infrared fluorophore DY-731 maleimide (acceptor) allowed us to establish donor/acceptor pairs sensitive to conformational motions. The interactions of all three double-labeled mutants with the allosteric activators α-naphthoflavone and testosterone resulted in an increase in the distance between the probes. A similar effect was elicited by cholesterol. These changes in distance vary from 1.3 to 8.5 Å, depending on the position of the donor/acceptor pair and the nature of the effector. In contrast, the changes in the interprobe distance caused by such substrates as bromocriptine or 1-pyrenebutanol were only marginal. Our results provide a decisive support to the paradigm of allosteric modulation of CYP3A4 and indicate that the conformational transition caused by allosteric effectors increases the spatial separation between the beta-domain of the enzyme (bearing residues Cys₆₄ and Cys₃₇₇) and the alpha-domain, where Cys₁₂₁ and Cys₄₆₈ are located.

Identifying cytochrome p450 functional networks and their allosteric regulatory elements

Cytochrome P450 (CYP) enzymes play key roles in drug metabolism and adverse drug-drug interactions. Despite tremendous efforts in the past decades, essential questions regarding the function and activity of CYPs remain unanswered. Here, we used a combination of sequence-based co-evolutionary analysis and structure-based anisotropic thermal diffusion (ATD) molecular dynamics simulations to detect allosteric networks of amino acid residues and characterize their biological and molecular functions. We investigated four CYP subfamilies (CYP1A, CYP2D, CYP2C, and CYP3A) that are involved in 90% of all metabolic drug transformations and identified four amino acid interaction networks associated with specific CYP functionalities, i.e., membrane binding, heme binding, catalytic activity, and dimerization. Interestingly, we did not detect any co-evolved substrate-binding network, suggesting that substrate recognition is specific for each subfamily. Analysis of the membrane binding networks revealed that different CYP proteins adopt different membrane-bound orientations, consistent with the differing substrate preference for each isoform. The catalytic networks were associated with conservation of catalytic function among CYP isoforms, whereas the dimerization network was specific to different CYP isoforms. We further applied low-temperature ATD simulations to verify proposed allosteric sites associated with the heme-binding network and their role in regulating metabolic fate. Our approach allowed for a broad characterization of CYP properties, such as membrane interactions, catalytic mechanisms, dimerization, and linking these to groups of residues that can serve as allosteric regulators. The presented combined co-evolutionary analysis and ATD simulation approach is also generally applicable to other biological systems where allostery plays a role.

Pivotal role of P450-P450 interactions in CYP3A4 allostery: the case of α-naphthoflavone

We investigated the relationship between oligomerization of CYP3A4 (cytochrome P450 3A4) and its response to ANF (α-naphthoflavone), a prototypical heterotropic activator. The addition of ANF resulted in over a 2-fold increase in the rate of CYP3A4-dependent debenzylation of 7-BFC [7-benzyloxy-4-(trifluoromethyl)coumarin] in HLM (human liver microsomes), but failed to produce activation in BD Supersomes or Baculosomes containing recombinant CYP3A4 and NADPH-CPR (cytochrome P450 reductase). However, incorporation of purified CYP3A4 into Supersomes containing only recombinant CPR reproduced the behaviour observed with HLM. The activation in this system was dependent on the surface density of the enzyme. Although no activation was detectable at an L/P (lipid/P450) ratio ≥750, it reached 225% at an L/P ratio of 140. To explore the relationship between this effect and CYP3A4 oligomerization, we probed P450-P450 interactions with a new technique that employs LRET (luminescence resonance energy transfer). The amplitude of LRET in mixed oligomers of the haem protein labelled with donor and acceptor fluorophores exhibited a sigmoidal dependence on the surface density of CYP3A4 in Supersomes™. The addition of ANF eliminated this sigmoidal character and increased the degree of oligomerization at low enzyme concentrations. Therefore the mechanisms of CYP3A4 allostery with ANF involve effector-dependent modulation of P450-P450 interactions.

Heterotropic activation of the midazolam hydroxylase activity of CYP3A by a positive allosteric modulator of mGlu5: in vitro to in vivo translation and potential impact on clinically relevant drug-drug interactions

Allosteric modulation of G protein-coupled receptors has gained considerable attention in the drug discovery arena because it opens avenues to achieve greater selectivity over orthosteric ligands. We recently identified a series of positive allosteric modulators (PAMs) of metabotropic glutamate receptor 5 (mGlu(5)) for the treatment of schizophrenia that exhibited robust heterotropic activation of CYP3A4 enzymatic activity. The prototypical compound from this series, 5-(4-fluorobenzyl)-2-((3-fluorophenoxy)methyl)-4,5,6,7-tetrahydropyrazolo[1,5-a]pyrazine (VU0448187), was found to activate CYP3A4 to >100% of its baseline intrinsic midazolam (MDZ) hydroxylase activity in vitro; activation was CYP3A substrate specific and mGlu(5) PAM dependent. Additional studies revealed the concentration-dependence of CYP3A activation by VU0448187 in multispecies hepatic and intestinal microsomes and hepatocytes, as well as a diminished effect observed in the presence of ketoconazole. Kinetic analyses of the effect of VU0448187 on MDZ metabolism in recombinant P450 or human liver microsomes resulted in a significant increase in V(max) (minimal change in K(m)) and required the presence of cytochrome b5. The atypical kinetics translated in vivo, as rats receiving an intraperitoneal administration of VU0448187 prior to MDZ treatment demonstrated a significant increase in circulating 1- and 4-hydroxy- midazolam (1-OH-MDZ, 4-OH-MDZ) levels compared with rats administered MDZ alone. The discovery of a potent substrate-selective activator of rodent CYP3A with an in vitro to in vivo translation serves to illuminate the impact of increasing intrinsic enzymatic activity of hepatic and extrahepatic CYP3A in rodents, and presents the basis to build models capable of framing the clinical relevance of substrate-dependent heterotropic activation.

Protein mechanics: how force regulates molecular function

BACKGROUND

Regulation of proteins is ubiquitous and vital for any organism. Protein activity can be altered chemically, by covalent modifications or non-covalent binding of co-factors. Mechanical forces are emerging as an additional way of regulating proteins, by inducing a conformational change or by partial unfolding.

SCOPE

We review some advances in experimental and theoretical techniques to study protein allostery driven by mechanical forces, as opposed to the more conventional ligand driven allostery. In this respect, we discuss recent single molecule pulling experiments as they have substantially augmented our view on the protein allostery by mechanical signals in recent years. Finally, we present a computational analysis technique, Force Distribution Analysis, that we developed to reveal allosteric pathways in proteins.

MAJOR CONCLUSIONS

Any kind of external perturbation, being it ligand binding or mechanical stretching, can be viewed as an external force acting on the macromolecule, rendering force-based experimental or computational techniques, a very general approach to the mechanics involved in protein allostery.

GENERAL SIGNIFICANCE

This unifying view might aid to decipher how complex allosteric protein machineries are regulated on the single molecular level.

Characterizing the membrane-bound state of cytochrome P450 3A4: structure, depth of insertion, and orientation

Cytochrome P450 3A4 (CYP3A4) is the most abundant membrane-associated isoform of the P450 family in humans and is responsible for biotransformation of more than 50% of drugs metabolized in the body. Despite the large number of crystallographic structures available for CYP3A4, no structural information for its membrane-bound state at an atomic level is available. In order to characterize binding, depth of insertion, membrane orientation, and lipid interactions of CYP3A4, we have employed a combined experimental and simulation approach in this study. Taking advantage of a novel membrane representation, highly mobile membrane mimetic (HMMM), with enhanced lipid mobility and dynamics, we have been able to capture spontaneous binding and insertion of the globular domain of the enzyme into the membrane in multiple independent, unbiased simulations. Despite different initial orientations and positions of the protein in solution, all the simulations converged into the same membrane-bound configuration with regard to both the depth of membrane insertion and the orientation of the enzyme on the surface of the membrane. In tandem, linear dichroism measurements performed on CYP3A4 bound to Nanodisc membranes were used to characterize the orientation of the enzyme in its membrane-bound form experimentally. The heme tilt angles measured experimentally are in close agreement with those calculated for the membrane-bound structures resulted from the simulations, thereby verifying the validity of the developed model. Membrane binding of the globular domain in CYP3A4, which appears to be independent of the presence of the transmembrane helix of the full-length enzyme, significantly reshapes the protein at the membrane interface, causing conformational changes relevant to access tunnels leading to the active site of the enzyme.

Structural diversity of eukaryotic membrane cytochrome p450s

X-ray crystal structures are available for 29 eukaryotic microsomal, chloroplast, or mitochondrial cytochrome P450s, including two non-monooxygenase P450s. These structures provide a basis for understanding structure-function relations that underlie their distinct catalytic activities. Moreover, structural plasticity has been characterized for individual P450s that aids in understanding substrate binding in P450s that mediate drug clearance.