Substrate Conformation Regulates Aromatic C-H Vs C-F Bond Activation in Heme-Dependent Tyrosine Hydroxylase.

A recently discovered heme-dependent enzyme tyrosine hydroxylase (TyrH) offers a green approach for functionalizing the high-strength C-H and C-F bonds in aromatic compounds. However, there is ambiguity regarding the nature of the oxidant (compound 0 or compound I) involved in activating these bonds. Herein, using comprehensive molecular dynamics (MD) simulations and hybrid quantum mechanical/molecular mechanical calculations, we reveal that it is compound I (Cpd I) that acts as the primary oxidant involved in the functionalization of both C-F and C-H bonds. The energy barrier for C-H and C-F activation using compound 0 (Cpd 0) as an oxidant was very high, indicating that Cpd 0 cannot be an oxidant. Consistent with the previous experimental finding, our simulation shows two different conformations of the substrate, where one orientation favors the C-H activation, while the other conformation prefers the C-F activation. As such, our mechanistic study shows that nature utilizes just one oxidant, that is, Cpd I, but it is the active site conformation that decides whether it selects C-F or C-H functionalization which may resemble involvement of two different oxidants.

How Can Static and Oscillating Electric Fields Serve in Decomposing Alzheimer's and Other Senile Plaques?

Alzheimer's disease is one of the most common neurodegenerative conditions, which are ascribed to extracellular accumulation of ß-amyloid peptides into plaques. This phenomenon seems to typify other related neurodegenerative diseases. The present study uses classical molecular-dynamics simulations to decipher the aggregation-disintegration behavior of ß-amyloid peptide plaques in the presence of static and oscillating oriented external electric fields (OEEFs). A long-term disintegration of such plaques is highly desirable since this may improve the prospects of therapeutic treatments of Alzheimer's disease and of other neurodegenerative diseases typified by senile plaques. Our study illustrates the spontaneous aggregation of the ß-amyloid, its prevention and breakdown when OEEF is applied, and the fate of the broken aggregate when the OEEF is removed. Notably, we demonstrate that the usage of an oscillating OEEF on ß-amyloid aggregates appears to lead to an irreversible disintegration. Insight is provided into the root causes of the various modes of aggregation, as well as into the different fates of OEEF-induced disintegration in oscillating vs static fields. Finally, our simulation results are compared to the well-established TTFields and the Deep Brain Stimulation (DBS) therapies, which are currently used options for treatments of Alzheimer's disease and other related neurodegenerative diseases.

Electro-Freezing of Supercooled Water Is Induced by Hydrated Al3+ and Mg2+ Ions: Experimental and Theoretical Studies.

This work reports that the octahedral hydrated Al3+ and Mg2+ ions operate within electrolytic cells as kosmotropic (long-range order-making) "ice makers" of supercooled water (SCW). 10-5 M solutions of hydrated Al3+ and Mg2+ ions each trigger, near the cathode (-20 ± 5 V), electro-freezing of SCW at -4 °C. The hydrated Al3+ ions do so with 100% efficiency, whereas the Mg2+ ions induce icing with 40% efficiency. In contrast, hydrated Na+ ions, under the same experimental conditions, do not induce icing differently than pure water. As such, our study shows that the role played by Al3+ and Mg2+ ions in water electro-freezing is impacted by two synchronous effects: (1) a geometric effect due to the octahedral packing of the coordinated water molecules around the metallic ions, and (2) the degree of polarization which these two ions induce and thereby acidify the coordinated water molecules, which in turn imparts them with an ice-like structure. Long-duration molecular dynamics (MD) simulations of the Al3+ and Mg2+ indeed reveal the formation of "ice-like" hexagons in the vicinity of these ions. Furthermore, the MD shows that these hexagons and the electric fields of the coordinate water molecules give rise to ultimate icing. As such, the MD simulations provide a rational explanation for the order-making properties of these ions during electro-freezing.

Design and Synthesis of Co-initiators via Base-Catalysed Sequential Conjugate Addition: Application in Photoinduced Radical Polymerisation Reaction.

Applications of photochemistry are becoming very popular in modern-day life due to its operational simplicity, environmentally friendly and economically sustainable nature in comparison to thermochemistry. In particular photoinduced radical polymerisation (PRP) reactions are finding more biological applications and especially in the areas of dental restoration processes, tissue engineering and artificial bone generation. A type-II photoinitiator and co-initiator-promoted PRP turned out to be a cost-effective protocol, and herein we report the design and synthesis of a new efficient co-initiator for a PRP reaction via a barrierless sequential conjugate addition reaction. Experimental mechanistic observations have been further complemented by computational data. Time for newly synthesised 1,2-benzenedithiol (DTH) based co-initiator promoted polymerisation of urethane dimethacrylate (UDMA, 70 %) and triethylene glycol dimethacrylate (TEGDMA, 30 %) in presence of 450 nm LED (15 W) under the aerobic conditions is 38 seconds. Polymeric material has high glass transition temperature, improved mechanical strength (860 BHN) and longer in-depth polymerisation (3 cm).

Tyr 118-Mediated Electron Transfer is Key to the Chlorite Decomposition in Heme-Dependent Chlorite Dismutase.

Chlorite dismutase (Cld) is a crucial enzyme that catalyzes the decomposition of chlorite ions into chloride ions (Cl-) and molecular oxygen (O2). Despite playing an important role in the detoxification of toxic chlorite ions, the mechanism of cleavage of the Cl-O bond by Cld remains highly debatable. The present study highlights the mechanism of such Cl-O bond cleavage in Cld using sophisticated computational tools such as hybrid quantum mechanical/molecular mechanical calculations and long-time scale molecular dynamics simulations. Here, we show that Cld forms a high spin ferric hexacoordinated substrate adduct in the presence of a chlorite ion, which subsequently reduces to a ferrous state. Our study shows a stepwise pathway with the homolytic cleavage of the Cl-O bond that produces a high spin Fe(III)-OH species and a diradicaloid species formed by the combination of a chlorine-based ClO• radical and a protein-based tyrosine118• radical. The findings provide significant insights into Cl-O bond cleavage and O2 formation which shows a crucial role of the tyrosine118 during the electron transfer process.

A single site mutation can induce functional promiscuity in homoserine kinase.

L-Homoserine kinase is crucial in the biosynthesis of L-threonine, L-isoleucine, and L-methionine, where it catalyzes ATP-dependent phosphorylation of L-homoserine (Hse) to yield L-homoserine phosphate as its native activity. However, a single site mutation of H138 → L shows the emergence of ATPase activity as a secondary function. However, a previous mechanistic study proposes direct involvement of ATP and the substrate without any catalytic base; therefore, how the mutation of H138 → L causes the secondary function remains an enigma. Using computational tools herein, we provide new insight into the catalytic mechanism of L-homoserine kinase, showing direct involvement of H138 as a catalytic base. We show that mutation of H138 → L opens a new water channel connecting ATP, which facilitates the ATPase activity and reduces the native activity. The proposed mechanism agrees with the experimental finding that an H138 → L mutation reduces the kinase activity but enhances a promiscuous function, i.e. ATPase activity. Since homoserine kinase catalyzes the biosynthesis of amino acids, we believe that an accurate mechanism could be significant for enzyme engineering to synthesize amino acid analogs.

Tuning Potency of Bioactive Molecules via Polymorphic Modifications: A Case Study.

Polymorphism in drugs and bioactive molecules is not uncommon, and it has remained as one of the critical issues in drug development processes. While improving physicochemical properties of bioactive molecules has been a prime focus of the pharmaceutical chemists, not much efforts have been put toward the improvement of their potency via polymorphic modifications. Here, we consider five cases of 5-arylidene-2-aminothiazolidinones derivatives, the known anticancer agents, and discover eights polymorphs in three out of the five cases. We perform systematic crystallization experiments and detailed crystal structure analysis of the eight polymorphs and two compounds, estimate both their energetic and thermal stabilities, and compare their solid state properties. We also compare in-solution properties, e.g., equilibrium solubility, intrinsic dissolution rate, and phase stability, of three polymorphs of one of the cases. Further, we study the extent of inhibition imposed by those eight polymorphs and seven bulk and crystal forms of the compounds on the proliferation of MCF7 breast cancer cells and also the extent of their binding to the isozyme γ-enolase. Furthermore, we perform MD simulations on the eight polymorphs and one compound to estimate and compare their binding affinity with γ-enolase. Our experimental and MD simulation analyses in general emphasize the importance of polymorphism in improving the biological potency of individual molecules.

Local Electric Fields Dictate Function: The Different Product Selectivities Observed for Fatty Acid Oxidation by Two Deceptively Very Similar P450-Peroxygenases OleT and BSß.

Cytochrome P450 peroxygenases use hydrogen peroxide to hydroxylate long-chain fatty acids by bypassing the use of O2 and a redox partner. Among the peroxygenases, P450OleT uniquely performs decarboxylation of fatty acids and production of terminal olefins. This route taken by P450OleT is intriguing, and its importance is augmented by the practical importance of olefin production. As such, this mechanistic choice merits elucidation. To address this puzzle, we use hybrid QM/MM calculations and MD simulations for the OleT enzyme as well as for the structurally analogous enzyme, P450BSß. The study of P450OleT reveals that the protonated His85 in the wild-type P450OleT plays a crucial role in steering decarboxylation activity by stabilizing the corresponding hydroxoiron(IV) intermediate (Cpd II). In contrast, for P450BSß in which Q85 replaces H85, the respective Cpd II species is unstable and it reacts readily with the substrate radical by rebound, producing hydroxylation products. As shown, this single-site difference creates in P450OleT a local electric field (LEF), which is significantly higher than that in P450BSß. In turn, these LEF differences are responsible for the different stabilities of the respective Cpd II/radical intermediates and hence for different functions of the two enzymes. P450BSß uses the common rebound mechanism and leads to hydroxylation, whereas P450OleT proceeds via decarboxylation and generates terminal olefins. Olefin production projects the power of a single residue to alter the LEF and the enzyme's function.

The Performance of Different Water Models on the Structure and Function of Cytochrome P450 Enzymes.

Modeling approaches and modern simulations to investigate the biomolecular structure and function rely on various methods. Since water molecules play a crucial role in all sorts of chemistry, the accurate modeling of water molecules is vital for such simulations. In cytochrome P450 (CYP450), in particular, water molecules play a key role in forming active oxidant that ultimately performs oxidation and metabolism. In the present study, we have highlighted the behavior of the three most widely used water models─TIP3P, SPC/E, and OPC─for three different CYP450 enzymes─CYP450BM3, CYP450OleT, and CYP450BSß─during MD simulations and QM/MM calculations. We studied the various properties, such as RMSD, RMSF, H-bond, water occupancy, and hydrogen atom transfer (HAT), using QM/MM calculations and compared them for all three water models. Our study shows that the stabilities of the enzyme complexes are well maintained in all three water models. However, the OPC water model performs well for the polar active sites, that is, in CYP450OleT and CYP450BSß, while the TIP3P water model is superior for the hydrophobic site, such as CYP450BM3.

Can the local electric field be a descriptor of catalytic activity? A case study on chorismate mutase.

The current theoretical perception of enzymatic activity is highly reliant on the determination of the activation energy of the reactions, which is often calculated using computationally demanding quantum mechanical calculations. With the ever-increasing use of bioengineering techniques that produce too many variants of the same enzyme, a fast and accurate way to study the relative efficiency of enzymes is currently in high demand. Here, we propose the local electric field (LEF) of the enzyme along the reaction axis as a descriptor for the enzymatic activity using the example of chorismate mutase in its native form and several variants (R90A, R90G, and R90K/C88S). The study shows a direct correlation between the calculated enzymatic EF and the enzymatic activity for all the complexes. MD simulations of the Michaelis complex and the transition state analog (TSA) show a stabilizing force on the TSA due to the enzymatic EF. QM/MM and QM-only DFT calculations in the presence of an external electric field (EEF) oriented along the reaction axis show that the electric field can interact with the dipole moment of the TS, thereby stabilizing it and thus lowering the activation energy.

Simulations reveal the key role of Arg15 in the promiscuous activity in the HisA enzyme.

The HisA enzyme catalyzes the first step of histidine biosynthesis via the Amadori rearrangement of the substrate ProFAR. Since it possesses the most conserved and ancient TIM-barrel fold, it provides an ideal framework for bioengineering of a new function from ancestral enzymes. In the present study, first, the catalytic mechanism of HisA biosynthesis was elucidated using hybrid Quantum Mechanical/Molecular Mechanical calculations, and thereafter, key residues contributing towards the promiscuity for TrpF activity were revealed using several MD simulations of a wild type enzyme and its variant with the native (ProFAR) and promiscuous (PRA) substrates. Our study reveals that the two loops (ßα)1 and (ßα)5 on the catalytic site of the HisA enzyme have incredible adaptability for the native and promiscuous substrates. The conformational interplay between these two loops is substrate driven and precise bioengineering targeting these loops is key to the emergence of new functions. Furthermore, the study reveals a key role of the Arg 15 residue which is close to the catalytic center of the enzyme in the bifunctionality of the HisA enzyme by increasing the loop flexibility. Therefore, our study provides crucial information for future bioengineering work to use the HisA enzyme as a scaffold for new enzymatic activity.

Drug Repurposing to Identify Nilotinib as a Potential SARS-CoV-2 Main Protease Inhibitor: Insights from a Computational and In Vitro Study.

COVID-19, an acute viral pneumonia, has emerged as a devastating pandemic. Drug repurposing allows researchers to find different indications of FDA-approved or investigational drugs. In this current study, a sequence of pharmacophore and molecular modeling-based screening against COVID-19 Mpro (PDB: 6LU7) suggested a subset of drugs, from the Drug Bank database, which may have antiviral activity. A total of 44 out of 8823 of the most promising virtual hits from the Drug Bank were subjected to molecular dynamics simulation experiments to explore the strength of their interactions with the SARS-CoV-2 Mpro active site. MD findings point toward three drugs (DB04020, DB12411, and DB11779) with very low relative free energies for SARS-CoV-2 Mpro with interactions at His41 and Met49. MD simulations identified an additional interaction with Glu166, which enhanced the binding affinity significantly. Therefore, Glu166 could be an interesting target for structure-based drug design. Quantitative structural-activity relationship analysis was performed on the 44 most promising hits from molecular docking-based virtual screening. Partial least square regression accurately predicted the values of independent drug candidates' binding energy with impressively high accuracy. Finally, the EC50 and CC50 of 10 drug candidates were measured against SARS-CoV-2 in cell culture. Nilotinib and bemcentinib had EC50 values of 2.6 and 1.1 µM, respectively. In summary, the results of our computer-aided drug design provide a roadmap for rational drug design of Mpro inhibitors and the discovery of certified medications as COVID-19 antiviral therapeutics.

Does Antibody Stabilize the Ligand Binding in GP120 of HIV-1 Envelope Protein? Evidence from MD Simulation.

CD4-mimetic HIV-1 entry inhibitors are small sized molecules which imitate similar conformational flexibility, in gp120, to the CD4 receptor. However, the mechanism of the conformational flexibility instigated by these small sized inhibitors is little known. Likewise, the effect of the antibody on the function of these inhibitors is also less studied. In this study, we present a thorough inspection of the mechanism of the conformational flexibility induced by a CD4-mimetic inhibitor, NBD-557, using Molecular Dynamics Simulations and free energy calculations. Our result shows the functional importance of Asn425 in substrate induced conformational dynamics in gp120. The MD simulations of Asn425Gly mutant provide a less dynamic gp120 in the presence of NBD-557 without incapacitating the binding enthalpy of NBD-557. The MD simulations of complexes with the antibody clearly show the enhanced affinity of NBD-557 due to the presence of the antibody, which is in good agreement with experimental Isothermal Titration Calorimetry results (Biochemistry2006, 45, 10973-10980).

Cytochrome P450-The Wonderful Nanomachine Revealed through Dynamic Simulations of the Catalytic Cycle.

This Account addresses the catalytic cycle of the enzyme cytochrome P450 (CYP450) as a prototypical biological machine with automatic features. CYP450 is a nanomachine that uses dioxygen and two reducing and two proton equivalents to oxidize a plethora of molecules (so-called substrates) as a means of supplying bio-organisms with essential molecules (e.g., brain neurotransmitters, sex hormones, etc.) and protecting biosystems against poisoning. An enticing property of CYP450s is that entrance of an oxidizable substrate into the active site initiates a series of events that constitute the catalytic cycle, which functions "automatically" in a regulated sequence of events culminating in the production of the oxidized substrates (e.g., hydroxylated, epoxidized, etc.), oftentimes with remarkable stereo- and regioselectivities. It is timely to demonstrate how theory uses molecular dynamics (MD) simulations and quantum-mechanical/molecular-mechanical (QM/MM) calculations to complement experiments and elucidate the choreography by which the protein regulates the catalytic cycle. CYP450 is a heme enzyme that contains a ferric ion (FeIII) coordinated by a porphyrin ligand, a water molecule, and a cysteinate ligand that is provided by a strategic residue of the encapsulating protein. While many of the individual steps are sufficiently well-understood, we shall provide here an overview of the factors that cause all of the steps to be sequentially coordinated. To this end, we use examples from three different CYP450 enzymes: the bacterial ones CYP450BM3 and CYP450CAM and the mammalian enzyme CYP4503A4. The treatment is limited to the catalytic cycle, as aspects of two-state reactivity were reviewed previously (e.g., Shaik , S. ; et al. Chem. Rev. 2005 , 105 , 2279 ). What are the principles that govern the seeming automatic feature? For example, how do substrate entrance and binding gate the enzyme? How does the reductase attachment to the enzyme affect the next steps? What triggers the attachment of the reductase? How does the electron transfer (ET) that converts FeIII to FeII occur? Is the ET coordinated with the entrance of O2 into the active site? What is the mechanism of the latter step? Since the entrance of the substrate expels the water molecules from the active site, how do water molecules re-enter to form a proton channel, which is necessary for creating the ultimate oxidant Compound I? How do mutations that disrupt the water channel nevertheless create a competent oxidant? By what means does the enzyme produce regio- and stereoselective oxidation products? What triggers the departure of the oxidized product, and how does the exit occur in a manner that generates the resting state ready for the next cycle? This Account shows that the entrance of the substrate triggers all of the ensuing events.

Iodine-Catalyzed Aerobic Diazenylation-Amination of Indole Derivatives.

A mild strategy for consecutive diazenylation and amination of indole moieties has been demonstrated. The functionalization occurs at C3 and C2 carbon atoms, respectively, at the indole scaffold in the presence of catalytic iodine and air at 40 °C in the 1,1,1,3,3,3-hexafluoroisopropanol (HFIP) solvent. It is noteworthy that the aromatic amines are generated in situ by the reaction of aryl hydrazine with iodine. In general, bright red products are obtained in moderate to good yield. Control reactions are conducted to establish the reaction mechanism.

Choreography of the Reductase and P450BM3 Domains Toward Electron Transfer Is Instigated by the Substrate.

The driving force for the electron transfer (ET) step in the catalytic cycle of cytochrome P450BM3 is investigated using three sets of 1 µs molecular dynamic simulations for the resting state of P450 in complex with the flavin (FMN) in the reductase domain. These sets involve the following: (i) substrate-free (SF), (ii) substrate (N-palmitoyl glycine, i.e., NPG)-bound (SB), and (iii) SB with the semiquinone radical anion (SQ-) of FMN. Starting from the X-ray structure of the SF heme domain, we observe that the α1-helix (of the reductase) and the C-helix (of the heme) undergo reorientation to a parallel orientation, which is the thermochemically stable form. The reorientation of the helices pushes away the FMN to a distance of 18.4 Å from the heme's center. When the substrate binds it causes the I-helix of the heme domain to kink and push the C-helix toward the α1-helix, thereby locking the latter two into a stabilized perpendicular conformation, wherein the FMN-heme distance is 12 Å. The distance drops further in the SQ- form, and upon QM/MM geometry optimization the two moieties approach 8.8 Å, which enhances the ET rate (by 104-106 fold) to the heme's Fe3+ ion. These motions are driven by hydrogen bond strengthening between the C- and the α1-helices. Finally, substrate binding leads to formation of an organized water chain connecting the FMN and heme moieties. The water channel assists the ET and couples it to the proton transfer steps that should activate O2 and create the oxo-iron active species.

Molecular Dynamics and QM/MM Calculations Predict the Substrate-Induced Gating of Cytochrome P450 BM3 and the Regio- and Stereoselectivity of Fatty Acid Hydroxylation.

Theory predicts herein enzymatic activity from scratch. We show that molecular dynamics (MD) simulations and quantum-mechanical/molecular mechanics (QM/MM) calculations of the fatty acid hydroxylase P450 BM3 predict the binding mechanism of the fatty acid substrate and its enantio/regioselective hydroxylation by the active species of the enzyme, Compound I. The MD simulations show that the substrate's entrance involves hydrogen-bonding interactions with Pro25, Glu43, and Leu188, which induce a huge conformational rearrangement that closes the substrate channel by pulling together the A helix and the ß1 sheet to the F/G loop. In turn, at the bottom of the substrate's channel, residue Phe87 controls the regioselectivity by causing the substrate's chain to curl up and juxtapose its CH2 positions ω-1/ω-2/ω-3 to Compound I while preventing access to the endmost position, ω-CH3. Phe87 also controls the stereoselectivity by the enantioselective steric blocking of the pro-S C-H bond, thus preferring R hydroxylation. Indeed, the MD simulations of the mutant Phe87Ala predict predominant ω hydroxylation. These findings, which go well beyond the X-ray structural data, demonstrate the predictive power of theory and its insight, which can potentially be used as a partner of experiment for eventual engineering of P450 BM3 with site-selective C-H functionalization capabilities.

Emergence of Function in P450-Proteins: A Combined Quantum Mechanical/Molecular Mechanical and Molecular Dynamics Study of the Reactive Species in the H2O2-Dependent Cytochrome P450SPα and Its Regio- and Enantioselective Hydroxylation of Fatty Acids.

This work uses combined quantum mechanical/molecular mechanical and molecular dynamics simulations to investigate the mechanism and selectivity of H2O2-dependent hydroxylation of fatty acids by the P450SPα class of enzymes. H2O2 is found to serve as the surrogate oxidant for generating the principal oxidant, Compound I (Cpd I), in a mechanism that involves homolytic O-O bond cleavage followed by H-abstraction from the Fe-OH moiety. Our results rule out a substrate-assisted heterolytic cleavage of H2O2 en route to Cpd I. We show, however, that substrate binding stabilizes the resultant Fe-H2O2 complex, which is crucial for the formation of Cpd I in the homolytic pathway. A network of hydrogen bonds locks the HO· radical, formed by the O-O homolysis, thus directing it to exclusively abstract the hydrogen atom from Fe-OH, thereby forming Cpd I, while preventing the autoxoidative reaction, with the porphyrin ligand, and the substrate oxidation. The so formed Cpd I subsequently hydroxylates fatty acids at their α-position with S-enantioselectivity. These selectivity patterns are controlled by the active site: substrate's binding by Arg241 determines the α-regioselectivity, while the Pro242 residue locks the prochiral α-CH2, thereby leading to hydroxylation of the pro-S C-H bond. Our study of the mutant Pro242Ala sheds light on potential modifications of the enzyme's active site in order to modify reaction selectivity. Comparisons of P450SPα to P450BM3 and to P450BSß reveal that function has evolved in these related metalloenzymes by strategically placing very few residues in the active site.

How do Enzymes Utilize Reactive OH Radicals? Lessons from Nonheme HppE and Fenton Systems.

The iron(IV)-oxo (ferryl) intermediate has been amply established as the principal oxidant in nonheme enzymes and the key player in C-H bond activations and functionalizations. In contrast to this status, our present QM/MM calculations of the mechanism of fosfomycin biosynthesis (a broad range antibiotic) by the nonheme HppE enzyme rule out the iron(IV)-oxo as the reactive species in the hydrogen abstraction (H-abstraction) step of the pro-R hydrogen from the (S)-2-hydroxypropylphosphonic substrate. Moreover, the study reveals that the ferryl species is bypassed in HppE, while the actual oxidant is an HO(•) radical hydrogen-bonded to a ferric-hydroxo complex, resulting via the homolytic dissociation of the hydrogen peroxide complex, Fe(II)-H2O2. The computed energy barrier of this pathway is 12.0 kcal/mol, in fair agreement with the experimental datum of 9.8 kcal/mol. An alternative mechanism involves the iron-complexed hydroxyl radical (Fe(III)-(HO(•))) that is obtained by protonation of the iron(IV)-oxo group via the O-H group of the substrate. The barrier for this pathway, 23.0 kcal/mol, is higher than the one in the first mechanism. In both mechanisms, the HO(•) radical is highly selective; its H-abstraction leading to the final cis-fosfomycin product. It appears that HppE is prone to usage of HO(•) radicals for C-H bond activation, because the ferryl oxidant requires a specific H-abstraction trajectory (∠FeOH ∼ 180°) that cannot be met for intramolecular H-abstraction. Thus, this work broadens the landscape of nonheme iron enzymes and makes a connection to Fenton chemistry, with implications on new potential biocatalysts that may harness hydroxyl radicals for C-H bond functionalizations.

Quantum mechanical/molecular mechanical calculated reactivity networks reveal how cytochrome P450cam and Its T252A mutant select their oxidation pathways.

Quantum mechanical/molecular mechanical calculations address the longstanding-question of a "second oxidant" in P450 enzymes wherein the proton-shuttle, which leads to formation of the "primary-oxidant" Compound I (Cpd I), was severed by mutating the crucial residue (in P450cam: Threonine-252-to-Alanine, hence T252A). Investigating the oxidant candidates Cpd I, ferric hydroperoxide, and ferric hydrogen peroxide (Fe(III)(O2H2)), and their reactions, generates reactivity networks which enable us to rule out a "second oxidant" and at the same time identify an additional coupling pathway that is responsible for the epoxidation of 5-methylenylcamphor by the T252A mutant. In this "second-coupling pathway", the reaction starts with the Fe(III)(O2H2) intermediate, which transforms to Cpd I via a O-O homolysis/H-abstraction mechanism. The persistence of Fe(III)(O2H2) and its oxidative reactivity are shown to be determined by interplay of substrate and protein. The substrate 5-methylenylcamphor prevents H2O2 release, while the protein controls the Fe(III)(O2H2) conversion to Cpd I by nailing-through hydrogen-bonding interactions-the conformation of the HO(•) radical produced during O-O homolysis. This conformation prevents HO(•) attack on the porphyrin's meso position, as in heme oxygenase, and prefers H-abstraction from Fe(IV)OH thereby generating H2O + Cpd I. Cpd I then performs substrate oxidations. Camphor cannot prevent H2O2 release and hence the T252A mutant does not oxidize camphor. This "second pathway" transpires also during H2O2 shunting of the cycle of wild-type P450cam, where the additional hydrogen-bonding with Thr252 prevents H2O2 release, and contributes to a successful Cpd I formation. The present results lead to a revised catalytic cycle of Cytochrome P450cam.