Mapping L1 ligase ribozyme conformational switch

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Prebiotic Factors Influencing the Activity of a Ligase Ribozyme

An RNA-lipid origin of life scenario provides a plausible route for compartmentalized replication of an informational polymer and subsequent division of the container. However, a full narrative to form such RNA protocells implies that catalytic RNA molecules, called ribozymes, can operate in the presence of self-assembled vesicles composed of prebiotically relevant constituents, such as fatty acids. Hereby, we subjected a newly engineered truncated variant of the L1 ligase ribozyme, named tL1, to various environmental conditions that may have prevailed on the early Earth with the objective to find a set of control parameters enabling both tL1-catalyzed ligation and formation of stable myristoleic acid (MA) vesicles. The separate and concurrent effects of temperature, concentrations of Mg²⁺, MA, polyethylene glycol and various solutes were investigated. The most favorable condition tested consists of 100 mM NaCl, 1 mM Mg²⁺, 5 mM MA, and 4 °C temperature, whereas the addition of Mg²⁺-chelating solutes, such as citrate, tRNAs, aspartic acid, and nucleoside triphosphates severely inhibits the reaction. These results further solidify the RNA-lipid world hypothesis and stress the importance of using a systems chemistry approach whereby a wide range of prebiotic factors interfacing with ribozymes are considered.

Multiscale methods for computational RNA enzymology

RNA catalysis is of fundamental importance to biology and yet remains ill-understood due to its complex nature. The multidimensional "problem space" of RNA catalysis includes both local and global conformational rearrangements, changes in the ion atmosphere around nucleic acids and metal ion binding, dependence on potentially correlated protonation states of key residues, and bond breaking/forming in the chemical steps of the reaction. The goal of this chapter is to summarize and apply multiscale modeling methods in an effort to target the different parts of the RNA catalysis problem space while also addressing the limitations and pitfalls of these methods. Classical molecular dynamics simulations, reference interaction site model calculations, constant pH molecular dynamics (CpHMD) simulations, Hamiltonian replica exchange molecular dynamics, and quantum mechanical/molecular mechanical simulations will be discussed in the context of the study of RNA backbone cleavage transesterification. This reaction is catalyzed by both RNA and protein enzymes, and here we examine the different mechanistic strategies taken by the hepatitis delta virus ribozyme and RNase A.

Efficient in silico exploration of RNA interhelical conformations using Euler angles and WExplore

HIV-1 TAR RNA is a two-helix bulge motif that plays a critical role in HIV viral replication and is an important drug target. However, efforts at designing TAR inhibitors have been challenged by its high degree of structural flexibility, which includes slow large-amplitude reorientations of its helices with respect to one another. Here, we use the recently introduced algorithm WExplore in combination with Euler angles to achieve unprecedented sampling of the TAR conformational ensemble. Our ensemble achieves similar agreement with experimental NMR data when compared with previous TAR computational studies, and is generated at a fraction of the computational cost. It clearly emerges from configuration space network analysis that the intermittent formation of the A22-U40 base pair acts as a reversible switch that enables sampling of interhelical conformations that would otherwise be topologically disallowed. We find that most previously determined ligand-bound structures are found in similar location in the network, and we use a sample-and-select approach to guide the construction of a set of novel conformations which can serve as the basis for future drug development efforts. Collectively, our findings demonstrate the utility of WExplore in combination with suitable order parameters as a method for exploring RNA conformational space.

The spatial structure of cell signaling systems

The spatial structure of the cell is highly organized at all levels: from small complexes and assemblies, to local nano- and microclusters, to global, micrometer scales across and between cells. We suggest that this multiscale spatial cell organization also organizes signaling and coordinates cellular behavior. We propose a new view of the spatial structure of cell signaling systems. This new view describes cell signaling in terms of dynamic allosteric interactions within and among distinct, spatially organized transient clusters. The clusters vary over time and space and are on length scales from nanometers to micrometers. When considered across these length scales, primary factors in the spatial organization are cell membrane domains and the actin cytoskeleton, both also highly dynamic. A key challenge is to understand the interplay across these multiple scales, link it to the physicochemical basis of the conformational behavior of single molecules and ultimately relate it to cellular function. Overall, our premise is that at these scales, cell signaling should be thought of not primarily as a sequence of diffusion-controlled molecular collisions, but instead transient, allostery-driven cluster re-forming interactions.

Native states of fast-folding proteins are kinetic traps

It has been suggested that the native state of a protein acts as a kinetic hub that can facilitate transitions between nonnative states. Using recently developed tools to quantify mediation probabilities ("hub scores"), we quantify hub-like behavior in atomic resolution trajectories for the first time. We use a data set of trajectory ensembles for 12 fast-folding proteins previously published by D. E. Shaw Research (Lindorff-Larsen, K.; et al. How Fast-Folding Proteins Fold. Science2011, 334, 517) with an aggregate simulation time of over 8.2 ms. We visualize the free-energy landscape of each molecule using configuration space networks, and show that dynamic quantities can be qualitatively understood from visual inspection of the networks. Modularity optimization is used to provide a parameter-free means of tessellating the network into a group of communities. Using hub scores, we find that the percentage of trajectories that are mediated by the native state is 31% when averaged over all molecules, and reaches a maximum of 52% for the homeodomain and chignolin. Furthermore, for these mediated transitions, we use Markov models to determine whether the native state acts as a facilitator for the transition, or as a trap (i.e., an off-pathway detour). Although instances of facilitation are found in 4 of the 12 molecules, we conclude that the native state acts primarily as a trap, which is consistent with the idea of a funnel-like landscape.