Chargaff's second parity rule and the kinetics of DNA replication

This paper presents the study of a DNA replication model grounded in the biochemical kinetics of DNA polymerases, which copy each DNA strand into a complementary strand, except for rare point-like mutations caused by nucleotide substitution errors. Numerical simulations of many successive replications, starting from an arbitrary initial DNA sequence, show that the fractions of mono- and oligonucleotides converge toward compliance with Chargaff’s second parity rule. The theoretical framework developed for this multireplication process demonstrates that the near-equalities of complementary nucleotide fractions arise from two key features: (1) the dominant role of base-pair complementarity in replication kinetics and (2) the low intrinsic error rate of DNA polymerases. Together, these two features yield a robust mechanistic basis for Chargaff’s second parity rule. These considerations explain the existence of deviations with respect to the predictions of models assuming no-strand-bias conditions.

💡 Research Summary

The paper investigates the origin of Chargaff’s second parity rule (CSPR), which states that within each single strand of double‑stranded DNA the frequencies of complementary bases are approximately equal (A≈T, C≈G). While Chargaff’s first rule follows directly from Watson‑Crick base pairing, the second rule has long been considered an emergent property of genome evolution, lacking a mechanistic explanation.

To address this gap, the authors construct a kinetic model of DNA replication grounded in experimentally measured parameters of several representative DNA polymerases: the B‑family polymerases Dpo1, Dpo3, Dpo4 from the archaeon Sulfolobus solfataricus; a D‑family polymerase from Thermococcus sp.; and the eukaryotic polymerase β from rat. For each enzyme they use pre‑steady‑state kinetic data to obtain the polymerization rate constants (k⁺ₘₙ) and dissociation constants (Kₘₙ) for all 16 possible base‑pairings, both correct (Watson‑Crick) and mismatched. The data reveal that correct pairings have rates orders of magnitude larger than mismatches, resulting in low intrinsic error probabilities (η≈10⁻⁴–10⁻³).

Replication is simulated using Gillespie’s stochastic kinetic Monte‑Carlo algorithm. An initial template of length L = 10⁶ nucleotides is generated as a Bernoulli chain with arbitrary base fractions. Each replication proceeds from the 3′→5′ direction of the template, adding nucleotides to the growing copy according to the kinetic rates. After a copy is completed, it becomes the template for the next round, allowing the authors to iterate the process tens of thousands of times. Simulations are performed under three nucleotide‑concentration regimes: (I) physiological resting‑cell levels, (II) dividing‑cell levels, and (III) saturating concentrations.

Across all polymerases, concentration conditions, and initial compositions, the simulations show a robust convergence of mononucleotide frequencies toward A ≈ T and C ≈ G. Even when the starting strand is highly biased (e.g., A = 40 %, T = 10 %, C = 30 %, G = 20 %), after many replication cycles the frequencies approach the averages (A₀+T₀)/2 and (C₀+G₀)/2, with residual differences on the order of the error probability η. Oligonucleotide frequencies (dinucleotides and trinucleotides) exhibit the same trend, confirming the extension of CSPR to short motifs.

The authors then develop an analytical framework. By treating the replication process as a Markov chain, they construct a transition matrix whose stationary distribution corresponds to the asymptotic nucleotide composition. The matrix’s dominant eigenvector is shown to be symmetric with respect to complementary bases, with asymmetries proportional to η. This demonstrates mathematically that (1) the dominance of Watson‑Crick complementarity in the kinetic rates forces the system toward strand symmetry, and (2) the low error rate ensures that deviations are minute. The analysis is first carried out for a “memoryless” model where the attachment/detachment rates depend only on the current template base; a later extension incorporates dependence on the previously incorporated nucleotide (nearest‑neighbor effects). Even with such context dependence, the convergence to CSPR persists, confirming the robustness of the mechanism.

The paper concludes that CSPR does not require any selective advantage, structural constraints, or special evolutionary pressures; it is an inevitable consequence of the fundamental chemistry of DNA polymerization—high fidelity and strong complementarity. Deviations observed in mitochondrial genomes, single‑stranded viruses, or organisms with unusually high polymerase error rates are naturally explained as departures from the low‑η regime. The authors suggest that future work could incorporate exonuclease proofreading, biased mutation spectra, or non‑equilibrium nucleotide pools to model the few known exceptions more precisely.

Overall, the study provides a compelling kinetic‑theoretic foundation for Chargaff’s second parity rule, linking microscopic enzymatic properties to a macroscopic genomic pattern observed across the tree of life.