Molecular basis for the increased affinity of an RNA recognition motif with re-engineered specificity: A molecular dynamics and enhanced sampling simulations study

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.

An overview of structural approaches to study therapeutic RNAs

RNAs provide considerable opportunities as therapeutic agent to expand the plethora of classical therapeutic targets, from extracellular and surface proteins to intracellular nucleic acids and its regulators, in a wide range of diseases. RNA versatility can be exploited to recognize cell types, perform cell therapy, and develop new vaccine classes. Therapeutic RNAs (aptamers, antisense nucleotides, siRNA, miRNA, mRNA and CRISPR-Cas9) can modulate or induce protein expression, inhibit molecular interactions, achieve genome editing as well as exon-skipping. A common RNA thread, which makes it very promising for therapeutic applications, is its structure, flexibility, and binding specificity. Moreover, RNA displays peculiar structural plasticity compared to proteins as well as to DNA. Here we summarize the recent advances and applications of therapeutic RNAs, and the experimental and computational methods to analyze their structure, by biophysical techniques (liquid-state NMR, scattering, reactivity, and computational simulations), with a focus on dynamic and flexibility aspects and to binding analysis. This will provide insights on the currently available RNA therapeutic applications and on the best techniques to evaluate its dynamics and reactivity.

In vitro evolution reveals non-cationic protein-RNA interaction mediated by metal ions

RNA–peptide/protein interactions have been of utmost importance to life since its earliest forms, reaching even before the last universal common ancestor (LUCA). However, the ancient molecular mechanisms behind this key biological interaction remain enigmatic because extant RNA–protein interactions rely heavily on positively charged and aromatic amino acids that were absent (or heavily under-represented) in the early pre-LUCA evolutionary period. Here, an RNA-binding variant of the ribosomal uL11 C-terminal domain was selected from an approximately 10¹⁰ library of partially randomized sequences, all composed of ten prebiotically plausible canonical amino acids. The selected variant binds to the cognate RNA with a similar overall affinity although it is less structured in the unbound form than the wild-type protein domain. The variant complex association and dissociation are both slower than for the wild-type, implying different mechanistic processes involved. The profile of the wild-type and mutant complex stabilities along with molecular dynamics simulations uncovers qualitative differences in the interaction modes. In the absence of positively charged and aromatic residues, the mutant uL11 domain uses ion bridging (K⁺/Mg²⁺) interactions between the RNA sugar-phosphate backbone and glutamic acid residues as an alternative source of stabilization. This study presents experimental support to provide a new perspective on how early protein–RNA interactions evolved, where the lack of aromatic/basic residues may have been compensated by acidic residues plus metal ions.

Spontaneous Membrane Nanodomain Formation in the Absence or Presence of the Neurotransmitter Serotonin

Detailed knowledge on the formation of biomembrane domains, their structure, composition, and physical characteristics is scarce. Despite its frequently discussed importance in signaling, e.g., in obtaining localized non-homogeneous receptor compositions in the plasma membrane, the nanometer size as well as the dynamic and transient nature of domains impede their experimental characterization. In turn, atomistic molecular dynamics (MD) simulations combine both, high spatial and high temporal resolution. Here, using microsecond atomistic MD simulations, we characterize the spontaneous and unbiased formation of nano-domains in a plasma membrane model containing phosphatidylcholine (POPC), palmitoyl-sphingomyelin (PSM), and cholesterol (Chol) in the presence or absence of the neurotransmitter serotonin at different temperatures. In the ternary mixture, highly ordered and highly disordered domains of similar composition coexist at 303 K. The distinction of domains by lipid acyl chain order gets lost at lower temperatures of 298 and 294 K, suggesting a phase transition at ambient temperature. By comparison of domain ordering and composition, we demonstrate how the domain-specific binding of the neurotransmitter serotonin results in a modified domain lipid composition and a substantial downward shift of the phase transition temperature. Our simulations thus suggest a novel mode of action of neurotransmitters possibly of importance in neuronal signal transmission.

Molecular Dynamics Simulations with NAMD2

X-ray diffraction crystallography is the primary technique to determine the three-dimensional structures of biomolecules. Although a robust method, X-ray crystallography is not able to access the dynamical behavior of macromolecules. To do so, we have to carry out molecular dynamics simulations taking as an initial system the three-dimensional structure obtained from experimental techniques or generated using homology modeling. In this chapter, we describe in detail a tutorial to carry out molecular dynamics simulations using the program NAMD2. We chose as a molecular system to simulate the structure of human cyclin-dependent kinase 2.

Determination of a Structural Ensemble Representing the Dynamics of a G-Quadruplex DNA

It is increasingly recognized that the structures and dynamics of G-quadruplex DNA molecules are dictated by their sequences and greatly affected by environmental factors. The core guanine tetrads (G-tetrads) coordinate cations and display a strong conformational rigidity compared with that of the connecting loops. Although long loops linking the G-tetrads are typically disfavored, when present, they provide a striking illustration of the dynamics of short, single-stranded DNA regions. In addition to their role in determining the stability of the G-quadruplex state, these loops are also interesting as potential drug targets. To characterize accurately the dynamics of this DNA state, we apply here the principles of structural ensemble determination developed in the past two decades for protein molecules to DNA molecules. We thus perform extensive molecular dynamics simulations restrained with nuclear magnetic resonance residual dipolar couplings to determine a structural ensemble of the human CEB25 minisatellite G-quadruplex, which contains a connecting loop of nine nucleotides. This structural ensemble displays a wide set of arrangements for the loop and a compact, well-defined G-quadruplex core. Our results show the importance of stacking interactions in the loop and strengthen the ability of the closing base pairs to confer a large thermodynamic stability to the G-quadruplex structure.

Disclosing the Impact of Carcinogenic SF3b Mutations on Pre-mRNA Recognition Via All-Atom Simulations

The spliceosome accurately promotes precursor messenger-RNA splicing by recognizing specific noncoding intronic tracts including the branch point sequence (BPS) and the 3'-splice-site (3'SS). Mutations of Hsh155 (yeast)/SF3B1 (human), which is a protein of the SF3b factor involved in BPS recognition and induces altered BPS binding and 3'SS selection, lead to mis-spliced mRNA transcripts. Although these mutations recur in hematologic malignancies, the mechanism by which they change gene expression remains unclear. In this study, multi-microsecond-long molecular-dynamics simulations of eighth distinct ∼700,000 atom models of the spliceosome Bact complex, and gene sequencing of SF3B1, disclose that these carcinogenic isoforms destabilize intron binding and/or affect the functional dynamics of Hsh155/SF3B1 only when binding non-consensus BPSs, as opposed to the non-pathogenic variants newly annotated here. This pinpoints a cross-talk between the distal Hsh155 mutation and BPS recognition sites. Our outcomes unprecedentedly contribute to elucidating the principles of pre-mRNA recognition, which provides critical insights on the mechanism underlying constitutive/alternative/aberrant splicing.