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

Substrate recognition by norovirus polymerase: microsecond molecular dynamics study

Molecular dynamics simulations of complexes between Norwalk virus RNA dependent RNA polymerase and its natural CTP and 2dCTP (both containing the O5'-C5'-C4'-O4' sequence of atoms bridging the triphosphate and sugar moiety) or modified coCTP (C5'-O5'-C4'-O4'), cocCTP (C5'-O5'-C4'-C4'') substrates were produced by means of CUDA programmable graphical processing units and the ACEMD software package. It enabled us to gain microsecond MD trajectories clearly showing that similar nucleoside triphosphates can bind surprisingly differently into the active site of the Norwalk virus RNA dependent RNA polymerase. It corresponds to their different modes of action (CTP-substrate, 2dCTP-poor substrate, coCTP-chain terminator, cocCTP-inhibitor). Moreover, extremely rare events-as repetitive pervasion of Arg182 into a potentially reaction promoting arrangement-were captured.

Structural model of the Y-Family DNA polymerase V/RecA mutasome

To synthesize past DNA damaged by chemicals or radiation, cells have lesion bypass DNA polymerases (DNAPs), most of which are in the Y-Family. One class of Y-Family DNAPs includes DNAP η in eukaryotes and DNAP V in bacteria, which have low fidelity when replicating undamaged DNA. In Escherchia coli, DNAP V is carefully regulated to insure it is active for lesion bypass only, and one mode of regulation involves interaction of the polymerase subunit (UmuC) and two regulatory subunits (UmuD') with a RecA-filament bound to ss-DNA. Taking a docking approach, ∼150,000 unique orientations involving UmuC, UmuD' and RecA were evaluated to generate models, one of which was judged best able to rationalize the following published findings. (1) In the UmuD'(2)C/RecA-filament model, R64-UmuC interacts with S117-RecA, which is known to be at the UmuC/RecA interface. (2) At the model's UmuC/RecA interface, UmuC has three basic amino acids (K59/R63/R64) that anchor it to RecA. No other Y-Family DNAP has three basic amino acids clustered in this region, making it a plausible site for UmuC to form its unique interaction with RecA. (3) In the model, residues N32/N33/D34 of UmuC form a second interface with RecA, which is consistent with published findings. (4) Active UmuD' is generated when 24 amino acids in the N-terminal tail of UmuD are proteolyzed, which occurs when UmuD(2)C binds the RecA-filament. When UmuD is included in an UmuD(2)C/RecA-filament model, plausible UmuD/RecA contacts guide the UmuD cleavage site (C24/G25) into the UmuD proteolysis active site (S60/K97). One contact involves E11-UmuD interacting with R243-RecA, where the latter is known to be important for UmuD cleavage. (5) The UmuD(2)C/RecA-filament model rationalizes published findings that at least some UmuD-to-UmuD' cleavage occurs intermolecularly. (6) Active DNAP V is known to be the heterotetramer UmuD'(2)C/RecA, a model of which can be generated by a simple rearrangement of the RecA monomer at the 3'-end of the RecA-filament. The rearranged UmuD'(2)C/RecA model rationalizes published findings about UmuD' residues in proximity to RecA. In summary, docking and molecular simulations are used to develop an UmuD'(2)C/RecA model, whose structure rationalizes much of the known properties of the active form of DNA polymerase V.

Why does the G117H mutation considerably improve the activity of human butyrylcholinesterase against sarin? Insights from quantum mechanical/molecular mechanical free energy calculations

Human butyrylcholinesterase (BChE) is recognized as the most promising bioscavenger for organophosphorus (OP) warfare nerve agents. The G117H mutant of human BChE has been identified as a potential catalytic bioscavenger with a remarkably improved activity against OP nerve agents such as sarin, but it still does not satisfy the clinical use. For further design of the higher-activity mutants against OP nerve agents, it is essential to understand how the G117H mutation improves the activity. The reaction mechanisms and the free energy profiles for spontaneous reactivation of wild-type BChE and its G117H mutant phosphorylated by sarin have been explored, in this study, by performing first-principles quantum mechanical/molecular mechanical free energy calculations, and the remarkable role of the G117H mutation on the activity has been elucidated. For both the wild-type and G117H mutant enzymes, H438 acts as a general base to initiate the spontaneous reactivation that consists of two reaction steps: the nucleophilic attack at the phosphorus by a water molecule and decomposition of the pentacoordinated phosphorus intermediate. The calculated overall free energy barriers, i.e., 30.2 and 23.9 kcal/mol for the wild type and G117H mutant, respectively, are in good agreement with available experimental kinetic data. On the basis of the calculated results, the mutated residue (H117 in the G117H mutant) cannot initiate the spontaneous reactivation as a general base. Instead, it skews the oxyanion hole and makes the phosphorus more open to the nucleophilic water molecule, resulting in a remarkable change in the rate-determining step and significantly improved catalytic activity of human BChE.

Preferred WMSA catalytic mechanism of the nucleotidyl transfer reaction in human DNA polymerase κ elucidates error-free bypass of a bulky DNA lesion

Human DNA Pol κ is a polymerase enzyme, specialized for near error-free bypass of certain bulky chemical lesions to DNA that are derived from environmental carcinogens present in tobacco smoke, automobile exhaust and cooked food. By employing ab initio QM/MM-MD (Quantum Mechanics/Molecular Mechanics-Molecular Dynamics) simulations with umbrella sampling, we have determined the entire free energy profile of the nucleotidyl transfer reaction catalyzed by Pol κ and provided detailed mechanistic insights. Our results show that a variant of the Water Mediated and Substrate Assisted (WMSA) mechanism that we previously deduced for Dpo4 and T7 DNA polymerases is preferred for Pol κ as well, suggesting its broad applicability. The hydrogen on the 3'-OH primer terminus is transferred through crystal and solvent waters to the γ-phosphate of the dNTP, followed by the associative nucleotidyl transfer reaction; this is facilitated by a proton transfer from the γ-phosphate to the α,β-bridging oxygen as pyrophosphate leaves, to neutralize the evolving negative charge. MD simulations show that the near error-free incorporation of dCTP opposite the major benzo[a]pyrene-derived dG lesion is compatible with the WMSA mechanism, allowing for an essentially undisturbed pentacovalent phosphorane transition state, and explaining the bypass of this lesion with little mutation by Pol κ.

Revelation of a catalytic calcium-binding site elucidates unusual metal dependence of a human apyrase

Human soluble calcium-activated nucleotidase 1 (hSCAN-1) represents a new family of apyrase enzymes that catalyze the hydrolysis of nucleotide di- and triphosphates, thereby modulating extracellular purinergic and pyrimidinergic signaling. Among well-characterized phosphoryl transfer enzymes, hSCAN-1 is unique not only in its unusual calcium-dependent activation, but also in its novel phosphate-binding motif. Its catalytic site does not utilize backbone amide groups to bind phosphate, as in the common P-loop, but contains a large cluster of acidic ionizable side chains. By employing a state-of-the-art computational approach, we have revealed a previously uncharacterized catalytic calcium-binding site in hSCAN-1, which elucidates the unusual calcium-dependence of its apyrase activity. In a high-order coordination shell, the newly identified calcium ion organizes the active site residues to mediate nucleotide binding, to orient the nucleophilic water, and to facilitate the phosphoryl transfer reaction. From ab initio QM/MM molecular dynamics simulations with umbrella sampling, we have characterized a reverse protonation catalytic mechanism for hSCAN-1 and determined its free energy reaction profile. Our results are consistent with available experimental studies and provide new detailed insight into the structure-function relationship of this novel calcium-activated phosphoryl transfer enzyme.

Watching DNA polymerase η make a phosphodiester bond

DNA synthesis has been extensively studied, but the chemical reaction itself has not been visualized. Here we follow the course of phosphodiester bond formation using time-resolved X-ray crystallography. Native human DNA polymerase η, DNA and dATP were co-crystallized at pH 6.0 without Mg(2+). The polymerization reaction was initiated by exposing crystals to 1 mM Mg(2+) at pH 7.0, and stopped by freezing at desired time points for structural analysis. The substrates and two Mg(2+) ions are aligned within 40 s, but the bond formation is not evident until 80 s. From 80 to 300 s structures show a mixture of decreasing substrate and increasing product of the nucleotidyl-transfer reaction. Transient electron densities indicate that deprotonation and an accompanying C2'-endo to C3'-endo conversion of the nucleophile 3'-OH are rate limiting. A third Mg(2+) ion, which arrives with the new bond and stabilizes the intermediate state, may be an unappreciated feature of the two-metal-ion mechanism.

Reaction Pathway and Free Energy Profile for Cocaine Hydrolase-Catalyzed Hydrolysis of (-)-Cocaine

Reaction pathway of (-)-cocaine hydrolysis catalyzed by our recently discovered most efficient cocaine hydrolase, which is the A199S/F227A/S287G/A328W/Y332G mutant of human butyrylcholinesterase (BChE), and the corresponding free energy profile have been studied by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Based on the QM/MM-FE results, the catalytic hydrolysis process consists of four major reaction steps, including the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by hydroxyl group of S198, dissociation of (-)-cocaine benzoyl ester, nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by water, and finally the dissociation between (-)-cocaine benzoyl group and S198 of the enzyme. The second reaction step is rate-determining. The calculated free energy barrier associated with the transition state for the rate-determining step is ~15.0 kcal/mol, which is in excellent agreement with the experimentally-derived activation free energy of ~14.7 kcal/mol. The mechanistic insights obtained from the present study will be valuable for rational design of more active cocaine hydrolase against (-)-cocaine. In particular, future efforts aiming at further increasing the catalytic activity of the enzyme against (-)-cocaine should focus on stabilization of the transition state for the second reaction step in which the benzoyl ester of (-)-cocaine dissociates.

Prechemistry versus preorganization in DNA replication fidelity

The molecular origin of nucleotide insertion catalysis and fidelity of DNA polymerases is explored by means of computational simulations. Special attention is paid to the examination of the validity of proposals that invoke prechemistry effects, checkpoints concepts, and dynamical effects. The simulations reproduce the observed fidelity in Pol β, starting with the relevant observed X-ray structures of the complex with the right (R) and wrong (W) nucleotides. The generation of free energy surfaces for the R and W systems also allowed us to analyze different proposals about the origin of the fidelity and to reach several important conclusions. It is found that the potential of mean force (PMF) obtained by proper sampling does not support QM/MM-based proposals of a large barrier before the prechemistry state. Furthermore, examination of dynamical proposals by the renormalization approach indicates that the motions from open to close configurations do not contribute to catalysis or fidelity. Finally we discuss and analyze the induced fit concept and show that, despite its importance, it does not explain fidelity. That is, the fidelity is apparently due to the change in the preorganization of the chemical site, as a result of the relaxation of the binding site upon binding of the incorrect nucleotide. Finally and importantly, since the issue is the barrier associated with the enzyme-substrate (ES)/DNA complex at the chemical transition state and not the path to this complex formation (unless this path involves rate determining steps), it is also not useful to invoke checkpoints while discussing fidelity.

An abridged transition state model to derive structure, dynamics, and energy components of DNA polymerase β fidelity

We show how a restricted reaction surface can be used to facilitate the calculation of biologically important contributions of active site geometries and dynamics to DNA polymerase fidelity. Our analysis, using human DNA polymerase beta (pol β), is performed within the framework of an electrostatic linear free energy response (EFER) model. The structure, dynamics, and energetics of pol β-DNA-dNTP interactions are computed between two points on the multidimensional reaction free energy surface. "Point 1" represents a ground state activation intermediate (GSA), which is obtained by deprotonating the terminal 3'OH group of the primer DNA strand. "Point 2" is the transition state (PTS) for the attack of the 3'O(-) (O(nuc)) on the P(α) atom of dNTP substrate, having the electron density of a dianionic phosphorane intermediate. Classical molecular dynamics simulations are used to compute the geometric and dynamic contributions to the formation of right and wrong O(nuc)-P chemical bonds. Matched dCTP·G and mismatched dATP·G base pairs are used to illustrate the analysis. Compared to the dCTP·G base pair, the dATP·G mismatch has fewer GSA configurations with short distances between O(nuc) and P(α) atoms and between the oxygen in the scissile P-O bond (O(lg)) and the nearest structural water. The thumb subdomain conformation of the GSA complex is more open for the mismatch, and the H-bonds in the mispair become more extended during the nucleophilic attack than in the correct pair. The electrostatic contributions of pol β and DNA residues to catalysis of the right and wrong P-O(nuc) bond formation are 5.3 and 3.1 kcal/mol, respectively, resulting in an 80-fold contribution to fidelity. The EFER calculations illustrate the considerable importance of Arg183 and an O(lg)-proximal water molecule to pol β fidelity.

Catalytic mechanism of RNA backbone cleavage by ribonuclease H from quantum mechanics/molecular mechanics simulations

We use quantum mechanics/molecular mechanics simulations to study the cleavage of the ribonucleic acid (RNA) backbone catalyzed by ribonuclease H. This protein is a prototypical member of a large family of enzymes that use two-metal catalysis to process nucleic acids. By combining Hamiltonian replica exchange with a finite-temperature string method, we calculate the free energy surface underlying the RNA-cleavage reaction and characterize its mechanism. We find that the reaction proceeds in two steps. In a first step, catalyzed primarily by magnesium ion A and its ligands, a water molecule attacks the scissile phosphate. Consistent with thiol-substitution experiments, a water proton is transferred to the downstream phosphate group. The transient phosphorane formed as a result of this nucleophilic attack decays by breaking the bond between the phosphate and the ribose oxygen. In the resulting intermediate, the dissociated but unprotonated leaving group forms an alkoxide coordinated to magnesium ion B. In a second step, the reaction is completed by protonation of the leaving group, with a neutral Asp132 as a likely proton donor. The overall reaction barrier of ∼15 kcal mol(-1), encountered in the first step, together with the cost of protonating Asp132, is consistent with the slow measured rate of ∼1-100/min. The two-step mechanism is also consistent with the bell-shaped pH dependence of the reaction rate. The nonmonotonic relative motion of the magnesium ions along the reaction pathway agrees with X-ray crystal structures. Proton-transfer reactions and changes in the metal ion coordination emerge as central factors in the RNA-cleavage reaction.

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.

Reaction mechanism for cocaine esterase-catalyzed hydrolyses of (+)- and (-)-cocaine: unexpected common rate-determining step

First-principles quantum mechanical/molecular mechanical free energy calculations have been performed to examine the catalytic mechanism for cocaine esterase (CocE)-catalyzed hydrolysis of (+)-cocaine in comparison with CocE-catalyzed hydrolysis of (-)-cocaine. It has been shown that the acylation of (+)-cocaine consists of nucleophilic attack of the hydroxyl group of Ser117 on the carbonyl carbon of (+)-cocaine benzoyl ester and the dissociation of (+)-cocaine benzoyl ester. The first reaction step of deacylation of (+)-cocaine, which is identical to that of (-)-cocaine, is rate-determining, indicating that CocE-catalyzed hydrolyses of (+)- and (-)-cocaine have a common rate-determining step. The computational results predict that the catalytic rate constant of CocE against (+)-cocaine should be the same as that of CocE against (-)-cocaine, in contrast with the remarkable difference between human butyrylcholinesterase-catalyzed hydrolyses of (+)- and (-)-cocaine. The prediction has been confirmed by experimental kinetic analysis on CocE-catalyzed hydrolysis of (+)-cocaine in comparison with CocE-catalyzed hydrolysis of (-)-cocaine. The determined common rate-determining step indicates that rational design of a high-activity mutant of CocE should be focused on the first reaction step of the deacylation. Furthermore, the obtained mechanistic insights into the detailed differences in the acylation between the (+)- and (-)-cocaine hydrolyses provide indirect clues for rational design of amino acid mutations that could more favorably stabilize the rate-determining transition state in the deacylation and, thus, improve the catalytic activity of CocE. This study provides a valuable mechanistic base for rational design of an improved esterase for therapeutic treatment of cocaine abuse.

Probing deuterium isotope effect on structure and solvation dynamics of human serum albumin

The deuterium isotopic effect on the structure and solvation dynamics of the protein, human serum albumin (HSA), has been studied by using circular dichroism (CD), femtosecond up-conversion, FRET, and single-molecule spectroscopy. The CD spectra suggest that D(2)O affects the structure of HSA, leading to a 20% decrease in the helical structure. The FRET study indicates that the distance of C153 from the lone tryptophan residue of HSA is quite similar (≈21 Å) in H(2)O and D(2)O, and hence, the location of the probe in the protein remains the same in the two solvents. The single-molecule study suggests that coumarin 153 (C153) binds almost exclusively (>96%) to one site of HSA. Solvation dynamics of C153 in HSA is found to be markedly retarded in D(2)O compared with H(2)O. In H(2)O, the solvation of C153 bound to HSA is found to be biexponential with one component of 7 ps (30%) and a long component of 350 ps (70%). In D(2)O, we detected a short component of 4 ps (41%) and a long component of 950 ps (59%). Thus, the ultraslow component of the solvation dynamics of C153 bound to HSA in D(2)O (950 ps) is 2.5-fold slower than that in H(2)O (350 ps). The marked deuterium isotope effect has been ascribed to water molecules confined in the protein environment and to a lesser extent to the structural modification of protein by D(2)O.

The glmS ribozyme: use of a small molecule coenzyme by a gene-regulatory RNA

The glmS ribozyme is the first known example of a natural ribozyme that has evolved to require binding of an exogenous small molecule for activity. In Gram-positive bacteria, this RNA domain is part of the messenger RNA (mRNA) encoding the essential enzyme that synthesizes glucosamine-6-phosphate (GlcN6P). When present at physiologic concentration, this small molecule binds to the glmS ribozyme and uncovers a latent self-cleavage activity that ultimately leads to degradation of the mRNA. Biochemical and structural studies reveal that the RNA adopts a rigid fold stabilized by three pseudoknots and the packing of a peripheral domain against the ribozyme core. GlcN6P binding to this pre-organized RNA does not induce conformational changes; rather, the small molecule functions as a coenzyme, providing a catalytically essential amine group to the active site. The ribozyme is not a passive player, however. Active site functional groups are essential for catalysis, even in the presence of GlcN6P. In addition to being a superb experimental system with which to analyze how RNA catalysts can exploit small molecule coenzymes to broaden their chemical versatility, the presence of the glmS ribozyme in numerous pathogenic bacteria make this RNA an attractive target for the development of new antibiotics and antibacterial strategies.

Architecture of y-family DNA polymerases relevant to translesion DNA synthesis as revealed in structural and molecular modeling studies

DNA adducts, which block replicative DNA polymerases (DNAPs), are often bypassed by lesion-bypass DNAPs, which are mostly in the Y-Family. Y-Family DNAPs can do non-mutagenic or mutagenic dNTP insertion, and understanding this difference is important, because mutations transform normal into tumorigenic cells. Y-Family DNAP architecture that dictates mechanism, as revealed in structural and modeling studies, is considered. Steps from adduct blockage of replicative DNAPs, to bypass by a lesion-bypass DNAP, to resumption of synthesis by a replicative DNAP are described. Catalytic steps and protein conformational changes are considered. One adduct is analyzed in greater detail: the major benzo[a]pyrene adduct (B[a]P-N²-dG), which is bypassed non-mutagenically (dCTP insertion) by Y-family DNAPs in the IV/κ-class and mutagenically (dATP insertion) by V/η-class Y-Family DNAPs. Important architectural differences between IV/κ-class versus V/η-class DNAPs are discussed, including insights gained by analyzing ~400 sequences each for bacterial DNAPs IV and V, along with sequences from eukaryotic DNAPs kappa, eta and iota. The little finger domains of Y-Family DNAPs do not show sequence conservation; however, their structures are remarkably similar due to the presence of a core of hydrophobic amino acids, whose exact identity is less important than the hydrophobic amino acid spacing.

Atomistic details of the associative phosphodiester cleavage in human ribonuclease H

During translation of the genetic information of DNA into proteins, mRNA is synthesized by RNA polymerase and after the transcription process degraded by RNase H. The endoribonuclease RNase H is a member of the nucleotidyl-transferase (NT) superfamily and is known to hydrolyze the phosphodiester bonds of RNA which is hybridized to DNA. Retroviral RNase H is part of the viral reverse transcriptase enzyme that is indispensable for the proliferation of retroviruses, such as HIV. Inhibitors of this enzyme could therefore provide new drugs against diseases like AIDS. In our study we investigated the molecular mechanism of RNA cleavage by human RNase H using a comprehensive high level DFT/B3LYP QM/MM theoretical method for the calculation of the stationary points and nudged elastic band (NEB) and free energy calculations to identify the transition state structures, the rate limiting step and the reaction barrier. Our calculations reveal that the catalytic mechanism proceeds in two steps and that the nature of the nucleophile is a water molecule. In the first step, the water attack on the scissile phosphorous is followed by a proton transfer from the water to the O2P oxygen and a trigonal bipyramidal pentacoordinated phosphorane is formed. Subsequently, in the second step the proton is shuttled to the O3' oxygen to generate the product state. During the reaction mechanism two Mg(2+) ions support the formation of a stable associated in-line S(N)2-type phosphorane intermediate. Our calculated energy barrier of 19.3 kcal mol(-1) is in excellent agreement with experimental findings (20.5 kcal mol(-1)). These results may contribute to the clarification and understanding of the RNase H reaction mechanism and of further enzymes from the RNase family.

On the role of water in peroxidase catalysis: a theoretical investigation of HRP compound I formation

We have investigated the dynamics of water molecules in the distal pocket of horseradish peroxidase to elucidate the role that they may play in the formation of the principal active species of the enzymatic cycle (compound I, Por(o+)-Fe(IV)=O) upon reaction of the resting Fe(III) state with hydrogen peroxide. The equilibrium molecular dynamics simulations show that, in accord with experimental evidence, the active site access channel is hydrated with an average of two to three water molecules within 5 A from the bound hydrogen peroxide. Although the channel is always hydrated, the specific conformations in which a water molecule bridges H(2)O(2) and the distal histidine, which were found (Derat; et al. J. Am. Chem. Soc. 2007, 129, 6346.) to display a low-energy barrier for the initial acid-base step of the reaction, occur with low probability but are relevant within the time scale of catalysis. Metadynamics simulations, which were used to reconstruct the free-energy landscape of water motion in the access channel, revealed that preferred interaction sites within the channel are separated by small energy barriers (<1.5 kcal/mol). Most importantly, water-bridged conformations lie on a shoulder just 1 kcal/mol above one local minimum and thus are easily accessible. Such an energy landscape appears as a requisite for the effectiveness of compound I formation, whereby the H-bonding pattern involving reactants and catalytic residues (including the intervening water molecule) has to rearrange to deliver the proton to the distal OH moiety of the hydrogen peroxide and thereby lead to heterolytic O-O cleavage. Our study provides an example of a system for which the "reactive configurations" (i.e., structures characterized by a low barrier for the chemical transformation) correspond to a minor population of the system and show how equilibrium molecular dynamics and free-energy calculations may conveniently be used to ascertain that such reactive conformations are indeed accessible to the system. Once again, the MD and QM/MM combination shows that a single water molecule acts as a biocatalyst in the cycle of HRP.

Influence of local sequence context on damaged base conformation in human DNA polymerase iota: molecular dynamics studies of nucleotide incorporation opposite a benzo[a]pyrene-derived adenine lesion

Human DNA polymerase iota is a lesion bypass polymerase of the Y family, capable of incorporating nucleotides opposite a variety of lesions in both near error-free and error-prone bypass. With undamaged templating purines polymerase iota normally favors Hoogsteen base pairing. Polymerase iota can incorporate nucleotides opposite a benzo[a]pyrene-derived adenine lesion (dA*); while mainly error-free, the identity of misincorporated bases is influenced by local sequence context. We performed molecular modeling and molecular dynamics simulations to elucidate the structural basis for lesion bypass. Our results suggest that hydrogen bonds between the benzo[a]pyrenyl moiety and nearby bases limit the movement of the templating base to maintain the anti glycosidic bond conformation in the binary complex in a 5'-CAGA*TT-3' sequence. This facilitates correct incorporation of dT via a Watson-Crick pair. In a 5'-TTTA*GA-3' sequence the lesion does not form these hydrogen bonds, permitting dA* to rotate around the glycosidic bond to syn and incorporate dT via a Hoogsteen pair. With syn dA*, there is also an opportunity for increased misincorporation of dGTP. These results expand our understanding of the versatility and flexibility of polymerase iota and its lesion bypass functions in humans.

Reaction pathway and free-energy barrier for reactivation of dimethylphosphoryl-inhibited human acetylcholinesterase

The dephosphorylation/reactivation mechanism and the corresponding free-energy profile of the dimethylphosphoryl-inhibited conjugate of human acetylcholinesterase (AChE) has been studied by performing first-principles quantum mechanical/molecular mechanical free-energy (QM/MM-FE) calculations. On the basis of the QM/MM-FE results, for the favorable reaction pathway, the entire dephosphorylation/reactivation process consists of three reaction steps, including the nucleophilic water attack on the P atom, the spatial reorganization of the dimethylphosphoryl group, and the dissociation between the dimethylphosphoryl group and Ser203 of AChE. The overall free-energy barrier for the entire dephosphorylation/reactivation reaction is found to be the free-energy change from the initial reactant to the transition state associated with the spatial reorganization step, and the calculated overall free-energy barrier (20.1 to 23.5 kcal/mol) is reasonably close to the experimentally derived activation free energy of 22.3 kcal/mol. In addition, key amino acid residues and their specific roles in the reaction process have been identified.

Reaction pathway and free-energy barrier for reactivation of dimethylphosphoryl-inhibited human acetylcholinesterase

The dephosphorylation/reactivation mechanism and the corresponding free-energy profile of the dimethylphosphoryl-inhibited conjugate of human acetylcholinesterase (AChE) has been studied by performing first-principles quantum mechanical/molecular mechanical free-energy (QM/MM-FE) calculations. On the basis of the QM/MM-FE results, for the favorable reaction pathway, the entire dephosphorylation/reactivation process consists of three reaction steps, including the nucleophilic water attack on the P atom, the spatial reorganization of the dimethylphosphoryl group, and the dissociation between the dimethylphosphoryl group and Ser203 of AChE. The overall free-energy barrier for the entire dephosphorylation/reactivation reaction is found to be the free-energy change from the initial reactant to the transition state associated with the spatial reorganization step, and the calculated overall free-energy barrier (20.1 to 23.5 kcal/mol) is reasonably close to the experimentally derived activation free energy of 22.3 kcal/mol. In addition, key amino acid residues and their specific roles in the reaction process have been identified.

Fundamental reaction mechanism and free energy profile for (-)-cocaine hydrolysis catalyzed by cocaine esterase

The fundamental reaction mechanism of cocaine esterase (CocE)-catalyzed hydrolysis of (-)-cocaine and the corresponding free energy profile have been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations. On the basis of the QM/MM-FE results, the entire hydrolysis reaction consists of four reaction steps, including the nucleophilic attack on the carbonyl carbon of (-)-cocaine benzoyl ester by the hydroxyl group of Ser117, dissociation of (-)-cocaine benzoyl ester, nucleophilic attack on the carbonyl carbon of (-)-cocaine benzoyl ester by water, and finally dissociation between the (-)-cocaine benzoyl group and Ser117 of CocE. The third reaction step involving the nucleophilic attack of a water molecule was found to be rate-determining, which is remarkably different from (-)-cocaine hydrolysis catalyzed by wild-type butyrylcholinesterase (BChE; where the formation of the prereactive BChE-(-)-cocaine complex is rate-determining) or its mutants containing Tyr332Gly or Tyr332Ala mutation (where the first chemical reaction step is rate-determining). Besides, the role of Asp259 in the catalytic triad of CocE does not follow the general concept of the "charge-relay system" for all serine esterases. The free energy barrier calculated for the rate-determining step of CocE-catalyzed hydrolysis of (-)-cocaine is 17.9 kcal/mol, which is in good agreement with the experimentally derived activation free energy of 16.2 kcal/mol. In the present study, where many sodium ions are present, the effects of counterions are found to be significant in determining the free energy barrier. The finding of the significant effects of counterions on the free energy barrier may also be valuable in guiding future mechanistic studies on other charged enzymes.

Impact of conformational heterogeneity of OxoG lesions and their pairing partners on bypass fidelity by Y family polymerases

7,8-Dihydro-8-oxoguanine (oxoG), the predominant oxidative DNA damage lesion, is processed differently by high-fidelity and Y-family lesion bypass polymerases. Although high-fidelity polymerases extend predominantly from an A base opposite an oxoG, the Y-family polymerases Dpo4 and human Pol eta preferentially extend from the oxoG*C base pair. We have determined crystal structures of extension Dpo4 ternary complexes with oxoG opposite C, A, G, or T and the next nascent base pair. We demonstrate that neither template backbone nor the architecture of the active site is perturbed by the oxoG(anti)*C and oxoG*A pairs. However, the latter manifest conformational heterogeneity, adopting both oxoG(syn)*A(anti) and oxoG(anti)*A(syn) alignment. Hence, the observed reduced primer extension from the dynamically flexible 3'-terminal primer base A is explained. Because of homology between Dpo4 and Pol eta, such a dynamic screening mechanism might be utilized by Dpo4 and Pol eta to regulate error-free versus error-prone bypass of oxoG and other lesions.

Artificial reaction coordinate "tunneling" in free-energy calculations: the catalytic reaction of RNase H

We describe a method for the systematic improvement of reaction coordinates in quantum mechanical/molecular mechanical (QM/MM) calculations of reaction free-energy profiles. In umbrella-sampling free-energy calculations, a biasing potential acting on a chosen reaction coordinate is used to sample the system in reactant, product, and transition states. Sharp, nearly discontinuous changes along the resulting reaction path are used to identify coordinates that are relevant for the reaction but not properly sampled. These degrees of freedom are then included in an extended reaction coordinate. The general formalism is illustrated for the catalytic cleavage of the RNA backbone of an RNA/DNA hybrid duplex by the RNase H enzyme of Bacillus halodurans. We find that in the initial attack of the phosphate diester by water, the oxygen-phosphorus distances alone are not sufficient as reaction coordinates, resulting in substantial hysteresis in the proton degrees of freedom and a barrier that is too low (approximately 10 kcal/mol). If the proton degrees of freedom are included in an extended reaction coordinate, we obtain a barrier of 21.6 kcal/mol consistent with the experimental rates. As the barrier is approached, the attacking water molecule transfers one of its protons to the O1P oxygen of the phosphate group. At the barrier top, the resulting hydroxide ion forms a penta-coordinated phosphate intermediate. The method used to identify important degrees of freedom, and the procedure to optimize the reaction coordinate are general and should be useful both in classical and in QM/MM free-energy calculations.

Visualizing sequence-governed nucleotide selectivities and mutagenic consequences through a replicative cycle: processing of a bulky carcinogen N2-dG lesion in a Y-family DNA polymerase

Understanding how DNA polymerases process lesions remains fundamental to determining the molecular origins of mutagenic translesion bypass. We have investigated how a benzo[a]pyrene-derived N(2)-dG adduct, 10S-(+)-trans-anti-[BP]-N(2)-dG ([BP]G*), is processed in Dpo4, the well-characterized Y-family bypass DNA polymerase. This polymerase has a slippage-prone spacious active site region. Experimental results in a 5'-C[BP]G*G-3' sequence context reveal significant selectivity for dGTP insertion that predominantly yields -1 deletion extension products. A less pronounced error-prone nonslippage pathway that leads to full extension products with insertion of A > C > G opposite the lesion is also observed. Molecular modeling and dynamics simulations follow the bypass of [BP]G* through an entire replication cycle for the first time in Dpo4, providing structural interpretations for the experimental observations. The preference for dGTP insertion is explained by a 5'-slippage pattern in which the unmodified G rather than G* is skipped, the incoming dGTP pairs with the C on the 5'-side of G*, and the -1 deletion is produced upon further primer extension which is more facile than nucleotide insertion. In addition, the simulations suggest that the [BP]G* may undergo an anti/syn conformational rearrangement during the stages of the replication cycle. In the minor nonslippage pathway, the nucleotide insertion preferences opposite the lesion are explained by relative distortions to the active site region. These structural insights, provided by the modeling and dynamics studies, augment kinetic and limited available crystallographic investigations with bulky lesions, by providing molecular explanations for lesion bypass activities over an entire replication cycle.

Amino acid architecture that influences dNTP insertion efficiency in Y-family DNA polymerase V of E. coli

Y-family DNA polymerases (DNAPs) are often required in cells to synthesize past DNA-containing lesions, such as [+ta]-B[a]P-N(2)-dG, which is the major adduct of the potent mutagen/carcinogen benzo[a]pyrene. The current model for the non-mutagenic pathway in Escherichia coli involves DNAP IV inserting deoxycytidine triphosphate opposite [+ta]-B[a]P-N(2)-dG and DNAP V doing the next step(s), extension. We are investigating what structural differences in these related Y-family DNAPs dictate their functional differences. X-ray structures of Y-family DNAPs reveal a number of interesting features in the vicinity of the active site, including (1) the "roof-amino acid" (roof-aa), which is the amino acid that lies above the nucleobase of the deoxynucleotide triphosphate (dNTP) and is expected to play a role in dNTP insertion efficiency, and (2) a cluster of three amino acids, including the roof-aa, which anchors the base of a loop, whose detailed structure dictates several important mechanistic functions. Since no X-ray structures existed for UmuC (the polymerase subunit of DNAP V) or DNAP IV, we previously built molecular models. Herein, we test the accuracy of our UmuC(V) model by investigating how amino acid replacement mutants affect lesion bypass efficiency. A ssM13 vector containing a single [+ta]-B[a]P-N(2)-dG is transformed into E. coli carrying mutations at I38, which is the roof-aa in our UmuC(V) model, and output progeny vector yield is monitored as a measure of the relative efficiency of the non-mutagenic pathway. Findings show that (1) the roof-aa is almost certainly I38, whose beta-carbon branching R-group is key for optimal activity, and (2) I38/A39/V29 form a hydrophobic cluster that anchors an important mechanistic loop, aa29-39. In addition, bypass efficiency is significantly lower both for the I38A mutation of the roof-aa and for the adjacent A39T mutation; however, the I38A/A39T double mutant is almost as active as wild-type UmuC(V), which probably reflects the following. Y-family DNAPs fall into several classes with respect to the [roof-aa/next amino acid]: one class has [isoleucine/alanine] and includes UmuC(V) and DNAP eta (from many species), while the second class has [alanine (or serine)/threonine] and includes DNAP IV, DNAP kappa (from many species), and Dpo4. Thus, the high activity of the I38A/A39T double mutant probably arises because UmuC(V) was converted from the V/eta class to the IV/kappa class with respect to the [roof-aa/next amino acid]. Structural and mechanistic aspects of these two classes of Y-family DNAPs are discussed.

Calculating solution redox free energies with ab initio quantum mechanical/molecular mechanical minimum free energy path method

A quantum mechanical/molecular mechanical minimum free energy path (QM/MM-MFEP) method was developed to calculate the redox free energies of large systems in solution with greatly enhanced efficiency for conformation sampling. The QM/MM-MFEP method describes the thermodynamics of a system on the potential of mean force surface of the solute degrees of freedom. The molecular dynamics (MD) sampling is only carried out with the QM subsystem fixed. It thus avoids "on-the-fly" QM calculations and thus overcomes the high computational cost in the direct QM/MM MD sampling. In the applications to two metal complexes in aqueous solution, the new QM/MM-MFEP method yielded redox free energies in good agreement with those calculated from the direct QM/MM MD method. Two larger biologically important redox molecules, lumichrome and riboflavin, were further investigated to demonstrate the efficiency of the method. The enhanced efficiency and uncompromised accuracy are especially significant for biochemical systems. The QM/MM-MFEP method thus provides an efficient approach to free energy simulation of complex electron transfer reactions.

Proton transfer in the mechanism of polyadenylate polymerase

PAP (polyadenylate polymerase) is the template-independent RNA polymerase responsible for synthesis of the 3' poly(A) tails of mRNA. To investigate the role of proton transfer in the catalytic mechanism of PAP, the pH dependence of the steady-state kinetic parameters of yeast PAP were determined for the forward (adenyl transfer) and reverse (pyrophosphorolysis) reactions. The results indicate that productive formation of an enzyme-RNA-MgATP complex is pH independent over a broad pH range, but that formation of an active enzyme-RNA-MgPPi complex is strongly pH dependent, consistent with the production of a proton on the enzyme in the forward reaction. The pH dependence of the maximum velocity of the forward reaction suggests two protonic species are involved in enzyme catalysis. Optimal enzyme activity requires one species to be protonated and the other deprotonated. The deuterium solvent isotope effect on Vmax is also consistent with proton transfer involved in catalysis of a rate-determining step. Finally, pKa calculations of PAP were performed by the MCCE (multiconformational continuum electrostatic) method. Together, the data support that the protonation of residues Lys215 and Tyr224 exhibit co-operativity that is important for MgATP2- and MgPPi2- binding/dissociation, and suggest these residues function in electrostatic, but not in general acid, catalysis.

Quantum mechanics/molecular mechanics investigation of the mechanism of phosphate transfer in human uridine-cytidine kinase 2

The mechanisms of enzyme-catalyzed phosphate transfer and hydrolysis reactions are of great interest due to their importance and abundance in biochemistry. The reaction may proceed in a stepwise fashion, with either a pentavalent phosphorane or a metaphosphate anion intermediate, or by a concerted SN2 mechanism. Despite much theoretical work focused on a few key enzymes, a consensus for the mechanism has not been reached, and examples of all three possibilities have been demonstrated. We have investigated the mechanism of human uridine-cytidine kinase 2 (UCK2, EC 2.7.1.48), which catalyzes the transfer of a phosphate group from ATP to the ribose 5'-hydroxyl of cytidine and uridine. UCK2 is normally expressed in human placenta, but is overexpressed in certain cancer cells, where it is responsible for activating a class of antitumor prodrugs. The UCK2 mechanism was investigated by generating a 2D potential energy surface as a function of the P-O bonds forming and breaking, with energies calculated using a quantum mechanics/molecular mechanics potential (B3LYP/6-31G(d):AMBER). The mechanism of phosphate transfer is shown to be concerted, and is accompanied by concerted proton transfer from the 5'-hydroxyl to a conserved active site aspartic acid that serves as a catalytic base. The calculated barrier for this reaction is 15.1 kcal/mol, in relatively good agreement with the experimental barrier of 17.5 kcal/mol. The interactions of the enzyme active site with the reactant, transition state, and product are examined for their implications on the design of anticancer prodrugs or positron emission tomography (PET) reporter probes for this enzyme.

Polymerase-tailored variations in the water-mediated and substrate-assisted mechanism for nucleotidyl transfer: insights from a study of T7 DNA polymerase

The nucleotidyl transfer reaction catalyzed by DNA polymerases is the critical step governing the accurate transfer of genetic information during DNA replication, and its malfunctioning can cause mutations leading to human diseases, including cancer. Here, utilizing ab initio quantum mechanical/molecular mechanical calculations with free-energy perturbation, we carried out an extensive investigation of the nucleotidyl transfer reaction mechanism in the well-characterized high-fidelity replicative DNA polymerase from phage T7. Our defined mechanism entails an initial concerted deprotonation of a conserved crystal water molecule with protonation of the gamma-phosphate of the deoxynucleotide triphosphate(dNTP) via a solvent water molecule, and then the proton on the primer 3'-terminus is transferred to the resulting hydroxide ion. Subsequently, the nucleophilic attack takes place, with the formation of a metastable pentacovalent phosphorane intermediate. Finally, the pyrophosphate leaves, facilitated by the relay of the proton on the gamma-phosphate to the alpha-beta bridging oxygen via solvent water. The computed activation free-energy barrier is consistent with kinetic data for the chemistry step with correct nucleotide incorporation in T7 DNA polymerase. This variant of the water-mediated and substrate-assisted mechanism has features tailored to the structure of the T7 DNA polymerase. However, a unifying theme in the water-mediated and substrate-assisted mechanism is the cycling through crystal and solvent water molecules of the proton originating from the primer 3'-terminus to the alpha-beta bridging oxygen of the deoxynucleotide triphosphate; this neutralizes the evolving negative charge as pyrophosphate leaves and restores the polymerase to its pre-chemistry state. These unifying features are likely requisite elements for nucleotidyl transfer reactions.

Reaction mechanism of the epsilon subunit of E. coli DNA polymerase III: insights into active site metal coordination and catalytically significant residues

The 28 kDa epsilon subunit of Escherichia coli DNA polymerase III is the exonucleotidic proofreader responsible for editing polymerase insertion errors. Here, we study the mechanism by which epsilon carries out the exonuclease activity. We performed quantum mechanics/molecular mechanics calculations on the N-terminal domain containing the exonuclease activity. Both the free-epsilon and a complex epsilon bound to a theta homologue (HOT) were studied. For the epsilon-HOT complex Mg(2+) or Mn(2+) were investigated as the essential divalent metal cofactors, while only Mg(2+) was used for free-epsilon. In all calculations a water molecule bound to the catalytic metal acts as the nucleophile for hydrolysis of the phosphate bond. Initially, a direct proton transfer to H162 is observed. Subsequently, the nucleophilic attack takes place followed by a second proton transfer to E14. Our results show that the reaction catalyzed with Mn(2+) is faster than that with Mg(2+), in agreement with experiment. In addition, the epsilon-HOT complex shows a slightly lower energy barrier compared to free-epsilon. In all cases the catalytic metal is observed to be pentacoordinated. Charge and frontier orbital analyses suggest that charge transfer may stabilize the pentacoordination. Energy decomposition analysis to study the contribution of each residue to catalysis suggests that there are several important residues. Among these, H98, D103, D129, and D146 have been implicated in catalysis by mutagenesis studies. Some of these residues were found to be structurally conserved on human TREX1, the exonuclease domains from E. coli DNA-Pol I, and the DNA polymerase of bacteriophage RB69.

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.

Accommodation of an N-(deoxyguanosin-8-yl)-2-acetylaminofluorene adduct in the active site of human DNA polymerase iota: Hoogsteen or Watson-Crick base pairing?

Bypass across DNA lesions by specialized polymerases is essential for maintenance of genomic stability. Human DNA polymerase iota (poliota) is a bypass polymerase of the Y family. Crystal structures of poliota suggest that Hoogsteen base pairing is employed to bypass minor groove DNA lesions, placing them on the spacious major groove side of the enzyme. Primer extension studies have shown that poliota is also capable of error-free nucleotide incorporation opposite the bulky major groove adduct N-(deoxyguanosin-8-yl)-2-acetylaminofluorene (dG-AAF). We present molecular dynamics simulations and free energy calculations suggesting that Watson-Crick base pairing could be employed in poliota for bypass of dG-AAF. In poliota with Hoogsteen-paired dG-AAF the bulky AAF moiety would reside on the cramped minor groove side of the template. The Hoogsteen-capable conformation distorts the active site, disrupting interactions necessary for error-free incorporation of dC opposite the lesion. Watson-Crick pairing places the AAF rings on the spacious major groove side, similar to the position of minor groove adducts observed with Hoogsteen pairing. Watson-Crick-paired structures show a well-ordered active site, with a near reaction-ready ternary complex. Thus our results suggest that poliota would utilize the same spacious region for lesion bypass of both major and minor groove adducts. Therefore, purine adducts with bulk on the minor groove side would use Hoogsteen pairing, while adducts with the bulky lesion on the major groove side would utilize Watson-Crick base pairing as indicated by our MD simulations for dG-AAF. This suggests the possibility of an expanded role for poliota in lesion bypass.

Highly dissociative and concerted mechanism for the nicotinamide cleavage reaction in Sir2Tm enzyme suggested by ab initio QM/MM molecular dynamics simulations

Sir2 enzymes catalyze the NAD+-dependent protein deacetylation and play critical roles in epigenetics, cell death, and lifespan regulation. In spite of a current flurry of experimental studies, the catalytic mechanism for this unique and important class of enzymes remains elusive. Employing on-the-fly Born-Oppenheimer molecular dynamics simulations with the B3LYP/6-31G(d) QM/MM potential and the umbrella sampling method, we have characterized the initial step of the Sir2Tm-catalyzed reaction, which is also the most controversial portion of its mechanism. Our results indicate that the nicotinamide cleavage reaction employs a highly dissociative and concerted displacement mechanism: the cleavage of the glycosidic bond is facilitated by the nucleophilic participation of the acetyl-lysine, and the dissociative transition state has a significant oxocarbenium ion character. During this step of the reaction, the Sir2Tm enzyme strongly stabilizes the covalent O-alkylamidate intermediate whereas its effect on the transition state is quite minimal. In addition, functional roles of key residues and motifs have been elucidated. This work further demonstrates the feasibility and applicability of the state-of-the-art ab initio QM/MM molecular dynamics approach in simulating enzyme reactions.

Y-Family DNA polymerases may use two different dNTP shapes for insertion: a hypothesis and its implications

Chemicals and radiation can damage DNA leading to the formation of adducts/lesions, which - if not removed by DNA repair pathways - usually block replicative DNA polymerases (DNAPs). To overcome such potentially lethal blockage, cells have lesion bypass DNAPs, which are often in the Y-Family and include several classes. One class includes human DNAP kappa and E. coli DNAP IV, and they insert dCTP in the non-mutagenic pathway opposite [+ta]-B[a]P-N(2)-dG, which is the major adduct formed by the environmental carcinogen benzo[a]pyrene. Another class includes hDNAP eta and ecDNAP V, and they insert dATP opposite [+ta]-B[a]P-N(2)-dG in the dominant G-->T mutagenic pathway. Herein we develop a hypothesis for why the IV/kappa-class preferentially does cellular dCTP insertion. On the minor groove side of the active site, Y-Family DNAPs have a cleft/hole that can be analyzed based on an analogy to a "chimney." Our models of DNAP IV show a large chimney opening from which the pyrene of [+ta]-B[a]P-N(2)-dG can protrude, which allows canonical adduct-dG:dCTP pairing. In contrast, our models of DNAP V have small chimney openings that forces adduct-dG downward in the active site such that canonical adduct-dG:dCTP pairing is not possible. Based on X-ray structures, sequence alignment and our modeled structures of Y-Family DNAPs, chimney opening size seems primarily controlled by one amino acid ("flue-handle"), which dictates whether nearby amino acids ("flue") plug the chimney or not. Based on this analysis, a correlation is apparent: the flue is closed in V/eta-class DNAPs giving small chimney openings, while the flue is open for the IV/kappa-class giving large chimney openings. Secondarily, a hypothesis is developed for why the V/eta-class might preferentially do cellular dATP insertion opposite [+ta]-B[a]P-N(2)-dG: the small chimney forces adduct-dG lower in the active site, possibly leading to catalysis using a non-canonical dNTP shape that permits syn-adenine:adduct-dG base pairing. In summary, a hypothesize is developed that the pyrene moiety of [+ta]-B[a]P-N(2)-dG protrudes from the large chimney opening of DNAP IV, thus permitting canonical dCTP:adduct-dG pairing, while the small chimney opening of DNAP V forces [+ta]-B[a]P-N(2)-dG lower down in the active site, in which syn-adenine can pair with adduct-dG via a non-canonical dNTP shape.

The N-clasp of human DNA polymerase kappa promotes blockage or error-free bypass of adenine- or guanine-benzo[a]pyrenyl lesions

DNA bypass polymerases are utilized to transit bulky DNA lesions during replication, but the process frequently causes mutations. The structural origins of mutagenic versus high fidelity replication in lesion bypass is therefore of fundamental interest. As model systems, we investigated the molecular basis of the experimentally observed essentially faithful bypass of the guanine 10S-(+)-trans-anti-benzo[a]pyrene-N(2)-dG adduct by the Y-family human DNA polymerase kappa, and the observed blockage of pol kappa produced by the adenine 10S-(+)-trans-anti-benzo[a]pyrene-N(2)-dA adduct. These lesions are derived from the most tumorigenic metabolite of the ubiquitous cancer-causing pollutant, benzo[a]pyrene. We compare our results for the dG adduct with our earlier studies for the pol kappa archaeal homolog Dpo4, which processes the same lesion in an error-prone manner. Molecular modeling, molecular mechanics calculations and molecular dynamics simulations were utilized. Our results show that the pol kappa N-clasp is a key structural feature that accounts for the dA adduct blockage and the near-error-free bypass of the dG lesion. Absence of the N-clasp in Dpo4 explains the error-prone processing of the same lesion by this enzyme. Thus, our studies elucidate structure-function relationships in the fidelity of lesion bypass.

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.

Catalytic mechanism of human DNA polymerase lambda with Mg2+ and Mn2+ from ab initio quantum mechanical/molecular mechanical studies

DNA polymerases play a crucial role in the cell cycle due to their involvement in genome replication and repair. Understanding the reaction mechanism by which these polymerases carry out their function can provide insights into these processes. Recently, the crystal structures of human DNA polymerase lambda (Pollambda) have been reported both for pre- and post-catalytic complexes [García-Díaz et al., DNA Repair 3 (2007), 1333]. Here we employ the pre-catalytic complex as a starting structure for the determination of the catalytic mechanism of Pollambda using ab initio quantum mechanical/molecular mechanical methods. The reaction path has been calculated using Mg(2+) and Mn(2+) as the catalytic metals. In both cases the reaction proceeds through a two-step mechanism where the 3'-OH of the primer sugar ring is deprotonated by one of the conserved Asp residues (D490) in the active site before the incorporation of the nucleotide to the nascent DNA chain. A significant charge transfer is observed between both metals and some residues in the active site as the reaction proceeds. The optimized reactant and product structures agree with the reported crystal structures. In addition, the calculated reaction barriers for both metals are close to experimentally estimated barriers. Energy decomposition analysis to explain individual residue contributions suggests that several amino acids surrounding the active site are important for catalysis. Some of these residues, including R420, R488 and E529, have been implicated in catalysis by previous mutagenesis experiments on the homologous residues on Polbeta. Furthermore, Pollambda residues R420 and E529 found to be important from the energy decomposition analysis, are homologous to residues R183 and E295 in Polbeta, both of which are linked to cancer. In addition, residues R386, E391, K422 and K472 appear to have an important role in catalysis and could be a potential target for mutagenesis experiments. There is partial conservation of these residues across the Pol X family of DNA polymerases.

Dineopentyl phosphate hydrolysis: evidence for stepwise water attack

Phosphate ester hydrolysis is ubiquitous in biology, playing a central role in energy production, signaling, biosynthesis, and the regulation of protein function among other things. Although the mechanism of action of the enzymes regulating this reaction has been the focus of intensive research in the past few decades, the correct description of this apparently simple reaction remains controversial. A clear understanding of the mechanism that takes place in solution is crucial to be able to evaluate whether proposals for the enzyme-catalyzed mechanisms are reasonable. For the pH-independent hydrolysis of phosphate diesters, several kinetically equivalent mechanisms are plausible, including hydroxide attack on the neutral phosphate. However, it is very difficult to measure the rate of this reaction directly by experimental methods, so it has been evaluated by examining the rate of hydrolysis of neutral phosphate triesters, where a methyl group has replaced a proton. This may not be an accurate model of the neutral phosphate diester and does not provide information about a reaction pathway that is concerted with nucleophilic attack to generate a similar phosphorane. We have carefully mapped out free energy surfaces for both hydroxide and water attack on the dineopentyl phosphate anion and for water attack on the neutral diester. In doing so, we have accurately reproduced existing experimental data and demonstrate that water attack proceeds through an associative mechanism with proton transfer to the phosphate to generate a phosphorane intermediate. Our data show that the substrate-as-base mechanism is viable for phosphate ester hydrolysis, which may have important implications for the studies of phosphate ester hydrolysis by enzymes.

The associative nature of adenylyl transfer catalyzed by T4 DNA ligase

DNA ligase seals nicks in dsDNA using chemical energy of the phosphoanhydride bond in ATP or NAD(+) and assistance of a divalent metal cofactor Mg(2+). Molecular details of ligase catalysis are essential for understanding the mechanism of metal-promoted phosphoryl transfer reactions in the living cell responsible for a wide range of processes, e.g., DNA replication and transcription, signaling and differentiation, energy coupling and metabolism. Here we report a single-turnover (31)P solid-state NMR study of adenylyl transfer catalyzed by DNA ligase from bacteriophage T4. Formation of a high-energy covalent ligase-nucleotide complex is triggered in situ by the photo release of caged Mg(2+), and sequentially formed intermediates are monitored by NMR. Analyses of reaction kinetics and chemical-shift changes indicate that the pentacoordinated phosphorane intermediate builds up to 35% of the total reacting species after 4-5 h of reaction. This is direct experimental evidence of the associative nature of adenylyl transfer catalyzed by DNA ligase. NMR spectroscopy in rotating solids is introduced as an analytical tool for recording molecular movies of reaction processes. Presented work pioneers a promising direction in structural studies of biochemical transformations.

Probing the structures of leishmanial farnesyl pyrophosphate synthases: homology modeling and docking studies

Leishmania donovani and Leishmania major farnesyl pyrophosphate synthase ( LdFPPS and LmFPPS) are potential targets for the development of antileishmanial therapy. The protein sequence for LdFPPS was recently elucidated in our laboratory. Highly refined homology models were generated using the protein sequences of LdFPPS and the closely related LmFPPS enzyme. A ligand-refined model of LmFPPS with a bound bisphosphonate ligand was generated using restraint-guided molecular mechanics followed by quantum mechanics/molecular mechanics refinement. The ligand-refined model of LmFPPS was further validated through extensive pose validation, enrichment, and other docking studies involving known bisphosphonate inhibitors. The model was able to explain the critical binding site interactions and site-directed mutagenesis data obtained from experimental studies on related FPPS enzymes. The ligand-refined model in conjunction with the validated docking protocol could be utilized in the future for structure-based virtual screening and rational drug design studies against these targets.

Lesion processing: high-fidelity versus lesion-bypass DNA polymerases

When a high-fidelity DNA polymerase encounters certain DNA-damage sites, its progress can be stalled and one or more lesion-bypass polymerases are recruited to transit the lesion. Here, we consider two representative types of lesions: (i) 7,8-dihydro-8-oxoguanine (8-oxoG), a small, highly prevalent lesion caused by oxidative damage; and (ii) bulky lesions derived from the environmental pre-carcinogen benzo[a]pyrene, in the high-fidelity DNA polymerase Bacillus fragment (BF) from Bacillus stearothermophilus and in the lesion-bypass DNA polymerase IV (Dpo4) from Sulfolobus solfataricus. The tight fit of the BF polymerase around the nascent base pair contrasts with the more spacious, solvent-exposed active site of Dpo4, and these differences in architecture result in distinctions in their respective functions: one-step versus stepwise polymerase translocation, mutagenic versus accurate bypass of 8-oxoG, and polymerase stalling versus mutagenic bypass at bulky benzo[a]pyrene-derived lesions.

Ab initio QM/MM studies of the phosphoryl transfer reaction catalyzed by PEP mutase suggest a dissociative metaphosphate transition state

The interconversion between phosphoenolpyruvate (PEP) and phosphonopyruvate (P-pyr) catalyzed by PEP mutase is investigated using an ab initio QM/MM method with the QM region treated at the B3LYP/6-31G* level of theory. Two-dimensional minimum energy path calculations were carried out for both the wild-type enzyme and the N122A mutant. The calculations suggest a dissociative transition state featuring metaphosphate and Mg(2+)-coordinating pyruvate enolate, stabilized by an extensive hydrogen bond network involving Asn122, Ser123, Arg159, His190, Ser46, and Leu48. It is also found that a substantial conformational change in the pyruvyl group is required for the interconversion.

On the interpretation of the observed linear free energy relationship in phosphate hydrolysis: a thorough computational study of phosphate diester hydrolysis in solution

The hydrolysis of phosphate esters is crucially important to biological systems, being involved in, among other things, signaling, energy transduction, biosynthesis, and the regulation of protein function. Despite this, there are many questions that remain unanswered in this important field, particularly with regard to the preferred mechanism of hydrolysis of phosphate esters, which can proceed through any of multiple pathways that are either associative or dissociative in nature. Previous comparisons of calculated and observed linear free energy relationships (LFERs) for phosphate monoester dianions with different leaving groups showed that the TS character gradually changes from associative to dissociative with the increasing acidity of the leaving group, while reproducing the experimental LFER. Here, we have generated ab initio potential energy surfaces for the hydrolysis of phosphate diesters in solution, with a variety of leaving groups. Once again, the reaction changes from a compact concerted pathway to one that is more expansive in character when the acidity of the leaving group increases. When such systems are examined in solution, it is essential to take into consideration the contribution of solute to the overall activation entropy, which remains a major computational challenge. The popular method of calculating the entropy using a quasi-harmonic approximation appears to markedly overestimate the configurational entropy for systems with multiple occupied energy wells. We introduce an improved restraint release approach for evaluating configurational entropies and apply this approach to our systems. We demonstrate that when this factor is taken into account, it is possible to reproduce the experimental LFER for this system with reasonable accuracy.