Multireference Ab Initio Quantum Mechanics/Molecular Mechanics Study on Intermediates in the Catalytic Cycle of Cytochrome P450cam

Exploring the Bonding Nature of Iron(IV)-Oxo Species through Valence Bond Calculations and Electron Density Analysis

Compound I (Cpd I) plays a pivotal role in substrate transformations within the catalytic cycle of cytochrome P450 enzymes (P450s). A key constituent of Cpd I is the iron(IV)-oxo unit, a structural motif also found in other heme enzymes and nonheme enzymes. In this study, we performed ab initio valence bond (VB) calculations, employing the valence bond self-consistent field (VBSCF) and breathing orbital valence bond (BOVB) methods, to unveil the bonding nature of this vital "Fe(IV)═O″ unit in bioinorganic chemistry. Comparisons were drawn with the triplet O2 molecule, which shares some electronic characteristics with iron(IV)-oxo. Additionally, Cpd I models of horseradish peroxidase (HRP) and catalase (CAT) were analyzed to assess the proximal ligand effect on the electronic structure of iron(IV)-oxo. Our VB analysis underscores the significant role of noncovalent resonance effects in shaping the iron(IV)-oxo bonding. The resonance stabilization within the π and σ frameworks occurs to comparable degrees, with additional stabilization resulting from resonance between VB structures from these frameworks. Furthermore, we elucidated the substantial influence of proximal and equatorial ligands in modulating the relative significance of different VB structures. Notably, in the presence of these ligands, iron(IV)-oxo is better described as iron(III)-oxyl or iron(II)-oxygen, displaying weak covalent character but enhanced by resonance effects. Although both species exhibit diradicaloid characters, resonance stabilization in iron(IV)-oxo is weaker than in O2. Further exploration using the Laplacian of electron density shows that, unlike O2, which exhibits a charge concentration region between its two oxygen atoms, iron(IV)-oxo species display a charge depletion region.

Multiscale Study on the Intramolecular C-S Bond Formation Catalyzed by P450 Monooxygenase CxnD Involved in the Biosynthesis of Chuangxinmycin: The Critical Roles of Noncrystal Water Molecule and Conformational Change

Cytochrome P450 monooxygenase CxnD catalyzes intramolecular C-S bond formation in the biosynthesis of chuangxinmycin, which is representative of the synthesis of sulfur-containing natural heterocyclic compounds. The intramolecular cyclization usually requires the activation of two reaction sites and a large conformational change; thus, illuminating its detailed reaction mechanism remains challengeable. Here, the reaction pathway of CxnD-catalyzed C-S bond formation was clarified by a series of calculations, including Gaussian accelerated molecular dynamics simulations and quantum mechanical-molecular mechanical calculations. Our results revealed that the C-S formation follows a diradical coupling mechanism. CxnD first employs Cpd I to abstract the hydrogen atom from the imino group of the indole ring, and then, the resulted Cpd II further extracts another hydrogen atom from the thiol group of the side chain to afford a diradical intermediate, in which a noncrystal water molecule entering into the active site after the formation of Cpd I was proved to play an indispensable role. Moreover, the diradical intermediate cannot directly perform the coupling reaction. It should first undergo a series of conformational changes leading to the proximity of two reaction sites. It is the flexibility of the active site of the enzyme and the side chain of the substrate that makes the diradical coupling to be successful.

Computational Characterization of the Reactivity of Compound I in Unspecific Peroxygenases

Unspecific peroxygenases (UPOs) are emerging as promising biocatalysts for selective oxyfunctionalization of unactivated C-H bonds. However, their potential in large-scale synthesis is currently constrained by suboptimal chemical selectivity. Improving the selectivity of UPOs requires a deep understanding of the molecular basis of their catalysis. Recent molecular simulations have sought to unravel UPO's selectivity and inform their design principles. However, most of these studies focused on substrate-binding poses. Few researchers have investigated how the reactivity of CpdI, the principal oxidizing intermediate in the catalytic cycle, influences selectivity in a realistic protein environment. Moreover, the influence of protein electrostatics on the reaction kinetics of CpdI has also been largely overlooked. To bridge this gap, we used multiscale simulations to interpret the regio- and enantioselective hydroxylation of the n-heptane substrate catalyzed by Agrocybe aegerita UPO (AaeUPO). We comprehensively characterized the energetics and kinetics of the hydrogen atom-transfer (HAT) step, initiated by CpdI, and the subsequent oxygen rebound step forming the product. Notably, our approach involved both free energy and potential energy evaluations in a quantum mechanics/molecular mechanics (QM/MM) setting, mitigating the dependence of results on the choice of initial conditions. These calculations illuminate the thermodynamics and kinetics of the HAT and oxygen rebound steps. Our findings highlight that both the conformational selection and the distinct chemical reactivity of different substrate hydrogen atoms together dictate the regio- and enantio-selectivity. Building on our previous study of CpdI's formation in AaeUPO, our results indicate that the HAT step is the rate-limiting step in the overall catalytic cycle. The subsequent oxygen rebound step is swift and retains the selectivity determined by the HAT step. We also pinpointed several polar and charged amino acid residues whose electrostatic potentials considerably influence the reaction barrier of the HAT step. Notably, the Glu196 residue is pivotal for both the CpdI's formation and participation in the HAT step. Our research offers in-depth insights into the catalytic cycle of AaeUPO, which will be instrumental in the rational design of UPOs with enhanced properties.

Reliably assessing the electronic structure of cytochrome P450 on today's classical computers and tomorrow's quantum computers

Chemical simulation is one of the most promising applications for future quantum computers. It is thought that quantum computers may enable accurate simulation for complex molecules that are otherwise impossible to simulate classically; that is, it displays quantum advantage. To better understand quantum advantage in chemical simulation, we explore what quantum and classical resources are required to simulate a series of pharmaceutically relevant molecules. Using classical methods, we show that reliable classical simulation of these molecules requires significant resources and therefore is a promising candidate for quantum simulation. We estimate the quantum resources, both in overall simulation time and the size. The insights from this study pave the way for future quantum simulation of complex molecules.

An accurate assessment of how quantum computers can be used for chemical simulation, especially their potential computational advantages, provides important context on how to deploy these future devices. To perform this assessment reliably, quantum resource estimates must be coupled with classical computations attempting to answer relevant chemical questions and to define the classical algorithms simulation frontier. Herein, we explore the quantum computation and classical computation resources required to assess the electronic structure of cytochrome P450 enzymes (CYPs) and thus define a classical–quantum advantage boundary. This is accomplished by analyzing the convergence of density matrix renormalization group plus n-electron valence state perturbation theory (DMRG+NEVPT2) and coupled-cluster singles doubles with noniterative triples [CCSD(T)] calculations for spin gaps in models of the CYP catalytic cycle that indicate multireference character. The quantum resources required to perform phase estimation using qubitized quantum walks are calculated for the same systems. Compilation into the surface code provides runtime estimates to compare directly to DMRG runtimes and to evaluate potential quantum advantage. Both classical and quantum resource estimates suggest that simulation of CYP models at scales large enough to balance dynamic and multiconfigurational electron correlation has the potential to be a quantum advantage problem and emphasizes the important interplay between classical computations and quantum algorithms development for chemical simulation.

Histidine versus Cysteine-Bearing Heme-Dependent Halogen Peroxidases: Parallels and Differences for Cl- Oxidation

The homodimeric myeloperoxidase (MPO) features a histidine as a proximal ligand and a sulfonium linkage covalently attaching the heme porphyrin ring to the protein. MPO is able to catalyze Cl^- oxidation with about the same efficiency as chloroperoxidase at pH 7.0. In this study, we seek to explore the parallels and differences between the histidine and cysteine heme-dependent halogen peroxidases. Transition states, reaction barriers, and relevant thermodynamic properties are calculated on protein models. Together with electronic structure calculations, it gives an overview of the reaction mechanisms and of the factors that determine the selectivity between one- and two-electron paths. Conclusions point to the innate oxidizing nature of MPO with the ester and sulfonium linkages hiking up the reactivity to enable chloride oxidation. The installation of a deprotonated imidazolate as a proximal ligand does not shift the equilibrium from one- to two-electron events without influencing the chemistry of the oxidation reaction.

Regio- and stereoselectivity in the CYP450BM3-catalyzed hydroxylation of complex terpenoids: a QM/MM study

The selective oxidation of C–H in artemisinin by P450_BM3 variants was disclosed by combining QM/MM and MD simulations.

Kinetic Isotope Effect Probes the Reactive Spin State, As Well As the Geometric Feature and Constitution of the Transition State during H-Abstraction by Heme Compound II Complexes

What do experimentally measured kinetic isotope effects (KIEs) tell us about H-abstraction reactions with multispin-state reactivity options? Using DFT calculations with tunneling corrections for experimentally studied H-abstraction reactions of porphyrin-Compound II species (Chem.-Eur. J. 2014, 20, 14437; Angew. Chem., Int. Ed. 2008, 47, 7321) with cyclohexane, dihydroanthracene (DHA), and xanthene (Xan), we show here that KIE is a selective probe that identifies the experimentally reactive spin state. At the same time, comparison of calculated and experimental KIE values permits us to determine the structural orientation of the transition states, as well as the presence/absence of an axial ligand, and the effect of porphyrin substituents. The studied compound II (Cpd II) species involve porphine, and porphyrin ligands with different meso-substituents, TPFPP (tetrakis(pentafluorophenyl)porphyrin dianion) and TMP (tetramesitylporphyrin dianion), with and without imidazole axial ligands. The DFT calculations reveal three potential pathways: quintet and triplet σ-pathways (⁵Hσ and ³Hσ) that possess linear transition state (TS) structures, and a triplet π -pathway (³Hπ) having a bent TS structure. Without an axial ligand, the ⁵Hσ pathways for these Cpd II complexes cross below the triplet states. The axial ligand raises the barriers for the quintet and triplet σ-pathways and quenches any chances for two-state reactivity, thus proceeding via the ³Hπ pathway. All of these pathways exhibit characteristic KIE values: very large for ³Hπ (48-200), small for ⁵Hσ (3-9), and intermediate for ³Hσ (23-51). The calculated KIEs for (TPFPP)Fe^IV═O without an axial ligand reveal that ³Hσ is the only pathway having a KIE that matches the experimental values, for the reactions with DHA and Xan (Angew. Chem., Int. Ed. 2008, 47, 7321). Indeed, theory shows that tunneling significantly lowers the ³Hσ barrier rendering it the sole reactive state for the reaction. A prediction is made for the reactivity and KIE of (TMP)Fe^IV═O complex, and a comparison is made with the analogous nonheme complexes.

Hydrogen Abstraction of Camphor Catalyzed by Cytochrome P450cam: A QM/MM Study

A combined quantum mechanics and molecular mechanics (QM/MM, QM = UB3LYP-D3, MM = AMBER) method is used to study the hydrogen abstraction reaction in P450_cam catalyzed hydroxylation of camphor in the quartet state. Compared to QM/MM calculations in the literature, this study uses larger basis sets for the most important atoms at the active site and QM/MM Hessian harmonic frequency calculations to determine the standard Gibbs free energy of activation and kinetic isotope effect. The QM/MM covalent boundary is treated with a capping hydrogen atom method, which is simple and robust. An energy barrier of 21.3 kcal/mol and a standard free energy of activation of 16.8 kcal/mol are obtained for this hydrogen abstraction reaction. These values are similar to those reported in the literature, suggesting that when a general protocol is followed, QM/MM results are reproducible. It is found that using a sufficiently large basis set is important to minimize basis set errors.

Insights into structure and redox potential of lignin peroxidase from QM/MM calculations

Redox potentials are computed for the active form (compound I) of lignin peroxidase (LiP) using a suitable QM/MM methodology (B3LYP/SDD/6-311G**//BP86/SVP:CHARMM). Allowing for dynamic conformational averaging, a potential of 0.67(33) V relative to ferrocenium/ferrocene is obtained for the active form with its oxoiron(iv) core. The computed redox potential is very sensitive to the charge distribution around the active site: protonation of titratable residues close to the metal center increases the redox potential, thereby rationalising the known pH dependence of LiP activity. A simple MM-charge deletion scheme is used to identify residues that are critical for the redox potential. Two mutant proteins are studied through homology modelling, E40Q and D183N, which are predicted to have an increased redox potential by 140 mV and 190 mV, respectively, relative to the wild type. These mutant proteins are thus promising targets for synthesis and further exploration toward a rational design of biocatalytic systems for oxidative degradation of lignin.

Theoretical Study of the Mechanism of Exemestane Hydroxylation Catalyzed by Human Aromatase Enzyme

Human aromatase (CYP19A1) aromatizes the androgens to form estrogens via a three-step oxidative process. The estrogens are necessary in humans, mainly in women, because of the role they play in sexual and reproductive development. However, these also are involved in the development and growth of hormone-dependent breast cancer. Therefore, inhibition of the enzyme aromatase, by means of drugs known as aromatase inhibitors, is the frontline therapy for these types of cancers. Exemestane is a suicidal third-generation inhibitor of aromatase, currently used in breast cancer treatment. In this study, the hydroxylation of exemestane catalyzed by aromatase has been studied by means of hybrid QM/MM methods. The Free Energy Perturbation calculations provided a free energy of activation for the hydrogen abstraction step (rate-limiting step) of 17 kcal/mol. The results reveal that the hydroxylation of exemestane is not the inhibition stage, suggesting a possible competitive mechanism between the inhibitor and the natural substrate androstenedione in the first catalytic subcycle of the enzyme. Furthermore, the analysis of the interaction energy for the substrate and the cofactor in the active site shows that the role of the enzymatic environment during this reaction consists of a transition state stabilization by means of electrostatic effects.

Quantum mechanical/molecular mechanical calculated reactivity networks reveal how cytochrome P450cam and Its T252A mutant select their oxidation pathways

Quantum mechanical/molecular mechanical calculations address the longstanding-question of a "second oxidant" in P450 enzymes wherein the proton-shuttle, which leads to formation of the "primary-oxidant" Compound I (Cpd I), was severed by mutating the crucial residue (in P450cam: Threonine-252-to-Alanine, hence T252A). Investigating the oxidant candidates Cpd I, ferric hydroperoxide, and ferric hydrogen peroxide (Fe(III)(O2H2)), and their reactions, generates reactivity networks which enable us to rule out a "second oxidant" and at the same time identify an additional coupling pathway that is responsible for the epoxidation of 5-methylenylcamphor by the T252A mutant. In this "second-coupling pathway", the reaction starts with the Fe(III)(O2H2) intermediate, which transforms to Cpd I via a O-O homolysis/H-abstraction mechanism. The persistence of Fe(III)(O2H2) and its oxidative reactivity are shown to be determined by interplay of substrate and protein. The substrate 5-methylenylcamphor prevents H2O2 release, while the protein controls the Fe(III)(O2H2) conversion to Cpd I by nailing-through hydrogen-bonding interactions-the conformation of the HO(•) radical produced during O-O homolysis. This conformation prevents HO(•) attack on the porphyrin's meso position, as in heme oxygenase, and prefers H-abstraction from Fe(IV)OH thereby generating H2O + Cpd I. Cpd I then performs substrate oxidations. Camphor cannot prevent H2O2 release and hence the T252A mutant does not oxidize camphor. This "second pathway" transpires also during H2O2 shunting of the cycle of wild-type P450cam, where the additional hydrogen-bonding with Thr252 prevents H2O2 release, and contributes to a successful Cpd I formation. The present results lead to a revised catalytic cycle of Cytochrome P450cam.

Theory uncovers an unusual mechanism of DNA repair of a lesioned adenine by AlkB enzymes

DNA-base lesions cause cancer and propagate into the genome. We use in-protein QM/MM calculations to study the repair of etheno-bridged adenine (εA) by the iron(IV)-oxo species of AlkB enzymes. Recent experimental investigations, using mass-spectrometry and in crystallo isolation, suggested that εA was repaired by formation of an epoxide (εA-ep) that further transforms to a glycol (εA-gl), ending finally in adenine and glyoxal. Theory reproduces the experimentally observed barrier for the rate-determining step and its pH dependence. However, as we show, the mass-spectrometrically identified species are side-byproducts unassociated with the repair mechanism. The repair is mediated by a zwitterionic species, of the same molecular mass as the epoxide, which transforms to an intermediate that matches the in crystallo trapped species in structure and mass, but is NOT the assumed εA-gl iron-glycol complex. Verifiable/falsifiable predictions, regarding the key protein residues, follow. The paper underscores the indispensable role of theory by providing atomistic descriptions of this vital mechanism, and guiding further experimental investigations.

Applications of density functional theory to iron-containing molecules of bioinorganic interest

The past decades have seen an explosive growth in the application of density functional theory (DFT) methods to molecular systems that are of interest in a variety of scientific fields. Owing to its balanced accuracy and efficiency, DFT plays particularly useful roles in the theoretical investigation of large molecules. Even for biological molecules such as proteins, DFT finds application in the form of, e.g., hybrid quantum mechanics and molecular mechanics (QM/MM), in which DFT may be used as a QM method to describe a higher prioritized region in the system, while a MM force field may be used to describe remaining atoms. Iron-containing molecules are particularly important targets of DFT calculations. From the viewpoint of chemistry, this is mainly because iron is abundant on earth, iron plays powerful (and often enigmatic) roles in enzyme catalysis, and iron thus has the great potential for biomimetic catalysis of chemically difficult transformations. In this paper, we present a brief overview of several recent applications of DFT to iron-containing non-heme synthetic complexes, heme-type cytochrome P450 enzymes, and non-heme iron enzymes, all of which are of particular interest in the field of bioinorganic chemistry. Emphasis will be placed on our own work.

EPR, resonance Raman, and DFT calculations on thiolate- and imidazole-bound iron(III) porphyrin complexes: role of the axial ligand in tuning the electronic structure

Iron(III) porphyrin complexes bearing covalently attached imidazole and thiolate axial ligands are investigated using resonance Raman, electron paramagnetic resonance, and cyclic voltammetry. The thiolate ligand stabilizes a low-spin ground state in solvent-bound six-coordinate species, weakens the Fe-N(pyr) bonds, and shifts the Fe(III/II) potential more negative by ~500 mV relative to an imidazole-bound species. Density functional theory calculations reproduce the experimental observation and indicate that the covalent charge donation from thiolate to iron reduces the Z(eff) on the iron. This increases the Fe(3d) orbital energies, which changes the bonding interaction present in these complexes significantly. In particular, the increase of the Fe(3d) energies activates an iron-to-porphyrin π*-back-bonding interaction not present in the imidazole-bound complex.

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.

Multireference and multiconfiguration ab initio methods in heme-related systems: what have we learned so far?

This work reviews the recent applications of ab initio multireference/multiconfiguration (MR/MC) electronic structure methods to heme-related systems, involving tetra-, penta-, and hexa-coordinate species, as well as the high-valent iron-oxo species. The current accuracy of these methods in the various systems is discussed, with special attention to potential sources of systematic errors. Thus, the review summarizes and tries to rationalize the key elements of MR/MC calculations, namely, the choice of the employed active space, especially the so-called double-shell effect that has already been recognized to be important in transition-metal-containing systems, and the impact of these elements on the spin-state energetics of heme species, as well as on the bonding mechanism of small molecules to the heme. It is shown that expansion of the MC wave function into one based on localized orbitals provides a compact and insightful view on some otherwise complex electronic structures. The effects of protein environment on the MR/MC results are summarized for the few available quantum mechanical/molecular mechanical (QM/MM) studies. Comparisons with corresponding DFT results are also made wherever available. Potential future directions are proposed.

Insights into the different dioxygen activation pathways of methane and toluene monooxygenase hydroxylases

The methane and toluene monooxygenase hydroxylases (MMOH and TMOH, respectively) have almost identical active sites, yet the physical and chemical properties of their oxygenated intermediates, designated P*, H(peroxo), Q, and Q* in MMOH and ToMOH(peroxo) in a subclass of TMOH, ToMOH, are substantially different. We review and compare the structural differences in the vicinity of the active sites of these enzymes and discuss which changes could give rise to the different behavior of H(peroxo) and Q. In particular, analysis of multiple crystal structures reveals that T213 in MMOH and the analogous T201 in TMOH, located in the immediate vicinity of the active site, have different rotatory configurations. We study the rotational energy profiles of these threonine residues with the use of molecular mechanics (MM) and quantum mechanics/molecular mechanics (QM/MM) computational methods and put forward a hypothesis according to which T213 and T201 play an important role in the formation of different types of peroxodiiron(III) species in MMOH and ToMOH. The hypothesis is indirectly supported by the QM/MM calculations of the peroxodiiron(III) models of ToMOH and the theoretically computed Mössbauer spectra. It also helps explain the formation of two distinct peroxodiiron(III) species in the T201S mutant of ToMOH. Additionally, a role for the ToMOD regulatory protein, which is essential for intermediate formation and protein functioning in the ToMO system, is advanced. We find that the low quadrupole splitting parameter in the Mössbauer spectrum observed for a ToMOH(peroxo) intermediate can be explained by protonation of the peroxo moiety, possibly stabilized by the T201 residue. Finally, similarities between the oxygen activation mechanisms of the monooxygenases and cytochrome P450 are discussed.

Elucidating the role of the proximal cysteine hydrogen-bonding network in ferric cytochrome P450cam and corresponding mutants using magnetic circular dichroism spectroscopy

Although extensive research has been performed on various cytochrome P450s, especially Cyt P450cam, there is much to be learned about the mechanism of how its functional unit, a heme b ligated by an axial cysteine, is finely tuned for catalysis by its second coordination sphere. Here we study how the hydrogen-bonding network affects the proximal cysteine and the Fe-S(Cys) bond in ferric Cyt P450cam. This is accomplished using low-temperature magnetic circular dichroism (MCD) spectroscopy on wild-type (wt) Cyt P450cam and on the mutants Q360P (pure ferric high-spin at low temperature) and L358P where the "Cys pocket" has been altered (by removing amino acids involved in the hydrogen-bonding network), and Y96W (pure ferric low-spin). The MCD spectrum of Q360P reveals fourteen electronic transitions between 15200 and 31050 cm(-1). Variable-temperature variable-field (VTVH) saturation curves were used to determine the polarizations of these electronic transitions with respect to in-plane (xy) and out-of-plane (z) polarization relative to the heme. The polarizations, oscillator strengths, and TD-DFT calculations were then used to assign the observed electronic transitions. In the lower energy region, prominent bands at 15909 and 16919 cm(-1) correspond to porphyrin (P) → Fe charge transfer (CT) transitions. The band at 17881 cm(-1) has distinct sulfur S(π) → Fe CT contributions. The Q band is observed as a pseudo A-term (derivative shape) at 18604 and 19539 cm(-1). In the case of the Soret band, the negative component of the expected pseudo A-term is split into two features due to mixing with another π → π* and potentially a P → Fe CT excited state. The resulting three features are observed at 23731, 24859, and 25618 cm(-1). Most importantly, the broad, prominent band at 28570 cm(-1) is assigned to the S(σ) → Fe CT transition, whose intensity is generated through a multitude of CT transitions with strong iron character. For wt, Q360P, and L358P, this band occurs at 28724, 28570, and 28620 cm(-1), respectively. The small shift of this feature upon altering the hydrogen bonds to the proximal cysteine indicates that the role of the Cys pocket is not primarily for electronic fine-tuning of the sulfur donor strength but is more for stabilizing the proximal thiolate against external reactants (NO, O(2), H(3)O(+)), and for properly positioning cysteine to coordinate to the iron center. This aspect is discussed in detail.

Color tuning in short wavelength-sensitive human and mouse visual pigments: ab initio quantum mechanics/molecular mechanics studies

We have investigated the protonation state and photoabsorption spectrum of Schiff-base (SB) nitrogen bound 11-cis-retinal in human blue and mouse UV cone visual pigments as well as in bovine rhodopsin by hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. We have employed both multireference (MRCISD+Q, MR-SORCI+Q, and MR-DDCI2+Q) and single reference (TD-B3LYP and RI-CC2) QM methods. The calculated ground-state and vertical excitation energies show that UV-sensitive pigments have deprotonated SB nitrogen, while violet-sensitive pigments have protonated SB nitrogen, in agreement with some indirect experimental evidence. A significant blue shift of the absorption maxima of violet-sensitive pigments relative to rhodopsins arises from the increase in bond length alternation of the polyene chain of 11-cis-retinal induced by polarizing fields of these pigments. The main counterion is Glu113 in both violet-sensitive vertebrate pigments and bovine rhodopsin. Neither Glu113 nor the remaining pigment has a significant influence on the first excitation energy of 11-cis-retinal in the UV-sensitive pigments that have deprotonated SB nitrogen. There is no charge transfer between the SB and beta-ionone terminals of 11-cis-retinal in the ground and first excited states.

Finite-temperature effects in enzymatic reactions — Insights from QM/MM free-energy simulations

We report potential-energy and free-energy data for three enzymatic reactions: carbon–halogen bond formation in fluorinase, hydrogen abstraction from camphor in cytochrome P450cam, and chorismate-to-prephenate Claisen rearrangement in chorismate mutase. The results were obtained by combined quantum mechanics/molecular mechanics (QM/MM) optimizations and two types of QM/MM free-energy simulations (free-energy perturbation and umbrella sampling) using semi-empirical or density-functional QM methods. Based on these results and our previously published free-energy data on electrophilic substitution in para-hydroxybenzoate hydroxylase, we discuss the importance of finite-temperature effects in the chemical step of enzyme reactions. We find that the entropic contribution to the activation barrier is generally rather small, usually of the order of 5 kJ mol–1 or less, consistent with the notion that enzymes bind and pre-organize the reactants in the active site. A somewhat larger entropic contribution is encountered in the case of chorismate mutase where the pericyclic transition state is intrinsically more rigid than the chorismate reactant (also in the enzyme). The present results suggest that barriers from QM/MM geometry optimization may often be close to free-energy barriers for the chemical step in enzymatic reactions.

QM/MM methods for biomolecular systems

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.