RNA recognition motifs (RRMs) are a key class of proteins that primarily bind single-stranded RNAs. In this study, we applied standard atomistic molecular dynamics simulations to obtain insights into the intricate binding dynamics between uridine-rich RNAs and TbRGG2 RRM using the recently developed OL3-Stafix AMBER force field, which improves the description of single-stranded RNA molecules. Complementing structural experiments that unveil a primary binding mode with a single uridine bound, our simulations uncover two supplementary binding modes in which adjacent nucleotides encroach upon the binding pocket. This leads to a unique molecular mechanism through which the TbRGG2 RRM is capable of rapidly transitioning the U-rich sequence. In contrast, the presence of non-native cytidines induces stalling and destabilization of the complex. By leveraging extensive equilibrium dynamics and a large variety of binding states, TbRGG2 RRM effectively expedites diffusion along the RNA substrate while ensuring robust selectivity for U-rich sequences despite featuring a solitary binding pocket. We further substantiate our description of the complex dynamics by simulating the fully spontaneous association process of U-rich sequences to the TbRGG2 RRM. Our study highlights the critical role of dynamics and auxiliary binding states in interface dynamics employed by RNA-binding proteins, which is not readily apparent in traditional structural studies but could represent a general type of binding strategy employed by many RNA-binding proteins.

• MeSH
• konformace nukleové kyseliny MeSH
• motiv rozpoznávající RNA MeSH
• proteiny vázající RNA metabolismus   chemie   MeSH
• RNA * chemie   metabolismus   MeSH
• simulace molekulární dynamiky * MeSH
• vazba proteinů MeSH
• vazebná místa MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• proteiny vázající RNA MeSH
• RNA * MeSH

Recognition of single-stranded RNA (ssRNA) by RNA recognition motif (RRM) domains is an important class of protein-RNA interactions. Many such complexes were characterized using nuclear magnetic resonance (NMR) and/or X-ray crystallography techniques, revealing ensemble-averaged pictures of the bound states. However, it is becoming widely accepted that better understanding of protein-RNA interactions would be obtained from ensemble descriptions. Indeed, earlier molecular dynamics simulations of bound states indicated visible dynamics at the RNA-RRM interfaces. Here, we report the first atomistic simulation study of spontaneous binding of short RNA sequences to RRM domains of HuR and SRSF1 proteins. Using a millisecond-scale aggregate ensemble of unbiased simulations, we were able to observe a few dozen binding events. HuR RRM3 utilizes a pre-binding state to navigate the RNA sequence to its partially disordered bound state and then to dynamically scan its different binding registers. SRSF1 RRM2 binding is more straightforward but still multiple-pathway. The present study necessitated development of a goal-specific force field modification, scaling down the intramolecular van der Waals interactions of the RNA which also improves description of the RNA-RRM bound state. Our study opens up a new avenue for large-scale atomistic investigations of binding landscapes of protein-RNA complexes, and future perspectives of such research are discussed.

• MeSH
• HuR protein metabolismus   MeSH
• motiv rozpoznávající RNA genetika   MeSH
• proteiny vázající RNA * metabolismus   MeSH
• RNA * chemie   MeSH
• RRM proteiny metabolismus   MeSH
• vazba proteinů MeSH
• vazebná místa MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• HuR protein MeSH
• proteiny vázající RNA * MeSH
• RNA * MeSH
• RRM proteiny MeSH

The conserved protein Hfq is a key factor in the RNA-mediated control of gene expression in most known bacteria. The transient intermediates Hfq forms with RNA support intricate and robust regulatory networks. In Pseudomonas, Hfq recognizes repeats of adenine-purine-any nucleotide (ARN) in target mRNAs via its distal binding side, and together with the catabolite repression control (Crc) protein, assembles into a translation-repression complex. Earlier experiments yielded static, ensemble-averaged structures of the complex, but details of its interface dynamics and assembly pathway remained elusive. Using explicit solvent atomistic molecular dynamics simulations, we modeled the extensive dynamics of the Hfq-RNA interface and found implications for the assembly of the complex. We predict that syn/anti flips of the adenine nucleotides in each ARN repeat contribute to a dynamic recognition mechanism between the Hfq distal side and mRNA targets. We identify a previously unknown binding pocket that can accept any nucleotide and propose that it may serve as a 'status quo' staging point, providing nonspecific binding affinity, until Crc engages the Hfq-RNA binary complex. The dynamical components of the Hfq-RNA recognition can speed up screening of the pool of the surrounding RNAs, participate in rapid accommodation of the RNA on the protein surface, and facilitate competition among different RNAs. The register of Crc in the ternary assembly could be defined by the recognition of a guanine-specific base-phosphate interaction between the first and last ARN repeats of the bound RNA. This dynamic substrate recognition provides structural rationale for the stepwise assembly of multicomponent ribonucleoprotein complexes nucleated by Hfq-RNA binding.

• Klíčová slova
• ARN repeats, Crc protein, Hfq protein, RNA metabolism, RNA-binding protein, dynamic recognition, molecular dynamics, protein–nucleic acid interaction,
• MeSH
• bakteriální RNA chemie   genetika   metabolismus   MeSH
• konformace nukleové kyseliny MeSH
• konformace proteinů MeSH
• nukleotidové motivy * MeSH
• protein hostitelského faktoru 1 chemie   genetika   metabolismus   MeSH
• Pseudomonas aeruginosa genetika   metabolismus   MeSH
• regulace genové exprese u bakterií * MeSH
• vazba proteinů MeSH
• vazebná místa MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• bakteriální RNA MeSH
• protein hostitelského faktoru 1 MeSH

Holliday junctions (HJs) are four-way DNA structures that occur in DNA repair by homologous recombination. Specialized nucleases, termed resolvases, remove (i.e., resolve) HJs. The bacterial protein RuvC is a canonical resolvase that introduces two symmetric cuts into the HJ. For complete resolution of the HJ, the two cuts need to be tightly coordinated. They are also specific for cognate DNA sequences. Using a combination of structural biology, biochemistry, and a computational approach, here we show that correct positioning of the substrate for cleavage requires conformational changes within the bound DNA. These changes involve rare high-energy states with protein-assisted base flipping that are readily accessible for the cognate DNA sequence but not for non-cognate sequences. These conformational changes and the relief of protein-induced structural tension of the DNA facilitate coordination between the two cuts. The unique DNA cleavage mechanism of RuvC demonstrates the importance of high-energy conformational states in nucleic acid readouts.

• MeSH
• arginin chemie   MeSH
• bakteriální proteiny chemie   MeSH
• biokatalýza MeSH
• DNA bakterií chemie   metabolismus   MeSH
• křížová struktura DNA chemie   MeSH
• párování bází MeSH
• sekvence nukleotidů MeSH
• simulace molekulární dynamiky MeSH
• Thermus thermophilus metabolismus   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• arginin MeSH
• bakteriální proteiny MeSH
• DNA bakterií MeSH
• křížová struktura DNA MeSH

Molecular dynamics (MD) simulations became a leading tool for investigation of structural dynamics of nucleic acids. Despite recent efforts to improve the empirical potentials (force fields, ffs), RNA ffs have persisting deficiencies, which hamper their utilization in quantitatively accurate simulations. Previous studies have shown that at least two salient problems contribute to difficulties in the description of free-energy landscapes of small RNA motifs: (i) excessive stabilization of the unfolded single-stranded RNA ensemble by intramolecular base-phosphate and sugar-phosphate interactions and (ii) destabilization of the native folded state by underestimation of stability of base pairing. Here, we introduce a general ff term (gHBfix) that can selectively fine-tune nonbonding interaction terms in RNA ffs, in particular, the H bonds. The gHBfix potential affects the pairwise interactions between all possible pairs of the specific atom types, while all other interactions remain intact; i.e., it is not a structure-based model. In order to probe the ability of the gHBfix potential to refine the ff nonbonded terms, we performed an extensive set of folding simulations of RNA tetranucleotides and tetraloops. On the basis of these data, we propose particular gHBfix parameters to modify the AMBER RNA ff. The suggested parametrization significantly improves the agreement between experimental data and the simulation conformational ensembles, although our current ff version still remains far from being flawless. While attempts to tune the RNA ffs by conventional reparametrizations of dihedral potentials or nonbonded terms can lead to major undesired side effects, as we demonstrate for some recently published ffs, gHBfix has a clear promising potential to improve the ff performance while avoiding introduction of major new imbalances.

• MeSH
• RNA chemie   MeSH
• simulace molekulární dynamiky * MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• RNA MeSH

The RNA recognition motif (RRM) is the most common RNA binding domain across eukaryotic proteins. It is therefore of great value to engineer its specificity to target RNAs of arbitrary sequence. This was recently achieved for the RRM in Rbfox protein, where four mutations R118D, E147R, N151S, and E152T were designed to target the precursor to the oncogenic miRNA 21. Here, we used a variety of molecular dynamics-based approaches to predict specific interactions at the binding interface. Overall, we have run approximately 50 microseconds of enhanced sampling and plain molecular dynamics simulations on the engineered complex as well as on the wild-type Rbfox·pre-miRNA 20b from which the mutated systems were designed. Comparison with the available NMR data on the wild type molecules (protein, RNA, and their complex) served to establish the accuracy of the calculations. Free energy calculations suggest that further improvements in affinity and selectivity are achieved by the S151T replacement.

• MeSH
• konformace nukleové kyseliny MeSH
• lidé MeSH
• mikro RNA chemie   genetika   metabolismus   MeSH
• molekulární modely MeSH
• motiv rozpoznávající RNA * genetika   MeSH
• nukleární magnetická rezonance biomolekulární MeSH
• proteinové inženýrství MeSH
• proteiny vázající RNA chemie   genetika   metabolismus   MeSH
• RNA chemie   metabolismus   MeSH
• sekvence aminokyselin MeSH
• simulace molekulární dynamiky MeSH
• stabilita RNA MeSH
• vazba proteinů MeSH
• vazebná místa genetika   MeSH
• výpočetní biologie MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• mikro RNA MeSH
• MIRN20b microRNA, human MeSH Prohlížeč
• MIRN21 microRNA, human MeSH Prohlížeč
• proteiny vázající RNA MeSH
• RNA MeSH

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.

• MeSH
• DNA chemie   MeSH
• katalýza MeSH
• konformace nukleové kyseliny * MeSH
• počítačová simulace MeSH
• RNA chemie   MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• přehledy MeSH
• Research Support, N.I.H., Extramural MeSH
• Názvy látek
• DNA MeSH
• RNA MeSH

The cyclooxygenase-2 is a pro-inflammatory and cancer marker, whose mRNA stability and translation is regulated by the CUG-binding protein 2 interacting with AU-rich sequences in the 3' untranslated region. Here, we present the solution NMR structure of CUG-binding protein 2 RRM3 in complex with 5'-UUUAA-3' originating from the COX-2 3'-UTR. We show that RRM3 uses the same binding surface and protein moieties to interact with AU- and UG-rich RNA motifs, binding with low and high affinity, respectively. Using NMR spectroscopy, isothermal titration calorimetry and molecular dynamics simulations, we demonstrate that distinct sub-states characterized by different aromatic side-chain conformations at the RNA-binding surface allow for high- or low-affinity binding with functional implications. This study highlights a mechanism for RNA discrimination possibly common to multiple RRMs as several prominent members display a similar rearrangement of aromatic residues upon binding their targets.The RNA Recognition Motif (RRM) is the most ubiquitous RNA binding domain. Here the authors combined NMR and molecular dynamics simulations and show that the RRM RNA binding surface exists in different states and that a conformational switch of aromatic side-chains fine-tunes sequence specific binding affinities.

• MeSH
• 3' nepřekládaná oblast MeSH
• aminokyselinové motivy MeSH
• CELF proteiny chemie   genetika   metabolismus   MeSH
• cyklooxygenasa 2 genetika   MeSH
• fenylalanin chemie   metabolismus   MeSH
• konformace proteinů MeSH
• magnetická rezonanční spektroskopie MeSH
• messenger RNA chemie   metabolismus   MeSH
• proteiny nervové tkáně chemie   genetika   metabolismus   MeSH
• simulace molekulární dynamiky MeSH
• substituce aminokyselin MeSH
• úseky bohaté na AU MeSH
• vazebná místa MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• 3' nepřekládaná oblast MeSH
• CELF proteiny MeSH
• CELF2 protein, human MeSH Prohlížeč
• cyklooxygenasa 2 MeSH
• fenylalanin MeSH
• messenger RNA MeSH
• proteiny nervové tkáně MeSH
• PTGS2 protein, human MeSH Prohlížeč

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.

• MeSH
• krystalografie rentgenová MeSH
• lidé MeSH
• motiv rozpoznávající RNA genetika   MeSH
• mutageneze cílená MeSH
• nukleární magnetická rezonance biomolekulární MeSH
• rekombinantní proteiny chemie   genetika   metabolismus   MeSH
• RNA metabolismus   MeSH
• sestřihové faktory chemie   genetika   metabolismus   MeSH
• simulace molekulární dynamiky MeSH
• substituce aminokyselin MeSH
• vazebná místa MeSH
• voda chemie   MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• RBFOX1 protein, human MeSH Prohlížeč
• rekombinantní proteiny MeSH
• RNA MeSH
• sestřihové faktory MeSH
• voda MeSH

Replication of human immunodeficiency virus 1 (HIV-1) involves conversion of its single-stranded RNA genome to double-stranded DNA, which is integrated into the genome of the host. This conversion is catalyzed by reverse transcriptase (RT), which possesses DNA polymerase and RNase H domains. The available crystal structures suggest that at any given time the RNA/DNA substrate interacts with only one active site of the two domains of HIV-1 RT. Unknown is whether a simultaneous interaction of the substrate with polymerase and RNase H active sites is possible. Therefore, the mechanism of the coordination of the two activities is not fully understood. We performed molecular dynamics simulations to obtain a conformation of the complex in which the unwound RNA/DNA substrate simultaneously interacts with the polymerase and RNase H active sites. When the RNA/DNA hybrid was immobilized at the polymerase active site, RNase H cleavage occurred, experimentally verifying that the substrate can simultaneously interact with both active sites. These findings demonstrate the existence of a transient conformation of the HIV-1 RT substrate complex, which is important for modulating and coordinating the enzymatic activities of HIV-1 RT.

• MeSH
• DNA chemie   metabolismus   MeSH
• HIV reverzní transkriptasa chemie   metabolismus   MeSH
• katalytická doména MeSH
• ribonukleasa H chemie   metabolismus   MeSH
• RNA chemie   metabolismus   MeSH
• simulace molekulární dynamiky MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• DNA MeSH
• HIV reverzní transkriptasa MeSH
• reverse transcriptase, Human immunodeficiency virus 1 MeSH Prohlížeč
• ribonukleasa H MeSH
• RNA MeSH