Rearrangement-Free Hydroxylation of Methylcubanes by a Cytochrome P450: The Case for Dynamical Coupling of C–H Abstraction and Rebound

Dynamical Origin of Rebound versus Dissociation Selectivity during Fe-Oxo-Mediated C-H Functionalization Reactions

The mechanism of catalytic C-H functionalization of alkanes by Fe-oxo complexes is often suggested to involve a hydrogen atom transfer (HAT) step with the formation of a radical-pair intermediate followed by diverging pathways for radical rebound, dissociation, or desaturation. Recently, we showed that in some Fe-oxo reactions, the radical pair is a nonstatistical-type intermediate and dynamic effects control rebound versus dissociation pathway selectivity. However, the effect of the solvent cage on the stability and lifetime of the radical-pair intermediate has never been analyzed. Moreover, because of the extreme complexity of motion that occurs during dynamics trajectories, the underlying physical origin of pathway selectivity has not yet been determined. For the reaction between [(TQA_Cl)FeIVO]+ and cyclohexane, here, we report explicit solvent trajectories and machine learning analysis on transition-state sampled features (e.g., vibrational, velocity, and geometric) that identified the transferring hydrogen atom kinetic energy as the most important factor controlling rebound versus nonrebound dynamics trajectories, which provides an explanation for our previously proposed dynamic matching effect in fast rebound trajectories that bypass the radical-pair intermediate. Manual control of the reaction trajectories confirmed the importance of this feature and provides a mechanism to enhance or diminish selectivity for the rebound pathway. This led to a general catalyst design principle and proof-of-principle catalyst design that showcases how to control rebound versus dissociation reaction pathway selectivity.

C-H Alkylation of Cubanes via Catalytic Generation of Cubyl Radicals

A catalytic method for the C-H alkylation of cubanes is described. Some hydrogen atom transfer catalysts enable the direct abstraction of a hydrogen atom from the C-H bond of cubanes, followed by conjugate addition of the generated cubyl radicals to electron-deficient alkenes. Synthetic applications of the functionalization method developed are also described.

Heme Spin Distribution in the Substrate-Free and Inhibited Novel CYP116B5hd: A Multifrequency Hyperfine Sublevel Correlation (HYSCORE) Study

The cytochrome P450 family consists of ubiquitous monooxygenases with the potential to perform a wide variety of catalytic applications. Among the members of this family, CYP116B5hd shows a very prominent resistance to peracid damage, a property that makes it a promising tool for fine chemical synthesis using the peroxide shunt. In this meticulous study, we use hyperfine spectroscopy with a multifrequency approach (X- and Q-band) to characterize in detail the electronic structure of the heme iron of CYP116B5hd in the resting state, which provides structural details about its active site. The hyperfine dipole-dipole interaction between the electron and proton nuclear spins allows for the locating of two different protons from the coordinated water and a beta proton from the cysteine axial ligand of heme iron with respect to the magnetic axes centered on the iron. Additionally, since new anti-cancer therapies target the inhibition of P450s, here we use the CYP116B5hd system-imidazole as a model for studying cytochrome P450 inhibition by an azo compound. The effects of the inhibition of protein by imidazole in the active-site geometry and electron spin distribution are presented. The binding of imidazole to CYP116B5hd results in an imidazole-nitrogen axial coordination and a low-spin heme FeIII. HYSCORE experiments were used to detect the hyperfine interactions. The combined interpretation of the gyromagnetic tensor and the hyperfine and quadrupole tensors of magnetic nuclei coupled to the iron electron spin allowed us to obtain a precise picture of the active-site geometry, including the orientation of the semi-occupied orbitals and magnetic axes, which coincide with the porphyrin N-Fe-N axes. The electronic structure of the iron does not seem to be affected by imidazole binding. Two different possible coordination geometries of the axial imidazole were observed. The angles between gx (coinciding with one of the N-Fe-N axes) and the projection of the imidazole plane on the heme were determined to be -60° and -25° for each of the two possibilities via measurement of the hyperfine structure of the axially coordinated 14N.

Computational Characterization of the Reactivity of Compound I in Unspecific Peroxygenases

Unspecific peroxygenases (UPOs) are emerging as promising biocatalysts for selective oxyfunctionalization of unactivated C-H bonds. However, their potential in large-scale synthesis is currently constrained by suboptimal chemical selectivity. Improving the selectivity of UPOs requires a deep understanding of the molecular basis of their catalysis. Recent molecular simulations have sought to unravel UPO's selectivity and inform their design principles. However, most of these studies focused on substrate-binding poses. Few researchers have investigated how the reactivity of CpdI, the principal oxidizing intermediate in the catalytic cycle, influences selectivity in a realistic protein environment. Moreover, the influence of protein electrostatics on the reaction kinetics of CpdI has also been largely overlooked. To bridge this gap, we used multiscale simulations to interpret the regio- and enantioselective hydroxylation of the n-heptane substrate catalyzed by Agrocybe aegerita UPO (AaeUPO). We comprehensively characterized the energetics and kinetics of the hydrogen atom-transfer (HAT) step, initiated by CpdI, and the subsequent oxygen rebound step forming the product. Notably, our approach involved both free energy and potential energy evaluations in a quantum mechanics/molecular mechanics (QM/MM) setting, mitigating the dependence of results on the choice of initial conditions. These calculations illuminate the thermodynamics and kinetics of the HAT and oxygen rebound steps. Our findings highlight that both the conformational selection and the distinct chemical reactivity of different substrate hydrogen atoms together dictate the regio- and enantio-selectivity. Building on our previous study of CpdI's formation in AaeUPO, our results indicate that the HAT step is the rate-limiting step in the overall catalytic cycle. The subsequent oxygen rebound step is swift and retains the selectivity determined by the HAT step. We also pinpointed several polar and charged amino acid residues whose electrostatic potentials considerably influence the reaction barrier of the HAT step. Notably, the Glu196 residue is pivotal for both the CpdI's formation and participation in the HAT step. Our research offers in-depth insights into the catalytic cycle of AaeUPO, which will be instrumental in the rational design of UPOs with enhanced properties.

Characterisation of the heme aqua-ligand coordination environment in an engineered peroxygenase cytochrome P450 variant

The cytochrome P450 enzymes (CYPs) are heme-thiolate monooxygenases that catalyse the insertion of an oxygen atom into the C-H bonds of organic molecules. In most CYPs, the activation of dioxygen by the heme is aided by an acid-alcohol pair of residues located in the I-helix of the enzyme. Mutation of the threonine residue of this acid-alcohol pair of CYP199A4, from the bacterium Rhodospeudomonas palustris HaA2, to a glutamate residue induces peroxygenase activity. In the X-ray crystal structures of this variant an interaction of the glutamate side chain and the distal aqua ligand of the heme was observed and this results in this ligand not being readily displaced in the peroxygenase mutant on the addition of substrate. Here we use a range of bulky hydrophobic and nitrogen donor containing ligands in an attempt to displace the distal aqua ligand of the T252E mutant of CYP199A4. Ligand binding was assessed by UV-visible absorbance spectroscopy, native mass spectrometry, electron paramagnetic resonance and X-ray crystallography. None of the ligands tested, even the nitrogen donor ligands which bind directly to the iron in the wild-type enzyme, resulted in displacement of the aqua ligand. Therefore, modification of the I-helix threonine residue to a glutamate residue results in a significant strengthening of the ferric distal aqua ligand. This ligand was not displaced using any of the ligands during this study and this provides a rationale as to why this mutant can shutdown the monooxygenase pathway of this enzyme and switch to peroxygenase activity.

Molecular Dynamics of Iron Porphyrin-Catalyzed C-H Hydroxylation of Ethylbenzene

Quasi-classical molecular dynamics (MD) simulations were carried out to study the mechanism of iron porphyrin-catalyzed hydroxylation of ethylbenzene. The hydrogen atom abstraction from ethylbenzene by iron-oxo species is the rate-determining step, which generates the radical pair of iron-hydroxo species and the benzylic radical. In the subsequent radical rebound step, the iron-hydroxo species and benzylic radical recombine to form the hydroxylated product, which is barrierless on the doublet energy surface. In the gas-phase quasi-classical MD study on the doublet energy surface, 45% of the reactive trajectories lead directly to the hydroxylated product, and this increases to 56% in implicit solvent model simulations. The percentage of reactive trajectories leading to the separated radical pair is 98-100% on high-spin (quartet/sextet) energy surfaces. The low-spin state reactivity dominates in the hydroxylation of ethylbenzene, which is dynamically both concerted and stepwise, since the time gap between C-H bond cleavage and C-O bond formation ranges from 41 to 619 fs. By contrast, the high-spin state catalysis is an energetically stepwise process, which has a negligible contribution to the formation of hydroxylation products.

Cytochrome P450-catalyzed oxidation of halogen-containing substrates

Cytochrome P450 (CYP) enzymes are heme-thiolate monooxygenases which catalyze the oxidation of aliphatic and aromatic C-H bonds and other reactions. The oxidation of halogens by cytochrome P450 enzymes has also been reported. Here we use CYP199A4, from the bacterium Rhodopseudomonas palustris strain HaA2, with a range of para-substituted benzoic acid ligands, which contain halogens, to assess if this enzyme can oxidize these species or if the presence of these electronegative atoms can alter the outcome of P450-catalyzed reactions. Despite binding to the enzyme, there was no detectable oxidation of any of the 4-halobenzoic acids. CYP199A4 was, however, able to efficiently catalyze the oxidation of both 4-chloromethyl- and 4-bromomethyl-benzoic acid to 4-formylbenzoic acid via hydroxylation of the α‑carbon. The 4-chloromethyl substrate bound in the enzyme active site in a similar manner to 4-ethylbenzoic acid. This places the benzylic α‑carbon hydrogens in an unfavorable position for abstraction indicating a degree of substrate mobility must be possible within the active site. CYP199A4 catalyzed oxidations of 4-(2'-haloethyl)benzoic acids yielding α-hydroxylation and desaturation metabolites. The α-hydroxylation product was the major metabolite. The desaturation pathway is significantly disfavored compared to 4-ethylbenzoic acid. This may be due to the electron-withdrawing halogen atom or a different positioning of the substrate within the active site. The latter was demonstrated by the X-ray crystal structures of CYP199A4 with these substrates. Overall, the presence of a halogen atom positioned close to the heme iron can alter the binding orientation and outcomes of enzyme-catalyzed oxidation.

In situ rapid synthesis of hydrogels based on a redox initiator and persistent free radicals

The development of fast and economical hydrogel manufacturing methods is crucial for expanding the application of hydrogels. However, the commonly used rapid initiation system is not conducive to the performance of hydrogels. Therefore, the research focuses on how to improve the preparation speed of hydrogels and avoid affecting the properties of hydrogels. Herein, a redox initiation system with nanoparticle-stabilized persistent free radicals was introduced to rapidly synthesize high-performance hydrogels at room temperature. A redox initiator composed of vitamin C and ammonium persulfate rapidly provides hydroxyl radicals at room temperature. Simultaneously, three-dimensional nanoparticles can stabilize free radicals and prolong their lifetime, thereby increasing the free radical concentration and accelerating the polymerization rate. And casein enabled the hydrogel to achieve impressive mechanical properties, adhesion, and electrical conductivity. This method greatly facilitates the rapid and economical synthesis of high-performance hydrogels and presents broad application prospects in the field of flexible electronics.

Direct Dynamics Trajectories Demonstrate Dynamic Matching and Nonstatistical Radical Pair Intermediates during Fe-Oxo-Mediated C-H Functionalization Reactions

The generally proposed mechanism for the reaction between non-heme Fe-oxo complexes and alkane C-H bonds involves a hydrogen atom transfer (HAT) reaction step with a radical pair intermediate that then has competitive radical rebound, dissociation, or desaturation pathways. Here, we report density functional theory-based quasiclassical direct dynamics trajectories that examine post-HAT reaction dynamics. Trajectories revealed that the radical pair intermediate can be a nonstatistical type intermediate without complete internal vibrational redistribution and post-HAT selectivity is generally determined by dynamic effects. Fast rebound trajectories occur through dynamic matching between the rotational motion of the newly formed Fe-OH bond and collision with the alkane radical, and all of this occurs through a nonsynchronous dynamically concerted process that circumvents the radical pair intermediate structure. For radical pair dissociation, trajectories proceeded to the radical pair intermediate for a very brief time, followed by complete dissociation. These trajectories provide a new viewpoint and model to understand the inherent reaction pathway selectivity for non-heme Fe-oxo-mediated C-H functionalization reactions.

Cage hydrocarbons as linkers in dimeric drug design: Case studies with trimethoprim and tedizolid

The looming threat of a "post-antibiotic era" has been caused by a rapid rise in antibacterial resistance and subsequent depletion of effective antibiotic agents in the clinic. An efficient strategy to address this shortfall lies in the reengineering of pre-existing and commercially available antibiotic drugs. This is exemplified by dimerization, a design concept in which two pharmacophores are covalently linked to form a new chemical entity. The cage hydrocarbons cubane (1), bicyclo[2.2.2]octane (BCO) (2), adamantane (3), and bicyclo[1.1.1]pentane (BCP) (4) present themselves as an attractive family of linkers in this regard. In this report, all four hydrocarbon cages were employed as linkers in a series of dimers based on the commercially available antibiotics trimethoprim and tedizolid. A detailed synthetic roadmap for the protection and deprotection of each pharmacophore is outlined. Several members of the trimethoprim series showed activity on par with that of their trimethoprim progenitor, although this was not the case for the tedizolid series. The design strategy outlined herein highlights the utility of the group as a platform for the rapid and modular construction of future novel antibiotics.

Investigating the Active Oxidants Involved in Cytochrome P450 Catalyzed Sulfoxidation Reactions

Cytochrome P450 (CYP) heme-thiolate monooxygenases catalyze the hydroxylation of the C-H bonds of organic molecules. This reaction is initiated by a ferryl-oxo heme radical cation (Cpd I). These enzymes can also catalyze sulfoxidation reactions and the ferric-hydroperoxy complex (Cpd 0) and the Fe(III)-H₂ O₂ complex have been proposed as alternative oxidants for this transformation. To investigate this, the oxidation of 4-alkylthiobenzoic acids and 4-methoxybenzoic acid by the CYP199A4 enzyme from Rhodopseudomonas palustris HaA2 was compared using both monooxygenase and peroxygenase pathways. By examining mutants at the mechanistically important, conserved acid alcohol-pair (D251N, T252A and T252E) the relative amounts of the reactive intermediates that would form in these reactions were disturbed. Substrate binding and X-ray crystal structures helped to understand changes in the activity and enabled an attempt to evaluate whether multiple oxidants can participate in these reactions. In peroxygenase reactions the T252E mutant had higher activity towards sulfoxidation than O-demethylation but in the monooxygenase reactions with the WT enzyme the activity of both reactions was similar. The peroxygenase activity of the T252A mutant was greater for sulfoxidation reactions than the WT enzyme, which is the reverse of the activity changes observed for O-demethylation. The monooxygenase activity and coupling efficiency of sulfoxidation and oxidative demethylation were reduced by similar degrees with the T252A mutant. These observations infer that while Cpd I is required for O-dealkylation, another oxidant may contribute to sulfoxidation. Based on the activity of the CYP199A4 mutants it is proposed that this is the Fe(III)-H₂ O₂ complex which would be more abundant in the peroxide-driven reactions.

Nonclassical Phenyl Bioisosteres as Effective Replacements in a Series of Novel Open-Source Antimalarials

The replacement of one chemical motif with another that is broadly similar is a common method in medicinal chemistry to modulate the physical and biological properties of a molecule (i.e., bioisosterism). In recent years, bioisosteres such as cubane and bicyclo[1.1.1]pentane (BCP) have been used as highly effective phenyl mimics. Herein, we show the successful incorporation of a range of phenyl bioisosteres during the open-source optimization of an antimalarial series. Cubane (19) and closo-carborane (23) analogues exhibited improved in vitro potency against Plasmodium falciparum compared to the parent phenyl compound; however, these changes resulted in a reduction in metabolic stability; unusually, enzyme-mediated oxidation was found to take place on the cubane core. A BCP analogue (22) was found to be equipotent to its parent phenyl compound and showed significantly improved metabolic properties. While these results demonstrate the utility of these atypical bioisosteres when used in a medicinal chemistry program, the search to find a suitable bioisostere may well require the preparation of many candidates, in our case, 32 compounds.

Structural insights into the role of the acid-alcohol pair of residues required for dioxygen activation in cytochrome P450 enzymes

The cytochrome P450 heme monooxygenases commonly use an acid-alcohol pair of residues, within the I-helix, to activate iron-bound dioxygen. This work aims to clarify conflicting reports on the importance of the alcohol functionality in this process. Mutants of the P450, CYP199A4 (CYP199A4_D251N and CYP199A4_T252A), were prepared, characterised and their crystal structures were solved. The acid residue of CYP199A4 is not part of a salt bridge network, a key feature of paradigmatic model system P450cam. Instead, there is a direct proton delivery network, via a chain of water molecules, extending to the surface. Nevertheless, CYP199A4_D251N dramatically reduced the activity of the enzyme consistent with a role in proton delivery. CYP199A4_T252A decreased the coupling efficiency of the enzyme with a concomitant increase in the hydrogen peroxide uncoupling pathway. However, the effect of this mutation was much less pronounced than reported with P450cam. Its crystal structures revealed fewer changes at the I-helix, compared to the P450cam system. The structural changes observed within the I-helix of P450cam during oxygen activation do not seem to be required in this P450. These differences are due to the presence of a second threonine residue at position 253, which is absent in P450cam. This threonine forms part of the hydrogen bonding network, resulting in subtle structural changes and is also present across the majority of the P450 superfamily. Overall, the results suggest that while the acid-alcohol pair is important for dioxygen activation this process and the method of proton delivery can differ across P450s.Graphic abstract.

Biophysical Techniques for Distinguishing Ligand Binding Modes in Cytochrome P450 Monooxygenases

The cytochrome P450 superfamily of heme monooxygenases catalyzes important chemical reactions across nature. The changes in the optical spectra of these enzymes, induced by the addition of substrates or inhibitors, are critical for assessing how these molecules bind to the P450, enhancing or inhibiting the catalytic cycle. Here we use the bacterial CYP199A4 enzyme (Uniprot entry Q2IUO2), from Rhodopseudomonas palustris HaA2, and a range of substituted benzoic acids to investigate different binding modes. 4-Methoxybenzoic acid elicits an archetypal type I spectral response due to a ≥95% switch from the low- to high-spin state with concomitant dissociation of the sixth aqua ligand. 4-(Pyridin-3-yl)- and 4-(pyridin-2-yl)benzoic acid induced different type II ultraviolet-visible (UV-vis) spectral responses in CYP199A4. The former induced a greater red shift in the Soret wavelength (424 nm vs 422 nm) along with a larger overall absorbance change and other differences in the α-, β-, and δ-bands. There were also variations in the ferrous UV-vis spectra of these two substrate-bound forms with a spectrum indicative of Fe-N bond formation with 4-(pyridin-3-yl)benzoic acid. The crystal structures of CYP199A4, with the pyridinyl compounds bound, revealed that while the nitrogen of 4-(pyridin-3-yl)benzoic acid is coordinated to the heme, with 4-(pyridin-2-yl)benzoic acid an aqua ligand remains. Continuous wave and pulse electron paramagnetic resonance data in frozen solution revealed that the substrates are bound in the active site in a form consistent with the crystal structures. The redox potential of each CYP199A4-substrate combination was measured, allowing correlation among binding modes, spectroscopic properties, and the observed biochemical activity.

Complementary and selective oxidation of hydrocarbon derivatives by two cytochrome P450 enzymes of the same family

The cytochrome P450 enzymes CYP101B1 and CYP101C1, from a Novosphingobium bacterium, can efficiently hydroxylate hydrocarbon derivatives containing a carbonyl moiety. Cyclic ketones (C9 to C15) were oxidised with contrasting yet high selectivity.