The mechanism of the translocation step in DNA replication by DNA polymerase I: a computer simulation analysis

Structures of the human leading strand Polε-PCNA holoenzyme

In eukaryotes, the leading strand DNA is synthesized by Polε and the lagging strand by Polδ. These replicative polymerases have higher processivity when paired with the DNA clamp PCNA. While the structure of the yeast Polε catalytic domain has been determined, how Polε interacts with PCNA is unknown in any eukaryote, human or yeast. Here we report two cryo-EM structures of human Polε-PCNA-DNA complex, one in an incoming nucleotide bound state and the other in a nucleotide exchange state. The structures reveal an unexpected three-point interface between the Polε catalytic domain and PCNA, with the conserved PIP (PCNA interacting peptide)-motif, the unique P-domain, and the thumb domain each interacting with a different protomer of the PCNA trimer. We propose that the multi-point interface prevents other PIP-containing factors from recruiting to PCNA while PCNA functions with Polε. Comparison of the two states reveals that the finger domain pivots around the [4Fe-4S] cluster-containing tip of the P-domain to regulate nucleotide exchange and incoming nucleotide binding.

Formation and Effects of Upstream DNA-RNA Base Pairing in Telomerase

Telomere elongation by telomerase consists of two types of translocation: duplex translocation during each repeat synthesis and template translocation at the end of repeat synthesis. Our replica exchange molecular dynamics simulations show that in addition to the Watson-Crick interactions in the active site, templating RNA can also form base pairs with the upstream regions of DNA, mostly with the second upstream DNA repeat with respect to the 3'-end. At the end of the repeat synthesis, dG10-P442 and dG11-N446 hydrogen bonds form. Then, active-site base pairs dissociate one by one, and the RNA bases reanneal with the complementary base on the upstream DNA repeat. For each dissociating base pair a new one forms, thus conserving the number of base pairs during template translocation.

Translocation pause of remdesivir-containing primer/template RNA duplex within SARS-CoV-2’s RNA polymerase complexes

The mechanism of remdesivir incorporation into the RNA primer by the RNA-dependent RNA polymerase (RdRp) of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) remains to be fully established at the molecular level. Here, we compare molecular dynamics (MD) simulations after incorporation of either remdesivir monophosphate (RMP) or adenosine monophosphate (AMP). We find that the Mg²⁺-pyrophosphate (PPi) binds more tightly to the polymerase when the added RMP is at the third primer position than in the AMP added complex. The increased affinity of Mg²⁺-PPi to the RMP-added primer/template (P/T) RNA duplex complex introduces a new hydrogen bond of a substituted cyano group in RMP with the K593 sidechain. The new interactions disrupt a switching mechanism of a hydrogen bond network that is essential for translocation of the P/T duplex product and for opening of a vacant NTP-binding site necessary for next primer extension. Furthermore, steric interactions between the sidechain of S861 and the 1′-cyano group of RMP at position i+3 hinders translocation of RMP to the i + 4 position, where i labels the insertion site. These findings are particularly valuable to guide the design of more effective inhibitors of SARS-CoV-2 RNA polymerase.

Transient State Kinetics of Plasmodium falciparum Apicoplast DNA Polymerase Suggests the Involvement of Accessory Factors for Efficient and Accurate DNA Synthesis

Plasmodium, the causative agent of malaria, belongs to the phylum Apicomplexa. Most apicomplexans, including Plasmodium, contain an essential nonphotosynthetic plastid called the apicoplast that harbors its own genome that is replicated by a dedicated organellar replisome. This replisome employs a single DNA polymerase (apPol), which is expected to perform both replicative and translesion synthesis. Unlike other replicative polymerases, no processivity factor for apPol has been identified. While preliminary structural and biochemical studies have provided an overall characterization of apPol, the kinetic mechanism of apPol's activity remains unknown. We have used transient state methods to determine the kinetics of replicative and translesion synthesis by apPol and show that apPol has low processivity and efficiency while copying undamaged DNA. Moreover, while apPol can bypass oxidatively damaged lesions, the bypass is error-prone. Taken together, our results raise the following question─how does a polymerase with low processivity, efficiency, and fidelity (for translesion synthesis) faithfully replicate the apicoplast organellar DNA within the hostile environment of the human host? We hypothesize that interactions with putative components of the apicoplast replisome and/or an as-yet-undiscovered processivity factor transform apPol into an efficient and accurate enzyme.

The Rate-limiting Step of DNA Synthesis by DNA Polymerase Occurs in the Fingers-closed Conformation

DNA polymerases maintain genomic integrity by copying DNA with high fidelity, part of which relies on the polymerase fingers opening-closing transition, a series of conformational changes during the DNA synthesis reaction cycle. Fingers opening and closing has been challenging to study, mainly due to the need to synchronise molecular ensembles. We previously studied fingers opening-closing on single polymerase-DNA complexes using single-molecule FRET; however, our work was limited to pre-chemistry reaction steps. Here, we advance our analysis to extensible substrates, and observe DNA polymerase (Pol) conformational changes across the entire DNA polymerisation reaction in real-time, gaining direct access to an elusive post-chemistry step rate-limiting for DNA synthesis. Our results showed that Pol adopts the fingers-closed conformation during polymerisation, and that the post-chemistry rate-limiting step occurs in the fingers-closed conformation. We found that fingers-opening in the Pol-DNA binary complex in the absence of polymerisation is slow (∼5.3 s-1), and comparable to the rate of fingers-opening after polymerisation (3.4 s-1); this indicates that the fingers-opening step itself could be largely responsible for the slow post-chemistry step, with the residual rate potentially accounted for by pyrophosphase release. We also observed that DNA chain-termination of the 3' end of the primer increases substantially the rate of fingers-opening in the Pol-DNA binary complex (5.3 → 29 s-1), demonstrating that the 3'-OH residue is important for the kinetics of fingers conformational changes. Our observations offer mechanistic insight and tools to offer mechanistic insight for all nucleic acid polymerases.

Controlled Trafficking of Multiple and Diverse Cations Prompts Nucleic Acid Hydrolysis

Recent in crystallo reaction intermediates have detailed how nucleic acid hydrolysis occurs in the RNA ribonuclease H1 (RNase H1), a fundamental metalloenzyme involved in maintaining the human genome. At odds with the previous characterization, these in crystallo structures unexpectedly captured multiple metal ions (K⁺ and Mg²⁺) transiently bound in the vicinity of the two-metal-ion active site. Using multi-microsecond all-atom molecular dynamics and free-energy simulations, we investigated the functional implications of the dynamic exchange of multiple K⁺ and Mg²⁺ ions at the RNase H1 reaction center. We found that such ions are timely positioned at non-overlapping locations near the active site, at different stages of catalysis, being crucial for both reactants' alignment and leaving group departure. We also found that this cation trafficking is tightly regulated by variations of the solution's ionic strength and is aided by two conserved second-shell residues, E188 and K196, suggesting a mechanism for the cations' recruitment during catalysis. These results indicate that controlled trafficking of multi-cation dynamics, opportunely prompted by second-shell residues, is functionally essential to the complex enzymatic machinery of the RNase H1. These findings revise the current knowledge on the RNase H1 catalysis and open new catalytic possibilities for other similar metalloenzymes including, but not limited to, CRISPR-Cas9, group II intron ribozyme and the human spliceosome.

Following replicative DNA synthesis by time-resolved X-ray crystallography

The mechanism of DNA synthesis has been inferred from static structures, but the absence of temporal information raises longstanding questions about the order of events in one of life's most central processes. Here we follow the reaction pathway of a replicative DNA polymerase using time-resolved X-ray crystallography to elucidate the order and transition between intermediates. In contrast to the canonical model, the structural changes observed in the time-lapsed images reveal a catalytic cycle in which translocation precedes catalysis. The translocation step appears to follow a push-pull mechanism where the O-O1 loop of the finger subdomain acts as a pawl to facilitate unidirectional movement along the template with conserved tyrosine residues 714 and 719 functioning as tandem gatekeepers of DNA synthesis. The structures capture the precise order of critical events that may be a general feature of enzymatic catalysis among replicative DNA polymerases.

DNA replication machinery: Insights from in vitro single-molecule approaches

The replisome is the multiprotein molecular machinery that replicates DNA. The replisome components work in precise coordination to unwind the double helix of the DNA and replicate the two strands simultaneously. The study of DNA replication using in vitro single-molecule approaches provides a novel quantitative understanding of the dynamics and mechanical principles that govern the operation of the replisome and its components. ‘Classical’ ensemble-averaging methods cannot obtain this information. Here we describe the main findings obtained with in vitro single-molecule methods on the performance of individual replisome components and reconstituted prokaryotic and eukaryotic replisomes. The emerging picture from these studies is that of stochastic, versatile and highly dynamic replisome machinery in which transient protein-protein and protein-DNA associations are responsible for robust DNA replication.

Programmed Allelic Mutagenesis of a DNA Polymerase with Single Amino Acid Resolution

Most DNA polymerase libraries sample unknown portions of mutational space and are constrained by the limitations of random mutagenesis. Here we describe a programmed allelic mutagenesis (PAM) strategy to comprehensively evaluate all possible single-point mutations in the entire catalytic domain of a replicative DNA polymerase. By applying the PAM strategy with ultrafast high-throughput screening, we show how DNA polymerases can be mapped for allelic mutations that exhibit enhanced activity for unnatural nucleic acid substrates. We suggest that comprehensive missense mutational scans may aid the discovery of specificity determining residues that are necessary for reprogramming the biological functions of natural DNA polymerases.

Synthesis of phosphoramidate-linked DNA by a modified DNA polymerase

Life on Earth depends on polymerases. These enzymes copy genetic information to produce the DNA and RNA strands at the core of the central dogma. Polymerases act by forming phosphodiester linkages to produce polynucleotide strands. While synthetic chemistry can generate a broad range of alternative genetic materials with unnatural linkages, polymerases have so far been limited to forming O-P bonds. Here, we show that, in fact, unnatural N-P bonds can also be formed by a modified DNA polymerase. This template-directed activity generates complementary strands linked by phosphoramidate (NP) esters, an alternative backbone linkage only known to exist in the laboratory. The emergence of NP-DNA polymerase activity implies the biochemical plausibility of alternative central dogmas for cellular life.

All known polymerases copy genetic material by catalyzing phosphodiester bond formation. This highly conserved activity proceeds by a common mechanism, such that incorporated nucleoside analogs terminate chain elongation if the resulting primer strand lacks a terminal hydroxyl group. Even conservatively substituted 3′-amino nucleotides generally act as chain terminators, and no enzymatic pathway for their polymerization has yet been found. Although 3′-amino nucleotides can be chemically coupled to yield stable oligonucleotides containing N3′→P5′ phosphoramidate (NP) bonds, no such internucleotide linkages are known to occur in nature. Here, we report that 3′-amino terminated primers are, in fact, slowly extended by the DNA polymerase from B. stearothermophilus in a template-directed manner. When its cofactor is Ca²⁺ rather than Mg²⁺, the reaction is fivefold faster, permitting multiple turnover NP bond formation to yield NP-DNA strands from the corresponding 3′-amino-2′,3′-dideoxynucleoside 5′-triphosphates. A single active site mutation further enhances the rate of NP-DNA synthesis by an additional 21-fold. We show that DNA-dependent NP-DNA polymerase activity depends on conserved active site residues and propose a likely mechanism for this activity based on a series of crystal structures of bound complexes. Our results significantly broaden the catalytic scope of polymerase activity and suggest the feasibility of a genetic transition between native nucleic acids and NP-DNA.

The Structural Basis for Processing of Unnatural Base Pairs by DNA Polymerases

Unnatural base pairs (UBPs) greatly increase the diversity of DNA and RNA, furthering their broad range of molecular biological and biotechnological approaches. Different candidates have been developed whereby alternative hydrogen-bonding patterns and hydrophobic and packing interactions have turned out to be the most promising base-pairing concepts to date. The key in many applications is the highly efficient and selective acceptance of artificial base pairs by DNA polymerases, which enables amplification of the modified DNA. In this Review, computational as well as experimental studies that were performed to characterize the pairing behavior of UBPs in free duplex DNA or bound to the active site of KlenTaq DNA polymerase are highlighted. The structural studies, on the one hand, elucidate how base pairs lacking hydrogen bonds are accepted by these enzymes and, on the other hand, highlight the influence of one or several consecutive UBPs on the structure of a DNA double helix. Understanding these concepts facilitates optimization of future UBPs for the manifold fields of applications.

Extensive free-energy simulations identify water as the base in nucleotide addition by DNA polymerase

Transphosphorylation of nucleotide triphosphates is the central reaction in DNA replication by DNA polymerase as well as many other biological processes. Despite its importance, the microscopic chemical mechanism of transphosphorylation of nucleotide triphosphates is, in most cases, unknown. Here we use extensive simulations of DNA polymerase η to test mechanistic hypotheses. We systematically survey the reactive space by calculating 2D free-energy surfaces for 10 different plausible mechanisms that have been proposed. We supplement these free-energy surfaces with calculations of pKa for a number of potentially acidic protons in different states relevant to the catalytic cycle. We find that among all of the conditions that we test, the smallest activation barrier occurs for a reaction where a Mg2+-coordinated water deprotonates the nucleophilic 3'-OH, and this deprotonation is concerted with the phosphoryl transfer. The presence of a third Mg2+ in the active site lowers the activation barrier for the water-as-base mechanism, as does protonation of the pyrophosphate leaving group, which is consistent with general acid catalysis. The results demonstrate the value of simulations, when used in conjunction with experimental data, to help establish a microscopic chemical mechanism in a complex environment.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Pyrophosphate Release in the Protein HIV Reverse Transcriptase

Enzymatic reactions usually occur in several steps: a step of substrate binding to the surface of the protein, a step of protein reorganization around the substrate and conduction of a chemical reaction, and a step of product release. The release of inorganic phosphate-PPi-from the matrix of the protein HIV reverse transcriptase is investigated computationally. Atomically detailed simulations with explicit solvent are analyzed to obtain the free energy profile, mean first passage time, and detailed molecular mechanisms of PPi escape. A challenge for the computations is of time scales. The experimental time scale of the process of interest is in milliseconds, and straightforward molecular dynamics simulations are in sub-microseconds. To overcome the time scale gap, we use the algorithm of Milestoning along a reaction coordinate to compute the overall free energy profile and rate. The methods of locally enhanced sampling and steered molecular dynamics determine plausible reaction coordinates. The observed molecular mechanism couples the transfer of the PPi to positively charged lysine side chains that are found on the exit pathway and to an exiting magnesium ion. In accord with experimental findings, the release rate is comparable to the chemical step, allowing for variations in substrate (DNA or RNA template) in which the release becomes rate determining.

Simulating the fidelity and the three Mg mechanism of pol η and clarifying the validity of transition state theory in enzyme catalysis

Pol η belongs to the important Y family of DNA polymerases that can catalyze translesion synthesis across sites of damaged DNA. This activity involves the reduced fidelity of Pol η for 8-oxo-7,8-dhyedro-2'-deoxoguanosin(8-oxoG). The fundamental interest in Pol η has grown recently with the demonstration of the importance of a 3rd Mg2+ ion. The current work explores both the fidelity of Pol η and the role of the 3rd metal ion, by using empirical valence bond (EVB) simulations. The simulations reproduce the observed trend in fidelity and shed a new light on the role of the 3rd metal ion. It is found that this ion does not lead to a major catalytic effect, but most probably plays an important role in reducing the product release barrier. Furthermore, it is concluded, in contrast to some implications, that the effect of this metal does not violate transition state theory, and the evaluation of the catalytic effect must conserve the molecular composition upon moving from the reactant to the transition state. Proteins 2017; 85:1446-1453. © 2017 Wiley Periodicals, Inc.

Structural Insights into the Post-Chemistry Steps of Nucleotide Incorporation Catalyzed by a DNA Polymerase

DNA polymerases are essential enzymes that faithfully and efficiently replicate genomic information.1-3 The mechanism of nucleotide incorporation by DNA polymerases has been extensively studied structurally and kinetically, but several key steps following phosphodiester bond formation remain structurally uncharacterized due to utilization of natural nucleotides. It is thought that the release of pyrophosphate (PP_i) triggers reverse conformational changes in a polymerase in order to complete a full catalytic cycle as well as prepare for DNA translocation and subsequent incorporation events. Here, by using the triphosphates of chain-terminating antiviral drugs lamivudine ((-)3TC-TP) and emtricitabine ((-)FTC-TP), we structurally reveal the correct sequence of post-chemistry steps during nucleotide incorporation by human DNA polymerase β (hPolβ) and provide a structural basis for PP_i release. These post-catalytic structures reveal hPolβ in an open conformation with PP_i bound in the active site, thereby strongly suggesting that the reverse conformational changes occur prior to PP_i release. The results also help to refine the role of the newly discovered third divalent metal ion for DNA polymerase-catalyzed nucleotide incorporation. Furthermore, a post-chemistry structure of hPolβ in the open conformation, following incorporation of (-)3TC-MP, with a second (-)3TC-TP molecule bound to the active site in the absence of PP_i, suggests that nucleotide binding stimulates PP_i dissociation and occurs before polymerase translocation. Our structural characterization defines the order of the elusive post-chemistry steps in the canonical mechanism of a DNA polymerase.

A Self-Activated Mechanism for Nucleic Acid Polymerization Catalyzed by DNA/RNA Polymerases

The enzymatic polymerization of DNA and RNA is the basis for genetic inheritance for all living organisms. It is catalyzed by the DNA/RNA polymerase (Pol) superfamily. Here, bioinformatics analysis reveals that the incoming nucleotide substrate always forms an H-bond between its 3'-OH and β-phosphate moieties upon formation of the Michaelis complex. This previously unrecognized H-bond implies a novel self-activated mechanism (SAM), which synergistically connects the in situ nucleophile formation with subsequent nucleotide addition and, importantly, nucleic acid translocation. Thus, SAM allows an elegant and efficient closed-loop sequence of chemical and physical steps for Pol catalysis. This is markedly different from previous mechanistic hypotheses. Our proposed mechanism is corroborated via ab initio QM/MM simulations on a specific Pol, the human DNA polymerase-η, an enzyme involved in repairing damaged DNA. The structural conservation of DNA and RNA Pols supports the possible extension of SAM to Pol enzymes from the three domains of life.

Divalent ions attenuate DNA synthesis by human DNA polymerase α by changing the structure of the template/primer or by perturbing the polymerase reaction

DNA polymerases (pols) are sophisticated protein machines operating in the replication, repair and recombination of genetic material in the complex environment of the cell. DNA pol reactions require at least two divalent metal ions for the phosphodiester bond formation. We explore two understudied roles of metals in pol transactions with emphasis on polα, a crucial enzyme in the initiation of DNA synthesis. We present evidence that the combination of many factors, including the structure of the template/primer, the identity of the metal, the metal turnover in the pol active site, and the influence of the concentration of nucleoside triphosphates, affect DNA pol synthesis. On the poly-dT70 template, the increase of Mg(2+) concentration within the range typically used for pol reactions led to the severe loss of the ability of pol to extend DNA primers and led to a decline in DNA product sizes when extending RNA primers, simulating the effect of "counting" of the number of nucleotides in nascent primers by polα. We suggest that a high Mg(2+) concentration promotes the dynamic formation of unconventional DNA structure(s), thus limiting the apparent processivity of the enzyme. Next, we found that Zn(2+) supported robust polα reactions when the concentration of nucleotides was above the concentration of ions; however, there was only one nucleotide incorporation by the Klenow fragment of DNA pol I. Zn(2+) drastically inhibited polα, but had no effect on Klenow, when Mg(2+) was also present. It is possible that Zn(2+) perturbs metal-mediated transactions in pol active site, for example affecting the step of pyrophosphate removal at the end of each pol cycle necessary for continuation of polymerization.

Cooperative motion of a key positively charged residue and metal ions for DNA replication catalyzed by human DNA Polymerase-η

Trans-lesion synthesis polymerases, like DNA Polymerase-η (Pol-η), are essential for cell survival. Pol-η bypasses ultraviolet-induced DNA damages via a two-metal-ion mechanism that assures DNA strand elongation, with formation of the leaving group pyrophosphate (PPi). Recent structural and kinetics studies have shown that Pol-η function depends on the highly flexible and conserved Arg61 and, intriguingly, on a transient third ion resolved at the catalytic site, as lately observed in other nucleic acid-processing metalloenzymes. How these conserved structural features facilitate DNA replication, however, is still poorly understood. Through extended molecular dynamics and free energy simulations, we unravel a highly cooperative and dynamic mechanism for DNA elongation and repair, which is here described by an equilibrium ensemble of structures that connect the reactants to the products in Pol-η catalysis. We reveal that specific conformations of Arg61 help facilitate the recruitment of the incoming base and favor the proper formation of a pre-reactive complex in Pol-η for efficient DNA editing. Also, we show that a third transient metal ion, which acts concertedly with Arg61, serves as an exit shuttle for the leaving PPi. Finally, we discuss how this effective and cooperative mechanism for DNA repair may be shared by other DNA-repairing polymerases.

Bacterial siderophore mimicking iron complexes as DNA targeting antimicrobials

Microbial secretion of siderophores for iron uptake can be employed as an efficient strategy to smuggle in bactericidal agents by conjugation to iron.

A Jump-from-Cavity Pyrophosphate Ion Release Assisted by a Key Lysine Residue in T7 RNA Polymerase Transcription Elongation

Pyrophosphate ion (PPi) release during transcription elongation is a signature step in each nucleotide addition cycle. The kinetics and energetics of the process as well as how it proceeds with substantial conformational changes of the polymerase complex determine the mechano-chemical coupling mechanism of the transcription elongation. Here we investigated detailed dynamics of the PPi release process in a single-subunit RNA polymerase (RNAP) from bacteriophage T7, implementing all-atom molecular dynamics (MD) simulations. We obtained a jump-from-cavity kinetic model of the PPi release utilizing extensive nanosecond MD simulations. We found that the PPi release in T7 RNAP is initiated by the PPi dissociation from two catalytic aspartic acids, followed by a comparatively slow jump-from-cavity activation process. Combining with a number of microsecond long MD simulations, we also found that the activation process is hindered by charged residue associations as well as by local steric and hydrogen bond interactions. On the other hand, the activation is greatly assisted by a highly flexible lysine residue Lys472 that swings its side chain to pull PPi out. The mechanism can apply in general to single subunit RNA and DNA polymerases with similar molecular structures and conserved key residues. Remarkably, the flexible lysine or arginine residue appears to be a universal module that assists the PPi release even in multi-subunit RNAPs with charge facilitated hopping mechanisms. We also noticed that the PPi release is not tightly coupled to opening motions of an O-helix on the fingers domain of T7 RNAP according to the microsecond MD simulations. Our study thus supports the Brownian ratchet scenario of the mechano-chemical coupling in the transcription elongation of the single-subunit polymerase.

A critical residue selectively recruits nucleotides for t7 RNA polymerase transcription fidelity control

Nucleotide selection is essential for fidelity control in gene replication and transcription. Recent work on T7 RNA polymerase suggested that a small posttranslocation free energy bias stabilizes Tyr(639) in the active site to aid nucleotide selection. However, it was not clear exactly how Tyr(639) assists the selection. Here we report a molecular-dynamics simulation study revealing atomistic detail of this critical selectivity. The study shows first that Tyr(639) blocks the active site at posttranslocation by marginally stacking to the end basepair of the DNA-RNA hybrid. The study then demonstrates that at the nucleotide preinsertion state, a cognate RNA nucleotide does not affect the local Tyr(639) stabilization, whereas a noncognate nucleotide substantially stabilizes Tyr(639) so that Tyr(639) keeps blocking the active site. As a result, further nucleotide insertion into the active site, which requires moving Tyr(639) out of the site, would be hindered for the noncognate nucleotide, but not for the cognate nucleotide. In particular, we note that water molecules assist the ribose recognition in the RNA nucleotide preinsertion, and help Tyr(639) stacking to the end basepair in the case of a DNA nucleotide. It was also seen that a base-mismatched nucleotide at preinsertion directly grabs Tyr(639) for the active site stabilization. We also find that in a mutant polymerase Y639F the strong stabilization of residue 639 in the active site cannot establish upon the DNA nucleotide preinsertion. The finding explains the reduced differentiation between ribo- and deoxyribonucleotides that has been recorded experimentally for the mutant polymerase.

The Closing Mechanism of DNA Polymerase I at Atomic Resolution

DNA polymerases must quickly and accurately distinguish between similar nucleic acids to form Watson-Crick base pairs and avoid DNA replication errors. Deoxynucleoside triphosphate (dNTP) binding to the DNA polymerase active site induces a large conformational change that is difficult to characterize experimentally on an atomic level. Here, we report an X-ray crystal structure of DNA polymerase I bound to DNA in the open conformation with a dNTP present in the active site. We use this structure to computationally simulate the open to closed transition of DNA polymerase in the presence of a Watson-Crick base pair. Our microsecond simulations allowed us to characterize the key steps involved in active site assembly, and propose the sequence of events involved in the prechemistry steps of DNA polymerase catalysis. They also reveal new features of the polymerase mechanism, such as a conserved histidine as a potential proton acceptor from the primer 3'-hydroxyl.

Mechano-chemical kinetics of DNA replication: identification of the translocation step of a replicative DNA polymerase

During DNA replication replicative polymerases move in discrete mechanical steps along the DNA template. To address how the chemical cycle is coupled to mechanical motion of the enzyme, here we use optical tweezers to study the translocation mechanism of individual bacteriophage Phi29 DNA polymerases during processive DNA replication. We determine the main kinetic parameters of the nucleotide incorporation cycle and their dependence on external load and nucleotide (dNTP) concentration. The data is inconsistent with power stroke models for translocation, instead supports a loose-coupling mechanism between chemical catalysis and mechanical translocation during DNA replication. According to this mechanism the DNA polymerase works by alternating between a dNTP/PPi-free state, which diffuses thermally between pre- and post-translocated states, and a dNTP/PPi-bound state where dNTP binding stabilizes the post-translocated state. We show how this thermal ratchet mechanism is used by the polymerase to generate work against large opposing loads (∼50 pN).

Molecular dynamics study of the opening mechanism for DNA polymerase I

During DNA replication, DNA polymerases follow an induced fit mechanism in order to rapidly distinguish between correct and incorrect dNTP substrates. The dynamics of this process are crucial to the overall effectiveness of catalysis. Although X-ray crystal structures of DNA polymerase I with substrate dNTPs have revealed key structural states along the catalytic pathway, solution fluorescence studies indicate that those key states are populated in the absence of substrate. Herein, we report the first atomistic simulations showing the conformational changes between the closed, open, and ajar conformations of DNA polymerase I in the binary (enzyme:DNA) state to better understand its dynamics. We have applied long time-scale, unbiased molecular dynamics to investigate the opening process of the fingers domain in the absence of substrate for B. stearothermophilis DNA polymerase in silico. These simulations are biologically and/or physiologically relevant as they shed light on the transitions between states in this important enzyme. All closed and ajar simulations successfully transitioned into the fully open conformation, which is known to be the dominant binary enzyme-DNA conformation from solution and crystallographic studies. Furthermore, we have detailed the key stages in the opening process starting from the open and ajar crystal structures, including the observation of a previously unknown key intermediate structure. Four backbone dihedrals were identified as important during the opening process, and their movements provide insight into the recognition of dNTP substrate molecules by the polymerase binary state. In addition to revealing the opening mechanism, this study also demonstrates our ability to study biological events of DNA polymerase using current computational methods without biasing the dynamics.

Investigations of α-helix↔β-sheet transition pathways in a miniprotein using the finite-temperature string method

A parallel implementation of the finite-temperature string method is described, which takes into account the invariance of coordinates with respect to rigid-body motions. The method is applied to the complex α-helix↔β-sheet transition in a β-hairpin miniprotein in implicit solvent, which exhibits much of the complexity of conformational changes in proteins. Two transition paths are considered, one derived from a linear interpolant between the endpoint structures and the other derived from a targeted dynamics simulation. Two methods for computing the conformational free energy (FE) along the string are compared, a restrained method, and a tessellation method introduced by E. Vanden-Eijnden and M. Venturoli [J. Chem. Phys. 130, 194103 (2009)]. It is found that obtaining meaningful free energy profiles using the present atom-based coordinates requires restricting sampling to a vicinity of the converged path, where the hyperplanar approximation to the isocommittor surface is sufficiently accurate. This sampling restriction can be easily achieved using restraints or constraints. The endpoint FE differences computed from the FE profiles are validated by comparison with previous calculations using a path-independent confinement method. The FE profiles are decomposed into the enthalpic and entropic contributions, and it is shown that the entropy difference contribution can be as large as 10 kcal/mol for intermediate regions along the path, compared to 15-20 kcal/mol for the enthalpy contribution. This result demonstrates that enthalpic barriers for transitions are offset by entropic contributions arising from the existence of different paths across a barrier. The possibility of using systematically coarse-grained representations of amino acids, in the spirit of multiple interaction site residue models, is proposed as a means to avoid ad hoc sampling restrictions to narrow transition tubes.

Exploring transition pathway and free-energy profile of large-scale protein conformational change by combining normal mode analysis and umbrella sampling molecular dynamics

Large-scale conformational changes of proteins are usually associated with the binding of ligands. Because the conformational changes are often related to the biological functions of proteins, understanding the molecular mechanisms of these motions and the effects of ligand binding becomes very necessary. In the present study, we use the combination of normal-mode analysis and umbrella sampling molecular dynamics simulation to delineate the atomically detailed conformational transition pathways and the associated free-energy landscapes for three well-known protein systems, viz., adenylate kinase (AdK), calmodulin (CaM), and p38α kinase in the absence and presence of respective ligands. For each protein under study, the transient conformations along the conformational transition pathway and thermodynamic observables are in agreement with experimentally and computationally determined ones. The calculated free-energy profiles reveal that AdK and CaM are intrinsically flexible in structures without obvious energy barrier, and their ligand binding shifts the equilibrium from the ligand-free to ligand-bound conformation (population shift mechanism). In contrast, the ligand binding to p38α leads to a large change in free-energy barrier (ΔΔG ≈ 7 kcal/mol), promoting the transition from DFG-in to DFG-out conformation (induced fit mechanism). Moreover, the effect of the protonation of D168 on the conformational change of p38α is also studied, which reduces the free-energy difference between the two functional states of p38α and thus further facilitates the conformational interconversion. Therefore, the present study suggests that the detailed mechanism of ligand binding and the associated conformational transition is not uniform for all kinds of proteins but correlated to their respective biological functions.

Structural insights into DNA replication without hydrogen bonds

The genetic alphabet is composed of two base pairs, and the development of a third, unnatural base pair would increase the genetic and chemical potential of DNA. d5SICS-dNaM is one of the most efficiently replicated unnatural base pairs identified to date, but its pairing is mediated by only hydrophobic and packing forces, and in free duplex DNA it forms a cross-strand intercalated structure that makes its efficient replication difficult to understand. Recent studies of the KlenTaq DNA polymerase revealed that the insertion of d5SICSTP opposite dNaM proceeds via a mutually induced-fit mechanism, where the presence of the triphosphate induces the polymerase to form the catalytically competent closed structure, which in turn induces the pairing nucleotides of the developing unnatural base pair to adopt a planar Watson-Crick-like structure. To understand the remaining steps of replication, we now report the characterization of the prechemistry complexes corresponding to the insertion of dNaMTP opposite d5SICS, as well as multiple postchemistry complexes in which the already formed unnatural base pair is positioned at the postinsertion site. Unlike with the insertion of d5SICSTP opposite dNaM, addition of dNaMTP does not fully induce the formation of the catalytically competent closed state. The data also reveal that once synthesized and translocated to the postinsertion position, the unnatural nucleobases again intercalate. Two modes of intercalation are observed, depending on the nature of the flanking nucleotides, and are each stabilized by different interactions with the polymerase, and each appear to reduce the affinity with which the next correct triphosphate binds. Thus, continued primer extension is limited by deintercalation and rearrangements with the polymerase active site that are required to populate the catalytically active, triphosphate bound conformation.

Cooperative dynamics of a DNA polymerase replicating complex

Engineered DNA polymerases continue to be the workhorses of many applications in biotechnology, medicine and nanotechnology. However, the dynamic interplay between the enzyme and the DNA remains unclear. In this study, we performed an extensive replica exchange with flexible tempering (REFT) molecular dynamics simulation of the ternary replicating complex of the archaeal family B DNA polymerase from the thermophile Thermococcus gorgonarius, right before the chemical step. The convoluted dynamics of the enzyme are reducible to rigid-body motions of six subdomains. Upon binding to the enzyme, the DNA double helix conformation changes from a twisted state to a partially untwisted state. The twisted state displays strong bending motion, whereby the DNA oscillates between a straight and a bent conformation. The dynamics of double-stranded DNA are strongly correlated with rotations of the thumb toward the palm, which suggests an assisting role of the enzyme during DNA translocation. In the complex, the primer-template duplex displays increased preference for the B-DNA conformation at the n-2 and n-3 dinucleotide steps. Interactions at the primer 3' end indicate that Thr541 and Asp540 are the acceptors of the first proton transfer in the chemical step, whereas in the translocation step both residues hold the primer 3' terminus in the vicinity of the priming site, which is crucial for high processivity.

Structures of intermediates along the catalytic cycle of terminal deoxynucleotidyltransferase: dynamical aspects of the two-metal ion mechanism

Terminal deoxynucleotidyltransferase (Tdt) is a non-templated eukaryotic DNA polymerase of the polX family that is responsible for the random addition of nucleotides at the V(D)J junctions of immunoglobulins and T-cell receptors. Here we describe a series of high-resolution X-ray structures that mimic the pre-catalytic state, the post-catalytic state and a competent state that can be transformed into the two other ones in crystallo via the addition of dAMPcPP and Zn(2+), respectively. We examined the effect of Mn(2+), Co(2+) and Zn(2+) because they all have a marked influence on the kinetics of the reaction. We demonstrate a dynamic role of divalent transition metal ions bound to site A: (i) Zn(2+) (or Co(2+)) in Metal A site changes coordination from octahedral to tetrahedral after the chemical step, which explains the known higher affinity of Tdt for the primer strand when these ions are present, and (ii) metal A has to leave to allow the translocation of the primer strand and to clear the active site, a typical feature for a ratchet-like mechanism. Except for Zn(2+), the sugar puckering of the primer strand 3' terminus changes from C2'-endo to C3'-endo during catalysis. In addition, our data are compatible with a scheme where metal A is the last component that binds to the active site to complete its productive assembly, as already inferred in human pol beta. The new structures have potential implications for modeling pol mu, a closely related polX implicated in the repair of DNA double-strand breaks, in a complex with a DNA synapsis.

Kinetic mechanism of translocation and dNTP binding in individual DNA polymerase complexes

Complexes formed between phi29 DNA polymerase (DNAP) and DNA fluctuate discretely between the pre-translocation and post-translocation states on the millisecond time scale. The translocation fluctuations can be observed in ionic current traces when individual complexes are captured atop the α-hemolysin nanopore in an electric field. The presence of complementary 2'-deoxynucleoside triphosphate (dNTP) shifts the equilibrium across the translocation step toward the post-translocation state. Here we have determined quantitatively the kinetic relationship between the phi29 DNAP translocation step and dNTP binding. We demonstrate that dNTP binds to phi29 DNAP-DNA complexes only after the transition from the pre-translocation state to the post-translocation state; dNTP binding rectifies the translocation but it does not directly drive the translocation. Based on the measured time traces of current amplitude, we developed a method for determining the forward and reverse translocation rates and the dNTP association and dissociation rates, individually at each dNTP concentration and each voltage. The translocation rates, and their response to force, match those determined for phi29 DNAP-DNA binary complexes and are unaffected by dNTP. The dNTP association and dissociation rates do not vary as a function of voltage, indicating that force does not distort the polymerase active site and that dNTP binding does not directly involve a displacement in the translocation direction. This combined experimental and theoretical approach and the results obtained provide a framework for separately evaluating the effects of biological variables on the translocation transitions and their effects on dNTP binding.

Structures of coxsackievirus, rhinovirus, and poliovirus polymerase elongation complexes solved by engineering RNA mediated crystal contacts

RNA-dependent RNA polymerases play a vital role in the growth of RNA viruses where they are responsible for genome replication, but do so with rather low fidelity that allows for the rapid adaptation to different host cell environments. These polymerases are also a target for antiviral drug development. However, both drug discovery efforts and our understanding of fidelity determinants have been hampered by a lack of detailed structural information about functional polymerase-RNA complexes and the structural changes that take place during the elongation cycle. Many of the molecular details associated with nucleotide selection and catalysis were revealed in our recent structure of the poliovirus polymerase-RNA complex solved by first purifying and then crystallizing stalled elongation complexes. In the work presented here we extend that basic methodology to determine nine new structures of poliovirus, coxsackievirus, and rhinovirus elongation complexes at 2.2-2.9 Å resolution. The structures highlight conserved features of picornaviral polymerases and the interactions they make with the template and product RNA strands, including a tight grip on eight basepairs of the nascent duplex, a fully pre-positioned templating nucleotide, and a conserved binding pocket for the +2 position template strand base. At the active site we see a pre-bound magnesium ion and there is conservation of a non-standard backbone conformation of the template strand in an interaction that may aid in triggering RNA translocation via contact with the conserved polymerase motif B. Moreover, by engineering plasticity into RNA-RNA contacts, we obtain crystal forms that are capable of multiple rounds of in-crystal catalysis and RNA translocation. Together, the data demonstrate that engineering flexible RNA contacts to promote crystal lattice formation is a versatile platform that can be used to solve the structures of viral RdRP elongation complexes and their catalytic cycle intermediates.

Molecular machines

Molecular machines (MM) are essential components of living cells. They conduct mechanical work, transport materials into and out of cells, assist in processing enzymatic reactions, and more. Their operations are frequently combined with significant conformational transitions. Computational studies of these conformational transitions and their coupling to molecular functions are discussed. It is argued that coarse descriptions of these molecules which are based on mass density and shape provide useful information on directions of action. It is further argued that MM are likely to have well focused and narrow reaction pathways. The proposal for such pathways is supported by evolutionary analyses of homologous machines. Finally, these observations are used to build atomically detailed models of these systems that are making the link from structure to functions (kinetics and thermodynamics). For that purpose enhanced sampling techniques are required.

Dynamics of the translocation step measured in individual DNA polymerase complexes

Complexes formed between the bacteriophage phi29 DNA polymerase (DNAP) and DNA fluctuate between the pre-translocation and post-translocation states on the millisecond time scale. These fluctuations can be directly observed with single-nucleotide precision in real-time ionic current traces when individual complexes are captured atop the α-hemolysin nanopore in an applied electric field. We recently quantified the equilibrium across the translocation step as a function of applied force (voltage), active-site proximal DNA sequences, and the binding of complementary dNTP. To gain insight into the mechanism of this step in the DNAP catalytic cycle, in this study, we have examined the stochastic dynamics of the translocation step. The survival probability of complexes in each of the two states decayed at a single exponential rate, indicating that the observed fluctuations are between two discrete states. We used a robust mathematical formulation based on the autocorrelation function to extract the forward and reverse rates of the transitions between the pre-translocation state and the post-translocation state from ionic current traces of captured phi29 DNAP-DNA binary complexes. We evaluated each transition rate as a function of applied voltage to examine the energy landscape of the phi29 DNAP translocation step. The analysis reveals that active-site proximal DNA sequences influence the depth of the pre-translocation and post-translocation state energy wells and affect the location of the transition state along the direction of the translocation.

Analysis and elimination of a bias in targeted molecular dynamics simulations of conformational transitions: application to calmodulin

The popular targeted molecular dynamics (TMD) method for generating transition paths in complex biomolecular systems is revisited. In a typical TMD transition path, the large-scale changes occur early and the small-scale changes tend to occur later. As a result, the order of events in the computed paths depends on the direction in which the simulations are performed. To identify the origin of this bias, and to propose a method in which the bias is absent, variants of TMD in the restraint formulation are introduced and applied to the complex open ↔ closed transition in the protein calmodulin. Due to the global best-fit rotation that is typically part of the TMD method, the simulated system is guided implicitly along the lowest-frequency normal modes, until the large spatial scales associated with these modes are near the target conformation. The remaining portion of the transition is described progressively by higher-frequency modes, which correspond to smaller-scale rearrangements. A straightforward modification of TMD that avoids the global best-fit rotation is the locally restrained TMD (LRTMD) method, in which the biasing potential is constructed from a number of TMD potentials, each acting on a small connected portion of the protein sequence. With a uniform distribution of these elements, transition paths that lack the length-scale bias are obtained. Trajectories generated by steered MD in dihedral angle space (DSMD), a method that avoids best-fit rotations altogether, also lack the length-scale bias. To examine the importance of the paths generated by TMD, LRTMD, and DSMD in the actual transition, we use the finite-temperature string method to compute the free energy profile associated with a transition tube around a path generated by each algorithm. The free energy barriers associated with the paths are comparable, suggesting that transitions can occur along each route with similar probabilities. This result indicates that a broad ensemble of paths needs to be calculated to obtain a full description of conformational changes in biomolecules. The breadth of the contributing ensemble suggests that energetic barriers for conformational transitions in proteins are offset by entropic contributions that arise from a large number of possible paths.

Active site opening and closure control translocation of multisubunit RNA polymerase

Multisubunit RNA polymerase (RNAP) is the central information-processing enzyme in all cellular life forms, yet its mechanism of translocation along the DNA molecule remains conjectural. Here, we report direct monitoring of bacterial RNAP translocation following the addition of a single nucleotide. Time-resolved measurements demonstrated that translocation is delayed relative to nucleotide incorporation and occurs shortly after or concurrently with pyrophosphate release. An investigation of translocation equilibrium suggested that the strength of interactions between RNA 3′ nucleotide and nucleophilic and substrate sites determines the translocation state of transcription elongation complexes, whereas active site opening and closure modulate the affinity of the substrate site, thereby favoring the post- and pre-translocated states, respectively. The RNAP translocation mechanism is exploited by the antibiotic tagetitoxin, which mimics pyrophosphate and induces backward translocation by closing the active site.

Trigger loop dynamics mediate the balance between the transcriptional fidelity and speed of RNA polymerase II

During transcription, RNA polymerase II (RNAPII) must select the correct nucleotide, catalyze its addition to the growing RNA transcript, and move stepwise along the DNA until a gene is fully transcribed. In all kingdoms of life, transcription must be finely tuned to ensure an appropriate balance between fidelity and speed. Here, we used an optical-trapping assay with high spatiotemporal resolution to probe directly the motion of individual RNAPII molecules as they pass through each of the enzymatic steps of transcript elongation. We report direct evidence that the RNAPII trigger loop, an evolutionarily conserved protein subdomain, serves as a master regulator of transcription, affecting each of the three main phases of elongation, namely: substrate selection, translocation, and catalysis. Global fits to the force-velocity relationships of RNAPII and its trigger loop mutants support a Brownian ratchet model for elongation, where the incoming NTP is able to bind in either the pre- or posttranslocated state, and movement between these two states is governed by the trigger loop. Comparison of the kinetics of pausing by WT and mutant RNAPII under conditions that promote base misincorporation indicate that the trigger loop governs fidelity in substrate selection and mismatch recognition, and thereby controls aspects of both transcriptional accuracy and rate.

A small post-translocation energy bias aids nucleotide selection in T7 RNA polymerase transcription

The RNA polymerase (RNAP) of bacteriophage T7 is a single subunit enzyme that can transcribe DNA to RNA in the absence of additional protein factors. In this work, we present a model of T7 RNAP translocation during elongation. Based on structural information and experimental data from single-molecule force measurements, we show that a small component of facilitated translocation or power stroke coexists with the Brownian-ratchet-driven motions, and plays a crucial role in nucleotide selection at pre-insertion. The facilitated translocation is carried out by the conserved Tyr(639) that moves its side chain into the active site, pushing aside the 3'-end of the RNA, and forming a locally stabilized post-translocation intermediate. Pre-insertion of an incoming nucleotide into this stabilized intermediate state ensures that Tyr(639) closely participates in selecting correct nucleotides. A similar translocation mechanism has been suggested for multi-subunit RNAPs involving the bridge-helix bending. Nevertheless, the bent bridge-helix sterically prohibits nucleotide binding in the post-transolocation intermediate analog; moreover, the analog is not stabilized unless an inhibitory protein factor binds to the enzyme. Using our scheme, we also compared the efficiencies of different strategies for nucleotide selection, and examined effects of facilitated translocation on forward tracking.

Dynamics of pyrophosphate ion release and its coupled trigger loop motion from closed to open state in RNA polymerase II

Pyrophosphate ion (PP(i)) release after nucleotide incorporation is a necessary step for RNA polymerase II (pol II) to enter the next nucleotide addition cycle during transcription elongation. However, the role of pol II residues in PP(i) release and the mechanistic relationship between PP(i) release and the conformational change of the trigger loop remain unclear. In this study, we constructed a Markov state model (MSM) from extensive all-atom molecular dynamics (MD) simulations in the explicit solvent to simulate the PP(i) release process along the pol II secondary channel. Our results show that the trigger loop has significantly larger intrinsic motion after catalysis and formation of PP(i), which in turn aids PP(i) release mainly through the hydrogen bonding between the trigger loop residue H1085 and the (Mg-PP(i))(2-) group. Once PP(i) leaves the active site, it adopts a hopping model through several highly conserved positively charged residues such as K752 and K619 to release from the pol II pore region of the secondary channel. These positive hopping sites form favorable interactions with PP(i) and generate four kinetically metastable states as identified by our MSM. Furthermore, our single-mutant simulations suggest that H1085 and K752 aid PP(i) exit from the active site after catalysis, whereas K619 facilitates its passage through the secondary channel. Finally, we suggest that PP(i) release could help the opening motion of the trigger loop, even though PP(i) release precedes full opening of the trigger loop due to faster PP(i) dynamics. Our simulations provide predictions to guide future experimental tests.

The structural basis of the kinetic mechanism of a gap-filling X-family DNA polymerase that binds Mg(2+)-dNTP before binding to DNA

DNA with single-nucleotide (1-nt) gaps can arise during various DNA processing events. These lesions are repaired by X-family DNA polymerases (PolXs) with high gap-filling activity. Some PolXs can bind productively to dNTPs in the absence of DNA and fill these 1-nt gaps. Although PolXs have a crucial role in efficient gap filling, currently, little is known of the kinetic and structural details of their productive dNTP binding. Here, we show that Thermus thermophilus HB8 PolX (ttPolX) had strong binding affinity for Mg(2+)-dNTPs in the absence of DNA and that it follows a Theorell-Chance (hit-and-run) mechanism with nucleotide binding first. Comparison of the intermediate crystal structures of ttPolX in a binary complex with dGTP and in a ternary complex with 1-nt gapped DNA and Mg(2+)-ddGTP revealed that the conformation of the incoming nucleotide depended on whether or not DNA was present. Furthermore, the Lys263 residue located between two guanosine conformations was essential to the strong binding affinity of the enzyme. The ability to bind to either syn-dNTP or anti-dNTP and the involvement of a Theorell-Chance mechanism are key aspects of the strong nucleotide-binding and efficient gap-filling activities of ttPolX.

Structural evidence for the rare tautomer hypothesis of spontaneous mutagenesis

Even though high-fidelity polymerases copy DNA with remarkable accuracy, some base-pair mismatches are incorporated at low frequency, leading to spontaneous mutagenesis. Using high-resolution X-ray crystallographic analysis of a DNA polymerase that catalyzes replication in crystals, we observe that a C • A mismatch can mimic the shape of cognate base pairs at the site of incorporation. This shape mimicry enables the mismatch to evade the error detection mechanisms of the polymerase, which would normally either prevent mismatch incorporation or promote its nucleolytic excision. Movement of a single proton on one of the mismatched bases alters the hydrogen-bonding pattern such that a base pair forms with an overall shape that is virtually indistinguishable from a canonical, Watson-Crick base pair in double-stranded DNA. These observations provide structural evidence for the rare tautomer hypothesis of spontaneous mutagenesis, a long-standing concept that has been difficult to demonstrate directly.

Biomolecularmodeling and simulation: a field coming of age

We assess the progress in biomolecular modeling and simulation, focusing on structure prediction and dynamics, by presenting the field’s history, metrics for its rise in popularity, early expressed expectations, and current significant applications. The increases in computational power combined with improvements in algorithms and force fields have led to considerable success, especially in protein folding, specificity of ligand/biomolecule interactions, and interpretation of complex experimental phenomena (e.g. NMR relaxation, protein-folding kinetics and multiple conformational states) through the generation of structural hypotheses and pathway mechanisms. Although far from a general automated tool, structure prediction is notable for proteins and RNA that preceded the experiment, especially by knowledge-based approaches. Thus, despite early unrealistic expectations and the realization that computer technology alone will not quickly bridge the gap between experimental and theoretical time frames, ongoing improvements to enhance the accuracy and scope of modeling and simulation are propelling the field onto a productive trajectory to become full partner with experiment and a field on its own right.

The structure of a high fidelity DNA polymerase bound to a mismatched nucleotide reveals an "ajar" intermediate conformation in the nucleotide selection mechanism

To achieve accurate DNA synthesis, DNA polymerases must rapidly sample and discriminate against incorrect nucleotides. Here we report the crystal structure of a high fidelity DNA polymerase I bound to DNA primer-template caught in the act of binding a mismatched (dG:dTTP) nucleoside triphosphate. The polymerase adopts a conformation in between the previously established "open" and "closed" states. In this "ajar" conformation, the template base has moved into the insertion site but misaligns an incorrect nucleotide relative to the primer terminus. The displacement of a conserved active site tyrosine in the insertion site by the template base is accommodated by a distinctive kink in the polymerase O helix, resulting in a partially open ternary complex. We suggest that the ajar conformation allows the template to probe incoming nucleotides for complementarity before closure of the enzyme around the substrate. Based on solution fluorescence, kinetics, and crystallographic analyses of wild-type and mutant polymerases reported here, we present a three-state reaction pathway in which nucleotides either pass through this intermediate conformation to the closed conformation and catalysis or are misaligned within the intermediate, leading to destabilization of the closed conformation.

In silico study of the inhibition of DNA polymerase by a novel catalpol derivative

In this work, a novel catalpol derivative (6,10,2',6'-tetraacetyl-O-catalpol), which was previously obtained by our group and shown experimentally to inhibit a type of Taq DNA polymerase, was studied in silico. Studies of the interaction of 6,10,2',6'-tetraacetyl-O-catalpol with the Klentaq fragment of the Taq DNA polymerase I from Thermus aquaticus helped to elucidate the mechanism of inhibition of the enzyme, and offered valuable information that can be used to propose substrate structural modifications aimed at increasing the binding affinity. Classical and semi-empirical methods were used to characterize the conformational preferences of this organic compound in solution. Using docking simulations, the most probable binding mode was found, and the stabilities of the docked solutions were tested in a series of molecular dynamics experiments. Results indicated that the mechanism of inhibition may be competitive, which agrees with previous binding experiments done with 6,10,2',6'-tetraacetyl-O-catalpol.