We identified over 1000 instances of water-nucleobase stacking contacts in a variety of RNA molecules from a non-redundant set of crystal structures with resolution ≤3.0 Å. Such contacts may be of either the lone pair-π (lp–π) or the OH–π type, in nature. The distribution of the distances of the water oxygen from the nucleobase plane peaks at 3.5 Å for A, G and C, and approximately at 3.1–3.2 Å for U. Quantum mechanics (QM) calculations confirm, as expected, that the optimal energy is reached at a shorter distance for the lp–π interaction as compared to the OH–π one (3.0 versus 3.5 Å). The preference of each nucleobase for either type of interaction closely correlates with its electrostatic potential map. Furthermore, QM calculations show that for all the nucleobases a favorable interaction, of either the lp–π or the OH–π type, can be established at virtually any position of the water molecule above the nucleobase skeleton, which is consistent with the uniform projection of the OW atoms over the nucleobases ring we observed in the experimental occurrences. Finally, molecular dynamics simulations of a model system for the characterization of water-nucleobase stacking contacts confirm the stability of these interactions also under dynamic conditions.

We employ density functional theory (DFT) and time-dependent DFT (TDDFT) calculations to investigate the structural, energetic and optical properties of a new computationally designed RNA alphabet, where the nucleobases, A,G, C, andU (ts-bases), have been derived by replacing sulfur with selenium in the previously reported tz-bases, based on the isothiazolo[4,3-d]pyrimidine heterocycle core. We find out that the modeled non-natural bases have minimal impact on the geometry and energetics of the classical Watson-Crick base pairs, thus potentially mimicking the natural bases in a RNA duplex in terms of H-bonding. In contrast, our calculations indicate that H-bonded base pairs involving the Hoogsteen edge of purines are destabilized as compared to their natural counterparts. We also focus on the photophysical properties of the non-natural bases and correlate their absorption/emission peaks to the strong impact of the modification on the energy of the lowest unoccupied molecular orbital. It is indeed stabilized by roughly 1.1-1.6 eV as compared to the natural analogues, resulting in a reduction of the gap between the highest occupied and the lowest unoccupied molecular orbital from 5.3-5.5 eV in the natural bases to 3.9-4.2 eV in the modified ones, with a consequent bathochromic shift in the absorption and emission spectra. Overall, our analysis clearly indicates that the newly modelled ts-bases are expected to exhibit better fluorescent properties as compared to the previously reported tz-bases, while retaining similar H-bonding properties. In addition, we show that a new RNA alphabet based on size-extended benzo-homologated ts-bases can also form stable Watson-Crick base pairs with the natural complementary nucleobases.

We present a systematic structural and energetic characterization of phosphate(OP)–nucleobase anion…π stacking interactions in RNAs. We observed OP–nucleobase stacking contacts in a variety of structural motifs other than regular helices and spanning broadly diverse sequence distances. Apart from the stacking between a phosphate and a guanine or a uracil two-residue upstream in specific U-turns, such interactions in RNA have been scarcely characterized to date. Our QM calculations showed an energy minimum at a distance between the OP atom and the nucleobase plane centroid slightly below 3 Å for all the nucleobases. By sliding the OP atom over the nucleobase plane we localized the optimal mutual positioning of the stacked moieties, corresponding to an energy minimum below -6 kcal•mol−1, for all the nucleobases, consistently with the projections of the OP atoms over the different π-rings we observed in experimental occurrences. We also found that the strength of the interaction clearly correlates with its electrostatic component, pointing to it as the most relevant contribution. Finally, as OP–uracil and OP–guanine interactions represent together 86% of the instances we detected, we also proved their stability under dynamic conditions in model systems simulated by state-of-the art DFT-MD calculations.

No abstract