Quantum Mechanical Cluster Models for Calculations on Enzymatic Reaction Mechanisms: Set-Up and Accuracy.

Enzymes turnover substrates into products with amazing efficiency and selectivity and as such have great potential for use in biotechnology and pharmaceutical applications. However, details of their catalytic cycles and the origins surrounding the regio- and chemoselectivity of enzymatic reaction processes remain unknown, which makes the engineering of enzymes and their use in biotechnology challenging. Computational modelling can assist experimental work in the field and establish the factors that influence the reaction rates and the product distributions. A popular approach in modelling is the use of quantum mechanical cluster models of enzymes that take the first- and second coordination sphere of the enzyme active site into consideration. These QM cluster models are widely applied but often the results are dependent on model choice and selection. Herein, we show that QM cluster models can produce highly accurate results that reproduce experimental product distributions and free energies of activation, regarded that large cluster models with >300 atoms are used. In this tutorial review, we give general guidelines on the set-up and applications of the QM cluster method and discuss its accuracy and reproducibility. Finally, several representative QM cluster model examples on metal-containing enzymes are presented, which highlight the strength of the approach.

Mechanism of the Oxidative Ring-Closure Reaction during Gliotoxin Biosynthesis by Cytochrome P450 GliF.

During gliotoxin biosynthesis in fungi, the cytochrome P450 GliF enzyme catalyzes an unusual C-N ring-closure step while also an aromatic ring is hydroxylated in the same reaction cycle, which may have relevance to drug synthesis reactions in biotechnology. However, as the details of the reaction mechanism are still controversial, no applications have been developed yet. To resolve the mechanism of gliotoxin biosynthesis and gain insight into the steps leading to ring-closure, we ran a combination of molecular dynamics and density functional theory calculations on the structure and reactivity of P450 GliF and tested a range of possible reaction mechanisms, pathways and models. The calculations show that, rather than hydrogen atom transfer from the substrate to Compound I, an initial proton transfer transition state is followed by a fast electron transfer en route to the radical intermediate, and hence a non-synchronous hydrogen atom abstraction takes place. The radical intermediate then reacts by OH rebound to the aromatic ring to form a biradical in the substrate that, through ring-closure between the radical centers, gives gliotoxin products. Interestingly, the structure and energetics of the reaction mechanisms appear little affected by the addition of polar groups to the model and hence we predict that the reaction can be catalyzed by other P450 isozymes that also bind the same substrate. Alternative pathways, such as a pathway starting with an electrophilic attack on the arene to form an epoxide, are high in energy and are ruled out.

What is the Origin of the Regioselective C3-Hydroxylation of L-Arg by the Nonheme Iron Enzyme Capreomycin C?

The nonheme iron dioxygenase capreomycin C (CmnC) hydroxylates a free L-arginine amino acid regio- and stereospe-cifically at the C3-position as part of the capreomycin antibiotics biosynthesis. Little is known on its structure, catalytic cycle and substrate specificity and, therefore, a comprehensive computational study was performed. A large QM cluster model of CmnC was created of 297 atoms and the mechanisms for C3-H, C4-H and C5-H hydroxylation and C3-C4 desaturation were investigated. All low-energy pathways correspond to radical reaction mechanisms with an initial hydrogen atom abstraction followed by OH rebound to form alcohol product complexes. The work is compared to alternative L-Arg hydroxylating nonheme iron dioxygenases and the differences in active site polarity are compared. We show that a tight hydrogen bonding network in the substrate binding pocket positions the substrate in an ideal orientation for C3-H activation, whereby the polar groups in the substrate binding pocket induce an electric field effect that guides the selectivity.

Catalytic divergencies in the mechanism of L-arginine hydroxylating nonheme iron enzymes.

Many enzymes in nature utilize a free arginine (L-Arg) amino acid to initiate the biosynthesis of natural products. Examples include nitric oxide synthases, which generate NO from L-Arg for blood pressure control, and various arginine hydroxylases involved in antibiotic biosynthesis. Among the groups of arginine hydroxylases, several enzymes utilize a nonheme iron(II) active site and let L-Arg react with dioxygen and α-ketoglutarate to perform either C3-hydroxylation, C4-hydroxylation, C5-hydroxylation, or C4-C5-desaturation. How these seemingly similar enzymes can react with high specificity and selectivity to form different products remains unknown. Over the past few years, our groups have investigated the mechanisms of L-Arg-activating nonheme iron dioxygenases, including the viomycin biosynthesis enzyme VioC, the naphthyridinomycin biosynthesis enzyme NapI, and the streptothricin biosynthesis enzyme OrfP, using computational approaches and applied molecular dynamics, quantum mechanics on cluster models, and quantum mechanics/molecular mechanics (QM/MM) approaches. These studies not only highlight the differences in substrate and oxidant binding and positioning but also emphasize on electronic and electrostatic differences in the substrate-binding pockets of the enzymes. In particular, due to charge differences in the active site structures, there are changes in the local electric field and electric dipole moment orientations that either strengthen or weaken specific substrate C-H bonds. The local field effects, therefore, influence and guide reaction selectivity and specificity and give the enzymes their unique reactivity patterns. Computational work using either QM/MM or density functional theory (DFT) on cluster models can provide valuable insights into catalytic reaction mechanisms and produce accurate and reliable data that can be used to engineer proteins and synthetic catalysts to perform novel reaction pathways.

Computational Study Into the Oxidative Ring-Closure Mechanism During the Biosynthesis of Deoxypodophyllotoxin.

The nonheme iron dioxygenase deoxypodophyllotoxin synthase performs an oxidative ring-closure reaction as part of natural product synthesis in plants. How the enzyme enables the oxidative ring-closure reaction of (-)-yatein and avoids substrate hydroxylation remains unknown. To gain insight into the reaction mechanism and understand the details of the pathways leading to products and by-products we performed a comprehensive computational study. The work shows that substrate is bound tightly into the substrate binding pocket with the C7'-H bond closest to the iron(IV)-oxo species. The reaction proceeds through a radical mechanism starting with hydrogen atom abstraction from the C7'-H position followed by ring-closure and a final hydrogen transfer to form iron(II)-water and deoxypodophyllotoxin. Alternative mechanisms including substrate hydroxylation and an electron transfer pathway were explored but found to be higher in energy. The mechanism is guided by electrostatic perturbations of charged residues in the second-coordination sphere that prevent alternative pathways.

QM/MM Study Into the Mechanism of Oxidative C=C Double Bond Cleavage by Lignostilbene-α,ß-Dioxygenase.

The enzymatic biosynthesis of fragrance molecules from lignin fragments is an important reaction in biotechnology for the sustainable production of fine chemicals. In this work we investigated the biosynthesis of vanillin from lignostilbene by a nonheme iron dioxygenase using QM/MM and tested several suggested proposals via either an epoxide or dioxetane intermediate. Binding of dioxygen to the active site of the protein results in the formation of an iron(II)-superoxo species with lignostilbene cation radical. The dioxygenase mechanism starts with electrophilic attack of the terminal oxygen atom of the superoxo group on the central C=C bond of lignostilbene, and the second-coordination sphere effects in the substrate binding pocket guide the reaction towards dioxetane formation. The computed mechanism is rationalized with thermochemical cycles and valence bond schemes that explain the electron transfer processes during the reaction mechanism. Particularly, the polarity of the protein and the local electric field and dipole moments enable a facile electron transfer and an exergonic dioxetane formation pathway.

Mechanism of CO2 Reduction to Methanol with H2 on an Iron(II)-scorpionate Catalyst.

CO2 utilization is an important process in the chemical industry with great environmental power. In this work we show how CO2 and H2 can be reacted to form methanol on an iron(II) center and highlight the bottlenecks for the reaction and what structural features of the catalyst are essential for efficient turnover. The calculations predict the reactions to proceed through three successive reaction cycles that start with heterolytic cleavage of H2 followed by sequential hydride and proton transfer processes. The H2 splitting process is an endergonic process and hence high pressures will be needed to overcome this step and trigger the hydrogenation reaction. Moreover, H2 cleavage into a hydride and proton requires a metal to bind hydride and a nearby source to bind the proton, such as an amide or pyrazolyl group, which the scorpionate ligand used here facilitates. As such the computations highlight the non-innocence of the ligand scaffold through proton shuttle from H2 to substrate as an important step in the reaction mechanism.

Insights into Cytochrome P450 Enzyme Catalyzed Defluorination of Aromatic Fluorides.

Density functional calculations establish a novel mechanism of aromatic defluorination by P450 Compound I. This is achieved via either an initial epoxide intermediate or through a 1,2-fluorine shift in an electrophilic intermediate, which highlights that the P450s can defluorinate fluoroarenes. However, in the absence of a proton donor a strong Fe-F bond can be obtained as shown from the calculations.

Reactivity Differences of Trigonal Pyramidal Nonheme Iron(IV)-Oxo and Iron(III)-Oxo Complexes: Experiment and Theory.

High-valent metal-oxo species play critical roles in enzymatic catalysis yet their properties are still poorly understood. In this work we report a combined experimental and computational study into biomimetic iron(IV)-oxo and iron(III)-oxo complexes with tight second-coordination sphere environments that restrict substrate access. The work shows that the second-coordination sphere slows the hydrogen atom abstraction step from toluene dramatically and the kinetics is zeroth order in substrate. However, the iron(II)-hydroxo that is formed has a low reduction potential and hence cannot do OH rebound favorably. The tolyl radical in solution then reacts further with alternative reaction partners. By contrast, the iron(IV)-oxo species reacts predominantly through OH rebound to form alcohol products. Our studies show that the oxidation state of the metal influences reactivities and selectivities with substrate dramatically and that enzymes will likely need an iron(IV) center to catalyze C-H hydroxylation reactions.

Caffeine Biodegradation by Cytochrome P450 1A2. What Determines the Product Distributions?

Caffeine is a natural compound found in plant seeds that after consumption by humans effects the central nervous system as well as the cardiovascular system. In general, the cytochrome P450 enzymes in the liver are involved in the biodegradation of caffeine, which gives paraxanthine, theobromine and theophylline products. There has been debate for many years why multiple products are obtained and how their distributions are determined. To this end we performed a high-level computational study using a combination of molecular dynamics and quantum mechanical approaches. A series of quantum chemical cluster models on the mechanism of caffeine activation by P450 model complexes give hydrogen atom abstraction barriers that predicts the correct ordering and statistical distribution of products. Our studies highlight that second-coordination sphere effects and thermochemical properties of the substrate determine the product distributions.

Melatonin Activation by Cytochrome P450 Isozymes: How Does CYP1A2 Compare to CYP1A1?

Cytochrome P450 enzymes are versatile enzymes found in most biosystems that catalyze mono-oxygenation reactions as a means of biosynthesis and biodegradation steps. In the liver, they metabolize xenobiotics, but there are a range of isozymes with differences in three-dimensional structure and protein chain. Consequently, the various P450 isozymes react with substrates differently and give varying product distributions. To understand how melatonin is activated by the P450s in the liver, we did a thorough molecular dynamics and quantum mechanics study on cytochrome P450 1A2 activation of melatonin forming 6-hydroxymelatonin and N-acetylserotonin products through aromatic hydroxylation and O-demethylation pathways, respectively. We started from crystal structure coordinates and docked substrate into the model, and obtained ten strong binding conformations with the substrate in the active site. Subsequently, for each of the ten substrate orientations, long (up to 1 µs) molecular dynamics simulations were run. We then analyzed the orientations of the substrate with respect to the heme for all snapshots. Interestingly, the shortest distance does not correspond to the group that is expected to be activated. However, the substrate positioning gives insight into the protein residues it interacts with. Thereafter, quantum chemical cluster models were created and the substrate hydroxylation pathways calculated with density functional theory. These relative barrier heights confirm the experimental product distributions and highlight why certain products are obtained. We make a detailed comparison with previous results on CYP1A1 and identify their reactivity differences with melatonin.

Second Coordination Sphere Effects on the Mechanistic Pathways for Dioxygen Activation by a Ferritin: Involvement of a Tyr Radical and the Identification of a Cation Binding Site.

Ferritins are ubiquitous diiron enzymes involved in iron(II) detoxification and oxidative stress responses and can act as metabolic iron stores. The overall reaction mechanisms of ferritin enzymes are still unclear, particularly concerning the role of the conserved, near catalytic center Tyr residue. Thus, we carried out a computational study of a ferritin using a large cluster model of well over 300 atoms including its first- and second-coordination sphere. The calculations reveal important insight into the structure and reactivity of ferritins. Specifically, the active site Tyr residue delivers a proton and electron in the catalytic cycle prior to iron(II) oxidation. In addition, the calculations highlight a likely cation binding site at Asp65 , which through long-range electrostatic interactions, influences the electronic configuration and charge distributions of the metal center. The results are consistent with experimental observations but reveal novel detail of early mechanistic steps that lead to an unusual mixed-valent iron(III)-iron(II) center.

Biodegradation of Herbicides by a Plant Nonheme Iron Dioxygenase: Mechanism and Selectivity of Substrate Analogues.

Aryloxyalkanoate dioxygenases are unique herbicide biodegrading nonheme iron enzymes found in plants and hence, from environmental and agricultural point of view they are important and valuable. However, they often are substrate specific and little is known on the details of the mechanism and the substrate scope. To this end, we created enzyme models and calculate the mechanism for 2,4-dichlorophenoxyacetic acid biodegradation and 2-methyl substituted analogues by density functional theory. The work shows that the substrate binding is tight and positions the aliphatic group close to the metal center to enable a chemoselective reaction mechanism to form the C2 -hydroxy products, whereas the aromatic hydroxylation barriers are well higher in energy. Subsequently, we investigated the metabolism of R- and S-methyl substituted inhibitors and show that these do not react as efficiently as 2,4-dichlorophenoxyacetic acid substrate due to stereochemical clashes in the active site and particularly for the R-isomer give high rebound barriers.

Electrostatic Perturbations in the Substrate-Binding Pocket of Taurine/α-Ketoglutarate Dioxygenase Determine its Selectivity.

Taurine/α-ketoglutarate dioxygenase is an important enzyme that takes part in the cysteine catabolism process in the human body and selectively hydroxylates taurine at the C1 -position. Recent computational studies showed that in the gas-phase the C2 -H bond of taurine is substantially weaker than the C1 -H bond, yet no evidence exists of 2-hydroxytaurine products. To this end, a detailed computational study on the selectivity patterns in TauD was performed. The calculations show that the second-coordination sphere and the protonation states of residues play a major role in guiding the enzyme to the right selectivity. Specifically, a single proton on an active site histidine residue can change the regioselectivity of the reaction through its electrostatic perturbations in the active site and effectively changes the C1 -H and C2 -H bond strengths of taurine. This is further emphasized by many polar and hydrogen bonding interactions of the protein cage in TauD with the substrate and the oxidant that weaken the pro-R C1 -H bond and triggers a chemoselective reaction process. The large cluster models reproduce the experimental free energy of activation excellently.

Structure and Functional Differences of Cysteine and 3-Mercaptopropionate Dioxygenases: A Computational Study.

Thiol dioxygenases are important enzymes for human health; they are involved in the detoxification and catabolism of toxic thiol-containing natural products such as cysteine. As such, these enzymes have relevance to the development of Alzheimer's and Parkinson's diseases in the brain. Recent crystal structure coordinates of cysteine and 3-mercaptopropionate dioxygenase (CDO and MDO) showed major differences in the second-coordination spheres of the two enzymes. To understand the difference in activity between these two analogous enzymes, we created large, active-site cluster models. We show that CDO and MDO have different iron(III)-superoxo-bound structures due to differences in ligand coordination. Furthermore, our studies show that the differences in the second-coordination sphere and particularly the position of a positively charged Arg residue results in changes in substrate positioning, mobility and enzymatic turnover. Furthermore, the substrate scope of MDO is explored with cysteinate and 2-mercaptosuccinic acid and their reactivity is predicted.