Spontaneous binding of single-stranded RNAs to RRM proteins visualized by unbiased atomistic simulations with a rescaled RNA force field

How Binding Site Flexibility Promotes RNA Scanning by TbRGG2 RRM: A Molecular Dynamics Simulation Study

RNA recognition motifs (RRMs) are a key class of proteins that primarily bind single-stranded RNAs. In this study, we applied standard atomistic molecular dynamics simulations to obtain insights into the intricate binding dynamics between uridine-rich RNAs and TbRGG2 RRM using the recently developed OL3-Stafix AMBER force field, which improves the description of single-stranded RNA molecules. Complementing structural experiments that unveil a primary binding mode with a single uridine bound, our simulations uncover two supplementary binding modes in which adjacent nucleotides encroach upon the binding pocket. This leads to a unique molecular mechanism through which the TbRGG2 RRM is capable of rapidly transitioning the U-rich sequence. In contrast, the presence of non-native cytidines induces stalling and destabilization of the complex. By leveraging extensive equilibrium dynamics and a large variety of binding states, TbRGG2 RRM effectively expedites diffusion along the RNA substrate while ensuring robust selectivity for U-rich sequences despite featuring a solitary binding pocket. We further substantiate our description of the complex dynamics by simulating the fully spontaneous association process of U-rich sequences to the TbRGG2 RRM. Our study highlights the critical role of dynamics and auxiliary binding states in interface dynamics employed by RNA-binding proteins, which is not readily apparent in traditional structural studies but could represent a general type of binding strategy employed by many RNA-binding proteins.

Transcriptome-scale analysis uncovers conserved residues in the hydrophobic core of the bacterial RNA chaperone Hfq required for small regulatory RNA stability

The RNA chaperone Hfq plays crucial roles in bacterial gene expression and is a major facilitator of small regulatory RNA (sRNA) action. The toroidal architecture of the Hfq hexamer presents three well-characterized surfaces that allow it to bind sRNAs to stabilize them and engage target transcripts. Hfq-interacting sRNAs are categorized into two classes based on the surfaces they use to bind Hfq. By characterizing a systematic alanine mutant library of Hfq to identify amino acid residues that impact survival of Escherichia coli experiencing nitrogen (N) starvation, we corroborated the important role of the three RNA-binding surfaces for Hfq function. We uncovered two, previously uncharacterized, conserved residues, V22 and G34, in the hydrophobic core of Hfq, to have a profound impact on Hfq's RNA-binding activity in vivo. Transcriptome-scale analysis revealed that V22A and G34A Hfq mutants cause widespread destabilization of both sRNA classes, to the same extent as seen in bacteria devoid of Hfq. However, the alanine substitutions at these residues resulted in only modest alteration in stability and structure of Hfq. We propose that V22 and G34 have impact on Hfq function, especially critical under cellular conditions when there is an increased demand for Hfq, such as N starvation.

Can We Ever Develop an Ideal RNA Force Field? Lessons Learned from Simulations of the UUCG RNA Tetraloop and Other Systems

Molecular dynamics (MD) simulations are an important and well-established tool for investigating RNA structural dynamics, but their accuracy relies heavily on the quality of the employed force field (ff). In this work, we present a comprehensive evaluation of widely used pair-additive and polarizable RNA ffs using the challenging UUCG tetraloop (TL) benchmark system. Extensive standard MD simulations, initiated from the NMR structure of the 14-mer UUCG TL, revealed that most ffs did not maintain the native state, instead favoring alternative loop conformations. Notably, three very recent variants of pair-additive ffs, OL3CP-gHBfix21, DES-Amber, and OL3R2.7, successfully preserved the native structure over a 10 × 20 μs time scale. To further assess these ffs, we performed enhanced sampling folding simulations of the shorter 8-mer UUCG TL, starting from the single-stranded conformation. Estimated folding free energies (ΔG°fold) varied significantly among these three ffs, with values of 0.0 ± 0.6, 2.4 ± 0.8, and 7.4 ± 0.2 kcal/mol for OL3CP-gHBfix21, DES-Amber, and OL3R2.7, respectively. The ΔG°fold value predicted by the OL3CP-gHBfix21 ff was closest to experimental estimates, ranging from -1.6 to -0.7 kcal/mol. In contrast, the higher ΔG°fold values obtained using DES-Amber and OL3R2.7 were unexpected, suggesting that key interactions are inaccurately described in the folded, unfolded, or misfolded ensembles. These discrepancies led us to further test DES-Amber and OL3R2.7 ffs on additional RNA and DNA systems, where further performance issues were observed. Our results emphasize the complexity of accurately modeling RNA dynamics and suggest that creating an RNA ff capable of reliably performing across a wide range of RNA systems remains extremely challenging. In conclusion, our study provides valuable insights into the capabilities of current RNA ffs and highlights key areas for future ff development.

The Drosophila RNA binding protein Hrp48 binds a specific RNA sequence of the msl-2 mRNA 3' UTR to regulate translation

Repression of msl-2 mRNA translation is essential for viability of Drosophila melanogaster females to prevent hypertranscription of both X chromosomes. This translational control event is coordinated by the female-specific protein Sex-lethal (Sxl) which recruits the RNA binding proteins Unr and Hrp48 to the 3' untranslated region (UTR) of the msl-2 transcript and represses translation initiation. The mechanism exerted by Hrp48 during translation repression and its interaction with msl-2 are not well understood. Here we investigate the RNA binding specificity and affinity of the tandem RNA recognition motifs of Hrp48. Using NMR spectroscopy, molecular dynamics simulations and isothermal titration calorimetry, we identified the exact region of msl-2 3' UTR recognized by Hrp48. Additional biophysical experiments and translation assays give further insights into complex formation of Hrp48, Unr, Sxl and RNA. Our results show that Hrp48 binds independent of Sxl and Unr downstream of the E and F binding sites of Sxl and Unr to msl-2.

On the specificity of the recognition of m6A-RNA by YTH reader domains

Most processes of life are the result of polyvalent interactions between macromolecules, often of heterogeneous types and sizes. Frequently, the times associated with these interactions are prohibitively long for interrogation using atomistic simulations. Here, we study the recognition of N6-methylated adenine (m6A) in RNA by the reader domain YTHDC1, a prototypical, cognate pair that challenges simulations through its composition and required timescales. Simulations of RNA pentanucleotides in water reveal that the unbound state can impact (un)binding kinetics in a manner that is both model- and sequence-dependent. This is important because there are two contributions to the specificity of the recognition of the Gm6AC motif: from the sequence adjacent to the central adenine and from its methylation. Next, we establish a reductionist model consisting of an RNA trinucleotide binding to the isolated reader domain in high salt. An adaptive sampling protocol allows us to quantitatively study the dissociation of this complex. Through joint analysis of a data set including both the cognate and control sequences (GAC, Am6AA, and AAA), we derive that both contributions to specificity, sequence, and methylation, are significant and in good agreement with experimental numbers. Analysis of the kinetics suggests that flexibility in both the RNA and the YTHDC1 recognition loop leads to many low-populated unbinding pathways. This multiple-pathway mechanism might be dominant for the binding of unstructured polymers, including RNA and peptides, to proteins when their association is driven by polyvalent, electrostatic interactions.

Human CSTF2 RNA Recognition Motif Domain Binds to a U-Rich RNA Sequence through a Multistep Binding Process

The RNA recognition motif (RRM) is a conserved and ubiquitous RNA-binding domain that plays essential roles in mRNA splicing, polyadenylation, transport, and stability. RRM domains exhibit remarkable diversity in binding partners, interacting with various sequences of single- and double-stranded RNA, despite their small size and compact fold. During pre-mRNA cleavage and polyadenylation, the RRM domain from CSTF2 recognizes U- or G/U-rich RNA sequences downstream from the cleavage and polyadenylation site to regulate the process. Given the importance of alternative cleavage and polyadenylation in increasing the diversity of mRNAs, the exact mechanism of binding of RNA to the RRM of CSTF2 remains unclear, particularly in the absence of a structure of this RRM bound to a native RNA substrate. Here, we performed a series of NMR titration and spin relaxation experiments, which were complemented by paramagnetic relaxation enhancement measurements and rigid-body docking, to characterize the interactions of the CSTF2 RRM with a U-rich ligand. Our results reveal a multistep binding process involving differences in ps-ns time scale dynamics and potential structural changes, particularly in the C-terminalα-helix. These results provide insights into how the CSTF2 RRM domain binds to U-rich RNA ligands and offer a greater understanding for the molecular basis of the regulation of pre-mRNA cleavage and polyadenylation.

RNA dynamics from experimental and computational approaches

Conformational dynamics is crucial for the biological function of RNA molecules and for their potential as therapeutic targets. This meeting report outlines key "take-home" messages that emerged from the presentations and discussions during the CECAM workshop "RNA dynamics from experimental and computational approaches" in Paris, June 26-28, 2023.

Dynamics of Protein-RNA Interfaces Using All-Atom Molecular Dynamics Simulations

Facing the current challenges posed by human health diseases requires the understanding of cell machinery at a molecular level. The interplay between proteins and RNA is key for any physiological phenomenon, as well protein-RNA interactions. To understand these interactions, many experimental techniques have been developed, spanning a very wide range of spatial and temporal resolutions. In particular, the knowledge of tridimensional structures of protein-RNA complexes provides structural, mechanical, and dynamical pieces of information essential to understand their functions. To get insights into the dynamics of protein-RNA complexes, we carried out all-atom molecular dynamics simulations in explicit solvent on nine different protein-RNA complexes with different functions and interface size by taking into account the bound and unbound forms. First, we characterized structural changes upon binding and, for the RNA part, the change in the puckering. Second, we extensively analyzed the interfaces, their dynamics and structural properties, and the structural waters involved in the binding, as well as the contacts mediated by them. Based on our analysis, the interfaces rearranged during the simulation time showing alternative and stable residue-residue contacts with respect to the experimental structure.

Frontiers and Challenges of Computing ncRNAs Biogenesis, Function and Modulation

Non-coding RNAs (ncRNAs), generated from nonprotein coding DNA sequences, constitute 98-99% of the human genome. Non-coding RNAs encompass diverse functional classes, including microRNAs, small interfering RNAs, PIWI-interacting RNAs, small nuclear RNAs, small nucleolar RNAs, and long non-coding RNAs. With critical involvement in gene expression and regulation across various biological and physiopathological contexts, such as neuronal disorders, immune responses, cardiovascular diseases, and cancer, non-coding RNAs are emerging as disease biomarkers and therapeutic targets. In this review, after providing an overview of non-coding RNAs' role in cell homeostasis, we illustrate the potential and the challenges of state-of-the-art computational methods exploited to study non-coding RNAs biogenesis, function, and modulation. This can be done by directly targeting them with small molecules or by altering their expression by targeting the cellular engines underlying their biosynthesis. Drawing from applications, also taken from our work, we showcase the significance and role of computer simulations in uncovering fundamental facets of ncRNA mechanisms and modulation. This information may set the basis to advance gene modulation tools and therapeutic strategies to address unmet medical needs.

Simple Adjustment of Intranucleotide Base-Phosphate Interaction in the OL3 AMBER Force Field Improves RNA Simulations

Molecular dynamics (MD) simulations represent an established tool to study RNA molecules. The outcome of MD studies depends, however, on the quality of the force field (ff). Here we suggest a correction for the widely used AMBER OL3 ff by adding a simple adjustment of the nonbonded parameters. The reparameterization of the Lennard-Jones potential for the -H8···O5'- and -H6···O5'- atom pairs addresses an intranucleotide steric clash occurring in the type 0 base-phosphate interaction (0BPh). The nonbonded fix (NBfix) modification of 0BPh interactions (NBfix0BPh modification) was tuned via a reweighting approach and subsequently tested using an extensive set of standard and enhanced sampling simulations of both unstructured and folded RNA motifs. The modification corrects minor but visible intranucleotide clash for the anti nucleobase conformation. We observed that structural ensembles of small RNA benchmark motifs simulated with the NBfix0BPh modification provide better agreement with experiments. No side effects of the modification were observed in standard simulations of larger structured RNA motifs. We suggest that the combination of OL3 RNA ff and NBfix0BPh modification is a viable option to improve RNA MD simulations.

Protein-based delivery systems for RNA delivery

RNA-based therapeutics have emerged as promising approaches to modulate gene expression and generate therapeutic proteins or antigens capable of inducing immune responses to treat a variety of diseases, such as infectious diseases, cancers, immunologic disorders, and genetic disorders. However, the efficient delivery of RNA molecules into cells poses significant challenges due to their large molecular weight, negative charge, and susceptibility to degradation by RNase enzymes. To overcome these obstacles, viral and non-viral vectors have been developed, including lipid nanoparticles, viral vectors, proteins, dendritic macromolecules, among others. Among these carriers, protein-based delivery systems have garnered considerable attention due to their potential to address specific issues associated with nanoparticle-based systems, such as liver accumulation and immunogenicity. This review provides an overview of currently marketed RNA drugs, underscores the significance of RNA delivery vector development, delineates the essential characteristics of an ideal RNA delivery vector, and introduces existing protein carriers for RNA delivery. By offering valuable insights, this review aims to serve as a reference for the future development of protein-based delivery vectors for RNA therapeutics.

Atomistic Picture of Opening-Closing Dynamics of DNA Holliday Junction Obtained by Molecular Simulations

Holliday junction (HJ) is a noncanonical four-way DNA structure with a prominent role in DNA repair, recombination, and DNA nanotechnology. By rearranging its four arms, HJ can adopt either closed or open state. With enzymes typically recognizing only a single state, acquiring detailed knowledge of the rearrangement process is an important step toward fully understanding the biological function of HJs. Here, we carried out standard all-atom molecular dynamics (MD) simulations of the spontaneous opening-closing transitions, which revealed complex conformational transitions of HJs with an involvement of previously unconsidered "half-closed" intermediates. Detailed free-energy landscapes of the transitions were obtained by sophisticated enhanced sampling simulations. Because the force field overstabilizes the closed conformation of HJs, we developed a system-specific modification which for the first time allows the observation of spontaneous opening-closing HJ transitions in unbiased MD simulations and opens the possibilities for more accurate HJ computational studies of biological processes and nanomaterials.