Protein environment facilitates O2 binding in non-heme iron enzyme. An insight from ONIOM calculations on isopenicillin N synthase (IPNS)

Binding of dioxygen to a non-heme enzyme has been modeled using the ONIOM combined quantum mechanical/molecular mechanical (QM/MM) method. For the present system, isopenicillin N synthase (IPNS), binding of dioxygen is stabilized by 8-10 kcal/mol for a QM:MM (B3LYP:Amber) protein model compared to a quantum mechanical model of the active site only. In the protein system, the free energy change of O2 binding is close to zero. Two major factors consistently stabilize O2 binding. The first effect, evaluated at the QM level, originates from a change in coordination geometry of the iron center. The active-site model artificially favors the deoxy state (O2 not bound) because it allows too-large rearrangements of the five-coordinate iron site. This error is corrected when the protein is included. The corresponding effect on binding energies is 3-6 kcal/mol, depending on the coordination mode of O2 (side-on or end-on). The second major factor that stabilizes O2 binding is van der Waals interactions between dioxygen and the surrounding enzyme. These interactions, 3-4 kcal/mol at the MM level, are neglected in models that include only the active site. Polarization of the active site by surrounding amino acids does not have a significant effect on the binding energy in the present system.     

Toward understanding macrocycle specificity of iron on the dioxygen-binding ability: a theoretical study

In an effort to examine the interaction between dioxygen and iron-macrocyclic complexes, and to understand how this interaction was affected by those different macrocyclic ligands, dioxygen binding with iron-porphyrin, iron-phthalocyanine, iron-dibenzotetraaza[14]annulene, and iron-salen complexes is investigated by means of quantum chemical calculations utilizing Density Functional Theory (DFT). Based on the analysis of factors influencing the corresponding dioxygen binding process, it showed that different macrocyclic ligands possess different O-O bond distances, and different electronic configurations for the bound O(2) and non-aromatic macrocyclic ligands favor dioxygen activation. Furthermore, the smaller the energy gap between the HOMO of iron-macrocyclic complexes and the LUMO of dioxygen, the more active the bound O(2) becomes, with a longer O-O bond distance and a shorter Fe-O bond length.     

Studies of iron(II) and iron(III) complexes with fac-N2O, cis-N2O2 and N2O3 donor ligands: models for the 2-His 1-carboxylate motif of non-heme iron monooxygenases

Enzymes in the oxygen-activating class of mononuclear non-heme iron oxygenases (MNOs) contain a highly conserved iron center facially ligated by two histidine nitrogen atoms and one carboxylate oxygen atom that leave one face of the metal center (three binding sites) open for coordination to cofactor, substrate, and/or dioxygen. A comparative family of [Fe(II/III)(N(2)O(n))(L)(4-n))](±x), n = 1-3, L = solvent or Cl(-), model complexes, based on a ligand series that supports a facially ligated N,N,O core that is then modified to contain either one or two additional carboxylate chelate arms, has been structurally and spectroscopically characterized. EPR studies demonstrate that the high-spin d(5) Fe(III)g = 4.3 signal becomes more symmetrical as the number of carboxylate ligands decreases across the series Fe(N(2)O(3)), Fe(N(2)O(2)), and Fe(N(2)O(1)), reflecting an increase in the E/D strain of these complexes as the number of exchangeable/solvent coordination sites increases, paralleling the enhanced distribution of electronic structures that contribute to the spectral line shape. The observed systematic variations in the Fe(II)-Fe(III) oxidation-reduction potentials illustrate the fundamental influence of differential carboxylate ligation. The trend towards lower reduction potential for the iron center across the [Fe(III)(N(2)O(1))Cl(3)](-), [Fe(III)(N(2)O(2))Cl(2)](-) and [Fe(III)(N(2)O(3))Cl](-) series is consistent with replacement of the chloride anions with the more strongly donating anionic O-donor carboxylate ligands that are expected to stabilize the oxidized ferric state. This electrochemical trend parallels the observed dioxygen sensitivity of the three ferrous complexes (Fe(II)(N(2)O(1)) < Fe(II)(N(2)O(2)) < Fe(II)(N(2)O(3))), which form μ-oxo bridged ferric species upon exposure to air or oxygen atom donor (OAD) molecules. The observed oxygen sensitivity is particularly interesting and discussed in the context of α-ketoglutarate-dependent MNO enzyme mechanisms.     

New force field for calcium binding sites in annexin-membrane complexes

For accurate classical molecular dynamics (MD) simulations of the calcium mediated bound complexes of annexin and membrane we have developed new force-field parameters correctly describing the interaction of the Ca ion with its environment. We have used quantum chemical calculations to investigate the potential energy surface experienced by the Ca ion within the three different binding sites found in domain 1 of annexin V (ANX V/1). Based on these calculations we were able to quantify the charge polarization of atoms within the binding sites, and to determine the geometry and force constants of harmonic restraints between the Ca ion and its coordinating oxygen atoms. Harmonic restraints were introduced to compensate for the deviations between the quantum mechanical potential energy surface and that of the classical force field. Our analysis has shown that using the refined force field for the Ca binding sites enables long-time MD simulations that conserve the initial structure of ANX V/1 significantly better than MD simulations using the standard force field.     

Computational analysis of the proton translocation from Asp96 to schiff base in bacteriorhodopsin

The potential energy change during the M --> N process in bacteriorhodopsin has been evaluated by ab initio quantum chemical and advanced quantum chemical calculations following molecular dynamics (MD) simulations. Many previous experimental studies have suggested that the proton transfer from Asp96 to the Schiff base occurs under the following two conditions: (1) the hydrogen bond between Thr46 and Asp96 breaks and Thr46 is detached from Asp96 and (2) a stable chain of four water molecules spans an area from Asp96 --> Schiff base. In this work, we successfully reproduced the proton-transfer process occurring under these two conditions by molecular dynamics and quantum chemical calculations. The quantum chemical computation revealed that the proton transfer from Asp96 to Shiff base occurs in two-step reactions via an intermediate in which an H(3)O(+) appears around Ala215. The activation energy for the proton transfer in the first reaction was calculated to be 9.7 kcal/mol, which enables fast and efficient proton pump action. Further QM/MM (quantum mechanical/molecular mechanical) and FMO (fragment molecular orbital) calculations revealed that the potential energy change during the proton transfer is tightly regulated by the composition and the geometry of the surrounding amino acid residues of bacteriorhodopsin. Here, we report in detail the Asp96 --> Schiff base proton translocation mechanism of bacteriorhodopsin. Additionally, we discuss the effectiveness of combining quantum chemical calculations with truncated cluster models followed by advanced quantum chemical calculations applied to a whole protein to elucidate its reaction mechanism.     

Quantum mechanical map for protein-ligand binding with application to beta-trypsin/benzamidine complex

We report full ab initio Hartree-Fock calculation to compute quantum mechanical interaction energies for beta-trypsin/benzamidine binding complex. In this study, the full quantum mechanical ab initio energy calculation for the entire protein complex with 3238 atoms is made possible by using a recently developed MFCC (molecular fractionation with conjugate caps) approach in which the protein molecule is decomposed into amino acid-based fragments that are properly capped. The present MFCC ab initio calculation enables us to obtain an "interaction spectrum" that provides detailed quantitative information on protein-ligand binding at the amino acid levels. These detailed information on individual residue-ligand interaction gives a quantitative molecular insight into our understanding of protein-ligand binding and provides a guidance to rational design of potential inhibitors of protein targets.     

Experimental and theoretical investigation of topological and energetic characteristics of Sb complexes reversibly binding molecular oxygen

The experimental distribution of electron density in Ph(3)(4,5-OMe-3,6-Bu(t)-Cat)Sb·MeCN (1*) and Ph(3)(4,5-N(2)C(4)H(6)-3,6-Bu(t)-Cat)Sb·MeOH (2*) complexes was studied. According to atoms in molecules theory, the Sb-C(Ph), Sb-O(catecholate), and Sb···N(O) bonds are intermediate, whereas the O-C and C-C bonds are covalent, respectively. The energy of the Sb···N(MeCN) and Sb···O(MeOH) bonds are 7.0 and 11.3 kcal/mol according to the Espinosa equation. Density functional theory and Hartree-Fock calculations were carried out for a series of catecholate and amidophenolate complexes of antimony(V). It was shown that such calculations reliably reproduce geometrical and topological parameters and therefore can be used for a criterion search of dioxygen reversible binding by the catecholate and amidophenolate complexes of antimony(V). It was found that the "critical" value of the HOMO energy vary in the range from -5.197 to -5.061 eV for reversible binding of dioxygen complexes. This can serve as a thermodynamic criterion to predict the possibility of the dioxygen reversible binding by the catecholate and amidophenolate complexes of Sb(V). The HOMO energies correlate with the conversion of the catecholate and amidophenolate complexes in corresponding spiroendoperoxide derivatives as well. The contribution of the atom orbitals of the carbon atoms in the five-membered metallocycle to HOMO in complexes with different substitutes in the 4- and 5-positions of the catecholate ligand allows predicting the place of dioxygen addition.     

Oxygen activation by nonheme iron(II) complexes: alpha-keto carboxylate versus carboxylate

Mononuclear iron(II) alpha-keto carboxylate and carboxylate compounds of the sterically hindered tridentate face-capping ligand Tp(Ph2) (Tp(Ph2) = hydrotris(3,5-diphenylpyrazol-1-yl)borate) were prepared as models for the active sites of nonheme iron oxygenases. The structures of an aliphatic alpha-keto carboxylate complex, [Fe(II)(Tp(Ph2))(O(2)CC(O)CH(3))], and the carboxylate complexes [Fe(II)(Tp(Ph2))(OBz)] and [Fe(II)(Tp(Ph2))(OAc)(3,5-Ph(2)pzH)] were determined by single-crystal X-ray diffraction, all of which have five-coordinate iron centers. Both the alpha-keto carboxylate and the carboxylate compounds react with dioxygen resulting in the hydroxylation of a single ortho phenyl position of the Tp(Ph2) ligand. The oxygenation products were characterized spectroscopically, and the structure of the octahedral iron(III) phenolate product [Fe(III)(Tp(Ph2))(OAc)(3,5-Ph(2)pzH)] was established by X-ray diffraction. The reaction of the alpha-keto carboxylate model compounds with oxygen to produce the phenolate product occurs with concomitant oxidative decarboxylation of the alpha-keto acid. Isotope labeling studies show that (18)O(2) ends up in the Tp(Ph2) phenolate oxygen and the carboxylate derived from the alpha-keto acid. The isotope incorporation mirrors the dioxygenase nature of the enzymatic systems. Parallel studies on the carboxylate complexes demonstrate that the oxygen in the hydroxylated ligand is also derived from molecular oxygen. The oxygenation of the benzoylformate complex is demonstrated to be first order in metal complex and dioxygen, with activation parameters DeltaH++ = 25 +/- 2 kJ mol(-1) and DeltaS++ = -179 +/- 6 J mol(-1) K(-1). The rate of appearance of the iron(III) phenolate product is sensitive to the nature of the substituent on the benzoylformate ligand, exhibiting a Hammett rho value of +1.3 indicative of a nucleophilic mechanism. The proposed reaction mechanism involves dioxygen binding to produce an iron(III) superoxide species, nucleophilic attack of the superoxide at the alpha-keto functionality, and oxidative decarboxylation of the adduct to afford the oxidizing species that attacks the Tp(Ph2) phenyl ring. Interestingly, the alpha-keto carboxylate complexes react 2 orders of magnitude faster than the carboxylate complexes, thus emphasizing the key role that the alpha-keto functionality plays in oxygen activation by alpha-keto acid-dependent iron enzymes.     

Complexes of thiomandelate and captopril mercaptocarboxylate inhibitors to metallo-beta-lactamase by polarizable molecular mechanics. Validation on model binding sites by quantum chemistry

Using the polarizable molecular mechanics method SIBFA, we have performed a search for the most stable binding modes of D- and L-thiomandelate to a 104-residue model of the metallo-beta-lactamase from B. fragilis, an enzyme involved in the acquired resistance of bacteria to antibiotics. Energy balances taking into account solvation effects computed with a continuum reaction field procedure indicated the D-isomer to be more stably bound than the L-one, conform to the experimental result. The most stably bound complex has the S(-) ligand bridging monodentately the two Zn(II) cations and one carboxylate O(-) H-bonded to the Asn193 side chain. We have validated the SIBFA energy results by performing additional SIBFA as well as quantum chemical (QC) calculations on small (88 atoms) model complexes extracted from the 104-residue complexes, which include the residues involved in inhibitor binding. Computations were done in parallel using uncorrelated (HF) as well as correlated (DFT, LMP2, MP2) computations, and the comparisons extended to corresponding captopril complexes (Antony et al., J Comput Chem 2002, 23, 1281). The magnitudes of the SIBFA intermolecular interaction energies were found to correctly reproduce their QC counterparts and their trends for a total of twenty complexes.     

X-ray crystal structures of manganese(II)-reconstituted and native toluene/o-xylene monooxygenase hydroxylase reveal rotamer shifts in conserved residues and an enhanced view of the protein interior

We report the X-ray crystal structures of native and manganese(II)-reconstituted toluene/o-xylene monooxygenase hydroxylase (ToMOH) from Pseudomonas stutzeri OX1 to 1.85 and 2.20 A resolution, respectively. The structures reveal that reduction of the dimetallic active site is accompanied by a carboxylate shift and alteration of the coordination environment for dioxygen binding and activation. A rotamer shift in a strategically placed asparagine 202 accompanies dimetallic center reduction and is proposed to influence protein component interactions. This rotamer shift is conserved between ToMOH and the corresponding residue in methane monooxygenase hydroxylase (MMOH). Previously unidentified hydrophobic pockets similar to those present in MMOH are assigned.     

Carboxylate–phenolate tautomerism in 5‐[(nitrophenyl)diazenyl]salicylate anions

Aryldiazenyl derivatives of salicylic acid and their salts are used as dyes. In these structures, the carboxylate groups are engaged in short contacts with the cations and in hydrogen bonds with water molecules, if present. If both O atoms of the carboxylate group take part in such interactions, the negative charge is delocalized over the two atoms. In the absence of hydrogen bonds and contacts with cations, the negative charge is localized on one of the O atoms. In the crystal structures of tetramethylammonium 2‐hydroxy‐5‐[(E)‐(4‐nitrophenyl)diazenyl]benzoate and tetramethylammonium 2‐hydroxy‐5‐[(E)‐(2‐nitrophenyl)diazenyl]benzoate, both C₄H₁₂N⁺·C₁₃H₈N₃O₅⁻, all the interactions between the cations and anions are weak, and their effect on the geometry of the anions is negligible. Under these conditions, the 2‐nitro‐substituted anion is an almost pure phenol–carboxylate tautomer, whereas in the 4‐nitro‐substituted anion, the phenolic H atom is shifted towards the carboxylate group, and thus the structure of this anion is intermediate between the phenol–carboxylate and phenolate–carboxylic acid tautomeric forms. The probable formation of such an intermediate form is supported by quantum chemical calculations. Being the characteristic feature of this form, a short distance between the phenolic and carboxylate O atoms is observed in the 4‐nitro‐substituted anion, as well as in the structures of some 3,5‐dinitrosalicylates reported in the literature.     

Quantum mechanical/molecular mechanical investigation of the mechanism of C-H hydroxylation of camphor by cytochrome P450cam: theory supports a two-state rebound mechanism

The stereospecific cytochrome P450-catalyzed hydroxylation of the C(5)-H((5-exo)) bond in camphor has been studied theoretically by a combined quantum mechanical/molecular mechanical (QM/MM) approach. Density functional theory is employed to treat the electronic structure of the active site (40-100 atoms), while the protein and solvent environment (ca. 24,000 atoms) is described by the CHARMM force field. The calculated energy profile of the hydrogen-abstraction oxygen-rebound mechanism indicates that the reaction takes place in two spin states (doublet and quartet), as has been suggested earlier on the basis of calculations on simpler models ("two-state reactivity"). While the reaction on the doublet potential energy surface is nonsynchronous, yet effectively concerted, the quartet pathway is truly stepwise, including formation of a distinct intermediate substrate radical and a hydroxo-iron complex. Comparative calculations in the gas phase demonstrate the effect of the protein environment on the geometry and relative stability of intermediates (in terms of spin states and redox electromers) through steric constraints and electronic polarization.     

Unraveling the mechanisms of O2 activation by size-selected gold clusters: transition from superoxo to peroxo chemisorption

The activation of dioxygen is a key step in CO oxidation catalyzed by gold nanoparticles. It is known that small gold cluster anions with even-numbered atoms can molecularly chemisorb O(2) via one-electron transfer from Au(n)(-) to O(2), whereas clusters with odd-numbered atoms are inert toward O(2). Here we report spectroscopic evidence of two modes of O(2) activation by the small even-sized Au(n)(-) clusters: superoxo and peroxo chemisorption. Photoelectron spectroscopy of O(2)Au(8)(-) revealed two distinct isomers, which can be converted from one to the other depending on the reaction time. Ab initio calculations show that there are two close-lying molecular O(2)-chemisorbed isomers for O(2)Au(8)(-): the lower energy isomer involves a peroxo-type binding of O(2) onto Au(8)(-), while the superoxo chemisorption is a slightly higher energy isomer. The computed detachment transitions of the superoxo and peroxo species are in good agreement with the experimental observation. The current work shows that there is a superoxo to peroxo chemisorption transition of O(2) on gold clusters at Au(8)(-): O(2)Au(n)(-) (n = 2, 4, 6) involves superoxo binding and n = 10, 12, 14, 18 involves peroxo binding, whereas the superoxo binding re-emerges at n = 20 due to the high symmetry tetrahedral structure of Au(20), which has a very low electron affinity. Hence, the two-dimensional (2D) Au(8)(-) is the smallest anionic gold nanoparticle that prefers peroxo binding with O(2). At Au(12)(-), although both 2D and 3D isomers coexist in the cluster beam, the 3D isomer prefers the peroxo binding with O(2).     

茂铁基稀土铈配合物的合成与结构

本文利用单羧基二茂铁Fc-COOH和1，1′-二羧基二茂铁HOOC-Fc-COOH作为配体分别合成了双核铈配合物Ce₂(FcCOO)₆ (1)和二维层状配位聚合物Ce₂(OOC-Fc-COO)₃(2)，测定了两种配合物的晶体结构。配合物1中的金属铈离子为九配位结构，分别与周围的羧基二茂铁上的氧原子和作为辅助配体的水分子配位，茂铁间的π-π相互作用将配合物1的二聚体单元连结在一起形成二维的网状结构。配合物2中的金属铈离子亦为九配位结构，分别与周围的羧基二茂铁上的氧原子，作为辅助配体的水分子和甲醇配位形成类似于配合物1的二聚体单元，1，1′-二羧基二茂铁HOOC-Fc-COOH作为桥基配体将二聚体单元连结在一起，形成二维网状的配位聚合物。     

Spectroscopic and computational studies of α-keto acid binding to Dke1: understanding the role of the facial triad and the reactivity of β-diketones

The O(2) activating mononuclear nonheme iron enzymes generally have a common facial triad (two histidine and one carboxylate (Asp or Glu) residue) ligating Fe(II) at the active site. Exceptions to this motif have recently been identified in nonheme enzymes, including a 3His triad in the diketone cleaving dioxygenase Dke1. This enzyme is used to explore the role of the facial triad in directing reactivity. A combination of spectroscopic studies (UV-vis absorption, MCD, and resonance Raman) and DFT calculations is used to define the nature of the binding of the α-keto acid, 4-hydroxyphenlpyruvate (HPP), to the active site in Dke1 and the origin of the atypical cleavage (C2-C3 instead of C1-C2) pattern exhibited by this enzyme in the reaction of α-keto acids with dioxygen. The reduced charge of the 3His triad induces α-keto acid binding as the enolate dianion, rather than the keto monoanion, found for α-keto acid binding to the 2His/1 carboxylate facial triad enzymes. The mechanistic insight from the reactivity of Dke1 with the α-keto acid substrate is then extended to understand the reaction mechanism of this enzyme with its native substrate, acac. This study defines a key role for the 2His/1 carboxylate facial triad in α-keto acid-dependent mononuclear nonheme iron enzymes in stabilizing the bound α-keto acid as a monoanion for its decarboxylation to provide the two additional electrons required for O(2) activation.     

Binding modes of oxalate in UO2(oxalate) in aqueous solution studied with first-principles molecular dynamics simulations. Implications for the chelate effect

Car-Parrinello molecular dynamics simulations are reported for aqueous UO(2)(H(2)O)(n)(C(2)O(4)) (n = 3, 4), calling special attention to the binding modes of oxalate and the thermodynamics of the so-called chelate effect. Based on free energies from thermodynamic integration (BLYP functional), the κ(1),κ(1')-binding mode of the oxalate (with one O atom from each carboxylate coordinating) is more stable than κ(2) (2 O atoms from the same carboxylate) and κ(1) forms by 23 and 39 kJ mol(-1), respectively. The free energy of binding a fourth water ligand to UO(2)(H(2)O)(3)(κ(1)-C(2)O(4)) is computed to be low, 12 kJ mol(-1). Changes of the hydration shell about oxalate during chelate opening are discussed. Composite enthalpies and free energies, obtained from both experiment and quantum-chemical modeling, are proposed for the formation of monodentate UO(2)(H(2)O)(4)(κ(1)-C(2)O(4)). These data suggest that the largest entropy change in the overall complex formation occurs at this stage, and that the subsequent chelate closure under water release is essentially enthalpy-driven.     

Single d-metal atoms on F(s) and F(s+) defects of MgO(001): a theoretical study across the periodic table

Single d-metal atoms on oxygen defects F(s) and F(s+) of the MgO(001) surface were studied theoretically. We employed an accurate density functional method combined with cluster models, embedded in an elastic polarizable environment, and we applied two gradient-corrected exchange-correlation functionals. In this way, we quantified how 17 metal atoms from groups 6-11 of the periodic table (Cu, Ag, Au; Ni, Pd, Pt; Co, Rh, Ir; Fe, Ru, Os; Mn, Re; and Cr, Mo, W) interact with terrace sites of MgO. We found bonding with F(s) and F(s+) defects to be in general stronger than that with O2- sites, except for Mn-, Re-, and Fe/F(s) complexes. In M/F(s) systems, electron density is accumulated on the metal center in a notable fashion. The binding energy on both kinds of O defects increases from 3d- to 4d- to 5d-atoms of a given group, at variance with the binding energy trend established earlier for the M/O2- complexes, 4d < 3d < 5d. Regarding the evolution of the binding energy along a period, group 7 atoms are slightly destabilized compared to their group 6 congeners in both the F(s) and F(s+) complexes; for later transition elements, the binding energy increases gradually up to group 10 and finally decreases again in group 11, most strongly on the F(s) site. This trend is governed by the negative charge on the adsorbed atoms. We discuss implications for an experimental detection of metal atoms on oxide supports based on computed core-level energies.     

Heterovalent trinuclear Co(III)-Co(II)-Co(III) complexes with N-(2-Carboxyethyl)salicylaldimines

The heterovalent trinuclear cobalt complexes [Co₂^IIIL₄ ⁱ · Co^II(H₂O)₄] · nXmY (L ⁱ are deprotonated Schiff bases derived from substituted salicylaldehydes and β-alanine; i = 1–3) were obtained and characterized. An X-ray diffraction study of the trinuclear cobalt complex with N-(2-carboxyethyl)salicylaldimine showed that the central Co(II) ion and the terminal Co(III) ions are linked by bridging carboxylate groups. Either terminal Co(III) atom is coordinated to two ligand molecules. They form an octahedral environment consisting of two azomethine N atoms, two phenolate O atoms, and two O atoms of two carboxylate groups. The central Co(II) atom is coordinated to four water molecules and to two O atoms of two bridging carboxylate ligands involved in the coordination sphere of the terminal Co(III) atoms.     

Parallel iterative reaction path optimization in ab initio quantum mechanical/molecular mechanical modeling of enzyme reactions

The determination of reaction paths for enzyme systems remains a great challenge for current computational methods. In this paper we present an efficient method for the determination of minimum energy reaction paths with the ab initio quantum mechanical/molecular mechanical approach. Our method is based on an adaptation of the path optimization procedure by Ayala and Schlegel for small molecules in gas phase, the iterative quantum mechanical/molecular mechanical (QM/MM) optimization method developed earlier in our laboratory and the introduction of a new metric defining the distance between different structures in the configuration space. In this method we represent the reaction path by a discrete set of structures. For each structure we partition the atoms into a core set that usually includes the QM subsystem and an environment set that usually includes the MM subsystem. These two sets are optimized iteratively: the core set is optimized to approximate the reaction path while the environment set is optimized to the corresponding energy minimum. In the optimization of the core set of atoms for the reaction path, we introduce a new metric to define the distances between the points on the reaction path, which excludes the soft degrees of freedom from the environment set and includes extra weights on coordinates describing chemical changes. Because the reaction path is represented by discrete structures and the optimization for each can be performed individually with very limited coupling, our method can be executed in a natural and efficient parallelization, with each processor handling one of the structures. We demonstrate the applicability and efficiency of our method by testing it on two systems previously studied by our group, triosephosphate isomerase and 4-oxalocrotonate tautomerase. In both cases the minimum energy paths for both enzymes agree with the previously reported paths.     

Mechanism for catechol ring cleavage by non-heme iron intradiol dioxygenases: a hybrid DFT study

The mechanism of the catalytic reaction of protocatechuate 3,4-dioxygenase (3,4-PCD), a representative intradiol dioxygenase, was studied with the hybrid density functional method B3LYP. First, a smaller model involving only the iron first-shell ligands (His460, His462, and Tyr408) and the substrates (catechol and dioxygen) was used to probe various a priori plausible reaction mechanisms. Then, an extended model involving also the most important second-shell groups (Arg457, Gln477, and Tyr479) was used for the refinement of the preselected mechanisms. The computational results suggest that the chemical reactions constituting the catalytic cycle of intradiol dioxygenases involve: (1) binding of the substrate as a dianion, in agreement with experimental suggestions, (2) binding of dioxygen to the metal aided by an electron transfer from the substrate to O(2), (3) formation of a bridging peroxo intermediate and its conformational change, which opens the coordination site trans to His462, (4) binding of a neutral XOH ligand (H(2)O or Tyr447) at the open site, (5) proton transfer from XOH to the neighboring peroxo ligand yielding the hydroperoxo intermediate, (6) a Criegee rearrangement leading to the anhydride intermediate, and (7) hydrolysis of the anhydride to the final acyclic product. One of the most important results obtained is that the Criegee mechanism requires an in-plane orientation of the four atoms (two oxygen and two carbon atoms) mainly involved in the reaction. This orientation yields a good overlap between the two sigma orbitals involved, C-C sigma and O-O sigma, allowing an efficient electron flow between them. Another interesting result is that under some conditions, a homolytic O-O bond cleavage might compete with the Criegee rearrangement. The role of the second-shell residues and the substituent effects are also discussed.