Abstract
DNA-RNA hybrid duplexes play essential roles during the reverse transcription of RNA viruses and DNA replication. The opening and conformation changes of individual base pairs are critical to their biological functions. However, the microscopic mechanisms governing base pair closing and opening at the atomic level remain poorly understood. In this study, we investigated the thermodynamic and kinetic parameters of the dA-rU base pair in a DNA-RNA hybrid duplex using 4 μs all-atom molecular dynamics (MD) simulations at different temperatures. Our results showed that the thermodynamic parameters of the dA-rU base pair aligned with the predictions of the nearest-neighbor model and were close to those of the AU base pair in RNA. The temperature dependence of the average lifetimes of both the open and the closed states, as well as the transition path times, were obtained. The free-energy barrier for a single base pair opening and closing arises from an increase in enthalpy due to the disruption of the base-stacking interactions and hydrogen bonding, along with an entropy loss attributed to the accompanying restrictions, such as torsional angle constraints and solvent viscosity.