Detection of substrate-dependent conformational changes in the P450 fold by nuclear magnetic resonance

Catalytic Site Constraints in the P450s Mediating Loganic Acid (7DLH) and Secologanic Acid Synthesis (SLAS) in Camptotheca

Terpene indole alkaloids (TIAs) are plant-derived natural products synthesized in low levels in medicinal plants such as Catharanthus roseus and Camptotheca acuminata. TIA pathways species utilize several CYP72A subfamily members to form loganic acid from 7-deoxyloganic acid (a simple hydroxylation) as well as secologanin and secologanic acid from loganin and loganic acid (a C-C bond scission). Divergences in the specificities of these P450s have allowed Camptotheca secologanic acid synthases (SLASs) to become bifunctional enzymes capable of performing both reactions. In contrast, Catharanthus 7-deoxyloganic acid hydroxylase (7DLH) and secologanin synthase (SLS) have remained monofunctional enzymes capable either of monooxygenation or C-C bond scission. Our in vitro reconstitutions have now demonstrated that Camptotheca also contains a monofunctional 7DLH capable only of hydroxylating 7-deoxyloganic acid. Mutageneses aimed at evaluating residues important for the tight specificity of Camptotheca 7DLH (CYP72A729) and the broad specificity of SLAS (CYP72A564) have identified several residues where reciprocal switches substantially affect their activities: Lys128His in 7DLH increases hydroxylation of 7-deoxyloganic acid, and His132Lys in SLAS decreases this hydroxylation and C-C bond scissions of loganic acid and loganin; Gly321Ser in 7DLH does not affect hydroxylation of 7-deoxyloganic acid, whereas Ser324Gly in SLAS significantly increases C-C bond scission of loganic acid; Asp332Glu in the acid-alcohol pair of 7DLH increases hydroxylation of 7-deoxyloganic acid, whereas Glu335Asp in SLAS completely eliminates both of its activities. These mutations that enhance or eliminate these respective activities have significant potential to aid engineering efforts aimed at increasing TIA production in cell cultures, microbial systems, and/or other plants.

Conformational heterogeneity suggests multiple substrate binding modes in CYP106A2

CYP106A2 (cytochrome P450meg) is a bacterial enzyme originally isolated from B. megaterium, and has been shown to hydroxylate a wide variety of substrates, including steroids. The regio- and stereochemistry of CYP106A2 hydroxylation has been shown to be dependent on a variety of factors, and hydroxylation often occurs at more than one site and/or with lack of stereospecificity for some substrates. Comprehensive backbone 15N, 1H and 13C resonance assignments based on multidimensional nuclear magnetic resonance (NMR) experiments performed with uniform and selective isotopically labeled CYP106A2 samples are reported herein, and broadening and splitting of resonances assigned to regions of the enzyme shown to be affected by substrate binding in other P450 enzymes indicate that substrate binding does not reduce structural heterogeneity as has been observed previously in P450 enzymes CYP101A1 and MycG. Paramagnetic relaxation enhancement (PRE) due to proximity between substrate protons and the heme iron were measured for three different substrates, and the relatively uniform nature of the PREs support the proposal that multiple substrate binding modes are occupied at saturating substrate concentrations.

Hydroxylation Regiochemistry Is Robust to Active Site Mutations in Cytochrome P450cam (CYP101A1)

Cytochrome P450_cam (CYP101A1) catalyzes the hydroxylation of d-camphor by molecular oxygen. The enzyme-catalyzed hydroxylation exhibits a high degree of regioselectivity and stereoselectivity, with a single major product, d-5-exo-hydroxycamphor, suggesting that the substrate is oriented to facilitate this specificity. In previous work, we used an elastic network model and perturbation response scanning to show that normal deformation modes of the enzyme structure are highly responsive not only to the presence of a substrate but also to the substrate orientation. This work examines the effects of mutations near the active site on substrate localization and orientation. The investigated mutations were designed to promote a change in substrate orientation and/or location that might give rise to different hydroxylation products, while maintaining the same carbon and oxygen atom balances as in the wild type (WT) enzyme. Computational experiments and parallel in vitro site-directed mutations of CYP101A1 were used to examine reaction products and enzyme activity. ¹H-¹⁵N TROSY-HSQC correlation maps were used to compare the computational results with detectable perturbations in the enzyme structure and dynamics. We found that all of the mutant enzymes retained the same regio- and stereospecificity of hydroxylation as the WT enzyme, with varying degrees of efficiency, which suggests that large portions of the enzyme have been subjected to evolutionary pressure to arrive at the appropriate sequence-structure combination for efficient 5-exo hydroxylation of camphor.

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.

Single mutations toggle the substrate selectivity of multifunctional Camptotheca secologanic acid synthases (CYP72As)

Terpene indole alkaloids (TIAs) are plant-derived specialized metabolites with widespread use in medicine. Species-specific pathways derive various TIAs from common intermediates, strictosidine or strictosidinic acid, produced by coupling tryptamine with secologanin or secologanic acid. The penultimate reaction in this pathway is catalyzed by either secologanin synthase (SLS) or secologanic acid synthase (SLAS) according to whether plants produce secologanin from loganin or secologanic acid from loganic acid. Previous work has identified SLSs and SLASs from different species, but the determinants of selectivity remain unclear. Here, combining molecular modeling, ancestral sequence reconstruction, and biochemical methodologies, we identified key residues that toggle SLS and SLAS selectivity in two CYP72A (cytochrome P450) subfamily enzymes from Camptotheca acuminata. We found that the positions of foremost importance are in substrate recognition sequence 1 (SRS1), where mutations to either of two adjacent histidine residues switched selectivity; His131Phe selects for and increases secologanin production whereas His132Asp selects for secologanic acid production. Furthermore, a change in SRS3 in the predicted substrate entry channel (Arg/Lys270Thr) and another in SRS4 at the start of the I-helix (Ser324Glu) decreased enzyme activity toward either substrate. We propose that the Camptotheca SLASs have maintained the broadened activities found in a common asterid ancestor, even as the Camptotheca lineage lost its ability to produce loganin while the campanulid and lamiid lineages specialized to produce secologanin by acquiring mutations in SRS1. The identification here of the residues essential for the broad substrate scope of SLASs presents opportunities for more tailored heterologous production of TIAs.

Updating the Paradigm: Redox Partner Binding and Conformational Dynamics in Cytochromes P450

This Account summarizes recent findings centered on the role that redox partner binding, allostery, and conformational dynamics plays in cytochrome P450 proton coupled electron transfer. P450s are one of Nature's largest enzyme families and it is not uncommon to find a P450 wherever substrate oxidation is required in the formation of essential molecules critical to the life of the organism or in xenobiotic detoxification. P450s can operate on a remarkably large range of substrates from the very small to the very large, yet the overall P450 three-dimensional structure is conserved. Given this conservation of structure, it is generally assumed that the basic catalytic mechanism is conserved. In nearly all P450s, the O₂ O-O bond must be cleaved heterolytically enabling one oxygen atom, the distal oxygen, to depart as water and leave behind a heme iron-linked O atom as the powerful oxidant that is used to activate the nearby substrate. For this process to proceed efficiently, externally supplied electrons and protons are required. Two protons must be added to the departing O atom while an electron is transferred from a redox partner that typically contains either a Fe₂S₂ or FMN redox center. The paradigm P450 used to unravel the details of these mechanisms has been the bacterial CYP101A1 or P450cam. P450cam is specific for its own Fe₂S₂ redox partner, putidaredoxin or Pdx, and it has long been postulated that Pdx plays an effector/allosteric role by possibly switching P450cam to an active conformation. Crystal structures, spectroscopic data, and direct binding experiments of the P450cam-Pdx complex provide some answers. Pdx shifts the conformation of P450cam to a more open state, a transition that is postulated to trigger the proton relay network required for O₂ activation. An essential part of this proton relay network is a highly conserved Asp (sometimes Glu) that is known to be critical for activity in a number of P450s. How this Asp and proton delivery networks are connected to redox partner binding is quite simple. In the closed state, this Asp is tied down by salt bridges, but these salt bridges are ruptured when Pdx binds, leaving the Asp free to serve its role in proton transfer. An alternative hypothesis suggests that a specific proton relay network is not really necessary. In this scenario, the Asp plays a structural role in the open/close transition and merely opening the active site access channel is sufficient to enable solvent protons in for O₂ protonation. Experiments designed to test these various hypotheses have revealed some surprises in both P450cam and other bacterial P450s. Molecular dynamics and crystallography show that P450cam can undergo rather significant conformational gymnastics that result in a large restructuring of the active site requiring multiple cis/trans proline isomerizations. It also has been found that X-ray driven substrate hydroxylation is a useful tool for better understanding the role that the essential Asp and surrounding residues play in catalysis. Here we summarize these recent results which provide a much more dynamic picture of P450 catalysis.

A cytochrome P450 CYP81AM1 from Tripterygium wilfordii catalyses the C-15 hydroxylation of dehydroabietic acid

A novel cytochrome P450 from Tripterygium wilfordii, CYP81AM1, specifically catalyses the C-15 hydroxylation of dehydroabietic acid. This is the first CYP450 to be found in plants with this function. Cytochrome P450 oxygenases (CYPs) play an important role in the post-modification in biosynthesis of plant bioactive terpenoids. Here, we found that CYP81AM1 can catalyze the formation of 15-hydroxydehydroabietic acid by in vitro enzymatic reactions and in vivo yeast feeding assays. This is the first study to show that CYP81 family enzymes are involved in the hydroxylation of abietane diterpenoids. At the same time, we found that CYP81AM1 could not catalyse abietatriene and dehydroabietinol, suggesting that the occurrence of the reaction may be related to the carboxyl group. Through molecular docking and site mutations, it was found that some amino acid sites (F104, K107) near the carboxyl group had an important effect on enzyme activity, also suggesting that the carboxyl group played an important role in the occurrence of the reaction.

19F-NMR reveals substrate specificity of CYP121A1 in Mycobacterium tuberculosis

Cytochromes P450 are versatile enzymes that function in endobiotic and xenobiotic metabolism and undergo meaningful structural changes that relate to their function. However, the way in which conformational changes inform the specific recognition of the substrate is often unknown. Here, we demonstrate the utility of fluorine (¹⁹F)-NMR spectroscopy to monitor structural changes in CYP121A1, an essential enzyme from Mycobacterium tuberculosis. CYP121A1 forms functional dimers that catalyze the phenol-coupling reaction of the dipeptide dicyclotyrosine. The thiol-reactive compound 3-bromo-1,1,1-trifluoroacetone was used to label an S171C mutation of the enzyme FG loop, which is located adjacent to the homodimer interface. Substrate titrations and inhibitor-bound ¹⁹F-NMR spectra indicate that ligand binding reduces conformational heterogeneity at the FG loop in both the dimer and in an engineered monomer of CYP121A1. However, only the dimer was found to promote a substrate-bound conformation that was preexisting in the substrate-free spectra, thus confirming a role for the dimer interface in dicyclotyrosine recognition. Moreover, ¹⁹F-NMR spectra in the presence of substrate analogs indicate the hydrogen-bonding feature of the dipeptide aromatic side chain as a dicyclotyrosine specificity criterion. This study demonstrates the utility of ¹⁹F-NMR as applied to a multimeric cytochrome P450, while also revealing mechanistic insights for an essential M. tuberculosis enzyme.

Partial Opening of Cytochrome P450cam (CYP101A1) Is Driven by Allostery and Putidaredoxin Binding

Cytochrome P450cam (CYP101A1) catalyzes the regio- and stereo-specific 5-exo-hydroxylation of camphor via a multistep catalytic cycle that involves two-electron transfer steps, with an absolute requirement that the second electron be donated by the ferrodoxin, putidaredoxin (Pdx). Whether P450cam, once camphor has bound to the active site and the substrate entry channel has closed, opens up upon Pdx binding, during the second electron transfer step, or it remains closed is still a matter of debate. A potential allosteric site for camphor binding has been identified and postulated to play a role in the binding of Pdx. Here, we have revisited paramagnetic NMR spectroscopy data and determined a heterogeneous ensemble of structures that explains the data, provides a complete representation of the P450cam/Pdx complex in solution, and reconciles alternative hypotheses. The allosteric camphor binding site is always present, and the conformational changes induced by camphor binding to this site facilitates Pdx binding. We also determined that the state to which Pdx binds comprises an ensemble of structures that have features of both the open and closed state. These results demonstrate that there is a finely balanced interaction between allosteric camphor binding and the binding of Pdx at high camphor concentrations.

Reconciling conformational heterogeneity and substrate recognition in cytochrome P450

Cytochrome P450, the ubiquitous metalloenzyme involved in detoxification of foreign components, has remained one of the most popular systems for substrate-recognition process. However, despite being known for its high substrate specificity, the mechanistic basis of substrate-binding by archetypal system cytochrome P450cam has remained at odds with the contrasting reports of multiple diverse crystallographic structures of its substrate-free form. Here, we address this issue by elucidating the probability of mutual dynamical transition to the other crystallographic pose of cytochrome P450cam and vice versa via unbiased all-atom computer simulation. A robust Markov state model, constructed using adaptively sampled 84-μs-long molecular dynamics simulation trajectories, maps the broad and heterogenous P450cam conformational landscape into five key substates. In particular, the Markov state model identifies an intermediate-assisted dynamic equilibrium between a pair of conformations of P450cam, in which the substrate-recognition sites remain "closed" and "open," respectively. However, the estimate of a significantly higher stationary population of closed conformation, coupled with faster rate of open → closed transition than its reverse process, dictates that the net conformational equilibrium would be swayed in favor of "closed" conformation. Together, the investigation quantitatively infers that although a potential substrate of cytochrome P450cam would, in principle, explore a diverse array of conformations of substrate-free protein, it would mostly encounter a "closed" or solvent-occluded conformation and hence would follow an induced-fit-based recognition process. Overall, the work reconciles multiple precedent crystallographic, spectroscopic investigations and establishes how a statistical elucidation of conformational heterogeneity in protein would provide crucial insights in the mechanism of potential substrate-recognition process.

Dynamics underlying hydroxylation selectivity of cytochrome P450cam

Structural heterogeneity and the dynamics of the complexes of enzymes with substrates can determine the selectivity of catalysis; however, fully characterizing how remains challenging as heterogeneity and dynamics can vary at the spatial level of an amino acid residue and involve rapid timescales. We demonstrate the nascent approach of site-specific two-dimensional infrared (IR) spectroscopy to investigate the archetypical cytochrome P450, P450cam, to better delineate the mechanism of the lower regioselectivity of hydroxylation of the substrate norcamphor in comparison to the native substrate camphor. Specific locations are targeted throughout the enzyme by selectively introducing cyano groups that have frequencies in a spectrally isolated region of the protein IR spectrum as local vibrational probes. Linear and two-dimensional IR spectroscopy were applied to measure the heterogeneity and dynamics at each probe and investigate how they differentiate camphor and norcamphor recognition. The IR data indicate that the norcamphor complex does not fully induce a large-scale conformational change to a closed state of the enzyme adopted in the camphor complex. Additionally, a probe directed at the bound substrate experiences rapidly interconverting states in the norcamphor complex that explain the hydroxylation product distribution. Altogether, the study reveals large- and small-scale structural heterogeneity and dynamics that could contribute to selectivity of a cytochrome P450 and illustrates the approach of site-selective IR spectroscopy to elucidate protein dynamics.

A dynamic understanding of cytochrome P450 structure and function through solution NMR

Many economically important biosyntheses incorporate regiospecific and stereospecific oxidations at unactivated carbons. Such oxidations are commonly catalyzed by cytochrome P450 monooxygenases, heme-containing enzymes that activate molecular oxygen while selectively binding and orienting the substrate for reaction. Despite the plethora of P450-catalyzed reactions, the P450 fold is highly conserved, and static structures are often insufficient for characterizing conformational states that contribute to specificity. High-resolution solution nuclear magnetic resonance (NMR) offers insights into dynamic processes and conformational changes that are required of a P450 in order to attain the combination of specificity and efficiency required for these reactions.

Linkage between Proximal and Distal Movements of P450cam Induced by Putidaredoxin

Putidaredoxin (Pdx) is the exclusive reductase and a structural effector for P450cam (CYP101A1). However, the mechanism of how Pdx modulates the conformational states of P450cam remains elusive. Here we report a putative communication pathway for the Pdx-induced conformational change in P450cam using results of double electron-electron resonance (DEER) spectroscopy and molecular dynamics simulations. Use of solution state DEER measurements allows us to observe subtle conformational changes in the internal helices in P450cam among closed, open, and P450cam-Pdx complex states. Molecular dynamics simulations and dynamic network analysis suggest that Pdx binding is coupled to small coordinated movements of several regions of P450cam, including helices C, B', I, G, and F. These changes provide a linkage between the Pdx binding site on the proximal side of the enzyme and helices F/G on the distal side and the site of the largest movement resulting from the Pdx-induced closed-to-open transition. This study provides a detailed rationale for how Pdx exerts its long-recognized effector function at the active site from its binding site on the opposite face of the enzyme.

Heterologous coexpression of the benzoate-para-hydroxylase CYP53B1 with different cytochrome P450 reductases in various yeasts

Cytochrome P450 monooxygenases (P450) are enzymes with high potential as biocatalysts for industrial applications. Their large-scale applications are, however, limited by instability and requirement for coproteins and/or expensive cofactors. These problems are largely overcome when whole cells are used as biocatalysts. We previously screened various yeast species heterologously expressing self-sufficient P450s for their potential as whole-cell biocatalysts. Most P450s are, however, not self-sufficient and consist of two or three protein component systems. Therefore, in the present study, we screened different yeast species for coexpression of P450 and P450-reductase (CPR) partners, using CYP53B1 from Rhodotorula minuta as an exemplary P450. The abilities of three different coexpressed CPR partners to support P450 activity were investigated, two from basidiomycetous origin and one from an ascomycete. The various P450-CPR combinations were cloned into strains of Saccharomyces cerevisiae, Kluyveromyces marxianus, Hansenula polymorpha, Yarrowia lipolytica and Arxula adeninivorans, using a broad-range yeast expression vector. The results obtained supported the previous finding that recombinant A. adeninivorans strains perform excellently as whole-cell biocatalysts. This study also demonstrated for the first time the P450 reductase activity of the CPRs from R. minuta and U. maydis. A very interesting observation was the variation in the supportive activity provided by the different reductase partners tested and demonstrated better P450 activity enhancement by a heterologous CPR compared to its natural partner CPR. This study highlights reductase selection as a critical variable for consideration in the pursuit of optimal P450-based catalytic systems. The usefulness of A. adeninivorans as both a host for recombinant P450s and whole-cell biocatalyst was emphasized, supporting earlier findings.

Crystal structure of bacterial CYP116B5 heme domain: New insights on class VII P450s structural flexibility and peroxygenase activity

Class VII cytochromes P450 are self-sufficient enzymes carrying a phthalate family oxygenase-like reductase domain and a P450 domain fused in a single polypeptide chain. The biocatalytic applications of CYP116B members are limited by the need of the NADPH cofactor and the lack of crystal structures as a starting point for protein engineering. Nevertheless, we demonstrated that the heme domain of CYP116B5 can use hydrogen peroxide as electron donor bypassing the need of NADPH. Here, we report the crystal structure of CYP116B5 heme domain in complex with histidine at 2.6 Å of resolution. The structure reveals the typical P450 fold and a closed conformation with an active site cavity of 284 ų in volume, accommodating a histidine molecule forming a hydrogen bond with the water molecule present as 6th heme iron ligand. MD simulations in the absence of any ligand revealed the opening of a tunnel connecting the active site to the protein surface through the movement of F-, G- and H-helices. A structural alignment with bacterial cytochromes P450 allowed the identification of amino acids in the proximal heme site potentially involved in peroxygenase activity. The availability of the crystal structure provides the bases for the structure-guided design of new biocatalysts.

Human cytochrome P450 enzymes bind drugs and other substrates mainly through conformational-selection modes

Cytochrome P450 (P450) enzymes are major catalysts involved in the oxidations of most drugs, steroids, carcinogens, fat-soluble vitamins, and natural products. The binding of substrates to some of the 57 human P450s and other mammalian P450s is more complex than a two-state system and has been proposed to involve mechanisms such as multiple ligand occupancy, induced-fit, and conformational-selection. Here, we used kinetic analysis of binding with multiple concentrations of substrates and computational modeling of these data to discern possible binding modes of several human P450s. We observed that P450 2D6 binds its ligand rolapitant in a mechanism involving conformational-selection. P450 4A11 bound the substrate lauric acid via conformational-selection, as did P450 2C8 with palmitic acid. Binding of the steroid progesterone to P450 21A2 was also best described by a conformational-selection model. Hexyl isonicotinate binding to P450 2E1 could be described by either a conformational-selection or an induced-fit model. Simulation of the binding of the ligands midazolam, bromocriptine, testosterone, and ketoconazole to P450 3A4 was consistent with an induced-fit or a conformational-selection model, but the concentration dependence of binding rates for varying both P450 3A4 and midazolam concentrations revealed discordance in the parameters, indicative of conformational-selection. Binding of the P450s 2C8, 2D6, 3A4, 4A11, and 21A2 was best described by conformational-selection, and P450 2E1 appeared to fit either mode. These findings highlight the complexity of human P450-substrate interactions and that conformational-selection is a dominant feature of many of these interactions.

What Your Crystal Structure Will Not Tell You about Enzyme Function

Enzyme function requires that enzyme structures be dynamic. Substrate binding, product release, and transition state stabilization typically involve different enzyme conformers. Furthermore, in multistep enzyme-catalyzed reactions, more than one enzyme conformation may be important for stabilizing different transition states. While X-ray crystallography provides the most detailed structural information of any current methodology, X-ray crystal structures of enzymes capture only those conformations that fit into the crystal lattice, which may or may not be relevant to function. Solution nuclear magnetic resonance (NMR) methods can provide an alternative approach to characterizing enzymes under nonperturbing and controllable conditions, allowing one to identify and localize dynamic processes that are important to function. However, many enzymes are too large for standard approaches to making sequential resonance assignments, a critical first step in analyzing and interpreting the wealth of information inherent in NMR spectra.

This Account describes our long-standing NMR-based research into structural and dynamic aspects of function in the cytochrome P450 monooxygenase superfamily. These heme-containing enzymes typically catalyze the oxidation of unactivated C–H and C=C bonds in a multitude of substrates, often with complete regio- and stereospecificity. Over 600 000 genes in GenBank have been assigned to P450s, yet all known P450 structures exhibit a highly conserved and unique fold. This combination of functional and structural conservation with a vast substrate clientele, each substrate having multiple possible sites for oxidation, makes the P450s a unique target for understanding the role of enzyme structure and dynamics in determining a particular substrate–product combination. P450s are large by solution NMR standards, requiring us to develop specialized approaches for making sequential resonance assignments and interpreting the spectral changes that occur as a function of changing conditions (e.g., oxidation and spin state changes, ligand, substrate or effector binding). Solution conformations are characterized by the fitting of residual dipolar couplings (RDCs) measured for sequence-specifically assigned amide N–H correlations to alignment tensors optimized in the course of restrained molecular dynamics (MD) simulations. The conformational ensembles obtained by such RDC-restrained simulations, which we call “soft annealing”, are then tested by site-directed mutation and spectroscopic and activity assays for relevance. These efforts have gained us insights into cryptic conformational changes associated with substrate and redox partner binding that were not suspected from crystal structures, but were shown by subsequent work to be relevant to function. Furthermore, it appears that many of these changes can be generalized to P450s besides those that we have characterized, providing guidance for enzyme engineering efforts. While past research was primarily directed at the more tractable prokaryotic P450s, our current efforts are aimed at medically relevant human enzymes, including CYP17A1, CYP2D6, and CYP3A4.

Lipid composition and macromolecular crowding effects on CYP2J2-mediated drug metabolism in nanodiscs

Lipid composition and macromolecular crowding are key external effectors of protein activity and stability whose role varies between different proteins. Therefore, it is imperative to study their effects on individual protein function. CYP2J2 is a membrane-bound cytochrome P450 in the heart involved in the metabolism of fatty acids and xenobiotics. In order to facilitate this metabolism, cytochrome P450 reductase (CPR), transfers electrons to CYP2J2 from NADPH. Herein, we use nanodiscs to show that lipid composition of the membrane bilayer affects substrate metabolism of the CYP2J2-CPR nanodisc (ND) system. Differential effects on both NADPH oxidation and substrate metabolism by CYP2J2-CPR are dependent on the lipid composition. For instance, sphingomyelin containing nanodiscs produced more secondary substrate metabolites than discs of other lipid compositions, implying a possible conformational change leading to processive metabolism. Furthermore, we demonstrate that macromolecular crowding plays a role in the lipid-solubilized CYP2J2-CPR system by increasing the K_m and decreasing the V_max , and effect that is size-dependent. Crowding also affects the CYP2J2-CPR-ND system by decreasing both the K_m and V_max for Dextran-based macromolecular crowding agents, implying an increase in substrate affinity but a lack of metabolism. Finally, protein denaturation studies show that crowding agents destabilize CYP2J2, while the multidomain protein CPR is stabilized. Overall, these studies are the first report on the role of the surrounding lipid environment and macromolecular crowding in modulating enzymatic function of CYP2J2-CPR membrane protein system.

Ligand and Redox Partner Binding Generates a New Conformational State in Cytochrome P450cam (CYP101A1)

It has become increasingly clear that cytochromes P450 can cycle back and forth between two extreme conformational states termed the closed and open states. In the well-studied cytochrome P450cam, the binding of its redox partner, putidaredoxin (Pdx), shifts P450cam toward the open state. Shifting to the open state is thought to be important in the formation of a proton relay network essential for O-O bond cleavage and formation of the active Fe(IV)═O intermediate. Another important intermediate is the oxy-P450cam complex when bound to Pdx. Trapping this intermediate in crystallo is challenging owing to its instability, but the CN^- complex is both stable and an excellent mimic of the O₂ complex. Here we present the P450cam-Pdx structure complexed with CN^-. CN^- results in large conformational changes including cis/trans isomerization of proline residues. Changes include large rearrangements of active-site residues and the formation of new active-site access channel that we have termed channel 2. The formation of channel 2 has also been observed in our previous molecular dynamics simulations wherein substrate binding to an allosteric site remote from the active site opens up channel 2. This new structure supports an extensive amount of previous work showing that distant regions of the structure are dynamically coupled and underscores the potentially important role that large conformational changes and dynamics play in P450 catalysis.

Entropic contribution to enhanced thermal stability in the thermostable P450 CYP119

The enhanced thermostability of thermophilic proteins with respect to their mesophilic counterparts is often attributed to the enthalpy effect, arising from strong interactions between protein residues. Intuitively, these strong interresidue interactions will rigidify the biomolecules. However, the present work utilizing neutron scattering and solution NMR spectroscopy measurements demonstrates a contrary example that the thermophilic cytochrome P450, CYP119, is much more flexible than its mesophilic counterpart, CYP101A1, something which is not apparent just from structural comparison of the two proteins. A mechanism to explain this apparent contradiction is that higher flexibility in the folded state of CYP119 increases its conformational entropy and thereby reduces the entropy gain during denaturation, which will increase the free energy needed for unfolding and thus stabilize the protein. This scenario is supported by thermodynamic data on the temperature dependence of unfolding free energy, which shows a significant entropic contribution to the thermostability of CYP119 and lends an added dimension to enhanced stability, previously attributed only to presence of aromatic stacking interactions and salt bridge networks. Our experimental data also support the notion that highly thermophilic P450s such as CYP119 may use a mechanism that partitions flexibility differently from mesophilic P450s between ligand binding and thermal stability.

Substrate-induced conformational change in cytochrome P450 OleP

The regulation of cytochrome P450 activity is often achieved by structural transitions induced by substrate binding. We describe the conformational transition experienced upon binding by the P450 OleP, an epoxygenase involved in oleandomycin biosynthesis. OleP bound to the substrate analog 6DEB crystallized in 2 forms: one with an ensemble of open and closed conformations in the asymmetric unit and another with only the closed conformation. Characterization of OleP-6DEB binding kinetics, also using the P450 inhibitor clotrimazole, unveiled a complex binding mechanism that involves slow conformational rearrangement with the accumulation of a spectroscopically detectable intermediate where 6DEB is bound to open OleP. Data reported herein provide structural snapshots of key precatalytic steps in the OleP reaction and explain how structural rearrangements induced by substrate binding regulate activity.-Parisi, G., Montemiglio, L. C., Giuffrè, A., Macone, A., Scaglione, A., Cerutti, G., Exertier, C., Savino, C., Vallone, B. Substrate-induced conformational change in cytochrome P450 OleP.

Some Surprising Implications of NMR-directed Simulations of Substrate Recognition and Binding by Cytochrome P450cam (CYP101A1)

Cytochrome P450_cam (CYP101A1) catalyzes the stereospecific 5-exo hydroxylation of d-camphor by molecular oxygen. Previously, residual dipolar couplings measured for backbone amide ¹H-¹⁵N correlations in both substrate-free and bound forms of CYP101A1 were used as restraints in soft annealing molecular dynamic simulations in order to identify average conformations of the enzyme with and without substrate bound. Multiple substrate-dependent conformational changes remote from the enzyme active site were identified, and site-directed mutagenesis and activity assays confirmed the importance of these changes in substrate recognition. The current work makes use of perturbation response scanning (PRS) and umbrella sampling molecular dynamic of the residual dipolar coupling-derived CYP101A1 structures to probe the roles of remote structural features in enforcing the regio- and stereospecific nature of the hydroxylation reaction catalyzed by CYP101A1. An improper dihedral angle Ψ was defined and used to maintain substrate orientation in the CYP101A1 active site, and it was observed that different values of Ψ result in different PRS response maps. Umbrella sampling methods show that the free energy of the system is sensitive to Ψ, and bound substrate forms an important mechanical link in the transmission of mechanical coupling through the enzyme structure. Finally, a qualitative approach to interpreting PRS maps in terms of the roles of secondary structural features is proposed.

Spectroscopic studies of the cytochrome P450 reaction mechanisms

The cytochrome P450 monooxygenases (P450s) are thiolate heme proteins that can, often under physiological conditions, catalyze many distinct oxidative transformations on a wide variety of molecules, including relatively simple alkanes or fatty acids, as well as more complex compounds such as steroids and exogenous pollutants. They perform such impressive chemistry utilizing a sophisticated catalytic cycle that involves a series of consecutive chemical transformations of heme prosthetic group. Each of these steps provides a unique spectral signature that reflects changes in oxidation or spin states, deformation of the porphyrin ring or alteration of dioxygen moieties. For a long time, the focus of cytochrome P450 research was to understand the underlying reaction mechanism of each enzymatic step, with the biggest challenge being identification and characterization of the powerful oxidizing intermediates. Spectroscopic methods, such as electronic absorption (UV-Vis), electron paramagnetic resonance (EPR), nuclear magnetic resonance (NMR), electron nuclear double resonance (ENDOR), Mössbauer, X-ray absorption (XAS), and resonance Raman (rR), have been useful tools in providing multifaceted and detailed mechanistic insights into the biophysics and biochemistry of these fascinating enzymes. The combination of spectroscopic techniques with novel approaches, such as cryoreduction and Nanodisc technology, allowed for generation, trapping and characterizing long sought transient intermediates, a task that has been difficult to achieve using other methods. Results obtained from the UV-Vis, rR and EPR spectroscopies are the main focus of this review, while the remaining spectroscopic techniques are briefly summarized. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Substrate recognition by two different P450s: Evidence for conserved roles in a common fold

Cytochrome P450 monooxygenases CYP101A1 and MycG catalyze regio- and stereospecific oxidations of their respective substrates, d-camphor and mycinamicin IV. Despite the low sequence homology between the two enzymes (29% identity) and differences in size and hydrophobicity of their substrates, the conformational changes that occur upon substrate binding in both enzymes as determined by solution NMR methods show some striking similarities. Many of the same secondary structural features in both enzymes are perturbed, suggesting the existence of a common mechanism for substrate binding and recognition in the P450 superfamily.

The Oxidation of Hydrophobic Aromatic Substrates by Using a Variant of the P450 Monooxygenase CYP101B1

The cytochrome P450 monooxygenase CYP101B1, from a Novosphingobium bacterium is able to bind and oxidise aromatic substrates but at a lower activity and efficiency than norisoprenoids and monoterpenoid esters. Histidine 85 of CYP101B1 aligns with tyrosine 96 of CYP101A1, which, in the latter enzyme forms the only hydrophilic interaction with its substrate, camphor. The histidine residue of CYP101B1 was mutated to phenylalanine with the aim of improving the activity of the enzyme for hydrophobic substrates. The H85F mutant lowered the binding affinity and activity of the enzyme for β-ionone and altered the oxidation selectivity. This variant also showed enhanced affinity and activity towards alkylbenzenes, styrenes and methylnaphthalenes. For example the rate of product formation for acenaphthene oxidation was improved sixfold to 245 nmol per nmol CYP per min. Certain disubstituted naphthalenes and substrates, such as phenylcyclohexane and biphenyls, were oxidised with lower activity by the H85F variant. Variants at H85 (A and G) designed to introduce additional space into the active site so as to accommodate these larger substrates did not improve the oxidation activity. As the H85F mutant of CYP101B1 improved the oxidation of hydrophobic substrates, this residue is likely to be in the substrate binding pocket or the access channel of the enzyme. The side chain of the histidine might interact with the carbonyl groups of the favoured norisoprenoid substrates of CYP101B1.

The Metal Drives the Chemistry: Dual Functions of Acireductone Dioxygenase

Acireductone dioxygenase (ARD) from the methionine salvage pathway (MSP) is a unique enzyme that exhibits dual chemistry determined solely by the identity of the divalent transition-metal ion (Fe²⁺ or Ni²⁺) in the active site. The Fe²⁺-containing isozyme catalyzes the on-pathway reaction using substrates 1,2-dihydroxy-3-keto-5-methylthiopent-1-ene (acireductone) and dioxygen to generate formate and the ketoacid precursor of methionine, 2-keto-4-methylthiobutyrate, whereas the Ni²⁺-containing isozyme catalyzes an off-pathway shunt with the same substrates, generating methylthiopropionate, carbon monoxide, and formate. The dual chemistry of ARD was originally discovered in the bacterium Klebsiella oxytoca, but it has recently been shown that mammalian ARD enzymes (mouse and human) are also capable of catalyzing metal-dependent dual chemistry in vitro. This is particularly interesting, since carbon monoxide, one of the products of off-pathway reaction, has been identified as an antiapoptotic molecule in mammals. In addition, several biochemical and genetic studies have indicated an inhibitory role of human ARD in cancer. This comprehensive review describes the biochemical and structural characterization of the ARD family, the proposed experimental and theoretical approaches to establishing mechanisms for the dual chemistry, insights into the mechanism based on comparison with structurally and functionally similar enzymes, and the applications of this research to the field of artificial metalloenzymes and synthetic biology.

Conformational Heterogeneity and the Affinity of Substrate Molecular Recognition by Cytochrome P450cam

The broad and variable substrate specificity of cytochrome P450 enzymes makes them a model system for studying the determinants of protein molecular recognition. The archetypal cytochrome P450cam (P450cam) is a relatively specific P450, a feature once attributed to the high rigidity of its active site. However, increasingly studies have provided evidence of the importance of conformational changes to P450cam activity. Here we used infrared (IR) spectroscopy to investigate the molecular recognition of P450cam. Toward this goal, and to assess the influence of a hydrogen bond (H-bond) between active site residue Y96 and substrates, two variants in which Y96 is replaced by a cyanophenyl (Y96CNF) or phenyl (Y96F) group were characterized in complexes with the substrates camphor, isoborneol, and camphane. These combinations allow for a comparison of complexes in which the moieties on both the protein and substrate can serve as a H-bond donor, acceptor, or neither. The IR spectra of heme-bound CO and the site-specifically incorporated CN of Y96CNF were analyzed to characterize the number and nature of environments in each protein, both in the free and bound states. Although the IR spectra do not support the idea that protein-substrate H-bonding is central to P450cam recognition, the data altogether suggest that the differing conformational heterogeneity in the active site of the P450cam variants and changes in heterogeneity upon binding of different substrates likely contribute to their variable affinities via a conformational selection mechanism. This study further extends our understanding of the molecular recognition of archetypal P450cam and demonstrates the application of IR spectroscopy combined with selective protein modification to delineate protein-ligand interactions.

Solution Conformations and Dynamics of Substrate-Bound Cytochrome P450 MycG

MycG is a P450 monooxygenase that catalyzes the sequential hydroxylation and epoxidation of mycinamicin IV (M-IV), the last two steps in the biosynthesis of mycinamicin II, a macrolide antibiotic isolated from Micromonospora griseorubida. The crystal structure of MycG with M-IV bound was previously determined but showed the bound substrate in an orientation that did not rationalize the observed regiochemistry of M-IV hydroxylation. Nuclear magnetic resonance paramagnetic relaxation enhancements provided evidence of an orientation of M-IV in the MycG active site more compatible with the observed chemistry, but substrate-induced changes in the enzyme structure were not characterized. We now describe the use of amide ¹H-¹⁵N residual dipolar couplings as experimental restraints in solvated "soft annealing" molecular dynamics simulations to generate solution structural ensembles of M-IV-bound MycG. Chemical shift perturbations, hydrogen-deuterium exchange, and ¹⁵N relaxation behavior provide insight into the dynamic and electronic perturbations in the MycG structure in response to M-IV binding. The solution and crystallographic structures are compared, and the possibility that the crystallographic orientation of bound M-IV represents an inhibitory mode is discussed.

The catalytic function of cytochrome P450 is entwined with its membrane-bound nature

Cytochrome P450, a family of monooxygenase enzymes, is organized as a catalytic metabolon, which requires enzymatic partners as well as environmental factors that tune its complex dynamic. P450 and its reducing counterparts—cytochrome P450-reductase and cytochrome b ₅—are membrane-bound proteins located in the cytosolic side of the endoplasmic reticulum. They are believed to dynamically associate to form functional complexes. Increasing experimental evidence signifies the role(s) played by both protein-protein and protein-lipid interactions in P450 catalytic function and efficiency. However, the biophysical challenges posed by their membrane-bound nature have severely limited high-resolution understanding of the molecular interfaces of these interactions. In this article, we provide an overview of the current knowledge on cytochrome P450, highlighting the environmental factors that are entwined with its metabolic function. Recent advances in structural biophysics are also discussed, setting up the bases for a new paradigm in the study of this important class of membrane-bound enzymes.

NADH reduction of nitroaromatics as a probe for residual ferric form high-spin in a cytochrome P450

The existence of a substrate-sensitive equilibrium between high spin (S=5/2) and low spin (S=1/2) ferric iron is a well-established phenomenon in the cytochrome P450 (CYP) superfamily, although its origins are still a subject of discussion. A series of mutations that strongly perturb the spin state equilibrium in the camphor hydroxylase CYP101A1 were recently described (Colthart et al., Sci. Rep. 6, 22035 (2016)). Wild type CYP101A1 as well as some CYP101A1 mutants are herein shown to be capable of catalyzing the reduction of nitroacetophenones by NADH to the corresponding anilino compounds (nitroreductase or NRase activity). The distinguishing characteristic between those mutants that catalyze the reduction and those that cannot appears to be the extent to which residual high spin form exists in the absence of the native substrate d-camphor, with those showing the largest spin state shifts upon camphor binding also exhibiting NRase activity. Optical and EPR spectroscopy was used to further examine these phenomena. These results suggest that reduction of nitroaromatics may provide a useful probe of residual high spin states in the CYP superfamily. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Triphenyltin degradation and proteomic response by an engineered Escherichia coli expressing cytochrome P450 enzyme

Although triphenyltin (TPT) degradation pathway has been determined, information about the enzyme and protein networks involved was severely limited. To this end, a cytochrome P450 hydroxylase (CYP450) gene from Bacillus thuringiensis was cloned and expressed in Escherichia coli BL21 (DE3), namely E. coli pET32a-CYP450, whose dosage at 1gL^-1 could degrade 54.6% TPT at 1mgL^-1 within 6 d through attacking the carbon-tin bonds of TPT by CYP450. Sequence analysis verified that the CYP450 gene had a 1214bp open reading frame, encoding a protein with 404 amino acids. Proteomic analysis determined that 60 proteins were significantly differentially regulated expression in E. coli pET32a-CYP450 after TPT degradation. The up-regulated proteins enriched in a network related to transport, cell division, biosynthesis of amino acids and secondary metabolites, and microbial metabolism in diverse environments. The current findings demonstrated for the first time that P450 received electrons transferring from NADH could effectively cleave carbon-metal bonds.

Recent Structural Insights into Cytochrome P450 Function

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

Site-Specific Characterization of Cytochrome P450cam Conformations by Infrared Spectroscopy

Conformational changes are central to protein function but challenging to characterize with both high spatial and temporal precision. The inherently fast time scale and small chromophores of infrared (IR) spectroscopy are well-suited for characterization of potentially rapidly fluctuating environments, and when frequency-resolved probes are incorporated to overcome spectral congestion, enable characterization of specific sites in proteins. We selectively incorporated p-cyanophenylalanine (CNF) as a vibrational probe at five distinct locations in the enzyme cytochrome P450cam and used IR spectroscopy to characterize the environments in substrate and/or ligand complexes reflecting those in the catalytic cycle. Molecular dynamics (MD) simulations were performed to provide a structural basis for spectral interpretation. Together the experimental and simulation data suggest that the CN frequencies are sensitive to both long-range influences, resulting from the particular location of a residue within the enzyme, as well as short-range influences from hydrogen bonding and packing interactions. The IR spectra demonstrate that the environments and effects of substrate and/or ligand binding are different at each position probed and also provide evidence that a single site can experience multiple environments. This study illustrates how IR spectroscopy, when combined with the spectral decongestion and spatial selectivity afforded by CNF incorporation, provides detailed information about protein structural changes that underlie function.