Free energy simulations of uncatalyzed DNA replication fidelity: structure and stability of T.G and dTTP.G terminal DNA mismatches flanked by a single dangling nucleotide

A reference system for DNA replication fidelity was studied by free energy perturbation (FEP) and linear interaction energy (LIE) methods. The studied system included a hydrated duplex DNA with the 5'-CG dangling end of the templating strand, and dCTP4-.Mg2+ or dTTP4-.Mg2+ inserted opposite the dangling G to form a correct (i.e., Watson-Crick) or incorrect (i.e., wobble) base pair, respectively. The average distance between the 3'-terminal oxygen of the primer strand and the alpha-phosphorus of dNTP was found to be 0.2 A shorter for the correct base pair than for the incorrect base pair. Binding of the incorrect dNTP was found to be disfavored by 0.4 kcal/mol relative to the correct dNTP. We estimated that improved binding and more near-attack configurations sampled by the correct base pair should translate in aqueous solution and in the absence of DNA polymerase into a six times faster rate for the incorporation of the correct dNTP into DNA. The accuracy of the calculated binding free energy difference was verified by examining the relative free energy for melting duplex DNA containing GC and GT terminal base pairs flanked by a 5' dangling C. The calculated LIE and FEP free energies of 1.7 and 1.1 kcal/mol, respectively, compared favorably with the experimental estimate of 1.4 kcal/mol obtained using the nearest neighbor parameters. To decompose the calculated free energies into additive electrostatic and van der Waals contributions and to provide a set of rigorous theoretical data for the parametrization of the LIE method, we suggested a variant of the FEP approach, for which we coined a binding-relevant free energy (BRFE) acronym. BRFE approach is characterized by its unique perturbation pathway and by its exclusion of the intramolecular energy of a rigid part of the ligand from the total potential energy.     

DNA polymerase beta catalysis: are different mechanisms possible?

DNA polymerases are crucial constituents of the complex cellular machinery for replicating and repairing DNA. Discerning mechanistic pathways of DNA polymerase on the atomic level is important for revealing the origin of fidelity discrimination. Mammalian DNA polymerase beta (pol beta), a small (39 kDa) member of the X-family, represents an excellent model system to investigate polymerase mechanisms. Here, we explore several feasible low-energy pathways of the nucleotide transfer reaction of pol beta for correct (according to Watson-Crick hydrogen bonding) G:C basepairing versus the incorrect G:G case within a consistent theoretical framework. We use mixed quantum mechanics/molecular mechanics (QM/MM) techniques in a constrained energy minimization protocol to effectively model not only the reactive core but also the influence of the rest of the enzymatic environment and explicit solvent on the reaction. The postulated pathways involve initial proton abstraction from the terminal DNA primer O3'H group, nucleophilic attack that extends the DNA primer chain, and elimination of pyrophosphate. In particular, we analyze several possible routes for the initial deprotonation step: (i) direct transfer to a phosphate oxygen O(Palpha) of the incoming nucleotide, (ii) direct transfer to an active site Asp group, and (iii) transfer to explicit water molecules. We find that the most probable initial step corresponds to step (iii), involving initial deprotonation to water, which is followed by proton migration to active site Asp residues, and finally to the leaving pyrophosphate group, with an activation energy of about 15 kcal/mol. We argue that initial deprotonation steps (i) and (ii) are less likely as they are at least 7 and 11 kcal/mol, respectively, higher in energy. Overall, the rate-determining step for both the correct and the incorrect nucleotide cases is the initial deprotonation in concert with nucleophilic attack at the phosphate center; however, the activation energy we obtain for the mismatched G:G case is 5 kcal/mol higher than that of the matched G:C complex, due to active site structural distortions. Taken together, our results support other reported mechanisms and help define a framework for interpreting nucleotide specificity differences across polymerase families, in terms of the concept of active site preorganization or the so-called "pre-chemistry avenue".     

Realistic simulations of proton transport along the gramicidin channel: demonstrating the importance of solvation effects

The nature of proton transduction (PTR) through a file of water molecules, along the gramicidin A (gA) channel, has long been considered as being highly relevant to PTR in biological systems. Previous attempts to model this process implied that the so-called Grotthuss mechanism and the corresponding orientation of the water file plays a major role. The present work reexamines the PTR in gA by combining a fully microscopic empirical valence bond (EVB) model and a recently developed simplified EVB-based model with Langevin dynamics (LD) simulations. The full model is used first to evaluate the free energy profile for a stepwise PTR process. The corresponding results are then used to construct the effective potential of the simplified EVB. This later model is then used in Langevin dynamics simulations, taking into account the correct physics of possible concerted motions and the effect of the solvent reorganization. The simulations reproduce the observed experimental trend and lead to a picture that is quite different from that assumed previously. It is found that the PTR in gA is controlled by the change in solvation energy of the transferred proton along the channel axis. Although the time dependent electrostatic fluctuations of the channel and water dipoles play their usual role in modulating the proton-transfer process (Proc. Natl. Acad. Sci. U.S.A. 1984, 81, 444), the PTR rate is mainly determined by the free energy profile. Furthermore, the energetics of the reorientation of the unprotonated water file do not appear to provide a consistent way of assessing the activation barrier for the PTR process. It seems to us that in the case of gA, and probably other systems with significant electrostatic barriers for the transfer of the proton charge, the PTR rate is controlled by the electrostatic barrier. This finding has clear consequences with regards to PTR processes in biological systems.     

Quantum mechanics/molecular mechanics investigation of the chemical reaction in Dpo4 reveals water-dependent pathways and requirements for active site reorganization

The nucleotidyl-transfer reaction coupled with the conformational transitions in DNA polymerases is critical for maintaining the fidelity and efficiency of DNA synthesis. We examine here the possible reaction pathways of a Y-family DNA polymerase, Sulfolobus solfataricus DNA polymerase IV (Dpo4), for the correct insertion of dCTP opposite 8-oxoguanine using the quantum mechanics/molecular mechanics (QM/MM) approach, both from a chemistry-competent state and a crystal closed state. The latter examination is important for understanding pre-chemistry barriers to interpret the entire enzyme mechanism, since the crystal closed state is not an ideal state for initiating the chemical reaction. The most favorable reaction path involves initial deprotonation of O3'H via two bridging water molecules to O1A, overcoming an overall potential energy barrier of approximately 20.0 kcal/mol. The proton on O1A-P(alpha) then migrates to the gamma-phosphate oxygen of the incoming nucleotide as O3' attacks P(alpha), and the P(alpha)-O3A bond breaks. The other possible pathway in which the O3'H proton is transferred directly to O1A on P(alpha) has an overall energy barrier of 25.0 kcal/mol. In both reaction paths, the rate-limiting step is the initial deprotonation, and the trigonal-bipyramidal configuration for P(alpha) occurs during the concerted bond formation (O3'-P(alpha)) and breaking (P(alpha)-O3A), indicating the associative nature of the chemical reaction. In contrast, the Dpo4/DNA complex with an imperfect active-site geometry corresponding to the crystal state must overcome a much higher activation energy barrier (29.0 kcal/mol) to achieve a tightly organized site due to hindered O3'H deprotonation stemming from larger distances and distorted conformation of the proton acceptors. This significant difference demonstrates that the pre-chemistry reorganization in Dpo4 costs approximately 4.0 to 9.0 kcal/mol depending on the primer terminus environment. Compared to the higher fidelity DNA polymerase beta from the X-family, Dpo4 has a higher chemical reaction barrier (20.0 vs 15.0 kcal/mol) due to the more solvent-exposed active site.     

Synthesis of Selenium‐Triphosphates (dNTPαSe) for More Specific DNA Polymerization

2′‐Deoxynucleoside 5′‐(alpha‐P‐seleno)‐triphosphates (dNTPαSe) have been conveniently synthesized using a protection‐free, one‐pot strategy. One of two diastereomers of each dNTPαSe can be efficiently recognized by DNA polymerases, while the other is neither a substrate nor an inhibitor. Furthermore, this Se‐atom modification can significantly inhibit non‐specific DNA polymerization caused by mis‐priming. Se–DNAs amplified with dNTPαSe via polymerase chain reaction have sequences identical to the corresponding native DNA. In conclusion, a simple strategy for more specific DNA polymerization has been established by replacing native dNTPs with dNTPαSe.     

Key role of active-site water molecules in bacteriorhodopsin proton-transfer reactions

The functional mechanism of the light-driven proton pump protein bacteriorhodopsin depends on the location of water molecules in the active site at various stages of the photocycle and on their roles in the proton-transfer steps. Here, free energy computations indicate that electrostatic interactions favor the presence of a cytoplasmic-side water molecule hydrogen bonding to the retinal Schiff base in the state preceding proton transfer from the retinal Schiff base to Asp85. However, the nonequilibrium nature of the pumping process means that the probability of occupancy of a water molecule in a given site depends both on the free energies of insertion of the water molecule in this and other sites during the preceding photocycle steps and on the kinetic accessibility of these sites on the time scale of the reaction steps. The presence of the cytoplasmic-side water molecule has a dramatic effect on the mechanism of proton transfer: the proton is channeled on the Thr89 side of the retinal, whereas the transfer on the Asp212 side is hindered. Reaction-path simulations and molecular dynamics simulations indicate that the presence of the cytoplasmic-side water molecule permits a low-energy bacteriorhodopsin conformer in which the water molecule bridges the twisted retinal Schiff base and the proton acceptor Asp85. From this low-energy conformer, proton transfer occurs via a concerted mechanism in which the water molecule participates as an intermediate proton carrier.     

Computational analysis of the proton translocation from Asp96 to schiff base in bacteriorhodopsin

The potential energy change during the M --> N process in bacteriorhodopsin has been evaluated by ab initio quantum chemical and advanced quantum chemical calculations following molecular dynamics (MD) simulations. Many previous experimental studies have suggested that the proton transfer from Asp96 to the Schiff base occurs under the following two conditions: (1) the hydrogen bond between Thr46 and Asp96 breaks and Thr46 is detached from Asp96 and (2) a stable chain of four water molecules spans an area from Asp96 --> Schiff base. In this work, we successfully reproduced the proton-transfer process occurring under these two conditions by molecular dynamics and quantum chemical calculations. The quantum chemical computation revealed that the proton transfer from Asp96 to Shiff base occurs in two-step reactions via an intermediate in which an H(3)O(+) appears around Ala215. The activation energy for the proton transfer in the first reaction was calculated to be 9.7 kcal/mol, which enables fast and efficient proton pump action. Further QM/MM (quantum mechanical/molecular mechanical) and FMO (fragment molecular orbital) calculations revealed that the potential energy change during the proton transfer is tightly regulated by the composition and the geometry of the surrounding amino acid residues of bacteriorhodopsin. Here, we report in detail the Asp96 --> Schiff base proton translocation mechanism of bacteriorhodopsin. Additionally, we discuss the effectiveness of combining quantum chemical calculations with truncated cluster models followed by advanced quantum chemical calculations applied to a whole protein to elucidate its reaction mechanism.     

Microscopic models for quantum mechanical calculations of chemical processes in solutions: LD/AMPAC and SCAAS/AMPAC calculations of solvation energies

Our previously developed approaches for integrating quantum mechanical molecular orbital methods with microscopic solvent models are refined and examined. These approaches consider the nonlinear solute–solvent coupling in a self-consistent way by incorporating the potential from the solvent dipoles in the solute Hamiltonian, while considering the polarization of the solvent by the potential from the solute charges. The solvent models used include the simplified Langevin Dipoles (LD) model and the much more expensive surface constrained All Atom Solvent (SCAAS) model, which is combined with a free energy pertubation (FEP) approach. Both methods are effectively integrated with the quantum mechanical AMPAC package and can be easily combined with other quantum mechanical programs. The advantages of the present approaches and their earlier versions over macroscopic reaction field models and supermolecular approaches are considered. A LD/MNDO study of solvated organic ions demonstrates that this model can yield reliable solvation energies, provided the quantum mechanical charges are scaled to have similar magnitudes to those obtained by high level ab initio methods. The incorporation of a field-dependent hydrophobic term in the LD free energy makes the present approach capable of evaluating the free energy of transfer of polar molecules from non polar solvents to aqueous solutions. The reliability of the LD approach is examined not only by evaluating a rather standard set of solvation energies of organic ions and polar molecules, but also by considering the stringent test case of sterically hindered hydrophobic ions. In this case, we compare the LD/MNDO solvation energies to the more rigorous FEP/SCAAS/MNDO solvation energies. Both methods are found to give similar results even in this challenging test case. The FEP/SCAAS/AMPAC method is incorporated into the current version of the program ENZYMIX. This option allows one to study chemical reactions in enzymes and in solutions using the MNDO and AM1 approximations. A special procedure that uses the EVB method as a reference potential for SCF MO calculations should help in improving the reliability of such studies.     

A quantum mechanical investigation of possible mechanisms for the nucleotidyl transfer reaction catalyzed by DNA polymerase beta

Several quantum mechanical (QM) and hybrid quantum/molecular mechanical (QM/MM) studies have been employed recently to analyze the nucleotidyl transfer reaction in DNA polymerase beta (pol beta). Our examination reveals strong dependence of the reported mechanism on the initial molecular model. Thus, we explore here several model systems by QM methods to investigate pol beta's possible pathway variations. Although our most favorable pathway involves a direct proton transfer from O3'(primer) to O2alpha(Palpha), we also discuss other initial proton-transfer steps--to an adjacent water, to triphosphate, or to aspartic units--and the stabilizing effect of crystallographic water molecules in the active site. Our favored reaction route has an energetically undemanding initial step of less than 1.0 kcal/mol (at the B3LYP/6-31G(d,p) level), and involves a slight rearrangement in the geometry of the active site. This is followed by two major steps: (1) direct proton transfer from O3'(primer) to O2alpha(Palpha) leading to the formation of a pentavalent, trigonal bipyramidal Palpha center, via an associative mechanism, at a cost of about 28 kcal/mol, and (2) breakage of the triphosphate unit (exothermic process, approximately 22 kcal/mol) that results in the full transfer of the nucleotide to the DNA and the formation of pyrophosphate. These energy values are expected to be lower in the physical system when full protein effects are incorporated. We also discuss variations from this dominant pathway, and their impact on the overall repair process. Our calculated barrier for the chemical reaction clearly indicates that chemistry is rate-limiting overall for correct nucleotide insertion in pol beta, in accord with other studies. Protonation studies on relevant intermediates suggest that, although protonation at a single aspartic residue may occur, the addition of a second proton to the system significantly disturbs the active site. We conclude that the active site rearrangement step necessary to attain a reaction-competent geometry is essential and closely related to the "pre-chemistry" avenue described recently as a key step in the overall kinetic cycle of DNA polymerases. Thus, our work emphasizes the many possible ways for DNA polymerase beta's chemical reaction to occur, determined by the active site environment and initial models.     

DFT study and Monte Carlo simulation on proton transfers of 2-amino-2-oxazoline, 2-amino-2-thiazoline, and 2-amino-2-imidazoline in the gas phase and in water

Density functional theory (DFT) and Monte Carlo (MC) simulation with free energy perturbation (FEP) techniques have been used to study the tautomeric proton transfer reaction of 2-amino-2-oxazoline, 2-amino-2-thiazoline, and 2-amino-2-imidazoline in the gas phase and in water. Two reaction pathways were considered: the direct and water-assisted transfers. The optimized structures and thermodynamic properties of stationary points for the title reaction system in the gas phase were calculated at the B3LYP/6-311+G(d, p) level of theory. The potential energy profiles along the minimum energy path in the gas phase and in water were obtained. The study of the solvent effect of water on the proton transfer of 2-amino-2-oxozoline, 2-amino-2-thiazoline, and 2-amino-2-imidazoline indicates that water as a solvent is favorable for the water-assisted process and slows down the rate of the direct transfer pathway.     

Prechemistry barriers and checkpoints do not contribute to fidelity and catalysis as long as they are not rate limiting

In the preceding article, “Perspective: Pre-chemistry conformational changes in DNA polymerase mechanisms” contributed by Schlick and coworkers as well as previous studies of these workers (Schlick et al. in Theor Chem Acc 131:1287, 2012; Radhakrishnan and Schlick in J Am Chem Soc 127:13245–13252, 2005; Radhakrishnan and Schlick in Biochem Biophys Res Commun 350:521–529, 2006; Radhakrishnan et al. in Biochemistry 45:15142–15156, 2006; Radhakrishnan and Schlick in Proc Natl Acad Sci USA 101:5970–5975, 2004) have argued that the conformational changes preceding the chemical step contribute to DNA synthesis and to the fidelity of DNA polymerases. In one of our previous investigations (Ram Prasad and Warshel in Proteins 79:2900–2919, 2011), we argued and showed that as long as the free energy barriers associated with any of the prechemistry steps are not rate limiting, they could not contribute to the catalysis and then to the fidelity. Though all our arguments are based on exact and well-defined scientific logics, Schlick and coworkers seem to overlook some of the clear conditions in these arguments and in particular the requirement that the chemical step is rate limiting in their arguments that the prechemistry barriers contribute to the catalysis. In fact, as long as the prechemistry steps are not rate limiting, we have shown that the enzymes cannot carry the memory of the previous steps. We also address other potential misunderstandings about several key issues; First, we clarify that it is misleading to relate the prechemistry proposal to the clear fact that the substrate-induced conformational changes determine the final preorganization (the issue is the height of the barrier of the enzyme substrate system and not the trivial fact that the enzyme has to change its structure when the substrate binds). Second, we address the presumed role of dynamical effects in enzyme catalysis and the assumption that any observable should be explored in studies of biological function even if they are not relevant to the given effect. Third, we clarify that the fidelity cannot be explained or quantified by invoking the induced fit or conformational selection effects but by evaluating the free energy contributions to the rate-limiting steps from the structures of the corresponding systems (that of course can reflect the induce fit structural changes). Overall, we put a major emphasis on clarifying what is the prechemistry proposal and thus on trying to force the reader to focus on the only real controversy. We of course dismiss any implication that our studies cannot explore mutational effects as we actually pioneered such computational studies and we clarify that in studies of chemical rates, the focus must be placed on evaluating the chemical barriers, rather than on irrelevant factors, but that the calculations of the chemical barriers must consider all the factors that determine this barrier (including metal ions) and also examine if needed different problematic proposals such as dynamical effects, tunneling, and prechemistry.     

Simulation of tunneling in enzyme catalysis by combining a biased propagation approach and the quantum classical path method: application to lipoxygenase

The ability of using wave function propagation approaches to simulate isotope effects in enzymes is explored, focusing on the large H/D kinetic isotope effect of soybean lipoxygenase-1 (SLO-1). The H/D kinetic isotope effect (KIE) is calculated as the ratio of the rate constants for hydrogen and deuterium transfer. The rate constants are calculated from the time course of the H and D nuclear wave functions. The propagations are done using one-dimensional proton potentials generated as sections from the full multidimensional surface of the reacting system in the protein. The sections are obtained during a classical empirical valence bond (EVB) molecular dynamics simulation of SLO-1. Since the propagations require an extremely long time for treating realistic activation barriers, it is essential to use an effective biasing approach. Thus, we develop here an approach that uses the classical quantum path (QCP) method to evaluate the quantum free energy change associated with the biasing potential. This approach provides an interesting alternative to full QCP simulations and to other current approaches for simulating isotope effects in proteins. In particular, this approach can be used to evaluate the quantum mechanical transmission factor or other dynamical effects, while still obtaining reliable quantized activation free energies due to the QCP correction.     

A Theoretical Investigation on N7(H)→N9(H)Tautomerization of Isolated and Monohydrated 2,6-Dithiopurine Combined with IPCM Solvent Model

The tautomerization reaction mechanism has been reported between N7(H) and N9(H) of isolated and monohydrated 2,6‐dithiopurine using B3LYP/6‐311+G(d,p). The isodensity polarized continuum model (IPCM) in the self‐consistent reaction field (SCRF) method is employed to account for the solvent effect of water on the tautomerization reaction activation energies. The results show that the two pathways P(1) (via the carbene intermediate I1) and P(2) (via the sp³‐hybrid intermediate I2) are found in intramolecular proton transfer, and each pathway is composed by two primary steps. The calculated activation energy barriers of the rate‐determining steps in isolated 2,6‐dithiopurine N7(H)→N9(H) tautomerism are 308.2 and 220.0 kJ·mol^?1 in the two pathways, respectively. Interestingly, in one‐water molecule catalyst, it dramatically lowers the N7(H)→N9(H) energy barriers by the concerted double proton transfer mechanism in P(1), favoring the formation of 2,6‐dithiopurine N9(H). However, the single proton transfer mechanism assisted with out‐of‐plane water in the first step of P(2) increases the activation energy barrier from 220.0 to 232.3 kJ·mol^?1, while the second step is the out‐of‐plane concerted double proton transfer mechanism, indicating that they will be less preferable for proton transfer. Additionally, the results also show that all the pathways are put into the aqueous solution, and the activation energy barriers have no significant changes. Therefore, the long‐range electrostatic effect of bulk solvent has no significant impact on proton transfer reactions and the interaction with explicit water molecules will significantly influence proton transfer reactions.     

Using the constrained DFT approach in generating diabatic surfaces and off diagonal empirical valence bond terms for modeling reactions in condensed phases

The empirical valence bond (EVB) model provides an extremely powerful way for modeling and analyzing chemical reactions in solutions and proteins. However, this model is based on the unverified assumption that the off diagonal elements of the EVB Hamiltonian do not change significantly upon transfer of the reacting system from one phase to another. This ad hoc assumption has been rationalized by its consistency with empirically observed linear free energy relationships, as well as by other qualitative considerations. Nevertheless, this assumption has not been rigorously established. The present work explores the validity of the above EVB key assumption by a rigorous numerical approach. This is done by exploiting the ability of the frozen density functional theory (FDFT) and the constrained density functional theory (CDFT) models to generate convenient diabatic states for QM/MM treatments, and thus to examine the relationship between the diabatic and adiabatic surfaces, as well as the corresponding effective off diagonal elements. It is found that, at least for the test case of S(N)()2 reactions, the off diagonal element does not change significantly upon moving from the gas phase to solutions and thus the EVB assumption is valid and extremely useful.     

A water-mediated and substrate-assisted catalytic mechanism for Sulfolobus solfataricus DNA polymerase IV

DNA polymerases are enzymes responsible for the synthesis of DNA from nucleotides. Understanding their molecular fundamentals is a prerequisite for elucidating their aberrant activities in diseases such as cancer. Here we have carried out ab initio quantum mechanical/molecular mechanical (QM/MM) studies on the nucleotidyl-transfer reaction catalyzed by the lesion-bypass DNA polymerase IV (Dpo4) from Sulfolobus solfataricus, with template guanine and Watson-Crick paired dCTP as the nascent base pair. The results suggested a novel water-mediated and substrate-assisted (WMSA) mechanism: the initial proton transfer to the alpha-phosphate of the substrate via a bridging crystal water molecule is the rate-limiting step, the nucleotidyl-transfer step is associative with a metastable pentacovalent phosphorane intermediate, and the pyrophosphate leaving is facilitated by a highly coordinated proton relay mechanism through mediation of water which neutralizes the evolving negative charge. The conserved carboxylates, which retain their liganding to the two Mg2+ ions during the reaction process, are found to be essential in stabilizing transition states. This WMSA mechanism takes specific advantage of the unique structural features of this low-fidelity lesion-bypass Y-family polymerase, which has a more spacious and solvent-exposed active site than replicative and repair polymerases.