Sequence-Dependent DNA Base Selection Fidelity: A Kinetics-based Model and its validation

DNA replication achieves error rates as low as 10-9-10-11 per base pair through the combined contribution of base selection, exonuclease proofreading, and mismatch repair. Among those processes, base selection determines the initial level of accuracy and is strongly modulated by local sequence context. Existing models address this dependence either by fitting individual rate constants for each sequence context or by invoking the global template properties, neither of which derives sequence dependence from the underlying thermodynamics and kinetics of base pair formation. Here we present a mechanism for sequence-dependent base selection fidelity, built from two physical properties: nearest-neighbor stacking thermodynamics and directional kinetic asymmetry. The model fits the experimentally observed mutation spectra from three mismatch repair-deficient organisms well (r=0.74, 0.70, and 0.63), and predicts that base-selection accuracy varies non-monotonically with temperature in a sequence-dependent manner. Our model, therefore, provides a framework that connects sequence-dependent thermodynamic and kinetic effects during nucleotide incorporation to experimentally observed mutation rates.