Elucidating the influence of RNA modifications and Magnesium ions on tRNA Phe conformational dynamics in S. cerevisiae : Insights from Replica Exchange Molecular Dynamics simulations.

Post-transcriptional modifications in RNA can significantly impact their structure and function. In particular, transfer RNAs (tRNAs) are heavily modified, with around 100 different naturally occurring nucleotide modifications contributing to codon bias and decoding efficiency. Here, we describe our efforts to investigate the impact of RNA modifications on the structure and stability of tRNA Phenylalanine (tRNA Phe ) from S. cerevisiae using molecular dynamics (MD) simulations. Through temperature replica exchange MD (T-REMD) studies, we explored the unfolding pathway to understand how RNA modifications influence the conformational dynamics of tRNA Phe , both in the presence and absence of magnesium ions (Mg 2+ ). We observe that modified nucleotides in key regions of the tRNA establish a complex network of hydrogen bonds and stacking interactions which is essential for tertiary structure stability of the tRNA. Furthermore, our simulations show that modifications facilitate the formation of ion binding sites on the tRNA. However, high concentrations of Mg 2+ ions can stabilize the tRNA tertiary structure in the absence of modifications. Our findings illuminate the intricate interactions between modifications, magnesium ions, and RNA structural stability.

Mutation of two intronic nucleotides alters RNA structure and dynamics inhibiting MBNL1 and RBFOX1 regulated splicing of the Insulin Receptor.

Alternative splicing (AS) of Exon 11 of the Insulin Receptor ( INSR ) is highly regulated and disrupted in several human disorders. To better understand INSR exon 11 AS regulation, splicing activity of an INSR exon 11 minigene reporter was measured across a gradient of the AS regulator muscleblind-like 1 protein (MBNL1). The RNA-binding protein Fox-1 (RBFOX1) was added to determine its impact on MBNL1-regulated splicing. The role of the RBFOX1 UGCAUG binding site within intron 11 was assessed across the MBNL1 gradient. Mutating the UGCAUG motif inhibited RBFOX1 regulation of exon 11 and had the unexpected effect of reducing MBNL1 regulation of this exon. Molecular dynamics simulations showed that exon 11 and the adjacent RNA adopts a dynamically stable conformation. Mutation of the RBFOX1 binding site altered RNA structure and dynamics, while a mutation that created an optimal MBNL1 binding site at the RBFOX1 site shifted the RNA back to wild type. An antisense oligonucleotide (ASO) was used to confirm the structure in this region of the pre-mRNA. This example of intronic mutations shifting pre-mRNA structure and dynamics to modulate splicing suggests RNA structure and dynamics should be taken into consideration for AS regulation and therapeutic interventions targeting pre-mRNA.

Aggregation in Aqueous Solutions of 2-(Tetrafluoro(trifluoromethyl)-λ6-sulfanyl-ethan-1-ol (CF3SF4-ethanol)): A Comparison with Aqueous Trifluoroethanol and Hexafluoroisopropanol Using Molecular Dynamics Simulations and Dynamic Light Scattering Experiments.

2-Tetrafluoro(trifluoromethyl)-λ6-sulfanylethan-1-ol (CF3SF4-ethanol) combines the polar hydrophobicity of tetrafluoro(trifluoromethyl)-λ6-sulfanyl (CF3SF4) group with the polarity of simple alcohols. The properties of aqueous solutions of the well-known fluorinated alcohols 2,2,2-trifluoroethanol (TFE) and 1,1,1,3,3,3-hexafluoroisopropanol (HFIP) were compared with those of aqueous solutions of the novel CF3SF4-ethanol. Those properties were computed using all atom molecular dynamics simulations with OPLS-compatible parameters. DFT ab initio calculations were used to accurately describe the nonsymmetrical, hypervalent sulfur in CF3SF4-ethanol. Although the molecular and conformational characteristics of CF3SF4-ethanol are like those of both TFE and HFIP, the greater hydrophobicity and lower polarity of CF3SF4-ethanol resulted in solution phase aggregation at a much lower concentration. The properties computed for TFE and HFIP in this work were consistent with published computational and experimental studies. CF3SF4-ethanol is predicted to be environmentally benign and hence an excellent green solvent candidate while possessing many of the same properties as TFE or HFIP.

Characterizing the Folding Transition-State Ensembles in the Energy Landscape of an RNA Tetraloop.

Molecular dynamics (MD) simulations have become increasingly powerful and can now describe the folding/unfolding of small biomolecules in atomic detail. However, a major challenge in MD simulations is to represent the complex energy landscape of biomolecules using a small number of reaction coordinates. In this study, we investigate the folding pathways of an RNA tetraloop, gcGCAAgc, using five classical MD simulations with a combined simulation time of approximately 120 µs. Our approach involves analyzing the tetraloop dynamics, including the folding transition state ensembles, using the energy landscape visualization method (ELViM). The ELViM is an approach that uses internal distances to compare any two conformations, allowing for a detailed description of the folding process without requiring root mean square alignment of structures. This method has previously been applied to describe the energy landscape of disordered ß-amyloid peptides and other proteins. The ELViM results in a non-linear projection of the multidimensional space, providing a comprehensive representation of the tetraloop's energy landscape. Our results reveal four distinct transition-state regions and establish the paths that lead to the folded tetraloop structure. This detailed analysis of the tetraloop's folding process has important implications for understanding RNA folding, and the ELViM approach can be used to study other biomolecules.

Toehold clipping: A mechanism for remote control of DNA strand displacement.

The ability to create stimuli-responsive DNA nanostructures has played a prominent role in dynamic DNA nanotechnology. Primary among these is the process of toehold-based strand displacement, where a nucleic acid molecule can act as a trigger to cause conformational changes in custom-designed DNA nanostructures. Here, we add another layer of control to strand displacement reactions through a 'toehold clipping' process. By designing DNA complexes with a photocleavable linker-containing toehold or an RNA toehold, we show that we can use light (UV) or enzyme (ribonuclease) to eliminate the toehold, thus preventing strand displacement reactions. We use molecular dynamics simulations to analyze the structural effects of incorporating a photocleavable linker in DNA complexes. Beyond simple DNA duplexes, we also demonstrate the toehold clipping process in a model DNA nanostructure, by designing a toehold containing double-bundle DNA tetrahedron that disassembles when an invading strand is added, but stays intact after the toehold clipping process even in the presence of the invading strand. This work is an example of combining multiple physical or molecular stimuli to provide additional remote control over DNA nanostructure reconfiguration, advances that hold potential use in biosensing, drug delivery or molecular computation.

Challenges with Simulating Modified RNA: Insights into Role and Reciprocity of Experimental and Computational Approaches.

RNA is critical to a broad spectrum of biological and viral processes. This functional diversity is a result of their dynamic nature; the variety of three-dimensional structures that they can fold into; and a host of post-transcriptional chemical modifications. While there are many experimental techniques to study the structural dynamics of biomolecules, molecular dynamics simulations (MDS) play a significant role in complementing experimental data and providing mechanistic insights. The accuracy of the results obtained from MDS is determined by the underlying physical models i.e., the force-fields, that steer the simulations. Though RNA force-fields have received a lot of attention in the last decade, they still lag compared to their protein counterparts. The chemical diversity imparted by the RNA modifications adds another layer of complexity to an already challenging problem. Insight into the effect of RNA modifications upon RNA folding and dynamics is lacking due to the insufficiency or absence of relevant experimental data. This review provides an overview of the state of MDS of modified RNA, focusing on the challenges in parameterization of RNA modifications as well as insights into relevant reference experiments necessary for their calibration.

Cotranscriptional RNA strand exchange underlies the gene regulation mechanism in a purine-sensing transcriptional riboswitch.

RNA folds cotranscriptionally to traverse out-of-equilibrium intermediate structures that are important for RNA function in the context of gene regulation. To investigate this process, here we study the structure and function of the Bacillus subtilis yxjA purine riboswitch, a transcriptional riboswitch that downregulates a nucleoside transporter in response to binding guanine. Although the aptamer and expression platform domain sequences of the yxjA riboswitch do not completely overlap, we hypothesized that a strand exchange process triggers its structural switching in response to ligand binding. In vivo fluorescence assays, structural chemical probing data and experimentally informed secondary structure modeling suggest the presence of a nascent intermediate central helix. The formation of this central helix in the absence of ligand appears to compete with both the aptamer's P1 helix and the expression platform's transcriptional terminator. All-atom molecular dynamics simulations support the hypothesis that ligand binding stabilizes the aptamer P1 helix against central helix strand invasion, thus allowing the terminator to form. These results present a potential model mechanism to explain how ligand binding can induce downstream conformational changes by influencing local strand displacement processes of intermediate folds that could be at play in multiple riboswitch classes.

Riboswitches have challenged our understanding of biological regulation for almost two decades. The ability of small molecules to bind to RNA and control gene expression offers another layer of regulation and the potential for direct action by compounds in the environment. While some riboswitches have been well studied, we lack a general understanding of how changes in RNA structure switch genetic expression from "On" to "Off". In this study, the authors propose an elegant "strand displacement" model to explain how the RNA structure shifts between "On" and "Off" states as the concentration of small molecule ligand changes. These observations help us to understand how riboswitches enable genetic decision-making. The data provide a possible general mechanism for understanding how the competition between different strand displacement outcomes can influence RNA folding. Understanding RNA folding pathways could advance the successful design of drugs that target RNA.

Structural and Binding Effects of Chemical Modifications on Thrombin Binding Aptamer (TBA).

The thrombin binding aptamer (TBA) is a promising nucleic acid-based anticoagulant. We studied the effects of chemical modifications, such as dendrimer Trebler and NHS carboxy group, on TBA with respect to its structures and thrombin binding affinity. The two dendrimer modifications were incorporated into the TBA at the 5' end and the NHS carboxy group was added into the thymine residues in the thrombin binding site of the TBA G-quadruplex (at T4, T13 and both T4/T13) using solid phase oligonucleotide synthesis. Circular dichroism (CD) spectroscopy confirmed that all of these modified TBA variants fold into a stable G-quadruplex. The binding affinity of TBA variants with thrombin was measured by surface plasmon resonance (SPR). The binding patterns and equilibrium dissociation constants (KD) of the modified TBAs are very similar to that of the native TBA. Molecular dynamics simulations studies indicate that the additional interactions or stability enhancement introduced by the modifications are minimized either by the disruption of TBA-thrombin interactions or destabilization elsewhere in the aptamer, providing a rational explanation for our experimental data. Overall, this study identifies potential positions on the TBA that can be modified without adversely affecting its structure and thrombin binding preference, which could be useful in the design and development of more functional TBA analogues.

Computationally reconstructing cotranscriptional RNA folding from experimental data reveals rearrangement of non-native folding intermediates.

The series of RNA folding events that occur during transcription can critically influence cellular RNA function. Here, we present reconstructing RNA dynamics from data (R2D2), a method to uncover details of cotranscriptional RNA folding. We model the folding of the Escherichia coli signal recognition particle (SRP) RNA and show that it requires specific local structural fluctuations within a key hairpin to engender efficient cotranscriptional conformational rearrangement into the functional structure. All-atom molecular dynamics simulations suggest that this rearrangement proceeds through an internal toehold-mediated strand-displacement mechanism, which can be disrupted with a point mutation that limits local structural fluctuations and rescued with compensating mutations that restore these fluctuations. Moreover, a cotranscriptional folding intermediate could be cleaved in vitro by recombinant E. coli RNase P, suggesting potential cotranscriptional processing. These results from experiment-guided multi-scale modeling demonstrate that even an RNA with a simple functional structure can undergo complex folding and processing during synthesis.

Terpene Chain Length Affects the Base Pairing Discrimination of S-geranyl-2-thiouridine in RNA Duplex.

Geranylation is a hydrophobic modification discovered in several bacteria tRNAs with the function of promoting codon bias during translation. However, why nature selects this C10-geranyl group remains a question. We conduct synthesis, UV-thermal denaturation, and molecular simulation studies in RNA duplexes and reveal possible reasons behind this natural selection. Among methyl-(C1), dimethylallyl-(C5), geranyl-(C10), and farnesyl-(C15) modified 2-thiouridines, only geranyl-group promotes U:G over U:A pair. Molecular simulation shows all the modified terpene groups point to the minor groove of RNA duplexes. The discrimination between U:G and U:A pairs derives from the difference in hydrogen bonding and interactions of the chain with the hydrophobic area in the minor groove. Geranyl group has perfect length to discriminate U:G and U:A pairs, whereas the others are either too long or too short to achieve the same behavior. This work indicates that geranyl group cannot be replaced by other terpene groups in promoting codon-specificity.

Hybrid DNA/RNA nanostructures with 2'-5' linkages.

Nucleic acid nanostructures with different chemical compositions have shown utility in biological applications as they provide additional assembly parameters and enhanced stability. The naturally occurring 2'-5' linkage in RNA is thought to be a prebiotic analogue and has potential use in antisense therapeutics. Here, we report the first instance of DNA/RNA motifs containing 2'-5' linkages. We synthesized and incorporated RNA strands with 2'-5' linkages into different DNA motifs with varying number of branch points (a duplex, four arm junction, double crossover motif and tensegrity triangle motif). Using experimental characterization and molecular dynamics simulations, we show that hybrid DNA/RNA nanostructures can accommodate interspersed 2'-5' linkages with relatively minor effect on the formation of these structures. Further, the modified nanostructures showed improved resistance to ribonuclease cleavage, indicating their potential use in the construction of robust drug delivery vehicles with prolonged stability in physiological conditions.

Base-pair conformational switch modulates miR-34a targeting of Sirt1 mRNA.

MicroRNAs (miRNAs) regulate the levels of translation of messenger RNAs (mRNAs). At present, the major parameter that can explain the selection of the target mRNA and the efficiency of translation repression is the base pairing between the 'seed' region of the miRNA and its counterpart mRNA1. Here we use R1ρ relaxation-dispersion nuclear magnetic resonance2 and molecular simulations3 to reveal a dynamic switch-based on the rearrangement of a single base pair in the miRNA-mRNA duplex-that elongates a weak five-base-pair seed to a complete seven-base-pair seed. This switch also causes coaxial stacking of the seed and supplementary helix fitting into human Argonaute 2 protein (Ago2), reminiscent of an active state in prokaryotic Ago4,5. Stabilizing this transient state leads to enhanced repression of the target mRNA in cells, revealing the importance of this miRNA-mRNA structure. Our observations tie together previous findings regarding the stepwise miRNA targeting process from an initial 'screening' state to an 'active' state, and unveil the role of the RNA duplex beyond the seed in Ago2.

CoSIMS: An Optimized Trajectory-Based Collision Simulator for Ion Mobility Spectrometry.

A new, multithreaded, trajectory method based software platform, CoSIMS, is revealed and compared to reference MOBCAL collision cross sections (CCS). CoSIMS employs various molecular mechanics algorithms to lessen the computational resources required to simulate thousands of buffer gas-ion collisions, including the neglect of London dispersion interactions at long distances and the removal of trajectories that insignificantly contribute to the total CCS via an ellipsoidal projection approximation. The showcased program is used to calculate the collision cross sections of carbon fullerenes, proteins, and DNA strands of various lengths, sizes, and molecular weights, and these are compared against the CCSs calculated by MOBCAL. Through this analysis, it is shown that the application of the aforementioned algorithms enables both faster and more reasonable CCS calculations than MOBCAL for highly elongated molecules such as nucleic acids; for all other molecules, CoSIMS is able to reproduce the CCSs generated by MOBCAL's trajectory method within a few percent. Overall, CoSIMS is able to calculate nearly identical CCSs as MOBCAL in nearly 2 orders of magnitude less CPU time due to the various numerical methods implemented into the software, even when run on a single CPU core.

A two-dimensional replica-exchange molecular dynamics method for simulating RNA folding using sparse experimental restraints.

We present a 2D replica exchange protocol incorporating secondary structure information to dramatically improve 3D RNA folding using molecular dynamics simulations. We show that incorporating base-pairing restraints into all-atom, explicit solvent simulations enables the accurate recapitulation of the global tertiary fold for 4 representative RNAs ranging in length from 24 to 68 nt. This method can potentially utilize base-pairing information from a wide variety of experimental inputs to predict complex RNA tertiary folds including pseudoknots, multi-loop junctions, and non-canonical interactions.

tRNA Modification Detection Using Graphene Nanopores: A Simulation Study.

There are over 100 enzyme-catalyzed modifications on transfer RNA (tRNA) molecules. The levels and identity of wobble uridine (U) modifications are affected by environmental conditions and diseased states, making wobble U detection a potential biomarker for exposures and pathological conditions. The current detection of RNA modifications requires working with nucleosides in bulk samples. Nanopore detection technology uses a single-molecule approach that has the potential to detect tRNA modifications. To evaluate the feasibility of this approach, we have performed all-atom molecular dynamics (MD) simulation studies of a five-layered graphene nanopore by localizing canonical and modified uridine nucleosides. We found that in a 1 M KCl solution with applied positive and negative biases not exceeding 2 V, nanopores can distinguish U from 5-carbonylmethyluridine (cm5U), 5-methoxycarbonylmethyluridine (mcm5U), 5-methoxycarbonylmethyl-2-thiouridine (mcm5s2U), and 5-methoxycarbonylmethyl-2'-O-methyluridine (mcm5Um) based on changes in the resistance of the nanopore. Specifically, we observed that in nanopores with dimensions less than 3 nm diameter, a localized mcm5Um and mcm5U modifications could be clearly distinguished from the canonical uridine, while the other modifications showed a modest yet detectable decrease in their respective nanopore conductance. We have compared the results between nanopores of various sizes to aid in the design, optimization, and fabrication of graphene nanopores devices for tRNA modification detection.

Advances in RNA molecular dynamics: a simulator's guide to RNA force fields.

Molecular simulations have become an essential tool for biochemical research. When they work properly, they are able to provide invaluable interpretations of experimental results and ultimately provide novel, experimentally testable predictions. Unfortunately, not all simulation models are created equal, and with inaccurate models it becomes unclear what is a bona fide prediction versus a simulation artifact. RNA models are still in their infancy compared to the many robust protein models that are widely in use, and for that reason the number of RNA force field revisions in recent years has been rapidly increasing. As there is no universally accepted 'best' RNA force field at the current time, RNA simulators must decide which one is most suited to their purposes, cognizant of its essential assumptions and their inherent strengths and weaknesses. Hopefully, armed with a better understanding of what goes inside the simulation 'black box,' RNA biochemists can devise novel experiments and provide crucial thermodynamic and structural data that will guide the development and testing of improved RNA models. WIREs RNA 2017, 8:e1396. doi: 10.1002/wrna.1396 For further resources related to this article, please visit the WIREs website.

Free-energy landscape of a hyperstable RNA tetraloop.

We report the characterization of the energy landscape and the folding/unfolding thermodynamics of a hyperstable RNA tetraloop obtained through high-performance molecular dynamics simulations at microsecond timescales. Sampling of the configurational landscape is conducted using temperature replica exchange molecular dynamics over three isochores at high, ambient, and negative pressures to determine the thermodynamic stability and the free-energy landscape of the tetraloop. The simulations reveal reversible folding/unfolding transitions of the tetraloop into the canonical A-RNA conformation and the presence of two alternative configurations, including a left-handed Z-RNA conformation and a compact purine Triplet. Increasing hydrostatic pressure shows a stabilizing effect on the A-RNA conformation and a destabilization of the left-handed Z-RNA. Our results provide a comprehensive description of the folded free-energy landscape of a hyperstable RNA tetraloop and highlight the significant advances of all-atom molecular dynamics in describing the unbiased folding of a simple RNA secondary structure motif.

Complex Thermodynamic Behavior of Single-Stranded Nucleic Acid Adsorption to Graphene Surfaces.

In just over a decade since its discovery, research on graphene has exploded due to a number of potential applications in electronics, materials, and medicine. In its water-soluble form of graphene oxide, the material has shown promise as a biosensor due to its preferential absorption of single-stranded polynucleotides and fluorescence quenching properties. The rational design of these biosensors, however, requires an improved understanding of the binding thermodynamics and ultimately a predictive model of sequence-specific binding. Toward these goals, here we directly measured the binding of nucleosides and oligonucleotides to graphene oxide nanoparticles using isothermal titration calorimetry and used the results to develop molecular models of graphene-nucleic acid interactions. We found individual nucleosides binding KD values lie in the submillimolar range with binding order of rG < rA < rC < dT < rU, while 5mer and 15mer oligonucleotides had markedly higher binding affinities in the range of micromolar and submicromolar KD values, respectively. The molecular models developed here are calibrated to quantitatively reproduce the above-mentioned experimental results. For oligonucleotides, our model predicts complex binding features such as double-stacked bases and a decrease in the fraction of graphene stacked bases with increasing oligonucleotide length until plateauing beyond ∼10-15 nucleotides. These experimental and computational results set the platform for informed design of graphene-based biosensors, further increasing their potential and application.

Using in-cell SHAPE-Seq and simulations to probe structure-function design principles of RNA transcriptional regulators.

Antisense RNA-mediated transcriptional regulators are powerful tools for controlling gene expression and creating synthetic gene networks. RNA transcriptional repressors derived from natural mechanisms called attenuators are particularly versatile, though their mechanistic complexity has made them difficult to engineer. Here we identify a new structure-function design principle for attenuators that enables the forward engineering of new RNA transcriptional repressors. Using in-cell SHAPE-Seq to characterize the structures of attenuator variants within Escherichia coli, we show that attenuator hairpins that facilitate interaction with antisense RNAs require interior loops for proper function. Molecular dynamics simulations of these attenuator variants suggest these interior loops impart structural flexibility. We further observe hairpin flexibility in the cellular structures of natural RNA mechanisms that use antisense RNA interactions to repress translation, confirming earlier results from in vitro studies. Finally, we design new transcriptional attenuators in silico using an interior loop as a structural requirement and show that they function as desired in vivo. This work establishes interior loops as an important structural element for designing synthetic RNA gene regulators. We anticipate that the coupling of experimental measurement of cellular RNA structure and function with computational modeling will enable rapid discovery of structure-function design principles for a diverse array of natural and synthetic RNA regulators.

High-resolution reversible folding of hyperstable RNA tetraloops using molecular dynamics simulations.

We report the de novo folding of three hyperstable RNA tetraloops to 1-3 Å rmsd from their experimentally determined structures using molecular dynamics simulations initialized in the unfolded state. RNA tetraloops with loop sequences UUCG, GCAA, or CUUG are hyperstable because of the formation of noncanonical loop-stabilizing interactions, and they are all faithfully reproduced to angstrom-level accuracy in replica exchange molecular dynamics simulations, including explicit solvent and ion molecules. This accuracy is accomplished using unique RNA parameters, in which biases that favor rigid, highly stacked conformations are corrected to accurately capture the inherent flexibility of ssRNA loops, accurate base stacking energetics, and purine syn-anti interconversions. In a departure from traditional quantum chemistrycentric approaches to force field optimization, our parameters are calibrated directly from thermodynamic and kinetic measurements of intra- and internucleotide structural transitions. The ability to recapitulate the signature noncanonical interactions of the three most abundant hyperstable stem loop motifs represents a significant milestone to the accurate prediction of RNA tertiary structure using unbiased all-atom molecular dynamics simulations.