Molecular dynamics simulation of the human U2B" protein complex with U2 snRNA hairpin IV in aqueous solution

Role of salt-bridging interactions in recognition of viral RNA by arginine-rich peptides

Interactions between RNA molecules and proteins are critical to many cellular processes and are implicated in various diseases. The RNA-peptide complexes are good model systems to probe the recognition mechanism of RNA by proteins. In this work, we report studies on the binding-unbinding process of a helical peptide from a viral RNA element using nonequilibrium molecular dynamics simulations. We explored the existence of various dissociation pathways with distinct free-energy profiles that reveal metastable states and distinct barriers to peptide dissociation. We also report the free-energy differences for each of the four pathways to be 96.47 ± 12.63, 96.1 ± 10.95, 91.83 ± 9.81, and 92 ± 11.32 kcal/mol. Based on the free-energy analysis, we further propose the preferred pathway and the mechanism of peptide dissociation. The preferred pathway is characterized by the formation of sequential hydrogen-bonding and salt-bridging interactions between several key arginine amino acids and the viral RNA nucleotides. Specifically, we identified one arginine amino acid (R8) of the peptide to play a significant role in the recognition mechanism of the peptide by the viral RNA molecule.

Atomistic Analysis of ToxN and ToxI Complex Unbinding Mechanism

ToxIN is a triangular structure formed by three protein toxins (ToxNs) and three specific noncoding RNA antitoxins (ToxIs). To respond to stimuli, ToxI is preferentially degraded, releasing the ToxN. Thus, the dynamic character is essential in the normal function interactions between ToxN and ToxI. Here, equilibrated molecular dynamics (MD) simulations were performed to study the stability of ToxN and ToxI. The results indicate that ToxI adjusts the conformation of 3' and 5' termini to bind to ToxN. Steered molecular dynamics (SMD) simulations combined with the recently developed thermodynamic integration in 3nD (TI3nD) method were carried out to investigate ToxN unbinding from the ToxIN complex. The potentials of mean force (PMFs) and atomistic pictures suggest the unbinding mechanism as follows: (1) dissociation of the 5' terminus from ToxN, (2) missing the interactions involved in the 3' terminus of ToxI without three nucleotides (G31, A32, and A33), (3) starting to unfold for ToxI, (4) leaving the binding package of ToxN for three nucleotides of ToxI, (5) unfolding of ToxI. This work provides information on the structure-function relationship at the atomistic level, which is helpful for designing new potent antibacterial drugs in the future.

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.

Synergy between NMR measurements and MD simulations of protein/RNA complexes: application to the RRMs, the most common RNA recognition motifs

RNA recognition motif (RRM) proteins represent an abundant class of proteins playing key roles in RNA biology. We present a joint atomistic molecular dynamics (MD) and experimental study of two RRM-containing proteins bound with their single-stranded target RNAs, namely the Fox-1 and SRSF1 complexes. The simulations are used in conjunction with NMR spectroscopy to interpret and expand the available structural data. We accumulate more than 50 μs of simulations and show that the MD method is robust enough to reliably describe the structural dynamics of the RRM-RNA complexes. The simulations predict unanticipated specific participation of Arg142 at the protein-RNA interface of the SRFS1 complex, which is subsequently confirmed by NMR and ITC measurements. Several segments of the protein-RNA interface may involve competition between dynamical local substates rather than firmly formed interactions, which is indirectly consistent with the primary NMR data. We demonstrate that the simulations can be used to interpret the NMR atomistic models and can provide qualified predictions. Finally, we propose a protocol for 'MD-adapted structure ensemble' as a way to integrate the simulation predictions and expand upon the deposited NMR structures. Unbiased μs-scale atomistic MD could become a technique routinely complementing the NMR measurements of protein-RNA complexes.

Molecular Dynamics: New Frontier in Personalized Medicine

The field of drug discovery has witnessed infinite development over the last decade with the demand for discovery of novel efficient lead compounds. Although the development of novel compounds in this field has seen large failure, a breakthrough in this area might be the establishment of personalized medicine. The trend of personalized medicine has shown stupendous growth being a hot topic after the successful completion of Human Genome Project and 1000 genomes pilot project. Genomic variant such as SNPs play a vital role with respect to inter individual's disease susceptibility and drug response. Hence, identification of such genetic variants has to be performed before administration of a drug. This process requires high-end techniques to understand the complexity of the molecules which might bring an insight to understand the compounds at their molecular level. To sustenance this, field of bioinformatics plays a crucial role in revealing the molecular mechanism of the mutation and thereby designing a drug for an individual in fast and affordable manner. High-end computational methods, such as molecular dynamics (MD) simulation has proved to be a constitutive approach to detecting the minor changes associated with an SNP for better understanding of the structural and functional relationship. The parameters used in molecular dynamic simulation elucidate different properties of a macromolecule, such as protein stability and flexibility. MD along with docking analysis can reveal the synergetic effect of an SNP in protein-ligand interaction and provides a foundation for designing a particular drug molecule for an individual. This compelling application of computational power and the advent of other technologies have paved a promising way toward personalized medicine. In this in-depth review, we tried to highlight the different wings of MD toward personalized medicine.

Hot Spot Occlusion from Bulk Water: a Comprehensive Study of the Complex between the Lysozyme HEL and the Antibody FVD1.3

Alanine scanning of protein-protein interfaces has shown that there are some residues in the protein-protein interfaces, responsible for most of the binding free energy, which are called hot spots. Hot spots tend to exist in densely packed central clusters, and a hypothesis has been proposed that considers that inaccessibility to the solvent must be a necessary condition to define a residue as a binding hot spot. This O-ring hypothesis is mainly based on the analysis of the accessible surface area (ASA) of 23 static, crystallographic structures of protein complexes. It is known, however, that protein flexibility allows for temporary exposures of buried interfacial groups, and even though the ASA provides a general trend of the propensity for hydration, protein/solvent-specific interactions or hydrogen bonding cannot be considered here. Therefore, a microscopic level, atomistic picture of hot spot solvation is needed to support the O-ring hypothesis. In this study, we began by applying a computational alanine-scanning mutagenesis technique, which reproduces the experimental results and allows for decomposing the binding free energy difference in its different energetic factors. Subsequently, we calculated the radial distribution function and residence times of the water molecules near the hot/warm spots to study the importance of the water environment around those energetically important amino acid residues. This study shows that within a flexible, dynamic protein framework, the warm/hot spot residues are, indeed, kept sheltered from the bulk solvent during the whole simulation, which allows a better interacting microenvironment.

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .

Cations and hydration in catalytic RNA: molecular dynamics of the hepatitis delta virus ribozyme

The hepatitis delta virus (HDV) ribozyme is an RNA enzyme from the human pathogenic HDV. Cations play a crucial role in self-cleavage of the HDV ribozyme, by promoting both folding and chemistry. Experimental studies have revealed limited but intriguing details on the location and structural and catalytic functions of metal ions. Here, we analyze a total of approximately 200 ns of explicit-solvent molecular dynamics simulations to provide a complementary atomistic view of the binding of monovalent and divalent cations as well as water molecules to reaction precursor and product forms of the HDV ribozyme. Our simulations find that an Mg2+ cation binds stably, by both inner- and outer-sphere contacts, to the electronegative catalytic pocket of the reaction precursor, in a position to potentially support chemistry. In contrast, protonation of the catalytically involved C75 in the precursor or artificial placement of this Mg2+ into the product structure result in its swift expulsion from the active site. These findings are consistent with a concerted reaction mechanism in which C75 and hydrated Mg2+ act as general base and acid, respectively. Monovalent cations bind to the active site and elsewhere assisted by structurally bridging long-residency water molecules, but are generally delocalized.

Loss of G-A base pairs is insufficient for achieving a large opening of U4 snRNA K-turn motif

Upon binding to the 15.5K protein, two tandem-sheared G–A base pairs are formed in the internal loop of the kink-turn motif of U4 snRNA (Kt-U4). We have reported that the folding of Kt-U4 is assisted by protein binding. Unstable interactions that contribute to a large opening of the free RNA (‘k–e motion’) were identified using locally enhanced sampling molecular dynamics simulations, results that agree with experiments. A detailed analysis of the simulations reveals that the k–e motion in Kt-U4 is triggered both by loss of G–A base pairs in the internal loop and backbone flexibility in the stems. Essential dynamics show that the loss of G–A base pairs is correlated along the first mode but anti-correlated along the third mode with the k–e motion. Moreover, when enhanced sampling was confined to the internal loop, the RNA adopted an alternative conformation characterized by a sharper kink, opening of G–A base pairs and modified stacking interactions. Thus, loss of G–A base pairs is insufficient for achieving a large opening of the free RNA. These findings, supported by previously published RNA structure probing experiments, suggest that G–A base pair formation occurs upon protein binding, thereby stabilizing a selective orientation of the stems.

Hinge-like motions in RNA kink-turns: the role of the second a-minor motif and nominally unpaired bases

Kink-turn (K-turn) motifs are asymmetric internal loops found at conserved positions in diverse RNAs, with sharp bends in phosphodiester backbones producing V-shaped structures. Explicit-solvent molecular dynamics simulations were carried out for three K-turns from 23S rRNA, i.e., Kt-38 located at the base of the A-site finger, Kt-42 located at the base of the L7/L12 stalk, and Kt-58 located in domain III, and for the K-turn of human U4 snRNA. The simulations reveal hinge-like K-turn motions on the nanosecond timescale. The first conserved A-minor interaction between the K-turn stems is entirely stable in all simulations. The angle between the helical arms of Kt-38 and Kt-42 is regulated by local variations of the second A-minor (type I) interaction between the stems. Its variability ranges from closed geometries to open ones stabilized by insertion of long-residency waters between adenine and cytosine. The simulated A-minor geometries fully agree with x-ray data. Kt-58 and Kt-U4 exhibit similar elbow-like motions caused by conformational change of the adenosine from the nominally unpaired region. Despite the observed substantial dynamics of K-turns, key tertiary interactions are stable and no sign of unfolding is seen. We suggest that some K-turns are flexible elements mediating large-scale ribosomal motions during the protein synthesis cycle.

Long-residency hydration, cation binding, and dynamics of loop E/helix IV rRNA-L25 protein complex

Molecular dynamics simulations of RNA-protein complex between Escherichia coli loop E/helix IV (LE/HeIV) rRNA and L25 protein reveal a qualitative agreement between the experimental and simulated structures. The major groove of LE is a prominent rRNA cation-binding site. Divalent cations rigidify the LE major groove geometry whereas in the absence of divalent cations LE extensively interacts with monovalent cations via inner-shell binding. The HeIV region shows bistability of its major groove explaining the observed differences between x-ray and NMR structures. In agreement with the experiments, the simulations suggest that helix-alpha1 of L25 is the least stable part of the protein. Inclusion of Mg2+ cations into the simulations causes perturbation of basepairing at the LE/HeIV junction, which does not, however, affect the protein binding. The rRNA-protein complex is mediated by a number of highly specific hydration sites with long-residing water molecules and two of them are bound throughout the entire 24-ns simulation. Long-residing water molecules are seen also outside the RNA-protein contact areas with water-binding times substantially enhanced compared to simulations of free RNA. Long-residency hydration sites thus represent important elements of the three-dimensional structure of rRNA.

The structure and dynamics of the U4/U6 snRNP: implications for pre-mRNA splicing and use as a model system to investigate the RNA-mediated effects of (5F)Ura

Pre-mRNA splicing is one of the most complex and intricate processes in eukaryotic cell biology. Over the past decade, my laboratory has been interested in determining the structures of RNA components of the spliceosome, and in investigating how the structure, stability and dynamics of these RNA components are perturbed by nucleoside analog substitution. In particular, we have investigated the U4/U6 snRNA complex as a model system for understanding the biophysical basis for the RNA-mediated effects of the widely-used anticancer drug 5-fluorouracil ((5F)Ura). In this review, our studies that have provided novel information concerning the structure of U4 snRNA and its interactions with U6 snRNA and the Sm (or common) snRNA binding proteins are summarized. These studies have also quantified the structural and thermodynamic consequences of (5F)Ura in this model system. Our work to date provides the foundation on which future studies investigating the biophysical basis for spliceosomal assembly and for clarifying the mechanisms of anticancer drugs targeted at nucleic acid-mediated processes will be developed.