The dynamics of t1 adenosine binding on human Argonaute 2: Understanding recognition with conformational selection

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.

Probing the conformational dynamics of an Ago-RNA complex in water/methanol solution

Argonaute (Ago) proteins mediate target recognition guiding miRNA to bind complementary mRNA primarily in the seed region. However, additional pairing can occur beyond the seed, forming a supplementary duplex that can contribute to the guide-target affinity. In order to shed light on the connection, between protein-RNA interactions and miRNA-mRNA seed and supplementary duplex mobility, we carried out molecular dynamics simulations at the microsecond time-scale using a different approach compared to the ones normally used. Until now, theoretical investigations with classical MD on Ago-RNA complexes have been focused primarily on pure water solvent, which mimics the natural environment of biological molecules. Here, we explored the conformational space of a human Ago2 (hAgo2) bound to the seed + supplementary miRNA-mRNA duplex, using the solvent environment as a molecular probe. MD simulations have been performed in a mixture of water/MeOH at a molar ratio of 70 : 30 as well as in pure water for comparison. Our findings revealed that the mixed solvent promotes protein RNA association, principally enhancing salt-linkages between basic amino acid side-chains and acidic phosphates of the sugar-phosphate backbone. The primary effect registered was the restriction of supplementary duplex flexibility and the stabilization of the miRNA 3' terminus. Interestingly, we observed that the influence of the solvent appears to have almost no impact on the conformation of the seed duplex.

RNA-protein complexes and force field polarizability

Molecular dynamic (MD) simulations offer a way to study biomolecular interactions and their dynamics at the atomistic level. There are only a few studies of RNA-protein complexes in MD simulations, and here we wanted to study how force fields differ when simulating RNA-protein complexes: 1) argonaute 2 with bound guide RNA and a target RNA, 2) CasPhi-2 bound to CRISPR RNA and 3) Retinoic acid-inducible gene I C268F variant in complex with double-stranded RNA. We tested three non-polarizable force fields: Amber protein force fields ff14SB and ff19SB with RNA force field OL3, and the all-atom OPLS4 force field. Due to the highly charged and polar nature of RNA, we also tested the polarizable AMOEBA force field and the ff19SB and OL3 force fields with a polarizable water model O3P. Our results show that the non-polarizable force fields lead to compact and stable complexes. The polarizability in the force field or in the water model allows significantly more movement from the complex, but in some cases, this results in the disintegration of the complex structure, especially if the protein contains longer loop regions. Thus, one should be cautious when running long-scale simulations with polarizability. As a conclusion, all the tested force fields can be used to simulate RNA-protein complexes and the choice of the optimal force field depends on the studied system and research question.