The intracellular environment is crowded and heterogeneous. Although the thermodynamic stability of nucleic acid duplexes is predictable in dilute solutions, methods of predicting such stability under specific intracellular conditions are not yet available. We recently showed that the nearest-neighbor model for self-complementary DNA is valid under molecular crowding condition of 40% polyethylene glycol with an average molecular weight of 200 (PEG 200) in 100 mM NaCl. Here, we determined nearest-neighbor parameters for DNA duplex formation under the same crowding condition to predict the thermodynamics of DNA duplexes in the intracellular environment. Preferential hydration of the nucleotides was found to be the key factor for nearest-neighbor parameters in the crowding condition. The determined parameters were shown to predict the thermodynamic parameters (∆H°, ∆S°, and ∆G°₃₇) and melting temperatures (T_m) of the DNA duplexes in the crowding condition with significant accuracy. Moreover, we proposed a general method for predicting the stability of short DNA duplexes in different cosolutes based on the relationship between duplex stability and the water activity of the cosolute solution. The method described herein would be valuable for investigating biological processes that occur under specific intracellular crowded conditions and for the application of DNA-based biotechnologies in crowded environments. Significance The intracellular environment is occupied by a high concentration of various types of biomolecules that lead to a situation known as molecular crowding, which affects important biological reactions such as replication and transcription by altering the stability of nucleic acids. Because intracellular crowded conditions are heterogeneous, a general method is needed to predict the stability of nucleic acids in such conditions. Here, we predicted the stability of DNA duplexes under crowded conditions induced by a large amount of a cosolute at a physiological salt concentration using the nearest-neighbor model. We also extended the parameters for predicting DNA stability under conditions imposed by different cosolutes that mimicked a localized intracellular crowded region.

The stability of Watson–Crick paired RNA/DNA hybrids is important for designing optimal oligonucleotides for ASO (Antisense Oligonucleotide) and CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)–Cas9 techniques. Previous nearest-neighbour (NN) parameters for predicting hybrid stability in a 1 M NaCl solution, however, may not be applicable for predicting stability at salt concentrations closer to physiological condition (e.g. ∼100 mM Na⁺ or K⁺ in the presence or absence of Mg²⁺). Herein, we report measured thermodynamic parameters of 38 RNA/DNA hybrids at 100 mM NaCl and derive new NN parameters to predict duplex stability. Predicted ΔG°₃₇ and T_m values based on the established NN parameters agreed well with the measured values with 2.9% and 1.1°C deviations, respectively. The new results can also be used to make precise predictions for duplexes formed in 100 mM KCl or 100 mM NaCl in the presence of 1 mM Mg²⁺, which can mimic an intracellular and extracellular salt condition, respectively. Comparisons of the predicted thermodynamic parameters with published data using ASO and CRISPR–Cas9 may allow designing shorter oligonucleotides for these techniques that will diminish the probability of non-specific binding and also improve the efficiency of target gene regulation.

During translation, intracellular mRNA folds co-transcriptionally and must refold following the passage of ribosome. The mRNAs can be entrapped in metastable structures during these folding events. In the present study, we evaluated the conformational dynamics of the kinetically favored, metastable, and hairpin-like structure, which disturbs the thermodynamically favored G-quadruplex structure, and its effect on co-transcriptional translation in prokaryotic cells. We found that nascent mRNA forms a metastable hairpin-like structure during co-transcriptional folding instead of the G-quadruplex structure. When the translation progressed co-transcriptionally before the metastable hairpin-like structure transition to the G-quadruplex, function of the G-quadruplex as a roadblock of the ribosome was sequestered. This suggested that kinetically formed RNA structures had a dominant effect on gene expression in prokaryotes. The results of this study indicate that it is critical to consider the conformational dynamics of RNA-folding to understand the contributions of the mRNA structures in controlling gene expression.

RNA performs various spatiotemporal functions in living cells. As the solution environments significantly affect the stability of RNA duplexes, a stability prediction of the RNA duplexes in diverse crowded conditions is required to understand and modulate gene expression in heterogeneously crowded intracellular conditions. Herein, we determined the nearest-neighbor (NN) parameters for RNA duplex formation when subjected to crowding conditions with an ionic concentration relevant to that found in cells. Determination of the individual contributions of excluded volume effect and water activity to each of the NN parameters in crowded environments enabled prediction of the thermodynamic parameters and their melting temperatures for plenty of tested RNA duplex formation in vitro and in cell with significant accuracy. The parameters reported herein will help predicting RNA duplex stability in different crowded environments, which will lead to an improved understanding of the stability-function relationship for RNAs in various cellular organelles with different molecular environments.

Climate change negatively affects crop yield, which hinders efforts to reach agricultural sustainability and food security. Here, we show that a previously unidentified allele of the nitrate transporter gene OsNRT2.3 is required to maintain high yield and high nitrogen use efficiency under high temperatures. We demonstrate that this tolerance to high temperatures in rice accessions harboring the HTNE-2 (high temperature resistant and nitrogen efficient-2) alleles from enhanced translation of the OsNRT2.3b mRNA isoform and the decreased abundance of a unique small RNA (sNRT2.3-1) derived from the 5′ untranslated region of OsNRT2.3 . sNRT2.3-1 binds to the OsNRT2.3a mRNA in a temperature-dependent manner. Our findings reveal that allelic variation in the 5′ untranslated region of OsNRT2.3 leads to an increase in OsNRT2.3b protein levels and higher yield during high-temperature stress. Our results also provide a breeding strategy to produce rice varieties with higher grain yield and lower N fertilizer input suitable for a sustainable agriculture that is resilient against climate change.