Effect of ethanol on the elasticities of double-stranded RNA and DNA revealed by magnetic tweezers and simulations.

The elasticities of double-stranded (ds) DNA and RNA, which are critical to their biological functions and applications in materials science, can be significantly modulated by solution conditions such as ions and temperature. However, there is still a lack of a comprehensive understanding of the role of solvents in the elasticities of dsRNA and dsDNA in a comparative way. In this work, we explored the effect of ethanol solvent on the elasticities of dsRNA and dsDNA by magnetic tweezers and all-atom molecular dynamics simulations. We found that the bending persistence lengths and contour lengths of dsRNA and dsDNA decrease monotonically with the increase in ethanol concentration. Furthermore, the addition of ethanol weakens the positive twist-stretch coupling of dsRNA, while promotes the negative twist-stretch coupling of dsDNA. Counter-intuitively, the lower dielectric environment of ethanol causes a significant re-distribution of counterions and enhanced ion neutralization, which overwhelms the enhanced repulsion along dsRNA/dsDNA, ultimately leading to the softening in bending for dsRNA and dsDNA. Moreover, for dsRNA, ethanol causes slight ion-clamping across the major groove, which weakens the major groove-mediated twist-stretch coupling, while for dsDNA, ethanol promotes the stretch-radius correlation due to enhanced ion binding and consequently enhances the helical radius-mediated twist-stretch coupling.

The origin of different bending stiffness between double-stranded RNA and DNA revealed by magnetic tweezers and simulations.

The subtle differences in the chemical structures of double-stranded (ds) RNA and DNA lead to significant variations in their biological roles and medical implications, largely due to their distinct biophysical properties, such as bending stiffness. Although it is well known that A-form dsRNA is stiffer than B-form dsDNA under physiological salt conditions, the underlying cause of this difference remains unclear. In this study, we employ high-precision magnetic-tweezer experiments along with molecular dynamics simulations and reveal that the relative bending stiffness between dsRNA and dsDNA is primarily determined by the structure- and salt-concentration-dependent ion distribution around their helical structures. At near-physiological salt conditions, dsRNA shows a sparser ion distribution surrounding its phosphate groups compared to dsDNA, causing its greater stiffness. However, at very high monovalent salt concentrations, phosphate groups in both dsRNA and dsDNA become fully neutralized by excess ions, resulting in a similar intrinsic bending persistence length of approximately 39 nm. This similarity in intrinsic bending stiffness of dsRNA and dsDNA is coupled to the analogous fluctuations in their total groove widths and further coupled to the similar fluctuation of base-pair inclination, despite their distinct A-form and B-form helical structures.

Multifunctional AuPt Nanoparticles for Synergistic Photothermal and Radiation Therapy.

Background: Photothermal therapy (PTT) has gained considerable interest as an emerging modality for cancer treatment in recent years. Radiation therapy (RT) has been widely used in the clinic as a traditional treatment method. However, RT and PTT treatments are limited by side effects and penetration depth, respectively. In addition, hypoxia within the tumor can lead to increased resistance to treatment. Methods: We synthesized multiple sizes of AuPt by modulating the reaction conditions. The smallest size of AuPt was selected and modified with folic acid (FA) for PTT and RT synergy therapy. Various methods including transmission electron microscope (TEM), X-ray photoelectron spectroscopy (XPS), X-ray diffraction (XRD), and Fourier transform infrared spectroscopy (FITR) are used to determine the structure and composition of AuPt-FA (AF). In addition, we researched the photothermal properties of AF with IR cameras and infrared lasers. Flow cytometry, colony formation assays, CCK8, and fluorescent staining for probing the treatment effect in vitro. Also, we explored the targeting of AF by TEM and In Vivo Imaging Systems (IVIS). In vivo experiments, we record changes in tumor volume and weight as well as staining of tumor sections (ROS, Ki67, and hematoxylin and eosin). Results: The AuPt with particle size of 16 nm endows it with remarkably high photothermal conversion efficiency (46.84%) and catalase activity compared to other sizes of AuPt (30 nm and 100 nm). AF alleviates hypoxia in the tumor microenvironment, leading to the production of more reactive oxygen species (ROS) during the treatment. In addition, the therapeutic effect was significantly enhanced by combining RT and PTT, with an apoptosis rate of 81.1% in vitro and an in vivo tumor volume reduction rate of 94.0% in vivo. Conclusion: These results demonstrate that AF potentiates the synergistic effect of PTT and RT and has the potential for clinical translation.

Modeling Coil-Globule-Helix Transition in Polymers by Self-Interacting Random Walks.

Random walks (RWs) have been important in statistical physics and can describe the statistical properties of various processes in physical, chemical, and biological systems. In this study, we have proposed a self-interacting random walk model in a continuous three-dimensional space, where the walker and its previous visits interact according to a realistic Lennard-Jones (LJ) potential uLJr=εr0/r12-2r0/r6. It is revealed that the model shows a novel globule-to-helix transition in addition to the well-known coil-to-globule collapse in its trajectory when the temperature decreases. The dependence of the structural transitions on the equilibrium distance r0 of the LJ potential and the temperature T were extensively investigated. The system showed many different structural properties, including globule-coil, helix-globule-coil, and line-coil transitions depending on the equilibrium distance r0 when the temperature T increases from low to high. We also obtained a correlation form of kBTc = λε for the relationship between the transition temperature Tc and the well depth ε, which is consistent with our numerical simulations. The implications of the random walk model on protein folding are also discussed. The present model provides a new way towards understanding the mechanism of helix formation in polymers like proteins.

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.

Universality in RNA and DNA deformations induced by salt, temperature change, stretching force, and protein binding.

Nucleic acid deformations play important roles in many biological processes. The physical understanding of nucleic acid deformation by environmental stimuli is limited due to the challenge in the precise measurement of RNA and DNA deformations and the complexity of interactions in RNA and DNA. Magnetic tweezers experiments provide an excellent opportunity to precisely measure DNA and RNA twist changes induced by environmental stimuli. In this work, we applied magnetic tweezers to measure double-stranded RNA twist changes induced by salt and temperature changes. We observed RNA unwinds when lowering salt concentration, or increasing temperature. Our molecular dynamics simulations revealed the mechanism: lowering salt concentration or increasing temperature enlarges RNA major groove width, which causes twist decrease through twist-groove coupling. Combining these results with previous results, we found some universality in RNA and DNA deformations induced by three different stimuli: salt change, temperature, and stretching force. For RNA, these stimuli first modify the major groove width, which is transduced into twist change through twist-groove coupling. For DNA, these stimuli first modify diameter, which is transduced into twist change through twist-diameter coupling. Twist-groove coupling and twist-diameter coupling appear to be utilized by protein binding to reduce DNA and RNA deformation energy cost upon protein binding.

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.

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.

5-Methyl-cytosine stabilizes DNA but hinders DNA hybridization revealed by magnetic tweezers and simulations.

5-Methyl-cytosine (5mC) is one of the most important DNA modifications and plays versatile biological roles. It is well known that 5mC stabilizes DNA duplexes. However, it remains unclear how 5mC affects the kinetics of DNA melting and hybridization. Here, we studied the kinetics of unzipping and rezipping using a 502-bp DNA hairpin by single-molecule magnetic tweezers. Under constant loading rates, 5mC increases the unzipping force but counterintuitively decreases the rezipping force at various salt and temperature conditions. Under constant forces, the non-methylated DNA hops between metastable states during unzipping and rezipping, which implies low energy barriers. Surprisingly, the 5mC DNA can't rezip after fully unzipping unless much lower forces are applied, where it rezips stochastically in a one-step manner, which implies 5mC kinetically hinders DNA hybridization and high energy barriers in DNA hybridization. All-atom molecular dynamics simulations reveal that the 5mC kinetically hinders DNA hybridization due to steric effects rather than electrostatic effects caused by the additional methyl groups of cytosines. Considering the possible high speed of DNA unzipping and zipping during replication and transcription, our findings provide new insights into the biological roles of 5mC.

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.

Twist-diameter coupling drives DNA twist changes with salt and temperature.

DNA deformations upon environmental changes, e.g., salt and temperature, play crucial roles in many biological processes and material applications. Here, our magnetic tweezers experiments observed that the increase in NaCl, KCl, or RbCl concentration leads to substantial DNA overwinding. Our simulations and theoretical calculation quantitatively explain the salt-induced twist change through the mechanism: More salt enhances the screening of interstrand electrostatic repulsion and hence reduces DNA diameter, which is transduced to twist increase through twist-diameter coupling. We determined that the coupling constant is 4.5 ± 0.8 kBT/(degrees∙nm) for one base pair. The coupling comes from the restraint of the contour length of DNA backbone. On the basis of this coupling constant and diameter-dependent DNA conformational entropy, we predict the temperature dependence of DNA twist Δωbp/ΔT ≈ -0.01 degree/°C, which agrees with our and previous experimental results. Our analysis suggests that twist-diameter coupling is a common driving force for salt- and temperature-induced DNA twist changes.

Multivalent Cations Reverse the Twist-Stretch Coupling of RNA.

When stretched, both DNA and RNA duplexes change their twist angles through twist-stretch coupling. The coupling is negative for DNA but positive for RNA, which is not yet completely understood. Here, our magnetic tweezers experiments show that the coupling of RNA reverses from positive to negative by multivalent cations. Combining with the previously reported tension-induced negative-to-positive coupling reversal of DNA, we propose a unified mechanism of the couplings of both RNA and DNA based on molecular dynamics simulations. Two deformation pathways are competing when stretched: shrinking the radius causes positive couplings but widening the major groove causes negative couplings. For RNA whose major groove is clamped by multivalent cations and canonical DNA, their radii shrink when stretched, thus exhibiting positive couplings. For elongated DNA whose radius already shrinks to the minimum and canonical RNA, their major grooves are widened when stretched, thus exhibiting negative couplings.

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.

An intratumoral injectable nanozyme hydrogel for hypoxia-resistant thermoradiotherapy.

Hypoxia in local tumors leads to the failure or resistance of radiotherapy (RT) and high-dose RT will cause systemic reactions and local radiation damage. As a non-chemotherapeutic intervention, photothermal therapy (PTT) can remove tumor tissues through thermal ablation as well as effectively improve the microenvironment of hypoxic cells. Therefore, the combined use of PTT and RT (thermoradiotherapy) has urgently become an efficient treatment. In this work, by encapsulating prussian blue (PB) nanoparticles in agarose hydrogel, we developed an injectable hybrid light-controlled hydrogel system as a PB reservoir and release controller (PRC) which can realize single injection and multiple treatments in vivo. Under the irradiation of 808 nm near-infrared (NIR) laser, PB nanoparticles convert laser energy into heat energy, causing degradation of agarose hydrogel and the release of PB nanoparticles. Due to the excellent photothermal properties of PB, photothermal treatment in the NIR Biological Windows can greatly enhance the sensitivity of tumor cells to RT. Meanwhile, PB nanoparticles can also be a nanozyme to drive the decomposition of endogenous hydrogen peroxide (H2O2), and then generate oxygen (O2) to improve the tumor hypoxic microenvironment, achieving the further enhancement of the radiation sensitivity. Notably, this study is the first design to utilize hydrogel for thermoradiotherapy. Both in vitro and in vivo experiments, the PRC demonstrated excellent effects of PTT-RT, good stability and biocompatibility, indicating our nanoplatform promote the development of anti-cancer combination thermoradiotherapy with greater clinical significance.

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.

Multivalent Ion-Mediated Attraction between Like-Charged Colloidal Particles: Nonmonotonic Dependence on the Particle Charge.

Ion-mediated effective interactions are important for the structure and stability of charged particles such as colloids and nucleic acids. It has been known that the intrinsic electrostatic repulsion between like-charged particles can be modulated into effective attraction by multivalent ions. In this work, we examined the dependence of multivalent ion-mediated attraction between like-charged colloidal particles on the particle charge in a wide range by extensive Monte Carlo simulations. Our calculations show that for both divalent and trivalent salts, the effective attraction between like-charged colloidal particles becomes stronger with the increase of the particle charge, whereas it gradually becomes weakened when the particle charge exceeds a "critical" value. Correspondingly, as the particle charge is increased, the driving force for such effective attraction transits from an attractive electrostatic force to an attractive depletion force, and the attraction weakening by high particle charges is attributed to the transition of electrostatic force from attraction to repulsion. Our analyses suggest that the attractive depletion force and the repulsive electrostatic force at high particle charges result from the Coulomb depletion which suppresses the counterion condensation in the limited region between two like-charged colloidal particles. Moreover, our extensive calculations indicate that the "critical" particle charge decreases apparently for larger ions and smaller colloidal particles due to stronger Coulomb depletion and decreases slightly at higher salt concentrations due to the slightly enhanced Coulomb depletion in the intervening space between colloidal particles. Encouragingly, we derived an analytical formula for the "critical" particle charge based on the Lindemann melting law.

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.

Cytosine Methylation Enhances DNA Condensation Revealed by Equilibrium Measurements Using Magnetic Tweezers.

CpG methylation of DNA is common in mammalian cells. In sperm, the DNA has the highest level of CpG methylation and is condensed into toroidal structures. How CpG methylation affects DNA structures and interactions is important to understand its biological roles but is largely unknown. Using an RNA-DNA-RNA structure, we observed the equilibrium hopping dynamics between the condensed and extended states of DNA in the presence of polyamines or polylysine peptide as a reduced model of histone tails. Combing with the measured DNA elasticities, we report that CpG methylation of each cytosine nucleotide substantially increases DNA-DNA attraction by up to 0.2 kBT. For the DNA with 57% GC content, the relative increase caused by CpG methylation is up to 32% for the spermine-induced DNA-DNA attraction and up to 9% for the polylysine-induced DNA-DNA attraction. These findings help us to evaluate the energetic contributions of CpG methylation in sperm development and chromatin regulation.

Opposite Effects of High-Valent Cations on the Elasticities of DNA and RNA Duplexes Revealed by Magnetic Tweezers.

We report that trivalent cobalt hexammine cations decrease the persistence length, stretching modulus, helical density, and size of plectonemes formed under torque of DNA but increase those of RNA. Divalent magnesium cations, however, decrease the persistence lengths, contour lengths, and sizes of plectonemes while increasing the helical densities of both DNA and RNA. The experimental results are explained by different binding modes of the cations on DNA and RNA in our all-atom molecular dynamics simulations. The significant variations of the helical densities and structures of DNA and RNA duplexes induced by high-valent cations may affect interactions of the duplexes with proteins.

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.