Effect of Geometry Optimizations on QM-Cluster and QM/MM Studies of Reaction Energies in Proteins

Multiscale Modelling of Lytic Polysaccharide Monooxygenases

Lytic polysaccharide monooxygenase (LPMO) enzymes have attracted considerable attention owing to their ability to enhance polysaccharide depolymerization, making them interesting with respect to production of biofuel from cellulose. LPMOs are metalloenzymes that contain a mononuclear copper active site, capable of activating dioxygen. However, many details of this activation are unclear. Some aspects of the mechanism have previously been investigated from a computational angle. Yet, either these studies have employed only molecular mechanics (MM), which are inaccurate for metal active sites, or they have described only the active site with quantum mechanics (QM) and neglected the effect of the protein. Here, we employ hybrid QM and MM (QM/MM) methods to investigate the first steps of the LPMO mechanism, which is reduction of Cu^II to Cu^I and the formation of a Cu^II-superoxide complex. In the latter complex, the superoxide can bind either in an equatorial or an axial position. For both steps, we obtain structures that are markedly different from previous suggestions, based on small QM-cluster calculations. Our calculations show that the equatorial isomer of the superoxide complex is over 60 kJ/mol more stable than the axial isomer because it is stabilized by interactions with a second-coordination-sphere glutamine residue, suggesting a possible role for this residue. The coordination of superoxide in this manner agrees with recent experimental suggestions.

Influence of Coupling and Embedding Schemes on QM Size Convergence in QM/MM Approaches for the Example of a Proton Transfer in DNA

The influence of embedding and coupling schemes on the convergence of the QM size in the QM/MM approach is investigated for the transfer of a proton in a DNA base pair. We find that the embedding scheme (mechanical or electrostatic) has a much greater impact on the convergence behavior than the coupling scheme (additive QM/MM or subtractive ONIOM). To achieve size convergence, QM regions with up to 6000 atoms are necessary for pure QM or mechanical embedding. In contrast, electrostatic embedding converges faster: for the example of the transfer of a proton between DNA base pairs, we recommend including at least five base pairs and 5 Å of solvent (including counterions) into the QM region, i.e., a total of 1150 atoms.

H2binding to the active site of [NiFe] hydrogenase studied by multiconfigurational and coupled-cluster methods

CCSD(T) and DMRG-CASPT2 calculations show that H₂prefers to bind to Ni rather than to Fe in [NiFe] hydrogenase.

Computational tools for the evaluation of laboratory-engineered biocatalysts

Biocatalysis is based on the application of natural catalysts for new purposes, for which enzymes were not designed. Although the first examples of biocatalysis were reported more than a century ago, biocatalysis was revolutionized after the discovery of an in vitro version of Darwinian evolution called Directed Evolution (DE). Despite the recent advances in the field, major challenges remain to be addressed. Currently, the best experimental approach consists of creating multiple mutations simultaneously while limiting the choices using statistical methods. Still, tens of thousands of variants need to be tested experimentally, and little information is available on how these mutations lead to enhanced enzyme proficiency. This review aims to provide a brief description of the available computational techniques to unveil the molecular basis of improved catalysis achieved by DE. An overview of the strengths and weaknesses of current computational strategies is explored with some recent representative examples. The understanding of how this powerful technique is able to obtain highly active variants is important for the future development of more robust computational methods to predict amino-acid changes needed for activity.

How Large Should the QM Region Be in QM/MM Calculations? The Case of Catechol O-Methyltransferase

Hybrid quantum mechanical-molecular mechanical (QM/MM) simulations are widely used in studies of enzymatic catalysis. Until recently, it has been cost prohibitive to determine the asymptotic limit of key energetic and structural properties with respect to increasingly large QM regions. Leveraging recent advances in electronic structure efficiency and accuracy, we investigate catalytic properties in catechol O-methyltransferase, a prototypical methyltransferase critical to human health. Using QM regions ranging in size from reactants-only (64 atoms) to nearly one-third of the entire protein (940 atoms), we show that properties such as the activation energy approach within chemical accuracy of the large-QM asymptotic limits rather slowly, requiring approximately 500-600 atoms if the QM residues are chosen simply by distance from the substrate. This slow approach to asymptotic limit is due to charge transfer from protein residues to the reacting substrates. Our large QM/MM calculations enable identification of charge separation for fragments in the transition state as a key component of enzymatic methyl transfer rate enhancement. We introduce charge shift analysis that reveals the minimum number of protein residues (approximately 11-16 residues or 200-300 atoms for COMT) needed for quantitative agreement with large-QM simulations. The identified residues are not those that would be typically selected using criteria such as chemical intuition or proximity. These results provide a recipe for a more careful determination of QM region sizes in future QM/MM studies of enzymes.

Pathways for Arene Oxidation in Non-Heme Diiron Enzymes: Lessons from Computational Studies on Benzoyl Coenzyme A Epoxidase

Oxygenation of aromatic rings using O₂ is catalyzed by several non-heme carboxylate-bridged diiron enzymes. In order to provide a general mechanistic description for these reactions, computational studies were carried out at the ONIOM(B3LYP/BP86/Amber) level on the non-heme diiron enzyme benzoyl coenzyme A epoxidase, BoxB. The calculations revealed four possible pathways for attacking the aromatic ring: (a) electrophilic (2e^-) attack by a bis(μ-oxo)-diiron(IV) species (Q pathway); (b) electrophilic (2e^-) attack via the σ* orbital of a μ-η²:η²-peroxo-diiron(III) intermediate (Pσ* pathway); (c) radical (1e^-) attack via the π*-orbital of a superoxo-diiron(II,III) species (Pπ* pathway); (d) radical (1e^-) attack of a partially quenched bis(μ-oxo)-diiron(IV) intermediate (Q' pathway). The results allowed earlier work of de Visser on olefin epoxidation by diiron complexes and QM-cluster studies of Liao and Siegbahn on BoxB to be put into a broader perspective. Parallels with epoxidation using organic peracids were also examined. Specifically for the BoxB enzyme, the Q pathway was found to be the most preferred, but the corresponding bis(μ-oxo)-diiron(IV) species is significantly destabilized and not expected to be directly observable. Epoxidation via the Pσ* pathway represents an energetically somewhat higher lying alternative; possible strategies for experimental discrimination are discussed. The selectivity toward epoxidation is shown to stem from a combination of inherent electronic properties of the thioacyl substituent and enzymatic constraints. Possible implications of the results for toluene monooxygenases are considered as well.

How Many Conformations of Enzymes Should Be Sampled for DFT/MM Calculations? A Case Study of Fluoroacetate Dehalogenase

The quantum mechanics/molecular mechanics (QM/MM) method (e.g., density functional theory (DFT)/MM) is important in elucidating enzymatic mechanisms. It is indispensable to study "multiple" conformations of enzymes to get unbiased energetic and structural results. One challenging problem, however, is to determine the minimum number of conformations for DFT/MM calculations. Here, we propose two convergence criteria, namely the Boltzmann-weighted average barrier and the disproportionate effect, to tentatively address this issue. The criteria were tested by defluorination reaction catalyzed by fluoroacetate dehalogenase. The results suggest that at least 20 conformations of enzymatic residues are required for convergence using DFT/MM calculations. We also tested the correlation of energy barriers between small QM regions and big QM regions. A roughly positive correlation was found. This kind of correlation has not been reported in the literature. The correlation inspires us to propose a protocol for more efficient sampling. This saves 50% of the computational cost in our current case.

Noncovalent interactions between the second coordination sphere and the active site of [NiFeSe] hydrogenase

QM/MM studies on seven truncated models of the oxidized as-isolated state of the [NiFeSe] Hases reveal the influence of the residues in the second coordination sphere on the active site.

Requirement for transient metal ions revealed through computational analysis for DNA polymerase going in reverse

DNA polymerases facilitate faithful insertion of nucleotides, a central reaction occurring during DNA replication and repair. DNA synthesis (forward reaction) is "balanced," as dictated by the chemical equilibrium by the reverse reaction of pyrophosphorolysis. Two closely spaced divalent metal ions (catalytic and nucleotide-binding metals) provide the scaffold for these reactions. The catalytic metal lowers the pKa of O3' of the growing primer terminus, and the nucleotide-binding metal facilitates substrate binding. Recent time-lapse crystallographic studies of DNA polymerases have identified an additional metal ion (product metal) associated with pyrophosphate formation, leading to the suggestion of its possible involvement in the reverse reaction. Here, we establish a rationale for a role of the product metal using quantum mechanical/molecular mechanical calculations of the reverse reaction in the confines of the DNA polymerase β active site. Additionally, site-directed mutagenesis identifies essential residues and metal-binding sites necessary for pyrophosphorolysis. The results indicate that the catalytic metal site must be occupied by a magnesium ion for pyrophosphorolysis to occur. Critically, the product metal site is occupied by a magnesium ion early in the pyrophosphorolysis reaction path but must be removed later. The proposed dynamic nature of the active site metal ions is consistent with crystallographic structures. The transition barrier for pyrophosphorolysis was estimated to be significantly higher than that for the forward reaction, consistent with kinetic activity measurements of the respective reactions. These observations provide a framework to understand how ions and active site changes could modulate the internal chemical equilibrium of a reaction that is central to genome stability.

The Importance of the MM Environment and the Selection of the QM Method in QM/MM Calculations: Applications to Enzymatic Reactions

In this chapter, we discuss the influence of an anisotropic protein environment on the reaction mechanisms of saccharopine reductase and uroporphyrinogen decarboxylase, respectively, via the use of a quantum mechanical and molecular mechanical (QM/MM) approach. In addition, we discuss the importance of selecting a suitable DFT functional to be used in a QM/MM study of a key intermediate in the mechanism of 8R-lipoxygenase, a nonheme iron enzyme. In the case of saccharopine reductase, while the enzyme utilizes a substrate-assisted catalytic pathway, it was found that only through treating the polarizing effect of the active site, via the use of an electronic embedding formalism, was agreement with experimental kinetic data obtained. Similarly, in the case of uroporphyrinogen decarboxylase, the effect of the protein environment on the catalytic mechanism was found to be such that the calculated rate-limiting barrier is in good agreement with related experimentally determined values for the first decarboxylation of the substrate. For 8R-lipoxygenase, it was found that the geometries and energies of the multicentered open-shell intermediate complexes formed during the mechanism are quite sensitive to the choice of the density functional theory method. Thus, while density functional theory has become the method of choice in QM/MM studies, care must be taken in the selection of a particular high-level method.

Converging nuclear magnetic shielding calculations with respect to basis and system size in protein systems

Ab initio chemical shielding calculations greatly facilitate the interpretation of nuclear magnetic resonance (NMR) chemical shifts in biological systems, but the large sizes of these systems requires approximations in the chemical models used to represent them. Achieving good convergence in the predicted chemical shieldings is necessary before one can unravel how other complex structural and dynamical factors affect the NMR measurements. Here, we investigate how to balance trade-offs between using a better basis set or a larger cluster model for predicting the chemical shieldings of the substrates in two representative examples of protein-substrate systems involving different domains in tryptophan synthase: the N-(4'-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F9) ligand which binds in the α active site, and the 2-aminophenol quinonoid intermediate formed in the β active site. We first demonstrate that a chemically intuitive three-layer, locally dense basis model that uses a large basis on the substrate, a medium triple-zeta basis to describe its hydrogen-bonding partners and/or surrounding van der Waals cavity, and a crude basis set for more distant atoms provides chemical shieldings in good agreement with much more expensive large basis calculations. Second, long-range quantum mechanical interactions are important, and one can accurately estimate them as a small-basis correction to larger-basis calculations on a smaller cluster. The combination of these approaches enables one to perform density functional theory NMR chemical shift calculations in protein systems that are well-converged with respect to both basis set and cluster size.

A base-independent repair mechanism for DNA glycosylase--no discrimination within the active site

The ubiquitous occurrence of DNA damages renders its repair machinery a crucial requirement for the genomic stability and the survival of living organisms. Deficiencies in DNA repair can lead to carcinogenesis, Alzheimer, or Diabetes II, where increased amounts of oxidized DNA bases have been found in patients. Despite the highest mutation frequency among oxidized DNA bases, the base-excision repair process of oxidized and ring-opened guanine, FapydG (2,6-diamino-4-hydroxy-5-formamidopyrimidine), remained unclear since it is difficult to study experimentally. We use newly-developed linear-scaling quantum-chemical methods (QM) allowing us to include up to 700 QM-atoms and achieving size convergence. Instead of the widely assumed base-protonated pathway we find a ribose-protonated repair mechanism which explains experimental observations and shows strong evidence for a base-independent repair process. Our results also imply that discrimination must occur during recognition, prior to the binding within the active site.

Use of Broken-Symmetry Density Functional Theory To Characterize the IspH Oxidized State: Implications for IspH Mechanism and Inhibition

With current therapies becoming less efficacious due to increased drug resistance, new inhibitors of both bacterial and malarial targets are desperately needed. The recently discovered methylerythritol phosphate (MEP) pathway for isoprenoid synthesis provides novel targets for the development of such drugs. Particular attention has focused on the IspH protein, the final enzyme in the MEP pathway, which uses its [4Fe-4S] cluster to catalyze the formation of the isoprenoid precursors IPP and DMAPP from HMBPP. IspH catalysis is achieved via a 2e-/2H+ reductive dehydroxylation of HMBPP; the mechanism by which catalysis is achieved, however, is highly controversial. The work presented herein provides the first step in assessing different routes to catalysis by using computational methods. By performing broken-symmetry density functional theory (BS-DFT) calculations that employ both the conductor-like screening solvation model (DFT/COSMO) and a finite-difference Poisson-Boltzmann self-consistent reaction field methodology (DFT/SCRF), we evaluate geometries, energies, and Mössbauer signatures of the different protonation states that may exist in the oxidized state of the IspH catalytic cycle. From DFT/SCRF computations performed on the oxidized state, we find a state where the substrate, HMBPP, coordinates the apical iron in the [4Fe-4S] cluster as an alcohol group (ROH) to be one of two, isoenergetic, lowest-energy states. In this state, the HMBPP pyrophosphate moiety and an adjacent glutamate residue (E126) are both fully deprotonated, making the active site highly anionic. Our findings that this low-energy state also matches the experimental geometry of the active site and that its computed isomer shifts agree with experiment validate the use of the DFT/SCRF method to assess relative energies along the IspH reaction pathway. Additional studies of IspH catalytic intermediates are currently being pursued.

The structural and energetic aspects of substrate binding and the mechanism of action of the DapE-encoded N-succinyl-l,l-diaminopimelic acid desuccinylase (DapE) investigated using a hybrid QM/MM method

Substrate binding and the mechanism of action of the DapE-encodedN-succinyl-l,l-diaminopimelic acid desuccinylase (DapE).

Accurate calculations of geometries and singlet–triplet energy differences for active-site models of [NiFe] hydrogenase

The relative stability of singlet and triplet state models of [NiFe] hydrogenase have been studied by advanced quantum-mechanical methods.