Molecular Dynamics Simulations Unveil the Basis of the Sequential Binding of RsmE to the Noncoding RNA RsmZ

Dynamics and Function of sRNA/mRNAs Under the Scrutiny of Computational Simulation Methods

Molecular dynamics simulations have proved extremely useful in investigating the functioning of proteins with atomic-scale resolution. Many applications to the study of RNA also exist, and their number increases by the day. However, implementing MD simulations for RNA molecules in solution faces challenges that the MD practitioner must be aware of for the appropriate use of this tool. In this chapter, we present the fundamentals of MD simulations, in general, and the peculiarities of RNA simulations, in particular. We discuss the strengths and limitations of the technique and provide examples of its application to elucidate small RNA's performance.

Characterization of epitranscriptome reader proteins experimentally and in silico: Current knowledge and future perspectives beyond the YTH domain

To date, over 150 chemical modifications to the four canonical RNA bases have been discovered, known collectively as the epitranscriptome. Many of these modifications have been implicated in a variety of cellular processes and disease states. Additional work has been done to identify proteins known as "readers" that selectively interact with RNAs that contain specific chemical modifications. Protein interactomes with N6-methyladenosine (m⁶A), N1-methyladenosine (m¹A), N5-methylcytosine (m⁵C), and 8-oxo-7,8-dihydroguanosine (8-oxoG) have been determined, mainly through experimental advances in proteomics techniques. However, relatively few proteins have been confirmed to bind directly to RNA containing these modifications. Furthermore, for many of these protein readers, the exact binding mechanisms as well as the exclusivity for recognition of modified RNA species remain elusive, leading to questions regarding their roles within different cellular processes. In the case of the YT-521B homology (YTH) family of proteins, both experimental and in silico techniques have been leveraged to provide valuable biophysical insights into the mechanisms of m⁶A recognition at atomic resolution. To date, the YTH family is one of the best characterized classes of readers. Here, we review current knowledge about epitranscriptome recognition of the YTH domain proteins from previously published experimental and computational studies. We additionally outline knowledge gaps for proteins beyond the well-studied human YTH domains and the current in silico techniques and resources that can enable investigation of protein interactions with modified RNA outside of the YTH-m⁶A context.

On the Uses of PCA to Characterise Molecular Dynamics Simulations of Biological Macromolecules: Basics and Tips for an Effective Use

Principal Component Analysis (PCA) is a procedure widely used to examine data collected from molecular dynamics simulations of biological macromolecules. It allows for greatly reducing the dimensionality of their configurational space, facilitating further qualitative and quantitative analysis. Its simplicity and relatively low computational cost explain its extended use. However, a judicious implementation of PCA requires the knowledge of its theoretical grounds as well as its weaknesses and capabilities. In this article, we review these issues and discuss several strategies developed over the last years to mitigate the main PCA flaws and enhance the reproducibility of its results.

Recognition and Binding of RsmE to an AGGAC Motif of RsmZ: Insights from Molecular Dynamics Simulations

CsrA/RsmE is a post-transcriptional regulator protein widely distributed in bacteria. It impedes the expression of target mRNAs by attaching their 5' untranslated region. The translation is restored by small, noncoding RNAs that sequester CsrA/RsmE acting as sponges. In both cases, the protein recognizes and attaches to specific AGGAX and AXGGAX motifs, where X refers to any nucleotide. RsmZ of Pseudomonas protegens is one of these small RNAs. The structures of some of its complexes with RsmE were disclosed a few years ago. We have used umbrella sampling simulations to force the unbinding of RsmE from the AGGAC motif located in the single-stranded region sited between stem loops 2 and 3 of RsmZ. The calculations unveiled the identity of the main residues and nucleotides involved in the process. They also showed that the region adopts a hairpin-like conformation during the initial stages of the binding. The ability to acquire this conformation requires that the region has a length of at least nine nucleotides. Besides, we performed standard molecular dynamics simulations of the isolated fragments, analyzed their typical conformations, and characterized their movements. This analysis revealed that the free molecules oscillate along specific collective coordinates that facilitate the initial stages of the binding. The results strongly suggest that the flexibility of the single-stranded region of RsmZ crucially affects the ability of its binding motif to catch RsmE.