Kinetic mechanisms for the sequence dependence of transcriptional errors

The fidelity of template-dependent mRNA synthesis during transcription elongation is the primary determinant of accurate gene expression and the maintenance of functional RNA transcripts. However, the mechanisms governing transcription fidelity remain incompletely understood. While previous studies have characterized how error rates vary with nucleotide identity at upstream and downstream positions from the incorporation site, the comprehensive microscopic explanation of this sequence dependence has not been elucidated. In this study, we develop a novel theoretical approach that integrates transcription proofreading mechanisms and inhomogenous DNA sequence effects. Using first-passage analysis validated by Monte Carlo simulations, we quantitatively characterize nucleotide-specific error rates during RNA polymerase II transcription. The model accurately reproduces experimental error rates and predicts kinetic parameters influencing transcriptional fidelity. Analysis reveals nucleotide incorporation rates follow the hierarchy U < C < G < A, consistent with independent experimental observations. Notably, our model not only explains how the error rates depend on the nature of the base immediately downstream (+1) but also predicts that the identity of the nucleotide at the second downstream position (+2) also plays an important role. Pyrimidines at position +2 contribute to lower error rates than purines, whereas the third downstream base (+3) has no effect. These previously unreported correlations are corroborated by bioinformatic analysis of existing datasets. In addition, using the BRCA1 gene as an example, we explore the physiological implications of sequence-dependent error rates, identifying an increased likelihood of premature stop codon errors. These findings clarify how DNA sequence context modulates nucleotide incorporation kinetics, advancing our understanding of transcriptional fidelity and its functional consequences.