Relating DNA base-pairing in aqueous media to DNA polymerase fidelity

Modifying the Basicity of the dNTP Leaving Group Modulates Precatalytic Conformational Changes of DNA Polymerase β

The catalytic function of DNA polymerase β (pol β) fulfills the gap-filling requirement of the base excision DNA repair pathway by incorporating a single nucleotide into a gapped DNA substrate resulting from the removal of damaged DNA bases. Most importantly, pol β can select the correct nucleotide from a pool of similarly structured nucleotides to incorporate into DNA in order to prevent the accumulation of mutations in the genome. Pol β is likely to employ various mechanisms for substrate selection. Here, we use dCTP analogues that have been modified at the β,γ-bridging group of the triphosphate moiety to monitor the effect of leaving group basicity of the incoming nucleotide on precatalytic conformational changes, which are important for catalysis and selectivity. It has been previously shown that there is a linear free energy relationship between leaving group pKa and the chemical transition state. Our results indicate that there is a similar relationship with the rate of a precatalytic conformational change, specifically, the closing of the fingers subdomain of pol β. In addition, by utilizing analogue β,γ-CHX stereoisomers, we identified that the orientation of the β,γ-bridging group relative to R183 is important for the rate of fingers closing, which directly influences chemistry.

A Deep Dive into DNA Base Pairing Interactions Under Water

Base pairing plays a pivotal role in DNA functions and replication fidelity. But while the complementarity between Watson-Crick matched bases is generally believed to arise from the different number of hydrogen bonds in G|C pairs versus A|T, the energetics of these interactions are heavily renormalized by the aqueous solvent. Employing large-scale Monte Carlo simulations, we have extracted the solvent contribution to the free energy for canonical and some noncanonical and stacked base pairs. For all of them, the solvent's contribution to the base pairing free energy is exclusively destabilizing. While the direct hydrogen bonding interactions in the G|C pair is much stronger than A|T, the thermodynamic resistance produced by the solvent also pushes back much stronger against G|C compared to A|T, generating an only ∼1 kcal/mol free energy difference between them. We have profiled the density of water molecules in the solvent adjacent to the bases and observed a "freezing" behavior where waters are recruited into the gap between the bases to compensate for the unsatisfied hydrogen bonds between them. A very small number of water molecules that are associated with the Watson-Crick donor/acceptor atoms turn out to be responsible for the majority of the solvent's thermodynamic resistance to base pairing. The absence or presence of these near-field waters can be used to enhance fidelity during DNA replication.

Elucidating the Determinants of Polymerase Specificity by Microfluidic-Based Deep Mutational Scanning

Engineering polymerases to synthesize artificial genetic polymers with unique backbone structures is limited by a general lack of understanding about the structural determinants that govern substrate specificity. Here, we report a high-throughput microfluidic-based approach for mapping sequence-function relationships that combines droplet-based optical polymerase sorting with deep mutational scanning. We applied this strategy to map the finger subdomain of a replicative DNA polymerase isolated from Thermococcus kodakarensis (Kod). The enrichment profile provides an unbiased view of the ability of each mutant to synthesize threose nucleic acid, which was used as a model non-natural genetic polymer. From a single round of sorting, we discovered two cases of positive epistasis and demonstrate the near inversion of substrate specificity from a double mutant variant. This effort indicates that polymerase specificity may be governed by a small number of highly specific residues that can be elucidated by deep mutational scanning without the need for iterative rounds of directed evolution.

A Transition-State Perspective on Y-Family DNA Polymerase η Fidelity in Comparison with X-Family DNA Polymerases λ and β

Deoxynucleotide misincorporation efficiencies can span a wide 10⁴-fold range, from ∼10^-2 to ∼10^-6, depending principally on polymerase (pol) identity and DNA sequence context. We have addressed DNA pol fidelity mechanisms from a transition-state (TS) perspective using our "tool-kit" of dATP- and dGTP-β,γ substrate analogues in which the pyrophosphate leaving group (p K_a4 = 8.9) has been replaced by a series of bisphosphonates covering a broad acidity range spanning p K_a4 values from 7.8 (CF₂) to 12.3 [C(CH₃)₂]. Here, we have used a linear free energy relationship (LFER) analysis, in the form of a Brønsted plot of log( k_pol) versus p K_a4, for Y-family error-prone pol η and X-family pols λ and β to determine the extent to which different electrostatic active site environments alter k_pol values. The apparent chemical rate constant ( k_pol) is the rate-determining step for the three pols. The pols each exhibit a distinct catalytic signature that differs for formation of right (A·T) and wrong (G·T) incorporations observed as changes in slopes and displacements of the Brønsted lines, in relation to a reference LFER. Common to this signature among all three pols is a split linear pattern in which the analogues containing two halogens show k_pol values that are systematically lower than would be predicted from their p K_a4 values measured in aqueous solution. We discuss how metal ions and active site amino acids are responsible for causing "effective" p K_a4 values that differ for dihalo and non-dihalo substrates as well as for individual R and S stereoisomers for CHF and CHCl.