Modeling Structure, Stability, and Flexibility of Double-Stranded RNAs in Salt Solutions

Supercharged coiled-coil protein with N-terminal decahistidine tag boosts siRNA complexation and delivery efficiency of a lipoproteoplex

Short interfering RNA (siRNA) therapeutics have soared in popularity due to their highly selective and potent targeting of faulty genes, providing a non-palliative approach to address diseases. Despite their potential, effective transfection of siRNA into cells requires the assistance of an accompanying vector. Vectors constructed from non-viral materials, while offering safer and non-cytotoxic profiles, often grapple with lackluster loading and delivery efficiencies, necessitating substantial milligram quantities of expensive siRNA to confer the desired downstream effects. We detail the recombinant synthesis of a diverse series of coiled-coil supercharged protein (CSP) biomaterials systematically designed to investigate the impact of two arginine point mutations (Q39R and N61R) and decahistidine tags on liposomal siRNA delivery. The most efficacious variant, N8, exhibits a twofold increase in its affinity to siRNA and achieves a twofold enhancement in transfection activity with minimal cytotoxicity in vitro. Subsequent analysis unveils the destabilizing effect of the Q39R and N61R supercharging mutations and the incorporation of C-terminal decahistidine tags on α-helical secondary structure. Cross-correlational regression analyses reveal that the amount of helical character in these mutants is key in N8's enhanced siRNA complexation and downstream delivery efficiency.

Atomistic simulations of RNA duplex thermal denaturation: Sequence- and forcefield-dependence

Double-stranded RNA is the end-product of template-based replication, and is also the functional state of some biological RNAs. Similarly to proteins and DNA, they can be denatured by temperature, with important physiological and technological implications. Here, we use an in silico strategy to probe the thermal denaturation of RNA duplexes. Following previous results that were obtained on a few different duplexes, and which nuanced the canonical 2-state picture of nucleic acid denaturation, we here specifically address three different aspects that greatly improve our description of the temperature-induced dsRNA separation. First, we investigate the effect of the spatial distribution of weak and strong base-pairs among the duplex sequence. We show that the deviations from the two-state dehybridization mechanism are more pronounced when a strong core is flanked with weak extremities, while duplexes with a weak core but strong extremities exhibit a two-state behavior, which can be explained by the key role played by base fraying. This was later verified by generating artificial hairpin or circular states containing one or two locked duplex extremities, which results in an important reinforcement of the entire HB structure of the duplex and higher melting temperatures. Finally, we demonstrate that our results are little sensitive to the employed combination of RNA and water forcefields. The trends in thermal stability among the different sequences as well as the observed unfolding mechanisms (and the deviations from a two-state scenario) remain the same regardless of the employed atomistic models. However, our study points to possible limitations of recent reparametrizations of the Amber RNA forcefield, which sometimes results in duplexes that readily denature under ambient conditions, in contradiction with available experimental results.

RNA 3D Structure Prediction: Progress and Perspective

Ribonucleic acid (RNA) molecules play vital roles in numerous important biological functions such as catalysis and gene regulation. The functions of RNAs are strongly coupled to their structures or proper structure changes, and RNA structure prediction has been paid much attention in the last two decades. Some computational models have been developed to predict RNA three-dimensional (3D) structures in silico, and these models are generally composed of predicting RNA 3D structure ensemble, evaluating near-native RNAs from the structure ensemble, and refining the identified RNAs. In this review, we will make a comprehensive overview of the recent advances in RNA 3D structure modeling, including structure ensemble prediction, evaluation, and refinement. Finally, we will emphasize some insights and perspectives in modeling RNA 3D structures.

Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding

DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.

Predicting 3D structures and stabilities for complex RNA pseudoknots in ion solutions

RNA pseudoknots are a kind of important tertiary motif, and the structures and stabilities of pseudoknots are generally critical to the biological functions of RNAs with the motifs. In this work, we have carefully refined our previously developed coarse-grained model with salt effect through involving a new coarse-grained force field and a replica-exchange Monte Carlo algorithm, and employed the model to predict structures and stabilities of complex RNA pseudoknots in ion solutions beyond minimal H-type pseudoknots. Compared with available experimental data, the newly refined model can successfully predict 3D structures from sequences for the complex RNA pseudoknots including SARS-CoV-2 programming-1 ribosomal frameshifting element and Zika virus xrRNA, and can reliably predict the thermal stabilities of RNA pseudoknots with various sequences and lengths over broad ranges of monovalent/divalent salts. In addition, for complex pseudoknots including SARS-CoV-2 frameshifting element, our analyses show that their thermally unfolding pathways are mainly dependent on the relative stabilities of unfolded intermediate states, in analogy to those of minimal H-type pseudoknots.

cgRNASP: coarse-grained statistical potentials with residue separation for RNA structure evaluation

Knowledge-based statistical potentials are very important for RNA 3-dimensional (3D) structure prediction and evaluation. In recent years, various coarse-grained (CG) and all-atom models have been developed for predicting RNA 3D structures, while there is still lack of reliable CG statistical potentials not only for CG structure evaluation but also for all-atom structure evaluation at high efficiency. In this work, we have developed a series of residue-separation-based CG statistical potentials at different CG levels for RNA 3D structure evaluation, namely cgRNASP, which is composed of long-ranged and short-ranged interactions by residue separation. Compared with the newly developed all-atom rsRNASP, the short-ranged interaction in cgRNASP was involved more subtly and completely. Our examinations show that, the performance of cgRNASP varies with CG levels and compared with rsRNASP, cgRNASP has similarly good performance for extensive types of test datasets and can have slightly better performance for the realistic dataset-RNA-Puzzles dataset. Furthermore, cgRNASP is strikingly more efficient than all-atom statistical potentials/scoring functions, and can be apparently superior to other all-atom statistical potentials and scoring functions trained from neural networks for the RNA-Puzzles dataset. cgRNASP is available at https://github.com/Tan-group/cgRNASP.

Ab initio predictions for 3D structure and stability of single- and double-stranded DNAs in ion solutions

The three-dimensional (3D) structure and stability of DNA are essential to understand/control their biological functions and aid the development of novel materials. In this work, we present a coarse-grained (CG) model for DNA based on the RNA CG model proposed by us, to predict 3D structures and stability for both dsDNA and ssDNA from the sequence. Combined with a Monte Carlo simulated annealing algorithm and CG force fields involving the sequence-dependent base-pairing/stacking interactions and an implicit electrostatic potential, the present model successfully folds 20 dsDNAs (≤52nt) and 20 ssDNAs (≤74nt) into the corresponding native-like structures just from their sequences, with an overall mean RMSD of 3.4Å from the experimental structures. For DNAs with various lengths and sequences, the present model can make reliable predictions on stability, e.g., for 27 dsDNAs with/without bulge/internal loops and 24 ssDNAs including pseudoknot, the mean deviation of predicted melting temperatures from the corresponding experimental data is only ~2.0°C. Furthermore, the model also quantificationally predicts the effects of monovalent or divalent ions on the structure stability of ssDNAs/dsDNAs.

To determine 3D structures and quantify stability of single- (ss) and double-stranded (ds) DNAs is essential to unveil the mechanisms of their functions and to further guide the production and development of novel materials. Although many DNA models have been proposed to reproduce the basic structural, mechanical, or thermodynamic properties of dsDNAs based on the secondary structure information or preset constraints, there are very few models can be used to investigate the ssDNA folding or dsDNA assembly from the sequence. Furthermore, due to the polyanionic nature of DNAs, metal ions (e.g., Na⁺ and Mg²⁺) in solutions can play an essential role in DNA folding and dynamics. Nevertheless, ab initio predictions for DNA folding in ion solutions are still an unresolved problem. In this work, we developed a novel coarse-grained model to predict 3D structures and thermodynamic stabilities for both ssDNAs and dsDNAs in monovalent/divalent ion solutions from their sequences. As compared with the extensive experimental data and available existing models, we showed that the present model can successfully fold simple DNAs into their native-like structures, and can also accurately reproduce the effects of sequence and monovalent/divalent ions on structure stability for ssDNAs including pseudoknot and dsDNAs with/without bulge/internal loops.

FebRNA: An automated fragment-ensemble-based model for building RNA 3D structures

Knowledge of RNA three-dimensional (3D) structures is critical to understanding the important biological functions of RNAs. Although various structure prediction models have been developed, the high-accuracy predictions of RNA 3D structures are still limited to the RNAs with short lengths or with simple topology. In this work, we proposed a new model, namely FebRNA, for building RNA 3D structures through fragment assembly based on coarse-grained (CG) fragment ensembles. Specifically, FebRNA is composed of four processes: establishing the library of different types of non-redundant CG fragment ensembles regardless of the sequences, building CG 3D structure ensemble through fragment assembly, identifying top-scored CG structures through a specific CG scoring function, and rebuilding the all-atom structures from the top-scored CG ones. Extensive examination against different types of RNA structures indicates that FebRNA consistently gives the reliable predictions on RNA 3D structures, including pseudoknots, three-way junctions, four-way and five-way junctions, and RNAs in the RNA-Puzzles. FebRNA is available on the Web site: https://github.com/Tan-group/FebRNA.

Stability Modelling of mRNA Vaccine Quality Based on Temperature Monitoring throughout the Distribution Chain

The vaccine distribution chains in several low- and middle-income countries are not adequate to facilitate the rapid delivery of high volumes of thermosensitive COVID-19 mRNA vaccines at the required low and ultra-low temperatures. COVID-19 mRNA vaccines are currently distributed along with temperature monitoring devices to track and identify deviations from predefined conditions throughout the distribution chain. These temperature readings can feed into computational models to quantify mRNA vaccine critical quality attributes (CQAs) and the remaining vaccine shelf life more accurately. Here, a kinetic modelling approach is proposed to quantify the stability-related CQAs and the remaining shelf life of mRNA vaccines. The CQA and shelf-life values can be computed based on the conditions under which the vaccines have been distributed from the manufacturing facilities via the distribution network to the vaccination centres. This approach helps to quantify the degree to which temperature excursions impact vaccine quality and can also reduce vaccine wastage. In addition, vaccine stock management can be improved due to the information obtained on the remaining shelf life of mRNA vaccines. This model-based quantification of mRNA vaccine quality and remaining shelf life can improve the deployment of COVID-19 mRNA vaccines to low- and middle-income countries.

rsRNASP: A residue-separation-based statistical potential for RNA 3D structure evaluation

Knowledge-based statistical potentials have been shown to be rather effective in protein 3-dimensional (3D) structure evaluation and prediction. Recently, several statistical potentials have been developed for RNA 3D structure evaluation, while their performances are either still at a low level for the test datasets from structure prediction models or dependent on the "black-box" process through neural networks. In this work, we have developed an all-atom distance-dependent statistical potential based on residue separation for RNA 3D structure evaluation, namely rsRNASP, which is composed of short- and long-ranged potentials distinguished by residue separation. The extensive examinations against available RNA test datasets show that rsRNASP has apparently higher performance than the existing statistical potentials for the realistic test datasets with large RNAs from structure prediction models, including the newly released RNA-Puzzles dataset, and is comparable to the existing top statistical potentials for the test datasets with small RNAs or near-native decoys. In addition, rsRNASP is superior to RNA3DCNN, a recently developed scoring function through 3D convolutional neural networks. rsRNASP and the relevant databases are available to the public.

Salt-Dependent RNA Pseudoknot Stability: Effect of Spatial Confinement

Macromolecules, such as RNAs, reside in crowded cell environments, which could strongly affect the folded structures and stability of RNAs. The emergence of RNA-driven phase separation in biology further stresses the potential functional roles of molecular crowding. In this work, we employed the coarse-grained model that was previously developed by us to predict 3D structures and stability of the mouse mammary tumor virus (MMTV) pseudoknot under different spatial confinements over a wide range of salt concentrations. The results show that spatial confinements can not only enhance the compactness and stability of MMTV pseudoknot structures but also weaken the dependence of the RNA structure compactness and stability on salt concentration. Based on our microscopic analyses, we found that the effect of spatial confinement on the salt-dependent RNA pseudoknot stability mainly comes through the spatial suppression of extended conformations, which are prevalent in the partially/fully unfolded states, especially at low ion concentrations. Furthermore, our comprehensive analyses revealed that the thermally unfolding pathway of the pseudoknot can be significantly modulated by spatial confinements, since the intermediate states with more extended conformations would loss favor when spatial confinements are introduced.

Salt dependent mesoscopic model for RNA at multiple strand concentrations

Mesoscopic models can be used for the description of the thermodynamic properties of RNA duplexes. With the use of experimental melting temperatures, its parametrization can provide important insights into its hydrogen bonds and stacking interactions as has been done for high sodium concentrations. However, the RNA parametrization for lower salt concentrations is still missing due to the limited amount of published melting temperature data. While the Peyrard-Bishop (PB) parametrization was found to be largely independent of strand concentrations, it requires that all temperatures are provided at the same strand concentrations. Here we adapted the PB model to handle multiple strand concentrations and in this way we were able to make use of an experimental set of temperatures to model the hydrogen bond and stacking interactions at low and intermediate sodium concentrations. For the parametrizations we make a distinction between terminal and internal base pairs, and the resulting potentials were qualitatively similar as we obtained previously for DNA. The main difference from DNA parameters, was the Morse potentials at low sodium concentrations for terminal r(AU) which is stronger than d(AT), suggesting higher hydrogen bond strength.

Ion-mediated interactions between like-charged polyelectrolytes with bending flexibility

Ion-mediated interactions between polyelectrolytes (PEs) are crucial to the properties of flexible biopolymers such as nucleic acids and proteins but the effect of PE flexibility on such interactions has not been explicitly addressed until now. In this work, the potentials of mean force (PMFs) between like-charged PEs with different bending flexibility have been investigated by Monte Carlo simulations and a cylindrical confinement around each PE was involved to model two PEs in an array. We found that in the absence of trivalent salt, the PMFs between like-charged PEs in an array are apparently repulsive while the bending flexibility can visibly decrease the repulsive PMFs. With the addition of high trivalent salt, the PMFs become significantly attractive whereas the attractive PMFs can be apparently weakened by the bending flexibility. Our analyses reveal that the effect of bending flexibility is attributed to the increased PE conformational space, which allows the PEs to fluctuate away to decrease the monovalent ion-mediated repulsion or to weaken the trivalent ion-mediated attraction through disrupting trivalent ion-bridging configuration. Additionally, our further calculations show that the effect of bending flexibility on the ion-mediated interactions is less apparent for PEs without cylindrical confinement.

Salt effect on thermodynamics and kinetics of a single RNA base pair

Due to the polyanionic nature of RNAs, the structural folding of RNAs are sensitive to solution salt conditions, while there is still lack of a deep understanding of the salt effect on the thermodynamics and kinetics of RNAs at a single base-pair level. In this work, the thermodynamic and the kinetic parameters for the base-pair AU closing/opening at different salt concentrations were calculated by 3-µsec all-atom molecular dynamics (MD) simulations at different temperatures. It was found that for the base-pair formation, the enthalpy change [Formula: see text] is nearly independent of salt concentration, while the entropy change [Formula: see text] exhibits a linear dependence on the logarithm of salt concentration, verifying the empirical assumption based on thermodynamic experiments. Our analyses revealed that such salt concentration dependence of the entropy change mainly results from the dependence of ion translational entropy change for the base pair closing/opening on salt concentration. Furthermore, the closing rate increases with the increasing of salt concentration, while the opening rate is nearly independent of salt concentration. Additionally, our analyses revealed that the free energy surface for describing the base-pair opening and closing dynamics becomes more rugged with the decrease of salt concentration.

Salt Dependence of A-Form RNA Duplexes: Structures and Implications

The biological functions of RNA range from gene regulation through catalysis and depend critically on its structure and flexibility. Conformational variations of flexible, non-base-paired components, including RNA hinges, bulges, or single-stranded tails, are well documented. Recent work has also identified variations in the structure of ubiquitous, base-paired duplexes found in almost all functional RNAs. Duplexes anchor the structures of folded RNAs, and their surface features are recognized by partner molecules. To date, no consistent picture has been obtained that describes the range of conformations assumed by RNA duplexes. Here, we apply wide angle, solution X-ray scattering (WAXS) to quantify these variations, by sampling length scales characteristic of helical geometries under different solution conditions. To identify the radius, helical rise, twist, and length of dsRNA helices, we exploit molecular dynamics generated structures, explicit solvent models, and ensemble optimization methods. Our results quantify the substantial and salt-dependent deviations of double-stranded (ds) RNA duplexes from the assumed canonical A-form conformation. Recent experiments underscore the need to properly describe the structures of RNA duplexes when interpreting the salt dependence of RNA conformations.

Structure folding of RNA kissing complexes in salt solutions: predicting 3D structure, stability, and folding pathway

RNA kissing complexes are essential for genomic RNA dimerization and regulation of gene expression, and their structures and stability are critical to their biological functions. In this work, we used our previously developed coarse-grained model with an implicit structure-based electrostatic potential to predict three-dimensional (3D) structures and stability of RNA kissing complexes in salt solutions. For extensive RNA kissing complexes, our model shows great reliability in predicting 3D structures from their sequences, and our additional predictions indicate that the model can capture the dependence of 3D structures of RNA kissing complexes on monovalent/divalent ion concentrations. Moreover, the comparisons with extensive experimental data show that the model can make reliable predictions on the stability for various RNA kissing complexes over wide ranges of monovalent/divalent ion concentrations. Notably, for RNA kissing complexes, our further analyses show the important contribution of coaxial stacking to the 3D structures and stronger stability than the corresponding kissing-interface duplexes at high salts. Furthermore, our comprehensive analyses for RNA kissing complexes reveal that the thermally folding pathway for a complex sequence is mainly determined by the relative stability of two possible folded states of kissing complex and extended duplex, which can be significantly modulated by its sequence.

3D structure stability of the HIV-1 TAR RNA in ion solutions: A coarse-grained model study

As an extremely common structural motif, RNA hairpins with bulge loops [e.g., the human immunodeficiency virus type 1 (HIV-1) transactivation response (TAR) RNA] can play essential roles in normal cellular processes by binding to proteins and small ligands, which could be very dependent on their three-dimensional (3D) structures and stability. Although the structures and conformational dynamics of the HIV-1 TAR RNA have been extensively studied, there are few investigations on the thermodynamic stability of the TAR RNA, especially in ion solutions, and the existing studies also have some divergence on the unfolding process of the RNA. Here, we employed our previously developed coarse-grained model with implicit salt to predict the 3D structure, stability, and unfolding pathway for the HIV-1 TAR RNA over a wide range of ion concentrations. As compared with the extensive experimental/theoretical results, the present model can give reliable predictions on the 3D structure stability of the TAR RNA from the sequence. Based on the predictions, our further comprehensive analyses on the stability of the TAR RNA as well as its variants revealed that the unfolding pathway of an RNA hairpin with a bulge loop is mainly determined by the relative stability between different states (folded state, intermediate state, and unfolded state) and the strength of the coaxial stacking between two stems in folded structures, both of which can be apparently modulated by the ion concentrations as well as the sequences.

What is the best reference state for building statistical potentials in RNA 3D structure evaluation?

Knowledge-based statistical potentials have been shown to be efficient in protein structure evaluation/prediction, and the core difference between various statistical potentials is attributed to the choice of reference states. However, for RNA 3D structure evaluation, a comprehensive examination on reference states is still lacking. In this work, we built six statistical potentials based on six reference states widely used in protein structure evaluation, including averaging, quasi-chemical approximation, atom-shuffled, finite-ideal-gas, spherical-noninteracting, and random-walk-chain reference states, and we examined the six reference states against three RNA test sets including six subsets. Our extensive examinations show that, overall, for identifying native structures and ranking decoy structures, the finite-ideal-gas and random-walk-chain reference states are slightly superior to others, while for identifying near-native structures, there is only a slight difference between these reference states. Our further analyses show that the performance of a statistical potential is apparently dependent on the quality of the training set. Furthermore, we found that the performance of a statistical potential is closely related to the origin of test sets, and for the three realistic test subsets, the six statistical potentials have overall unsatisfactory performance. This work presents a comprehensive examination on the existing reference states and statistical potentials for RNA 3D structure evaluation.

Structural Flexibility of DNA-RNA Hybrid Duplex: Stretching and Twist-Stretch Coupling

DNA-RNA hybrid (DRH) duplexes play essential roles during the replication of DNA and the reverse transcription of RNA viruses, and their flexibility is important for their biological functions. Recent experiments indicated that A-form RNA and B-form DNA have a strikingly different flexibility in stretching and twist-stretch coupling, and the structural flexibility of DRH duplex is of great interest, especially in stretching and twist-stretch coupling. In this work, we performed microsecond all-atom molecular dynamics simulations with new AMBER force fields to characterize the structural flexibility of DRH duplex in stretching and twist-stretch coupling. We have calculated all the helical parameters, stretch modulus, and twist-stretch coupling parameters for the DRH duplex. First, our analyses on structure suggest that the DRH duplex exhibits an intermediate conformation between A- and B-forms and closer to A-form, which can be attributed to the stronger rigidity of the RNA strand than the DNA strand. Second, our calculations show that the DRH duplex has the stretch modulus of 834 ± 34 pN and a very weak twist-stretch coupling. Our quantitative analyses indicate that, compared with DNA and RNA duplexes, the different flexibility of the DRH duplex in stretching and twist-stretch coupling is mainly attributed to its apparently different basepair inclination in the helical structure.