Delicate conformational balance of the redox enzyme cytochrome P450cam

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.

Paramagnetic Chemical Probes for Studying Biological Macromolecules

Paramagnetic chemical probes have been used in electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) spectroscopy for more than four decades. Recent years witnessed a great increase in the variety of probes for the study of biological macromolecules (proteins, nucleic acids, and oligosaccharides). This Review aims to provide a comprehensive overview of the existing paramagnetic chemical probes, including chemical synthetic approaches, functional properties, and selected applications. Recent developments have seen, in particular, a rapid expansion of the range of lanthanoid probes with anisotropic magnetic susceptibilities for the generation of structural restraints based on residual dipolar couplings and pseudocontact shifts in solution and solid state NMR spectroscopy, mostly for protein studies. Also many new isotropic paramagnetic probes, suitable for NMR measurements of paramagnetic relaxation enhancements, as well as EPR spectroscopic studies (in particular double resonance techniques) have been developed and employed to investigate biological macromolecules. Notwithstanding the large number of reported probes, only few have found broad application and further development of probes for dedicated applications is foreseen.

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.

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.

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.

The Photocatalyzed Thiol-ene reaction: A New Tag to Yield Fast, Selective and reversible Paramagnetic Tagging of Proteins

Paramagnetic restraints have been used in biomolecular NMR for the last three decades to elucidate and refine biomolecular structures, but also to characterize protein-ligand interactions. A common technique to generate such restraints in proteins, which do not naturally contain a (paramagnetic) metal, consists in the attachment to the protein of a lanthanide-binding-tag (LBT). In order to design such LBTs, it is important to consider the efficiency and stability of the conjugation, the geometry of the complex (conformational exchanges and coordination) and the chemical inertness of the ligand. Here we describe a photo-catalyzed thiol-ene reaction for the cysteine-selective paramagnetic tagging of proteins. As a model, we designed an LBT with a vinyl-pyridine moiety which was used to attach our tag to the protein GB1 in fast and irreversible fashion. Our tag T1 yields magnetic susceptibility tensors of significant size with different lanthanides and has been characterized using NMR and relaxometry measurements.

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.

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.

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.

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.

Amino Acid Selective 13C Labeling and 13C Scrambling Profile Analysis of Protein α and Side-Chain Carbons in Escherichia coli Utilized for Protein Nuclear Magnetic Resonance

Amino acid selective isotope labeling is an important nuclear magnetic resonance technique, especially for larger proteins, providing strong bases for the unambiguous resonance assignments and information concerning the structure, dynamics, and intermolecular interactions. Amino acid selective ¹⁵N labeling suffers from isotope dilution caused by metabolic interconversion of the amino acids, resulting in isotope scrambling within the target protein. Carbonyl ¹³C atoms experience less isotope scrambling than the main-chain ¹⁵N atoms do. However, little is known about the side-chain ¹³C atoms. Here, the ¹³C scrambling profiles of the Cα and side-chain carbons were investigated for ¹⁵N scrambling-prone amino acids, such as Leu, Ile, Tyr, Phe, Thr, Val, and Ala. The level of isotope scrambling was substantially lower in ¹³Cα and ¹³C side-chain labeling than in ¹⁵N labeling. We utilized this reduced scrambling-prone character of ¹³C as a simple and efficient method for amino acid selective ¹³C labeling using an Escherichia coli cold-shock expression system and high-cell density fermentation. Using this method, the ¹³C labeling efficiency was >80% for Leu and Ile, ∼60% for Tyr and Phe, ∼50% for Thr, ∼40% for Val, and 30-40% for Ala. ¹H-¹⁵N heteronuclear single-quantum coherence signals of the ¹⁵N scrambling-prone amino acid were also easily filtered using ¹⁵N-{¹³Cα} spin-echo difference experiments. Our method could be applied to the assignment of the 55 kDa protein.

Thermodynamics of camphor migration in cytochrome P450cam by atomistic simulations

Understanding the mechanisms of ligand binding to enzymes is of paramount importance for the design of new drugs. Here, we report on the use of a novel biased molecular dynamics (MD) methodology to study the mechanism of camphor binding to cytochrome P450cam. Microsecond-long MD simulations allowed us to observe reaction coordinates characterizing ligand diffusion from the active site of cytochrome P450cam to solvent via three egress routes. These atomistic simulations were used to estimate thermodynamic quantities along the reaction coordinates and indicate diverse binding configurations. The results suggest that the diffusion of camphor along the pathway near the substrate recognition site (SRS) is thermodynamically preferred. In addition, we show that the diffusion near the SRS is triggered by a transition from a heterogeneous collection of closed ligand-bound conformers to the basin comprising the open conformations of cytochrome P450cam. The conformational change accompanying this switch is characterized by the retraction of the F and G helices and the disorder of the B' helix. These results are corroborated by experimental studies and provide detailed insight into ligand binding and conformational behavior of the cytochrome family. The presented methodology is general and can be applied to other ligand-protein systems.

Structure restraints from heteronuclear pseudocontact shifts generated by lanthanide tags at two different sites

Pseudocontact shifts (PCS) encode long-range information on 3D structures of protein backbones and side-chains. The level of structural detail that can be obtained increases with the number of different sites tagged with a paramagnetic metal ion to generate PCSs. Here we show that PCSs from two different sites can suffice to determine the structure of polypeptide chains and their location and orientation relative to the magnetic susceptibility tensor χ, provided that PCSs are available for ¹H as well as heteronuclear spins. In addition, PCSs from two different sites are shown to provide detailed structural information on the conformation of methyl group-bearing amino-acid side-chains. A previously published ensemble structure of ubiquitin is shown to explain the magnetic susceptibility and alignment tensors slightly better than structures that try to explain the experimental data by a single conformation, illustrating the potential of PCSs as a tool to investigate small conformational changes.

Identification of productive and futile encounters in an electron transfer protein complex

Well-defined, stereospecific states in protein complexes are often in exchange with an ensemble of more dynamic orientations: the encounter states. The structure of the stereospecific complex between cytochrome P450cam and putidaredoxin was solved recently by X-ray diffraction as well as paramagnetic NMR spectroscopy. Other than the stereospecific complex, the NMR data clearly show the presence of additional states in the complex in solution. In these encounter states, populated for a small percentage of the time, putidaredoxin assumes multiple orientations and samples a large part of the surface of cytochrome P450cam. To characterize the nature of the encounter states, an extensive paramagnetic NMR dataset has been analyzed using the Maximum Occurrence of Regions methodology. The analysis reveals the location and maximal spatial extent of the additional states needed to fully explain the NMR data. Under the assumption of sparsity of the size of the conformational ensemble, several minor states can be located quite precisely. The distribution of these minor states correlates with the electrostatic potential map around cytochrome P450cam. Whereas some minor states are on isolated positively charged patches, others are connected to the stereospecific site via positively charged paths. The existence of electrostatically favorable pathways between the stereospecific interaction site and the different minor states or lack thereof suggests a means to discriminate between productive and futile encounter states.

Effector Roles of Putidaredoxin on Cytochrome P450cam Conformational States

In this study, the effector role of Pdx (putidaredoxin) on cytochrome P450cam conformation is refined by attaching two different spin labels, MTSL or BSL (bifunctional spin-label) onto the F or G helices and using DEER (double electron-electron resonance) to measure the distance between labels. Recent EPR and crystallographic studies have observed that oxidized Pdx induces substrate-bound P450cam to change from the closed to the open state. However, this change was not observed by DEER in the reduced Pdx complex with carbon-monoxide-bound P450cam (Fe(2+)CO). In addition, recent NMR studies have failed to observe a change in P450cam conformation upon binding Pdx. Hence, resolving these issues is important for a full understanding the effector role of Pdx. Here we show that oxidized Pdx induces camphor-bound P450cam to shift from the closed to the open conformation when labeled on either the F or G helices with MTSL. BSL at these sites can either narrow the distance distribution widths dramatically or alter the extent of the conformational change. In addition, we report DEER spectra on a mixed oxidation state containing oxidized Pdx and ferrous CO-bound P450cam, showing that P450cam remains closed. This indicates that CO binding to the heme prevents P450cam from opening, overriding the influence exerted by Pdx binding. Finally, we report the open form P450cam crystal structure with substrate bound, which suggests that crystal packing effects may prevent conformational conversion. Using multiple labeling approaches, DEER provides a unique perspective to resolve how the conformation of P450cam depends on Pdx and ligand states.

Conformational selectivity in cytochrome P450 redox partner interactions

The heme iron of cytochromes P450 must be reduced to bind and activate molecular oxygen for substrate oxidation. Reducing equivalents are derived from a redox partner, which requires the formation of a protein-protein complex. A subject of increasing discussion is the role that redox partner binding plays, if any, in favoring significant structural changes in the P450s that are required for activity. Many P450s now have been shown to experience large open and closed motions. Several structural and spectral studies indicate that the well-studied P450cam adopts the open conformation when its redox partner, putidaredoxin (Pdx), binds, whereas recent NMR studies indicate that this view is incorrect. Given the relevance of this discrepancy to P450 chemistry, it is important to determine whether Pdx favors the open or closed form of P450cam. Here, we have used both computational and experimental isothermal titration calorimetry studies that unequivocally show Pdx favors binding to the open form of P450cam. Analyses of molecular-dynamic trajectories also provide insights into intermediate conformational states that could be relevant to catalysis.

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.

Structural Study of a Flexible Active Site Loop in Human Indoleamine 2,3-Dioxygenase and Its Functional Implications

Human indoleamine 2,3-dioxygenase catalyzes the oxidative cleavage of tryptophan to N-formyl kynurenine, the initial and rate-limiting step in the kynurenine pathway. Additionally, this enzyme has been identified as a possible target for cancer therapy. A 20-amino acid protein segment (the JK loop), which connects the J and K helices, was not resolved in the reported hIDO crystal structure. Previous studies have shown that this loop undergoes structural rearrangement upon substrate binding. In this work, we apply a combination of replica exchange molecular dynamics simulations and site-directed mutagenesis experiments to characterize the structure and dynamics of this protein region. Our simulations show that the JK loop can be divided into two regions: the first region (JK loop(C)) displays specific and well-defined conformations and is within hydrogen bonding distance of the substrate, while the second region (JK loop(N)) is highly disordered and exposed to the solvent. The peculiar flexible nature of JK loop(N) suggests that it may function as a target for post-translational modifications and/or a mediator for protein-protein interactions. In contrast, hydrogen bonding interactions are observed between the substrate and Thr379 in the highly conserved "GTGG" motif of JK loop(C), thereby anchoring JK loop(C) in a closed conformation, which secures the appropriate substrate binding mode for catalysis. Site-directed mutagenesis experiments confirm the key role of this residue, highlighting the importance of the JK loop(C) conformation in regulating the enzymatic activity. Furthermore, the existence of the partially and totally open conformations in the substrate-free form suggests a role of JK loop(C) in controlling substrate and product dynamics.

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

Cytochrome P450 monooxygenases typically catalyze the insertion of one atom of oxygen from O2 into unactivated carbon-hydrogen and carbon-carbon bonds, with concomitant reduction of the other oxygen atom to H2O by NAD(P)H. Comparison of the average structures of the camphor hydroxylase cytochrome P450(cam) (CYP101) obtained from residual dipolar coupling (RDC)-restrained molecular dynamics (MD) in the presence and absence of substrate camphor shows structural displacements resulting from the essential collapse of the active site upon substrate removal. This collapse has conformational consequences that extend across the protein structure, none of which were observed in analogous crystallographic structures. Mutations were made to test the involvement of the observed conformational changes in substrate binding and recognition. All of the mutations performed based upon the NMR-detected perturbations, even those remote from the active site, resulted in modified substrate selectivity, enzyme efficiency and/or haem iron spin state. The results demonstrate that solution NMR can provide insights into enzyme structure-function relationships that are difficult to obtain by other methods.

Determination of pseudocontact shifts of low-populated excited states by NMR chemical exchange saturation transfer

Angular and distance restraints for low populated excited conformations are studied using PCS–CEST NMR spectroscopy.

An encodable lanthanide binding tag with reduced size and flexibility for measuring residual dipolar couplings and pseudocontact shifts in large proteins

Metal ions serve important roles in structural biology applications from long-range perturbations seen in magnetic resonance experiments to electron-dense signatures in X-ray crystallography data; however, the metal ion must be secured in a molecular framework to achieve the maximum benefit. Polypeptide-based lanthanide-binding tags (LBTs) represent one option that can be directly encoded within a recombinant protein expression construct. However, LBTs often exhibit significant mobility relative to the target molecule. Here we report the characterization of improved LBTs sequences for insertion into a protein loop. These LBTs were inserted to connect two parallel alpha helices of an immunoglobulin G (IgG)-binding Z domain platform. Variants A and B bound Tb(3+) with high affinity (0.70 and 0.13 μM, respectively) and displayed restricted LBT motion. Compared to the parent construct, the metal-bound A experienced a 2.5-fold reduction in tag motion as measured by magnetic field-induced residual dipolar couplings and was further studied in a 72.2 kDa complex with the human IgG1 fragment crystallizable (IgG1 Fc) glycoprotein. The appearance of both pseudo-contact shifts (-0.221 to 0.081 ppm) and residual dipolar couplings (-7.6 to 14.3 Hz) of IgG1 Fc resonances in the IgG1 Fc:(variant A:Tb(3+))2 complex indicated structural restriction of the LBT with respect to the Fc. These studies highlight the applicability of improved LBT sequences with reduced mobility to probe the structure of macromolecular systems.

Enantiomeric two-armed lanthanide-binding tags for complementary effects in paramagnetic NMR spectroscopy

A new pair of two-armed lanthanide-binding tags provide distinct sets of structural restraints when attached to the same site of a protein.

Single-armed phenylsulfonated pyridine derivative of DOTA is rigid and stable paramagnetic tag in protein analysis

Single-armed DOTA-like phenylsulfonated pyridine derivatives are rigid and stable paramagnetic tags for site-specific labelling of proteins. The respective protein conjugates yield valuable long-range structural restraints for proteins.

Pulse EPR-enabled interpretation of scarce pseudocontact shifts induced by lanthanide binding tags

Pseudocontact shifts (PCS) induced by tags loaded with paramagnetic lanthanide ions provide powerful long-range structure information, provided the location of the metal ion relative to the target protein is known. Usually, the metal position is determined by fitting the magnetic susceptibility anisotropy (Δχ) tensor to the 3D structure of the protein in an 8-parameter fit, which requires a large set of PCSs to be reliable. In an alternative approach, we used multiple Gd(3+)-Gd(3+) distances measured by double electron-electron resonance (DEER) experiments to define the metal position, allowing Δχ-tensor determinations from more robust 5-parameter fits that can be performed with a relatively sparse set of PCSs. Using this approach with the 32 kDa E. coli aspartate/glutamate binding protein (DEBP), we demonstrate a structural transition between substrate-bound and substrate-free DEBP, supported by PCSs generated by C3-Tm(3+) and C3-Tb(3+) tags attached to a genetically encoded p-azidophenylalanine residue. The significance of small PCSs was magnified by considering the difference between the chemical shifts measured with Tb(3+) and Tm(3+) rather than involving a diamagnetic reference. The integrative sparse data approach developed in this work makes poorly soluble proteins of limited stability amenable to structural studies in solution, without having to rely on cysteine mutations for tag attachment.

Quantum Mechanics/Molecular Mechanics Modeling of Enzymatic Processes: Caveats and Breakthroughs

Nature has developed large groups of enzymatic catalysts with the aim to transfer substrates into useful products, which enables biosystems to perform all their natural functions. As such, all biochemical processes in our body (we drink, we eat, we breath, we sleep, etc.) are governed by enzymes. One of the problems associated with research on biocatalysts is that they react so fast that details of their reaction mechanisms cannot be obtained with experimental work. In recent years, major advances in computational hardware and software have been made and now large (bio)chemical systems can be studied using accurate computational techniques. One such technique is the quantum mechanics/molecular mechanics (QM/MM) technique, which has gained major momentum in recent years. Unfortunately, it is not a black-box method that is easily applied, but requires careful set-up procedures. In this work we give an overview on the technical difficulties and caveats of QM/MM and discuss work-protocols developed in our groups for running successful QM/MM calculations.