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.

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 (Fe III ) 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 CYP450 BM3 and CYP450 CAM and the mammalian enzyme CYP450 3A4 . The treatment is limited to the catalytic cycle, as aspects of two-state reactivity were reviewed previously (e.g., 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 Fe III to Fe II occur? Is the ET coordinated with the entrance of O 2 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.

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.

Glycosyltransferases (GTs) catalyze the formation of glycosidic bonds in carbohydrates and glycoconjugates, with various outcomes depending not only on the acceptor molecules they bind but also on the type of glycosidic bond they form (C−O, C−N, C−S, or C−C). Here we show that the glucosyltransferase UGT1 from the indigo plant Polygonum tinctorium catalyzes either N-, O-, or S-glycosylation with similar rates. We solve the structure of the enzyme in complex with its donor and acceptor substrates and elucidate the molecular basis of N-, O-, and S-specificities using experimental mutagenesis and QM/MM simulations, revealing distinct mechanisms for N-, O-, and S-glycosylation. We also show that the active site can be engineered to increase or favor one of the three glycosylation activities over another. These results will foster the design of more active and specific enzyme variants for production of glycosides.

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.