Simulations of RNA base pairs in a nanodroplet reveal solvation-dependent stability

A Hitchhiker's guide to RNA-RNA structure and interaction prediction tools

RNA biology has risen to prominence after a remarkable discovery of diverse functions of noncoding RNA (ncRNA). Most untranslated transcripts often exert their regulatory functions into RNA-RNA complexes via base pairing with complementary sequences in other RNAs. An interplay between RNAs is essential, as it possesses various functional roles in human cells, including genetic translation, RNA splicing, editing, ribosomal RNA maturation, RNA degradation and the regulation of metabolic pathways/riboswitches. Moreover, the pervasive transcription of the human genome allows for the discovery of novel genomic functions via RNA interactome investigation. The advancement of experimental procedures has resulted in an explosion of documented data, necessitating the development of efficient and precise computational tools and algorithms. This review provides an extensive update on RNA-RNA interaction (RRI) analysis via thermodynamic- and comparative-based RNA secondary structure prediction (RSP) and RNA-RNA interaction prediction (RIP) tools and their general functions. We also highlighted the current knowledge of RRIs and the limitations of RNA interactome mapping via experimental data. Then, the gap between RSP and RIP, the importance of RNA homologues, the relationship between pseudoknots, and RNA folding thermodynamics are discussed. It is hoped that these emerging prediction tools will deepen the understanding of RNA-associated interactions in human diseases and hasten treatment processes.

Research of Water Molecules Cluster Structuring during Haberlea rhodopensis Friv. Hydration

Gesneriaceae plant family is comprised of resurrection species, namely Boea hygrometrica and Paraboea rufescens, that are native to the Southeast Asia and Haberlea rhodopensis, Ramonda myconi, and Ramonda serbica, which are mainly found in the Balkan Peninsula. Haberlea rhodopensis is known to be able to survive extreme and prolonged dehydration. Study was carried out after the dried plant Haberlea rhodopensis Friv. had been hydrated and had reached its fresh state. Two juice samples were collected from the plant blossom: The first sample was prepared with 1% filtered water through a patented EVOdrop device. Then the sample was saturated with hydrogen with EVOdrop booster to a concentration of 1.2 ppm, pH = 7.3, ORP = -390 mV. This first sample was prepared with filtered tap water from Sofia, Bulgaria. The second sample, which was a control one, was developed with tap water from Sofia, Bulgaria, consisting of 1% solutions of Haberlea rhodopensis. A study revealed that during the drying process in H. rhodopensis the number of free water molecules decreases, and water dimers are formed. The aim of our study was to determine the number of water molecules in clusters in 1% solutions of hydrated H. rhodopensis plants. Results were analyzed according to the two types of water used in the experiment. Th EVOdrop device is equipped with an ultranano membrane and rotating jet nozzle to create a vortex water and saturation thanks to a second device EVObooster to obtain hydrogen-rich water. In the current study Hydrogen-rich water is referred to as Hydrogen EVOdrop Water (HEW). Research was conducted using the following methods-spectral methods non-equilibrium energy spectrum (NES) and differential non-equilibrium energy spectrum (DNES), mathematical models, and study of the distribution of water molecules in water clusters. In a licensed Eurotest Laboratory, the research of tap water before and after flowing through the EVOdrop device was proven. Studies have been carried out on the structuring of water molecule clusters after change of hydrogen bond energies. The restructuring comes with rearrangement of water molecules by the energy levels of hydrogen bonds. Local extrema can be observed in the spectrum with largest amount of water molecules. The structural changes were tested using the NES and DNES spectral methods. The conducted research proved that the application of EVOdrop device and EVObooster changes the parameters of water to benefit hydration and health.

Temperature-Induced Change of Water Structure in Aqueous Solutions of Some Kosmotropic and Chaotropic Salts

The hydrogen bond structure of water was examined by comparing the temperature dependent OH-stretching bands of water and aqueous NaClO₄, KClO₄, Na₂SO₄, and K₂SO₄ solutions. Results called attention to the role of cations on top of the importance of anions determining the emerging structure of a multi-layered system consisting single water rings or multi-ring water-clusters.

Use of a redox probe for an electrochemical RNA-ligand binding assay in microliter droplets

In this work, we report an affordable, sensitive, fast and user-friendly electroanalytical method for monitoring the binding between unlabeled RNA and small compounds in microliter-size droplets using a redox-probe and disposable miniaturized screen-printed electrochemical cells.

Water evaporation and conformational changes from partially solvated ubiquitin

Using molecular dynamics simulation, we study the evaporation of water molecules off partially solvated ubiquitin. The evaporation and cooling rates are determined for a molecule at the initial temperature of 300 K. The cooling rate is found to be around 3 K/ns, and decreases with water temperature in the course of the evaporation. The conformation changes are monitored by studying a variety of intermediate partially solvated ubiquitin structures. We find that ubiquitin shrinks with decreasing hydration shell and exposes more of its hydrophilic surface area to the surrounding.

Computational and single-molecule force studies of a macro domain protein reveal a key molecular determinant for mechanical stability

Resolving molecular determinants of mechanical stability of proteins is crucial in the rational design of advanced biomaterials for use in biomedical and nanotechnological applications. Here we present an interdisciplinary study combining bioinformatics screening, steered molecular dynamics simulations, protein engineering, and single-molecule force spectroscopy that explores the mechanical properties of a macro domain protein with mixed alpha + beta topology. The unique architecture is defined by a single seven-stranded beta-sheet in the core of the protein flanked by five alpha-helices. Unlike mechanically stable proteins studied thus far, the macro domain provides the distinct advantage of having the key load-bearing hydrogen bonds (H bonds) buried in the hydrophobic core protected from water attacks. This feature allows direct measurement of the force required to break apart the load-bearing H bonds under locally hydrophobic conditions. Steered molecular dynamics simulations predicted extremely high mechanical stability of the macro domain by using constant velocity and constant force methods. Single-molecule force spectroscopy experiments confirm the exceptional mechanical strength of the macro domain, measuring a rupture force as high as 570 pN. Furthermore, through selective deletion of shielding peptide segments, we examined the same key H bonds under hydrophilic environments in which the beta-strands are exposed to solvent and verify that the high mechanical stability of the macro domain results from excellent shielding of the load-bearing H bonds from competing water. Our study reveals that shielding water accessibility to the load-bearing strands is a critical molecular determinant for enhancing the mechanical stability of proteins.

Assessing the performance of implicit solvation models at a nucleic acid surface

Implicit solvation models are popular alternatives to explicit solvent methods due to their ability to "pre-average" solvent behavior and thus reduce the need for computationally-expensive sampling. Previously, we have demonstrated that Poisson-Boltzmann models for polar solvation and integral-based models for nonpolar solvation can reproduce explicit solvation forces in a low-charge density protein system. In the present work, we examine the ability of these continuum models to describe solvation forces at the surface of a RNA hairpin. While these models do not completely describe all of the details of solvent behavior at this highly-charged biomolecular interface, they do provide a reasonable description of average solvation forces and therefore show significant promise for developing more robust implicit descriptions of solvent around nucleic acid systems for use in biomolecular simulation and modeling. Additionally, we observe fairly good transferability in the nonpolar model parameters optimized for protein systems, suggesting its robustness for modeling general nonpolar solvation phenomena in biomolecular systems.

Microhydration of the guanine-cytosine (GC) base pair in the neutral and anionic radical states: a density functional study

A density functional study of the effects of microhydration on the guanine-cytosine (GC) base pair and its anion radical is presented. Geometries of the GC base pair in the presence of 6 and 11 water molecules were fully optimized in the neutral (GC-nH2O) and anion radical [(GC-nH2O)*-] (n = 6 and 11) states using the B3LYP method and the 6-31+G** basis set. Further, vibrational frequency analysis at the same level of theory (B3LYP/6-31+G**) was also performed to ensure the existence of local minima in these hydrated structures. It was found that water molecules surrounding the GC base pair have significant effects on the geometry of the GC base pair and promote nonplanarity in the GC base pair. The calculated structures were found to be in good agreement with those observed experimentally and obtained in molecular dynamics (MD) simulation studies. The water molecules in neutral GC-nH2O complexes lie near the ring plane of the GC base pair where they undergo hydrogen bonding with both GC and each other. However, in the GC anion radical complexes (GC-nH2O, n = 6, 11), the water molecules are displaced substantially from the GC ring plane. For GC-11H2O*-, a water molecule is hydrogen-bonded with the C6 atom of the cytosine base. We found that the hydration shell initially destabilizes the GC base pair toward electron capture as a transient anion. Energetically unstable diffuse states in the hydration shell are suggested to provide an intermediate state for the excess electron before molecular reorganization of the water molecules and the base pair results in a stable anion formation. The singly occupied molecular orbital (SOMO) in the anion radical complexes clearly shows that an excess electron localizes into a pi orbital of cytosine. The zero-point-energy (ZPE-) corrected adiabatic electron affinities (AEAs) of the GC-6H2O and GC-11H2O complexes, at the B3LYP/6-31+G** level of theory, were found to be 0.74 and 0.95 eV, respectively. However, the incorporation of bulk water as a solvent using the polarized continuum model (PCM) increases the EAs of these complexes to 1.77 eV.