Insights into the behaviour of phosphorylated DNA breaks from molecular dynamic simulations

Single-stranded breaks (SSBs) are the most frequent DNA lesions threatening genomic integrity-understanding how DNA sensor proteins recognize certain SSB types is crucial for studies of the DNA repair pathways. During repair of damaged DNA the final SSB that is to be ligated contains a 5'-phosphorylated end. The present work employed molecular simulation (MD) of DNA with a phosphorylated break in solution to address multiple questions regarding the dynamics of the break site. How does the 5'-phosphate group behave before it initiates a connection with other biomolecules? What is the conformation of the SSB site when it is likely to be recognized by DNA repair factors once the DNA repair response is triggered? And how is the structure and dynamics of DNA affected by the presence of a break? For this purpose, a series of MD simulations of 20 base pair DNAs, each with either a pyrimidine-based or purine-based break, were completed at a combined length of over 20,000 ns simulation time and compared with intact DNA of the same sequence. An analysis of the DNA forms, translational and orientational helical parameters, local break site stiffness, bending angles, 5'-phosphate group orientation dynamics, and the effects of the protonation state of the break site phosphate group provides insights into the mechanism for the break site recognition.