A short guide for molecular dynamics simulations of RNA systems

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.

Structural, Dynamical, and Entropic Differences between SARS-CoV and SARS-CoV-2 s2m Elements Using Molecular Dynamics Simulations

The functional role of the highly conserved stem-loop II motif (s2m) in SARS-CoV and SARS-CoV-2 in the viral lifecycle remains enigmatic and an intense area of research. Structure and dynamics of the s2m are key to establishing a structure-function connection, yet a full set of atomistic resolution coordinates is not available for SARS-CoV-2 s2m. Our work constructs three-dimensional coordinates consistent with NMR solution phase data for SARS-CoV-2 s2m and provides a comparative analysis with its counterpart SARS-CoV s2m. We employed initial coordinates based on PDB ID 1XJR for SARS-CoV s2m and two models for SARS-CoV-2 s2m: one based on 1XJR in which we introduced the mutations present in SARS-CoV-2 s2m and the second based on the available SARS-CoV-2 NMR NOE data supplemented with knowledge-based methods. For each of the three systems, 3.5 μs molecular dynamics simulations were used to sample the structure and dynamics, and principal component analysis (PCA) reduced the ensembles to hierarchal conformational substates for detailed analysis. Dilute solution simulations of SARS-CoV s2m demonstrate that the GNRA-like terminal pentaloop is rigidly defined by base stacking uniquely positioned for possible kissing dimer formation. However, the SARS-CoV-2 s2m simulation did not retain the reported crystallographic SARS-CoV motifs and the terminal loop expands to a highly dynamic "nonaloop." Increased flexibility and structural disorganization are observed for the larger terminal loop, where an entropic penalty is computed to explain the experimentally observed reduction in kissing complex formation. Overall, both SARS-CoV and SARS-CoV-2 s2m elements have a similarly pronounced L-shape due to different motif interactions. Our study establishes the atomistic three-dimensional structure and uncovers dynamic differences that arise from s2m sequence changes, which sets the stage for the interrogation of different mechanistic pathways of suspected biological function.

Protein-Based Virtual Screening Tools Applied for RNA-Ligand Docking Identify New Binders of the preQ1-Riboswitch

Targeting RNA with small molecules is an emerging field. While several ligands for different RNA targets are reported, structure-based virtual screenings (VSs) against RNAs are still rare. Here, we elucidated the general capabilities of protein-based docking programs to reproduce native binding modes of small-molecule RNA ligands and to discriminate known binders from decoys by the scoring function. The programs were found to perform similar compared to the RNA-based docking tool rDOCK, and the challenges faced during docking, namely, protomer and tautomer selection, target dynamics, and explicit solvent, do not largely differ from challenges in conventional protein-ligand docking. A prospective VS with the Bacillus subtilis preQ₁-riboswitch aptamer domain performed with FRED, HYBRID, and FlexX followed by microscale thermophoresis assays identified six active compounds out of 23 tested VS hits with potencies between 29.5 nM and 11.0 μM. The hits were selected not solely based on their docking score but for resembling key interactions of the native ligand. Therefore, this study demonstrates the general feasibility to perform structure-based VSs against RNA targets, while at the same time it highlights pitfalls and their potential solutions when executing RNA-ligand docking.

The effect of mutation in the stem of the MicroROSE thermometer on its thermosensing ability: insights from molecular dynamics simulation studies

Mutation induced thermosensing ability of MicroROSE thermometer.

Understanding the mechanistic basis of non-coding RNA through molecular dynamics simulations

Noncoding RNA (ncRNA) has a key role in regulating gene expression, mediating fundamental processes and diseases via a variety of yet unknown mechanisms. Here, we review recent applications of conventional and enhanced Molecular Dynamics (MD) simulations methods to address the mechanistic function of large biomolecular systems that are tightly involved in the ncRNA function and that are of key importance in life sciences. This compendium focuses of three biomolecular systems, namely the CRISPR-Cas9 genome editing machinery, group II intron ribozyme and the ribonucleoprotein complex of the spliceosome, which edit and process ncRNA. We show how the application of a novel accelerated MD simulations method has been key in disclosing the conformational transitions underlying RNA binding in the CRISPR-Cas9 complex, suggesting a mechanism for RNA recruitment and clarifying the conformational changes required for attaining genome editing. As well, we discuss the use of mixed quantum-classical MD simulations in deciphering the catalytic mechanism of RNA splicing as operated by group II intron ribozyme, one of the largest ncRNA structures crystallized so far. Finally, we debate the future challenges and opportunities in the field, discussing the recent application of MD simulations for unraveling the functional biophysics of the spliceosome, a multi-mega Dalton complex of proteins and small nuclear RNAs that performs RNA splicing in humans. This showcase of applications highlights the current talent of MD simulations to dissect atomic-level details of complex biomolecular systems instrumental for the design of finely engineered genome editing machines. As well, this review aims at inspiring future investigations of several other ncRNA regulatory systems, such as micro and small interfering RNAs, which achieve their function and specificity using RNA-based recognition and targeting strategies.

Nucleobase carbonyl groups are poor Mg2+ inner-sphere binders but excellent monovalent ion binders-a critical PDB survey

Precise knowledge of Mg²⁺ inner-sphere binding site properties is vital for understanding the structure and function of nucleic acid systems. Unfortunately, the PDB, which represents the main source of Mg²⁺ binding sites, contains a substantial number of assignment issues that blur our understanding of the functions of these ions. Here, following a previous study devoted to Mg²⁺ binding to nucleobase nitrogens, we surveyed nucleic acid X-ray structures from the PDB with resolutions ≤2.9 Å to classify the Mg²⁺ inner-sphere binding patterns to nucleotide carbonyl, ribose hydroxyl, cyclic ether, and phosphodiester oxygen atoms. From this classification, we derived a set of "prior-knowledge" nucleobase Mg²⁺ binding sites. We report that crystallographic examples of trustworthy nucleobase Mg²⁺ binding sites are fewer than expected since many of those are associated with misidentified Na⁺ or K⁺ We also emphasize that binding of Na⁺ and K⁺ to nucleic acids is much more frequent than anticipated. Overall, we provide evidence derived from X-ray structures that nucleobases are poor inner-sphere binders for Mg²⁺ but good binders for monovalent ions. Based on strict stereochemical criteria, we propose an extended set of guidelines designed to help in the assignment and validation of ions directly contacting nucleobase and ribose atoms. These guidelines should help in the interpretation of X-ray and cryo-EM solvent density maps. When borderline Mg²⁺ stereochemistry is observed, alternative placement of Na⁺, K⁺, or Ca²⁺ must be considered. We also critically examine the use of lanthanides (Yb³⁺, Tb³⁺) as Mg²⁺ substitutes in crystallography experiments.

Dynamic coarse-graining fills the gap between atomistic simulations and experimental investigations of mechanical unfolding

We present a dynamic coarse-graining technique that allows one to simulate the mechanical unfolding of biomolecules or molecular complexes on experimentally relevant time scales. It is based on Markov state models (MSMs), which we construct from molecular dynamics simulations using the pulling coordinate as an order parameter. We obtain a sequence of MSMs as a function of the discretized pulling coordinate, and the pulling process is modeled by switching among the MSMs according to the protocol applied to unfold the complex. This way we cover seven orders of magnitude in pulling speed. In the region of rapid pulling, we additionally perform steered molecular dynamics simulations and find excellent agreement between the results of the fully atomistic and the dynamically coarse-grained simulations. Our technique allows the determination of the rates of mechanical unfolding in a dynamical range from approximately 10^-8/ns to 1/ns thus reaching experimentally accessible time regimes without abandoning atomistic resolution.

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.

Structural study of the Fox-1 RRM protein hydration reveals a role for key water molecules in RRM-RNA recognition

The Fox-1 RNA recognition motif (RRM) domain is an important member of the RRM protein family. We report a 1.8 Å X-ray structure of the free Fox-1 containing six distinct monomers. We use this and the nuclear magnetic resonance (NMR) structure of the Fox-1 protein/RNA complex for molecular dynamics (MD) analyses of the structured hydration. The individual monomers of the X-ray structure show diverse hydration patterns, however, MD excellently reproduces the most occupied hydration sites. Simulations of the protein/RNA complex show hydration consistent with the isolated protein complemented by hydration sites specific to the protein/RNA interface. MD predicts intricate hydration sites with water-binding times extending up to hundreds of nanoseconds. We characterize two of them using NMR spectroscopy, RNA binding with switchSENSE and free-energy calculations of mutant proteins. Both hydration sites are experimentally confirmed and their abolishment reduces the binding free-energy. A quantitative agreement between theory and experiment is achieved for the S155A substitution but not for the S122A mutant. The S155 hydration site is evolutionarily conserved within the RRM domains. In conclusion, MD is an effective tool for predicting and interpreting the hydration patterns of protein/RNA complexes. Hydration is not easily detectable in NMR experiments but can affect stability of protein/RNA complexes.

Folding of guanine quadruplex molecules-funnel-like mechanism or kinetic partitioning? An overview from MD simulation studies

BACKGROUND

Guanine quadruplexes (GQs) play vital roles in many cellular processes and are of much interest as drug targets. In contrast to the availability of many structural studies, there is still limited knowledge on GQ folding.

SCOPE OF REVIEW

We review recent molecular dynamics (MD) simulation studies of the folding of GQs, with an emphasis paid to the human telomeric DNA GQ. We explain the basic principles and limitations of all types of MD methods used to study unfolding and folding in a way accessible to non-specialists. We discuss the potential role of G-hairpin, G-triplex and alternative GQ intermediates in the folding process. We argue that, in general, folding of GQs is fundamentally different from funneled folding of small fast-folding proteins, and can be best described by a kinetic partitioning (KP) mechanism. KP is a competition between at least two (but often many) well-separated and structurally different conformational ensembles.

MAJOR CONCLUSIONS

The KP mechanism is the only plausible way to explain experiments reporting long time-scales of GQ folding and the existence of long-lived sub-states. A significant part of the natural partitioning of the free energy landscape of GQs comes from the ability of the GQ-forming sequences to populate a large number of syn-anti patterns in their G-tracts. The extreme complexity of the KP of GQs typically prevents an appropriate description of the folding landscape using just a few order parameters or collective variables.

GENERAL SIGNIFICANCE

We reconcile available computational and experimental studies of GQ folding and formulate basic principles characterizing GQ folding landscapes. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.

How to understand atomistic molecular dynamics simulations of RNA and protein-RNA complexes?

We provide a critical assessment of explicit-solvent atomistic molecular dynamics (MD) simulations of RNA and protein/RNA complexes, written primarily for non-specialists with an emphasis to explain the limitations of MD. MD simulations can be likened to hypothetical single-molecule experiments starting from single atomistic conformations and investigating genuine thermal sampling of the biomolecules. The main advantage of MD is the unlimited temporal and spatial resolution of positions of all atoms in the simulated systems. Fundamental limitations are the short physical time-scale of simulations, which can be partially alleviated by enhanced-sampling techniques, and the highly approximate atomistic force fields describing the simulated molecules. The applicability and present limitations of MD are demonstrated on studies of tetranucleotides, tetraloops, ribozymes, riboswitches and protein/RNA complexes. Wisely applied simulations respecting the approximations of the model can successfully complement structural and biochemical experiments. WIREs RNA 2017, 8:e1405. doi: 10.1002/wrna.1405 For further resources related to this article, please visit the WIREs website.

Small molecules targeting viral RNA

Highly conserved noncoding RNA (ncRNA) elements in viral genomes and transcripts offer new opportunities to expand the repertoire of drug targets for the development of antiinfective therapy. Ligands binding to ncRNA architectures are able to affect interactions, structural stability or conformational changes and thereby block processes essential for viral replication. Proof of concept for targeting functional RNA by small molecule inhibitors has been demonstrated for multiple viruses with RNA genomes. Strategies to identify antiviral compounds as inhibitors of ncRNA are increasingly emphasizing consideration of drug‐like properties of candidate molecules emerging from screening and ligand design. Recent efforts of antiviral lead discovery for RNA targets have provided drug‐like small molecules that inhibit viral replication and include inhibitors of human immunodeficiency virus (HIV), hepatitis C virus (HCV), severe respiratory syndrome coronavirus (SARS CoV), and influenza A virus. While target selectivity remains a challenge for the discovery of useful RNA‐binding compounds, a better understanding is emerging of properties that define RNA targets amenable for inhibition by small molecule ligands. Insight from successful approaches of targeting viral ncRNA in HIV, HCV, SARS CoV, and influenza A will provide a basis for the future exploration of RNA targets for therapeutic intervention in other viral pathogens which create urgent, unmet medical needs. Viruses for which targeting ncRNA components in the genome or transcripts may be promising include insect‐borne flaviviruses (Dengue, Zika, and West Nile) and filoviruses (Ebola and Marburg). WIREs RNA 2016, 7:726–743. doi: 10.1002/wrna.1373

This article is categorized under: 1

RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems

2

RNA Interactions with Proteins and Other Molecules > Small Molecule–RNA Interactions

3

Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs

Elucidating the Relation between Internal Motions and Dihedral Angles in an RNA Hairpin Using Molecular Dynamics

Molecular dynamics simulations were performed to characterize the internal motions of the ribonucleic acid apical stem loop of human hepatitis B virus. The NMR relaxation rates calculated directly from the trajectory are in good agreement with the experiment. Calculated order parameters follow the experimental pattern. Order parameters lower than 0.8 are observed for nucleotides that are weakly hydrogen bonded to their base pair partner, unpaired, or part of the loop. These residues show slow decay of the internal correlation functions of their base and sugar C-H vectors. Concerted motions around backbone dihedral angles influence the amplitude of motion of the sugar and base C-H vectors. The order parameters for base C-H vectors are also affected by the fluctuation of the glycosidic dihedral angle.

Ab initio RNA folding

RNA molecules are essential cellular machines performing a wide variety of functions for which a specific three-dimensional structure is required. Over the last several years, the experimental determination of RNA structures through x-ray crystallography and NMR seems to have reached a plateau in the number of structures resolved each year, but as more and more RNA sequences are being discovered, the need for structure prediction tools to complement experimental data is strong. Theoretical approaches to RNA folding have been developed since the late nineties, when the first algorithms for secondary structure prediction appeared. Over the last 10 years a number of prediction methods for 3D structures have been developed, first based on bioinformatics and data-mining, and more recently based on a coarse-grained physical representation of the systems. In this review we are going to present the challenges of RNA structure prediction and the main ideas behind bioinformatic approaches and physics-based approaches. We will focus on the description of the more recent physics-based phenomenological models and on how they are built to include the specificity of the interactions of RNA bases, whose role is critical in folding. Through examples from different models, we will point out the strengths of physics-based approaches, which are able not only to predict equilibrium structures, but also to investigate dynamical and thermodynamical behavior, and the open challenges to include more key interactions ruling RNA folding.

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.

Molecular dynamic simulations of protein/RNA complexes: CRISPR/Csy4 endoribonuclease

BACKGROUND

Many prokaryotic genomes comprise Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) offering defense against foreign nucleic acids. These immune systems are conditioned by the production of small CRISPR-derived RNAs matured from long RNA precursors. This often requires a Csy4 endoribonuclease cleaving the RNA 3'-end.

METHODS

We report extended explicit solvent molecular dynamic (MD) simulations of Csy4/RNA complex in precursor and product states, based on X-ray structures of product and inactivated precursor (55 simulations; ~3.7μs in total).

RESULTS

The simulations identify double-protonated His29 and deprotonated terminal phosphate as the likely dominant protonation states consistent with the product structure. We revealed potential substates consistent with Ser148 and His29 acting as the general base and acid, respectively. The Ser148 could be straightforwardly deprotonated through solvent and could without further structural rearrangements deprotonate the nucleophile, contrasting similar studies investigating the general base role of nucleobases in ribozymes. We could not locate geometries consistent with His29 acting as general base. However, we caution that the X-ray structures do not always capture the catalytically active geometries and then the reactive structures may be unreachable by the simulation technique.

CONCLUSIONS

We identified potential catalytic arrangement of the Csy4/RNA complex but we also report limitations of the simulation technique. Even for the dominant protonation state we could not achieve full agreement between the simulations and the structural data.

GENERAL SIGNIFICANCE

Potential catalytic arrangement of the Csy4/RNA complex is found. Further, we provide unique insights into limitations of simulations of protein/RNA complexes, namely, the influence of the starting experimental structures and force field limitations. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Thermostability of in vitro evolved Bacillus subtilis lipase A: a network and dynamics perspective

Proteins in thermophilic organisms remain stable and function optimally at high temperatures. Owing to their important applicability in many industrial processes, such thermostable proteins have been studied extensively, and several structural factors attributed to their enhanced stability. How these factors render the emergent property of thermostability to proteins, even in situations where no significant changes occur in their three-dimensional structures in comparison to their mesophilic counter-parts, has remained an intriguing question. In this study we treat Lipase A from Bacillus subtilis and its six thermostable mutants in a unified manner and address the problem with a combined complex network-based analysis and molecular dynamic studies to find commonality in their properties. The Protein Contact Networks (PCN) of the wild-type and six mutant Lipase A structures developed at a mesoscopic scale were analyzed at global network and local node (residue) level using network parameters and community structure analysis. The comparative PCN analysis of all proteins pointed towards important role of specific residues in the enhanced thermostability. Network analysis results were corroborated with finer-scale molecular dynamics simulations at both room and high temperatures. Our results show that this combined approach at two scales can uncover small but important changes in the local conformations that add up to stabilize the protein structure in thermostable mutants, even when overall conformation differences among them are negligible. Our analysis not only supports the experimentally determined stabilizing factors, but also unveils the important role of contacts, distributed throughout the protein, that lead to thermostability. We propose that this combined mesoscopic-network and fine-grained molecular dynamics approach is a convenient and useful scheme not only to study allosteric changes leading to protein stability in the face of negligible over-all conformational changes due to mutations, but also in other molecular networks where change in function does not accompany significant change in the network structure.

Molecular Dynamics Simulations of Nucleic Acids. From Tetranucleotides to the Ribosome

We present a brief overview of explicit solvent molecular dynamics (MD) simulations of nucleic acids. We explain physical chemistry limitations of the simulations, namely, the molecular mechanics (MM) force field (FF) approximation and limited time scale. Further, we discuss relations and differences between simulations and experiments, compare standard and enhanced sampling simulations, discuss the role of starting structures, comment on different versions of nucleic acid FFs, and relate MM computations with contemporary quantum chemistry. Despite its limitations, we show that MD is a powerful technique for studying the structural dynamics of nucleic acids with a fast growing potential that substantially complements experimental results and aids their interpretation.

A novel implicit solvent model for simulating the molecular dynamics of RNA

Although molecular dynamics simulations can be accelerated by more than an order of magnitude by implicitly describing the influence of the solvent with a continuum model, most currently available implicit solvent simulations cannot robustly simulate the structure and dynamics of nucleic acids. The difficulties become exacerbated especially for RNAs, suggesting the presence of serious physical flaws in the prior continuum models for the influence of the solvent and counter ions on the nucleic acids. We present a novel, to our knowledge, implicit solvent model for simulating nucleic acids by combining the Langevin-Debye model and the Poisson-Boltzmann equation to provide a better estimate of the electrostatic screening of both the water and counter ions. Tests of the model involve comparisons of implicit and explicit solvent simulations for three RNA targets with 20, 29, and 75 nucleotides. The model provides reasonable agreement with explicit solvent simulations, and directions for future improvement are noted.

Prediction of the interaction between spermidine and the G-G mismatch containing acceptor stem in tRNA(Ile): molecular modeling, density functional theory, and molecular dynamics study

Polyamines, putrescine, spermidine (SPD), and spermine are closely linked to cell growth, and highly regulate the levels of transcription, translation and protein turnover. We propose that SPD stimulates the formation of Ile-tRNA(Ile) by inducing a selective structural change of the G-G mismatch containing acceptor stem in tRNA(Ile). Here, we provide insight into how SPD recognizes and stabilizes the G-G mismatch containing acceptor stem in tRNA(Ile) with molecular modeling (MM), density functional theory (DFT) calculations, and molecular dynamics (MD) simulations. The results of the MM and DFT calculations indicate that the negatively charged region of the G-G mismatch containing acceptor stem in tRNA(Ile) is preferentially recognized by positively charged SPD. In addition, MD simulations indicate that all of the positively charged amino groups of SPD under physiological conditions (N1(NH3(+)), N5(NH2(+)), and N10(NH3(+)) could form hydrogen bonds with tRNA(Ile) and trigger the SPD-induced stabilization and structural change of the G-G mismatch containing acceptor stem in tRNA(Ile). Thus, this approach should be useful for determining the preferential binding site and appropriate binding mode of polyamines on tRNA(Ile).

The impact of a ligand binding on strand migration in the SAM-I riboswitch

Riboswitches sense cellular concentrations of small molecules and use this information to adjust synthesis rates of related metabolites. Riboswitches include an aptamer domain to detect the ligand and an expression platform to control gene expression. Previous structural studies of riboswitches largely focused on aptamers, truncating the expression domain to suppress conformational switching. To link ligand/aptamer binding to conformational switching, we constructed models of an S-adenosyl methionine (SAM)-I riboswitch RNA segment incorporating elements of the expression platform, allowing formation of an antiterminator (AT) helix. Using Anton, a computer specially developed for long timescale Molecular Dynamics (MD), we simulated an extended (three microseconds) MD trajectory with SAM bound to a modeled riboswitch RNA segment. Remarkably, we observed a strand migration, converting three base pairs from an antiterminator (AT) helix, characteristic of the transcription ON state, to a P1 helix, characteristic of the OFF state. This conformational switching towards the OFF state is observed only in the presence of SAM. Among seven extended trajectories with three starting structures, the presence of SAM enhances the trend towards the OFF state for two out of three starting structures tested. Our simulation provides a visual demonstration of how a small molecule (<500 MW) binding to a limited surface can trigger a large scale conformational rearrangement in a 40 kDa RNA by perturbing the Free Energy Landscape. Such a mechanism can explain minimal requirements for SAM binding and transcription termination for SAM-I riboswitches previously reported experimentally.

Influence of the solvent on the stability of bis(terpyridine) structures on graphite

The effect of solvation on the adsorption of organic molecules on graphite at room temperature has been addressed with force-field molecular dynamics simulations. As a model system, the solvation of a bis(terpyridine) isomer in water and 1,2,4-trichlorobenzene was studied with an explicit solvation model. The inclusion of solvation has a noticeable effect on adsorption energies. Although the results of the various considered force fields differ quite significantly, they all agree that the adsorption of BTP from the TCB solvent is almost thermoneutral. The substrate simply acts as a template to allow a planar arrangement of the network, which is stabilized by the intermolecular interaction. Using an atomic thermodynamics approach, the order of the stability of various network structures as a function of the chemical potential is derived yielding a sequence in agreement with the experiment.

Investigating the malleability of RNA aptamers

Aptamers are short, single-stranded nucleic acids with structures that frequently change upon ligand binding and are sensitive to the ionic environment. To achieve facile application of aptamers in controlling cellular activities, a better understanding is needed of aptamer ligand binding parameters, structures, intramolecular mobilities and how these structures adapt to different ionic environments with consequent effects on their ligand binding characteristics. Here we discuss the integration of biochemical analysis with NMR spectroscopy and computational modeling to explore the relation between ligand binding and structural malleability of some well-studied aptamers. Several methods for determining aptamer binding affinity and specificity are discussed, including isothermal titration calorimetry, steady state fluorescence of 2-aminopurine substituted aptamers, and dye displacement assays. Also considered are aspects of molecular dynamics simulations specific to aptamers including adding ions and simulating aptamer structure in the absence of ligand when NMR spectroscopy or X-ray crystallography structures of the unoccupied aptamer are not available. We focus specifically on RNA aptamers that bind small molecule ligands as would be applied in sensors or integrated into riboswitches such as to measure the products of metabolic activity.

Structural principles of RNA catalysis in a 2'-5' lariat-forming ribozyme

RNA-catalyzed lariat formation is present in both eukaryotes and prokaryotes. To date we lack structural insights into the catalytic mechanism of lariat-forming ribozymes. Here, we study an artificial 2'-5' AG1 lariat-forming ribozyme that shares the sequence specificity of lariat formation with the pre-mRNA splicing reaction. Using NMR, we solve the structure of the inactive state of the ribozyme in the absence of magnesium. The reaction center 5'-guanosine appears to be part of a helix with an exceptionally widened major groove, while the lariat-forming A48 is looped out at the apex of a pseudoknot. The model of the active state built by mutational analysis, molecular modeling, and small-angle X-ray scattering suggests that A48 is recognized by a conserved adenosine, juxtaposed to the 5'-guanosine in one base-pair step distance, while the G1-N7 coordinates a magnesium ion essential for the activation of the nucleophile. Our findings offer implications for lariat formation in RNA enzymes including the mechanism of the recognition of the branch-site adenosine.

Bridging the gap between theory and experiment to derive a detailed understanding of hammerhead ribozyme catalysis

Herein we summarize our progress toward the understanding of hammerhead ribozyme (HHR) catalysis through a multiscale simulation strategy. Simulation results collectively paint a picture of HHR catalysis: HHR first folds to form an electronegative active site pocket to recruit a threshold occupation of cationic charges, either a Mg(2+) ion or multiple monovalent cations. Catalytically active conformations that have good in-line fitness are supported by specific metal ion coordination patterns that involve either a bridging Mg(2+) ion or multiple Na(+) ions, one of which is also in a bridging coordination pattern. In the case of a single Mg(2+) ion bound in the active site, the Mg(2+) ion undergoes a migration that is coupled with deprotonation of the nucleophile (C17:O2'). As the reaction proceeds, the Mg(2+) ion stabilizes the accumulating charge of the leaving group and significantly increases the general acid ability of G8:O2'. Further computational mutagenesis simulations suggest that the disruptions due to mutations may severely impact HHR catalysis at different stages of the reaction. Catalytic mechanisms supported by the simulation results are consistent with available structural and biochemical experiments, and together they advance our understanding of HHR catalysis.

Simulating DNA by molecular dynamics: aims, methods, and validation

The structure and dynamics of the B-DNA double helix involves subtle sequence-dependent effects which are decisive for its function, but difficult to characterize. These structural and dynamic effects can be addressed by simulations of DNA sequences in explicit solvent. Here, we present and discuss the state-of-art of B-DNA molecular dynamics simulations with the major force fields in use today. We explain why a critical analysis of the MD trajectories is required to assess their reliability, and estimate the value and limitations of these models. Overall, simulations of DNA bear great promise towards deciphering the structural and physical subtleties of this biopolymer, where much remains to be understood.

The role of salt concentration and magnesium binding in HIV-1 subtype-A and subtype-B kissing loop monomer structures

The subtype-B monomers of the human immunodeficiency virus type-1 (HIV-1) have experimentally been shown to dimerize at high salt concentration or in the presence of magnesium, while the dimerization of the subtype-A monomers requires magnesium binding at the G273 or G274 phosphate groups regardless of salt concentration. We used explicit solvent molecular dynamics (MD) simulations to investigate the conformational changes in subtype-A and -B monomers in different salt concentrations, and we found that our MD simulation results are consistent with those of experiments. At low salt concentration, hairpin loop structures of both subtypes were deformed and bases in the hairpin loop were turned inside. At high salt concentrations, the subtype-B monomer maintained the hairpin loop shape and most bases in the hairpin loop pointed out, while the subtype-A monomer showed a severe deformation. We also found that the flanking bases in the subtype-B stabilize the hairpin loop, while the flanking base G273 in the subtype-A caused a significant deformation. However, a bound magnesium ion at the G273 or G274 phosphate groups controlled the behavior of the G273 base and prevented the subtype-A monomer from deformation. We also applied restraints to both subtypes to examine the role of high salt concentration or magnesium binding. While restraints were applied, both subtypes at 0 M salt concentration maintained their shapes. However, when restraints were removed, they deformed significantly. Therefore, we suggest that the dimerization of both subtypes requires the proper conformation of the monomers which is induced by the appropriate salt strength and magnesium ion binding.

Molecular dynamics re-refinement of two different small RNA loop structures using the original NMR data suggest a common structure

Restrained molecular dynamics simulations are a robust, though perhaps underused, tool for the end-stage refinement of biomolecular structures. We demonstrate their utility-using modern simulation protocols, optimized force fields, and inclusion of explicit solvent and mobile counterions-by re-investigating the solution structures of two RNA hairpins that had previously been refined using conventional techniques. The structures, both domain 5 group II intron ribozymes from yeast ai5γ and Pylaiella littoralis, share a nearly identical primary sequence yet the published 3D structures appear quite different. Relatively long restrained MD simulations using the original NMR restraint data identified the presence of a small set of violated distance restraints in one structure and a possibly incorrect trapped bulge nucleotide conformation in the other structure. The removal of problematic distance restraints and the addition of a heating step yielded representative ensembles with very similar 3D structures and much lower pairwise RMSD values. Analysis of ion density during the restrained simulations helped to explain chemical shift perturbation data published previously. These results suggest that restrained MD simulations, with proper caution, can be used to "update" older structures or aid in the refinement of new structures that lack sufficient experimental data to produce a high quality result. Notable cautions include the need for sufficient sampling, awareness of potential force field bias (such as small angle deviations with the current AMBER force fields), and a proper balance between the various restraint weights.

The structural role of Mg2+ ions in a class I RNA polymerase ribozyme: a molecular simulation study

According to the RNA world hypothesis, self-replicating ribozymes, storing the genetic information and being able to perform catalysis, were the constituents of the first living organisms. In particular, RNA polymerase ribozymes, similar to current proteinaceous enzymatic polymerases, may have been able to promote the synthesis of RNA strands in a primitive world. Polymerase catalysis is usually assisted by Mg(2+) ions, but it is not always trivial to find out experimentally the number of Mg(2+) ions placed in the active site as well as the identity and the number of their coordination ligands. Here, we addressed this issue in an artificial class I ligase ribozyme. On the basis of a recently solved crystal structure, we constructed computational models of reactant and product states of this ribozyme, considering monometallic and bimetallic species. Our models were relaxed by force field based molecular dynamics (MD) simulations and mixed quantum-classical (QM/MM) MD. The structural and dynamical properties of these models were consistent with experimental data and were validated by a comparison with the catalytic sites of proteinaceous DNA and RNA polymerases. Consistently with enzymatic polymerases, our results suggest that class I RNA ligases most probably contain two magnesium ions in the active site and they may, therefore, catalyze the junction of two RNA strands via "a two Mg(2+) ions" mechanism.

Flexibility analysis of biomacromolecules with application to computer-aided drug design

Flexibility characteristics of biomacromolecules can be efficiently determined down to the atomic level by a graph-theoretical technique as implemented in the FIRST (Floppy Inclusion and Rigid Substructure Topology) and ProFlex software packages. The method has been successfully applied to a series of protein and nucleic acid structures. Here, we describe practical guidelines for setting up and performing a flexibility analysis, discuss current bottlenecks of the approach, and provide sample applications as to how this technique can support computer-aided drug design approaches.

Characterization of metal ion-nucleic acid interactions in solution

Metal ions are inextricably involved with nucleic acids due to their polyanionic nature. In order to understand the structure and function of RNAs and DNAs, one needs to have detailed pictures on the structural, thermodynamic, and kinetic properties of metal ion interactions with these biomacromolecules. In this review we first compile the physicochemical properties of metal ions found and used in combination with nucleic acids in solution. The main part then describes the various methods developed over the past decades to investigate metal ion binding by nucleic acids in solution. This includes for example hydrolytic and radical cleavage experiments, mutational approaches, as well as kinetic isotope effects. In addition, spectroscopic techniques like EPR, lanthanide(III) luminescence, IR and Raman as well as various NMR methods are summarized. Aside from gaining knowledge about the thermodynamic properties on the metal ion-nucleic acid interactions, especially NMR can be used to extract information on the kinetics of ligand exchange rates of the metal ions applied. The final section deals with the influence of anions, buffers, and the solvent permittivity on the binding equilibria between metal ions and nucleic acids. Little is known on some of these aspects, but it is clear that these three factors have a large influence on the interaction between metal ions and nucleic acids.

Interactions of cations with RNA loop-loop complexes

RNA loop-loop interactions are essential in many biological processes, including initiation of RNA folding into complex tertiary shapes, promotion of dimerization, and viral replication. In this article, we examine interactions of metal ions with five RNA loop-loop complexes of unique biological significance using explicit-solvent molecular-dynamics simulations. These simulations revealed the presence of solvent-accessible tunnels through the major groove of loop-loop interactions that attract and retain cations. Ion dynamics inside these loop-loop complexes were distinctly different from the dynamics of the counterion cloud surrounding RNA and depend on the number of basepairs between loops, purine sequence symmetry, and presence of unpaired nucleotides. The cationic uptake by kissing loops depends on the number of basepairs between loops. It is interesting that loop-loop complexes with similar functionality showed similarities in cation dynamics despite differences in sequence and loop size.

Understanding the origins of bacterial resistance to aminoglycosides through molecular dynamics mutational study of the ribosomal A-site

Paromomycin is an aminoglycosidic antibiotic that targets the RNA of the bacterial small ribosomal subunit. It binds in the A-site, which is one of the three tRNA binding sites, and affects translational fidelity by stabilizing two adenines (A1492 and A1493) in the flipped-out state. Experiments have shown that various mutations in the A-site result in bacterial resistance to aminoglycosides. In this study, we performed multiple molecular dynamics simulations of the mutated A-site RNA fragment in explicit solvent to analyze changes in the physicochemical features of the A-site that were introduced by substitutions of specific bases. The simulations were conducted for free RNA and in complex with paromomycin. We found that the specific mutations affect the shape and dynamics of the binding cleft as well as significantly alter its electrostatic properties. The most pronounced changes were observed in the U1406C∶U1495A mutant, where important hydrogen bonds between the RNA and paromomycin were disrupted. The present study aims to clarify the underlying physicochemical mechanisms of bacterial resistance to aminoglycosides due to target mutations.

A coarse-grained model of DNA with explicit solvation by water and ions

Solvation by water and ions has been shown to be vitally important for biological molecules, yet fully atomistic simulations of large biomolecules remain a challenge due to their high computational cost. The effect of solvation is the most pronounced in polyelectrolytes, of which DNA is a paradigmatic example. Coarse-grained (CG) representations have been developed to model the essential physics of the DNA molecule, yet almost without exception, these models replace the water and ions by implicit solvation in order to significantly reduce the computational expense. This work introduces the first coarse-grained model of DNA solvated explicitly with water and ions. To this end, we combined two established CG models; the recently developed mW-ion model [DeMille, R. C.; Molinero, V. J. Chem. Phys. 2009, 131, 034107], which reproduces the structure of aqueous ionic solutions without electrostatic interactions, was coupled to the three-sites-per-nucleotide (3SPN) CG model of DNA [Knotts, T. A., IV; et al. J. Chem. Phys. 2007, 126, 084901]. Using atomistic simulations of d(CGCGAATTCGCG)(2) as a reference, we optimized the coarse-grained interactions between DNA and solvent to reproduce the solvation structure of water and ions around CG DNA. The resulting coarse-grained model of DNA explicitly solvated by ions and water (mW/3SPN-DNA) exhibits base-pair specificity and ion-condensation effects and it is 2 orders of magnitude computationally more efficient than atomistic models. We describe the parametrization strategy and offer insight into how other CG models may be combined with a coarse-grained solvent model such as mW-ion.

Toward a consensus view of duplex RNA flexibility

The structure and flexibility of the RNA duplex has been studied using extended molecular dynamics simulations on four diverse 18-mer oligonucleotides designed to contain many copies of the 10 unique dinucleotide steps in different sequence environments. Simulations were performed using the two most popular force fields for nucleic acids simulations (AMBER and CHARMM) in their latest versions, trying to arrive to a consensus picture of the RNA flexibility. Contrary to what was found for DNA duplex (DNA(2)), no clear convergence is found for the RNA duplex (RNA(2)), but one of the force field seems to agree better with experimental data. MD simulations performed with this force field were used to fully characterize, for the first time to our knowledge, the sequence-dependent elastic properties of RNA duplexes at different levels of resolutions. The flexibility pattern of RNA(2) shows similarities with DNA(2), but also surprising differences, which help us to understand the different biological functions of both molecules. A full mesoscopic model of RNA duplex at different resolution levels is derived to be used for genome-wide description of the flexibility of double-helical fragments of RNA.

Structural dynamics of the box C/D RNA kink-turn and its complex with proteins: the role of the A-minor 0 interaction, long-residency water bridges, and structural ion-binding sites revealed by molecular simulations

Kink-turns (K-turns) are recurrent elbow-like RNA motifs that participate in protein-assisted RNA folding and contribute to RNA dynamics. We carried out a set of molecular dynamics (MD) simulations using parm99 and parmbsc0 force fields to investigate structural dynamics of the box C/D RNA and its complexes with two proteins: native archaeal L7ae protein and human 15.5 kDa protein, originally bound to very similar structure of U4 snRNA. The box C/D RNA forms K-turn with A-minor 0 tertiary interaction between its canonical (C) and noncanonical (NC) stems. The local K-turn architecture is thus different from the previously studied ribosomal K-turns 38 and 42 having A-minor I interaction. The simulations reveal visible structural dynamics of this tertiary interaction involving altogether six substates which substantially contribute to the elbow-like flexibility of the K-turn. The interaction can even temporarily shift to the A-minor I type pattern; however, this is associated with distortion of the G/A base pair in the NC-stem of the K-turn. The simulations show reduction of the K-turn flexibility upon protein binding. The protein interacts with the apex of the K-turn and with the NC-stem. The protein-RNA interface includes long-residency hydration sites. We have also found long-residency hydration sites and major ion-binding sites associated with the K-turn itself. The overall topology of the K-turn remains stable in all simulations. However, in simulations of free K-turn, we observed instability of the key C16(O2')-A7(N1) H-bond, which is a signature interaction of K-turns and which was visibly more stable in simulations of K-turns possessing A-minor I interaction. It may reflect either some imbalance of the force field or it may be a correct indication of early stages of unfolding since this K-turn requires protein binding for its stabilization. Interestingly, the 16(O2')-7(N1) H- bond is usually not fully lost since it is replaced by a water bridge with a tightly bound water, which is adenine-specific similarly as the original interaction. The 16(O2')-7(N1) H-bond is stabilized by protein binding. The stabilizing effect is more visible with the human 15.5 kDa protein, which is attributed to valine to arginine substitution in the binding site. The behavior of the A-minor interaction is force-field-dependent because the parmbsc0 force field attenuates the A-minor fluctuations compared to parm99 simulations. Behavior of other regions of the box C/D RNA is not sensitive to the force field choice. Simulation with net-neutralizing Na(+) and 0.2 M excess salt conditions appear in all aspects equivalent. The simulations show loss of a hairpin tetraloop, which is not part of the K-turn. This was attributed to force field limitations.

Molecular recognition of RNA: challenges for modelling interactions and plasticity

There is growing interest in molecular recognition processes of RNA because of RNA's widespread involvement in biological processes. Computational approaches are increasingly used for analysing and predicting binding to RNA, fuelled by encouraging progress in developing simulation, free energy and docking methods for nucleic acids. These developments take into account challenges regarding the energetics of RNA-ligand binding, RNA plasticity, and the presence of water molecules and ions in the binding interface. Accordingly, we will detail advances in force field and scoring function development for molecular dynamics (MD) simulations, free energy computations and docking calculations of nucleic acid complexes. Furthermore, we present methods that can detect moving parts within RNA structures based on graph-theoretical approaches or normal mode analysis (NMA). As an example of the successful use of these developments, we will discuss recent structure-based drug design approaches that focus on the bacterial ribosomal A-site RNA as a drug target.

Binding of the bacteriophage P22 N-peptide to the boxB RNA motif studied by molecular dynamics simulations

Protein-RNA interactions are important for many cellular processes. The Nut-utilization site (N)-protein of bacteriophages contains an N-terminal arginine-rich motif that undergoes a folding transition upon binding to the boxB RNA hairpin loop target structure. Molecular dynamics simulations were used to investigate the dynamics of the P22 N-peptide-boxB complex and to elucidate the energetic contributions to binding. In addition, the free-energy changes of RNA and peptide conformational adaptation to the bound forms, as well as the role of strongly bound water molecules at the peptide-RNA interface, were studied. The influence of peptide amino acid substitutions and the salt dependence of interaction were investigated and showed good agreement with available experimental results. Several tightly bound water molecules were found at the RNA-binding interface in both the presence and absence of N-peptide. Explicit consideration of the waters resulted in shifts of calculated contributions during the energetic analysis, but overall similar binding energy contributions were found. Of interest, it was found that the electrostatic field of the RNA has a favorable influence on the coil-to-alpha-helix transition of the N-peptide already outside of the peptide-binding site. This result may have important implications for understanding peptide-RNA complex formation, which often involves coupled folding and association processes. It indicates that electrostatic interactions near RNA molecules can lead to a shift in the equilibrium toward the bound form of an interacting partner before it enters the binding pocket.

Elbow flexibility of the kt38 RNA kink-turn motif investigated by free-energy molecular dynamics simulations

Kink-turns (K-turns) are common structural motifs that can introduce sharp kinks into double-stranded RNA, and have been proposed to mediate large-scale motions in the ribosome. K-turns consist of a bulge loop region flanked by trans sugar-Hoogsteen G:A pairs, and the sharp kink conformation is stabilized by A-minor interactions (adenine contacting a G:C basepair in the minor groove). Umbrella-sampling molecular dynamics simulations were used to disrupt an A-minor interaction in the ribosomal kt38 turn and to calculate the associated free-energy change. Coupling of umbrella sampling with replica exchanges between neighboring umbrella-sampling intervals could further improve the convergence of the free-energy calculations. The simulations revealed a coupled A-minor disruption and global opening of the K-turn motif, and allowed us to characterize several intermediate A-minor conformations. The calculated free-energy profile indicated a meta-stable, semi-open structure of slightly higher free energy ( approximately 1 kcal mol(-1)), separated by a small free-energy barrier ( approximately 1.5 kcal mol(-1)) from the closed (highly kinked) form. Both K-turn states are stabilized by distinct variants of the A-minor interaction. Further opening of the K-turn structure required significantly larger free-energy changes. The semi-open form had a reduced kink angle compatible with experimental data on K-turn solution structures, and opening was coupled to a continuous global unwinding of the K-turn motif. The range of free-energy changes associated with kt38 opening and unwinding are compatible with the idea that K-turns may facilitate biologically relevant motions during large-scale ribosome dynamics.