Orchestration of cooperative events in DNA synthesis and repair mechanism unraveled by transition path sampling of DNA polymerase beta's closing

In silico studies of the African swine fever virus DNA polymerase X support an induced-fit mechanism

The African swine fever virus DNA polymerase X (pol X), a member of the X family of DNA polymerases, is thought to be involved in base excision repair. Kinetics data indicate that pol X catalyzes DNA polymerization with low fidelity, suggesting a role in viral mutagenesis. Though pol X lacks the fingers domain that binds the DNA in other members of the X family, it binds DNA tightly. To help interpret details of this interaction, molecular dynamics simulations of free pol X at different salt concentrations and of pol X bound to gapped DNA, in the presence and in the absence of the incoming nucleotide, are performed. Anchors for the simulations are two NMR structures of pol X without DNA and a model of one NMR structure plus DNA and incoming nucleotide. Our results show that, in its free form, pol X can exist in two stable conformations that interconvert to one another depending on the salt concentration. When gapped double stranded DNA is introduced near the active site, pol X prefers an open conformation, regardless of the salt concentration. Finally, under physiological conditions, in the presence of both gapped DNA and correct incoming nucleotide, and two divalent ions, the thumb subdomain of pol X undergoes a large conformational change, closing upon the DNA. These results predict for pol X a substrate-induced conformational change triggered by the presence of DNA and the correct incoming nucleotide in the active site, as in DNA polymerase beta. The simulations also suggest specific experiments (e.g., for mutants Phe-102Ala, Val-120Gly, and Lys-85Val that may reveal crucial DNA binding and active-site organization roles) to further elucidate the fidelity mechanism of pol X.

Fidelity discrimination in DNA polymerase beta: differing closing profiles for a mismatched (G:A) versus matched (G:C) base pair

Understanding fidelity-the faithful replication or repair of DNA by polymerases-requires tracking of the structural and energetic changes involved, including the elusive transient intermediates, for nucleotide incorporation at the template/primer DNA junction. We report, using path sampling simulations and a reaction network model, strikingly different transition states in DNA polymerase beta's conformational closing for correct dCTP versus incorrect dATP incoming nucleotide opposite a template G. The cascade of transition states leads to differing active-site assembly processes toward the "two-metal-ion catalysis" geometry. We demonstrate that these context-specific pathways imply different selection processes: while active-site assembly occurs more rapidly with the correct nucleotide and leads to primer extension, the enzyme remains open longer, has a more transient closed state, and forms product more slowly when an incorrect nucleotide is present. Our results also suggest that the rate-limiting step in pol beta's conformational closing is not identical to that for overall nucleotide insertion and that the rate-limiting step in the overall nucleotide incorporation process for matched as well as mismatched systems occurs after the closing conformational change.

Computer simulations of protein functions: searching for the molecular origin of the replication fidelity of DNA polymerases

The use of computers to simulate the functions of complex biological macromolecules is essential to achieve a microscopic description of biological processes and to model and interpret experimental data. Here we apply theoretical computational approaches to investigate the fidelity of T7 DNA polymerase, divided into discrete steps that include contributions from substrate binding, pK(a) shifts, and rate constants for the PO bond-breaking and bond-making processes. We begin by defining the discrimination between right and wrong nucleotides in terms of the free energy landscape for the dNMP incorporation reaction. We then use the linear response approximation and the empirical valence bond methods to obtain converging results for the contribution of the binding and chemical steps to the overall fidelity. These approaches are successful in reproducing general trends in the observed polymerase incorporation fidelity. The calculations demonstrate the potential for further integration of theoretical and experimental studies to analyze high- and low-fidelity DNA polymerases.

Molecular dynamics and protein function

A fundamental appreciation for how biological macromolecules work requires knowledge of structure and dynamics. Molecular dynamics simulations provide powerful tools for the exploration of the conformational energy landscape accessible to these molecules, and the rapid increase in computational power coupled with improvements in methodology makes this an exciting time for the application of simulation to structural biology. In this Perspective we survey two areas, protein folding and enzymatic catalysis, in which simulations have contributed to a general understanding of mechanism. We also describe results for the F(1) ATPase molecular motor and the Src family of signaling proteins as examples of applications of simulations to specific biological systems.

Structural insights into DNA polymerase beta deterrents for misincorporation support an induced-fit mechanism for fidelity

DNA polymerases generally select the correct nucleotide from a pool of structurally similar molecules to preserve Watson-Crick base-pairing rules. We report the structure of DNA polymerase beta with DNA mismatches situated in the polymerase active site. This was achieved by using nicked product DNA that traps the mispair (template-primer, A-C or T-C) in the nascent base pair binding pocket. The structure of each mispair complex indicates that the bases do not form hydrogen bonds with one another, but form a staggered arrangement where the bases of the mispair partially overlap. This prevents closure/opening of the N subdomain that is believed to be required for catalytic cycling. The partially open conformation of the N subdomain results in distinct hydrogen bonding networks that are unique for each mispair. These structures define diverse molecular aspects of misinsertion that are consistent with the induced-fit model for substrate specificity.

Conformational Transition Pathway of Polymerase β/DNA upon Binding Correct Incoming Substrate

The closing conformational transition of wild-type polymerase beta bound to DNA template/primer before the chemical step (nucleotidyl transfer reaction) is simulated using the stochastic difference equation (in length version, "SDEL") algorithm that approximates long-time dynamics. The order of the events and the intermediate states during pol beta's closing pathway are identified and compared to a separate study of pol beta using transition path sampling (TPS) (Radhakrishnan, R.; Schlick, T. Proc. Natl. Acad. Sci. USA 2004, 101, 5970-5975). Results highlight the cooperative and subtle conformational changes in the pol beta active site upon binding the correct substrate that may help explain DNA replication and repair fidelity. These changes involve key residues that differentiate the open from the closed conformation (Asp192, Arg258, Phe272), as well as residues contacting the DNA template/primer strand near the active site (Tyr271, Arg283, Thr292, Tyr296) and residues contacting the beta and gamma phosphates of the incoming nucleotide (Ser180, Arg183, Gly189). This study compliments experimental observations by providing detailed atomistic views of the intermediates along the polymerase closing pathway and by suggesting additional key residues that regulate events prior to or during the chemical reaction. We also show general agreement between two sampling methods (the stochastic difference equation and transition path sampling) and identify methodological challenges involved in the former method relevant to large-scale biomolecular applications. Specifically, SDEL is very quick relative to TPS for obtaining an approximate path of medium resolution and providing qualitative information on the sequence of events; however, associated free energies are likely very costly to obtain because this will require both successful further refinement of the path segments close to the bottlenecks and large computational time.

Highly organized but pliant active site of DNA polymerase beta: compensatory mechanisms in mutant enzymes revealed by dynamics simulations and energy analyses

To link conformational transitions noted for DNA polymerases with kinetic results describing catalytic efficiency and fidelity, we investigate the role of key DNA polymerase beta residues on subdomain motion through simulations of five single-residue mutants: Arg-283-Ala, Tyr-271-Ala, Asp-276-Val, Arg-258-Lys, and Arg-258-Ala. Since a movement toward a closed state was only observed for R258A, we suggest that Arg(258) is crucial in modulating motion preceding chemistry. Analyses of protein/DNA interactions in the mutant active site indicate distinctive hydrogen bonding and van der Waals patterns arising from compensatory structural adjustments. By comparing closed mutant complexes with the wild-type enzyme, we interpret experimentally derived nucleotide binding affinities in molecular terms: R283A (decreased), Y271A (increased), D276V (increased), and R258A (decreased). Thus, compensatory interactions (e.g., in Y271A with adjacent residues Phe(272), Asn(279), and Arg(283)) increase the overall binding affinity for the incoming nucleotide although direct interactions may decrease. Together with energetic analyses, we predict that R258G might increase the rate of nucleotide insertion and maintain enzyme fidelity as R258A; D276L might increase the nucleotide binding affinity more than D276V; and R283A/K280A might decrease the nucleotide binding affinity and increase misinsertion more than R283A. The combined observations regarding key roles of specific residues (e.g., Arg(258)) and compensatory interactions echo the dual nature of polymerase active site, namely versatility (to accommodate various basepairs) and specificity (for preserving fidelity) and underscore an organized but pliant active site essential to enzyme function.

Exploring the helix-coil transition via all-atom equilibrium ensemble simulations

The ensemble folding of two 21-residue alpha-helical peptides has been studied using all-atom simulations under several variants of the AMBER potential in explicit solvent using a global distributed computing network. Our extensive sampling, orders of magnitude greater than the experimental folding time, results in complete convergence to ensemble equilibrium. This allows for a quantitative assessment of these potentials, including a new variant of the AMBER-99 force field, denoted AMBER-99 phi, which shows improved agreement with experimental kinetic and thermodynamic measurements. From bulk analysis of the simulated AMBER-99 phi equilibrium, we find that the folding landscape is pseudo-two-state, with complexity arising from the broad, shallow character of the "native" and "unfolded" regions of the phase space. Each of these macrostates allows for configurational diffusion among a diverse ensemble of conformational microstates with greatly varying helical content and molecular size. Indeed, the observed structural dynamics are better represented as a conformational diffusion than as a simple exponential process, and equilibrium transition rates spanning several orders of magnitude are reported. After multiple nucleation steps, on average, helix formation proceeds via a kinetic "alignment" phase in which two or more short, low-entropy helical segments form a more ideal, single-helix structure.

In silico evidence for DNA polymerase-beta's substrate-induced conformational change

Structural information for mammalian DNA pol-beta combined with molecular and essential dynamics studies have provided atomistically detailed views of functionally important conformational rearrangements that occur during DNA repair and replication. This conformational closing before the chemical reaction is explored in this work as a function of the bound substrate. Anchors for our study are available in crystallographic structures of the DNA pol-beta in "open" (polymerase bound to gapped DNA) and "closed" (polymerase bound to gapped DNA and substrate, dCTP) forms; these different states have long been used to deduce that a large-scale conformational change may help the polymerase choose the correct nucleotide, and hence monitor DNA synthesis fidelity, through an "induced-fit" mechanism. However, the existence of open states with bound substrate and closed states without substrates suggest that substrate-induced conformational closing may be more subtle. Our dynamics simulations of two pol-beta/DNA systems (with/without substrates at the active site) reveal the large-scale closing motions of the thumb and 8-kDa subdomains in the presence of the correct substrate--leading to nearly perfect rearrangement of residues in the active site for the subsequent chemical step of nucleotidyl transfer--in contrast to an opening trend when the substrate is absent, leading to complete disassembly of the active site residues. These studies thus provide in silico evidence for the substrate-induced conformational rearrangements, as widely assumed based on a variety of crystallographic open and closed complexes. Further details gleaned from essential dynamics analyses clarify functionally relevant global motions of the polymerase-beta/DNA complex as required to prepare the system for the chemical reaction of nucleotide extension.

Simultaneous DNA Binding, Bending, and Base Flipping

We measured the kinetics of DNA bending by M.EcoRI using DNA labeled at both 5'-ends and observed changes in fluorescence resonance energy transfer. Although known to bend its cognate DNA site, energy transfer is decreased upon enzyme binding. This unanticipated effect is shown to be robust because we observe the identical decrease with different dye pairs, when the dye pairs are placed on the respective 3'-ends, the effect is cofactor- and protein-dependent, and the effect is observed with duplexes ranging from 14 through 17 base pairs. The same labeled DNA shows the anticipated increased energy transfer with EcoRV endonuclease, which also bends this sequence, and no change in energy transfer with EcoRI endonuclease, which leaves this sequence unbent. We interpret these results as evidence for an increased end-to-end distance resulting from M.EcoRI binding, mediated by a mechanism novel for DNA methyltransferases, combining DNA bending and an overall expansion of the DNA duplex. The M.EcoRI protein sequence is poorly accommodated into well defined classes of DNA methyltransferases, both at the level of individual motifs and overall alignment. Interestingly, M.EcoRI has an intercalation motif observed in the FPG DNA glycosylase family of repair enzymes. Enzyme-dependent changes in anisotropy and fluorescence resonance energy transfer have similar rate constants, which are similar to the previously determined rate constant for base flipping; thus, the three processes are nearly coincidental. Similar fluorescence resonance energy transfer experiments following AdoMet-dependent catalysis show that the unbending transition determines the steady state product release kinetics.

Biomolecular free energy profiles by a shooting/umbrella sampling protocol, "BOLAS"

We develop an efficient technique for computing free energies corresponding to conformational transitions in complex systems by combining a Monte Carlo ensemble of trajectories generated by the shooting algorithm with umbrella sampling. Motivated by the transition path sampling method, our scheme "BOLAS" (named after a cowboy's lasso) preserves microscopic reversibility and leads to the correct equilibrium distribution. This makes possible computation of free energy profiles along complex reaction coordinates for biomolecular systems with a lower systematic error compared to traditional, force-biased umbrella sampling protocols. We demonstrate the validity of BOLAS for a bistable potential, and illustrate the method's scope with an application to the sugar repuckering transition in a solvated deoxyadenosine molecule.

Critical role of magnesium ions in DNA polymerase beta's closing and active site assembly

To dissect the effects of the nucleotide-binding and catalytic metal ions on DNA polymerase mechanisms for DNA repair and synthesis, aside from the chemical reaction, we investigate their roles in the conformational transitions between closed and open states and assembly/disassembly of the active site of polymerase beta/DNA complexes before and after the chemical reaction of nucleotide incorporation. Using dynamics simulations, we find that closing before chemical reaction requires both divalent metal ions in the active site while opening after the chemical reaction is triggered by release of the catalytic metal ion. The critical closing is stabilized by the interaction of the incoming nucleotide with conserved catalytic residues (Asp190, Asp192, Asp256) and the two functional magnesium ions; without the catalytic ion, other protein residues (Arg180, Arg183, Gly189) coordinate the incomer's triphosphate group through the nucleotide-binding ion. Because we also note microionic heterogeneity near the active site, Mg(2+) and Na(+) ions can diffuse into the active site relatively rapidly, we suggest that the binding of the catalytic ion itself is not a rate-limiting conformational or overall step. However, geometric adjustments associated with functional ions and proper positioning in the active site, including subtle but systematic motions of protein side chains (e.g., Arg258), define slow or rate-limiting conformational steps that may guide fidelity mechanisms. These sequential rearrangements are likely sensitively affected when an incorrect nucleotide approaches the active site. Our suggestion that subtle and slow adjustments of the nucleotide-binding and catalytic magnesium ions help guide polymerase selection for the correct nucleotide extends descriptions of polymerase pathways and underscores the importance of the delicate conformational events both before and after the chemical reaction to polymerase efficiency and fidelity mechanisms.