A pseudobond parametrization for improved electrostatics in quantum mechanical/molecular mechanical simulations of enzymes

A Polarizable QM/MM Model That Combines the State-Averaged CASSCF and AMOEBA Force Field for Photoreactions in Proteins

This study presents the polarizable quantum mechanics/molecular mechanics (QM/MM) embedding of the state-averaged complete active space self-consistent field (SA-CASSCF) in the atomic multipole optimized energetics for biomolecular applications (AMOEBA) force field for the purpose of studying photoreactions in protein environments. We describe two extensions of our previous work that combine SA-CASSCF with AMOEBA water models, allowing it to be generalized to AMOEBA models for proteins and other macromolecules. First, we discuss how our QM/MM model accounts for the discrepancy between the direct and polarization electric fields that arises in the AMOEBA description of intramolecular polarization. A second improvement is the incorporation of link atom schemes to treat instances in which the QM/MM boundary goes through covalent bonds. A single-link atom scheme and double-link atom scheme are considered in this work, and we will discuss how electrostatic interaction, van der Waals interaction, and various kinds of valence terms are treated across the boundary. To test the accuracy of the link atom scheme, we will compare QM/MM with full QM calculations and study how the errors in ground state properties, excited state properties, and excitation energies change when tuning the parameters in the link atom scheme. We will also test the new SA-CASSCF/AMOEBA method on an elementary reaction step in NanoLuc, an artificial bioluminescence luciferase. We will show how the reaction mechanism is different when calculated in the gas phase, in polarizable continuum medium (PCM), versus in protein AMOEBA models.

PyDFT-QMMM: A modular, extensible software framework for DFT-based QM/MM molecular dynamics

PyDFT-QMMM is a Python-based package for performing hybrid quantum mechanics/molecular mechanics (QM/MM) simulations at the density functional level of theory. The program is designed to treat short-range and long-range interactions through user-specified combinations of electrostatic and mechanical embedding procedures within periodic simulation domains, providing necessary interfaces to external quantum chemistry and molecular dynamics software. To enable direct embedding of long-range electrostatics in periodic systems, we have derived and implemented force terms for our previously described QM/MM/PME approach [Pederson and McDaniel, J. Chem. Phys. 156, 174105 (2022)]. Communication with external software packages Psi4 and OpenMM is facilitated through Python application programming interfaces (APIs). The core library contains basic utilities for running QM/MM molecular dynamics simulations, and plug-in entry-points are provided for users to implement custom energy/force calculation and integration routines, within an extensible architecture. The user interacts with PyDFT-QMMM primarily through its Python API, allowing for complex workflow development with Python scripting, for example, interfacing with PLUMED for free energy simulations. We provide benchmarks of forces and energy conservation for the QM/MM/PME and alternative QM/MM electrostatic embedding approaches. We further demonstrate a simple example use case for water solute in a water solvent system, for which radial distribution functions are computed from 100 ps QM/MM simulations; in this example, we highlight how the solvation structure is sensitive to different basis-set choices due to under- or over-polarization of the QM water molecule's electron density.

Effects of the Y432S Cancer-Associated Variant on the Reaction Mechanism of Human DNA Polymerase κ

Human DNA polymerases are vital for genetic information management. Their function involves catalyzing the synthesis of DNA strands with unparalleled accuracy, which ensures the fidelity and stability of the human genomic blueprint. Several disease-associated mutations and their functional impact on DNA polymerases have been reported. One particular polymerase, human DNA polymerase kappa (Pol κ), has been reported to be susceptible to several cancer-associated mutations. The Y432S mutation in Pol κ, associated with various cancers, is of interest due to its impact on polymerization activity and markedly reduced thermal stability. Here, we have used computational simulations to investigate the functional consequences of the Y432S using classical molecular dynamics (MD) and coupled quantum mechanics/molecular mechanics (QM/MM) methods. Our findings suggest that Y432S induces structural alterations in domains responsible for nucleotide addition and ternary complex stabilization while retaining structural features consistent with possible catalysis in the active site. Calculations of the minimum energy path associated with the reaction mechanism of the wild type (WT) and Y432S Pol κ indicate that, while both enzymes are catalytically competent (in terms of energetics and the active site's geometries), the cancer mutation results in an endoergic reaction and an increase in the catalytic barrier. Interactions with a third magnesium ion and environmental effects on nonbonded interactions, particularly involving key residues, contribute to the kinetic and thermodynamic distinctions between the WT and mutant during the catalytic reaction. The energetics and electronic findings suggest that active site residues favor the catalytic reaction with dCTP3- over dCTP4-.

Impact of an Ionic Liquid Solution on Horseradish Peroxidase Activity

Horseradish peroxidase (HRP) is an enzyme that oxidizes pollutants from wastewater. A previous report indicated that peroxidases can have an enhancement in initial enzymatic activity in an aqueous solution of 0.26 M 1-ethyl-3-methylimidazolium ethyl sulfate ([EMIm][EtSO4]) at neutral pH. However, the atomistic details remain elusive. In the enzymatic landscape of HRP, compound II (Cpd II) plays a key role and involves a histidine (H42) residue. Cpd II exists as oxoferryl (2a) or hydroxoferryl (2b(FeIV)) forms, where 2a is the predominantly observed form in experimental studies. Intriguingly, the ferric 2b(FeIII) form seen in synthetic complexes has not been observed in HRP. Here, we have investigated the structure and dynamics of HRP in pure water and aqueous [EMIm][EtSO4] (0.26 M), as well as the reaction mechanism of 2a to 2b conversion using polarizable molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations. When HRP is solvated in aq [EMIm][EtSO4], the catalytic water displaces, and H42 directly orients over the ferryl moiety, allowing a direct proton transfer (PT) with a significant energy barrier reduction. Conversely, in neat water, the reaction of 2a to 2b follows the previously reported mechanism. We further investigated the deprotonated form of H42. Analysis of the electric fields at the active site indicates that the aq [EMIm][EtSO4] medium facilitates the reaction by providing a more favorable environment compared with the system solvated in neat water. Overall, the atomic level supports the previous experimental observations and underscores the importance of favorable electric fields in the active site to promote catalysis.

Distal Mutations in the β-Clamp of DNA Polymerase III* Disrupt DNA Orientation and Affect Exonuclease Activity

DNA polymerases are responsible for the replication and repair of DNA found in all DNA-based organisms. DNA Polymerase III is the main replicative polymerase of E. coli and is composed of over 10 proteins. A subset of these proteins (Pol III*) includes the polymerase (α), exonuclease (ϵ), clamp (β), and accessory protein (θ). Mutations of residues in, or around the active site of the catalytic subunits (α and ϵ), can have a significant impact on catalysis. However, the effects of distal mutations in noncatalytic subunits on the activity of catalytic subunits are less well-characterized. Here, we investigate the effects of two Pol III* variants, β-L82E/L82'E and β-L82D/L82'D, on the proofreading reaction catalyzed by ϵ. MD simulations reveal major changes in the dynamics of Pol III*, which extend throughout the complex. These changes are mostly induced by a shift in the position of the DNA substrate inside the β-clamp, although no major structural changes are observed in the protein complex. Quantum mechanics/molecular mechanics (QM/MM) calculations indicate that the β-L82D/L82'D variant has reduced catalytic proficiency due to highly endoergic reaction energies resulting from structural changes in the active site and differences in the electric field at the active site arising from the protein and substrate. Conversely, the β-L82E/L82'E variant is predicted to maintain proofreading activity, exhibiting a similar reaction barrier for nucleotide excision compared with the WT system. However, significant differences in the reaction mechanism are obtained due to the changes induced by the mutations on the β-clamp.

Computational investigations of selected enzymes from two iron and α-ketoglutarate-dependent families

DNA alkylation is used as the key epigenetic mark in eukaryotes, however, most alkylation in DNA can result in deleterious effects. Therefore, this process needs to be tightly regulated. The enzymes of the AlkB and Ten-Eleven Translocation (TET) families are members of the Fe and alpha-ketoglutarate-dependent superfamily of enzymes that are tasked with dealkylating DNA and RNA in cells. Members of these families span all species and are an integral part of transcriptional regulation. While both families catalyze oxidative dealkylation of various bases, each has specific preference for alkylated base type as well as distinct catalytic mechanisms. This perspective aims to provide an overview of computational work carried out to investigate several members of these enzyme families including AlkB, ALKB Homolog 2, ALKB Homolog 3 and Ten-Eleven Translocate 2. Insights into structural details, mutagenesis studies, reaction path analysis, electronic structure features in the active site, and substrate preferences are presented and discussed.

Development and application of quantum mechanics/molecular mechanics methods with advanced polarizable potentials

Quantum mechanics/molecular mechanics (QM/MM) simulations are a popular approach to study various features of large systems. A common application of QM/MM calculations is in the investigation of reaction mechanisms in condensed-phase and biological systems. The combination of QM and MM methods to represent a system gives rise to several challenges that need to be addressed. The increase in computational speed has allowed the expanded use of more complicated and accurate methods for both QM and MM simulations. Here, we review some approaches that address several common challenges encountered in QM/MM simulations with advanced polarizable potentials, from methods to account for boundary across covalent bonds and long-range effects, to polarization and advanced embedding potentials.

Avoiding Electron Spill-Out in QM/MM Calculations on Excited States with Simple Pseudopotentials

QM/MM calculations of electronic excitations with diffuse basis sets often have large errors due to spill-out of electrons from the quantum subsystem. The Pauli repulsion of the electrons by the environment has to be included to avoid this. We propose transferable atomic all-electron pseudopotentials that can readily be combined with most MM force fields to avoid electron spill-out. QM/MM excitation energies computed with time-dependent Hartree-Fock and the algebraic diagrammatic construction through second-order are benchmarked against supermolecular calculations to validate these new pseudopotentials. The QM/MM calculations with pseudopotentials give accurate results that are stable with augmentation of the basis set with diffuse functions. We show that the largest contribution to residual deviations from full QM calculations is caused by the missing London dispersion interaction.

Ground State Destabilization in Uracil DNA Glycosylase: Let's Not Forget "Tautomeric Strain" in Substrates

Enzymes like uracil DNA glycosylase (UDG) can achieve ground state destabilization, by polarizing substrates to mimic rare tautomers. On the basis of computed nucleus independent chemical shifts, NICS(1)_zz, and harmonic oscillator model of electron delocalization (HOMED) analyses, of quantum mechanics (QM) and quantum mechanics/molecular mechanics (QM/MM) models of the UDG active site, uracil is strongly polarized when bound to UDG and resembles a tautomer >12 kcal/mol higher in energy. Natural resonance theory (NRT) analyses identified a dominant O2 imidate resonance form for residue bound 1-methyl-uracil. This "tautomeric strain" raises the energy of uracil, making uracilate a better than expected leaving group. Computed gas-phase S_N2 reactions of free and hydrogen bonded 1-methyl-uracil demonstrate the relationship between the degree of polarization in uracil and the leaving group ability of uracilate.

Towards large scale hybrid QM/MM dynamics of complex systems with advanced point dipole polarizable embeddings

In this work, we present a general route to hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) Molecular Dynamics for complex systems using a polarizable embedding. We extend the capabilities of our hybrid framework, combining the Gaussian and Tinker/Tinker-HP packages in the context of the AMOEBA polarizable force field to treat large (bio)systems where the QM and the MM subsystems are covalently bound, adopting pseudopotentials at the boundaries between the two regions. We discuss in detail the implementation and demonstrate the global energy conservation of our QM/MM Born-Oppenheimer molecular dynamics approach using Density Functional Theory. Finally, the approach is assessed on the electronic absorption properties of a 16 500 atom complex encompassing an organic dye embedded in a DNA matrix in solution, extending the hybrid method to a time-dependent Density Functional Theory approach. The results obtained comparing different partitions between the quantum and the classical subsystems also suggest that large QM portions are not necessary if accurate polarizable force fields are used in a variational formulation of the embedding, properly including the QM/MM mutual polarization.

LICHEM 1.1: Recent Improvements and New Capabilities

The QM/MM method has become a useful tool to investigate various properties of complex systems. We previously introduced the layered interacting chemical models (LICHEM) package to enable QM/MM simulations with advanced potentials by combining various (unmodified) QM and MM codes ( J. Comp. Chem., 2016, 37, 1019). LICHEM provides several capabilities such as the ability to use polarizable force fields, such as AMOEBA, for the MM environment. Here, we describe an updated version of LICHEM (v1.1), which includes several new functionalities including a new method to account for long-range electrostatic effects in QM/MM (QM/MM-LREC), a new implementation for QM/MM with the Gaussian electrostatic model (GEM), and new capabilities for path optimizations using the quadratic string model (QSM) coupled with restrained MM environment optimization.

Insight into wild-type and T1372E TET2-mediated 5hmC oxidation using ab initio QM/MM calculations

Ten-eleven translocation 2 (TET2) is an Fe/α-ketoglutarate (α-KG) dependent enzyme that dealkylates 5-methylcytosine (5mC). The reaction mechanism involves a series of three sequential oxidations that convert 5mC to 5-hydroxy-methylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC). Our previous biochemical and computational studies uncovered an active site scaffold that is required for wild-type (WT) stepwise oxidation (Nat. Chem. Bio., 13, 181). We showed that the mutation of a single residue, T1372 to some amino acids, such as Glu, can impact the iterative oxidation steps and stop the oxidation of 5hmC to 5fC/caC. However, the source of the stalling at the first oxidation step by some mutant TET proteins still remains unclear. Here, we studied the catalytic mechanism of oxidation of 5hmC to 5fC by WT and T1372E TET2 using an ab initio quantum mechanical/molecular mechanical (QM/MM) approach. Our results suggest that the rate limiting step for WT TET2 involves a hydrogen atom abstraction from the hydroxyl group of 5hmC by the ferryl moiety in the WT. By contrast, our calculations for the T1372E mutant indicate that the rate limiting step for this variant corresponds to a second proton abstraction and the calculated barrier is almost twice as large as for WT TET2. Our results suggest that the large barrier for the 5hmC to 5fC oxidation in this mutant is due (at least in part) to the unfavorable orientation of the substrate in the active site. Combined electron localization function (ELF) and non-covalent interaction (NCI) analyses provide a qualitative description of the evolution of the electronic structure of the active site along the reaction path. Energy decomposition analysis (EDA) has been performed on the WT to investigate the impact of each MM residue on catalytic activity.

Accurate Quantum Mechanical/Molecular Mechanical Calculations of Reduction Potentials in Azurin Variants

Understanding the regulation mechanism and molecular determinants of the reduction potential of metalloprotein is a major challenge. An ab initio quantum mechanical/molecular mechanical (QM/MM) method combining the minimum free energy path (MFEP) and fractional number of electron (FNE) approaches has been developed in our group to simulate the redox processes of large systems. The FNE scheme provides an efficient unique description for the redox process, while the MFEP method provides improved conformational sampling on complex environments such as protein in the QM/MM calculations. The reduction potentials of wild-type and seven mutants of azurin, a type 1 copper metalloprotein, were simulated with the QM/MM-MFEP+FNE approach in this paper. A range of 350 mV for the variations of the reduction potentials of these azurin proteins was reproduced faithfully with relative errors around 20 mV. The correlation between structural interactions and reduction potentials observed in simulations provides in-depth insight into the regulation of reduction potentials, which potentially can also be very useful to the engineering of metalloprotein-based electrocatalysts in artificial photosynthesis. The excellent accuracy and efficiency of the QM/MM-MFEP+FNE approach demonstrate the potential for simulations of many electron transfer processes in condensed phases and biochemical systems.

Tinker-HP: a massively parallel molecular dynamics package for multiscale simulations of large complex systems with advanced point dipole polarizable force fields

We present Tinker-HP, a massively MPI parallel package dedicated to classical molecular dynamics (MD) and to multiscale simulations, using advanced polarizable force fields (PFF) encompassing distributed multipoles electrostatics. Tinker-HP is an evolution of the popular Tinker package code that conserves its simplicity of use and its reference double precision implementation for CPUs. Grounded on interdisciplinary efforts with applied mathematics, Tinker-HP allows for long polarizable MD simulations on large systems up to millions of atoms. We detail in the paper the newly developed extension of massively parallel 3D spatial decomposition to point dipole polarizable models as well as their coupling to efficient Krylov iterative and non-iterative polarization solvers. The design of the code allows the use of various computer systems ranging from laboratory workstations to modern petascale supercomputers with thousands of cores. Tinker-HP proposes therefore the first high-performance scalable CPU computing environment for the development of next generation point dipole PFFs and for production simulations. Strategies linking Tinker-HP to Quantum Mechanics (QM) in the framework of multiscale polarizable self-consistent QM/MD simulations are also provided. The possibilities, performances and scalability of the software are demonstrated via benchmarks calculations using the polarizable AMOEBA force field on systems ranging from large water boxes of increasing size and ionic liquids to (very) large biosystems encompassing several proteins as well as the complete satellite tobacco mosaic virus and ribosome structures. For small systems, Tinker-HP appears to be competitive with the Tinker-OpenMM GPU implementation of Tinker. As the system size grows, Tinker-HP remains operational thanks to its access to distributed memory and takes advantage of its new algorithmic enabling for stable long timescale polarizable simulations. Overall, a several thousand-fold acceleration over a single-core computation is observed for the largest systems. The extension of the present CPU implementation of Tinker-HP to other computational platforms is discussed.

LICHEM: A QM/MM program for simulations with multipolar and polarizable force fields

We introduce an initial implementation of the LICHEM software package. LICHEM can interface with Gaussian, PSI4, NWChem, TINKER, and TINKER-HP to enable QM/MM calculations using multipolar/polarizable force fields. LICHEM extracts forces and energies from unmodified QM and MM software packages to perform geometry optimizations, single-point energy calculations, or Monte Carlo simulations. When the QM and MM regions are connected by covalent bonds, the pseudo-bond approach is employed to smoothly transition between the QM region and the polarizable force field. A series of water clusters and small peptides have been employed to test our initial implementation. The results obtained from these test systems show the capabilities of the new software and highlight the importance of including explicit polarization. © 2016 Wiley Periodicals, Inc.

Alternative Pathway for the Reaction Catalyzed by DNA Dealkylase AlkB from Ab Initio QM/MM Calculations

AlkB is the title enzyme of a family of DNA dealkylases that catalyze the direct oxidative dealkylation of nucleobases. The conventional mechanism for the dealkylation of N1-methyl adenine (1-meA) catalyzed by AlkB after the formation of FeIV-oxo is comprised by a reorientation of the oxo moiety, hydrogen abstraction, OH rebound from the Fe atom to the methyl adduct, and the dissociation of the resulting methoxide to obtain the repaired adenine base and formaldehyde. An alternative pathway with hydroxide as a ligand bound to the iron atom is proposed and investigated by QM/MM simulations. The results show OH- has a small impact on the barriers for the hydrogen abstraction and OH rebound steps. The effects of the enzyme and the OH- ligand on the hydrogen abstraction by the FeIV-oxo moiety are discussed in detail. The new OH rebound step is coupled with a proton transfer to the OH- ligand and results in a novel zwitterion intermediate. This zwitterion structure can also be characterized as Fe-O-C complex and facilitates the formation of formaldehyde. In contrast, for the pathway with H2O bound to iron, the hydroxyl product of the OH rebound step first needs to unbind from the metal center before transferring a proton to Glu136 or other residue/substrate. The consistency between our theoretical results and experimental findings is discussed. This study provides new insights into the oxidative repair mechanism of DNA repair by nonheme FeII and α-ketoglutarate (α-KG) dependent dioxygenases and a possible explanation for the substrate preference of AlkB.

Pseudobond parameters for QM/MM studies involving nucleosides, nucleotides, and their analogs

In biological systems involving nucleosides, nucleotides, or their respective analogs, the ribose sugar moiety is the most common reaction site, for example, during DNA replication and repair. However, nucleic bases, which comprise a sizable portion of nucleotide molecules, are usually unreactive during such processes. In quantum mechanical∕molecular simulations of nucleic acid reactivity, it may therefore be advantageous to describe specific ribosyl or ribosyl phosphate groups quantum mechanically and their respective nucleic bases with a molecular mechanics potential function. Here, we have extended the pseudobond approach to enable quantum mechanical∕molecular mechanical simulations involving nucleotides, nucleosides, and their analogs in which the interface between the two subsystems is located between the sugar and the base, namely, the C(sp(3))-N(sp(2)) bond. The pseudobond parameters were optimized on a training set of 10 molecules representing several nucleotide and nucleoside bases and analogs, and they were then tested on a larger test set of 20 diverse molecules. Particular emphasis was placed on providing accurate geometries and electrostatic properties, including electrostatic potential, natural bond orbital (NBO) and atoms in molecules (AIM) charges and AIM first moments. We also tested the optimized parameters on five nucleotide and nucleoside analogues of pharmaceutical relevance and a small polypeptide (triglycine). Accuracy was maintained for these systems, which highlights the generality and transferability of the pseudobond approach.

Ab initio QM/MM calculations show an intersystem crossing in the hydrogen abstraction step in dealkylation catalyzed by AlkB

AlkB is a bacterial enzyme that catalyzes the dealkylation of alkylated DNA bases. The rate-limiting step is known to be the abstraction of an H atom from the alkyl group on the damaged base by a Fe(IV)-oxo species in the active site. We have used hybrid ab initio quantum mechanical/molecular mechanical methods to study this step in AlkB. Instead of forming an Fe(III)-oxyl radical from Fe(IV)-oxo near the C-H activation transition state, the reactant is found to be an Fe(III)-oxyl with an intermediate-spin Fe (S = 3/2) ferromagnetically coupled to the oxyl radical, which we explore in detail using molecular orbital and quantum topological analyses. The minimum energy pathway remains on the quintet surface, but there is a transition between (IS)Fe(III)-oxyl and the state with a high-spin Fe (S = 5/2) antiferromagnetically coupled to the oxyl radical. These findings provide clarity for the evolution of the well-known π and σ channels on the quintet surface in the enzyme environment. Additionally, an energy decomposition analysis reveals nine catalytically important residues for the C-H activation step, some of which are conserved in two human homologues. These conserved residues are proposed as targets for experimental mutagenesis studies.

Liquid water simulations with the density fragment interaction approach

We reformulate the density fragment interaction (DFI) approach [Fujimoto and Yang, J. Chem. Phys., 2008, 129, 054102.] to achieve linear-scaling quantum mechanical calculations for large molecular systems. Two key approximations are developed to improve the efficiency of the DFI approach and thus enable the calculations for large molecules: the electrostatic interactions between fragments are computed efficiently by means of polarizable electrostatic-potential-fitted atomic charges; and frozen fragment pseudopotentials, similar to the effective fragment potentials that can be fitted from interactions between small molecules, are employed to take into account the Pauli repulsion effect among fragments. Our reformulated and parallelized DFI method demonstrates excellent parallel performance based on the benchmarks for the system of 256 water molecules. Molecular dynamics simulations for the structural properties of liquid water also show a qualitatively good agreement with experimental measurements including the heat capacity, binding energy per water molecule, and the radial distribution functions of atomic pairs of O-O, O-H, and H-H. With this approach, large-scale quantum mechanical simulations for water and other liquids become feasible.

Catalytic mechanism of 4-oxalocrotonate tautomerase: significances of protein-protein interactions on proton transfer pathways

4-Oxalocrotonate tautomerase (4-OT), a member of tautomerase superfamily, is an essential enzyme in the degradative metabolism pathway occurring in the Krebs cycle. The proton transfer process catalyzed by 4-OT has been explored previously using both experimental and theoretical methods; however, the elaborate catalytic mechanism of 4-OT still remains unsettled. By combining classical molecular mechanics with quantum mechanics, our results demonstrate that the native hexametric 4-OT enzyme, including six protein monomers, must be employed to simulate the proton transfer process in 4-OT due to protein-protein steric and electrostatic interactions. As a consequence, only three out of the six active sites in the 4-OT hexamer are observed to be occupied by three 2-oxo-4-hexenedioates (2o4hex), i.e., half-of-the-sites occupation. This agrees with experimental observations on negative cooperative effect between two adjacent substrates. Two sequential proton transfers occur: one proton from the C3 position of 2o4hex is initially transferred to the nitrogen atom of the general base, Pro1. Subsequently, the same proton is shuttled back to the position C5 of 2o4hex to complete the proton transfer process in 4-OT. During the catalytic reaction, conformational changes (i.e., 1-carboxyl group rotation) of 2o4hex may occur in the 4-OT dimer model but cannot proceed in the hexametric structure. We further explained that the docking process of 2o4hex can influence the specific reactant conformations and an alternative substrate (2-hydroxymuconate) may serve as reactant under a different reaction mechanism than 2o4hex.

Correlation between electron localization and metal ion mutagenicity in DNA synthesis from QM/MM calculations

DNA polymerases require two divalent metal ions in the active site for catalysis. Mg(2+) has been confirmed to be the most probable cation utilized by most polymerases in vivo. Other metal ions are either potent mutagens or inhibitors. We used structural and topological analyses based on ab initio QM/MM calculations to study human DNA polymerase λ (Polλ) with different metals in the active site. Our results indicate a slightly longer O3'-Pα distance (∼3.6 Å) for most inhibitor cations compared to the natural and mutagenic metals (∼3.3-3.4 Å). Optimization with a larger basis set for the previously reported transition state (TS) structures (Cisneros et al., DNA Repair, 2008, 7, 1824.) gives barriers of 17.4 kcal mol(-1) and 15.1 kcal mol(-1) for the Mg(2+) and Mn(2+) catalyzed reactions respectively. Relying on the key relation between the topological signature of a metal cation and its selectivity within biological systems (de Courcy et al., J. Chem. Theor. Comput., 2010, 6, 1048.) we have performed electron localization function (ELF) topological analyses. These analyses show that all inhibitor and mutagenic metals considered, except Na(+), present a "split" of the outer-shell density of the metal. This "splitting" is not observed for the non-mutagenic Mg(2+) metal. Population and multipole analyses on the ELF basins reveal that the electronic dipolar and quadrupolar polarization is significantly different with Mg(2+) compared to all other cations. Our results shed light at the atomic level on the subtle differences between Mg(2+), mutagenic, and inhibitor metals in DNA polymerases. These results provide a correlation between the electronic distribution of the cations in the active site and the possible consequences on DNA synthesis.

Autocatalytic intramolecular isopeptide bond formation in gram-positive bacterial pili: a QM/MM simulation

Many gram-positive pathogens possess external pili or fimbriae with which they adhere to host cells during the infection process. Unusual dual intramolecular isopeptide bonds between Asn and Lys side chains within the N-terminal and C-terminal domains of the pilus subunits have been observed initially in the Streptococcus pyogenes pilin subunit Spy0128 and subsequently in GBS52 from Streptococcus agalactiae, in the BcpA major pilin of Bacillus cereus and in the RrgB pilin of Streptococcus pneumoniae, among others. Within each pilin subunit, intramolecular isopeptide bonds serve to stabilize the protein. These bonds provide a means to withstand large external mechanical forces, as well as possibly assisting in supporting a conformation favored for pilin subunit polymerization via sortase transpeptidases. Genome-wide analyses of pili-containing gram-positive bacteria are known or suspected to contain isopeptide bonds in pilin subunits. For the autocatalytic formation of isopeptide cross-links, a conservation of three amino acids including Asn, Lys, and a catalytically important acidic Glu (or Asp) residue are responsible. However, the chemical mechanism of how isopeptide bonds form within pilin remains poorly understood. Although it is possible that several mechanistic paths could lead to isopeptide bond formation in pili, the requirement of a conserved glutamate and highly organized positioning of residues within the hydrophobic environment of the active site were found in numerous pilin crystal structures such as Spy0128 and RrgB. This suggests a mechanism involving direct coupling of lysine side chain amine to the asparagine carboxamide mediated by critical acid/base or hydrogen bonding interactions with the catalytic glutamate residue. From this mechanistic perspective, we used the QM/MM minimum free-energy path method to examine the reaction details of forming the isopeptide bonds with Spy0128 as a model pilin, specifically focusing on the role of the glutamate in catalysis. It was determined that the reaction mechanism likely consists of two major steps: the nucleophilic attack on Cγ by nitrogen in the unprotonated Lys ε-amino group and, then two concerted proton transfers occur during the formation of the intramolecular isopeptide bond to subsequently release ammonia. More importantly, within the dual active sites of Spy0128, Glu(117), and Glu(258) residues function as crucial catalysts for each isopeptide bond formation, respectively, by relaying two proton transfers. This work also suggests that domain-domain interactions within Spy0128 may modulate the reactivity of residues within each active site. Our results may hopefully shed light on the molecular mechanisms of pilin biogenesis in gram-positive bacteria.

Development and application of ab initio QM/MM methods for mechanistic simulation of reactions in solution and in enzymes

Determining the free energies and mechanisms of chemical reactions in solution and enzymes is a major challenge. For such complex reaction processes, combined quantum mechanics/molecular mechanics (QM/MM) method is the most effective simulation method to provide an accurate and efficient theoretical description of the molecular system. The computational costs of ab initio QM methods, however, have limited the application of ab initio QM/MM methods. Recent advances in ab initio QM/MM methods allowed the accurate simulation of the free energies for reactions in solution and in enzymes and thus paved the way for broader application of the ab initio QM/MM methods. We review here the theoretical developments and applications of the ab initio QM/MM methods, focusing on the determination of reaction path and the free energies of the reaction processes in solution and enzymes.