Temperature-dependent twist of double-stranded RNA probed by magnetic tweezers experiments and molecular dynamics simulations

Abstract

RNA plays critical roles in the transmission and regulation of genetic information and is increasingly used in biomedical and biotechnological applications. Functional RNAs contain extended double-stranded regions and the structure of double-stranded RNA (dsRNA) has been revealed at high-resolution. However, the dependence of the properties of the RNA double helix on environmental effects, notably temperature, is still poorly understood. Here, we use single-molecule magnetic tweezers measurements to determine the dependence of the dsRNA twist on temperature. We find that dsRNA unwinds with increasing temperature, even more than DNA, with {Delta}TwRNA = -14.4 {+/-} 0.7 {o}/({degrees}C{middle dot}kbp), compared to {Delta}TwDNA = -11.0 {+/-} 1.2 {o}/({degrees}C{middle dot}kbp). All-atom molecular dynamics (MD) simulations using a range of nucleic acid force fields, ion parameters, and water models correctly predict that dsRNA unwinds with rising temperature, but significantly underestimate the magnitude of the effect. These MD data, together with additional MD simulations involving DNA and DNA-RNA hybrid duplexes, reveal a linear correlation between twist temperature decrease and the helical rise, in line with DNA but at variance with RNA experimental data. We speculate that this discrepancy might be caused by some unknown bias in the RNA force fields tested, or by as yet undiscovered transient alternative structures in the RNA duplex. Our results provide a baseline to model more complex RNA assemblies and to test and develop new parameterizations for RNA simulations. They may also inspire physical models of temperature-dependent dsRNA structure.