This article describes the software suite GROMACS (Groningen MAchine for Chemical Simulation) that was developed at the University of Groningen, The Netherlands, in the early 1990s. The software, written in ANSI C, originates from a parallel hardware project, and is well suited for parallelization on processor clusters. By careful optimization of neighbor searching and of inner loop performance, GROMACS is a very fast program for molecular dynamics simulation. It does not have a force field of its own, but is compatible with GROMOS, OPLS, AMBER, and ENCAD force fields. In addition, it can handle polarizable shell models and flexible constraints. The program is versatile, as force routines can be added by the user, tabulated functions can be specified, and analyses can be easily customized. Nonequilibrium dynamics and free energy determinations are incorporated. Interfaces with popular quantum-chemical packages (MOPAC, GAMES-UK, GAUSSIAN) are provided to perform mixed MM/QM simulations. The package includes about 100 utility and analysis programs. GROMACS is in the public domain and distributed (with source code and documentation) under the GNU General Public License. It is maintained by a group of developers from the Universities of Groningen, Uppsala, and Stockholm, and the Max Planck Institute for Polymer Research in Mainz. Its Web site is http://www.gromacs.org.

/pdf/gromacs-fast-flexible-and-free-15aazpel5u.pdf

GROMACS: Fast, flexible, and free

We present the theoretical and technical foundations of the Amsterdam Density Functional (ADF) program with a survey of the characteristics of the code (numerical integration, density fitting for the Coulomb potential, and STO basis functions). Recent developments enhance the efficiency of ADF (e.g., parallelization, near order-N scaling, QM/MM) and its functionality (e.g., NMR chemical shifts, COSMO solvent effects, ZORA relativistic method, excitation energies, frequency-dependent (hyper)polarizabilities, atomic VDD charges). In the Applications section we discuss the physical model of the electronic structure and the chemical bond, i.e., the Kohn–Sham molecular orbital (MO) theory, and illustrate the power of the Kohn–Sham MO model in conjunction with the ADF-typical fragment approach to quantitatively understand and predict chemical phenomena. We review the “Activation-strain TS interaction” (ATS) model of chemical reactivity as a conceptual framework for understanding how activation barriers of various types of (competing) reaction mechanisms arise and how they may be controlled, for example, in organic chemistry or homogeneous catalysis. Finally, we include a brief discussion of exemplary applications in the field of biochemistry (structure and bonding of DNA) and of time-dependent density functional theory (TDDFT) to indicate how this development further reinforces the ADF tools for the analysis of chemical phenomena. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 931–967, 2001

Chemistry with ADF

▪ Abstract Proton-coupled electron transfer (PCET) is an important mechanism for charge transfer in a wide variety of systems including biology- and materials-oriented venues. We review several are...

/pdf/proton-coupled-electron-transfer-3oyvp0w9qd.pdf

Proton-Coupled Electron Transfer

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

QM/MM Methods for Biomolecular Systems

ion phase that leads to an alkyl radical coordinated to the iron-hydroxo complex by a weak OH---C hydrogen bond, labeled as CI; (ii) an alkyl (or OH) rotation phase whereby the alkyl group achieves a favorable orientation for rebound; and (iii) a rebound phase that leads to C-O bond making and the ferric-alcohol complexes, 4,2P. The two profiles remain close in energy throughout the first two phases and then bifurcate. Whereas the HS state exhibits a significant barrier and a genuine TS for rebound, in the LS state, once the right orientation of the alkyl group is achieved, the LS rebound proceeds in a virtually barrier-free fashion to the alcohol. As such, alkane hydroxylation proceeds by TSR,70-72,120 in which the HS mechanism is truly stepwise with a finite lifetime for the radical intermediate, whereas the LS mechanism is effectively concerted with an ultrashort lifetime for the radical intermediate. Subsequent studies of ethane and camphor hydroxylation by the Yoshizawa group117,181-183 arrived at basically the same conclusion, that the mechanism is typified by TSR. The differences between the results of Shaik et al.130,173,177-180 and Yoshizawa et al.117,181-183 were rationalized recently71,72 and shown to arise owing to technical problems and the choice of the mercaptide ligand,117,181-183 which is a powerful electron donor and is too far from the representation of cysteine in the protein environment. The most recent study of camphor hydroxylation, which was done at a higher quality,117 converged to the picture reported by Shaik et al.130,173,177-180 and shows a stepwise HS process with a barrier of more than 3 kcal/mol for C-O bond formation by rebound of the camphoryl radical vis-à-vis an effectively concerted LS process for which this barrier is 0.7 kcal mol-1 and is the rotational barrier for reaching the rebound position. By referring to Figure 21, it is possible to rationalize the clock data of Newcomb in a simple manner. The apparent lifetimes are based on the assumption that there is a single state that leads to the reaction, such that the radical lifetime can be quantitated from the rate constant of free radical rearrangement and the ratio of rearranged to unrearranged alcohol product. However, in TSR, the rearranged (R) product is formed only/mainly on the HS surface, while the unrearranged (U) product is formed mainly on Figure 20. Formal descriptions of iron(III)-peroxo, iron(III)-hydroperoxo, and iron(V)-oxo species with indication of the negative charges. The roles “electrophile” or “nucleophile” are assigned according to the charge type. Reprinted with permission from ref 7. Copyright 2000 Springer-Verlag Heidelberg. 3964 Chemical Reviews, 2004, Vol. 104, No. 9 Meunier et al.

Mechanism of oxidation reactions catalyzed by cytochrome p450 enzymes.

Theoretical perspective on the structure and mechanism of cytochrome P450 enzymes.

Two-state reactivity mechanisms of hydroxylation and epoxidation by cytochrome P-450 revealed by theory.

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.; Schroder, 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.

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

Epoxidation (C=C) vis-a-vis allylic hydroxylation (C-H) reactions of propene with a model compound I (Cpd I) of the enzyme cytochrome P450 were studied using B3LYP density functional theory. Potential energy profiles and kinetic isotope effects (KIE) were calculated. The interactions in the protein pocket were mimicked by adding two external NH- - -S hydrogen bonds to the thiolate ligand and by introducing a nonpolar medium (with a dielectric constant, epsilon = 5.7) that can exert a polarization effect on the reacting species. A two-state reactivity (TSR) with high-spin (HS) and low-spin (LS) states were located for both processes (Ogliaro, F.; Harris, N.; Cohen, S.; Filatov, M.; de Visser, S. P.; Shaik, S. J. Am. Chem. Soc. 2000, 122, 8977-8989. de Visser, S. P.; Ogilaro, F.; Harris, N.; Shaik, S. J. Am. Chem. Soc. 2001, 123, 3037-3047). The HS processes were found to be stepwise, whereas the LS processes were characterized as nonsynchronous but effectively concerted pathways. The computed KIE for C-H hydroxylation with and without tunneling corrections are large (>7), and they support the assignment of the corresponding transition states as hydrogen-abstraction species (Groves, J. T.; Han, Y.-Z. In Cytochrome P450: Structures, Mechanism and Biochemistry, 2nd ed.; Ortiz de Montellano, P. R., Ed.; Plenum Press: New York, 1995; Chapter 1; pp 3-48). In the gas phase, epoxidation is energetically favorable by 3.4 kcal mol(-1). Inclusion of zero-point energy reduces this difference but still predicts C=C/C-H > 1. Environmental effects were found to have major impact on the C=C/C-H ratio as well as on the stereoselectivity of the processes. Thus, two NH- - -S hydrogen bonds away from the reaction center reverse the regioselectivity and prefer hydroxylation, namely, C=C/C-H <1. The polarity of the medium further accentuates the trend and leads to a change by 2 orders of magnitude in the regioselectivity, C=C/C-H << 1. Furthermore, since the environmental interactions prefer the LS over the HS reactions, both hydroxylation and epoxidation processes are rendered more stereoselective, again by 2 orders of magnitude. It follows, therefore, that Cpd I is a chameleon oxidant (Ogliaro, F.; Cohen, S.; de Visser, S. P.; Shaik, S. J. Am. Chem. Soc. 2000, 122, 12892-12893; Ogliaro, F.; de Visser, S. P.; Cohen, S.; Kaneti, J.; Shaik, S. Chembiochem. 2001, 2, 848-851; Ogliaro, F.; de Visser, S. P.; Groves, J. T.; Shaik, S. Angew. Chem., Int. Ed. 2001, 40, 2874-2878) that tunes its reactivity and selectivity patterns in response to the protein environment in which it is accommodated. A valence bond (VB) model, akin to "redox mesomerism" (Bernadou, J.; Fabiano, A.-S.; Robert, A.; Meunier, B. J. Am. Chem. Soc. 1994, 116, 9375-9376), is constructed and enables the description of a chameleon transition state. It shows that the good donor ability of the thiolate ligand and the acceptor ability of the iron porphyrin create mixed-valent situations that endow the transition state with a great sensitivity to external perturbations as in the protein pocket. The model is used to discuss the computed results and to relate them to experimental findings.

What factors affect the regioselectivity of oxidation by cytochrome P450? A DFT study of allylic hydroxylation and double bond epoxidation in a model reaction

Sam P. de Visser

Papers

Mechanism of oxidation reactions catalyzed by cytochrome p450 enzymes.

Theoretical perspective on the structure and mechanism of cytochrome P450 enzymes.

Two-state reactivity mechanisms of hydroxylation and epoxidation by cytochrome P-450 revealed by theory.

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

What factors affect the regioselectivity of oxidation by cytochrome P450? A DFT study of allylic hydroxylation and double bond epoxidation in a model reaction