Mechanistic Insights into a Stibene Cleavage Oxygenase NOV1 from Quantum Mechanical/Molecular Mechanical Calculations

NOV1, a stilbene cleavage oxygenase, catalyzes the cleavage of the central double bond of stilbenes to two phenolic aldehydes, using a 4-His Fe(II) center and dioxygen. Herein, we use in-protein quantum mechanical/molecular mechanical (QM/MM) calculations to elucidate the reaction mechanism of the central double bond cleavage of phytoalexin resveratrol by NOV1. Our results showed that the oxygen molecule prefers to bind to the iron center in a side-on fashion, as suggested from the experiment. The quintet Fe−O₂ complex with the side-on superoxo antiferromagnetic coupled to the resveratrol radical is identified as the reactive oxygen species. The QM/MM results support the dioxygenase mechanism involving a dioxetane intermediate with a rate-limiting barrier of 10.0 kcal mol⁻¹. The alternative pathway through an epoxide intermediate is ruled out due to a larger rate-limiting barrier (26.8 kcal mol⁻¹). These findings provide important insight into the catalytic mechanism of carotenoid cleavage oxygenases and also the dioxygen activation of non-heme enzymes.     

A proton-shuttle mechanism mediated by the porphyrin in benzene hydroxylation by cytochrome p450 enzymes

Benzene hydroxylation is a fundamental process in chemical catalysis. In nature, this reaction is catalyzed by the enzyme cytochrome P450 via oxygen transfer in a still debated mechanism of considerable complexity. The paper uses hybrid density functional calculations to elucidate the mechanisms by which benzene is converted to phenol, benzene oxide, and ketone, by the active species of the enzyme, the high-valent iron-oxo porphyrin species. The effects of the protein polarity and hydrogen-bonding donation to the active species are mimicked, as before (Ogliaro, F.; Cohen, S.; de Visser, S. P.; Shaik, S. J. Am. Chem. Soc. 2000, 122, 12892-12893). It is verified that the reaction does not proceed either by hydrogen abstraction or by initial electron transfer (Ortiz de Montellano, P. R. In Cytochrome P450: Structure, Mechanism and Biochemistry, 2nd ed.; Ortiz de Montellano, P. R., Ed.; Plenum Press: New York, 1995; Chapter 8, pp 245-303). In accord with the latest experimental conclusions, the theoretical calculations show that the reactivity is an interplay of electrophilic and radicalar pathways, which involve an initial attack on the pi-system of the benzene to produce sigma-complexes (Korzekwa, K. R.; Swinney, D. C.; Trager, W. T. Biochemistry 1989, 28, 9019-9027). The dominant reaction channel is electrophilic and proceeds via the cationic sigma-complex,( 2)3, that involves an internal ion pair made from a cationic benzene moiety and an anionic iron porphyrin. The minor channel proceeds by intermediacy of the radical sigma-complex, (2)2, in which the benzene moiety is radicalar and the iron-porphyrin moiety is neutral. Ring closure in these intermediates produces the benzene oxide product ((2)4), which does not rearrange to phenol ((2)7) or cyclohexenone ((2)6). While such a rearrangement can occur post-enzymatically under physiological conditions by acid catalysis, the computations reveal a novel mechanism whereby the active species of the enzyme catalyzes directly the production of phenol and cyclohexenone. This enzymatic mechanism involves proton shuttles mediated by the porphyrin ring through the N-protonated intermediate, (2)5, which relays the proton either to the oxygen atom to form phenol ((2)7) or to the ortho-carbon atom to produce cyclohexenone product ((2)6). The formation of the phenol via this proton-shuttle mechanism will be competitive with the nonenzymatic conversion of benzene oxide to phenol by external acid catalysis. With the assumption that (2)5 is not fully thermalized, this novel mechanism would account also for the observation that there is a partial skeletal retention of the original hydrogen of the activated C-H bond, due to migration of the hydrogen from the site of hydroxylation to the adjacent carbon (so-called "NIH shift" (Jerina, D. M.; Daly, J. W. Science 1974, 185, 573-582)). Thus, in general, the computationally discovered mechanism of a porphyrin proton shuttle suggests thatthere is an enzymatic pathway that converts benzene directly to a phenol and ketone, in addition to nonenzymatic production of these species by conversion of arene oxide to phenol and ketone. The potential generality of protonated porphyrin intermediates in P450 chemistry is discussed in the light of the H/D exchange observed during some olefin epoxidation reactions (Groves, J. T.; Avaria-Neisser, G. E.; Fish, K. M.; Imachi, M.; Kuczkowski, R. J. Am. Chem. Soc. 1986, 108, 3837-3838) and the general observation of heme alkylation products (Kunze, K. L.; Mangold, B. L. K.; Wheeler, C.; Beilan, H. S.; Ortiz de Montellano, P. R. J. Biol. Chem. 1983, 258, 4202-4207). The competition, similarities, and differences between benzene oxidation viz. olefin epoxidation and alkanyl C-H hydroxylation are discussed, and comparison is made with relevant experimental and computational data. The dominance of low-spin reactivity in benzene hydroxylation viz. two-state reactivity (Shaik, S.; de Visser, S. P.; Ogliaro, F.; Schwarz, H.; Schr?der, D. Curr. Opin. Chem. Biol. 2002, 6, 556-567) in olefin epoxidation and alkane hydroxylation is traced to the loss of benzene resonance energy during the bond activation step.     

Electronic structure of compound I in human isoforms of cytochrome P450 from QM/MM modeling

Human cytochromes P450 play a vital role in drug metabolism. The key step in substrate oxidation involves hydrogen atom abstraction or C=C bond addition by the oxygen atom of the Compound I intermediate. The latter has three unpaired electrons, two on the Fe-O center and one shared between the porphyrin ring and the proximal cysteinyl sulfur atom. Changes in its electronic structure have been suggested to affect reactivity. The electronic and geometric structure of Compound I in three important human subfamilies of cytochrome P450 (P450, 2C, 2B, and 3A) that are major contributors to drug metabolism is characterized here using combined quantum mechanical/molecular mechanical (QM/MM) calculations at the B3LYP:CHARMM27 level. Compound I is remarkably similar in all isoforms, with the third unpaired electron located mainly on the porphyrin ring, and this prediction is not very sensitive to details of the QM/MM methodology, such as the DFT functional, the basis set, or the size of the QM region. The presence of substrate also has no effect. The main source of variability in spin density on the cysteinyl sulfur (from 26 to 50%) is the details of the system setup, such as the starting protein geometry used for QM/MM minimization. This conformational effect is larger than the differences between human isoforms, which are therefore not distinguishable on electronic grounds, so it is unlikely that observed large differences in substrate selectivity can be explained to a large extent in these terms.     

Systematic QM/MM investigation of factors that affect the cytochrome P450-catalyzed hydrogen abstraction of camphor

The hydrogen abstraction reaction of camphor in cytochrome P450(cam) has been investigated in the native enzyme environment by combined quantum mechanical/molecular mechanical (QM/MM) calculations and in the gas phase by density functional calculations. This work has been motivated by contradictory published QM/MM results. In an attempt to pinpoint the origin of these discrepancies, we have systematically studied the factors that may affect the computed barriers, including the QM/MM setup, the optimization procedures, and the choice of QM region, basis set, and protonation states. It is found that the ChemShell and QSite programs used in the published QM/MM calculations yield similar results at given geometries, and that the discrepancies mainly arise from two technical issues (optimization protocols and initial system preparation) that need to be well controlled in QM/MM work. In the course of these systematic investigations, new mechanistic insights have been gained. The crystallographic water 903 placed near the oxo atom of Compound I lowers the hydrogen abstraction barrier by ca. 4 kcal/mol, and thus acts as a catalyst for this reaction. Spin density may appear at the A-propionate side chain of the heme if the carboxylate group is not properly screened, which might be expected to happen during protein dynamics, but not in static equilibrium situations. There is no clear correlation between the computed A-propionate spin density and the hydrogen abstraction barrier, and hence, no support for a previously proposed side-chain mediated transition state stabilization mechanism. Standard QM/MM optimizations yield an A-propionate environment close to the X-ray structure only for protonated Asp297, and not for deprotonated Asp297, but the computed barriers are similar in both cases. An X-ray like A-propionate environment can also be obtained when deprotonated Asp297 is included in the QM region and His355 is singly protonated, but this Compound II-type species with a closed-shell porphyrin ring has a higher hydrogen abstraction barrier and should thus not be mechanistically relevant.     

The elusive oxidant species of cytochrome P450 enzymes: characterization by combined quantum mechanical/molecular mechanical (QM/MM) calculations

The primary oxidant of cytochrome P450 enzymes, Compound I, is hard to detect experimentally; in the case of cytochrome P450(cam), this intermediate does not accumulate in solution during the catalytic cycle even at temperatures as low as 200 K (ref 4). Theory can play an important role in characterizing such elusive species. We present here combined quantum mechanical/molecular mechanical (QM/MM) calculations of Compound I of cytochrome P450(cam) in the full enzyme environment as well as density functional studies of the isolated QM region. The calculations assign the ground state of the species, quantify the effect of polarization and hydrogen bonding on its properties, and show that the protein environment and its specific hydrogen bonding to the cysteinate ligand are crucial for sustaining the Fe-S bond and for preventing the full oxidation of the sulfur.     

Can ferric-superoxide act as a potential oxidant in P450(cam)? QM/MM investigation of hydroxylation, epoxidation, and sulfoxidation

In view of recent reports of high reactivity of ferric-superoxide species in heme and nonheme systems (Morokuma et al. J. Am. Chem. Soc. 2010, 132, 11993-12005; Que et al. Inorg. Chem. 2010, 49, 3618-3628; Nam et al. J. Am. Chem. Soc. 2010, 132, 5958-5959; J. Am. Chem. Soc. 2010, 132, 10668-10670), we use herein combined quantum mechanics/molecular mechanics (QM/MM) methods to explore the potential reactivity of P450(cam) ferric-superoxide toward hydroxylation, epoxidation, and sulfoxidation. The calculations demonstrate that P450 ferric-superoxide is a sluggish oxidant compared with the high-valent oxoiron porphyrin cation-radical species. As such, unlike heme enzymes with a histidine axial ligand, the P450 superoxo species does not function as an oxidant in P450(cam). The origin of this different behavior of the superoxo species of P450 vis-a?-vis other heme enzymes like tryptophan 2, 3-dioxygenase (TDO) is traced to the ability of the latter superoxo species to make a stronger FeOO-X (X = H,C) bond and to stabilize the corresponding bond-activation transition states by resonance with charge-transfer configurations. By contrast, the negatively charged thiolate ligand in the P450 superoxo species minimizes the mixing of charge transfer configurations in the transition state and raises the reaction barrier. However, as we demonstrate, an external electric field oriented along the Fe-O axis with a direction pointing from Fe toward O will quench Cpd I formation by slowing the reduction of ferric-superoxide and will simultaneously lower the barriers for oxidation by the latter species, thereby enabling observation of superoxo chemistry in P450. Other options for nascent superoxo reactivity in P450 are discussed.     

QM/MM modeling of compound I active species in cytochrome P450, cytochrome C peroxidase, and ascorbate peroxidase

QM/MM calculations provide a means for predicting the electronic structure of the metal center in metalloproteins. Two heme peroxidases, Cytochrome c Peroxidase (CcP) and Ascorbate Peroxidase (APX), have a structurally very similar active site, yet have active intermediates with very different electronic structures. We review our recent QM/MM calculations on these systems, and present new computational data. Our results are in good agreement with experiment, and suggest that the difference in electronic structure is due to a large number of small differences in structure from one protein to another. We also discuss recent QM/MM calculations on the active species of cytochrome P450, in which a similar sensitivity of the electronic structure to the environment is found. However, this does not appear to explain different catalytic profiles of the different drug-metabolizing isoforms of this class of enzyme.     

A new understanding on how heme metabolism occurs in heme oxygenase: water-assisted oxo mechanism

Heme metabolism by heme oxygenase (HO) is investigated with quantum mechanical/molecular mechanical (QM/MM) calculations. A mechanism assisted by water is proposed: (1) an iron-oxo species and a water molecule are generated by the heterolytic cleavage of the O-O bond of an iron-hydroperoxo species in a similar way to P450-mediated reactions, (2) a hydrogen atom abstraction by the iron-oxo species from the generated water molecule and the C-O bond formation between the water molecule and the α-meso carbon take place simultaneously. The water molecule is hydrogen-bonded to the oxo ligand and to the water cluster in the active site of HO. The water cluster can control the position of the generated water molecule to ensure the regioselective oxidation of heme at the α-meso position, at the same time, can facilitate the oxidation by stabilizing a positive charge on the water molecule in the transition state. A key difference between HO and P450 is observed in the structure of the active site; Thr252 in P450 blocks the access of the water molecule to the α-meso position, and can thus suppress the undesired heme oxidation for P450.     

QM/MM study of the taxadiene synthase mechanism

Combined quantum mechanics/molecular mechanics (QM/MM) calculations were used to investigate the reaction mechanism of taxadiene synthase (TXS). TXS catalyzes the cyclization of geranylgeranyl diphosphate (GGPP) to taxadiene (T) and four minor cyclic products. All these products originate from the deprotonation of carbocation intermediates. The reaction profiles for the conversion of GGPP to T as well as to minor products were calculated for different configurations of relevant TXS carbocation complexes. The QM region was treated at the M06-2X/TZVP level, while the CHARMM27 force field was used to describe the MM region. The QM/MM calculations suggest a reaction pathway for the conversion of GGPP to T, which slightly differs from previous proposals regarding the number of reaction steps and the conformation of the carbocations. The QM/MM results also indicate that the formation of minor products via water-assisted deprotonation of the carbocations is highly exothermic, by about −7 to −23 kcal/mol. Curiously, however, the computed barriers and reaction energies indicate that the formation of some of the minor products is more facile than the formation of T. Thus, the present QM/MM calculations provide detailed insights into possible reaction pathways and into the origin of the promiscuity of TXS, but they do not reproduce the product distribution observed experimentally. © 2019 Wiley Periodicals, Inc.     

Quantum mechanically derived AMBER-compatible heme parameters for various states of the cytochrome P450 catalytic cycle

Molecular mechanics (MM) methods are computationally affordable tools for screening chemical libraries of novel compounds for sites of P450 metabolism. One challenge for MM methods has been the absence of a consistent and transferable set of parameters for the heme within the P450 active site. Experimental data indicate that mammalian P450 enzymes vary greatly in the size, architecture, and plasticity of their active sites. Thus, obtaining X-ray-based geometries for the development of accurate MM parameters for the major classes of hepatic P450 remains a daunting task. Our previous work with preliminary gas-phase quantum mechanics (QM)-derived atomic partial charges greatly improved the accuracy of docking studies of raloxifene to CYP3A4. We have therefore developed and tested a consistent set of transferable MM parameters based on gas-phase QM calculations of two model systems of the heme-a truncated (T-HM) and a full (F-HM) for four states of the P450 catalytic cycle. Our results indicate that the use of the atomic partial charges from the F-HM further improves the accuracy of docked predictions for raloxifene to CYP3A4. Different patterns for substrate docking are also observed depending on the choice of heme model and state. Newly parameterized heme models are tested in implicit and explicitly solvated MD simulations in the absence and presence of enzyme structures, for CYP3A4, and appear to be stable on the nanosecond simulation timescale. The new force field for the various heme states may aid the community for simulations of P450 enzymes and other heme-containing enzymes.     

Multistate reactivity in styrene epoxidation by compound I of cytochrome p450: mechanisms of products and side products formation

Density functional theoretical calculations are used to elucidate the epoxidation mechanism of styrene with a cytochrome P450 model Compound I, and the formation of side products. The reaction features multistate reactivity (MSR) with different spin states (doublet and quartet) and different electromeric situations having carbon radicals and cations, as well as iron(III) and iron(IV) oxidation states. The mechanisms involve state-specific product formation, as follows: a) The low-spin pathways lead to epoxide formation in effectively concerted mechanisms. b) The high-spin pathways have finite barriers for ring-closure and may have a sufficiently long lifetime to undergo rearrangement and lead to side products. c) The high-spin radical intermediate, (4)2(rad)-IV, has a ring closure barrier as small as the C--C rotation barrier. This intermediate will therefore lose stereochemistry and lead to a mixture of cis and trans epoxides. The barriers for the production of aldehyde and suicidal complexes are too high for this intermediate. d) The high-spin radical intermediate, (4)2(rad)-III, has a substantial ring closure barrier and may survive long enough time to lead to suicidal, phenacetaldehyde and 2-hydroxostyrene side products. e) The phenacetaldehyde and 2-hydroxostyrene products both originate from crossover from the (4)2(rad)-III radical intermediate to the cationic state, (4)2(cat,z(2) ). The process involves an N-protonated porphyrin intermediate that re-shuttles the proton back to the substrate to form either phenacetaldehyde or 2-hydroxostyrene products. This resembles the internally mediated NIH-shift observed during benzene hydroxylation.     

Quantum mechanical/molecular mechanical study on the mechanisms of compound I formation in the catalytic cycle of chloroperoxidase: an overview on heme enzymes

The formation of Compound I (Cpd I), the active species of the enzyme chloroperoxidase (CPO), was studied using QM/MM calculation. Starting from the substrate complex with hydrogen peroxide, FeIII-HOOH, we examined two alternative mechanisms on the three lowest spin-state surfaces. The calculations showed that the preferred pathway involves heterolytic O-O cleavage that proceeds via the iron hydroperoxide species, i.e., Compound 0 (Cpd 0), on the doublet-state surface. This process is effectively concerted, with a barrier of 12.4 kcal/mol, and is catalyzed by protonation of the distal OH group of Cpd 0. By comparison, the path that involves a direct O-O cleavage from FeIII-HOOH is less favored. A proton coupled electron transfer (PCET) feature was found to play an important role in the mechanism nascent from Cpd 0. Initially, the O-O cleavage progresses in a homolytic sense, but as soon as the proton is transferred to the distal OH, it triggers an electron transfer from the heme-oxo moiety to form water and Cpd I. This study enables us to generalize the mechanisms of O-O activation, elucidated so far by QM/MM calculations, for other heme enzymes, e.g., cytochrome P450cam, horseradish peroxidase (HRP), nitric oxide synthase (NOS), and heme oxygenase (HO). Much like for CPO, in the cases of P450 and HRP, the PCET lowers the barrier below the purely homolytic cleavage alternative (in our case, the homolytic mechanism is calculated directly from FeIII-HOOH). By contrast, the absence of PCET in HO, along with the robust water cluster, prefers a homolytic cleavage mechanism.     

Multireference ab initio quantum mechanics/molecular mechanics study on intermediates in the catalytic cycle of cytochrome P450(cam)

We have investigated the elusive reactive species of cytochrome P450(cam) (Compound I), the hydroxo complex formed during camphor hydroxylation, and the ferric hydroperoxo complex (Compound 0) by combined quantum mechanical/molecular mechanical (QM/MM) calculations, employing both density functional theory (DFT) and correlated ab initio methods. The first two intermediates appear multiconfigurational in character, especially in the doublet state and less so in the quartet state. DFT(B3LYP)/MM calculations reproduce the relative energies from correlated ab initio QM/MM treatments quite well, except for the splitting of the lowest A(1u)-A(2u) radical states. The inclusion of dynamic correlation is crucial for the proper ab initio treatment of these intermediates.     

Amino acid analogues bind to carbon nanotube via π-π interactions: comparison of molecular mechanical and quantum mechanical calculations

Understanding the interaction between carbon nanotubes (CNTs) and biomolecules is essential to the CNT-based nanotechnology and biotechnology. Some recent experiments have suggested that the π-π stacking interactions between protein's aromatic residues and CNTs might play a key role in their binding, which raises interest in large scale modeling of protein-CNT complexes and associated π-π interactions at atomic detail. However, there is concern on the accuracy of classical fixed-charge molecular force fields due to their classical treatments and lack of polarizability. Here, we study the binding of three aromatic residue analogues (mimicking phenylalanine, tyrosine, and tryptophan) and benzene to a single-walled CNT, and compare the molecular mechanical (MM) calculations using three popular fixed-charge force fields (OPLSAA, AMBER, and CHARMM), with quantum mechanical (QM) calculations using the density-functional tight-binding method with the inclusion of dispersion correction (DFTB-D). Two typical configurations commonly found in π-π interactions are used, one with the aromatic rings parallel to the CNT surface (flat), and the other perpendicular (edge). Our calculations reveal that compared to the QM results the MM approaches can appropriately reproduce the strength of π-π interactions for both configurations, and more importantly, the energy difference between them, indicating that the various contributions to π-π interactions have been implicitly included in the van der Waals parameters of the standard MM force fields. Meanwhile, these MM models are less accurate in predicting the exact structural binding patterns (matching surface), meaning there are still rooms to be improved. In addition, we have provided a comprehensive and reliable QM picture for the π-π interactions of aromatic molecules with CNTs in gas phase, which might be used as a benchmark for future force field developments.     

Water-assisted oxo mechanism for heme metabolism

A mechanism of heme metabolism by heme oxygenase (HO) is discussed from B3LYP density functional theory calculations. The concerted OH group attack to the alpha-carbon by the iron-hydroperoxo species is investigated using a model with full protoporphyrin IX to confirm our previous conclusion that this species does not have sufficient oxidizing power for heme oxidation (J. Am. Chem. Soc. 2004, 126, 3672). Calculated activation energies and structures of the intermediates and transition state for this process remain unchanged from those for a small model with porphine in the previous study, which shows that the inclusion of the side chain of the porphyrin ring is not essential in describing the OH group transfer. The activation barrier for a direct oxo attack to the alpha-carbon by an iron-oxo model is calculated to be 49.8 kcal/mol, the barrier height of which looks very high for the enzymatic reaction under physiological conditions. This large activation energy is due to a highly bent porphyrin structure in the transition state. However, a bridging water molecule plays an important role in reducing the porphyrin distortion in the transition state, resulting in a remarkable decrease of the activation barrier to 13.9 kcal/mol. A whole-enzyme model with about 4000 atoms is constructed to elucidate functions of the protein environment in this enzymatic reaction using QM/MM calculations. The key water molecule is fixed in the protein environment to ensure the low-barrier and regioselective heme oxidation. A water-assisted oxo mechanism of heme oxidation by heme oxygenase is proposed from these calculational results.     

Geometry optimization based on linear response free energy with quantum mechanical/molecular mechanical method: applications to Menshutkin-type and Claisen rearrangement reactions in aqueous solution

The authors present a method based on a linear response theory that allows one to optimize the geometries of quantum mechanical/molecular mechanical (QM/MM) systems on the free energy surfaces. Two different forms of linear response free energy functionals are introduced, and electronic wave functions of the QM region, as well as the responses of electrostatic and Lennard-Jones potentials between QM and MM regions, are self-consistently determined. The covariant matrix relating the QM charge distribution to the MM response is evaluated by molecular dynamics (MD) simulation of the MM system. The free energy gradients with respect to the QM atomic coordinates are also calculated using the MD trajectory results. They apply the present method to calculate the free energy profiles of Menshutkin-type reaction of NH3 with CH3Cl and Claisen rearrangement of allyl vinyl ether in aqueous solution. For the Menshutkin reaction, the free energy profile calculated with the modified linear response free energy functional is in good agreement with that by the free energy perturbation calculations. They examine the nonequilibrium solvation effect on the transmission coefficient and the kinetic isotope effect for the Claisen rearrangement.     

Quantum mechanical/molecular mechanical study on the mechanism of the enzymatic Baeyer-Villiger reaction

We report a combined quantum mechanical/molecular mechanical (QM/MM) study on the mechanism of the enzymatic Baeyer-Villiger reaction catalyzed by cyclohexanone monooxygenase (CHMO). In QM/MM geometry optimizations and reaction path calculations, density functional theory (B3LYP/TZVP) is used to describe the QM region consisting of the substrate (cyclohexanone), the isoalloxazine ring of C4a-peroxyflavin, the side chain of Arg-329, and the nicotinamide ring and the adjacent ribose of NADP(+), while the remainder of the enzyme is represented by the CHARMM force field. QM/MM molecular dynamics simulations and free energy calculations at the semiempirical OM3/CHARMM level employ the same QM/MM partitioning. According to the QM/MM calculations, the enzyme-reactant complex contains an anionic deprotonated C4a-peroxyflavin that is stabilized by strong hydrogen bonds with the Arg-329 residue and the NADP(+) cofactor. The CHMO-catalyzed reaction proceeds via a Criegee intermediate having pronounced anionic character. The initial addition reaction has to overcome an energy barrier of about 9 kcal/mol. The formed Criegee intermediate occupies a shallow minimum on the QM/MM potential energy surface and can undergo fragmentation to the lactone product by surmounting a second energy barrier of about 7 kcal/mol. The transition state for the latter migration step is the highest point on the QM/MM energy profile. Gas-phase reoptimizations of the QM region lead to higher barriers and confirm the crucial role of the Arg-329 residue and the NADP(+) cofactor for the catalytic efficiency of CHMO. QM/MM calculations for the CHMO-catalyzed oxidation of 4-methylcyclohexanone reproduce and rationalize the experimentally observed (S)-enantioselectivity for this substrate, which is governed by the conformational preferences of the corresponding Criegee intermediate and the subsequent transition state for the migration step.     

Performance assessment of semiempirical molecular orbital methods in describing halogen bonding: quantum mechanical and quantum mechanical/molecular mechanical-molecular dynamics study

The performance of semiempirical molecular-orbital methods--MNDO, MNDO-d, AM1, RM1, PM3 and PM6--in describing halogen bonding was evaluated, and the results were compared with molecular mechanical (MM) and quantum mechanical (QM) data. Three types of performance were assessed: (1) geometrical optimizations and binding energy calculations for 27 halogen-containing molecules complexed with various Lewis bases (Two of the tested methods, AM1 and RM1, gave results that agree with the QM data.); (2) charge distribution calculations for halobenzene molecules, determined by calculating the solvation free energies of the molecules relative to benzene in explicit and implicit generalized Born (GB) solvents (None of the methods gave results that agree with the experimental data.); and (3) appropriateness of the semiempirical methods in the hybrid quantum-mechanical/molecular-mechanical (QM/MM) scheme, investigated by studying the molecular inhibition of CK2 protein by eight halobenzimidazole and -benzotriazole derivatives using hybrid QM/MM molecular-dynamics (MD) simulations with the inhibitor described at the QM level by the AM1 method and the rest of the system described at the MM level. The pure MM approach with inclusion of an extra point of positive charge on the halogen atom approach gave better results than the hybrid QM/MM approach involving the AM1 method. Also, in comparison with the pure MM-GBSA (generalized Born surface area) binding energies and experimental data, the calculated QM/MM-GBSA binding energies of the inhibitors were improved by replacing the G(GB,QM/MM) solvation term with the corresponding G(GB,MM) term.