Substrate routes to the buried active site may vary among cytochromes P450: mutagenesis of the F-G region in P450 2B1

Multiresolution molecular dynamics simulations reveal the interplay between conformational variability and functional interactions in membrane-bound cytochrome P450 2B4

Cytochrome P450 2B4 (CYP 2B4) is one of the best-characterized CYPs and serves as a key model system for understanding the mechanisms of microsomal class II CYPs, which metabolize most known drugs. The highly flexible nature of CYP 2B4 is apparent from crystal structures that show the active site with either a wide open or a closed heme binding cavity. Here, we investigated the conformational ensemble of the full-length CYP 2B4 in a phospholipid bilayer, using multiresolution molecular dynamics (MD) simulations. Coarse-grained MD simulations revealed two predominant orientations of CYP 2B4's globular domain with respect to the bilayer. Their refinement by atomistic resolution MD showed adaptation of the enzyme's interaction with the lipid bilayer, leading to open configurations that facilitate ligand access to the heme binding cavity. CAVER analysis of enzyme tunnels, AquaDuct analysis of water routes, and Random Acceleration Molecular Dynamics simulations of ligand dissociation support the conformation-dependent passage of molecules between the active site and the protein surroundings. Furthermore, simulation of the re-entry of the inhibitor bifonazole into the open conformation of CYP 2B4 resulted in binding at a transient hydrophobic pocket within the active site cavity that may play a role in substrate binding or allosteric regulation. Together, these results show how the open conformation of CYP 2B4 facilitates the binding of substrates from and release of products to the membrane, whereas the closed conformation prolongs the residence time of substrates or inhibitors and selectively allows the passage of smaller reactants via the solvent and water channels.

Exploring the biotransformation of N-(2-hydroxyphenyl)-2-propylpentanamide (an aryl valproic acid derivative) by CYP2C11, using in silico predictions and in vitro studies

OBJECTIVES

N-(2-hydroxyphenyl)-2-propylpentanamide (HO-AAVPA), a derivative of valproic acid (VPA), has been proposed as a potential anticancer agent due to its improved antiproliferative effects in some cancer cell lines. Although there is evidence that VPA is metabolized by cytochrome P450 2C11 rat isoform, HO-AAVPA CYP-mediated metabolism has not yet been fully explored. Therefore, in this work, the biotransformation of HO-AAVPA by CYP2C11 was investigated.

METHODS

Kinetic parameters and spectral interaction between HO-AAVPA and CYP were evaluated using rat liver microsomes. The participation of CYP2C11 in metabolism of HO-AAVPA was confirmed by cimetidine (CIM) inhibition assay. Docking and molecular dynamics simulations coupled to MMGBSA methods were used in theoretical study.

KEY FINDINGS

HO-AAVPA is metabolized by CYP enzymes (K_M  = 38.94 µm), yielding a hydroxylated metabolite according to its HPLC retention time (5.4 min) and MS analysis (252.2 m/z). In addition, CIM inhibition in rat liver microsomes (K_i  = 59.23 µm) confirmed that CYP2C11 is mainly involved in HO-AAVPA metabolism. Furthermore, HO-AAVPA interacts with CYP2C11 as a type I ligand. HO-AAVPA is stabilized at the CYP2C11 ligand recognition site through a map of interactions similar to other typical CYP2C11 substrates.

CONCLUSION

Therefore, rat liver CYP2C11 isoform is able to metabolize HO-AAVPA.

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.

ChannelsDB: database of biomacromolecular tunnels and pores

ChannelsDB (http://ncbr.muni.cz/ChannelsDB) is a database providing information about the positions, geometry and physicochemical properties of channels (pores and tunnels) found within biomacromolecular structures deposited in the Protein Data Bank. Channels were deposited from two sources; from literature using manual deposition and from a software tool automatically detecting tunnels leading to the enzymatic active sites and selected cofactors, and transmembrane pores. The database stores information about geometrical features (e.g. length and radius profile along a channel) and physicochemical properties involving polarity, hydrophobicity, hydropathy, charge and mutability. The stored data are interlinked with available UniProt annotation data mapping known mutation effects to channel-lining residues. All structures with channels are displayed in a clear interactive manner, further facilitating data manipulation and interpretation. As such, ChannelsDB provides an invaluable resource for research related to deciphering the biological function of biomacromolecular channels.

Atomic-level description of protein-lipid interactions using an accelerated membrane model

Peripheral membrane proteins are structurally diverse proteins that are involved in fundamental cellular processes. Their activity of these proteins is frequently modulated through their interaction with cellular membranes, and as a result techniques to study the interfacial interaction between peripheral proteins and the membrane are in high demand. Due to the fluid nature of the membrane and the reversibility of protein-membrane interactions, the experimental study of these systems remains a challenging task. Molecular dynamics simulations offer a suitable approach to study protein-lipid interactions; however, the slow dynamics of the lipids often prevents sufficient sampling of specific membrane-protein interactions in atomistic simulations. To increase lipid dynamics while preserving the atomistic detail of protein-lipid interactions, in the highly mobile membrane-mimetic (HMMM) model the membrane core is replaced by an organic solvent, while short-tailed lipids provide a nearly complete representation of natural lipids at the organic solvent/water interface. Here, we present a brief introduction and a summary of recent applications of the HMMM to study different membrane proteins, complementing the experimental characterization of the presented systems, and we offer a perspective of future applications of the HMMM to study other classes of membrane proteins. This article is part of a Special Issue entitled: Membrane proteins edited by J.C. Gumbart and Sergei Noskov.

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.

Efficient biosynthesis of rare natural product scopolamine using E. coli cells expressing a S14P/K97A mutant of hyoscyamine 6β-hydroxylase AaH6H

Hyoscyamine 6β-hydroxylase (H6H, EC 1.14.11.11), an α-ketoglutarate dependent dioxygenase catalyzes the hydroxylation of (-)-hyoscyamine and the subsequent epoxidation of 6β-hydroxyhyoscyamine to form scopolamine, a valuable natural alkaloid. In this study, random mutagenesis and site-directed saturation mutagenesis were used to enhance the hydroxylation activity of H6H from Anisodus acutangulus (AaH6H). A double mutant, AaH6HM1 (S14P/K97A), showed a 3.4-fold improved hydroxylation activity compared with the wild-type enzyme, and the in vivo epoxidation activity was also improved by 2.3 times. After 34h cultivation of Escherichia coli cells harboring Aah6hm1 in a 5-L bioreactor with a working volume of 3L, scopolamine was produced via a single-enzyme-mediated two-step transformation from 500mgL(-1) (-)-hyoscyamine in 97% conversion, and 1.068g of the product were isolated, corresponding to a space-time yield of 251mgL(-1)d(-1). This study shows that the protein engineering of some key enzymes is a promising and effective way for improving the production of rare natural products such as scopolamine.

Directed evolution of cytochrome P450 enzymes for biocatalysis: exploiting the catalytic versatility of enzymes with relaxed substrate specificity

Cytochrome P450 enzymes are renowned for their ability to insert oxygen into an enormous variety of compounds with a high degree of chemo- and regio-selectivity under mild conditions. This property has been exploited in Nature for an enormous variety of physiological functions, and representatives of this ancient enzyme family have been identified in all kingdoms of life. The catalytic versatility of P450s makes them well suited for repurposing for the synthesis of fine chemicals such as drugs. Although these enzymes have not evolved in Nature to perform the reactions required for modern chemical industries, many P450s show relaxed substrate specificity and exhibit some degree of activity towards non-natural substrates of relevance to applications such as drug development. Directed evolution and other protein engineering methods can be used to improve upon this low level of activity and convert these promiscuous generalist enzymes into specialists capable of mediating reactions of interest with exquisite regio- and stereo-selectivity. Although there are some notable successes in exploiting P450s from natural sources in metabolic engineering, and P450s have been proven repeatedly to be excellent material for engineering, there are few examples to date of practical application of engineered P450s. The purpose of the present review is to illustrate the progress that has been made in altering properties of P450s such as substrate range, cofactor preference and stability, and outline some of the remaining challenges that must be overcome for industrial application of these powerful biocatalysts.

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

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

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

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

Quantum mechanically derived AMBER-compatible heme parameters for various states of the cytochrome P450 catalytic cycle

Molecular mechanics (MM) methods are computationally affordable tools for screening chemical libraries of novel compounds for sites of P450 metabolism. One challenge for MM methods has been the absence of a consistent and transferable set of parameters for the heme within the P450 active site. Experimental data indicate that mammalian P450 enzymes vary greatly in the size, architecture, and plasticity of their active sites. Thus, obtaining X-ray-based geometries for the development of accurate MM parameters for the major classes of hepatic P450 remains a daunting task. Our previous work with preliminary gas-phase quantum mechanics (QM)-derived atomic partial charges greatly improved the accuracy of docking studies of raloxifene to CYP3A4. We have therefore developed and tested a consistent set of transferable MM parameters based on gas-phase QM calculations of two model systems of the heme-a truncated (T-HM) and a full (F-HM) for four states of the P450 catalytic cycle. Our results indicate that the use of the atomic partial charges from the F-HM further improves the accuracy of docked predictions for raloxifene to CYP3A4. Different patterns for substrate docking are also observed depending on the choice of heme model and state. Newly parameterized heme models are tested in implicit and explicitly solvated MD simulations in the absence and presence of enzyme structures, for CYP3A4, and appear to be stable on the nanosecond simulation timescale. The new force field for the various heme states may aid the community for simulations of P450 enzymes and other heme-containing enzymes.

Membrane position of ibuprofen agrees with suggested access path entrance to cytochrome P450 2C9 active site

Cytochrome P450 2C9 (CYP2C9) is a membrane-anchored human microsomal protein involved in the drug metabolism in liver. CYP2C9 consists of an N-terminal transmembrane anchor and a catalytic cytoplasmic domain. While the structure of the catalytic domain is well-known from X-ray experiments, the complete structure and its incorporation into the membrane remains unsolved. We constructed an atomistic model of complete CYP2C9 in a dioleoylphosphatidylcholine membrane and evolved it by molecular dynamics simulations in explicit water on a 100+ ns time-scale. The model agrees well with known experimental data about membrane positioning of cytochromes P450. The entry to the substrate access channel is proposed to be facing the membrane interior while the exit of the product egress channel is situated above the interface pointing toward the water phase. The positions of openings of the substrate access and product egress channels correspond to free energy minima of CYP2C9 substrate ibuprofen and its metabolite in the membrane, respectively.

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.

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.

Probing possible egress channels for multiple ligands in human CYP3A4: a molecular modeling study

Human cytochrome P450 (CYP) 3A4 extensively contributes to metabolize 50% of the marketed drugs. Recently, a CYP3A4 structure with two molecules of ketoconazole (2KT) was identified. However, channels for egresses of these inhibitors are unexplored. Thus, we applied molecular dynamics simulations followed by channel analyses. Two simulations of empty and 2KT-bound CYP3A4 results revealed the multiple ligand-induced conformational changes in channel forming regions, which appear to be important for the regulation of channels. In addition, we observed that the channel-3 entrance is closed due to the large structural deviation of the key residues from Phe-cluster. F215 and F220 are known as entrance blockers of channel-2 in metyrapone-bound CYP3A4. Currently, F220 blocks the channel-3 along with F213 and F241. Therefore, it suggested that channel-1 and 2 could potentially serve as egress routes for 2KT. It is also supported by the results from MOLAxis analyses, in which the frequency of channel occurrence and bottleneck radius during simulation favor channel-1 and 2. Several bottleneck residues of these channels may have critical roles in 2KT egresses, especially S119. Our modeling study for multiple ligand-channeling of CYP3A4 could be very helpful to gain new insights into channel selectivity of CYP3A4.

Pathways and mechanisms for product release in the engineered haloalkane dehalogenases explored using classical and random acceleration molecular dynamics simulations

Eight mutants of the DhaA haloalkane dehalogenase carrying mutations at the residues lining two tunnels, previously observed by protein X-ray crystallography, were constructed and biochemically characterized. The mutants showed distinct catalytic efficiencies with the halogenated substrate 1,2,3-trichloropropane. Release pathways for the two dehalogenation products, 2,3-dichloropropane-1-ol and the chloride ion, and exchange pathways for water molecules, were studied using classical and random acceleration molecular dynamics simulations. Five different pathways, denoted p1, p2a, p2b, p2c, and p3, were identified. The individual pathways showed differing selectivity for the products: the chloride ion releases solely through p1, whereas the alcohol releases through all five pathways. Water molecules play a crucial role for release of both products by breakage of their hydrogen-bonding interactions with the active-site residues and shielding the charged chloride ion during its passage through a hydrophobic tunnel. Exchange of the chloride ions, the alcohol product, and the waters between the buried active site and the bulk solvent can be realized by three different mechanisms: (i) passage through a permanent tunnel, (ii) passage through a transient tunnel, and (iii) migration through a protein matrix. We demonstrate that the accessibility of the pathways and the mechanisms of ligand exchange were modified by mutations. Insertion of bulky aromatic residues in the tunnel corresponding to pathway p1 leads to reduced accessibility to the ligands and a change in mechanism of opening from permanent to transient. We propose that engineering the accessibility of tunnels and the mechanisms of ligand exchange is a powerful strategy for modification of the functional properties of enzymes with buried active sites.

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.

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.

Possible pathway(s) of metyrapone egress from the active site of cytochrome P450 3A4: a molecular dynamics simulation

To identify a possible pathway(s) for metyrapone egress from the active site of P450 3A4, a 5-ns conventional molecular dynamics simulation followed by steered molecular dynamics simulations was performed on the complex with metyrapone. The steered molecular dynamics simulations showed that metyrapone egress via channel 1, threading through the B-C loop, only required a relatively small rupture force and small displacement of residues, whereas egress via the third channel, between helix I and helices F' and G', required a relatively large force and perturbation of helices I, B', and C. The conventional dynamics simulation indicated that channel 2, located between the beta1 sheet, B-B' loop, and F'-G' region, is closed because of the movement of residues in the mouth of this channel. The findings suggest that channel 1 can be used for metyrapone egress, whereas both channel 2 and channel 3 have a low probability of serving as an exit channel for metyrapone. In addition, residues F108 and I120 appear to act as two gatekeepers to prevent the inhibitor from leaving the active site. These results are in agreement with previous site-directed mutagenesis experiments.

Re-engineering cytochrome P450 2B11dH for enhanced metabolism of several substrates including the anti-cancer prodrugs cyclophosphamide and ifosfamide

Based on recent directed evolution of P450 2B1, six P450 2B11 mutants at three positions were created in an N-terminal modified construct termed P450 2B11dH and characterized for enzyme catalysis using five substrates. Mutant I209A demonstrated a 3.2-fold enhanced k(cat)/K(m) for 7-ethoxy-4-trifluoromethylcourmarin O-deethylation, largely due to a dramatic decrease in K(m) (0.72 microM vs. 18 microM). I209A also demonstrated enhanced selectivity for testosterone 16beta-hydroxylation over 16alpha-hydroxylation. In contrast, V183L showed a 4-fold increased k(cat) for 7-benzyloxyresorufin debenzylation and a 4.7-fold increased k(cat)/K(m) for testosterone 16alpha-hydroxylation. V183L also displayed a 1.7-fold higher k(cat)/K(m) than P450 2B11dH with the anti-cancer prodrugs cyclophosphamide and ifosfamide, resulting from a approximately 4-fold decrease in K(m). Introduction of the V183L mutation into full-length P450 2B11 did not enhance the k(cat)/K(m). Overall, the re-engineered P450 2B11dH enzymes exhibited enhanced catalytic efficiency with several substrates including the anti-cancer prodrugs.

Investigation of the role of cytochrome P450 2B4 active site residues in substrate metabolism based on crystal structures of the ligand-bound enzyme

Based on the X-ray crystal structures of 4-(4-chlorophenyl)imidazole (4-CPI)- and bifonazole (BIF)-bound P450 2B4, eight active site mutants at six positions were created in an N-terminal modified construct termed 2B4dH and characterized for enzyme inhibition and catalysis. I363A showed a >4-fold decrease in differential inhibition by BIF and 4-CPI (IC(50,BIF)/IC(50,4-CPI)). F296A, T302A, I363A, V367A, and V477A showed a 2-fold decreased k(cat) for 7-ethoxy-4-trifluoromethylcoumarin O-deethylation, whereas V367A and V477F showed an altered K(m). T302A, V367L, and V477A showed >4-fold decrease in total testosterone hydroxylation, whereas I363A, V367A, and V477F showed altered stereo- and regioselectivity. Interestingly, I363A showed a 150-fold enhanced k(cat)/K(m) with testosterone, and yielded a new metabolite. Furthermore, testosterone docking into three-dimensional models of selected mutants based on the 4-CPI-bound structure suggested a re-positioning of residues 363 and 477 to yield products. In conclusion, our results suggest that the 4-CPI-bound 2B4dH/H226Y crystal structure is an appropriate model for predicting enzyme catalysis.

Structure-function analysis of cytochromes P450 2B

In the last 4 years, breakthroughs were made in the field of P450 2B (CYP2B) structure-function through determination of one ligand-free and two inhibitor-bound X-ray crystal structures of CYP2B4, which revealed many of the structural features required for binding ligands of different size and shape. Large conformational changes of several plastic regions of CYP2B4 can dramatically reshape the active site of the enzyme to fit the size and shape of the bound ligand without perturbing the overall P450 fold. Solution biophysical studies using isothermal titration calorimetry (ITC) have revealed the large difference in the thermodynamic parameters of CYP2B4 in binding inhibitors of different ring chemistry and side chains. Other studies have revealed that the effects of site-specific mutations on steady-state kinetic parameters and mechanism-based inactivation are often substrate dependent. These findings agree with the structural data that the enzymes adopt different conformations to bind various ligands. Thus, the substrate specificity of an individual enzyme is determined not only by active site residues but also non-active site residues that modulate conformational changes that are important for substrate access and rearrangement of the active site to accommodate the bound substrate.

Alternative splicing within the human cytochrome P450 superfamily with an emphasis on the brain: the convolution continues

The human cytochrome P450 (CYP) superfamily of enzymes regulate hepatic phase 1 drug metabolism and subsequently play a significant role in pharmacokinetics, drug discovery and drug development. Alternative splicing of the cytochrome CYP gene transcripts enhances gene diversity and may play a role in transcriptional regulation of certain CYP proteins. Tissue-specific alternative splicing of CYPs is significant for its potential to add greater dimension to differential drug metabolism in hepatic and extrahepatic tissues, such as the brain, and to our understanding of the CYP family. This review provides an overview of tissue-specific splicing patterns, splicing types, regulation and the functional diversities between liver and splice variant CYP proteins and further explores the relevance of tissue-specific alternative splicing of CYPs in the nervous system.

Conformational flexibility of mammalian cytochrome P450 2B4 in binding imidazole inhibitors with different ring chemistry and side chains. Solution thermodynamics and molecular modeling

Recent x-ray structures of cytochrome P450 2B4 (CYP2B4) reveal an open form that undergoes a large-scale structural transition to a closed form upon binding to 4-(4-chlorophenyl)imidazole (4-CPI). Here, we report for the first time a complete solution thermodynamic study using isothermal titration calorimetry supported by spectroscopic studies to elucidate the conformational flexibility of CYP2B4 in binding imidazole inhibitors with different ring chemistry and side chains: 4-CPI, 1-benzylimidazole (1-BI), 1-CPI, 4-phenylimidazole (4-PI), 1-(2-(benzyloxy)ethyl)imidazole (BEI), and 1-PI. Each of the inhibitors induced type II spectral changes, and IC50 values for enzyme inhibition ranged from 0.1 to 2.4 microM, following the order 1-BI < 4-CPI < 1-CPI < 4-PI < BEI < 1-PI. Calorimetric titrations using monomeric enzyme yielded a 1:1 binding stoichiometry, with the associated KD values ranging from 0.3 to 4.8 microM and following the same rank order as the IC50 values. Changes in enthalpy at 25 degrees C ranged from -6.5 to -8.8 kcal mol(-1). The largest difference in binding entropy (+5.9 versus -4.1 cal mol(-1) K(-1)) was observed between 4-CPI and BEI, respectively, with a 2-fold difference in heat capacity changes (-604 versus -331 cal mol(-1) K(-1)), which is inferred to result from the reduction of apolar surface area of the enzyme ensuing from a conformational change upon 4-CPI binding. Accessibility to acrylamide of the only tryptophan (Trp121), which is located in helix C, was greatly decreased only in protein bound to 4-CPI. Steric restrictions hindered the perfect docking of only BEI to the closed conformation of the enzyme. The thermodynamic signature obtained for structurally similar inhibitors suggests remarkable plasticity of CYP2B4.

Use of directed evolution of mammalian cytochromes P450 for investigating the molecular basis of enzyme function and generating novel biocatalysts

Directed evolution has been successfully applied to the design of industrial biocatalysts for enhanced catalytic efficiency and stability, and for examining the molecular basis of enzyme function. Xenobiotic-metabolizing mammalian cytochromes P450 with their catalytic versatility and broad substrate specificity offer the possibility of widespread applications in industrial synthesis, medicine, and bioremediation. However, the requirement for NADPH-cytochrome P450 reductase, often cytochrome b5, and an expensive cofactor, NADPH, complicates the design of mammalian P450 enzymes as biocatalysts. Recently, Guengerich and colleagues have successfully performed directed evolution of P450s 1A2 and 2A6 initially by using colony-based colorimetric and genotoxicity screening assays, respectively, followed by in vitro fluorescence-based activity screening assays. More recently, our laboratory has developed a fluorescence-based in vitro activity screening assay system for enhanced catalytic activity of P450s 2B1 and 3A4. The studies indicate an important role of amino acid residues outside of the active site, which would be difficult to target by other methods. The approach can now be expanded to design these as well as new P450s using more targeted substrates of environmental, industrial, and medical importance.

Possible pathway(s) of testosterone egress from the active site of cytochrome P450 2B1: a steered molecular dynamics simulation

To probe the possible substrate exit channel(s) in cytochrome P450 (P450) 2B1 and to clarify the role of residues previously identified by site-directed mutagenesis, a homology model was constructed based on the X-ray crystal structure of a P450 2B4-inhibitor complex. Testosterone was docked into the active site of P450 2B1 and was then pulled out through three putative channels using steered molecular dynamics simulations. The results indicated that of the three channels, the "solvent channel," lined by helices E, F, and I and the beta3 hairpin, required the largest rupture force and backbone motion, which rendered it unlikely as an exit route. The relatively small rupture forces and backbone motions for the other two channels suggested them as possible candidates for testosterone passage. The opening of channel 1, located between helices G and I and the B'-C loop, is characterized by rotation of the aromatic ring of Phe297 together with a bending of the B'-C loop. The opening of channel 2, penetrating through the B'-C loop/B' helix, is achieved by an expansion of this region and a small displacement of the backbone. Interestingly, during the egress of testosterone along channel 1, Phe297 and Phe108 appear to act as two clamps to stabilize testosterone binding and prevent it from leaving the active site. Phe115 acts as a gatekeeper for channel 2. These results are in agreement with previous site-directed mutagenesis experiments.

Directed Evolution of Mammalian Cytochrome P450 2B1

Cytochrome P450 2B1 has been subjected to directed evolution to investigate the role of amino acid residues outside of the active site and to engineer novel, more active P450 catalysts. A high throughput screening system was developed to measure H(2)O(2)-supported oxidation of the marker fluorogenic substrate 7-ethoxy-4-trifluoromethylcoumarin (7-EFC). Random mutagenesis by error-prone polymerase chain reaction and activity screening were optimized using the L209A mutant of P450 2B1 in an N-terminally modified construct with a C-terminal His tag (P450 2B1dH). Two rounds of mutagenesis and screening and one subcloning step yielded the P450 2B1 quadruple mutant V183L/F202L/L209A/S334P, which demonstrated a 6-fold higher k(cat) than L209A. Further random or site-directed mutagenesis did not improve the activity. When assayed in an NADPH-supported reconstituted system, V183L/L209A demonstrated lower 7-EFC oxidation than L209A. Therefore, F202L/L209A/S334P was generated, which showed a 2.5-fold higher k(cat)/K(m) for NADPH-dependent 7-EFC oxidation than L209A. F202L/L209A/S334P also showed enhanced catalytic efficiency with 7-benzyloxyresorufin, benzphetamine, and testosterone, and a 10-fold increase in stereoselectivity for testosterone 16alpha-versus 16beta-hydroxylation compared with 2B1dH. Enhanced catalytic efficiency of F202L/L209A/S334P was also retained in the full-length P450 2B1 background with 7-EFC and testosterone as substrates. Finally, the individual mutants were tested for metabolism of the anti-cancer prodrugs cyclophosphamide and ifosfamide. Several of the mutants showed increased metabolism via the therapeutically beneficial 4-hydroxylation pathway, with L209A/S334P showing 2.8-fold enhancement of k(cat)/K(m) with cyclophosphamide and V183L/L209A showing 3.5-fold enhancement with ifosfamide. Directed evolution can thus be used to enhance P450 2B1 catalytic efficiency across a panel of substrates and to identify functionally important residues distant from the active site.

Functional role of residues in the helix B' region of cytochrome P450 2B1

Comparison of several recently determined X-ray crystal structures of mammalian cytochrome P450 family 2 enzymes suggests considerable movement of helix B' when ligands bind. To investigate the functional role of helix B' in P450 2B1, residues 100-109 were substituted with alanine and phenylalanine. Kinetic properties were examined with the typical 2B substrates 7-benzyloxyresorufin, 7-ethoxy-4-trifluoromethylcoumarin, benzphetamine, and testosterone. Several mutants showed 2- to 3-fold changes in k(cat) values and significant differences in catalytic efficiencies among the substrates examined, consistent with structural information suggesting that the helix B' region can adopt multiple conformations with different contact residues depending on the substrate. Homology modeling of P450 2B1 was performed based on an inhibitor-bound P450 2B4 structure, and the docking analyses were consistent with experimental results. The findings suggest that residues in the helix B' region affect regio- and stereoselective oxidation in P450 family 2 enzymes as well as substrate entry.

Cytochrome P450: what have we learned and what are the future issues?

The cytochrome P450 (P450) field came out of interest in the metabolism of drugs, carcinogens, and steroids, which remain major focal points. Over the years we have come to understand the P450 system components, the multiplicity of P450s, and many aspects of the regulation of the genes and also the catalytic mechanism. Many crystal structures are now becoming available. The significance of P450 in in vivo metabolism is appreciated, particularly in the context of pharmacogenetics. Current scientific issues involve posttranslational modification, gene regulation, component interactions, structures of P450 complexed with ligands, details of high-valent oxygen chemistry, the nature and influence of rate-limiting steps, greater details about some reaction steps, cooperativity, and the relevance of P450 variations to cancer risk. Some emerging research areas involve new methods of analysis of ligand interactions, roles of conformational changes linked to individual reaction steps, functions of orphan P450s, "molecular breeding" of new P450 functions and enhanced activity, and the utilization of P450s in chemical synthesis.

Structure of mammalian cytochrome P450 2B4 complexed with 4-(4-chlorophenyl)imidazole at 1.9-A resolution: insight into the range of P450 conformations and the coordination of redox partner binding

A 1.9-A molecular structure of the microsomal cytochrome P450 2B4 with the specific inhibitor 4-(4-chlorophenyl)imidazole (CPI) in the active site was determined by x-ray crystallography. In contrast to the previous experimentally determined 2B4 structure, this complex adopted a closed conformation similar to that observed for the mammalian 2C enzymes. The differences between the open and closed structures of 2B4 were primarily limited to the lid domain of helices F through G, helices B' and C, the N terminus of helix I, and the beta(4) region. These large-scale conformational changes were generally due to the relocation of conserved structural elements toward each other with remarkably little remodeling at the secondary structure level. For example, the F' and G' helices were maintained with a sharp turn between them but are placed to form the exterior ceiling of the active site in the CPI complex. CPI was closely surrounded by residues from substrate recognition sites 1, 4, 5, and 6 to form a small, isolated hydrophobic cavity. The switch from open to closed conformation dramatically relocated helix C to a more proximal position. As a result, heme binding interactions were altered, and the putative NADPH-cytochrome P450 reductase binding site was reformed. This suggests a structural mechanism whereby ligand-induced conformational changes may coordinate catalytic activity. Comparison of the 2B4/CPI complex with the open 2B4 structure yields insights into the dynamics involved in substrate access, tight inhibitor binding, and coordination of substrate and redox partner binding.

Mutagenesis and molecular dynamics suggest structural and functional roles for residues in the N-terminal portion of the cytochrome P450 2B1 I helix

To investigate their potential roles in ligand access, binding, and subsequent metabolism, residues in the N-terminal portion of the cytochrome P450 2B1 I helix were mutated to alanine and phenylalanine. Of the 18 mutants from E286 to S294 only 7 yielded holoprotein in an Escherichia coli expression system. Substitutions at positions 289, 290, 292, and 294 caused >/= 2-fold changes in kcat and/or Km for two or more of the 2B1 substrates examined, testosterone, 7-ethoxy-4-trifluoromethylcoumarin, 7-benzyloxyresorufin, and benzphetamine. I290 substitutions had the largest effects on steady-state parameters for three substrates and increased benzphetamine affinity. Steered molecular dynamics simulations of testosterone egress along the I helix identified hydrophobic interactions with I290, L293, and S294 and water bridges to E286 and S294. Sensitivity of holoprotein formation to substitution and effects on substrate binding and metabolism suggest structural and functional roles for residues in the N-terminus of the cytochrome P450 2B1 I helix.

A rational approach to Re-engineer cytochrome P450 2B1 regioselectivity based on the crystal structure of cytochrome P450 2C5

The regioselectivity for progesterone hydroxylation by cytochrome P450 2B1 was re-engineered based on the x-ray crystal structure of cytochrome P450 2C5. 2B1 is a high K(m) progesterone 16alpha-hydroxylase, whereas 2C5 is a low K(m) progesterone 21-hydroxylase. Initially, nine individual 2B1 active-site residues were changed to the corresponding 2C5 residues, and the mutants were purified from an Escherichia coli expression system and assayed for progesterone hydroxylation. At 150 microm progesterone, I114A, F297G, and V363L showed 5-15% of the 21-hydroxylase activity of 2C5, whereas F206V showed high activity for an unknown product and a 13-fold decrease in K(m). Therefore, a quadruple mutant, I114A/F206V/F297G/V363L (Q), was constructed that showed 60% of 2C5 progesterone 21-hydroxylase activity and 57% regioselectivity. Based on their 2C5-like testosterone hydroxylation profiles, S294D and I477F alone and in combination were added to the quadruple mutant. All three mutants showed enhanced regioselectivity (70%) for progesterone 21-hydroxylation, whereas only Q/I477F had a higher k(cat). Finally, the remaining three single mutants, V103I, V367L, and G478V, were added to Q/I477F and Q/S294D/I477F, yielding seven additional multiple mutants. Among these, Q/V103I/S294D/I477F showed the highest k(cat) (3-fold higher than that of 2C5) and 80% regioselectivity for progesterone 21-hydroxylation. Docking of progesterone into a three-dimensional model of this mutant indicated that 21-hydroxylation is favored. In conclusion, a systematic approach to convert P450 regioselectivity across subfamilies suggests that active-site residues are mainly responsible for regioselectivity differences between 2B1 and 2C5 and validates the reliability of 2B1 models based on the crystal structure of 2C5.