Interactions of cations with RNA loop-loop complexes

Structure folding of RNA kissing complexes in salt solutions: predicting 3D structure, stability, and folding pathway

RNA kissing complexes are essential for genomic RNA dimerization and regulation of gene expression, and their structures and stability are critical to their biological functions. In this work, we used our previously developed coarse-grained model with an implicit structure-based electrostatic potential to predict three-dimensional (3D) structures and stability of RNA kissing complexes in salt solutions. For extensive RNA kissing complexes, our model shows great reliability in predicting 3D structures from their sequences, and our additional predictions indicate that the model can capture the dependence of 3D structures of RNA kissing complexes on monovalent/divalent ion concentrations. Moreover, the comparisons with extensive experimental data show that the model can make reliable predictions on the stability for various RNA kissing complexes over wide ranges of monovalent/divalent ion concentrations. Notably, for RNA kissing complexes, our further analyses show the important contribution of coaxial stacking to the 3D structures and stronger stability than the corresponding kissing-interface duplexes at high salts. Furthermore, our comprehensive analyses for RNA kissing complexes reveal that the thermally folding pathway for a complex sequence is mainly determined by the relative stability of two possible folded states of kissing complex and extended duplex, which can be significantly modulated by its sequence.

RNA as a Complex Polymer with Coupled Dynamics of Ions and Water in the Outer Solvation Sphere

We unravel the internal and collective modes of a widely studied 58-nucleotide rRNA fragment in solvent using atomically detailed molecular dynamics simulations. The variation of lifetimes for water hydrogen bonds with nucleotide groups indicates heterogeneity of water dynamics on the RNA surface. The time scales of interactions of the discrete water molecules with RNA nucleotides extend from several hundred picoseconds to a few nanoseconds. We determine all of the association sites and the spatial distribution of residence times for Mg²⁺, K⁺, and water molecules in those sites. We provide insights into the population of Mg²⁺ and K⁺ ions and water molecules in the outer sphere and how their fluctuations are intricately linked with the kinetics of the 58-mer. We find that many of the long-lived Mg²⁺ sites identified from the simulations agree with the locations of ions in the X-ray structure. We determine the excess ion atmosphere around the rRNA and compare it with experimental data. We investigate the collective behavior of RNA, ions, and water, by performing a joint principle component analysis for the Cartesian coordinates of the RNA phosphorus atoms and for the occupation counts of the association sites. Our results indicate that the 58-mer system is a complex polymer, composed of RNA that is encased by a fluctuating network of associated counterions, co-ions, and water.

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.

Electrostatics Explains the Position-Dependent Effect of G⋅U Wobble Base Pairs on the Affinity of RNA Kissing Complexes

In the RNA realm, non-Watson-Crick base pairs are abundant and can affect both the RNA 3D structure and its function. Here, we investigated the formation of RNA kissing complexes in which the loop-loop interaction is modulated by non-Watson-Crick pairs. Mass spectrometry, surface plasmon resonance, and UV-melting experiments show that the G⋅U wobble base pair favors kissing complex formation only when placed at specific positions. We tried to rationalize this effect by molecular modeling, including molecular mechanics Poisson-Boltzmann surface area (MMPBSA) thermodynamics calculations and PBSA calculations of the electrostatic potential surfaces. Modeling reveals that the G⋅U stabilization is due to a specific electrostatic environment defined by the base pairs of the entire loop-loop region. The loop is not symmetric, and therefore the identity and position of each base pair matters. Predicting and visualizing the electrostatic environment created by a given sequence can help to design specific kissing complexes with high affinity, for potential therapeutic, nanotechnology or analytical applications.

Extended molecular dynamics of a c-kit promoter quadruplex

The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.

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.

Bioinformatics and molecular dynamics simulation study of L1 stalk non-canonical rRNA elements: kink-turns, loops, and tetraloops

The L1 stalk is a prominent mobile element of the large ribosomal subunit. We explore the structure and dynamics of its non-canonical rRNA elements, which include two kink-turns, an internal loop, and a tetraloop. We use bioinformatics to identify the L1 stalk RNA conservation patterns and carry out over 11.5 μs of MD simulations for a set of systems ranging from isolated RNA building blocks up to complexes of L1 stalk rRNA with the L1 protein and tRNA fragment. We show that the L1 stalk tetraloop has an unusual GNNA or UNNG conservation pattern deviating from major GNRA and YNMG RNA tetraloop families. We suggest that this deviation is related to a highly conserved tertiary contact within the L1 stalk. The available X-ray structures contain only UCCG tetraloops which in addition differ in orientation (anti vs syn) of the guanine. Our analysis suggests that the anti orientation might be a mis-refinement, although even the anti interaction would be compatible with the sequence pattern and observed tertiary interaction. Alternatively, the anti conformation may be a real substate whose population could be pH-dependent, since the guanine syn orientation requires protonation of cytosine in the tertiary contact. In absence of structural data, we use molecular modeling to explore the GCCA tetraloop that is dominant in bacteria and suggest that the GCCA tetraloop is structurally similar to the YNMG tetraloop. Kink-turn Kt-77 is unusual due to its 11-nucleotide bulge. The simulations indicate that the long bulge is a stalk-specific eight-nucleotide insertion into consensual kink-turn only subtly modifying its structural dynamics. We discuss a possible evolutionary role of helix H78 and a mechanism of L1 stalk interaction with tRNA. We also assess the simulation methodology. The simulations provide a good description of the studied systems with the latest bsc0χOL3 force field showing improved performance. Still, even bsc0χOL3 is unable to fully stabilize an essential sugar-edge H-bond between the bulge and non-canonical stem of the kink-turn. Inclusion of Mg(2+) ions may deteriorate the simulations. On the other hand, monovalent ions can in simulations readily occupy experimental Mg(2+) binding sites.

Effect of initial ion positions on the interactions of monovalent and divalent ions with a DNA duplex as revealed with atomistic molecular dynamics simulations

Monovalent (Na(+)) and divalent (Mg(2+)) ion distributions around the Dickerson-Drew dodecamer were studied by atomistic molecular dynamics (MD) simulations with AMBER molecular modeling software. Different initial placements of ions were tried and the resulting effects on the ion distributions around DNA were investigated. For monovalent ions, results were found to be nearly independent of initial cation coordinates. However, Mg(2+) ions demonstrated a strong initial coordinate dependent behavior. While some divalent ions initially placed near the DNA formed essentially permanent direct coordination complexes with electronegative DNA atoms, Mg(2+) ions initially placed further away from the duplex formed a full, nonexchanging, octahedral first solvation shell. These fully solvated cations were still capable of binding with DNA with events lasting up to 20 ns, and in comparison were bound much longer than Na(+) ions. Force field parameters were also investigated with modest and little differences arising from ion (ions94 and ions08) and nucleic acid description (ff99, ff99bsc0, and ff10), respectively. Based on known Mg(2+) ion solvation structure, we conclude that in most cases Mg(2+) ions retain their first solvation shell, making only solvent-mediated contacts with DNA duplex. The proper way to simulate Mg(2+) ions around DNA duplex, therefore, should begin with ions placed in the bulk water.

Weakly charged cationic nanoparticles induce DNA bending and strand separation

Weakly charged cationic nanoparticles cause structural changes including local denaturing and compaction to DNA under mild conditions. The charged ligands bind to the phosphate backbone of DNA and the uncharged ligands penetrate the helix and disrupt base pairing. Mobility shifts in electrophoresis, molecular dynamics, and UV-vis spectrophotometry give clues to the details of the interactions.

Magnesium Ion-Water Coordination and Exchange in Biomolecular Simulations

Magnesium ions have an important role in the structure and folding mechanism of ribonucleic acid systems. To properly simulate these biophysical processes, the applied molecular models should reproduce, among other things, the kinetic properties of the ions in water solution. Here, we have studied the kinetics of the binding of magnesium ions with water molecules and nucleic acid systems using molecular dynamics simulation in detail. We have validated the parameters used in biomolecular force fields, such as AMBER and CHARMM, for Mg(2+) ions and also for the biologically relevant ions Na(+), K(+), and Ca(2+) together with three different water models (TIP3P, SPC/E, and TIP5P). The results show that Mg(2+) ions have a slower exchange rate than Na(+), K(+), and Ca(2+) in agreement with the experimental trend, but the simulated value underestimates the experimentally observed Mg(2+)-water exchange rate by several orders of magnitude, irrespective of the force field and water model. A new set of parameters for Mg(2+) was developed to reproduce the experimental kinetic data. This set also leads to better reproduction of structural data than existing models. We have applied the new parameter set to Mg(2+) binding with a monophosphate model system and with the purine riboswitch, add A-riboswitch. In line with the Mg(2+)-water results, the newly developed parameters show a better description of the structure and kinetics of the Mg(2+)-phosphate binding than all other models. The characterization of the ion binding to the riboswitch system shows that the new parameter set does not affect the global structure of the ribonucleic acid system or the number of ions involved in direct or indirect binding. A slight decrease in the number of water-bridged contacts between A-riboswitch and the Mg(2+) ion is observed. The results support the ability of the newly developed parameters to improve the kinetic description of the Mg(2+) and phosphate ions and their applicability in nucleic acid simulation.

A-minor tertiary interactions in RNA kink-turns. Molecular dynamics and quantum chemical analysis

The RNA kink-turn is an important recurrent RNA motif, an internal loop with characteristic consensus sequence forming highly conserved three-dimensional structure. Functional arrangement of RNA kink-turns shows a sharp bend in the phosphodiester backbone. Among other signature interactions, kink-turns form A-minor interaction between their two stems. Most kink-turns possess extended A-minor I (A-I) interaction where adenine of the second A•G base pair of the NC-stem interacts with the first canonical pair of the C-stem (i.e., the receptor pair) via trans-sugar-edge/sugar-edge (tSS) and cis-sugar-edge/sugar-edge (cSS) interactions. The remaining kink-turns have less compact A-minor 0 (A-0) interaction with just one tSS contact. We show that kink-turns with A-I in ribosomal X-ray structures keep G═C receptor base pair during evolution while the inverted pair (C═G) is not realized. In contrast, kink-turns with A-0 in the observed structures alternate G═C and C═G base pairs in sequences. We carried out an extended set (~5 μs) of explicit-solvent molecular dynamics simulations of kink-turns to rationalize this structural/evolutionary pattern. The simulations were done using a net-neutral Na(+) cation atmosphere (with ~0.25 M cation concentration) supplemented by simulations with either excess salt KCl atmosphere or inclusion of Mg(2+). The results do not seem to depend on the treatment of ions. The simulations started with X-ray structures of several kink-turns while we tested the response of the simulated system to base substitutions, modest structural perturbations and constraints. The trends seen in the simulations reveal that the A-I/G═C arrangement is preferred over all three other structures. The A-I/C═G triple appears structurally entirely unstable, consistent with the covariation patterns seen during the evolution. The A-0 arrangements tend to shift toward the A-I pattern in simulations, which suggests that formation of the A-0 interaction is likely supported by the surrounding protein and RNA molecules. A-0 may also be stabilized by additional kink-turn nucleotides not belonging to the kink-turn consensus, as shown for the kink-turn from ribosomal Helix 15. Quantum-chemical calculations on all four A-minor triples suggest that there is a different balance of electrostatic and dispersion stabilization in the A-I/G═C and A-I/C═G triples, which may explain different behavior of these otherwise isosteric triples in the context of kink-turns.