Abstract
Multiple RNA strands can interact in solution and assume a large variety of configurations dictated by their potential for base pairing. Although duplex formation from two complementary oligonucleotides has been studied in detail, we still lack a systematic characterization of the behavior of higher order complexes. Here, we focus on the thermodynamic and kinetic effects of an upstream oligonucleotide on the binding of a downstream oligonucleotide to a common template, as we vary the sequence and structure of the contact interface. We show that coaxial stacking in RNA is well correlated with but much more stabilizing than helix propagation over an analogous intact double helix step (median ΔΔG°(37 °C) ≈ 1.7 kcal/mol). Consequently, approximating coaxial stacking in RNA with the helix propagation term leads to large discrepancies between predictions and our experimentally determined melting temperatures, with an offset of ≈10 °C. Our kinetic study reveals that the hybridization of the downstream probe oligonucleotide is impaired (lower k(on)) by the presence of the upstream oligonucleotide, with the thermodynamic stabilization coming entirely from an extended lifetime (lower k(off)) of the bound downstream oligonucleotide, which can increase from seconds to months. Surprisingly, we show that the effect of nicks is dependent on the length of the stacking oligonucleotides, and we discuss the binding of ultrashort (1-4 nt) oligonucleotides that are relevant in the context of the origin of life. The thermodynamic and kinetic data obtained in this work allow for the prediction of the formation and stability of higher-order multistranded complexes.