Putidaredoxin Binds to the Same Site on Cytochrome P450cam in the Open and Closed Conformation

Applications of Molecular Dynamics Simulations in Drug Discovery

In the current drug development process, molecular dynamics (MD) simulations have proven to be very useful. This chapter provides an overview of the current applications of MD simulations in drug discovery, from detecting protein druggable sites and validating drug docking outcomes to exploring protein conformations and investigating the influence of mutations on its structure and functions. In addition, this chapter emphasizes various strategies to improve the conformational sampling efficiency in molecular dynamics simulations. With a growing computer power and developments in the production of force fields and MD techniques, the importance of MD simulations in helping the drug development process is projected to rise significantly in the future.

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.

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.

Kinetic Evidence for an Induced Fit Mechanism in the Binding of the Substrate Camphor by Cytochrome P450cam

Bacterial cytochrome P450 (P450) 101A1 (P450cam) has served as a prototype among the P450 enzymes and has high catalytic activity towards its cognate substrate, camphor. X-ray crystallography and NMR and IR spectroscopy have demonstrated the existence of multiple conformations of many P450s, including P450cam. Kinetic studies have indicated that substrate binding to several P450s is dominated by a conformational selection process, in which the substrate binds an individual conformer(s) of the unliganded enzyme. P450cam was found to differ in that binding of the substrate camphor is dominated by an induced fit mechanism, in which the enzyme binds camphor and then changes conformation, as evidenced by the equivalence of binding eigenvalues observed when varying both camphor and P450cam concentrations. The accessory protein putidaredoxin had no effect on substrate binding. Estimation of the rate of dissociation of the P450cam·camphor complex (15 s-1) and fitting of the data yield a minimal kinetic mechanism in which camphor binds (1.5 × 107 M-1 s-1) and the initial P450cam•camphor complex undergoes a reversible equilibrium (k forward 112 s-1, k reverse 28 s-1) to a final complex. This induced fit mechanism differs from those reported for several mammalian P450s and bacterial P450BM-3, indicative of the diversity of how P450s recognize multiple substrates. However, similar behavior was not observed with the alternate substrates (+)-α-pinene and 2-adamantanone, which probably utilize a conformational selection process.

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.

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.

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.

Mechanisms of Cytochrome P450-Catalyzed Oxidations

Enzymes are complex biological catalysts and are critical to life. Most oxidations of chemicals are catalyzed by cytochrome P450 (P450, CYP) enzymes, which generally utilize mixed-function oxidase stoichiometry, utilizing pyridine nucleotides as electron donors: NAD(P)H + O₂ + R → NAD(P)⁺ + RO + H₂O (where R is a carbon substrate and RO is an oxidized product). The catalysis of oxidations is largely understood in the context of the heme iron-oxygen complex generally referred to as Compound I, formally FeO³⁺, whose basis was in peroxidase chemistry. Many X-ray crystal structures of P450s are now available (≥ 822 structures from ≥146 different P450s) and have helped in understanding catalytic specificity. In addition to hydroxylations, P450s catalyze more complex oxidations, including C-C bond formation and cleavage. Enzymes derived from P450s by directed evolution can even catalyze more unusual reactions, e.g. cyclopropanation. Current P450 questions under investigation include the potential role of the intermediate Compound 0 (formally Fe^III-O₂ ^-) in catalysis of some reactions, the roles of high- and low-spin forms of Compound I, the mechanism of desaturation, the roles of open and closed structures of P450s in catalysis, the extent of processivity in multi-step oxidations, and the role of the accessory protein cytochrome b ₅. More global questions include exactly how structure drives function, prediction of catalysis, and roles of multiple protein conformations.

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.

Conformational Changes in Cytochrome P450cam and the Effector Role of Putidaredoxin

The cytochromes P450 form an enormous family of over 20 000 enzyme variants found in all branches of life. They catalyze the O2 dependent monooxygenation of a wide range of substrates in reactions important to drug metabolism, biosynthesis and energy utilization. Understanding how they function is important for biomedical science and requires a full description of their notorious propensity for specificity and promiscuity. The bacterial P450cam is an unusual example, having the most well characterized chemical mechanism of all of the forms. It also undergoes an increasingly well characterized structural change upon substrate binding, which may be similar to to that displayed by some, but not all forms of P450. Finally, P450cam is one of the rare forms that have a strict requirement for a particular electron donor, putidaredoxin (pdx). Pdx provides the required electrons for enzyme turnover, but it also induces specific changes in the enzyme to allow enzyme turnover, long known as its effector role. This review summarizes recent crystallographic and double electron–electron resonance studies that have revealed the effects of substrate and pdx binding on the structure of P450cam. We describe an emerging idea for how pdx exerts its effector function by inducing a conformational change in the enzyme. This change then propagates to the active site to enable cleavage of the ferric–hydroperoxy bond during catalysis, and appears to provide a very elegant approach for P450cam to attain both high efficiency and protection from oxidative damage.

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.

A Structural Model of a P450-Ferredoxin Complex from Orientation-Selective Double Electron-Electron Resonance Spectroscopy

Cytochrome P450 (CYP) monooxygenases catalyze the oxidation of chemically inert carbon-hydrogen bonds in diverse endogenous and exogenous organic compounds by atmospheric oxygen. This C-H bond oxy-functionalization activity has huge potential in biotechnological applications. Class I CYPs receive the two electrons required for oxygen activation from NAD(P)H via a ferredoxin reductase and ferredoxin. The interaction of Class I CYPs with their cognate ferredoxin is specific. In order to reconstitute the activity of diverse CYPs, structural characterization of CYP-ferredoxin complexes is necessary, but little structural information is available. Here we report a structural model of such a complex (CYP199A2-HaPux) in frozen solution derived from distance and orientation restraints gathered by the EPR technique of orientation-selective double electron-electron resonance (os-DEER). The long-lived oscillations in the os-DEER spectra were well modeled by a single orientation of the CYP199A2-HaPux complex. The structure is different from the two known Class I CYP-Fdx structures: CYP11A1-Adx and CYP101A1-Pdx. At the protein interface, HaPux residues in the [Fe₂S₂] cluster-binding loop and the α3 helix and the C-terminus residue interact with CYP199A2 residues in the proximal loop and the C helix. These residue contacts are consistent with biochemical data on CYP199A2-ferredoxin binding and electron transfer. Electron-tunneling calculations indicate an efficient electron-transfer pathway from the [Fe₂S₂] cluster to the heme. This new structural model of a CYP-Fdx complex provides the basis for tailoring CYP enzymes for which the cognate ferredoxin is not known, to accept electrons from HaPux and display monooxygenase activity.

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

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

A Distal Loop Controls Product Release and Chemo- and Regioselectivity in Cytochrome P450 Decarboxylases

Cytochrome P450 OleT utilizes hydrogen peroxide (H₂O₂) to catalyze the decarboxylation or hydroxylation of fatty acid (FA) substrates. Both reactions are initiated through the abstraction of a substrate hydrogen atom by the high-valent iron-oxo intermediate known as Compound I. Here, we specifically probe the influence of substrate coordination on OleT reaction partitioning through the combined use of fluorescent and electron paramagnetic resonance (EPR)-active FA probes and mutagenesis of a structurally disordered F-G loop that is distal from the heme-iron active site. Both probes are efficiently metabolized by OleT and efficiently trigger the formation of Compound I. Transient fluorescence and EPR reveal a slow product release step, mediated by the F-G loop, that limits OleT turnover. A single-amino acid change or excision of the loop reveals that this region establishes critical interactions to anchor FA substrates in place. The stabilization afforded by the F-G loop is essential for regulating regiospecific C-H abstraction and allowing for efficient decarboxylation to occur. These results highlight a regulatory strategy whereby the fate of activated oxygen species can be controlled at distances far removed from the site of chemistry.