Effect of Redox Partner Binding on Cytochrome P450 Conformational Dynamics

3site Multisubstrate-Bound State of Cytochrome P450cam

We identified a multisubstrate-bound state, hereby referred as a 3site state, in cytochrome P450cam via integrating molecular dynamics simulation with nuclear magnetic resonance (NMR) pseudocontact shift measurements. The 3site state is a result of simultaneous binding of three camphor molecules in three locations around P450cam: (a) in a well-established "catalytic" site near heme, (b) in a kink-separated "waiting" site along channel-1, and (c) in a previously reported "allosteric" site at E, F, G, and H helical junctions. These three spatially distinct binding modes in the 3site state mutually communicate with each other via homotropic allostery and act cooperatively to render P450cam functional. The 3site state shows a significantly superior fit with NMR pseudo contact shift (PCS) data with a Q-score of 0.045 than previously known bound states and consists of D251 free of salt-bridges with K178 and R186, rendering the enzyme functionally primed. To date, none of the reported cocomplex of P450cam with its redox partner putidaredoxin (pdx) has been able to match solution NMR data and controversial pdx-induced opening of P450cam's channel-1 remains a matter of recurrent discourse. In this regard, inclusion of pdx to the 3site state is able to perfectly fit the NMR PCS measurement with a Q-score of 0.08 and disfavors the pdx-induced opening of channel-1, reconciling previously unexplained remarkably fast hydroxylation kinetics with a koff of 10.2 s-1. Together, our findings hint that previous experimental observations may have inadvertently captured the 3site state as an in vitro solution state, instead of the catalytic state alone, and provided a distinct departure from the conventional understanding of cytochrome P450.

Solar-powered P450 catalysis: Engineering electron transfer pathways from photosynthesis to P450s

With the increasing focus on green chemistry, biocatalysis is becoming more widely used in the pharmaceutical and other chemical industries for sustainable production of high value and structurally complex chemicals. Cytochrome P450 monooxygenases (P450s) are attractive biocatalysts for industrial application due to their ability to transform a huge range of substrates in a stereo- and regiospecific manner. However, despite their appeal, the industrial application of P450s is limited by their dependence on costly reduced nicotinamide adenine dinucleotide phosphate (NADPH) and one or more auxiliary redox partner proteins. Coupling P450s to the photosynthetic machinery of a plant allows photosynthetically-generated electrons to be used to drive catalysis, overcoming this cofactor dependency. Thus, photosynthetic organisms could serve as photobioreactors with the capability to produce value-added chemicals using only light, water, CO2 and an appropriate chemical as substrate for the reaction/s of choice, yielding new opportunities for producing commodity and high-value chemicals in a carbon-negative and sustainable manner. This review will discuss recent progress in using photosynthesis for light-driven P450 biocatalysis and explore the potential for further development of such systems.

Redox partner recognition and selectivity of cytochrome P450lin (CYP111A1)

The strict requirement of cytochrome P450cam for its native ferredoxin redox partner, putidaredoxin (Pdx), is not exhibited by any other known cytochrome P450 (CYP) system and the molecular details of redox partner selectivity are still not completely understood. We therefore examined the selectivity of a related Pseudomonas cytochrome P450, P450lin, by testing its activity with non-native redox partners. We found that P450lin could utilize Arx, the native redox partner of CYP101D1, to enable turnover of its substrate, linalool, while Pdx showed limited activity. Arx exhibited a higher sequence similarity to P450lins native redox partner, linredoxin (Ldx) than Pdx, including several residues that are believed to be at the interface of the two proteins, based on the P450cam-Pdx complex structure. We therefore mutated Pdx to resemble Ldx and Arx and found that a double mutant, D38L/∆106, displayed higher activity than Arx. In addition, Pdx D38L/∆106 does not induce a low-spin shift in linalool bound P450lin but does destabilize the P450lin-oxycomplex. Together our results suggest that P450lin and its redox partners may form a similar interface to P450cam-Pdx, but the interactions that allow for productive turnover are different.

Cooperative Substrate Binding Controls Catalysis in Bacterial Cytochrome P450terp (CYP108A1)

Despite being one of the most well-studied aspects of cytochrome P450 chemistry, important questions remain regarding the nature and ubiquity of allosteric regulation of catalysis. The crystal structure of a bacterial P450, P450terp, in the presence of substrate reveals two binding sites, one above the heme in position for regioselective hydroxylation and another in the substrate access channel. Unlike many bacterial P450s, P450terp does not exhibit an open to closed conformational change when substrate binds; instead, P450terp uses the second substrate molecule to hold the first substrate molecule in position for catalysis. Spectral titrations clearly show that substrate binding to P450terp is cooperative with a Hill coefficient of 1.4 and is supported by isothermal titration calorimetry. The importance of the allosteric site was explored by a series of mutations that weaken the second site and that help hold the first substrate in position for proper catalysis. We further measured the coupling efficiency of both the wild-type (WT) enzyme and the mutant enzymes. While the WT enzyme exhibits 97% efficiency, each of the variants showed lower catalytic efficiency. Additionally, the variants show decreased spin shifts upon binding of substrate. These results are the first clear example of positive homotropic allostery in a class 1 bacterial P450 with its natural substrate. Combined with our recent results from P450cam showing complex substrate allostery and conformational dynamics, our present study with P450terp indicates that bacterial P450s may not be as simple as once thought and share complex substrate binding properties usually associated with only mammalian P450s.

Unravelling the role of transient redox partner complexes in P450 electron transfer mechanics

The molecular evolution of cytochromes P450 and associated redox-driven oxidative catalysis remains a mystery in biology. It is widely believed that sterol 14α-demethylase (CYP51), an essential enzyme of sterol biosynthesis, is the ancestor of the whole P450 superfamily given its conservation across species in different biological kingdoms. Herein we have utilized X-ray crystallography, molecular dynamics simulations, phylogenetics and electron transfer measurements to interrogate the nature of P450-redox partner binding using the naturally occurring fusion protein, CYP51-ferredoxin found in the sterol-producing bacterium Methylococcus capsulatus. Our data advocates that the electron transfer mechanics in the M. capsulatus CYP51-ferredoxin fusion protein involves an ensemble of ferredoxin molecules in various orientations and the interactions are transient. Close proximity of ferredoxin, however, is required to complete the substrate-induced large-scale structural switch in the P450 domain that enables proton-coupled electron transfer and subsequent oxygen scission and catalysis. These results have fundamental implications regarding the early evolution of electron transfer proteins and for the redox reactions in the early steps of sterol biosynthesis. They also shed new light on redox protein mechanics and the subsequent diversification of the P450 electron transfer machinery in nature.

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.

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.

Active Site Hydrogen Bonding Induced in Cytochrome P450cam by Effector Putidaredoxin

Cytochrome P450s are diverse and powerful catalysts that can activate molecular oxygen to oxidize a wide variety of substrates. Catalysis relies on effective uptake of two electrons and two protons. For cytochrome P450cam, an archetypal member of the superfamily, the second electron must be supplied by the redox partner putidaredoxin (Pdx). Pdx also plays an effector role beyond electron transfer, but after decades the mechanism remains under investigation. We applied infrared spectroscopy to heme-ligated CN^- to examine the influence of Pdx binding. The results indicate that Pdx induces the population of a conformation wherein the CN^- ligand forms a strong hydrogen bond to a solvent water molecule, experimentally corroborating the formation of a proposed proton delivery network. Further, characterization of T252A P450cam implicates the side chain of Thr252 in regulating the population equilibrium of hydrogen-bonded states within the P450cam/Pdx complex, which could underlie its role in directing activated oxygen toward product formation and preventing reaction uncoupling through peroxide release.

Pervasive cooperative mutational effects on multiple catalytic enzyme traits emerge via long-range conformational dynamics

Multidimensional fitness landscapes provide insights into the molecular basis of laboratory and natural evolution. To date, such efforts usually focus on limited protein families and a single enzyme trait, with little concern about the relationship between protein epistasis and conformational dynamics. Here, we report a multiparametric fitness landscape for a cytochrome P450 monooxygenase that was engineered for the regio- and stereoselective hydroxylation of a steroid. We develop a computational program to automatically quantify non-additive effects among all possible mutational pathways, finding pervasive cooperative signs and magnitude epistasis on multiple catalytic traits. By using quantum mechanics and molecular dynamics simulations, we show that these effects are modulated by long-range interactions in loops, helices and β-strands that gate the substrate access channel allowing for optimal catalysis. Our work highlights the importance of conformational dynamics on epistasis in an enzyme involved in secondary metabolism and offers insights for engineering P450s.

Connecting conformational dynamics and epistasis has so far been limited to a few proteins and a single fitness trait. Here, the authors provide evidence of positive epistasis on multiple catalytic traits in the evolution and dynamics of engineered cytochrome P450 monooxygenase, offering insights for in silico protein design.

Unexpected Differences between Two Closely Related Bacterial P450 Camphor Monooxygenases

The bacterial cytochrome P450cam catalyzes the oxidation of camphor to 5-exo-hydroxycamphor as the first step in the oxidative assimilation of camphor as a carbon/energy source. CYP101D1 is another bacterial P450 that catalyzes the same reaction. A third P450 (P450tcu) has recently been discovered that has ≈86% sequence identity to P450cam as well as very similar enzymatic properties. P450tcu, however, exhibits three unusual features not found in P450cam. First, we observe product in at least two orientations in the X-ray structure that indicates that, unlike the case for P450cam, X-ray-generated reducing equivalents can drive substrate hydroxylation in crystallo. We postulate, on the basis of molecular dynamics simulations, that greater flexibility in P450tcu enables easier access of protons to the active site and, together with X-ray driven reduction, results in O₂ activation and substrate hydroxylation. Second, the characteristic low-spin to high-spin transition when camphor binds occurs immediately with P450cam but is very slow in P450tcu. Third, isothermal titration calorimetry shows that in P450cam substrate binding is entropically driven with a ΔH of >0 while in P450tcu with a ΔH of <0 with a more modest change in -TΔS. These results indicate that despite nearly identical structures and enzymatic properties, these two P450s exhibit quite different properties most likely related to differences in conformational dynamics.

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.

Hole Hopping through Cytochrome P450

High-potential iron-oxo species are intermediates in the catalytic cycles of oxygenase enzymes. They can cause heme degradation and irreversible oxidation of nearby amino acids. We have proposed that there are protective mechanisms in which hole hopping from oxidized hemes through tryptophan/tyrosine chains generates a surface-exposed amino-acid oxidant that could be rapidly disarmed by reaction with cellular reductants. In investigations of cytochrome P450_BM3, we identified Trp96 as a critical residue that could play such a protective role. This Trp is cation-π paired with Arg398 in 81% of mammalian P450s. Here we report on the effect of the Trp/Arg cation-π interaction on Trp96 formal potentials as well as on electronic coupling strengths between Trp96 and the heme both for wild type cytochrome P450 and selected mutants. Mutation of Arg398 to His, which decreases the Trp96 formal potential, increases Trp-heme electronic coupling; however, surprisingly, the rate of phototriggered electron transfer from a Ru-sensitizer (through Trp96) to the P450_BM3 heme was unaffected by the Arg398His mutation. We conclude that Trp96 has moved away from Arg398, suggesting that the protective mechanism for P450s with this Trp-Arg pair is conformationally gated.

Dual Character of Reactive Oxygen, Nitrogen, and Halogen Species: Endogenous Sources, Interconversions and Neutralization

Oxidative stress resulting from accumulation of reactive oxygen, nitrogen, and halogen species (ROS, RNS, and RHS, respectively) causes the damage of cells and biomolecules. However, over the long evolutionary time, living organisms have developed the mechanisms for adaptation to oxidative stress conditions including the activity of the antioxidant system (AOS), which maintains low intracellular levels of RONS (ROS and RNS) and RHS. Moreover, living organisms have adapted to use low concentrations of these electrophiles for the regulation of cell functions through the reversible post-translational chemical modifications of redox-sensitive amino acid residues in intracellular effectors of signal transduction pathways (protein kinases and protein phosphatases), transcription factors, etc. An important fine-tuning mechanism that ensures involvement of RONS and RHS in the regulation of physiological processes is interconversion between different reactive species. This review focuses on the complex networks of interacting RONS and RHS types and their endogenous sources, such as NOX family of NADPH oxidases, complexes I and III of the mitochondrial electron transport chain, NO synthases, cytochrome P450-containing monooxygenase system, xanthine oxidoreductase, and myeloperoxidases. We highlight that kinetic parameters of reactions involving RONS and RHS determine the effects of these reactive species on cell functions. We also describe the functioning of enzymatic and non-enzymatic AOS components and the mechanisms of RONS and RHS scavenging under physiological conditions. We believe that analysis of interactions between RONS and relationships between different endogenous sources of these compounds will contribute to better understanding of their role in the maintenance of cell redox homeostasis as well as initiation and progression of diseases.

Fungal P450 monooxygenases - the diversity in catalysis and their promising roles in biocontrol activity

The fungal P_450s catalyze vital monooxygenation reactions in primary and secondary metabolism, which may lead to the production of diverse secondary metabolites. Many of these, such as from the family of trichothecenes, involve in biocontrol activities. The diversified nature of fungal P₄₅₀ monooxygenases makes their host organisms adoptable to various ecological niches. The available genome data analysis provided an insight into the activity and mechanisms of the fungal P_450s. However, still more structural and functional studies are needed to elucidate the details of its catalytic mechanism, and the advance studies are also required to decipher further about their dynamic role in various aspects of trichothecene oxygenations. This mini review will provide updated information on different fungal P₄₅₀ monooxygenases, their genetic diversity, and their role in catalyzing various biochemical reactions leading to the production of plant growth promoting secondary metabolites.

Proton relay network in P450cam formed upon docking of putidaredoxin

Cytochromes P450 are versatile heme-based enzymes responsible for vital life processes. Of these, P450cam (substrate camphor) has been most studied. Despite this, precise mechanisms of the key O─O cleavage step remain partly elusive to date; effects observed in various enzyme mutants remain partly unexplained. We have carried out extended (to 1000 ns) MM-MD and follow-on quantum mechanics/molecular mechanics computations, both on the well-studied FeOO state and on Cpd(0) (compound 0). Our simulations include (all camphor-bound): (a) WT (wild type), FeOO state. (b) WT, Cpd(0). (c) Pdx (Putidaredoxin, redox partner of P450)-docked-WT, FeOO state. (d) Pdx-docked WT, Cpd(0). (e) Pdx-docked T252A mutant, Cpd(0). Among our key findings: (a) Effect of Pdx docking appears to go far beyond that indicated in prior studies: it leads to specific alterations in secondary structure that create the crucial proton relay network. (b) Specific proton relay networks we identify are: FeOO(H)⋯T252⋯nH ₂ O⋯D251 in WT; FeOO(H)⋯nH ₂ O⋯D251 in T252A mutant; both occur with Pdx docking. (c) Direct interaction of D251 with -FeOOH is, respectively, rare/frequent in WT/T252A mutant. (d) In WT, T252 is in the proton relay network. (e) Positioning of camphor appears significant: when camphor is part of H-bonding network, second protonation appears to be facilitated.

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.

An Intermediate Conformational State of Cytochrome P450cam-CN in Complex with Putidaredoxin

Cytochrome P450cam is an archetypal example of the vast family of heme monooxygenases and serves as a model for an enzyme that is highly specific for both its substrate and reductase. During catalysis, it undergoes significant conformational changes of the F and G helices upon binding its substrate and redox partner, putidaredoxin (Pdx). Recent studies have shown that Pdx binding to the closed camphor-bound form of ferric P450cam results in its conversion to a fully open state. However, during catalytic turnover, it remains unclear whether this same conformational change also occurs or whether it is coupled to the formation of the critical compound I intermediate. Here, we have examined P450cam bound simultaneously by camphor, CN^-, and Pdx as a mimic of the catalytically competent ferrous oxy-P450cam-Pdx state. The combined use of double electron-electron resonance and molecular dynamics showed direct observation of intermediate conformational states of the enzyme upon CN^- and subsequent Pdx binding. This state is coupled to the movement of the I helix and residues at the active site, including Arg-186, Asp-251, and Thr-252. These movements enable occupation of a water molecule that has been implicated in proton delivery and peroxy bond cleavage to give compound I. These findings provide a detailed understanding of how the Pdx-induced conformational change may sequentially promote compound I formation followed by product release, while retaining stereoselective hydroxylation of the substrate of this highly specific monooxygenase.

Use of Computational Biochemistry for Elucidating Molecular Mechanisms of Nitric Oxide Synthase

Nitric oxide (NO) is an essential signaling molecule in the regulation of multiple cellular processes. It is endogenously synthesized by NO synthase (NOS) as the product of L-arginine oxidation to L-citrulline, requiring NADPH, molecular oxygen, and a pterin cofactor. Two NOS isoforms are constitutively present in cells, nNOS and eNOS, and a third is inducible (iNOS). Despite their biological relevance, the details of their complex structural features and reactivity mechanisms are still unclear. In this review, we summarized the contribution of computational biochemistry to research on NOS molecular mechanisms. We described in detail its use in studying aspects of structure, dynamics and reactivity. We also focus on the numerous outstanding questions in the field that could benefit from more extensive computational investigations.

Site-Specific 1D and 2D IR Spectroscopy to Characterize the Conformations and Dynamics of Protein Molecular Recognition

Proteins exist as ensembles of interconverting states on a complex energy landscape. A complete, molecular-level understanding of their function requires knowledge of the populated states and thus the experimental tools to characterize them. Infrared (IR) spectroscopy has an inherently fast time scale that can capture all states and their dynamics with, in principle, bond-specific spatial resolution, and 2D IR methods that provide richer information are becoming more routine. Although application of IR spectroscopy for investigation of proteins is challenged by spectral congestion, the issue can be overcome by site-specific introduction of amino acid side chains that have IR probe groups with frequency-resolved absorptions, which furthermore enables selective characterization of different locations in proteins. Here, we briefly introduce the biophysical methods and summarize the current progress toward the study of proteins. We then describe our efforts to apply site-specific 1D and 2D IR spectroscopy toward elucidation of protein conformations and dynamics to investigate their involvement in protein molecular recognition, in particular mediated by dynamic complexes: plastocyanin and its binding partner cytochrome f, cytochrome P450s and substrates or redox partners, and Src homology 3 domains and proline-rich peptide motifs. We highlight the advantages of frequency-resolved probes to characterize specific, local sites in proteins and uncover variation among different locations, as well as the advantage of the fast time scale of IR spectroscopy to detect rapidly interconverting states. In addition, we illustrate the greater insight provided by 2D methods and discuss potential routes for further advancement of the field of biomolecular 2D IR spectroscopy.

Mapping the Substrate Recognition Pathway in Cytochrome P450

Cytochrome P450s are ubiquitous metalloenzymes involved in the metabolism and detoxification of foreign components via catalysis of the hydroxylation reactions of a vast array of organic substrates. However, the mechanism underlying the pharmaceutically critical process of substrate access to the catalytic center of cytochrome P450 is a long-standing puzzle, further complicated by the crystallographic evidence of a closed catalytic center in both substrate-free and substrate-bound cytochrome P450. Here, we address a crucial question whether the conformational heterogeneity prevalent in cytochrome P450 translates to heterogeneous pathways for substrate access to the catalytic center of these metalloenzymes. By atomistically capturing the full process of spontaneous substrate association from bulk solvent to the occluded catalytic center of an archetypal system P450cam in multi-microsecond-long continuous unbiased molecular dynamics simulations, we here demonstrate that the substrate recognition in P450cam always occurs through a single well-defined dominant pathway. The simulated final bound pose resulting from these unguided simulations is in striking resemblance with the crystallographic bound pose. Each individual binding trajectory reveals that the substrate, initially placed at random locations in bulk solvent, spontaneously lands on a single key channel on the protein-surface of P450cam and resides there for an uncharacteristically long period, before correctly identifying the occluded target-binding cavity. Surprisingly, the passage of substrate to the closed catalytic center is not accompanied by any large-scale opening in protein. Rather, the unbiased simulated trajectories (∼57 μs) and underlying Markov state model, in combination with free-energy analysis, unequivocally show that the substrate recognition process in P450cam needs a substrate-induced side-chain displacement coupled with a complex array of dynamical interconversions of multiple metastable substrate conformations. Further, the work reconciles multiple precedent experimental and theoretical observations on P450cam and establishes a comprehensive view of substrate-recognition in cytochrome P450 that only occurs via substrate-induced structural rearrangements.

Conformational Change Induced by Putidaredoxin Binding to Ferrous CO-ligated Cytochrome P450cam Characterized by 2D IR Spectroscopy

The importance of conformational dynamics to protein function is now well-appreciated. An outstanding question is whether they are involved in the effector role played by putidaredoxin (Pdx) in its reduction of the O₂ complex of cytochrome P450cam (P450cam), an archetypical member of the cytochrome P450 superfamily. Recent studies have reported that binding of Pdx induces a conformational change from a closed to an open state of ferric P450cam, but a similar conformational change does not appear to occur for the ferrous, CO-ligated enzyme. To better understand the effector role of Pdx when binding the ferrous, CO-ligated P450cam, we applied 2D IR spectroscopy to compare the conformations and dynamics of the wild-type (wt) enzyme in the absence and presence of Pdx, as well as of L358P P450cam (L358P), which has served as a putative model for the Pdx complex. The CO vibrations of the Pdx complex and L358P report population of two conformational states in which the CO experiences distinct environments. The dynamics among the CO frequencies indicate that the energy landscape of substates within one conformation are reflective of the closed state of P450cam, and for the other conformation, differ from the free wt enzyme, but are equivalent between the Pdx complex and L358P. The two states co-populated by the Pdx complex are postulated to reflect a loosely bound encounter complex and a more tightly bound state, as is commonly observed for the dynamic complexes of redox partners. Significantly, this study shows that the binding of Pdx to ferrous, CO-ligated P450cam does perturb the conformational ensemble in a way that might underlie the effector role of Pdx.

Molecular Dynamics Simulation for All

The impact of molecular dynamics (MD) simulations in molecular biology and drug discovery has expanded dramatically in recent years. These simulations capture the behavior of proteins and other biomolecules in full atomic detail and at very fine temporal resolution. Major improvements in simulation speed, accuracy, and accessibility, together with the proliferation of experimental structural data, have increased the appeal of biomolecular simulation to experimentalists-a trend particularly noticeable in, although certainly not limited to, neuroscience. Simulations have proven valuable in deciphering functional mechanisms of proteins and other biomolecules, in uncovering the structural basis for disease, and in the design and optimization of small molecules, peptides, and proteins. Here we describe, in practical terms, the types of information MD simulations can provide and the ways in which they typically motivate further experimental work.

Effect of redox partner binding on CYP101D1 conformational dynamics

We have compared the thermodynamics of substrate and redox partner binding of P450cam to its close homologue, CYP101D1, using isothermal titration calorimetry (ITC). CYP101D1 binds camphor about 10-fold more weakly than P450cam which is consistent with the inability of camphor to cause a complete low- to high-spin shift in CYP101D1. Even so molecular dynamics simulations show that camphor is very stable in the CYP101D1 active site similar to P450cam. ITC data on the binding of the CYP101D1 ferredoxin redox partner (abbreviated Arx) shows that the substrate-bound closed state of CYP101D1 binds Arx more tightly than the substrate-free open form. This is just the opposite to P450cam where Pdx (ferredoxin redox partner of P450cam) favors binding to the P450cam open state. In addition, CYP101D1-Arx binding has a large negative ΔS while the P450cam-Pdx has a much smaller ΔS indicating that interactions at the docking interface are different. The most obvious difference is that PDXD38 which forms an important ion pair with P450camR112 at the center of the interface is ArxL39 in Arx. This suggests that Arx may adopt a different orientation than Pdx in order to optimize nonpolar interactions with ArxL39.

Ligand Access Channels in Cytochrome P450 Enzymes: A Review

Quantitative structure-activity relationships may bring invaluable information on structural elements of both enzymes and substrates that, together, govern substrate specificity. Buried active sites in cytochrome P450 enzymes are connected to the solvent by a network of channels exiting at the distal surface of the protein. This review presents different in silico tools that were developed to uncover such channels in P450 crystal structures. It also lists some of the experimental evidence that actually suggest that these predicted channels might indeed play a critical role in modulating P450 functions. Amino acid residues at the entrance of the channels may participate to a first global ligand recognition of ligands by P450 enzymes before they reach the buried active site. Moreover, different P450 enzymes show different networks of predicted channels. The plasticity of P450 structures is also important to take into account when looking at how channels might play their role.

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.