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.

Functional characterization of cytochromes P450 2B from the desert woodrat Neotoma lepida

Mammalian detoxification processes have been the focus of intense research, but little is known about how wild herbivores process plant secondary compounds, many of which have medicinal value or are drugs. cDNA sequences that code for three enzymes of the cytochrome P450 (CYP) 2B subfamily, here termed 2B35, 2B36, and 2B37 have been recently identified from a wild rodent, the desert woodrat (Malenke et al., 2012). Two variant clones of each enzyme were engineered to increase protein solubility and to facilitate purification, as reported for CYP2B enzymes from multiple species. When expressed in Escherichia coli each of the woodrat proteins gave the characteristic maximum at 450nm in a reduced carbon monoxide difference spectrum but generally expressed at lower levels than rat CYP2B1. Two enzymes, 2B36 and 2B37, showed dealkylation activity with the model substrates 7-ethoxy-4-(trifluoromethyl)coumarin and 7-benzyloxyresorufin, whereas 2B35 was inactive. Binding of the monoterpene (+)-α-pinene produced a Type I shift in the absorbance spectrum of each enzyme. Mutation of 2B37 at residues 114, 262, or 480, key residues governing ligand interactions with other CYP2B enzymes, did not significantly change expression levels or produce the expected functional changes. In summary, two catalytic and one ligand-binding assay are sufficient to distinguish among CYP2B35, 2B36, and 2B37. Differences in functional profiles between 2B36 and 2B37 are partially explained by changes in substrate recognition site residue 114, but not 480. The results advance our understanding of the mechanisms of detoxification in wild mammalian herbivores and highlight the complexity of this system.

An in vivo assay for elucidating the importance of cytochromes P450 for the ability of a wild mammalian herbivore (Neotoma lepida) to consume toxic plants

An in vivo assay using the cytochrome P450 (P450) suicide inhibitor 1-aminobenzotriazole (ABT) and 24-h food intake was developed to determine the relative importance of P450s in two populations of Neotoma lepida with respect to biotransformation of plant secondary compounds in the animals' natural diets. The efficacy of ABT as a P450 inhibitor was first validated using hypnotic-state assays with and without pretreatment with ABT. Pretreatment with 100 mg/kg ABT by gavage increased hexobarbital sleep times 3-4-fold in both populations, showing effective inhibition of P450s in woodrats. Next, the Great Basin population was fed a terpene-rich juniper diet, and the Mojave population was fed a phenolic-rich creosote diet, with rabbit chow serving as the control diet in each group. Treatment with ABT inhibited food intake in the Great Basin population fed the juniper diet to a greater extent (35%) than the Great Basin population fed the control diet (19%) or the Mojave population fed the creosote diet (16%). The food intake of the Mojave population fed the control diet was not significantly inhibited by ABT. The findings suggest that the biotransformation of terpenes in juniper relies more heavily on P450s than that of phenolics in creosote. This assay provides an inexpensive and noninvasive method to explore the relative importance of P450s in the biotransformation strategies of wild mammalian herbivores.

Structural and thermodynamic basis of (+)-α-pinene binding to human cytochrome P450 2B6

Despite recent advances in atomic-level understanding of drug and inhibitor interactions with human cytochromes P450, the decades-old questions of chemical and structural determinants of hydrocarbon binding are still unanswered. (+)-α-Pinene is a monoterpene hydrocarbon that is widely distributed in the environment and a potent P450 2B inhibitor. Therefore, a combined biophysical and structural analysis of human P450 2B6 interactions with (+)-α-pinene was undertaken to elucidate the basis of the very high affinity binding. Binding of (+)-α-pinene to the P450 active site was demonstrated by a Type I spectral shift. Thermodynamics of ligand binding were explored using isothermal titration calorimetry and compared to those of P450 2A6, which is much less flexible than 2B6 based on comparison of multiple X-ray crystal structures. Consistent with expectation, entropy is the major driving force for hydrocarbon binding to P450 2A6, as evidenced by the calorimetric results. However, formation of the 2B6-(+)-α-pinene complex has a significant enthalpic component. A 2.0 Å resolution crystal structure of this enzyme-ligand complex reveals that the highly plastic 2B6 utilizes previously unrecognized rearrangements of protein motifs. The results indicate that the specific components of enthalpic contribution to ligand binding are closely tied to the degree of enzyme flexibility.

A structural snapshot of CYP2B4 in complex with paroxetine provides insights into ligand binding and clusters of conformational states

An X-ray crystal structure of CYP2B4 in complex with the drug paroxetine [(3S,4R)-3-[(2H-1,3-benzodioxol-5-yloxy)methyl]-4-(4-fluorophenyl)piperidine] was solved at 2.14 Å resolution. The structure revealed a conformation intermediate to that of the recently solved complex with amlodipine and that of the more compact complex with 4-(4-chlorophenyl)imidazole in terms of the placement of the F-G cassette. Moreover, comparison of the new structure with 15 previously solved structures of CYP2B4 revealed some new insights into the determinants of active-site size and shape. The 2B4-paroxetine structure is nearly superimposable on a previously solved closed structure in a ligand-free state. Despite the overall conformational similarity among multiple closed structures, the active-site cavity volume of the paroxetine complex is enlarged. Further analysis of the accessible space and binding pocket near the heme reveals a new subchamber that resulted from the movement of secondary structural elements and rearrangements of active-site side chains. Overall, the results from the comparison of all 16 structures of CYP2B4 demonstrate a cluster of protein conformations that were observed in the presence or absence of various ligands.

CYP261 enzymes from deep sea bacteria: a clue to conformational heterogeneity in cytochromes P450

We have explored the adaptation of the cytochromes P450 (P450) of deep-sea bacteria to high hydrostatic pressures. Strict conservation of the protein fold and functional importance of protein-bound water make P450 a unique subject for the studies of high-pressure adaptation. Earlier, we expressed and purified a fatty-acid binding P450 from the deep-sea bacteria Photobacterium profundum SS9 (CYP261C1). Here, we report purification and initial characterization of its mesophilic ortholog from the shallow-water P. profundum 3TCK (CYP261C2), as well as another piezophilic enzyme, CYP261D1, from deep-sea Moritella sp. PE36. Comparison of the three enzymes revealed a striking peculiarity of the piezophilic enzymes. Both CYP261C1 and CYP261D1 possess an apparent pressure-induced conformational toggle actuated at the pressures commensurate with the physiological pressure of habitation of the host bacteria. Furthermore, in contrast to CYP261C2, the piezophilic CYP261 enzymes may be chromatographically separated into two fractions with different properties, and different thermodynamic parameters of spin equilibrium in particular. According to our concept, the changes in the energy landscape that evolved in pressure-tolerant enzymes must stabilize the less-hydrated, closed conformers, which may be transient in the catalytic mechanisms of nonpiezophilic enzymes. The studies of enzymes of piezophiles should help unravel the mechanisms that control water access during the catalytic cycle.

Potent mechanism-based inactivation of cytochrome P450 2B4 by 9-ethynylphenanthrene: implications for allosteric modulation of cytochrome P450 catalysis

The mechanism-based inactivation of cytochrome P450 2B4 (CYP2B4) by 9-ethynylphenanthrene (9EP) has been investigated. The partition ratio and k(inact) are 0.2 and 0.25 min(-1), respectively. Intriguingly, the inactivation exhibits sigmoidal kinetics with a Hill coefficient of 2.5 and an S(50) of 4.5 μM indicative of homotropic cooperativity. Enzyme inactivation led to an increase in mass of the apo-CYP2B4 by 218 Da as determined by electrospray ionization liquid chromatography and mass spectrometry, consistent with covalent protein modification. The modified CYP2B4 was purified to homogeneity and its structure determined by X-ray crystallography. The structure showed that 9EP is covalently attached to Oγ of Thr 302 via an ester bond, which is consistent with the increased mass of the protein. The presence of the bulky phenanthrenyl ring resulted in inward rotations of Phe 297 and Phe 206, leading to a compact active site. Thus, binding of another molecule of 9EP in the active site is prohibited. However, results from the quenching of 9EP fluorescence by unmodified or 9EP-modified CYP2B4 revealed at least two binding sites with distinct affinities, with the low-affinity site being the catalytic site and the high-affinity site on the protein periphery. Computer-aided docking and molecular dynamics simulations with one or two ligands bound revealed that the high-affinity site is situated at the entrance of a substrate access channel surrounded by the F' helix, β1-β2 loop, and β4 loop and functions as an allosteric site to enhance the efficiency of activation of the acetylenic group of 9EP and subsequent covalent modification of Thr 302.

X-ray crystal structure of the cytochrome P450 2B4 active site mutant F297A in complex with clopidogrel: insights into compensatory rearrangements of the binding pocket

Prior X-ray crystal structures of cytochrome P450 2B4 revealed the pivotal role of rearrangement of the side chains of residues F206 and F297 in the active site in accommodating various inhibitors or substrates. To explore the role of these residues, 2B4 F206A and F297A were created by site-directed mutagenesis and characterized functionally. The structure of F297A with clopidogrel demonstrated the reorientation of the ligand such that the methyl ester group is oriented toward the heme, whereas the thiophene moiety now extends to the additional void in the F297A mutant. Most interestingly, movement of the I helix and several amino acid side chains within the active site was observed in apparent response to the altered binding orientation. Results of flexible docking using the 2B4 wild type or the F297A-virtual mutant positioned either the thiophene or chlorophenyl group closer to heme. However, docking of clopidogrel using the real F297A mutant or a virtual mutant with the I-helix re-positioned oriented clopidogrel preferentially with either the methyl ester or the chlorophenyl group closest to heme. The study provides insight into how the altered active site adapts to accommodate and interact with the substrate in a distinct orientation while maintaining the overall closed protein conformation.

Stereoselective formation and metabolism of 4-hydroxy-retinoic Acid enantiomers by cytochrome p450 enzymes

All-trans-retinoic acid (atRA), the major active metabolite of vitamin A, plays a role in many biological processes, including maintenance of epithelia, immunity, and fertility and regulation of apoptosis and cell differentiation. atRA is metabolized mainly by CYP26A1, but other P450 enzymes such as CYP2C8 and CYP3As also contribute to atRA 4-hydroxylation. Although the primary metabolite of atRA, 4-OH-RA, possesses a chiral center, the stereochemical course of atRA 4-hydroxylation has not been studied previously. (4S)- and (4R)-OH-RA enantiomers were synthesized and separated by chiral column HPLC. CYP26A1 was found to form predominantly (4S)-OH-RA. This stereoselectivity was rationalized via docking of atRA in the active site of a CYP26A1 homology model. The docked structure showed a well defined niche for atRA within the active site and a specific orientation of the β-ionone ring above the plane of the heme consistent with stereoselective abstraction of the hydrogen atom from the pro-(S)-position. In contrast to CYP26A1, CYP3A4 formed the 4-OH-RA enantiomers in a 1:1 ratio and CYP3A5 preferentially formed (4R)-OH-RA. Interestingly, CYP3A7 and CYP2C8 preferentially formed (4S)-OH-RA from atRA. Both (4S)- and (4R)-OH-RA were substrates of CYP26A1 but (4S)-OH-RA was cleared 3-fold faster than (4R)-OH-RA. In addition, 4-oxo-RA was formed from (4R)-OH-RA but not from (4S)-OH-RA by CYP26A1. Overall, these findings show that (4S)-OH-RA is preferred over (4R)-OH-RA by the enzymes regulating atRA homeostasis. The stereoselectivity observed in CYP26A1 function will aid in better understanding of the active site features of the enzyme and the disposition of biologically active retinoids.

Conformational adaptation of human cytochrome P450 2B6 and rabbit cytochrome P450 2B4 revealed upon binding multiple amlodipine molecules

Structures of human cytochrome P450 2B6 and rabbit cytochrome P450 2B4 in complex with two molecules of the calcium channel blocker amlodipine have been determined by X-ray crystallography. The presence of two drug molecules suggests clear substrate access channels in each P450. According to a previously established nomenclature, amlodipine molecules were trapped in access pathway 2f in P450 2B6 and in pathway 2a or 2f in P450 2B4. These pathways overlap for part of the length and then diverge as they extend toward the protein surface. A previously described solvent channel was also found in each enzyme. The results indicate that key residues located on the surface and at the entrance of the substrate access channels in each of these P450s may play a crucial role in guiding substrate entry. In addition, the region of P450 2B6 and 2B4 involving helices B', F, F', and G' and part of helix G is substantially more open in the amlodipine complexes than in the corresponding 4-(4-chlorophenyl)imidazole complexes. The increased active site volume observed results from the major retraction of helices F, F', and B' and the β4 sheet region located close to the binding cavity to accommodate amlodipine. These structures demonstrate novel insight into distinct conformational states not observed with previous P450 2B structures and provide clear evidence of the substrate access channels in two drug-metabolizing P450s. In addition, the structures exhibit the versatility that can be exploited via in silico studies with other P450 2B6 ligands as large as raloxifene and itraconazole.

Plasticity of CYP2B enzymes: structural and solution biophysical methods

In the past three years, major advances in understanding cytochrome P450 2B (CYP2B) structure-function relationships have been made through determination of multiple ligand-bound and one ligand-free X-ray crystal structure of CYP2B4 and one ligand-bound X-ray crystal structure of CYP2B6. These structures have provided insight into the features that provide the high degree of plasticity of the enzymes. A combination of a phenylalanine cluster that allows for concerted movement of helices F through G and a conserved set of electrostatic interactions involving Arg(262) facilitates movement of this region to accommodate binding of ligands of various sizes without perturbing most of the P450 fold. Integrating solution based techniques such as NMR or deuterium exchange mass spectrometry (DXMS) with computational methods including molecular docking has provided further insight into enzyme behavior upon ligand binding. In addition, extended molecular dynamics simulations have provided a link between an open and a closed conformation of ligand-free CYP2B4 found in crystal structures. Other studies revealed the utility of rational engineering in improving stability of P450s to facilitate structural studies. The solution and computational results combined with the X-ray crystal structures yield a comprehensive picture of how these enzymes adopt different conformations to bind various ligands.

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

The mechanisms of ligand binding and allostery in the major human drug-metabolizing enzyme cytochrome P450 3A4 (CYP3A4) were explored with fluorescence resonance energy transfer (FRET) using a laser dye, fluorol-7GA (F7GA), as a model substrate. Incorporation into the enzyme of a thiol-reactive FRET probe, pyrene iodoacetamide, allowed us to monitor the binding by FRET from the pyrene donor to the F7GA acceptor. Cooperativity of the interactions detected by FRET indicates that the enzyme possesses at least two F7GA-binding sites that have different FRET efficiencies and are therefore widely separated. To probe spatial localization of these sites, we studied FRET in a series of mutants bearing pyrene iodoacetamide at different positions, and we measured the distances from each of the sites to the donor. Our results demonstrate the presence of a high affinity binding site at the enzyme periphery. Analysis of the set of measured distances complemented with molecular modeling and docking allowed us to pinpoint the most probable peripheral site. It is located in the vicinity of residues 217-220, similar to the position of the progesterone molecule bound at the distal surface of the CYP3A4 in a prior x-ray crystal structure. Peripheral binding of F7GA causes a substantial spin shift and serves as a prerequisite for the binding in the active site. This is the first indication of functionally important ligand binding outside of the active site in cytochromes P450. The findings strongly suggest that the mechanisms of CYP3A4 cooperativity involve a conformational transition triggered by an allosteric ligand.

Investigation by site-directed mutagenesis of the role of cytochrome P450 2B4 non-active-site residues in protein-ligand interactions based on crystal structures of the ligand-bound enzyme

Residues located outside the active site of cytochromes P450 2B have exhibited importance in ligand binding, structural stability and drug metabolism. However, contributions of non-active-site residues to the plasticity of these enzymes are not known. Thus, a systematic investigation was undertaken of unique residue-residue interactions found in crystal structures of P450 2B4 in complex with 4-(4-chlorophenyl)imidazole (4-CPI), a closed conformation, or in complex with bifonazole, an expanded conformation. Nineteen mutants distributed over 11 sites were constructed, expressed in Escherichia coli and purified. Most mutants showed significantly decreased expression, especially in the case of interactions found in the 4-CPI structure. Six mutants (H172A, H172F, H172Q, L437A, E474D and E474Q) were chosen for detailed functional analysis. Among these, the K(s) of H172F for bifonazole was ∼ 20 times higher than for wild-type 2B4, and the K(s) of L437A for 4-CPI was ∼ 50 times higher than for wild-type, leading to significantly altered inhibitor selectivity. Enzyme function was tested with the substrates 7-ethoxy-4-(trifluoromethyl)coumarin, 7-methoxy-4-(trifluoromethyl)coumarin and 7-benzyloxyresorufin (7-BR). H172F was inactive with all three substrates, and L437A did not turn over 7-BR. Furthermore, H172A, H172Q, E474D and E474Q showed large changes in k(cat)/K(M) for each of the three substrates, in some cases up to 50-fold. Concurrent molecular dynamics simulations yielded distances between some of the residues in these putative interaction pairs that are not consistent with contact. The results indicate that small changes in the protein scaffold lead to large differences in solution behavior and enzyme function.

The structural basis for homotropic and heterotropic cooperativity of midazolam metabolism by human cytochrome P450 3A4

Human cytochrome P450 3A4 (CYP3A4) metabolizes a significant portion of clinically relevant drugs and often exhibits complex steady-state kinetics that can involve homotropic and heterotropic cooperativity between bound ligands. In previous studies, the hydroxylation of the sedative midazolam (MDZ) exhibited homotropic cooperativity via a decrease in the ratio of 1'-OH-MDZ to 4-OH-MDZ at higher drug concentrations. In this study, MDZ exhibited heterotropic cooperativity with the antiepileptic drug carbamazepine (CBZ) with characteristic decreases in the 1'-OH-MDZ to 4-OH-MDZ ratios. To unravel the structural basis of MDZ cooperativity, we probed MDZ and CBZ bound to CYP3A4 using longitudinal T(1) nuclear magnetic resonance (NMR) relaxation and molecular docking with AutoDock 4.2. The distances calculated from longitudinal T(1) NMR relaxation were used during simulated annealing to constrain the molecules to the substrate-free X-ray crystal structure of CYP3A4. These simulations revealed that either two MDZ molecules or an MDZ molecule and a CBZ molecule assume a stacked configuration within the CYP3A4 active site. In either case, the proton at position 4 of the MDZ molecule was closer to the heme than the protons of the 1'-CH(3) group. In contrast, molecular docking of a single molecule of MDZ revealed that the molecule was preferentially oriented with the 1'-CH(3) position closer to the heme than position 4. This study provides the first detailed molecular analysis of heterotropic and homotropic cooperativity of a human cytochrome P450 from an NMR-based model. Cooperativity of ligand binding through direct interaction between stacked molecules may represent a common motif for homotropic and heterotropic cooperativity.

Structures of cytochrome P450 2B6 bound to 4-benzylpyridine and 4-(4-nitrobenzyl)pyridine: insight into inhibitor binding and rearrangement of active site side chains

The biochemical, biophysical, and structural analysis of the cytochrome P450 2B subfamily of enzymes has provided a wealth of information regarding conformational plasticity and substrate recognition. The recent X-ray crystal structure of the drug-metabolizing P450 2B6 in complex with 4-(4-chlorophenyl)imidazole (4-CPI) yielded the first atomic view of this human enzyme. However, knowledge of the structural basis of P450 2B6 specificity and inhibition has remained limited. In this study, structures of P450 2B6 were determined in complex with the potent inhibitors 4-benzylpyridine (4-BP) and 4-(4-nitrobenzyl)pyridine (4-NBP). Comparison of the present structures with the previous P450 2B6-4-CPI complex showed that reorientation of side chains of the active site residue Phe206 on the F-helix and Phe297 on the I-helix was necessary to accommodate the inhibitors. However, P450 2B6 does not require any major side chain rearrangement to bind 4-NBP compared with 4-BP, and the enzyme provides no hydrogen-bonding partners for the polar nitro group of 4-NBP within the hydrophobic active site. In addition, on the basis of these new structures, substitution of residue 172 with histidine as observed in the single nucleotide polymorphism Q172H and in P450 2B4 may contribute to a hydrogen bonding network connecting the E- and I-helices, thereby stabilizing active site residues on the I-helix. These results provide insight into the role of active site side chains upon inhibitor binding and indicate that the recognition of the benzylpyridines in the closed conformation structure of P450 2B6 is based solely on hydrophobicity, size, and shape.

Institutional Profile: University of California San Diego Pharmacogenomics Education Program (PharmGenEd™): bridging the gap between science and practice

Clinical application of evidence-based pharmacogenomics information has the potential to help healthcare professionals provide safe and effective medication management to patients. However, there is a gap between the advances of pharmacogenomics discovery and the health professionals' knowledge regarding pharmacogenomics testing and therapeutic uses. Furthermore, pharmacogenomics education materials for healthcare professionals have not been readily available or accessible. Pharmacogenomics Education Program (PharmGenEd™) is an evidence-based pharmacogenomics education program developed at the University of California San Diego Skaggs School of Pharmacy and Pharmaceutical Sciences and the School of Medicine (CA, USA), with funding support from the Centers for Disease Control and Prevention. Program components include continuing education modules, train-the-trainer materials and shared curriculum modules based on therapeutic topics, and virtual communities with online resources.

Structural analysis of mammalian cytochrome P450 2B4 covalently bound to the mechanism-based inactivator tert-butylphenylacetylene: insight into partial enzymatic activity

A combined structural and computational analysis of rabbit cytochrome P450 2B4 covalently bound to the mechanism-based inactivator tert-butylphenylacetylene (tBPA) has yielded insight into how the enzyme retains partial activity. Since conjugation to tBPA modifies a highly conserved active site residue, the residual activity of tBPA-labeled 2B4 observed in previous studies was puzzling. Here we describe the first crystal structures of a modified mammalian P450, which show an oxygenated metabolite of tBPA conjugated to Thr 302 of helix I. These results are consistent with previous studies that identified Thr 302 as the site of conjugation. In each structure, the core of 2B4 remains unchanged, but the arrangement of plastic regions differs. This results in one structure that is compact and closed. In this conformation, tBPA points toward helix B', making a 31° angle with the heme plane. This conformation is in agreement with previously performed in silico experiments. However, dimerization of 2B4 in the other structure, which is caused by movement of the B/C loop and helices F through G, alters the position of tBPA. In this case, tBPA lies almost parallel to the heme plane due to the presence of helix F' of the opposite monomer entering the active site to stabilize the dimer. However, docking experiments using this open form show that tBPA is able to rotate upward to give testosterone and 7-ethoxy-4-trifluoromethylcoumarin access to the heme, which could explain the previously observed partial activity.