The specific folding pattern and function of RNA molecules lies in various weak interactions, in addition to the strong base-base pairing and stacking. One of these relatively weak interactions, characterized by the stacking of the O4′ atom of a ribose on top of the heterocycle ring of a nucleobase, has been known to occur but has largely been ignored in the description of RNA structures. We identified 2015 ribose–base stacking interactions in a high-resolution set of non-redundant RNA crystal structures. They are widespread in structured RNA molecules and are located in structural motifs other than regular stems. Over 50% of them involve an adenine, as we found ribose-adenine contacts to be recurring elements in A-minor motifs. Fewer than 50% of the interactions involve a ribose and a base of neighboring residues, while approximately 30% of them involve a ribose and a nucleobase at least four residues apart. Some of them establish inter-domain or inter-molecular contacts and often implicate functionally relevant nucleotides. In vacuo ribose-nucleobase stacking interaction energies were calculated by quantum mechanics methods. Finally, we found that lone pair–π stacking interactions also occur between ribose and aromatic amino acids in RNA–protein complexes.

Base pairs involving protonated nucleobases play important roles in mediating global macromolecular conformational changes and in facilitation of catalysis in a variety of functional RNA molecules. Here we present our attempts at understanding the role of such base pairs by detecting possible protonated base pairs in the available RNA crystal structures using BPFind software, in their specific structural contexts, and by the characterization of their geometries, interaction energies, and stabilities using advanced quantum chemical computations. We report occurrences of 18 distinct protonated base pair combinations from a representative data set of RNA crystal structures and propose a theoretical model for one putative base pair combination. Optimization of base pair geometries was carried out at the B3LYP/cc-pVTZ level, and the BSSE corrected interaction energies were calculated at the MP2/aug-cc-pVDZ level of theory. The geometries for each of the base pairs were characterized in terms of H-bonding patterns observed, rmsd values observed on optimization, and base pair geometrical parameters. In addition, the intermolecular interaction in these complexes was also analyzed using Morokuma energy decomposition. The gas phase interaction energies of the base pairs range from -24 to -49 kcal/mol and reveal the dominance of Hartree-Fock component of interaction energy constituting 73% to 98% of the total interaction energy values. On the basis of our combined bioinformatics and quantum chemical analysis of different protonated base pairs, we suggest resolution of structural ambiguities and correlate their geometric and energetic features with their structural and functional roles. In addition, we also examine the suitability of specific base pairs as key elements in molecular switches and as nucleators for higher order structures such as base triplets and quartets.

The G:C reverse Watson–Crick (W:W trans) base pair, also known as Levitt base pair in the context of tRNAs, is a structurally and functionally important base pair that contributes to tertiary interactions joining distant domains in functional RNA molecules and also participates in metabolite binding in riboswitches. We previously indicated that the isolated G:C W:W trans base pair is a rather unstable geometry, and that dicationic metal binding to the Guanine base or posttranscriptional modification of the Guanine can increase its stability. Herein, we extend our survey and report on other H-bonding interactions that can increase the stability of this base pair. To this aim, we performed a bioinformatics search of the PDB to locate all the occurencies of G:C trans base pairs. Interestingly, 66% of the G:C trans base pairs in the PDB are engaged in additional H-bonding interactions with other bases, the RNA backbone or structured water molecules. High level quantum mechanical calculations on a data set of representative crystal structures were performed to shed light on the structural stability and energetics of the various crystallographic motifs. This analysis was extended to the binding of the preQ1 metabolite to a preQ1-II riboswitch.

No abstract

No abstract

Posttranscriptional modifications greatly enhance the chemical information of RNA molecules, contributing to explain the diversity of their structures and functions. A significant fraction of RNA experimental structures available to date present modified nucleobases, with half of them being involved in H-bonding interactions with other bases, i.e. ‘modified base pairs’. Herein we present a systematic investigation of modified base pairs, in the context of experimental RNA structures. To this end, we first compiled an atlas of experimentally observed modified base pairs, for which we recorded occurrences and structural context. Then, for each base pair, we selected a representative for subsequent quantum mechanics calculations, to find out its optimal geometry and interaction energy. Our structural analyses show that most of the modified base pairs are non Watson–Crick like and are involved in RNA tertiary structure motifs. In addition, quantum mechanics calculations quantify and provide a rationale for the impact of the different modifications on the geometry and stability of the base pairs they participate in.

We present theoretical characterization of fluorescent non-natural nucleobases, (tz)A, (tz)G, (tz)C, and (tz)U, derived from the isothiazolo[4,3-d]pyrimidine heterocycle. Consistent with the experimental evidence, our calculations show that the non-natural bases have minimal impact on the geometry and stability of the classical Watson-Crick base pairs, allowing them to accurately mimic 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 are destabilized relative to their natural counterparts. Analysis of the photophysical properties of the non-natural bases allowed us to correlate their absorption/emission peaks to the strong impact of the modification on the energy of the lowest unoccupied molecular orbital, LUMO, which is stabilized by roughly 1.0-1.2 eV relative to the natural analogues, while the highest occupied molecular orbital, HOMO, is not substantially affected. As a result, the HOMO-LUMO gap is reduced from 5.3-5.5 eV in the natural bases to 4.0-4.4 eV in the modified ones, with a consequent bathochromic shift in the absorption and emission spectra.

Researchers have suggested that the reason behind infertility is pernicious effect of broad spectrum pesticides on non target, beneficial microorganism of soil. Here, studying the chelating effect of selective organophosphate and carbamate pesticides with essential metal ions, at all possible combinations of three different pH (4 ± 0.05, 7 ± 0.05 and 9 ± 0.05) and three different temperatures (15 ± 0.5°C, 30 ± 0.5°C and 45 ± 0.5°C), shows very fast rate of reaction which further increases with increase of pH and temperature. Carbonyl oxygen of carbamate and phosphate oxygen of organophosphate were found to be common ligating sites among all the complexes. Formed metal complexes were found to be highly stable and water insoluble on interaction with essential metal ions in solvent medium as well as over silica. Density functional theory (DFT) calculations not only reinforced the experimental observations, but, after a wide computational conformational analysis, unraveled the nature of the high stable undesired species that consist of pesticides complexed by metal ions from the soil. All in all, apart from the direct toxicity of pesticides, the indirect effect by means of complexation of free metal ions impoverishes the soil.