Structure-function analysis of cytochromes P450 2B

Multidisciplinary Insights into the Structure-Function Relationship of the CYP2B6 Active Site

Cytochrome P450 2B6 (CYP2B6) is a highly polymorphic human enzyme involved in the metabolism of many clinically relevant drugs, environmental toxins, and endogenous molecules with disparate structures. Over the last 20-plus years, in silico and in vitro studies of CYP2B6 using various ligands have provided foundational information regarding the substrate specificity and structure-function relationship of this enzyme. Approaches such as homology modeling, X-ray crystallography, molecular docking, and kinetic activity assays coupled with CYP2B6 mutagenesis have done much to characterize this originally neglected monooxygenase. However, a complete understanding of the structural details that make new chemical entities substrates of CYP2B6 is still lacking. Surprisingly little in vitro data has been obtained about the structure-function relationship of amino acids identified to be in the CYP2B6 active site. Since much attention has already been devoted to elucidating the function of CYP2B6 allelic variants, here we review the salient findings of in silico and in vitro studies of the CYP2B6 structure-function relationship with a deliberate focus on the active site. In addition to summarizing these complementary approaches to studying structure-function relationships, we note gaps/challenges in existing data such as the need for more CYP2B6 crystal structures, molecular docking results with various ligands, and data coupling CYP2B6 active site mutagenesis with kinetic parameter measurement under standard expression conditions. Harnessing in silico and in vitro techniques in tandem to understand the CYP2B6 structure-function relationship will likely offer further insights into CYP2B6-mediated metabolism. SIGNIFICANCE STATEMENT: The apparent importance of cytochrome P450 2B6 (CYP2B6) in the metabolism of various xenobiotics and endogenous molecules has grown since its discovery with many in silico and in vitro studies offering a partial description of its structure-function relationship. Determining the structure-function relationship of CYP2B6 is difficult but may be aided by thorough biochemical investigations of the CYP2B6 active site that provide a more complete pharmacological understanding of this important enzyme.

Computer-Aided (In Silico) Modeling of Cytochrome P450-Mediated Food–Drug Interactions (FDI)

Modifications of the activity of Cytochrome 450 (CYP) enzymes by compounds in food might impair medical treatments. These CYP-mediated food–drug interactions (FDI) play a major role in drug clearance in the intestine and liver. Inter-individual variation in both CYP expression and structure is an important determinant of FDI. Traditional targeted approaches have highlighted a limited number of dietary inhibitors and single-nucleotide variations (SNVs), each determining personal CYP activity and inhibition. These approaches are costly in time, money and labor. Here, we review computational tools and databases that are already available and are relevant to predicting CYP-mediated FDIs. Computer-aided approaches such as protein–ligand interaction modeling and the virtual screening of big data narrow down hundreds of thousands of items in databanks to a few putative targets, to which the research resources could be further directed. Structure-based methods are used to explore the structural nature of the interaction between compounds and CYP enzymes. However, while collections of chemical, biochemical and genetic data are available today and call for the implementation of big-data approaches, ligand-based machine-learning approaches for virtual screening are still scarcely used for FDI studies. This review of CYP-mediated FDIs promises to attract scientists and the general public.

Dissecting the Structural Plasticity and Dynamics of Cytochrome P450 2B4 by Molecular Dynamics Simulations

The plasticity of cytochromes P450 (P450s) is known to contribute significantly to their catalytic capacity of metabolizing various substrates. Although numerous studies have been performed, factors governing the plasticity and dynamics of P450s are still not fully understood. In this study, taking CYP2B4 as an example, we dissect the protein plasticity and dynamics in different environments. CYP2B4 is featured by a high degree of plasticity, which exhibits open, closed, and intermediate states. By analyzing the CYP2B4 crystal structures, we identified the structural features for the closed, open, and intermediate states. Interestingly, formation of the dimer structure was found in the open and intermediate states. The subsequent molecular dynamics (MD) simulations of the open structure in water confirmed the importance of the dimer form in stabilizing the open conformations. MD simulations of the closed and open structures in the membrane environment and the free energies for opening the F-G cassette obtained from the umbrella sampling calculations indicate that the membrane environment is important for stabilizing the F-G cassette. The dynamical network analysis indicates that Asp105 on the B-C loop plays an important role in transiting the structure from the open to the intermediate state. Our results thus unveil the mechanisms of dimer formation and open-to-intermediate transition for CYP2B4 in the water and membrane environments.

Cytochrome P450 2B diversity in a dietary specialist—the red tree vole (Arborimus longicaudus)

Although herbivores rely on liver enzymes to biotransform plant secondary metabolites ingested in plant-based diets, only a few enzymes from a handful of species have been characterized at the genomic level. In this study, we examined cytochrome P450 2B (CYP2B) sequence diversity and gene copy number in a conifer specialist, the red tree vole (Arborimus longicaudus). We fed captive individuals exclusively Douglas-fir (Pseudotsuga menziesii) foliage, cloned and sequenced their liver CYP2B cDNA, and estimated CYP2B gene copy number. We identified 21 unique CYP2B nucleotide sequences, and 20 unique CYP2B amino acid sequences. Gene copy number of CYP2B was estimated at 7.7 copies per haploid genome. We compared red tree vole CYP2B with CYP2B sequences of a generalist, the prairie vole (Microtus ochrogaster), found in GenBank. Our study revealed that the CYP2B enzymes of red tree voles possess unique sequences compared to CYP2B enzymes of other herbivorous species. The unique combination of amino acid residues at key substrate recognition sites of CYP2B enzymes may underlie the ability of the red tree vole to specialize on a highly toxic diet of Douglas-fir.

Use of Phenoxyaniline Analogues To Generate Biochemical Insights into the Interactio n of Polybrominated Diphenyl Ether with CYP2B Enzymes

Human hepatic cytochromes P450 (CYP) are integral to xenobiotic metabolism. CYP2B6 is a major catalyst of biotransformation of environmental toxicants, including polybrominated diphenyl ethers (PBDEs). CYP2B substrates tend to contain halogen atoms, but the biochemical basis for this selectivity and for species specific determinants of metabolism has not been identified. Spectral binding titrations and inhibition studies were performed to investigate interactions of rat CYP2B1, rabbit CYP2B4, and CYP2B6 with a series of phenoxyaniline (POA) congeners that are analogues of PBDEs. For most congeners, there was a <3-fold difference between the spectral binding constants (K_S) and IC₅₀ values. In contrast, large discrepancies between these values were observed for POA and 3-chloro-4-phenoxyaniline. CYP2B1 was the enzyme most sensitive to POA congeners, so the Val-363 residue from that enzyme was introduced into CYP2B4 or CYP2B6. This substitution partially altered the protein-ligand interaction profiles to make them more similar to that of CYP2B1. Addition of cytochrome P450 oxidoreductase (POR) to titrations of CYP2B6 with POA or 2'4'5'TCPOA decreased the affinity of both ligands for the enzyme. Addition of cytochrome b₅ to a recombinant enzyme system containing POR and CYP2B6 increased the POA IC₅₀ value and decreased the 2'4'5'TCPOA IC₅₀ value. Overall, the inconsistency between K_S and IC₅₀ values for POA versus 2'4'5'TCPOA is largely due to the effects of redox partner binding. These results provide insight into the biochemical basis of binding of diphenyl ethers to human CYP2B6 and changes in CYP2B6-mediated metabolism that are dependent on POA congener and redox partner identity.

A CYP2B6-humanized mouse model and its potential applications

CYP2B6 is a human microsomal cytochrome P450 enzyme with broad substrate selectivity. CYP2B6 is the only functional member of the human CYP2B gene subfamily, which differs from the situation in rodents, such as mouse, where multiple functional Cyp2b genes are expressed. Recent studies with Cyp2b knockout or knockdown mouse models have yielded insights into the in vivo roles of mouse CYP2B enzymes in drug disposition and xenobiotic toxicity. A CYP2B6-humanized mouse model (CYP2A13/2B6/2F1-transgenic/Cyp2abfgs-null), which expresses human CYP2B6 in the liver, and human CYP2A13 and CYP2F1 in the respiratory tract, but not any of the mouse Cyp2b genes, has also been established. In the CYP2B6-humanized mouse, the CYP2B6 transgene is expressed primarily in the liver, where it was found to be active toward prototype CYP2B6 substrate drugs. The regulatory elements of the CYP2B6 transgene appear to be compatible with mouse nuclear receptors that mediate CYP2B induction. Therefore, the CYP2B6-humanized mouse is a valuable animal model for studying the impact of CYP2B6 expression or induction on drug metabolism, drug efficacy, drug-drug interaction, and drug/xenobiotic toxicity. In this mini-review, we provide a brief background on CYP2B6 and the Cyp2b-knockout and CYP2B6-humanized mice, and discuss the potential applications and limitations of the current models.

Roles of Residues F206 and V367 in Human CYP2B6: Effects of Mutations on Androgen Hydroxylation, Mechanism-Based Inactivation, and Reversible Inhibition

The crystal structures of human CYP2B6 indicate that Phe206 and Val367 are in close proximity to the substrate binding site and suggest that both residues may play important roles in substrate metabolism and inhibitor binding. To test this hypothesis, we investigated the effects of mutating these residues to Ala on the regiospecificity of CYP2B6 for the metabolism of testosterone and androstenedione. For testosterone metabolism, 16β-OH-testosterone formation by the F206A mutant was <5% of the wild type (WT), whereas the V367A mutant exhibited a doubling of 16α-OH-testosterone formation with a 50% decrease in 16β-OH-testosterone formation compared with the WT. Significant alterations in the regiospecificity for androstenedione metabolism were also observed. To investigate the roles of these two residues in the metabolic activation of mechanism-based inactivators, tert-butylphenylacetylene (BPA) and bergamottin (BG) were used to test the susceptibility to inactivation. Although the rates of inactivation of both mutants by BG were not significantly decreased compared with the WT, the efficiency of inactivation by BPA of both mutants was more than an order of magnitude lower. Our results demonstrate that Phe206 plays a crucial role in determining the specificity of CYP2B6 for the 16β-hydroxylation of testosterone and androstenedione and that it also plays an important role in BG binding and mechanism-based inactivation by BPA. In addition, Val367 dramatically enhances the catalytic activity of CYP2B6 toward androstenedione and plays an important role in mechanism-based inactivation by BPA. The results presented here show the important roles of Phe206 and Val367 in interactions of CYP2B6 with substrates and inactivators/inhibitors and are consistent with the crystal structures.

Effect of detergent binding on cytochrome P450 2B4 structure as analyzed by X-ray crystallography and deuterium-exchange mass spectrometry

Multiple crystal structures of CYP2B4 have demonstrated the binding of the detergent 5-cyclohexyl-1-pentyl-β-D-maltoside (CYMAL-5) in a peripheral pocket located adjacent to the active site. To explore the consequences of detergent binding, X-ray crystal structures of the peripheral pocket mutant CYP2B4 F202W were solved in the presence of hexaethylene glycol monooctyl ether (C8E6) and CYMAL-5. The structure in the presence of CYMAL-5 illustrated a closed conformation indistinguishable from the previously solved wild-type. In contrast, the F202W structure in the presence of C8E6 revealed a detergent molecule that coordinated the heme-iron and extended to the protein surface through the substrate access channel 2f. Despite the overall structural similarity of these detergent complexes, remarkable differences were observed in the A, A', and H helices, the F-G cassette, the C-D and β4 loop region. Hydrogen-deuterium exchange mass spectrometry (DXMS) was employed to probe these differences and to test the effect of detergents in solution. The presence of either detergent increased the H/D exchange rate across the plastic regions, and the results obtained by DXMS in solution were consistent in general with the relevant structural snapshots. The study provides insight into effect of detergent binding and the interpretation of associated conformational dynamics of CYP2B4.

Coumarin Derivatives as Substrate Probes of Mammalian Cytochromes P450 2B4 and 2B6: Assessing the Importance of 7-Alkoxy Chain Length, Halogen Substitution, and Non-Active Site Mutations

Using a combined structural and biochemical approach, the functional importance of a recently described peripheral pocket bounded by the E-, F-, G-, and I-helices in CYP2B4 and 2B6 was probed. Three series of 4-substituted-7-alkoxycoumarin derivatives with -H, -CH3, or -CF3 at the 4 position of the coumarin core were used initially to monitor functional differences between CYP2B4 and 2B6. 7-Ethoxy-4-(trifluoromethyl)coumarin (7-EFC) displayed the highest catalytic efficiency among these substrates. Mutants were made to alter side-chain polarity (V/E194Q) or bulk (F/Y244W) to alter access to the peripheral pocket. Modest increases in catalytic efficiency of 7-EFC O-deethylation by the mutants were magnified considerably by chlorination or bromination of the substrate ethoxy chain. A structure of CYP2B6 Y244W in complex with (+)-α-pinene was solved at 2.2 Å and showed no CYMAL-5 in the peripheral pocket. A ligand free structure of CYP2B4 F244W was solved at 3.0 Å with CYMAL-5 in the peripheral pocket. In both instances, comparison of the respective wild-type and mutant CYP2B enzymes revealed that CYMAL-5 occupancy of the peripheral pocket had little effect on the topology of active site residue side-chains, despite the fact that the peripheral pocket and active site are located on opposite sides of the I-helix. Analysis of available CYP2B structures suggest that the effect of the amino acid substitutions within the peripheral pocket derive from altered interactions between the F and G helices.

Structure-Function Analysis of Mammalian CYP2B Enzymes Using 7-Substituted Coumarin Derivatives as Probes: Utility of Crystal Structures and Molecular Modeling in Understanding Xenobiotic Metabolism

Crystal structures of CYP2B35 and CYP2B37 from the desert woodrat were solved in complex with 4-(4-chlorophenyl)imidazole (4-CPI). The closed conformation of CYP2B35 contained two molecules of 4-CPI within the active site, whereas the CYP2B37 structure demonstrated an open conformation with three 4-CPI molecules, one within the active site and the other two in the substrate access channel. To probe structure-function relationships of CYP2B35, CYP2B37, and the related CYP2B36, we tested the O-dealkylation of three series of related substrates-namely, 7-alkoxycoumarins, 7-alkoxy-4-(trifluoromethyl)coumarins, and 7-alkoxy-4-methylcoumarins-with a C1-C7 side chain. CYP2B35 showed the highest catalytic efficiency (kcat/KM) with 7-heptoxycoumarin as a substrate, followed by 7-hexoxycoumarin. In contrast, CYP2B37 showed the highest catalytic efficiency with 7-ethoxy-4-(trifluoromethyl)coumarin (7-EFC), followed by 7-methoxy-4-(trifluoromethyl)coumarin (7-MFC). CYP2B35 had no dealkylation activity with 7-MFC or 7-EFC. Furthermore, the new CYP2B-4-CPI-bound structures were used as templates for docking the 7-substituted coumarin derivatives, which revealed orientations consistent with the functional studies. In addition, the observation of multiple -Cl and -NH-π interactions of 4-CPI with the aromatic side chains in the CYP2B35 and CYP2B37 structures provides insight into the influence of such functional groups on CYP2B ligand binding affinity and specificity. To conclude, structural, computational, and functional analysis revealed striking differences between the active sites of CYP2B35 and CYP2B37 that will aid in the elucidation of new structure-activity relationships.

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.

Plasticity of specialized metabolism as mediated by dynamic metabolons

The formation of specialized metabolites enables plants to respond to biotic and abiotic stresses, but requires the sequential action of multiple enzymes. To facilitate swift production and to avoid leakage of potentially toxic and labile intermediates, many of the biosynthetic pathways are thought to organize in multienzyme clusters termed metabolons. Dynamic assembly and disassembly enable the plant to rapidly switch the product profile and thereby prioritize its resources. The lifetime of metabolons is largely unknown mainly due to technological limitations. This review focuses on the factors that facilitate and stimulate the dynamic assembly of metabolons, including microenvironments, noncatalytic proteins, and allosteric regulation. Understanding how plants organize carbon fluxes within their metabolic grids would enable targeted bioengineering of high-value specialized metabolites.

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.

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.

Novel CYP2B6 enzyme variants in a Rwandese population: functional characterization and assessment of in silico prediction tools

Cytochrome P450 CYP2B6 is a highly polymorphic enzyme that metabolizes numerous drugs, pesticides, and environmental toxins. Sequence analysis of a Rwandese population identified eight functionally uncharacterized nonsynonymous variants c.329G>T (p.G110V), c.341T>C (p.I114T), c.444G>T (p.E148D), c.548T>G (p.V183G), c.637T>C (p.F213L), c.758G>A (p.R253H), c.835G>C (p.A279P), and c.1459C>A (p.R487S), and five novel alleles termed CYP2B6*33 to CYP2B6*37 were assigned. Recombinant expression in COS-1 cells and functional characterization using the antidepressant bupropion and the antiretroviral efavirenz (EFV) as substrates demonstrated complete or almost complete loss-of-function for variants p.G110V, p.I114T, p.V183G, and p.F213L, whereas p.E148D, p.R253H, p.A279P, and p.R487S variants were functional. The data were used to assess the predictive power of eight online available functional prediction programs for amino-acid changes. Although none of the programs correctly predicted the functionality of all variants, substrate docking simulation analyses indicated similar conformational changes by all four deleterious mutations within the enzyme's active site, thus explaining lack of enzymatic function for both substrates. Because low-activity alleles of CYP2B6 are associated with impaired EFV metabolism and adverse drug response, these results are of potential utility for personalized treatment strategies in HIV/AIDS therapy.

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.

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.

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.

Thr302 is the site for the covalent modification of human cytochrome P450 2B6 leading to mechanism-based inactivation by tert-butylphenylacetylene

The mechanism-based inactivation of human CYP2B6 by tert-butylphenylacetylene (BPA) in the reconstituted system was investigated. The inactivation of CYP2B6 by BPA is time-, concentration-, and NADPH-dependent and exhibits a K(I) of 2.8 μM, a k(inact) of 0.7 min(-1), and a t(1/2) of 1 min. The partition ratio is ∼5. Unlike CYP2B1 and CYP2B4, in addition to the formation of an apoprotein adduct and a glutathione conjugate, a small heme adduct was observed when CYP2B6 was incubated with BPA. The mass increase of the adducted apoprotein and GSH conjugate is 174 Da, equivalent to the mass of one molecule of BPA plus one oxygen atom. To identify the adducted residue, BPA-inactivated CYP2B6 was digested with trypsin, and the digest was then analyzed by liquid chromatography-tandem mass spectrometry. A mass shift of 174 Da was used for the SEQUEST database search, and the identity of the modified residue was confirmed by MS/MS fragmentation of the modified peptide. Two residues, Lys274 and Thr302, were identified as having been modified. Further mutagenesis studies have demonstrated that the residue that is modified to result in inactivation is Thr302, not Lys274. Docking studies show that in the enzyme-substrate complex, Thr302 is in close contact with the triple bond of BPA with a distance of 3.8 Å between the terminal carbon of BPA and the oxygen in the hydroxyl group of Thr302. In conclusion, Thr302 of CYP2B6 is covalently modified by a reactive metabolite of BPA, and this modification is responsible for the mechanism-based inactivation.

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

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

Analysis of Cytochrome P450 Conserved Sequence Motifs between Helices E and H: Prediction of Critical Motifs and Residues in Enzyme Functions

Rational approaches have been extensively used to investigate the role of active site residues in cytochrome P450 (CYP) functions. However, recent studies using random mutagenesis suggest an important role for non-active site residues in CYP functions. Meta-analysis of the random mutants showed that 75% of the functionally important non-active site residues are present in 20% of the entire protein between helices E and H (E-H) and conserved sequence motif (CSM) between 7 and 11. The CSM approach was developed recently to investigate the functional role of non-active site residues in CYP2B4. Furthermore, we identified and analyzed the CSM in multiple CYP families and subfamilies in the E-H region. Results from CSM analysis showed that CSM 7, 8, 10, and 11 are conserved in CYP1, CYP2, and CYP3 families, while CSM 9 is conserved only in CYP2 family. Analysis of different CYP2 subfamilies showed that CYP2B and CYP2C have similar characteristics in the CSM, while the characteristics of CYP2A and CYP2D subfamilies are different. Finally, we analyzed CSM 7, 8, 10, and 11, which are common in all the CYP families/subfamilies analyzed, in fifteen important drug-metabolizing CYPs. The results showed that while CSM 8 is most conserved among these CYPs, CSM 7, 9, and 10 have significant variations. We suggest that CSM8 has a common role in all the CYPs that have been analyzed, while CSM 7, 10, and 11 may have relatively specific role within the subfamily. We further suggest that these CSM play important role in opening and closing of the substrate access/egress channel by modulating the flexible/plastic region of the protein. Thus, site-directed mutagenesis of these CSM can be used to study structure-function and dynamic/plasticity-function relationships and to design CYP biocatalysts.

The TonB dimeric crystal structures do not exist in vivo

The TonB system energizes transport of nutrients across the outer membrane of Escherichia coli using cytoplasmic membrane proton motive force (PMF) for energy. Integral cytoplasmic membrane proteins ExbB and ExbD appear to harvest PMF and transduce it to TonB. The carboxy terminus of TonB then physically interacts with outer membrane transporters to allow translocation of ligands into the periplasmic space. The structure of the TonB carboxy terminus (residues ~150 to 239) has been solved several times with similar results. Our previous results hinted that in vitro structures might not mimic the dimeric conformations that characterize TonB in vivo. To test structural predictions and to identify irreplaceable residues, the entire carboxy terminus of TonB was scanned with Cys substitutions. TonB I232C and N233C, predicted to efficiently form disulfide-linked dimers in the crystal structures, did not do so. In contrast, Cys substitutions positioned at large distances from one another in the crystal structures efficiently formed dimers. Cys scanning identified seven functionally important residues. However, no single residue was irreplaceable. The phenotypes conferred by changes of the seven residues depended on both the specific assay used and the residue substituted. All seven residues were synergistic with one another. The buried nature of the residues in the structures was also inconsistent with these properties. Taken together, these results indicate that the solved dimeric crystal structures of TonB do not exist. The most likely explanation for the aberrant structures is that they were obtained in the absence of the TonB transmembrane domain, ExbB, ExbD, and/or the PMF.

The TonB system of Gram-negative bacteria is an attractive target for development of novel antibiotics because of its importance in iron acquisition and virulence. Logically, therefore, the structure of TonB must be accurately understood. TonB functions as a dimer in vivo, and two different but similar crystal structures of the dimeric carboxy-terminal ~90 amino acids gave rise to mechanistic models. Here we demonstrate that the crystal structures, and therefore the models based on them, are not biologically relevant. The idiosyncratic phenotypes conferred by substitutions at the only seven functionally important residues in the carboxy terminus suggest that similar to interaction of cytochromes P450 with numerous substrates, these residues allow TonB to differentially interact with different outer membrane transporters. Taken together, data suggest that TonB is maintained poised between order and disorder by ExbB, ExbD, and the proton motive force (PMF) before energy transduction to the outer membrane transporters.

Plasticity of cytochrome P450 2B4 as investigated by hydrogen-deuterium exchange mass spectrometry and X-ray crystallography

Crystal structures of the xenobiotic metabolizing cytochrome P450 2B4 have demonstrated markedly different conformations in the presence of imidazole inhibitors or in the absence of ligand. However, knowledge of the plasticity of the enzyme in solution has remained scant. Thus, hydrogen-deuterium exchange mass spectrometry (DXMS) was utilized to probe the conformations of ligand-free P450 2B4 and the complex with 4-(4-chlorophenyl)imidazole (4-CPI) or 1-biphenyl-4-methyl-1H-imidazole (1-PBI). The results of DXMS indicate that the binding of 4-CPI slowed the hydrogen-deuterium exchange rate over the B'- and C-helices and portions of the F-G-helix cassette compared with P450 2B4 in the absence of ligands. In contrast, there was little difference between the ligand-free and 1-PBI-bound exchange sets. In addition, DXMS suggests that the ligand-free P450 2B4 is predominantly open in solution. Interestingly, a new high resolution structure of ligand-free P450 2B4 was obtained in a closed conformation very similar to the 4-CPI complex. Molecular dynamics simulations performed with the closed ligand-free structure as the starting point were used to probe the energetically accessible conformations of P450 2B4. The simulations were found to equilibrate to a conformation resembling the 1-PBI-bound P450 2B4 crystal structure. The results indicate that conformational changes observed in available crystal structures of the promiscuous xenobiotic metabolizing cytochrome P450 2B4 are consistent with its solution structural behavior.

Structures of cytochrome P450 2B4 complexed with the antiplatelet drugs ticlopidine and clopidogrel

Prior X-ray crystal structures of rabbit cytochrome P450 2B4 (2B4) in complexes with various imidazoles have demonstrated markedly different enzyme conformations depending on the size of the inhibitor occupying the active site. In this study, structures of 2B4 were determined with the antiplatelet drugs clopidogrel and ticlopidine, which were expected to have greater freedom of movement in the binding pocket. Ticlopidine could be modeled into the electron density maps in two distinct orientations, both of which are consistent with metabolic data gathered with other mammalian P450 enzymes. Results of ligand docking and heme-induced NMR relaxation of drug protons showed that ticlopidine was preferentially oriented with the chlorophenyl group closest to the heme. Because of its stereocenter, clopidogrel was easier to fit in the electron density and exhibited a single orientation, which points the chlorophenyl ring toward the heme. The C(α) traces of both complexes aligned very well with each other and revealed a compact, closed structure that resembles the conformation observed in two previously determined 2B4 structures with the small molecule inhibitors 4-(4-chlorophenyl)imidazole and 1-(4-chlorophenyl)imidazole. The 2B4 active site is able to accommodate small ligands by moving only a small number of side chains, suggesting that ligand reorientation is energetically favored over protein conformational changes for binding of these similarly sized molecules. Adjusting both protein conformation and ligand orientation in the active site gives 2B4 the flexibility to bind to the widest range of molecules, while also being energetically favorable.

Engineering cytochrome P450 biocatalysts for biotechnology, medicine and bioremediation

IMPORTANCE OF THE FIELD

Cytochrome P450 enzymes comprise a superfamily of heme monooxygenases that are of considerable interest for the: i) synthesis of novel drugs and drug metabolites; ii) targeted cancer gene therapy; iii) biosensor design; and iv) bioremediation. However, their applications are limited because cytochrome P450, especially mammalian P450 enzymes, show a low turnover rate and stability, and require a complex source of electrons through cytochrome P450 reductase and NADPH.

AREAS COVERED IN THIS REVIEW

In this review, we discuss the recent progress towards the use of P450 enzymes in a variety of the above-mentioned applications. We also present alternate and cost-effective ways to perform P450-mediated reaction, especially using peroxides. Furthermore, we expand upon the current progress in P450 engineering approaches describing several recent examples that are utilized to enhance heterologous expression, stability, catalytic efficiency and utilization of alternate oxidants.

WHAT THE READER WILL GAIN

The review provides a comprehensive knowledge in the design of P450 biocatalysts for potentially practical purposes. Finally, we provide a prospective on the future aspects of P450 engineering and its applications in biotechnology, medicine and bioremediation.

TAKE HOME MESSAGE

Because of its wide applications, academic and pharmaceutical researchers, environmental scientists and healthcare providers are expected to gain current knowledge and future prospects of the practical use of P450 biocatalysts.

Covalent modification of Thr302 in cytochrome P450 2B1 by the mechanism-based inactivator 4-tert-butylphenylacetylene

The mechanism of inactivation of cytochrome P450 2B1 (CYP2B1) by 4-tert-butylphenylacetylene (BPA) has been characterized previously to be caused by the covalent binding of a reactive intermediate to the apoprotein rather than heme destruction (J Pharmacol Exp Ther 331:392-403, 2009). The identification of a BPA-glutathione conjugate and the increase in the mass of the BPA-adducted apoprotein have indicated that the mass of adduct is 174 Da, equivalent to the mass of BPA plus one oxygen atom. To identify the adducted residue, BPA-inactivated CYP2B1 was digested with trypsin, and the digest was then analyzed by using capillary liquid chromatography with a LTQ linear ion trap mass spectrometer as the detector. A mass shift of 174 Da was used for a SEQUEST database search. The tandem mass spectrometry fragmentation of the modified peptide and the identity of modified residue were determined. The results revealed a mass increase of 174 Da for the peptide sequence (296)FFAGTSSTTLR(308) in the I-helix of CYP2B1 and that the site of adduction formation is Thr302. Homology modeling and ligand docking studies showed that BPA binds in close proximity to both the heme iron and Thr302 with the distances being 2.96 and 3.42 A, respectively. The identification of Thr302 in the CYP2B1 active site as the site of covalent modification leading to inactivation by BPA supports previous hypotheses that this conserved Thr residue may play a crucial role for various functions in P450s.

Crystal structure of a cytochrome P450 2B6 genetic variant in complex with the inhibitor 4-(4-chlorophenyl)imidazole at 2.0-A resolution

The structure of the K262R genetic variant of human cytochrome P450 2B6 in complex with the inhibitor 4-(4-chlorophenyl)imidazole (4-CPI) has been determined using X-ray crystallography to 2.0-A resolution. Production of diffraction quality crystals was enabled through a combination of protein engineering, chaperone coexpression, modifications to the purification protocol, and the use of unique facial amphiphiles during crystallization. The 2B6-4-CPI complex is virtually identical to the rabbit 2B4 structure bound to the same inhibitor with respect to the arrangement of secondary structural elements and the placement of active site residues. The structure supports prior P450 2B6 homology models based on other mammalian cytochromes P450 and is consistent with the limited site-directed mutagenesis studies on 2B6 and extensive studies on P450 2B4 and 2B1. Although the K262R genetic variant shows unaltered binding of 4-CPI, altered binding affinity, kinetics, and/or product profiles have been previously shown with several other ligands. On the basis of new P450 2B6 crystal structure and previous 2B4 structures, substitutions at residue 262 affect a hydrogen-bonding network connecting the G and H helices, where subtle differences could be transduced to the active site. Docking experiments indicate that the closed protein conformation allows smaller ligands such as ticlopidine to bind to the 2B6 active site in the expected orientation. However, it is unknown whether 2B6 undergoes structural reorganization to accommodate bulkier molecules, as previously inferred from multiple P450 2B4 crystal structures.

Rational engineering of cytochromes P450 2B6 and 2B11 for enhanced stability: Insights into structural importance of residue 334

Rational mutagenesis was used to improve the thermal stability of human cytochrome P450 2B6 and canine P450 2B11. Comparison of the amino acid sequences revealed seven sites that are conserved between the stable 2B1 and 2B4 but different from those found in the less stable 2B6 and 2B11. P334S was the only mutant that showed increased heterologous expression levels and thermal stability in both 2B6 and 2B11. The mechanism of this effect was explored with pressure-perturbation spectroscopy. Compressibility of the heme pocket in variants of all four CYP2B enzymes containing proline at position 334 are characterized by lower compressibility than their more stable serine 334 counterpart. Therefore, the stabilizing effect of P334S is associated with increased conformational flexibility in the region of the heme pocket. Improved stability of P334S 2B6 and 2B11 may facilitate the studies of these enzymes by X-ray crystallography and biophysical techniques.

Mechanism-based inactivation of CYP2B1 and its F-helix mutant by two tert-butyl acetylenic compounds: covalent modification of prosthetic heme versus apoprotein

The mechanism-based inactivation of cytochrome CYP2B1 [wild type (WT)] and its Thr205 to Ala mutant (T205A) by tert-butylphenylacetylene (BPA) and tert-butyl 1-methyl-2-propynyl ether (BMP) in the reconstituted system was investigated. The inactivation of WT by BPA exhibited a k(inact)/K(I) value of 1343 min(-1)mM(-1) and a partition ratio of 1. The inactivation of WT by BMP exhibited a k(inact)/K(I) value of 33 min(-1)mM(-1) and a partition ratio of 10. Liquid chromatography/tandem mass spectrometry analysis (LC/MS/MS) of the WT revealed 1) inactivation by BPA resulted in the formation of a protein adduct with a mass increase equivalent to the mass of BPA plus one oxygen atom, and 2) inactivation by BMP resulted in the formation of multiple heme adducts that all exhibited a mass increase equivalent to BMP plus one oxygen atom. LC/MS/MS analysis indicated the formation of glutathione (GSH) conjugates by the reaction of GSH with the ethynyl moiety of BMP or BPA with the oxygen being added to the internal or terminal carbon. For the inactivation of T205A by BPA and BMP, the k(inact)/K(I) values were suppressed by 100- and 4-fold, respectively, and the partition ratios were increased 9- and 3.5-fold, respectively. Only one major heme adduct was detected following the inactivation of the T205A by BMP. These results show that the Thr205 in the F-helix plays an important role in the efficiency of the mechanism-based inactivation of CYP2B1 by BPA and BMP. Homology modeling and substrate docking studies were presented to facilitate the interpretation of the experimental results.

Crystal structures of cytochrome P450 2B4 in complex with the inhibitor 1-biphenyl-4-methyl-1H-imidazole: ligand-induced structural response through alpha-helical repositioning

Two different ligand occupancy structures of cytochrome P450 2B4 (CYP2B4) in complex with 1-biphenyl-4-methyl-1H-imidazole (1-PBI) have been determined by X-ray crystallography. 1-PBI belongs to a series of tight binding, imidazole-based CYP2B4 inhibitors. 1-PBI binding to CYP2B4 yields a type II spectrum with a K(s) value of 0.23 microM and inhibits enzyme activity with an IC(50) value of 0.035 microM. Previous CYP2B4 structures have shown a large degree of structural movement in response to ligand size. With two phenyl rings, 1-PBI is larger than 1-(4-chlorophenyl)imidazole (1-CPI) and 4-(4-chlorophenyl)imidazole (4-CPI) but smaller than bifonazole, which is branched and contains three phenyl rings. The CYP2B4-1-PBI complex is a structural intermediate to the closed CPI and the open bifonazole structures. The B/C-loop reorganizes itself to include two short partial helices while closing one side of the active site. The F-G-helix cassette pivots over the I-helix in direct response to the size of the ligand in the active site. A cluster of Phe residues at the fulcrum of this pivot point allows for dramatic repositioning of the cassette with only a relatively small amount of secondary structure rearrangement. Comparisons of ligand-bound CYP2B4 structures reveal trends in plastic region mobility that could allow for predictions of their position in future structures based on ligand shape and size.

Cooperative properties of cytochromes P450

Cytochromes P450 form a large and important class of heme monooxygenases with a broad spectrum of substrates and corresponding functions, from steroid hormone biosynthesis to the metabolism of xenobiotics. Despite decades of study, the molecular mechanisms responsible for the complex non-Michaelis behavior observed with many members of this superfamily during metabolism, often termed 'cooperativity', remain to be fully elucidated. Although there is evidence that oligomerization may play an important role in defining the observed cooperativity, some monomeric cytochromes P450, particularly those involved in xenobiotic metabolism, also display this behavior due to their ability to simultaneously bind several substrate molecules. As a result, formation of distinct enzyme-substrate complexes with different stoichiometry and functional properties can give rise to homotropic and heterotropic cooperative behavior. This review aims to summarize the current understanding of cooperativity in cytochromes P450, with a focus on the nature of cooperative effects in monomeric enzymes.

Allosteric P450 mechanisms: multiple binding sites, multiple conformers or both?

According to the initial hypothesis on the mechanisms of cooperativity in drug-metabolizing cytochromes P450, a loose fit of a single substrate molecule in the P450 active site results in a requirement for the binding of multiple ligand molecules for efficient catalysis. Although simultaneous occupancy of the active site by multiple ligands is now well established, there is increasing evidence that the mechanistic basis of cooperativity also involves an important ligand-induced conformational transition. Moreover, recent studies demonstrate that the conformational heterogeneity of the enzyme is stabilized by ligand-dependent interactions of several P450 molecules. Application of the concept of an oligomeric allosteric enzyme to microsomal cytochromes P450 in combination with a general paradigm of multiple ligand occupancy of the active site provides an excellent explanation for complex manifestations of the atypical kinetic behavior of the enzyme.

A survey of the year 2007 literature on applications of isothermal titration calorimetry

Elucidation of the energetic principles of binding affinity and specificity is a central task in many branches of current sciences: biology, medicine, pharmacology, chemistry, material sciences, etc. In biomedical research, integral approaches combining structural information with in-solution biophysical data have proved to be a powerful way toward understanding the physical basis of vital cellular phenomena. Isothermal titration calorimetry (ITC) is a valuable experimental tool facilitating quantification of the thermodynamic parameters that characterize recognition processes involving biomacromolecules. The method provides access to all relevant thermodynamic information by performing a few experiments. In particular, ITC experiments allow to by-pass tedious and (rarely precise) procedures aimed at determining the changes in enthalpy and entropy upon binding by van't Hoff analysis. Notwithstanding limitations, ITC has now the reputation of being the "gold standard" and ITC data are widely used to validate theoretical predictions of thermodynamic parameters, as well as to benchmark the results of novel binding assays. In this paper, we discuss several publications from 2007 reporting ITC results. The focus is on applications in biologically oriented fields. We do not intend a comprehensive coverage of all newly accumulated information. Rather, we emphasize work which has captured our attention with originality and far-reaching analysis, or else has provided ideas for expanding the potential of the method.

Allosteric transitions in cytochrome P450eryF explored with pressure-perturbation spectroscopy, lifetime FRET, and a novel fluorescent substrate, Fluorol-7GA

To establish a direct method for monitoring substrate binding in cytochrome P450eryF applicable at elevated hydrostatic pressures, we introduce a laser dye Fluorol-7GA (F7GA) as a novel fluorescent ligand. The high intensity of fluorescence and the reasonable resolution of the excitation band from the absorbance bands of P450 allowed us to establish highly sensitive binding assays compatible with pressure perturbation. The interactions of F7GA with P450eryF cause an ample spin shift revealing cooperative binding ( S50 = 8.2 +/- 1.3 microM; n = 2.3 +/- 0.1). Fluorescence resonance energy transfer (FRET) experiments suggest the presence of at least two substrate binding sites with apparent K D values in the ranges of 0.1-0.3 and 6-9 microM. Similar to that observed earlier with CYP3A4, increasing hydrostatic pressure does not cause either a complete dissociation of the substrate complexes or a displacement of the spin equilibrium toward the low-spin state. Rather, increased pressure enhances the cooperativity of the F7GA-induced spin shift, so that the Hill coefficient approaches 3 at 2 kbar. Lifetime FRET experiments revealed an important increase in the affinity of the enzyme for F7GA at elevated pressures, suggesting that the binding of the ligand induces a conformational transition associated with an important increase in the level of protein hydration. This transition largely attenuates the solvent accessibility of the heme pocket and causes an unusual stability of the high-spin, substrate-bound enzyme at elevated pressures.

Modification of serine 360 by a reactive intermediate of 17-alpha-ethynylestradiol results in mechanism-based inactivation of cytochrome P450s 2B1 and 2B6

17-alpha-Ethynylestradiol (17EE) is a mechanism-based inactivator of P450 2B1 and P450 2B6 in the reconstituted monooxygenase system. The loss in enzymatic activity was due to the binding of a reactive intermediate of 17EE to the apoprotein. P450 2B1 and P450 2B6 were inactivated by 17EE and digested with trypsin. The peptides obtained following digestion with trypsin of 17EE-inactivated P450 2B1 and P450 2B6 were separated by liquid chromatography and analyzed by ESI-MS. Adducted peptides exhibiting an increase in mass consistent with the addition of the mass of the reactive intermediate of 17EE were identified for each enzyme. Analysis of these modified peptides by ESI-MS/MS and precursor ion scanning facilitated the identification of the Ser360 in both enzymes as a site that had been adducted by a reactive intermediate of 17EE. A P450 2B1 mutant where Ser360 was replaced by alanine was constructed, expressed, and purified. Activity and inactivation studies indicated that mutation of the Ser360 residue to alanine did not prevent inactivation of the mutant enzyme by 17EE. These observations suggest that Ser360 is not critical for the catalytic function of these P450s. Spectral binding studies of the 17EE-inactivated P450 2B1 and P450 2B6 indicated that modification of the enzymes by the reactive intermediate of 17EE resulted in an enzyme that was no longer capable of binding substrates. These results suggest that the inactivation by 17EE may be due to modification of an amino acid residue in the substrate access channel near the point of entry into the active site.

Characterization of the microsomal cytochrome P450 2B4 O2 activation intermediates by cryoreduction and electron paramagnetic resonance

The oxy-ferrous complex of cytochrome P450 2B4 (2B4) has been prepared at -40 degrees C with and without bound substrate [butylated hydroxytoluene (BHT)] and radiolytically one-electron cryoreduced at 77 K. Electron paramagnetic resonance (EPR) shows that in both cases the observed product of cryoreduction is the hydroperoxo-ferriheme species, indicating that the microsomal P450 contains an efficient distal-pocket proton-delivery network. In the absence of substrate, two distinct hydroperoxo-ferriheme signals are observed, reflecting the presence of two major conformational substates in the oxy-ferrous precursor. Only one species is observed when BHT is bound, indicating a more ordered active site. BHT binding also changes the g-tensor components of the hydroperoxo-ferric 2B4 intermediate, indicating that the substrate modulates the properties of this intermediate. Step annealing the cryoreduced ternary 2B4 complex at >or=175 K causes the loss of hydroperoxo-ferric 2B4 and the parallel appearance of high-spin ferric 2B4; liquid chromatography-tandem mass spectroscopy (LC-MS/MS) analysis shows that in this process BHT is quantitatively converted to two products, hydroxymethyl BHT (1) and 3-hydroxy- tert-butyl BHT (2). This implies that the hydroperoxo-ferric 2B4 prepared by cryoreduction is catalytically active and that the high-spin state observed after annealing contains an enzyme-bound product of BHT monooxygenation. The ratio of products generated during cryoreduction and annealing (6.2/1) is significantly different from the ratio (2.5/1) at ambient temperature. These findings suggest that substrate is held more rigidly relative to the oxidizing species at low temperatures and/or that dissociation of FeOOH is inhibited at low temperature. As in experiments under ambient conditions, product formation is not observed with the inactive F429H 2B4 mutant.

Identification and analysis of conserved sequence motifs in cytochrome P450 family 2. Functional and structural role of a motif 187RFDYKD192 in CYP2B enzymes

Using a multiple alignment of 175 cytochrome P450 (CYP) family 2 sequences, 20 conserved sequence motifs (CSMs) were identified with the program PCPMer. Functional importance of the CSM in CYP2B enzymes was assessed from available data on site-directed mutants and genetic variants. These analyses suggested an important role of the CSM 8, which corresponds to(187)RFDYKD(192) in CYP2B4. Further analysis showed that residues 187, 188, 190, and 192 have a very high rank order of conservation compared with 189 and 191. Therefore, eight mutants (R187A, R187K, F188A, D189A, Y190A, K191A, D192A, and a negative control K186A) were made in an N-terminal truncated and modified form of CYP2B4 with an internal mutation, which is termed 2B4dH/H226Y. Function was examined with the substrates 7-methoxy-4-(trifluoromethyl)coumarin (7-MFC), 7-ethoxy-4-(trifluoromethyl)coumarin (7-EFC), 7-benzyloxy-4-(trifluoromethyl)coumarin (7-BFC), and testosterone and with the inhibitors 4-(4-chlorophenyl)imidazole (4-CPI) and bifonazole (BIF). Compared with the template and K186A, the mutants R187A, R187K, F188A, Y190A, and D192A showed > or =2-fold altered substrate specificity, k(cat), K(m), and/or k(cat)/K(m) for 7-MFC and 7-EFC and 3- to 6-fold decreases in differential inhibition (IC(50,BIF)/IC(50,4-CPI)). Subsequently, these mutants displayed 5-12 degrees C decreases in thermal stability (T(m)) and 2-8 degrees C decreases in catalytic tolerance to temperature (T(50)) compared with the template and K186A. Furthermore, when R187A and D192A were introduced in CYP2B1dH, the P450 expression and thermal stability were decreased. In addition, R187A showed increased activity with 7-EFC and decreased IC(50,BIF)/IC(50,4-CPI) compared with 2B1dH. Analysis of long range residue-residue interactions in the CYP2B4 crystal structures indicated strong hydrogen bonds involving Glu(149)-Asn(177)-Arg(187)-Tyr(190) and Asp(192)-Val(194), which were significantly-reduced/abolished by the Arg(187)-->Ala and Asp(192)-->Alasubstitutions, respectively.

The effect of mutation of F87 on the properties of CYP102A1-CYP4C7 chimeras: altered regiospecificity and substrate selectivity

CYP102A1 is a highly active water-soluble bacterial monooxygenase that contains both substrate-binding heme and diflavin reductase subunits, all in a single polypeptide that has been called a "self-sufficient enzyme." Several years ago we developed a procedure called "scanning chimeragenesis," where we focused on residues 73-82 of CYP102A1, which contact approximately 40% of the substrates palmitoleic acid and N-palmitoylglycine [Murataliev et al. (2004) Biochemistry 43:1771-1780]. These residues were replaced with the homologous residues of CYP4C7. In the current work, that study has been expanded to include residue 87. Phenylalanine 87 of wild-type CYP102A1 was replaced with the homologous residue of CYP4C7, leucine, as well as with alanine. The full-sized chimeric proteins C(73-78, F87L), C(73-78, F87A), C(75-80, F87L), C(75-80, F87A), C(78-82, F87L) and C(78-82, F87A) have been purified and characterized. Wild-type CYP102A1 is most active toward fatty acids (both lauric and palmitic acids produce omega-1, omega-2, and omega-3 hydroxylated fatty acids), but it also catalyzes the oxidation of farnesol to three products (2, 3- and 10,11-epoxyfarnesols and 9-hydroxyfarnesol). All of the F87-mutant chimeric proteins show dramatic decreases in activities with the natural CYP102A1 substrates. In contrast, C(78-82, F87A) and C(78-82, F87L) have markedly increased activities with farnesol, with the latter showing a 5.7-fold increase in catalytic activity as compared to wild-type CYP102A1. C(78-82, F87L) produces 10,11-epoxyfarnesol as the single primary metabolite. The results show that chimeragenesis involving only the second half of SRS-1 plus F87 is sufficient to change the substrate selectivity of CYP102A1 from fatty acids to farnesol and to produce a single primary product.

Resonance Raman studies of cytochrome P450 2B4 in its interactions with substrates and redox partners

Resonance Raman studies of P450 2B4 are reported for the substrate-free form and when bound to the substrates, benzphetamine (BZ) or butylated hydroxytoluene (BHT), the latter representing a substrate capable of inducing an especially effective conversion to the high-spin state. In addition to studies of the ferric resting state, spectra are acquired for the ferrous CO ligated form. Importantly, for the first time, the RR technique is effectively applied to interrogate the changes in active site structure induced by binding of cytochrome P450 reductase (CPR) and Mn(III) cytochrome b 5 (Mn cyt b 5); the manganese derivative of cyt b 5 was employed to avoid spectroscopic interferences. The results, consistent with early work on mammalian P450s, demonstrate that substrate structure has minimal effects on heme structure or the FeCO fragment of the ferrous CO derivatives. Similarly, the data indicate that the protein is flexible and that substrate binding does not exert significant strain on the heme peripheral groups, in contrast to P450 cam, where substantial effects on heme peripheral groups are seen. However, significant differences are observed in the RR spectra of P450 2B4 when bound with the different redox partners, indicating that the heme structure is clearly sensitive to perturbations near the proximal heme binding site. The most substantial changes are displacements of the peripheral vinyl groups toward planarity with the heme macrocycle by cyt b 5 but away from planarity by CPR. These changes can have an impact on heme reduction potential. Most interestingly, these RR results support an earlier observation that the combination of benzphetamine and cyt b 5 binding produce a synergy leading to unique active site structural changes when both are bound.

Structures of prostacyclin synthase and its complexes with substrate analog and inhibitor reveal a ligand-specific heme conformation change

Prostacyclin synthase (PGIS) is a cytochrome P450 (P450) enzyme that catalyzes production of prostacyclin from prostaglandin H(2). PGIS is unusual in that it catalyzes an isomerization rather than a monooxygenation, which is typical of P450 enzymes. To understand the structural basis for prostacyclin biosynthesis in greater detail, we have determined the crystal structures of ligand-free, inhibitor (minoxidil)-bound and substrate analog U51605-bound PGIS. These structures demonstrate a stereo-specific substrate binding and suggest features of the enzyme that facilitate isomerization. Unlike most microsomal P450s, where large substrate-induced conformational changes take place at the distal side of the heme, conformational changes in PGIS are observed at the proximal side and in the heme itself. The conserved and extensive heme propionate-protein interactions seen in all other P450s, which are largely absent in the ligand-free PGIS, are recovered upon U51605 binding accompanied by water exclusion from the active site. In contrast, when minoxidil binds, the propionate-protein interactions are not recovered and water molecules are largely retained. These findings suggest that PGIS represents a divergent evolution of the P450 family, in which a heme barrier has evolved to ensure strict binding specificity for prostaglandin H(2), leading to a radical-mediated isomerization with high product fidelity. The U51605-bound structure also provides a view of the substrate entrance and product exit channels.

Improvement of cyclophosphamide activation by CYP2B6 mutants: from in silico to ex vivo

Cyclophosphamide (CPA) is a chemotherapeutic agent that is primarily activated in the liver by cytochrome P4502B6 (CYP2B6) and then transported to the tumor via blood flow. To prevent deleterious secondary effects, P450-based gene-directed enzyme prodrug therapy (GDEPT) consists of expressing CYP2B6 in tumor cells before CPA treatment. Given the relatively low affinity of CYP2B6 for CPA, the aim of our work was to modify CYP2B6 to increase its catalytic efficiency (V(max)/K(m)) to metabolize CPA into 4'-OH CPA. A molecular model of CYP2B6 was built, and four residues in close contact with the substrate were subjected to mutagenesis. Canine CYP2B11 exhibiting a particularly low K(m) to CPA, the amino acids exclusively present in the CYP2B11 substrate recognition sequences were substituted in human CYP2B6. All mutants (n = 26) were expressed in Saccharomyces cerevisiae and their enzymatic constants (K(m), V(max)) evaluated using CPA as substrate. Five mutants exhibited a 2- to 3-fold higher catalytic efficiency than wild-type CYP2B6. A double mutant, comprising the two most effective mutations, showed a 4-fold increase in K(m)/V(max). Molecular dynamic simulations of several mutants were found to be consistent with the observed modifications in catalytic efficiency. Finally, expression of the CYP2B6 114V/477W double mutant, contrary to wt CYP2B6, allowed switching of a resistant human head and neck cancer cell line (A-253) into a sensitive cell line toward CPA. Thus, we were able to obtain a new efficient CYP2B6 mutant able to metabolize CPA, an important step in the GDEPT strategy for human cancer treatment.

Mechanism-based inactivation of human cytochromes p450s: experimental characterization, reactive intermediates, and clinical implications

The P450 type cytochromes are responsible for the metabolism of a wide variety of xenobiotics and endogenous compounds. Although P450-catalyzed reactions are generally thought to lead to detoxication of xenobiotics, the reactions can also produce reactive intermediates that can react with cellular macromolecules leading to toxicity or that can react with the P450s that form them leading to irreversible (i.e., mechanism-based) inactivation. This perspective describes the fundamentals of mechanism-based inactivation as it pertains to P450 enzymes. The experimental approaches used to characterize mechanism-based inactivators are discussed, and the criteria required for a compound to be classified as a mechanism-based inactivator are outlined. The kinetic scheme for mechanism-based inactivation and the calculation of the relevant kinetic constants that describe a particular inactivation event are presented. The structural aspects and important functional groups of several classes of molecules that have been found to impart mechanism-based inactivation upon metabolism by P450s such as acetylenes, thiol-containing compounds that include isothiocyanates, thiazolidinediones, and thiophenes, arylamines, quinones, furanocoumarins, and cyclic tertiary amines are described. Emphasis throughout this perspective is placed on more recent findings with human P450s where the site of modification, whether it be the apoprotein or the heme moiety, and, at least in part, the identity of the reactive intermediate responsible for the loss in P450 activity are known or inferred. Recent advances in trapping procedures as well as new methods for identification of reactive intermediates are presented. A variety of clinically important drugs that act as mechanism-based inactivators of P450s are discussed. The irreversible inactivation of human P450s by these drugs has the potential for causing serious drug-drug interactions that may have severe toxicological effects. The clinical significance of inactivating human P450s for improving drug efficacy as well as drug safety is discussed along with the potential for exploiting mechanism-based inactivators of P450s for therapeutic benefits.