Computational Study of the Force Dependence of Phosphoryl Transfer during DNA Synthesis by a High Fidelity Polymerase

High fidelity polymerases are efficient catalysts of phosphodiester bond formation during DNA replication or repair. We interpret molecular dynamics simulations of a polymerase bound to its substrate DNA and incoming nucleotide using a quasiharmonic model to study the effect of external forces applied to the bound DNA on the kinetics of phosphoryl transfer. The origin of the force dependence is shown to be an intriguing coupling between slow, delocalized polymerase-DNA modes and fast catalytic site motions. Using noncognate DNA substrates we show that the force dependence is context specific

The effect of oxidatively damaged DNA on the active site pre-organization during nucleotide incorporation in a high fidelity polymerase from \u3cem\u3eBacillus stearothermophilus\u3c/em\u3e

We study the effect of the oxidative lesion 8-oxoguanine (8oxoG) on the pre-organization of the active site for DNA replication in the closed (active) state of the Bacillus Fragment (BF), a Klenow analog from Bacillus stearothermophilus. Our molecular dynamics and free energy simulations of explicitly solvated model ternary complexes of BF bound to correct dCTP/incorrect dATP opposite guanine (G) and 8oxoG bases in DNA suggest that the lesion introduces structural and energetic changes at the catalytic site to favor dATP insertion. Despite the formation of a stable Watson-Crick pairing in the 8oxoG:dCTP system, the catalytic geometry is severely distorted to possibly slow down catalysis. Indeed, our calculated free energy landscapes associated with active site pre-organization suggest additional barriers to assemble an efficient catalytic site, which need to be overcome during dCTP incorporation opposite 8oxoG relative to that opposite undamaged G. In contrast, the catalytic geometry for the Hoogsteen pairing in the 8oxoG:dATP system is highly organized and poised for efficient nucleotide incorporation via the twometal- ion catalyzed phosphoryl transfer mechanism. However, the free energy calculations suggest that the catalytic geometry during dATP incorporation opposite 8oxoG is considerably less plastic than that during dCTP incorporation opposite G despite a very similar, well organized catalytic site for both systems. A correlation analysis of the dynamics trajectories suggests the presence of significant coupling between motions of the polymerase fingers and the primary distance for nucleophilic attack (i.e., between the terminal primer O3´ and the dNTP Pα atoms) during correct dCTP incorporation opposite undamaged G. This coupling is shown to be disrupted during nucleotide incorporation by the polymerase with oxidatively damaged DNA/dNTP substrates. We also suggest that the lesion affects DNA interactions with key polymerase residues, thereby affecting the enzymes ability to discriminate against noncomplementary DNA/dNTP substrates. Taken together, our results provide a unified structural, energetic, and dynamic platform to rationalize experimentally observed relative nucleotide incorporation rates for correct dCTP/incorrect dATP insertion opposite an undamaged/oxidatively damaged template G by BF

Calculation of free energies in fluid membranes subject to heterogeneous curvature fields

We present a computational methodology for incorporating thermal effects and calculating relative free energies for elastic fluid membranes subject to spatially dependent intrinsic curvature fields using the method of thermodynamic integration. Based on a simple model for the intrinsic curvature imposed only in a localized region of the membrane, we employ thermodynamic integration to calculate the free-energy change as a function of increasing strength of the intrinsic curvature field and a thermodynamic cycle to compute free-energy changes for different sizes of the localized region. By explicitly computing the free-energy changes and by quantifying the loss of entropy accompanied with increasing membrane deformation, we show that the membrane stiffness increases with increasing intrinsic field, thereby, renormalizing the membrane bending rigidity. The second main conclusion of this work is that the entropy of the membrane decreases with increasing size of the localized region subject to the curvature field. Our results help to quantify the free-energy change when a planar membrane deforms under the influence of curvature-inducing proteins at a finite temperature