Small Details Matter: The 2′-Hydroxyl as a Conformational Switch in RNA

modXNA: A Modular Approach to Parametrization of Modified Nucleic Acids for Use with Amber Force Fields

Modified nucleic acids have surged as a popular therapeutic route, emphasizing the importance of nucleic acid research in drug discovery and development. Beyond well-known RNA vaccines, antisense oligonucleotides and aptamers can incorporate various modified nucleic acids to target specific biomolecules for various therapeutic activities. Molecular dynamics simulations can accelerate the design and development of these systems with noncanonical nucleic acids by observing intricate dynamic properties and relative stability on the all-atom level. However, modeling these modified systems is challenging due to the time and resources required to parametrize components outside default force field parameters. Here, we present modXNA, a tool to derive and build modified nucleotides for use with Amber force fields. Several nucleic acid systems varying in size and number of modification sites were used to evaluate the accuracy of modXNA parameters, and results indicate the dynamics and structure are preserved throughout the simulations. We detail the protocol for quantum mechanics charge derivation and describe a workflow for implementing modXNA in Amber molecular dynamics simulations, which includes updates and added features to CPPTRAJ.

Hydration effect and molecular geometry conformation as critical factors affecting the longevity stability of G4-structure-based supramolecular hydrogels

Nucleoside-derived supramolecular hydrogels based on G4-structures have been extensively developed in the biomedical sector and recognized for superior excellent biocompatibility and biodegradability. However, limited longevity and stability present a significant challenge. Chemical modifications in the molecular structure have been shown to enhance the longevity stability of G4-structure-based supramolecular hydrogels, but the precise way in which the molecular structure impacts the stability of the G4-structures and consequently affects the properties of the hydrogel remains to be elucidated. This issue represents a notable challenge in the field, which restricts their further applications to some extent. In this study, single crystals of Gd, αGd and αGd* were cultivated and compared with G. Notably, before this study, the single crystal structures of all natural nucleosides, with the exception of Gd, had been determined. The investigation into the molecular structure and supramolecular self-assembly properties of four guanosine analogs at the atomic scale revealed that the formation of G-quartets is critical for their ability to form hydrogels. The stability of the sugar ring geometry conformation (an intrinsic factor) and the disorder and strength of the hydration effect (extrinsic factors) are vital for maintaining the stability of the G4-structures. The rapid cooling changes the molecular geometry conformation, and the organic solvent changes the hydration effect, which can improve the longevity stability of G4-structure-based supramolecular hydrogels instead of chemical modifications. Consequently, the lifespan of the hydrogels was extended from 2 h to over one week. This advancement is expected to offer significant insights for future research in designing and developing G4-structure-based supramolecular hydrogels.

Rational design and synthesis of triazene-amonafide derivatives as novel potential antitumor agents causing oxidative damage towards DNA through intercalation mode

In this work, we rationally designed and synthesized two novel triazene-amonafide derivatives 2-(2-(diisopropylamino)ethyl)-5-(3,3-dimethyltriaz-1-en-1-yl)-1H-benzo[de]isoquinoline-1,3(2H)-dione (D-11) and 5-(3,3-diethyltriaz-1-en-1-yl)-2-(2-(diisopropylamino)ethyl)-1H-benzo[de]isoquinoline-1,3(2H)-dione (D-12) as potential antitumor agents. The DNA damage induced by the intercalation mode of D-11 (D-12) towards DNA was electrochemically detected through the construction of efficient biosensors. The consecutive processes of reversible redox of naphthylimide ring and irreversible oxidation of triazene moiety were elucidated on the surface of glassy carbon electrode (GCE) by CV, SWV, and DPV methods. Electrochemical biosensors were obtained through the immobilization of ctDNA, G-quadruplexes, poly(dG), and poly(dA), respectively, on the clean surface of GCE. After the incubation of biosensors with D-11 or D-12, the peaks of dGuo and dAdo decreased prominently, and the peak of 8-oxoGua appeared at +0.50 V, suggesting that the interaction between D-11 (D-12) and DNA could result in the oxidative damage of guanine. Unexpected, the as-prepared DNA biosensor possessed satisfactory anti-interference property and good practicability in real samples. UV-vis and fluorescence spectra, and gel electrophoresis assays were employed to further confirm the intercalation mode of D-11 (D-12) towards DNA base pairs. Moreover, D-11 was proved to exhibit stronger anti-proliferation activity than mitionafide and amonafide against both A549 and HeLa cell lines.

The pseudotorsional space of RNA

The characterization of the conformational landscape of the RNA backbone is rather complex due to the ability of RNA to assume a large variety of conformations. These backbone conformations can be depicted by pseudotorsional angles linking RNA backbone atoms, from which Ramachandran-like plots can be built. We explore here different definitions of these pseudotorsional angles, finding that the most accurate ones are the traditional η (eta) and θ (theta) angles, which represent the relative position of RNA backbone atoms P and C4'. We explore the distribution of η - θ in known experimental structures, comparing the pseudotorsional space generated with structures determined exclusively by one experimental technique. We found that the complete picture only appears when combining data from different sources. The maps provide a quite comprehensive representation of the RNA accessible space, which can be used in RNA-structural predictions. Finally, our results highlight that protein interactions lead to significant changes in the population of the η - θ space, pointing toward the role of induced-fit mechanisms in protein-RNA recognition.

Synthesis and properties of fully-modified 4'-selenoRNA, an endonuclease-resistant RNA analog

A large number of chemically modified oligonucleotides (ONs) have been developed for RNA-based technologies. In most modified RNAs, the characteristic 2'-hydroxyl (2'-OH) groups are removed to enhance both nuclease resistance and hybridization ability. However, the importance of the 2'-OH group in RNA structure and function is well known. Here, we report the synthesis and properties of 4'-selenoRNA in which all four nucleoside units retain the 2'-OH groups but contain a selenium atom instead of an oxygen atom at the 4'-position of the furanose ring. 4'-SelenoRNA has enhanced ability to form duplexes with RNA, and high endonuclease resistance despite the presence of the 2'-OH groups. X-ray crystallography analysis showed that the 4'-selenoRNA duplex adopts an A-conformation, similar to natural RNA, although one 4'-selenocytidine residue has unusual South-type sugar puckering. Furthermore, preliminary studies using 4'-seleno-modified siRNAs suggest that 4'-selenoRNA may be applicable to RNA interference technology. Collectively, our results raise the possibility of a new class of modified RNA in which 2'-OH groups do not need to be masked.

Modification at the C2'-O-Position with 2-Methylbenzothiophene Induces Unique Structural Changes and Thermal Transitions on Duplexes of RNA and DNA

Oligonucleotides can be chemically modified for a variety of applications that include their use as biomaterials, in therapeutics, or as tools to understand biochemical processes, among others. This work focuses on the functionalization of oligonucleotides of RNA and DNA (12- or 14-nucleotides long) with methylbenzothiophene (BT), at the C2'-O-position, which led to unique structural features. Circular dichroism (CD) analyses showed that positioning the BT units on one strand led to significant thermal destabilization, while duplexes where each strand contained 4-BT rings formed a distinct arrangement with cooperativity/interactions among the modifications (evidenced from the appearance of a band with positive ellipticity at 235 nm). Interestingly, the structural arrays displayed increased duplex stabilization (>10 °C higher than the canonical analogue) as a function of [Na⁺] with an unexpected structural rearrangement at temperatures above 50 °C. Density functional theory-polarizable continuum model (DFT-PCM) calculations were carried out, and the analyses were in agreement with induced structural changes as a function of salt content. A model was proposed where the hydrophobic surface allows for an internal nucleobase rearrangement into a more thermodynamically stable structure, before undergoing full denaturation, with increased heat. While this behavior is not common, B- to Z-form duplex transitions can occur and are dependent on parameters that were probed in this work, i.e., temperature, nature of modification, or ionic content. To take advantage of this phenomenon, we probed the ability of the modified duplexes to be recognized by Zα (an RNA binding protein that targets Z-form RNA) via electrophoretic analysis and CD. Interestingly, the protein did not bind to canonical duplexes of DNA or RNA; however, it recognized the modified duplexes, in a [monovalent/divalent salt] dependent manner. Overall, the findings describe methodology to attain unique structural motifs of modified duplexes of DNA or RNA, and control their behavior as a function of salt concentration. While their affinity to RNA binding proteins, and the corresponding mechanism of action, requires further exploration, the tunable properties can be of potential use to study this, and other, types of modifications. The novel arrays that formed, under the conditions described herein, provide a useful way to explore the structure and behavior of modified oligonucleotides, in general.

Cryo neutron crystallography demonstrates influence of RNA 2'-OH orientation on conformation, sugar pucker and water structure

The ribose 2′-hydroxyl is the key chemical difference between RNA and DNA and primary source of their divergent structural and functional characteristics. Macromolecular X-ray diffraction experiments typically do not reveal the positions of hydrogen atoms. Thus, standard crystallography cannot determine 2′-OH orientation (H2′-C2′-O2′-HO2′ torsion angle) and its potential roles in sculpting the RNA backbone and the expansive fold space. Here, we report the first neutron crystal structure of an RNA, the Escherichia coli rRNA Sarcin-Ricin Loop (SRL). 2′-OD orientations were established for all 27 residues and revealed O-D bonds pointing toward backbone (O3′, 13 observations), nucleobase (11) or sugar (3). Most riboses in the SRL stem region show a 2′-OD backbone-orientation. GAGA-tetraloop riboses display a 2′-OD base-orientation. An atypical C2′-endo sugar pucker is strictly correlated with a 2′-OD sugar-orientation. Neutrons reveal the strong preference of the 2′-OH to donate in H-bonds and that 2′-OH orientation affects both backbone geometry and ribose pucker. We discuss 2′-OH and water molecule orientations in the SRL neutron structure and compare with results from a solution phase 10 μs MD simulation. We demonstrate that joint cryo-neutron/X-ray crystallography offers an all-in-one approach to determine the complete structural properties of RNA, i.e. geometry, conformation, protonation state and hydration structure.

Beyond ribose and phosphate: Selected nucleic acid modifications for structure-function investigations and therapeutic applications

Over the past 25 years, the acceleration of achievements in the development of oligonucleotide-based therapeutics has resulted in numerous new drugs making it to the market for the treatment of various diseases. Oligonucleotides with alterations to their scaffold, prepared with modified nucleosides and solid-phase synthesis, have yielded molecules with interesting biophysical properties that bind to their targets and are tolerated by the cellular machinery to elicit a therapeutic outcome. Structural techniques, such as crystallography, have provided insights to rationalize numerous properties including binding affinity, nuclease stability, and trends observed in the gene silencing. In this review, we discuss the chemistry, biophysical, and structural properties of a number of chemically modified oligonucleotides that have been explored for gene silencing.

Intrinsic Contributions of 2'-Hydroxyl to the Hydration of Nucleosides at the Monomeric Level

Although many reports have revealed structural differences between DNA and RNA at the polymeric level, there are no comparative studies with 2'-deoxyribonucleoside and ribonucleoside to explore the role of the 2'-OH group at the monomeric level under the same conditions. Inspired by this, herein, the intrinsic contributions of the 2'-OH group in the nucleoside have been systematically investigated by directly solving the single-crystal structures of 2'-deoxy-2-aminoadenosine (1), 2-aminoadenosine (2), and 2-aminoarabinofuranosyladenine (3) in water. The 2'-OH group not only influenced the conformation and base-pair pattern of the single-nucleoside molecule, but also played a fundamental role in the entire supramolecular structure. Interestingly, compound 1, which did not contain the 2'-OH group, displayed strong hydration, whereas 2 and 3 (with the 2'-OH group in the opposite direction) exhibited no hydration, which was completely different from that observed in nucleic acids. Meanwhile, compound 1 trapped water molecules to form unique trihydrol moieties, which further served as the backbone to construct the simplest double-chain DNA-like structures. To this end, to investigate the effect of the biological environment on these unique structures, the solvent was changed from water to phosphate-buffered saline (PBS). Surprisingly, such a subtle adjustment led to entirely different superstructures, consisting of 2D lamellar structures in water and 3D porous structures in PBS. These large morphological differences could be attributed to delicate ion hydration, which was also confirmed through variable-temperature X-ray analysis, SEM, and intermolecular interaction energy calculations. In summary, this study comprehensively investigated the intrinsic contributions of 2'-hydroxyl to the hydration of nucleosides at the monomeric level; this is helpful to further understand the differences in DNA/RNA and the impact of their surrounding environment.

DFFR: A New Method for High-Throughput Recalibration of Automatic Force-Fields for Drugs

We present drug force-field recalibration (DFFR), a new method for refining of automatic force-fields used to represent small drugs in docking and molecular dynamics simulations. The method is based on fine-tuning of torsional terms to obtain ensembles that reproduce observables derived from reference data. DFFR is fast and flexible and can be easily automatized for a high-throughput regime, making it useful in drug-design projects. We tested the performance of the method in a few model systems and also in a variety of druglike molecules using reference data derived from: (i) density functional theory coupled to a self-consistent reaction field (DFT/SCRF) calculations on highly populated conformers and (ii) enhanced sampling quantum mechanical/molecular mechanics (QM/MM) where the drug is reproduced at the QM level, while the solvent is represented by classical force-fields. Extension of the method to include other sources of reference data is discussed.

Coming Together: RNAs and Proteins Assemble under the Single-Molecule Fluorescence Microscope

RNAs, across their numerous classes, often work in concert with proteins in RNA-protein complexes (RNPs) to execute critical cellular functions. Ensemble-averaging methods have been instrumental in revealing many important aspects of these RNA-protein interactions, yet are insufficiently sensitive to much of the dynamics at the heart of RNP function. Single-molecule fluorescence microscopy (SMFM) offers complementary, versatile tools to probe RNP conformational and compositional changes in detail. In this review, we first outline the basic principles of SMFM as applied to RNPs, describing key considerations for labeling, imaging, and quantitative analysis. We then sample applications of in vitro and in vivo single-molecule visualization using the case studies of pre-messenger RNA (mRNA) splicing and RNA silencing, respectively. After discussing specific insights single-molecule fluorescence methods have yielded, we briefly review recent developments in the field and highlight areas of anticipated growth.

Polarizable force field for RNA based on the classical drude oscillator

RNA molecules are highly dynamic and capable of adopting a wide range of complex, folded structures. The factors driving the folding and dynamics of these structures are dependent on a balance of base pairing, hydration, base stacking, ion interactions, and the conformational sampling of the 2'-hydroxyl group in the ribose sugar. The representation of these features is a challenge for empirical force fields used in molecular dynamics simulations. Toward meeting this challenge, the inclusion of explicit electronic polarization is important in accurately modeling RNA structure. In this work, we present a polarizable force field for RNA based on the classical Drude oscillator model, which represents electronic degrees of freedom via negatively charged particles attached to their parent atoms by harmonic springs. Beginning with parametrization against quantum mechanical base stacking interaction energy and conformational energy data, we have extended the Drude-2017 nucleic acid force field to include RNA. The conformational sampling of a range of RNA sequences were used to validate the force field, including canonical A-form RNA duplexes, stem-loops, and complex tertiary folds that bind multiple Mg2+ ions. Overall, the Drude-2017 RNA force field reproduces important properties of these structures, including the conformational sampling of the 2'-hydroxyl and key interactions with Mg2+ ions. © 2018 Wiley Periodicals, Inc.

Effect of Ribose Conformation on RNA Cleavage via Internal Transesterification

RNA cleavage via internal transesterification is a fundamental reaction involved in RNA processing and metabolism, and the regulation thereof. Herein, the influence of ribose conformation on this reaction was investigated with conformationally constrained ribonucleotides. RNA cleavage rates were found to decrease in the order South-constrained ribonucleotide > native ribonucleotide ≫ North-constrained counterpart, indicating that the ribose conformation plays an important role in modulating RNA cleavage via internal transesterification.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Modeling of canonical and C2′-O-thiophenylmethyl modified hexamers of RNA. Insights into the nature of structural changes and thermal stability

Modification of the C2′-O-position with thiophenylmethyl groups on both strands leads to thermal stabilization of the duplex. Predicting the effects that modifications will have on structure of RNA is of importance in the development of new RNA technologies.

H. R, Thamleena A. H, Joseph A, Ben A, V. U. K. Roles of different amino-acid residues towards binding and selective transport of K+ through KcsA K+-ion channel

Each amino acid in the selectivity filter plays a distinct role towards binding and transport of K⁺ ion through KcsA.

Translational roles of the C75 2'OH in an in vitro tRNA transcript at the ribosomal A, P and E sites

Aminoacyl-tRNAs containing a deoxy substitution in the penultimate nucleotide (C75 2′OH → 2′H) have been widely used in translation for incorporation of unnatural amino acids (AAs). However, this supposedly innocuous modification surprisingly increased peptidyl-tRNA^Ala_ugc drop off in biochemical assays of successive incorporations. Here we predict the function of this tRNA 2′OH in the ribosomal A, P and E sites using recent co-crystal structures of ribosomes and tRNA substrates and test these structure-function models by systematic kinetics analyses. Unexpectedly, the C75 2′H did not affect A- to P-site translocation nor peptidyl donor activity of tRNA^Ala_ugc. Rather, the peptidyl acceptor activity of the A-site Ala-tRNA^Ala_ugc and the translocation of the P-site deacylated tRNA^Ala_ugc to the E site were impeded. Delivery by EF-Tu was not significantly affected. This broadens our view of the roles of 2′OH groups in tRNAs in translation.