Conformational selectivity in cytochrome P450 redox partner interactions

Molecular dynamics simulations assisted investigation of phytochemicals as potential lead candidates against anti-apoptotic Bcl-B protein

Due to the multifarious nature of cancer, finding a single definitive cure for this dreadful disease remains an elusive challenge. The dysregulation of the apoptotic pathway or programmed cell death, governed by the Bcl-2 family of proteins plays a crucial role in cancer development and progression. Bcl-B stands out as a unique anti-apoptotic protein from the Bcl-2 family that selectively binds to Bax which inhibits its pro-apoptotic function. Although several inhibitors are reported for Bcl-2 family proteins, no specific inhibitors are available against the anti-apoptotic Bcl-B protein. This study aims to address this research gap by using virtual screening of an in-house library of phytochemicals from seven anti-cancer medicinal plants to identify lead molecules against Bcl-B protein. Through pharmacokinetic analysis and molecular docking studies, we identified three lead candidates (Enterolactone, Piperine, and Protopine) based on appreciable drug-likeliness, ADME properties, and binding affinity values. The identified molecules also exhibited specific interactions with critical amino acid residues of the binding cleft, highlighting their potential as lead candidates. Finally, molecular dynamics simulations and MM/PBSA based binding free energy analysis revealed that Enterolactone (CID_114739) and Piperine (CID_638024) molecules were on par with Obatoclax (CID_11404337), which is a known inhibitor of the Bcl-2 family proteins.Communicated by Ramaswamy H. Sarma.

Interactions of human mitochondrial Ferredoxin 1 (Adrenodoxin) by NMR; modulation by cytochrome P450 substrate and by truncation of the C-terminal tail

Human Ferredoxin 1, also referred to as Adrenodoxin (Adx), is the sole electron carrier supporting the function of all seven mitochondrial cytochrome P450 (CYP) enzymes. Adx utilizes conserved negatively charged residues along its α-helix3 to interact with either the proximal surface of CYP enzymes or the binding surface of Adrendodoxin Reductase (AdR). However, in the oxidized state, Adx assumes a monomer-homodimer equilibrium that requires the presence of its unstructured C-terminal tail. Crystallographic structures of full-length human Adx dimers indicate that part of the binding surface necessary for its interactions with CYPs or with AdR is partially occluded by the dimer interface. In this study, protein NMR spectroscopy was used to interrogate the interactions between full-length (2-124) or truncated monomeric (2-108) human Adx and human CYP24A1 (with and without its vitamin-D substrate) as well as interactions with AdR. Here, monomeric Adx induced a similar pattern of peak broadening as that induced by addition of CYP24A1 substrate, consistent with a 1:1 Adx:CYP interaction as the functional complex. Additionally, removal of the C-terminal tail appears to enhance the interaction with AdR, despite removal of some of the AdR contacts in the tail region. This finding was also supported by an NMR competition assay. These findings suggest that the Adx dimers do not undergo meaningful interactions with either CYP or AdR, but may instead be responsible for regulating access to monomeric Adx. These conclusions are discussed in the context of a revised model of the Adx electron shuttle mechanism.

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.

Structural and functional characterization of 1-deoxy-D-xylulose-5-phosphate synthase (DXS) from Acinetobacter baumannii: identification of promising lead molecules from virtual screening, molecular docking and molecular dynamics simulations

The advent of multi drug resistance and extensive-drug resistance among various pathogens has caused a rise in nosocomial infections, which in turn has led to rising hospital-acquired infection-related mortality rates. Amongst them, carbapenem-resistant Acinetobacter baumannii is one of the most notorious bacterial species, categorized as a Priority 1: Critical pathogen by the WHO. Therefore, the discovery and development of novel antibiotics, as well as the identification of potential inhibitors, have become the need of the hour. In this study, we have employed computational methods to explore and identify molecules capable of inhibiting enzymes essential in the methylerythritol 4-phosphate (MEP) biosynthetic pathway. The high throughput virtual screening of small molecules (Enamine Advanced Collection (AC) library) against the highly conserved substrate-binding site of the DXS target protein provided us with a total of 1000 molecules. The top four potential candidate molecules, namely-Z3353989070, Z3353989049, Z2295848528, and Z1685501455, alongside fluoropyruvate (control), a known inhibitor of DXS, was chosen for a molecular dynamic simulation study. The molecular dynamic simulation trajectories suggested high structural and thermodynamical stability and strong binding affinity of all the DXS-ligand complexes. Moreover, the MM/PBSA-based binding free energy calculations also exhibited strong interactions of the selected ligand molecules with DXS. In conclusion, we have found that all four molecules displayed better results and stronger binding affinity than the control. In the end, based on all the above-mentioned criteria, we have proposed Z3353989049 to be the promising lead candidate against DXS from A. baumannii.Communicated by Ramaswamy H. Sarma.

Three pairs of surrogate redox partners comparison for Class I cytochrome P450 enzyme activity reconstitution

Most P450s require redox partners for the electron transfer during catalysis. However, little information is available on cognate redox partners for P450s, which greatly limits P450 function exploration and practical application. Thus, the stategy of building various hybrid P450 catalytic systems with surrogate redox partner has often adopted to engineer P450 biocatalysts. In this study, we compare three pairs of frequently-used surrogate redox partner SelFdx1499/SelFdR0978, Adx/AdR and Pdx/PdR and in terms of their electron transfer properties. The three selected bacterial Class I P450s include PikC, P450sca-2 and CYP-sb21, which are responsible for production of high-value-added products. Here we show that SelFdx1499/SelFdR0978 is the most promising redox partner compared to Adx/AdR and Pdx/PdR. The results provide insights into the domination for P450-redox partner interactions in modulating the catalytic activity of P450s. This study not only produces a more active biocatalyst but also suggests a general chose for a universal reductase which would facilitate engineering of P450 catalyst.

Aiming for an efficient Class I cytochrome P450 catalytic system, three pairs of surrogate redox partners for biocatalyst applications are tested.

A computational essential dynamics approach to investigate structural influences of ligand binding on Papain like protease from SARS-CoV-2

Papain like protease (PLpro) is a cysteine protease from the coronaviridae family of viruses. Coronaviruses possess a positive sense, single-strand RNA, leading to the translation of two viral polypeptides containing viral structural, non-structural and accessory proteins. PLpro is responsible for the cleavage of nsp1–3 from the viral polypeptide. PLpro also possesses deubiquitinating and deISGlyating activity, which sequesters the virus from the host's immune system. This indispensable attribute of PLpro makes it a protein of interest as a drug target. The present study aims to analyze the structural influences of ligand binding on PLpro. First, PLpro was screened against the ZINC-in-trials library, from which four lead compounds were identified based on estimated binding affinity and interaction patterns. Next, based on molecular docking results, ZINC000000596945, ZINC000064033452 and VIR251 (control molecule) were subjected to molecular dynamics simulation. The study evaluated global and essential dynamics analyses utilising principal component analyses, dynamic cross-correlation matrix, free energy landscape and time-dependant essential dynamics to predict the structural changes observed in PLpro upon ligand binding in a simulated environment. The MM/PBSA-based binding free energy calculations of the two selected molecules, ZINC000000596945 (−41.23 ± 3.70 kcal/mol) and ZINC000064033452 (−25.10 ± 2.65 kcal/mol), displayed significant values which delineate them as potential inhibitors of PLpro from SARS-CoV-2.

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.

Structural Dynamics of Cytochrome P450 3A4 in the Presence of Substrates and Cytochrome P450 Reductase

Cytochrome P450 3A4 (CYP3A4) is the most important drug-metabolizing enzyme in humans and has been associated with harmful drug interactions. The activity of CYP3A4 is known to be modulated by several compounds and by the electron transfer partner, cytochrome P450 reductase (CPR). The underlying mechanism of these effects, however, is poorly understood. We have used hydrogen-deuterium exchange mass spectrometry to investigate the impact of binding of CPR and of three different substrates (7-benzyloxy-4-trifluoromethyl-coumarin, testosterone, and progesterone) on the conformational dynamics of CYP3A4. Here, we report that interaction of CYP3A4 with substrates or with the oxidized or reduced forms of CPR leads to a global rigidification of the CYP3A4 structure. This was evident from the suppression of deuterium exchange in several regions of CYP3A4, including regions known to be involved in protein-protein interactions (helix C) and substrate binding and specificity (helices B' and E, and loop K/β1). Furthermore, the bimodal isotopic distributions observed for some CYP3A4-derived peptides were drastically impacted upon binding to CPR and/or substrates, suggesting the existence of stable CYP3A4 conformational populations that are perturbed by ligand/CPR binding. The results have implications for understanding the mechanisms of ligand binding, allostery, and catalysis in CYP enzymes.

Self-Sufficient Class VII Cytochromes P450: From Full-Length Structure to Synthetic Biology Applications

Members of class VII cytochromes P450 are catalytically self-sufficient enzymes containing a phthalate dioxygenase reductase-like domain fused to the P450 catalytic domain. Among these, CYP116B46 is the first enzyme for which the 3D structure of the whole polypeptide chain has been solved, shedding light on the interaction between its domains, which is crucial for catalysis. Most of these enzymes have been isolated from extremophiles or detoxifying bacteria that can carry out regio- and enantioselective oxidation of compounds of biotechnological interest. Protein engineering has generated mutants that can perform challenging organic reactions such as the anti-Markovnikov alkene oxidation. This potential, combined with the detailed 3D structure, forms the basis for further directed evolution studies aimed at widening their biotechnological exploitation.

A new twist of rubredoxin function in M. tuberculosis

Electron transfer mediated by metalloproteins drives many biological processes. Rubredoxins are a ubiquitous [1Fe-0S] class of electron carriers that play an important role in bacterial adaptation to changing environmental conditions. In Mycobacterium tuberculosis, oxidative and acidic stresses as well as iron starvation induce rubredoxins expression. However, their functions during M. tuberculosis infection are unknown. In the present work, we show that rubredoxin B (RubB) is able to efficiently shuttle electrons from cognate reductases, FprA and FdR to support catalytic activity of cytochrome P450s, CYP124, CYP125, and CYP142, which are important for bacterial viability and pathogenicity. We solved the crystal structure of RubB and characterized the interaction between RubB and CYPs using site-directed mutagenesis. Mutations that not only neutralize single charge but also change the specific residues on the surface of RubB did not dramatically decrease activity of studied CYPs. Together with isothermal calorimetry (ITC) experiments, the obtained results suggest that interactions are transient and not highly specific. The redox potential of RubB is -264 mV vs. Ag/AgCl and the measured extinction coefficients are 9931 M^-1cm^-1 and 8371 M^-1cm^-1 at 380 nm and 490 nm, respectively. Characteristic parameters of RubB along with the discovered function might be useful for biotechnological applications. Our findings suggest that a switch from ferredoxins to rubredoxins might be crucial for M. tuberculosis to support CYPs activity during the infection.

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.

Molecular and Structural Evolution of Cytochrome P450 Aromatase

Aromatase is the cytochrome P450 enzyme converting androgens into estrogen in the last phase of steroidogenesis. As estrogens are crucial in reproductive biology, aromatase is found in vertebrates and the invertebrates of the genus Branchiostoma, where it carries out the aromatization reaction of the A-ring of androgens that produces estrogens. Here, we investigate the molecular evolution of this unique and highly substrate-selective enzyme by means of structural, sequence alignment, and homology modeling, shedding light on its key role in species conservation. The alignments led to the identification of a core structure that, together with key and unique amino acids located in the active site and the substrate recognition sites, has been well conserved during evolution. Structural analysis shows what their roles are and the reason why they have been preserved. Moreover, the residues involved in the interaction with the redox partner and some phosphorylation sites appeared late during evolution. These data reveal how highly substrate-selective cytochrome P450 has evolved, indicating that the driving forces for evolution have been the optimization of the interaction with the redox partner and the introduction of phosphorylation sites that give the possibility of modulating its activity in a rapid way.

Effector role of cytochrome P450 reductase for androstenedione binding to human aromatase

Cytochromes P450 constitute a large superfamily of monooxygenases involved in many metabolic pathways. Most of them are not self-sufficient and need a reductase protein to provide the electrons necessary for catalysis. It was shown that the redox partner plays a role in the modulation of the structure and function of some bacterial P450 enzymes. Here, the effect of NADPH-cytochrome reductase (CPR) on human aromatase (Aro) is studied for what concerns its role in substrate binding. Pre-steady-state kinetic experiments indicate that both the substrate binding rates and the percentage of spin shift detected for aromatase are increased when CPR is present. Moreover, aromatase binds the substrate through a conformational selection mechanism, suggesting a possible effector role of CPR. The thermodynamic parameters for the formation of the CPR-Aro complex were studied by isothermal titration calorimetry. The dissociation constant of the complex formation is 4.5 folds lower for substrate-free compared to the substrate-bound enzyme. The enthalpy change observed when the CPR-Aro complex forms in the absence of the substrate are higher than in its presence, indicating that more interactions are formed/broken in the former case. Taken together, our data confirm that CPR has a role in promoting aromatase conformation optimal for substrate binding.

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.

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.

Conformational turn triggers regio-selectivity in the bioactivation of thiophene-contained compounds mediated by cytochrome P450

In the present work, we performed Density Functional Theory calculations to explore the bioactivation mechanism of thiophene-containing molecules mediated by P450s. For this purpose, relatively large size compounds, 2,5-diaminothiophene derivatives were selected particularly for this investigation. Here we found the successive regio-selectivity triggered by conformational turn played a significant role in the occurrence of bioactivation. 2,5-Diaminothiophene was oxidized to a 2,5-diimine thiophene-reactive intermediate by Compound I (Cpd I) through successive activations of two N-H bonds (H3-N11 and H1-N6). This reaction exhibited three special characteristics: (1) self-controlled regio-selectivity during the oxidation process. There was a large scale of conformational turn in the abstraction of the first H atom which triggers the selection of the second H for abstraction. (2) Proton-shuttle mechanism. In high spin (HS) state, proton-shuttle mechanism was observed for the abstraction of the second H atom. (3) Spin-selective manner. In protein environment, the energy barrier in HS state was much lower than that in low spin state. The novel proposed bioactivation mechanism of 2,5-diaminothiophene compounds can help us in rational design of thiophene-contained drugs avoiding the occurrence of bioactivation.

Development of a Novel AOP for Cyp2F2-Mediated Lung Cancer in Mice

Traditional methods for carcinogenicity testing rely heavily on the rodent bioassay as the standard for identification of tumorigenic risk. As such, identification of species-specific outcomes and/or metabolism are a frequent argument for regulatory exemption. One example is the association of tumor formation in the mouse lung after exposure to Cyp2F2 ligands. The adverse outcome pathway (AOP) framework offers a theoretical platform to address issues of species specificity that is consistent, transparent, and capable of integrating data from new approach methodologies as well as traditional data streams. A central premise of the AOP concept is that pathway progression from the molecular initiating event (MIE) implies a definable “response-response” (R-R) relationship between each key event (KE) that drives the pathway towards a specific adverse outcome (AO). This article describes an AOP for lung cancer in the mouse from an MIE of Cyp2F2-specific reactive metabolite formation, advancing through KE that include protein and/or nucleic acid adducts, diminished Club Cell 10 kDa (CC10) protein expression, hyperplasia of CC10 deficient Club cells, and culminating in the AO of mixed-cell tumor formation in the distal airways. This tumor formation is independent of route of exposure and our AOP construct is based on overlapping mechanistic events for naphthalene, styrene, ethyl benzene, isoniazid, and fluensulfone in the mouse. This AOP is intended to accelerate the explication of an apparent mouse-specific outcome and serve as a starting point for a quantitative analysis of mouse-human differences in susceptibility to the tumorigenic effects of Cyp2F2 ligands.

Specificity of the Redox Complex between Cytochrome P450 24A1 and Adrenodoxin Relies on Carbon-25 Hydroxylation of Vitamin-D Substrate

Metabolic deactivation of 1,25(OH)₂D3 is initiated by modification of the vitamin-D side chain, as carried out by the mitochondrial cytochrome P450 24A1 (CYP24A1). In addition to its role in vitamin-D metabolism, CYP24A1 is involved in catabolism of vitamin-D analogs, thereby reducing their efficacy. CYP24A1 function relies on electron transfer from the soluble ferredoxin protein adrenodoxin (Adx). Recent structural evidence suggests that regioselectivity of the CYP24A1 reaction may correlate with distinct modes of Adx recognition. Here we used nuclear magnetic resonance (NMR) spectroscopy to monitor the structure of ¹⁵N-labeled full-length Adx from rat while forming the complex with rat CYP24A1 in the ligand-free state or bound to either 1,25(OH)₂D3 or the vitamin-D supplement 1α(OH)D3. Although both vitamin-D ligands were found to induce a reduction in overall NMR peak broadening, thereby suggesting ligand-induced disruption of the complex, a crosslinking analysis suggested that ligand does not have a significant effect on the relative association affinities of the redox complexes. However, a key finding is that, whereas the presence of primary CYP24A1 substrate was found to induce NMR peak broadening focused on the putative recognition site α-helix 3 of rat adrenodoxin, the interaction in the presence of 1α(OH)D3, which is lacking the carbon-25 hydroxyl, results in disruption of the NMR peak broadening pattern, thus indicating a ligand-induced nonspecific protein interaction. These findings provide a structural basis for the poor substrate turnover of side-chain-modified vitamin-D analogs, while also confirming that specificity of the CYP24A1-ligand interaction influences specificity of CYP24A1-Adx recognition. SIGNIFICANCE STATEMENT: Mitochondrial cytochrome P450 enzymes, such as CYP24A1 responsible for catabolizing vitamin-D and its analogs, rely on a protein-protein interaction with a ferredoxin in order to receive delivery of the electrons required for catalysis. In this study, we demonstrate that this protein interaction is influenced by the enzyme-ligand interaction that precedes it. Specifically, vitamin-D missing carbon-25 hydroxylation binds the enzyme active site with high affinity but results in a loss of P450-ferredoxin binding specificity.

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.

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.

Cytochrome-P450-Induced Ordering of Microsomal Membranes Modulates Affinity for Drugs

Although membrane environment is known to boost drug metabolism by mammalian cytochrome P450s, the factors that stabilize the structural folding and enhance protein function are unclear. In this study, we use peptide-based lipid nanodiscs to "trap" the lipid boundaries of microsomal cytochrome P450 2B4. We report the first evidence that CYP2B4 is able to induce the formation of raft domains in a biomimetic compound of the endoplasmic reticulum. NMR experiments were used to identify and quantitatively determine the lipids present in nanodiscs. A combination of biophysical experiments and molecular dynamics simulations revealed a sphingomyelin binding region in CYP2B4. The protein-induced lipid raft formation increased the thermal stability of P450 and dramatically altered ligand binding kinetics of the hydrophilic ligand BHT. These results unveil membrane/protein dynamics that contribute to the delicate mechanism of redox catalysis in lipid membrane.

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.

The cytochrome P450 24A1 interaction with adrenodoxin relies on multiple recognition sites that vary among species

Mitochondrial cytochromes P450 (P450s) are responsible for important metabolic reactions, including steps involved in steroid and vitamin D metabolism. The mitochondrial P450 24A1 (CYP24A1) is responsible for deactivation of the bioactive form of vitamin D, 1,25(OH)₂D3. Its function relies on formation of a P450-redox partner complex with the ferredoxin and electron donor adrenodoxin (Adx). However, very little is known about how the Adx-CYP24A1 complex forms. In this study, we report the results of solution NMR in which we monitor isotopically labeled full-length Adx as it binds CYP24A1 in complex with the P450 inhibitor clotrimazole. The NMR titration data suggested a mode for P450-Adx interactions in which formation of the complex relies on contributions from multiple recognition sites on the Adx core domain, some of which have not previously been reported. To evaluate differences among CYP24A1-Adx complexes from different mammalian species and displaying distinct regioselectivity for 1,25(OH)₂D3, all bound spectra were acquired in parallel for human (carbon-23 and -24 hydroxylase), rat (carbon-24 hydroxylase), and opossum (carbon-23 hydroxylase) CYP24A1 isoforms. Binding data from a series of single and double charge-neutralizing substitutions of Adx confirmed that species-specific CYP24A1 isoforms differ in binding to Adx, providing evidence that variations in redox partner interactions correlate with P450 regioselectivity. In summary, these findings reveal that CYP24A1-Adx interactions rely on several recognition sites and that variations in CYP24A1 isoforms modulate formation of the complex, thus providing insight into the variable and complex nature of mitochondrial P450-Adx interactions.

Effect of Redox Partner Binding on Cytochrome P450 Conformational Dynamics

Previous crystal structures of cytochrome P450cam complexed with its redox partner, putidaredoxin (Pdx), shows that P450cam adopts the open conformation. It has been hypothesized that the Pdx-induced shift toward the open state frees the essential Asp251 from salt bridges with Arg186 and Lys178 so that Asp251 can participate in a proton relay network required for O₂ activation. This in part explains why P450cam has such a strict requirement for Pdx. One problem with this view is that looser substrate-protein interactions in the open state may not be compatible with the observed regio- and stereoselective hydroxylation. In the present study, molecular dynamics simulations show that Pdx binding favors a conformation that stabilizes the active site and decreases camphor mobility yet retains a partially open conformation compatible with the required proton relay network. The R186A mutant which frees Asp251 in the absence of Pdx retains good enzyme activity, and the crystal structure shows that product, 5-exo-hydroxycamphor, is bound. This indicates that rupture of the Asp251-Arg186 relaxes selectivity with respect to source of electrons and enables X-ray generated reducing equivalents to support substrate hydroxylation. These combined computational and experimental results are consistent with the proposed role of Pdx in assisting the release of Asp251 from ion pairs so that it can participate in proton-coupled electron transfer.

Structural insights into the mechanism of the drastic changes in enzymatic activity of the cytochrome P450 vitamin D3 hydroxylase (CYP107BR1) caused by a mutation distant from the active site

Cytochromes P450 (P450s) are haem-containing enzymes that catalyze medically and industrially important oxidative reactions, and many P450s have been subjected to directed evolution and site-directed mutagenesis to improve their activity and substrate specificity. Nonetheless, in most cases the mechanism that leads to drastic changes in specific activity after the introduction of an amino-acid substitution distant from the active-site pocket is unclear. Here, two crystal structures of inactive mutants of the P450 vitamin D3 hydroxylase (Vdh), Vdh-F106V and Vdh-L348M, which were obtained in the course of protein-engineering experiments on Vdh, are reported. The overall structures of these mutants show an open conformation similar to that of wild-type Vdh (Vdh-WT), whereas a rearrangement of the common main-chain hydrogen bonds is observed in the CD-loop (residues 102-106), resulting in a more compactly folded CD-loop relative to that of Vdh-WT. The previously reported structures of Vdh-WT and of the highly active Vdh-T107A and Vdh-K1 mutants have a more stretched CD-loop, with partial formation of 310-helix-type hydrogen bonds, both in the open and closed states. Molecular-dynamics simulations also showed that the frequency of the 310-helix is significantly reduced in Vdh-F106V and Vdh-L348M. The closed conformation is crucial for substrate and ferredoxin binding to initiate the catalytic reaction of Vdh. Therefore, it is implied that the small local structural changes observed in this study might disrupt the conformational transition from the open to the closed state, thereby leading to a complete loss of vitamin D3 hydroxylase activity.

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.

A short history of heme dioxygenases: rise, fall and rise again

It is well established that there are two different classes of enzymes—tryptophan 2,3-dioxygenase (TDO) and indoleamine 2,3-dioxygenase (IDO)—that catalyse the O₂-dependent oxidation of l-tryptophan to N-formylkynurenine. But it was not always so. This perspective presents a short history of the early TDO and IDO literature, the people that were involved in creating it, and the legacy that this left for the future.

A Comparative Analysis of the Effector Role of Redox Partner Binding in Bacterial P450s

The camphor monooxygenase, cytochrome P450cam, exhibits a strict requirement for its own redox partner, putidaredoxin (Pdx), a two-iron-sulfur ferredoxin. The closest homologue to P450cam, CYP101D1, is structurally very similar, uses a similar redox partner, and exhibits nearly identical enzymatic properties in the monooxygenation of camphor to give the same single 5-exo-hydroxy camphor product. However, CYP101D1 does not strictly require its own ferredoxin (Arx) for activity because Pdx can support CYP101D1 catalysis but Arx cannot support P450cam catalysis. We have further examined the differences between these two P450s by determining the effect of spin equilibrium, redox properties, and stability of oxygen complexes. We find that Arx shifts the spin state equilibrium toward high-spin, which is the opposite of the effect of Pdx on P450cam. In both P450s, redox partner binding destabilizes the oxy-P450 complex but this effect is much weaker with CYP101D1. In addition, resonance Raman data show that structural perturbations observed in P450cam upon addition of Pdx are absent in CYP101D1. These data indicate that Arx does not play the same effector role in catalysis as Pdx does with P450cam. The most relevant structural difference between these two P450s centers on a catalytically important Asp residue required for proton-coupled electron transfer. We postulate that with P450cam larger Pdx-assisted motions are required to free this Asp for catalysis while the smaller number of restrictions in CYP101D1 precludes the need for redox partner-assisted structural changes.