ABSTRACTGenetic information is encoded in the DNA double helix which, in its physiological milieu, is characterized by the iconical Watson-Crick nucleobase pairing. Recent NMR relaxation experiments revealed the transient presence of an alternative, Hoogsteen base pairing pattern in naked DNA duplexes and estimated its relative stability and lifetime. In contrast, HG transitions in RNA were not observed. Understanding Hoogsteen (HG) base pairing is important because the underlying "breathing" can modulate significantly DNA/RNA recognition by proteins. However, a detailed mechanistic insight into the transition pathways and kinetics is still missing. We performed enhanced sampling simulation (with combined metadynamics and adaptive force bias method) and Markov State modeling to obtain accurate free energy, kinetics and the intermediates in the transition pathway between WC and HG base pair for both naked B-DNA and A-RNA duplexes. The Markov state model constructed from our unbiased MD simulation data revealed previously unknown complex extra-helical intermediates in this seemingly simple process of base pair conformation switching in B-DNA. Extending our calculation to A-RNA, for which HG base pair is not observed experimentally, resulted in relatively unstable single hydrogen bonded distorted Hoogsteen like base pair. Unlike B-DNA the transition pathway primarily involved base paired and intra-helical intermediates with transition timescales much higher than that of B-DNA. The seemingly obvious flip-over reaction coordinate, i.e., the glycosidic torsion angle is unable to resolve the intermediates; so a multidimensional picture, involving backbone dihedral angles and distance between atoms participating in hydrogen bonds, is required to gain insight into the molecular mechanism.SIGNIFICANCEFormation of unconventional Hoogsteen (HG) base pairing is an important problem in DNA biophysics owing to its key role in facilitating the binding of DNA repairing enzymes, proteins and drugs to damaged DNA. X-ray crystallography and NMR relaxation experiments revealed the presence of HG base pair in naked DNA duplex and protein-DNA complex but no HG base pair was observed in RNA. Molecular dynamics simulations could reproduce the experimental free energy cost of HG base pairing in DNA although a detailed mechanistic insight is still missing. We performed enhanced sampling simulation and Markov state modeling to obtain accurate free energy, kinetics and the intermediates in the transition pathway between WC and HG base pair for both B-DNA and A-RNA.
Hoogsteen Base Pairing in DNA vs RNA: Thermodynamics and Kinetics from Enhanced Sampling Simulation and Markov State Modeling
