Prediction of 3D RNA Structures from Sequence Using Energy Landscapes of RNA Dimers: Application to RNA Tetraloops.

Access to the three-dimensional structure of RNA enables an ability to gain a more profound understanding of its biological mechanisms, as well as the ability to design RNA-targeting drugs, which can take advantage of the unique chemical environment imposed by a folded RNA structure. Due to the dynamic and structurally complex properties of RNA, both experimental and traditional computational methods have difficulty in determining RNA's 3D structure. Herein, we introduce TAPERSS (Theoretical Analyses, Prediction, and Evaluation of RNA Structures from Sequence), a physics-based fragment assembly method for predicting 3D RNA structures from sequence. Using a fragment library created using discrete path sampling calculations of RNA dinucleoside monophosphates, TAPERSS can sample the physics-based energy landscapes of any RNA sequence with relatively low computational complexity. We have benchmarked TAPERSS on 21 RNA tetraloops, using a combinatorial algorithm as a proof-of-concept. We show that TAPERSS was successfully able to predict the apo-state structures of all 21 RNA hairpins, with 16 of those structures also having low predicted energies as well. We demonstrate that TAPERSS performs most accurately on GNRA-like tetraloops with mostly stacked loop-nucleotides, while having limited success with more dynamic UNCG and CUYG tetraloops, most likely due to the influence of the RNA force field used to create the fragment library. Moreover, we show that TAPERSS can successfully predict the majority of the experimental non-apo states, highlighting its potential in anticipating biologically significant yet unobserved states. This holds great promise for future applications in drug design and related studies. With discussed improvements and implementation of more efficient sampling algorithms, we believe TAPERSS may serve as a useful tool for a physics-based conformational sampling of large RNA structures.

Conformational dynamics of RNA G4C2 and G2C4 repeat expansions causing ALS/FTD using NMR and molecular dynamics studies.

G4C2 and G2C4 repeat expansions in chromosome 9 open reading frame 72 (C9orf72) are the most common cause of genetically defined amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), or c9ALS/FTD. The gene is bidirectionally transcribed, producing G4C2 repeats [r(G4C2)exp] and G2C4 repeats [r(G2C4)exp]. The c9ALS/FTD repeat expansions are highly structured, and structural studies showed that r(G4C2)exp predominantly folds into a hairpin with a periodic array of 1 × 1 G/G internal loops and a G-quadruplex. A small molecule probe revealed that r(G4C2)exp also adopts a hairpin structure with 2 × 2 GG/GG internal loops. We studied the conformational dynamics adopted by 2 × 2 GG/GG loops using temperature replica exchange molecular dynamics (T-REMD) and further characterized the structure and underlying dynamics using traditional 2D NMR techniques. These studies showed that the loop's closing base pairs influence both structure and dynamics, particularly the configuration adopted around the glycosidic bond. Interestingly, r(G2C4) repeats, which fold into an array of 2 × 2 CC/CC internal loops, are not as dynamic. Collectively, these studies emphasize the unique sensitivity of r(G4C2)exp to small changes in stacking interactions, which is not observed in r(G2C4)exp, providing important considerations for further principles in structure-based drug design.

Dynamic docking of small molecules targeting RNA CUG repeats causing myotonic dystrophy type 1.

Expansion of RNA CUG repeats causes myotonic dystrophy type 1 (DM1). Once transcribed, the expanded CUG repeats strongly attract muscleblind-like 1 (MBNL1) proteins and disturb their functions in cells. Because of its unique structural form, expanded RNA CUG repeats are prospective drug targets, where small molecules can be utilized to target RNA CUG repeats to inhibit MBNL1 binding and ameliorate DM1-associated defects. In this contribution, we developed two physics-based dynamic docking approaches (DynaD and DynaD/Auto) and applied them to nine small molecules known to specifically target RNA CUG repeats. While DynaD uses a distance-based reaction coordinate to study the binding phenomenon, DynaD/Auto combines results of umbrella sampling calculations performed on 1 × 1 UU internal loops and AutoDock calculations to efficiently sample the energy landscape of binding. Predictions are compared with experimental data, displaying a positive correlation with correlation coefficient (R) values of 0.70 and 0.81 for DynaD and DynaD/Auto, respectively. Furthermore, we found that the best correlation was achieved with MM/3D-RISM calculations, highlighting the importance of solvation in binding calculations. Moreover, we detected that DynaD/Auto performed better than DynaD because of the use of prior knowledge about the binding site arising from umbrella sampling calculations. Finally, we developed dendrograms to present how bound states are connected to each other in a binding process. Results are exciting, as DynaD and DynaD/Auto will allow researchers to utilize two novel physics-based and computer-aided drug-design methodologies to perform in silico calculations on drug-like molecules aiming to target complex RNA loops.

Evaluating Geometric Definitions of Stacking for RNA Dinucleoside Monophosphates Using Molecular Mechanics Calculations.

RNA modulation via small molecules is a novel approach in pharmacotherapies, where the determination of the structural properties of RNA motifs is considered a promising way to develop drugs capable of targeting RNA structures to control diseases. However, due to the complexity and dynamic nature of RNA molecules, the determination of RNA structures using experimental approaches is not always feasible, and computational models employing force fields can provide important insight. The quality of the force field will determine how well the predictions are compared to experimental observables. Stacking in nucleic acids is one such structural property, originating mainly from London dispersion forces, which are quantum mechanical and are included in molecular mechanics force fields through nonbonded interactions. Geometric descriptions are utilized to decide if two residues are stacked and hence to calculate the stacking free energies for RNA dinucleoside monophosphates (DNMPs) through statistical mechanics for comparison with experimental thermodynamics data. Here, we benchmark four different stacking definitions using molecular dynamics (MD) trajectories for 16 RNA DNMPs produced by two different force fields (RNA-IL and ff99OL3) and show that our stacking definition better correlates with the experimental thermodynamics data. While predictions within an accuracy of 0.2 kcal/mol at 300 K were observed in RNA CC, CU, UC, AG, GA, and GG, stacked states of purine-pyrimidine and pyrimidine-purine DNMPs, respectively, were typically underpredicted and overpredicted. Additionally, population distributions of RNA UU DNMPs were poorly predicted by both force fields, implying a requirement for further force field revisions. We further discuss the differences predicted by each RNA force field. Finally, we show that discrete path sampling (DPS) calculations can provide valuable information and complement the MD simulations. We propose the use of experimental thermodynamics data for RNA DNMPs as benchmarks for testing RNA force fields.

Computational Investigation of Bending Properties of RNA AUUCU, CCUG, CAG, and CUG Repeat Expansions Associated With Neuromuscular Disorders.

Expansions of RNA AUUCU, CCUG, CAG, and CUG repeats cause spinocerebellar ataxia type 10, myotonic dystrophy type 2, Huntington's disease, and myotonic dystrophy type 1, respectively. By performing extensive molecular dynamic simulations, we investigated the bending propensities and conformational landscapes adopted by 3×3, 2×2, and 1×1 internal loops observed in RNA AUUCU, CCUG, CAG, and CUG repeat expansions using model systems having biologically relevant repeat sizes. We show that the conformational variability experienced by these loops is more complex than previous reports where a variety of unconventional hydrogen bonds are formed. At the global scale, strong bending propensity was observed in r(AUUCU)10, r(CCUG)15, r(CAG)20, and r(CUG)20, and, to a lesser extent, in r(AUUCU)4, r(CCUG)10, r(CAG)10, and r(CUG)10. Furthermore, RNA CAG repeats exhibit a tendency toward bent states with more than 50% of observed conformations having bending angles greater than 50°, while RNA CUG repeats display relatively linear-like conformations with extremely bent conformations accounting for less than 25% of the observed structures. Conformations experienced by RNA AUUCU repeats are a combination of strongly bent and kinked structures. The bent states in RNA CCUG repeats mostly fall into the moderately bent category with a marginal ensemble experiencing extreme bending. The general pattern observed in all the bent structures indicates the collapse of the major groove width as the mechanical trigger for bending, which is caused by alteration of base pair step parameters at multiple locations along the RNA due to local distortions at the loop sites. Overextension is also observed in all the RNA repeats that is attributed to widening of the major groove width as well as undertwisting phenomenon. This information and the rich structural repository could be applied for structure based small molecule design targeting disease-causing RNAs. The bending propensities of these constructs, at the global level, could also have implications on how expanded RNA repeats interact with proteins.

A Small Molecule Exploits Hidden Structural Features within the RNA Repeat Expansion That Causes c9ALS/FTD and Rescues Pathological Hallmarks.

The hexanucleotide repeat expansion GGGGCC [r(G4C2)exp] within intron 1 of C9orf72 causes genetically defined amyotrophic lateral sclerosis and frontotemporal dementia, collectively named c9ALS/FTD. , the repeat expansion causes neurodegeneration via deleterious phenotypes stemming from r(G4C2)exp RNA gain- and loss-of-function mechanisms. The r(G4C2)exp RNA folds into both a hairpin structure with repeating 1 × 1 nucleotide GG internal loops and a G-quadruplex structure. Here, we report the identification of a small molecule (CB253) that selectively binds the hairpin form of r(G4C2)exp. Interestingly, the small molecule binds to a previously unobserved conformation in which the RNA forms 2 × 2 nucleotide GG internal loops, as revealed by a series of binding and structural studies. NMR and molecular dynamics simulations suggest that the r(G4C2)exp hairpin interconverts between 1 × 1 and 2 × 2 internal loops through the process of strand slippage. We provide experimental evidence that CB253 binding indeed shifts the equilibrium toward the 2 × 2 GG internal loop conformation, inhibiting mechanisms that drive c9ALS/FTD pathobiology, such as repeat-associated non-ATG translation formation of stress granules and defective nucleocytoplasmic transport in various cellular models of c9ALS/FTD.

Ribonuclease recruitment using a small molecule reduced c9ALS/FTD r(G4C2) repeat expansion in vitro and in vivo ALS models.

The most common cause of amyotrophic lateral sclerosis and frontotemporal dementia (c9ALS/FTD) is an expanded G4C2 RNA repeat [r(G4C2)exp] in chromosome 9 open reading frame 72 (C9orf72), which elicits pathology through several mechanisms. Here, we developed and characterized a small molecule for targeted degradation of r(G4C2)exp. The compound was able to selectively bind r(G4C2)exp's structure and to assemble an endogenous nuclease onto the target, provoking removal of the transcript by native RNA quality control mechanisms. In c9ALS patient­derived spinal neurons, the compound selectively degraded the mutant C9orf72 allele with limited off-targets and reduced quantities of toxic dipeptide repeat proteins (DPRs) translated from r(G4C2)exp. In vivo work in a rodent model showed that abundance of both the mutant allele harboring the repeat expansion and DPRs were selectively reduced by this compound. These results demonstrate that targeted small-molecule degradation of r(G4C2)exp is a strategy for mitigating c9ALS/FTD-associated pathologies and studying disease-associated pathways in preclinical models.

Massively Parallel Optimization of the Linker Domain in Small Molecule Dimers Targeting a Toxic r(CUG) Repeat Expansion.

RNA contributes to disease pathobiology and is an important therapeutic target. The downstream biology of disease-causing RNAs can be short-circuited with small molecules that recognize structured regions. The discovery and optimization of small molecules interacting with RNA is, however, challenging. Herein, we demonstrate a massively parallel one-bead-one-compound methodology, employed to optimize the linker region of a dimeric compound that binds the toxic r(CUG) repeat expansion [r(CUG)exp] causative of myotonic dystrophy type 1 (DM1). Indeed, affinity selection on a 331,776-member library allowed the discovery of a compound with enhanced potency both in vitro (10-fold) and in DM1-patient-derived myotubes (5-fold). Molecular dynamics simulations revealed additional interactions between the optimized linker and the RNA, resulting in ca. 10 kcal/mol lower binding free energy. The compound was conjugated to a cleavage module, which directly cleaved the transcript harboring the r(CUG)exp and alleviated disease-associated defects.

A Small Molecule that Binds an RNA Repeat Expansion Stimulates Its Decay via the Exosome Complex.

Many diseases are caused by toxic RNA repeats. Herein, we designed a lead small molecule that binds the structure of the r(CUG) repeat expansion [r(CUG)exp] that causes myotonic dystrophy type 1 (DM1) and Fuchs endothelial corneal dystrophy (FECD) and rescues disease biology in patient-derived cells and in vivo. Interestingly, the compound's downstream effects are different in the two diseases, owing to the location of the repeat expansion. In DM1, r(CUG)exp is harbored in the 3' untranslated region, and the compound has no effect on the mRNA's abundance. In FECD, however, r(CUG)exp is located in an intron, and the small molecule facilitates excision of the intron, which is then degraded by the RNA exosome complex. Thus, structure-specific, RNA-targeting small molecules can act disease specifically to affect biology, either by disabling the gain-of-function mechanism (DM1) or by stimulating quality control pathways to rid a disease-affected cell of a toxic RNA (FECD).

A general fragment-based approach to identify and optimize bioactive ligands targeting RNA.

RNAs have important functions that are dictated by their structure. Indeed, small molecules that interact with RNA structures can perturb function, serving as chemical probes and lead medicines. Here we describe the development of a fragment-based approach to discover and optimize bioactive small molecules targeting RNA. We extended the target validation method chemical cross-linking and isolation by pull-down (Chem-CLIP) to identify and map the binding sites of low molecular weight fragments that engage RNA or Chem-CLIP fragment mapping (Chem-CLIP-Frag-Map). Using Chem-CLIP-Frag-Map, we identified several fragments that bind the precursor to oncogenic microRNA-21 (pre-miR-21). Assembly of these fragments provided a specific bioactive compound with improved potency that inhibits pre-miR-21 processing, reducing mature miR-21 levels. The compound exerted selective effects on the transcriptome and selectively mitigated a miR-21-associated invasive phenotype in triple-negative breast cancer cells. The Chem-CLIP-Frag-Map approach should prove general to expedite the identification and optimization of small molecules that bind RNA targets.

Structural Features of Small Molecules Targeting the RNA Repeat Expansion That Causes Genetically Defined ALS/FTD.

Genetically defined amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), collectively named c9ALS/FTD, are triggered by hexanucleotide GGGGCC repeat expansions [r(G4C2)exp] within the C9orf72 gene. In these diseases, neuronal loss occurs through an interplay of deleterious phenotypes, including r(G4C2)exp RNA gain-of-function mechanisms. Herein, we identified a benzimidazole derivative, CB096, that specifically binds to a repeating 1 × 1 GG internal loop structure, 5'CGG/3'GGC, that is formed when r(G4C2)exp folds. Structure-activity relationship (SAR) studies and molecular dynamics (MD) simulations were used to define the molecular interactions formed between CB096 and r(G4C2)exp that results in the rescue of disease-associated pathways. Overall, this study reveals a unique structural feature within r(G4C2)exp that can be exploited for the development of lead medicines and chemical probes.

Design of a small molecule that stimulates vascular endothelial growth factor A enabled by screening RNA fold-small molecule interactions.

Vascular endothelial growth factor A (VEGFA) stimulates angiogenesis in human endothelial cells, and increasing its expression is a potential treatment for heart failure. Here, we report the design of a small molecule (TGP-377) that specifically and potently enhances VEGFA expression by the targeting of a non-coding microRNA that regulates its expression. A selection-based screen, named two-dimensional combinatorial screening, revealed preferences in small-molecule chemotypes that bind RNA and preferences in the RNA motifs that bind small molecules. The screening program increased the dataset of known RNA motif-small molecule binding partners by 20-fold. Analysis of this dataset against the RNA-mediated pathways that regulate VEGFA defined that the microRNA-377 precursor, which represses Vegfa messenger RNA translation, is druggable in a selective manner. We designed TGP-377 to potently and specifically upregulate VEGFA in human umbilical vein endothelial cells. These studies illustrate the power of two-dimensional combinatorial screening to define molecular recognition events between 'undruggable' biomolecules and small molecules, and the ability of sequence-based design to deliver efficacious structure-specific compounds.

Macrocyclization of a Ligand Targeting a Toxic RNA Dramatically Improves Potency.

RNA molecules both contribute to and are causative of many human diseases. One method to perturb RNA function is to target its structure with small molecules. However, discovering bioactive ligands for RNA targets is challenging. Here, we show that the bioactivity of a linear dimeric ligand that inactivates the RNA trinucleotide repeat expansion that causes myotonic dystrophy type 1 [DM1; r(CUG)exp ] can be improved by macrocyclization. Indeed, the macrocyclic compound is ten times more potent than the linear compound for improving DM1-associated defects in cells, including in patient-derived myotubes (muscle cells). This enhancement in potency is due to the macrocycle's increased affinity and selectively for the target, which inhibit r(CUG)exp 's toxic interaction with muscleblind-like 1 (MBNL1), and its superior cell permeability. Macrocyclization could prove to be an effective way to enhance the bioactivity of modularly assembled ligands targeting RNA.

Small-molecule targeted recruitment of a nuclease to cleave an oncogenic RNA in a mouse model of metastatic cancer.

As the area of small molecules interacting with RNA advances, general routes to provide bioactive compounds are needed as ligands can bind RNA avidly to sites that will not affect function. Small-molecule targeted RNA degradation will thus provide a general route to affect RNA biology. A non-oligonucleotide-containing compound was designed from sequence to target the precursor to oncogenic microRNA-21 (pre-miR-21) for enzymatic destruction with selectivity that can exceed that for protein-targeted medicines. The compound specifically binds the target and contains a heterocycle that recruits and activates a ribonuclease to pre-miR-21 to substoichiometrically effect its cleavage and subsequently impede metastasis of breast cancer to lung in a mouse model. Transcriptomic and proteomic analyses demonstrate that the compound is potent and selective, specifically modulating oncogenic pathways. Thus, small molecules can be designed from sequence to have all of the functional repertoire of oligonucleotides, including inducing enzymatic degradation, and to selectively and potently modulate RNA function in vivo.

Improvement of RNA Simulations with Torsional Revisions of the AMBER Force Field.

Our current knowledge on the unique roles of RNA in cells makes it vital to investigate the properties of RNA systems using computational methods because of the potential pharmaceutical applications. With the continuous advancement of computer technology, it is now possible to study RNA folding. Molecular mechanics calculations are useful in discovering the structural and thermodynamic properties of RNA systems. Yet, the predictions depend on the quality of the RNA force field, which is a set of parameters describing the potential energy of the system. Torsional parameters are one of the terms in a force field that can be revised using physics-based approaches. This chapter focuses on improvements provided by revisions of torsional parameters of the AMBER (Assisted Model Building with Energy Refinement) RNA force field. The theory behind torsional revisions and re-parameterization of several RNA torsions is briefly described. Applications of the revised torsional parameters to study RNA nucleosides, single-stranded RNA tetramers, and RNA repeat expansions are described in detail. It is concluded that RNA force fields require constant revisions and should be benchmarked against diverse RNA systems such as single strands and internal loops in order to test their qualities.

The Hairpin Form of r(G4C2)exp in c9ALS/FTD Is Repeat-Associated Non-ATG Translated and a Target for Bioactive Small Molecules.

The most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) is an expanded G4C2 repeat [(G4C2)exp] in C9ORF72. ALS/FTD-associated toxicity has been traced to the RNA transcribed from the repeat expansion [r(G4C2)exp], which sequesters RNA-binding proteins (RBPs) and undergoes repeat-associated non-ATG (RAN) translation to generate toxic dipeptide repeats. Using in vitro and cell-based assays, we identified a small molecule (4) that selectively bound r(G4C2)exp, prevented sequestration of an RBP, and inhibited RAN translation. Indeed, biophysical characterization showed that 4 selectively bound the hairpin form of r(G4C2)exp, and nuclear magnetic resonance spectroscopy studies and molecular dynamics simulations defined this molecular recognition event. Cellular imaging revealed that 4 localized to r(G4C2)exp cytoplasmic foci, the putative sites of RAN translation. Collectively, these studies highlight that the hairpin structure of r(G4C2)exp is a therapeutically relevant target and small molecules that bind it can ameliorate c9ALS/FTD-associated toxicity.

Computational Investigation of RNA A-Bulges Related to the Microtubule-Associated Protein Tau Causing Frontotemporal Dementia and Parkinsonism.

Mutations in the human tau gene result in alternative splicing of the tau protein, which causes frontotemporal dementia and Parkinsonism. One disease mechanism is linked to the stability of a hairpin within the microtubule-associated protein tau (MAPT) mRNA, which contains an A-bulge. Here we employ computational methods to investigate the structural and thermodynamic properties of several A-bulge RNAs with different closing base-pairs. We find that the current amber RNA force field has a preference to overstabilize base-triple over stacked states, even though some of the A-bulges are known to prefer stacked states according to NMR studies. We further determined that if the neighboring base-pairs of A-bulges are AU, this situation can lead to base slippage. However, when the 3'-side of the A-bulge has an UA base-pair, the stacked state is stabilized by an extra interaction that is not observed in the other sequences. We suggest that these A-bulge RNA systems could be used as benchmarks to improve the current RNA force fields.

Structure and Dynamics of Electron Injection and Charge Recombination in i-Motif DNA Conjugates.

The dynamics of electron injection have been investigated in intramolecular i-motif conjugates possessing stilbenediether (Sd) and perylenediimide (PDI) chromophores separated by either four or six hemiprotonated cytosine C-C base pairs assembled with synthetic loops. Circular dichroism spectra are consistent with the formation of i-motif structures in the absence or presence of Sd and PDI chromophores. The fluorescence of the Sd chromophore is essentially completely quenched by neighboring C-C base pairs, consistent with their function as an electron donor and electron acceptor, respectively. However, the fluorescence of the PDI chromophore is only partially quenched. The dynamics of electron injection from singlet Sd into the i-motif and subsequent charge recombination has been determined by femtosecond transient absorption (fsTA) spectroscopy and compared with the results for electron injection and charge recombination in Sd-linked hairpins possessing cytosine-guanine (C-G) or 5-fluorouracil-adenine (F-A) base pairs. While charge injection is ultrafast (<0.8 ps) for the i-motifs, charge transport across the i-motif C-C base pairs to the PDI electron trap is not observed. The absence of electron transport is related to the structure of the stacked C-C base pairs in the i-motif.

Structure and Dynamics of RNA Repeat Expansions That Cause Huntington's Disease and Myotonic Dystrophy Type 1.

RNA repeat expansions cause a host of incurable, genetically defined diseases. The most common class of RNA repeats consists of trinucleotide repeats. These long, repeating transcripts fold into hairpins containing 1 × 1 internal loops that can mediate disease via a variety of mechanism(s) in which RNA is the central player. Two of these disorders are Huntington's disease and myotonic dystrophy type 1, which are caused by r(CAG) and r(CUG) repeats, respectively. We report the structures of two RNA constructs containing three copies of a r(CAG) [r(3×CAG)] or r(CUG) [r(3×CUG)] motif that were modeled with nuclear magnetic resonance spectroscopy and simulated annealing with restrained molecular dynamics. The 1 × 1 internal loops of r(3×CAG) are stabilized by one-hydrogen bond (cis Watson-Crick/Watson-Crick) AA pairs, while those of r(3×CUG) prefer one- or two-hydrogen bond (cis Watson-Crick/Watson-Crick) UU pairs. Assigned chemical shifts for the residues depended on the identity of neighbors or next nearest neighbors. Additional insights into the dynamics of these RNA constructs were gained by molecular dynamics simulations and a discrete path sampling method. Results indicate that the global structures of the RNA are A-form and that the loop regions are dynamic. The results will be useful for understanding the dynamic trajectory of these RNA repeats but also may aid in the development of therapeutics.

Improving Computational Predictions of Single-Stranded RNA Tetramers with Revised α/γ Torsional Parameters for the Amber Force Field.

With current advancements in RNA based therapeutics, it is becoming crucial to utilize theoretical and computational methods to describe properly the physical properties of RNA molecules. NMR and X-ray crystallography are two powerful techniques for investigating structural properties. However, if the RNA molecules are complex or dynamic, these methods might not be adequate. For computational approaches, the quality of the force field will determine accuracy of our predictions. In this contribution, we revise the α/γ torsional parameters of RNA for amber force field using a model system representing an RNA dimer backbone. Combined with revised χ torsional parameters, previously shown to improve computational predictions, we benchmarked the revised force field on five single-stranded RNA (ssRNA) tetramers, three RNA dodecamer duplexes, and an RNA hairpin. A total of 60 µs of molecular dynamics (MD) simulations were run. We also employ the discrete path sampling (DPS) approach to compare the predictions for the revised amber force field with those for amber10. Our results indicate that the unphysical states observed with amber10 in ssRNA MD simulations are suppressed for the revised amber force field. In line with NMR experimental observations, incorporation of the revised α/γ and χ torsional parameters leads to A-form-like conformational states as the most favorable ssRNA tetramer conformations. Furthermore, the revised force field maintains the A-form geometry in regular RNA duplexes. Our revised amber force field for RNA should therefore improve structural and thermodynamic predictions for challenging RNA systems.