Search

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.

The RNA chaperone Hfq plays crucial roles in bacterial gene expression and is a major facilitator of small regulatory RNA (sRNA) action. The toroidal architecture of the Hfq hexamer presents three well-characterized surfaces that allow it to bind sRNAs to stabilize them and engage target transcripts. Hfq-interacting sRNAs are categorized into two classes based on the surfaces they use to bind Hfq. By characterizing a systematic alanine mutant library of Hfq to identify amino acid residues that impact survival of Escherichia coli experiencing nitrogen (N) starvation, we corroborated the important role of the three RNA-binding surfaces for Hfq function. We uncovered two, previously uncharacterized, conserved residues, V22 and G34, in the hydrophobic core of Hfq, to have a profound impact on Hfq's RNA-binding activity in vivo. Transcriptome-scale analysis revealed that V22A and G34A Hfq mutants cause widespread destabilization of both sRNA classes, to the same extent as seen in bacteria devoid of Hfq. However, the alanine substitutions at these residues resulted in only modest alteration in stability and structure of Hfq. We propose that V22 and G34 have impact on Hfq function, especially critical under cellular conditions when there is an increased demand for Hfq, such as N starvation., (© The Author(s) 2025. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Molecular dynamics (MD) simulations are an important and well-established tool for investigating RNA structural dynamics, but their accuracy relies heavily on the quality of the employed force field ( ff ). In this work, we present a comprehensive evaluation of widely used pair-additive and polarizable RNA ff s using the challenging UUCG tetraloop (TL) benchmark system. Extensive standard MD simulations, initiated from the NMR structure of the 14-mer UUCG TL, revealed that most ff s did not maintain the native state, instead favoring alternative loop conformations. Notably, three very recent variants of pair-additive ff s, OL3 CP -gHBfix21, DES-Amber, and OL3 R2.7 , successfully preserved the native structure over a 10 × 20 μs time scale. To further assess these ff s, we performed enhanced sampling folding simulations of the shorter 8-mer UUCG TL, starting from the single-stranded conformation. Estimated folding free energies (Δ G ° fold ) varied significantly among these three ff s, with values of 0.0 ± 0.6, 2.4 ± 0.8, and 7.4 ± 0.2 kcal/mol for OL3 CP -gHBfix21, DES-Amber, and OL3 R2.7 , respectively. The Δ G ° fold value predicted by the OL3 CP -gHBfix21 ff was closest to experimental estimates, ranging from -1.6 to -0.7 kcal/mol. In contrast, the higher Δ G ° fold values obtained using DES-Amber and OL3 R2.7 were unexpected, suggesting that key interactions are inaccurately described in the folded, unfolded, or misfolded ensembles. These discrepancies led us to further test DES-Amber and OL3 R2.7 ff s on additional RNA and DNA systems, where further performance issues were observed. Our results emphasize the complexity of accurately modeling RNA dynamics and suggest that creating an RNA ff capable of reliably performing across a wide range of RNA systems remains extremely challenging. In conclusion, our study provides valuable insights into the capabilities of current RNA ff s and highlights key areas for future ff development.

Repression of msl-2 mRNA translation is essential for viability of Drosophila melanogaster females to prevent hypertranscription of both X chromosomes. This translational control event is coordinated by the female-specific protein Sex-lethal (Sxl) which recruits the RNA binding proteins Unr and Hrp48 to the 3' untranslated region (UTR) of the msl-2 transcript and represses translation initiation. The mechanism exerted by Hrp48 during translation repression and its interaction with msl-2 are not well understood. Here we investigate the RNA binding specificity and affinity of the tandem RNA recognition motifs of Hrp48. Using NMR spectroscopy, molecular dynamics simulations and isothermal titration calorimetry, we identified the exact region of msl-2 3' UTR recognized by Hrp48. Additional biophysical experiments and translation assays give further insights into complex formation of Hrp48, Unr, Sxl and RNA. Our results show that Hrp48 binds independent of Sxl and Unr downstream of the E and F binding sites of Sxl and Unr to msl-2., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

The N-terminal RNA recognition motif domain (RRM1) of polypyrimidine tract binding protein (PTB) forms an additional C-terminal helix α3, which docks to one edge of the β-sheet upon binding to a stem-loop RNA containing a UCUUU pentaloop. Importantly, α3 does not contact the RNA. The α3 helix therefore represents an allosteric means to regulate the conformation of adjacent domains in PTB upon binding structured RNAs. Here we investigate the process of dynamic adaptation by stem-loop RNA and RRM1 using NMR and MD in order to obtain mechanistic insights on how this allostery is achieved. Relaxation data and NMR structure determination of the free protein show that α3 is partially ordered and interacts with the domain transiently. Stem-loop RNA binding quenches fast time scale dynamics and α3 becomes ordered, however microsecond dynamics at the protein-RNA interface is observed. MD shows how RRM1 binding to the stem-loop RNA is coupled to the stabilization of the C-terminal helix and helps to transduce differences in RNA loop sequence into changes in α3 length and order. IRES assays of full length PTB and a mutant with altered dynamics in the α3 region show that this dynamic allostery influences PTB function in cultured HEK293T cells., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

N6-Methyladenosine (m 6 A) is a prevalent RNA post-transcriptional modification that plays crucial roles in RNA stability, structural dynamics, and interactions with proteins. The YT521-B (YTH) family of proteins, which are notable m 6 A readers, functions through its highly conserved YTH domain. Recent structural investigations and molecular dynamics (MD) simulations have shed light on the mechanism of recognition of m 6 A by the YTHDC1 protein. Despite advancements, using MD to predict the stabilization induced by m 6 A on the free energy of binding between RNA and YTH proteins remains challenging due to inaccuracy of the employed force field and limited sampling. For instance, simulations often fail to sufficiently capture the hydration dynamics of the binding pocket. This study addresses these challenges through an innovative methodology that integrates metadynamics, alchemical simulations, and force-field refinement. Importantly, our research identifies hydration of the binding pocket as giving only a minor contribution to the binding free energy and emphasizes the critical importance of precisely tuning force-field parameters to experimental data. By employing a fitting strategy built on alchemical calculations, we refine the m 6 A partial charge parameters, thereby enabling the simultaneous reproduction of N6 methylation on both the protein binding free energy and the thermodynamic stability of nine RNA duplexes. Our findings underscore the sensitivity of binding free energies to partial charges, highlighting the necessity for thorough parametrization and validation against experimental observations across a range of structural contexts.

Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The CHARMM force field accurately describes the C3'-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations and also highlight potential pitfalls and best practices for obtaining optimal performance.

The single-stranded RNA genome of SARS-CoV-2 is highly structured. Numerous helical stem-loop structures interrupted by mismatch motifs are present in the functionally important 5'- and 3'-UTRs. These mismatches modulate local helical geometries and feature unusual arrays of hydrogen bonding donor and acceptor groups. However, their conformational and dynamical properties cannot be directly inferred from chemical probing and are difficult to predict theoretically. A mismatch motif (SL1-motif) consisting of three consecutive U•U base pairs is located in stem-loop 1 of the 3'-UTR. We combined NMR-spectroscopy and MD-simulations to investigate its structure and dynamics. All three U•U base pairs feature two direct hydrogen bonds and are as stable as Watson-Crick A:U base pairs. Plasmodium falciparum 25S rRNA contains a triple U•U mismatch motif (Pf-motif) differing from SL1-motif only with respect to the orientation of the two closing base pairs. Interestingly, while the geometry of the outer two U•U mismatches was identical in both motifs the preferred orientation of the central U•U mismatch was different. MD simulations and potassium ion titrations revealed that the potassium ion-binding mode to the major groove is connected to the different preferred geometries of the central base pair in the two motifs., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Guanine quadruplex (GQ) is a noncanonical nucleic acid structure formed by guanine-rich DNA and RNA sequences. Folding of GQs is a complex process, where several aspects remain elusive, despite being important for understanding structure formation and biological functions of GQs. Pulling experiments are a common tool for acquiring insights into the folding landscape of GQs. Herein, we applied a computational pulling strategy─steered molecular dynamics (SMD) simulations─in combination with standard molecular dynamics (MD) simulations to explore the unfolding landscapes of tetrameric parallel GQs. We identified anisotropic properties of elastic conformational changes, unfolding transitions, and GQ mechanical stabilities. Using a special set of structural parameters, we found that the vertical component of pulling force (perpendicular to the average G-quartet plane) plays a significant role in disrupting GQ structures and weakening their mechanical stabilities. We demonstrated that the magnitude of the vertical force component depends on the pulling anchor positions and the number of G-quartets. Typical unfolding transitions for tetrameric parallel GQs involve base unzipping, opening of the G-stem, strand slippage, and rotation to cross-like structures. The unzipping was detected as the first and dominant unfolding event, and it usually started at the 3'-end. Furthermore, results from both SMD and standard MD simulations indicate that partial spiral conformations serve as a transient ensemble during the (un)folding of GQs.

Molecular dynamics (MD) simulations represent an established tool to study RNA molecules. The outcome of MD studies depends, however, on the quality of the force field ( ff ). Here we suggest a correction for the widely used AMBER OL3 ff by adding a simple adjustment of the nonbonded parameters. The reparameterization of the Lennard-Jones potential for the -H8···O5'- and -H6···O5'- atom pairs addresses an intranucleotide steric clash occurring in the type 0 base-phosphate interaction (0BPh). The nonbonded fix (NBfix) modification of 0BPh interactions (NBfix 0BPh modification) was tuned via a reweighting approach and subsequently tested using an extensive set of standard and enhanced sampling simulations of both unstructured and folded RNA motifs. The modification corrects minor but visible intranucleotide clash for the anti nucleobase conformation. We observed that structural ensembles of small RNA benchmark motifs simulated with the NBfix 0BPh modification provide better agreement with experiments. No side effects of the modification were observed in standard simulations of larger structured RNA motifs. We suggest that the combination of OL3 RNA ff and NBfix 0BPh modification is a viable option to improve RNA MD simulations.

Guanine quadruplexes (GQs) are non-canonical nucleic acid structures involved in many biological processes. GQs formed in single-stranded regions often need to be unwound by cellular machinery, so their mechanochemical properties are important. Here, we performed steered molecular dynamics simulations of human telomeric GQs to study their unfolding. We examined four pulling regimes, including a very slow setup with pulling velocity and force load accessible to high-speed atomic force microscopy. We identified multiple factors affecting the unfolding mechanism, i.e.,: (i) the more the direction of force was perpendicular to the GQ channel axis (determined by GQ topology), the more the base unzipping mechanism happened, (ii) the more parallel the direction of force was, GQ opening and cross-like GQs were more likely to occur, (iii) strand slippage mechanism was possible for GQs with an all- anti pattern in a strand, and (iv) slower pulling velocity led to richer structural dynamics with sampling of more intermediates and partial refolding events. We also identified that a GQ may eventually unfold after a force drop under forces smaller than those that the GQ withstood before the drop. Finally, we found out that different unfolding intermediates could have very similar chain end-to-end distances, which reveals some limitations of structural interpretations of single-molecule spectroscopic data.

Pseudouridine is the most frequently naturally occurring RNA modification, found in all classes of biologically functional RNAs. Compared to uridine, pseudouridine contains an additional hydrogen bond donor group and is therefore widely regarded as a structure stabilizing modification. However, the effects of pseudouridine modifications on the structure and dynamics of RNAs have so far only been investigated in a limited number of different structural contexts. Here, we introduced pseudouridine modifications into the U-turn motif and the adjacent U:U closing base pair of the neomycin-sensing riboswitch (NSR)-an extensively characterized model system for RNA structure, ligand binding, and dynamics. We show that the effects of replacing specific uridines with pseudouridines on RNA dynamics crucially depend on the exact location of the replacement site and can range from destabilizing to locally or even globally stabilizing. By using a combination of NMR spectroscopy, MD simulations and QM calculations, we rationalize the observed effects on a structural and dynamical level. Our results will help to better understand and predict the consequences of pseudouridine modifications on the structure and function of biologically important RNAs., (© 2023 Vögele et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Holliday junction (HJ) is a noncanonical four-way DNA structure with a prominent role in DNA repair, recombination, and DNA nanotechnology. By rearranging its four arms, HJ can adopt either closed or open state. With enzymes typically recognizing only a single state, acquiring detailed knowledge of the rearrangement process is an important step toward fully understanding the biological function of HJs. Here, we carried out standard all-atom molecular dynamics (MD) simulations of the spontaneous opening-closing transitions, which revealed complex conformational transitions of HJs with an involvement of previously unconsidered "half-closed" intermediates. Detailed free-energy landscapes of the transitions were obtained by sophisticated enhanced sampling simulations. Because the force field overstabilizes the closed conformation of HJs, we developed a system-specific modification which for the first time allows the observation of spontaneous opening-closing HJ transitions in unbiased MD simulations and opens the possibilities for more accurate HJ computational studies of biological processes and nanomaterials.

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., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

RNA-protein complexes use diverse binding strategies, ranging from structurally well-defined interfaces to completely disordered regions. Experimental characterization of flexible segments is challenging and can be aided by atomistic molecular dynamics (MD) simulations. Here, we used an extended set of microsecond-scale MD trajectories (400 μs in total) to study two FUS-RNA constructs previously characterized by nuclear magnetic resonance (NMR) spectroscopy. The FUS protein contains a well-structured RNA recognition motif domain followed by a presumably disordered RGG tail that binds RNA stem-loop hairpins. Our simulations not only provide several suggestions complementing the experiments but also reveal major methodological difficulties in studies of such complex RNA-protein interfaces. Despite efforts to stabilize the binding via system-specific force-field adjustments, we have observed progressive distortions of the RNA-protein interface inconsistent with experimental data. We propose that the dynamics is so rich that its converged description is not achievable even upon stabilizing the system. Still, after careful analysis of the trajectories, we have made several suggestions regarding the binding. We identify substates in the RNA loops, which can explain the NMR data. The RGG tail localized in the minor groove remains disordered, sampling countless transient interactions with the RNA. There are long-range couplings among the different elements contributing to the recognition, which can lead to allosteric communication throughout the system. Overall, the RNA-FUS systems form dynamical ensembles that cannot be fully represented by single static structures. Thus, albeit imperfect, MD simulations represent a viable tool to investigate dynamic RNA-protein complexes.

The capability of current force fields to reproduce RNA structural dynamics is limited. Several methods have been developed to take advantage of experimental data in order to enforce agreement with experiments. Here, we extend an existing framework which allows arbitrarily chosen force-field correction terms to be fitted by quantification of the discrepancy between observables back-calculated from simulation and corresponding experiments. We apply a robust regularization protocol to avoid overfitting and additionally introduce and compare a number of different regularization strategies, namely, L1, L2, Kish size, relative Kish size, and relative entropy penalties. The training set includes a GACC tetramer as well as more challenging systems, namely, gcGAGAgc and gcUUCGgc RNA tetraloops. Specific intramolecular hydrogen bonds in the AMBER RNA force field are corrected with automatically determined parameters that we call gHBfix opt . A validation involving a separate simulation of a system present in the training set (gcUUCGgc) and new systems not seen during training (CAAU and UUUU tetramers) displays improvements regarding the native population of the tetraloop as well as good agreement with NMR experiments for tetramers when using the new parameters. Then, we simulate folded RNAs (a kink-turn and L1 stalk rRNA) including hydrogen bond types not sufficiently present in the training set. This allows a final modification of the parameter set which is named gHBfix21 and is suggested to be applicable to a wider range of RNA systems.

The programmable synthesis of rationally engineered crystal architectures for the precise arrangement of molecular species is a foundational goal in nanotechnology, and DNA has become one of the most prominent molecules for the construction of these materials. In particular, branched DNA junctions have been used as the central building block for the assembly of 3D lattices. Here, crystallography is used to probe the effect of all 36 immobile Holliday junction sequences on self-assembling DNA crystals. Contrary to the established paradigm in the field, most junctions yield crystals, with some enhancing the resolution or resulting in unique crystal symmetries. Unexpectedly, even the sequence adjacent to the junction has a significant effect on the crystal assemblies. Six of the immobile junction sequences are completely resistant to crystallization and thus deemed "fatal," and molecular dynamics simulations reveal that these junctions invariably lack two discrete ion binding sites that are pivotal for crystal formation. The structures and dynamics detailed here could be used to inform future designs of both crystals and DNA nanostructures more broadly, and have potential implications for the molecular engineering of applied nanoelectronics, nanophotonics, and catalysis within the crystalline context., (© 2022. The Author(s).)

The YTH domain of YTHDC1 belongs to a class of protein "readers", recognizing the N6-methyladenosine (m 6 A) chemical modification in mRNA. Static ensemble-averaged structures revealed details of N6-methyl recognition via a conserved aromatic cage. Here, we performed molecular dynamics (MD) simulations along with nuclear magnetic resonance (NMR) and isothermal titration calorimetry (ITC) to examine how dynamics and solvent interactions contribute to the m 6 A recognition and negative selectivity toward an unmethylated substrate. The structured water molecules surrounding the bound RNA and the methylated substrate's ability to exclude bulk water molecules contribute to the YTH domain's preference for m 6 A. Intrusions of bulk water deep into the binding pocket disrupt binding of unmethylated adenosine. The YTHDC1's preference for the 5'-Gm 6 A-3' motif is partially facilitated by a network of water-mediated interactions between the 2-amino group of the guanosine and residues in the m 6 A binding pocket. The 5'-Im 6 A-3' (where I is inosine) motif can be recognized too, but disruption of the water network lowers affinity. The D479A mutant also disrupts the water network and destabilizes m 6 A binding. Our interdisciplinary study of the YTHDC1 protein-RNA complex reveals an unusual physical mechanism by which solvent interactions contribute toward m 6 A recognition.

Phosphorothioates (PTs) are important chemical modifications of the RNA backbone where a single nonbridging oxygen of the phosphate is replaced with a sulfur atom. PT can stabilize RNAs by protecting them from hydrolysis and is commonly used as a tool to explore their function. It is, however, unclear what basic physical effects PT has on RNA stability and electronic structure. Here, we present molecular dynamics (MD) simulations, quantum mechanical (QM) calculations, and NMR spectroscopy measurements, exploring the effects of PT modifications in the structural context of the neomycin-sensing riboswitch (NSR). The NSR is the smallest biologically functional riboswitch with a well-defined structure stabilized by a U-turn motif. Three of the signature interactions of the U-turn: an H-bond, an anion-π interaction, and a potassium binding site; are formed by RNA phosphates, making the NSR an ideal model for studying how PT affects RNA structure and dynamics. By comparing with high-level QM calculations, we reveal the distinct physical properties of the individual interactions facilitated by the PT. The sulfur substitution, besides weakening the direct H-bond interaction, reduces the directionality of H-bonding while increasing its dispersion and induction components. It also reduces the induction and increases the dispersion component of the anion-π stacking. The sulfur force-field parameters commonly employed in the literature do not reflect these distinctions, leading to the unsatisfactory description of PT in simulations of the NSR. We show that it is not possible to accurately describe the PT interactions using one universal set of van der Waals sulfur parameters and provide suggestions for improving the force-field performance.

The human prototypical SR protein SRSF1 is an oncoprotein that contains two RRMs and plays a pivotal role in RNA metabolism. We determined the structure of the RRM1 bound to RNA and found that the domain binds preferentially to a CN motif (N is for any nucleotide). Based on this solution structure, we engineered a protein containing a single glutamate to asparagine mutation (E87N), which gains the ability to bind to uridines and thereby activates SMN exon7 inclusion, a strategy that is used to cure spinal muscular atrophy. Finally, we revealed that the flexible inter-RRM linker of SRSF1 allows RRM1 to bind RNA on both sides of RRM2 binding site. Besides revealing an unexpected bimodal mode of interaction of SRSF1 with RNA, which will be of interest to design new therapeutic strategies, this study brings a new perspective on the mode of action of SRSF1 in cells.

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., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: The chromodomain (CD) of HP1 proteins is an established H3K9 me3 reader that also binds H1, EHMT2 and H3K23 lysine-methylated targets. Structural experiments have provided atomistic pictures of its recognition of the conserved ARK me3 S/T motif, but structural dynamics' contribution to the recognition may have been masked by ensemble averaging., Methods: We acquired ~350 μs of explicit solvent molecular dynamics (MD) simulations of the CD domain interacting with several peptides using the latest AMBER force fields., Results: The simulations reproduced the experimentally observed static binding patterns well but also revealed visible structural dynamics at the interfaces. While the buried K 0 me3 and A -2 target residues are tightly bound, several flanking sidechains sample diverse sites on the CD surface. Different amino acid positions of the targets can substitute for each other by forming mutually replaceable interactions with CD, thereby explaining the lack of strict requirement for cationic H3 target residues at the -3 position. The Q -4 residue of H3 targets further stabilizes the binding. The recognition pattern of the H3K23 ATK me3 A motif, for which no structure is available, is predicted., Conclusions: The CD reads a longer target segment than previously thought, ranging from positions -7 to +3. The CD anionic clamp can be neutralized not only by the -3 and -1 residues, but also by -7, -6, -5 and +3 residues., General Significance: Structural dynamics, not immediately apparent from the structural data, contribute to molecular recognition between the HP1 CD domain and its targets. Mutual replaceability of target residues increases target sequence flexibility., (Copyright © 2020. Published by Elsevier B.V.)

Explicit solvent atomistic molecular dynamics (MD) simulations represent an established technique to study structural dynamics of RNA molecules and an important complement for diverse experimental methods. However, performance of molecular mechanical (MM) force fields (ff's) remains far from satisfactory even after decades of development, as apparent from a problematic structural description of some important RNA motifs. Actually, some of the smallest RNA molecules belong to the most challenging systems for MD simulations and, among them, the UUCG tetraloop is saliently difficult. We report a detailed analysis of UUCG MD simulations, depicting the sequence of events leading to the loss of the UUCG native state during MD simulations. The total amount of MD simulation data analyzed in this work is close to 1.3 ms. We identify molecular interactions, backbone conformations, and substates that are involved in the process. Then, we unravel specific ff deficiencies using diverse quantum mechanical/molecular mechanical (QM/MM) and QM calculations. Comparison between the MM and QM methods shows discrepancies in the description of the 5'-flanking phosphate moiety and both signature sugar-base interactions. Our work indicates that poor behavior of the UUCG tetraloop in simulations is a complex issue that cannot be attributed to one dominant and straightforwardly correctable factor. Instead, there is a concerted effect of multiple ff inaccuracies that are coupled and amplifying each other. We attempted to improve the simulation behavior by some carefully tailored interventions, but the results were still far from satisfactory, underlying the difficulties in development of accurate nucleic acid ff's.

Human heterochromatin protein 1 (HP1) is a key factor in heterochromatin formation and maintenance. Its chromo-shadow domain (CSD) homodimerizes, and the HP1 dimer acts as a hub, transiently interacting with diverse binding partner (BP) proteins. We analyze atomistic details of interactions of the HP1γ(CSD) dimer with one of its targets, the histone H3 N-terminal tail, using molecular dynamics (MD) simulations. The goal is to complement the available X-ray crystallography data and unravel potential dynamic effects in the molecular recognition. Our results suggest that HP1(CSD)-BP recognition involves structural dynamics of both partners, including structural communication between adjacent binding pockets that may fine-tune the sequence recognition. For example, HP1 Trp174 sidechain substates may help in distinguishing residues bound in the conserved HP1(CSD) ±2 hydrophobic pockets. Further, there is intricate competition between the binding of negatively charged HP1 C-terminal extension and solvent anions near the ±2 hydrophobic pockets, which is also influenced by the BP sequence. Phosphorylated H3 Y41 can interact with the same site. We also analyze the ability of several pair-additive force fields to describe the protein-protein interface. ff14SB and ff99SB-ILDN* provide the closest correspondence with the crystallographic model. The ff15ipq local dynamics are somewhat less consistent with details of the experimental structure, while larger perturbations of the interface commonly occur in CHARMM36m simulations. The balance of some interactions, mainly around the anion binding site, also depends on the ion parameters. Some differences between the simulated and experimental structures are attributable to crystal packing.

The family of heterochromatin protein 1 (HP1) isoforms is essential for chromatin packaging, regulation of gene expression, and repair of damaged DNA. Here we document that γ-radiation reduced the number of HP1α-positive foci, but not HP1β and HP1γ foci, located in the vicinity of the fibrillarin-positive region of the nucleolus. The additional analysis confirmed that γ-radiation has the ability to significantly decrease the level of HP1α in rDNA promoter and rDNA encoding 28S rRNA. By mass spectrometry, we showed that treatment by γ-rays enhanced the HP1β serine 88 phosphorylation (S88ph), but other analyzed modifications of HP1β, including S161ph/Y163ph, S171ph, and S174ph, were not changed in cells exposed to γ-rays or treated by the HDAC inhibitor (HDACi). Interestingly, a combination of HDACi and γ-radiation increased the level of HP1α and HP1γ. The level of HP1β remained identical before and after the HDACi/γ-rays treatment, but HDACi strengthened HP1β interaction with the KRAB-associated protein 1 (KAP1) protein. Conversely, HP1γ did not interact with KAP1, although approximately 40% of HP1γ foci co-localized with accumulated KAP1. Especially HP1γ foci at the periphery of nucleoli were mostly absent of KAP1. Together, DNA damage changed the morphology, levels, and interaction properties of HP1 isoforms. Also, γ-irradiation-induced hyperphosphorylation of the HP1β protein; thus, HP1β-S88ph could be considered as an important marker of DNA damage., Competing Interests: The authors declare no conflict of interest.

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.

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.

Human antigen R (HuR) is a key regulator of cellular mRNAs containing adenylate/uridylate-rich elements (AU-rich elements; AREs). These are a major class of cis elements within 3' untranslated regions, targeting these mRNAs for rapid degradation. HuR contains three RNA recognition motifs (RRMs): a tandem RRM1 and 2, followed by a flexible linker and a C-terminal RRM3. While RRM1 and 2 are structurally characterized, little is known about RRM3. Here we present a 1.9-Å-resolution crystal structure of RRM3 bound to different ARE motifs. This structure together with biophysical methods and cell-culture assays revealed the mechanism of RRM3 ARE recognition and dimerization. While multiple RNA motifs can be bound, recognition of the canonical AUUUA pentameric motif is possible by binding to two registers. Additionally, RRM3 forms homodimers to increase its RNA binding affinity. Finally, although HuR stabilizes ARE-containing RNAs, we found that RRM3 counteracts this effect, as shown in a cell-based ARE reporter assay and by qPCR with native HuR mRNA targets containing multiple AUUUA motifs, possibly by competing with RRM12., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

Understanding the RNA binding specificity of protein is of primary interest to decipher their function in the cell. Here, we review the methodology used to solve the structures of protein-RNA complexes using solution-state NMR spectroscopy: from sample preparation to structure calculation procedures. We also describe how molecular dynamics simulations can help providing additional information on the role of key amino acid side chains and of water molecules in protein-RNA recognition., (© 2019 Elsevier Inc. All rights reserved.)

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., Competing Interests: The authors have declared that no competing interests exist.

The NMR solution structures of human telomeric (Htel) G-quadruplexes (GQs) are characterized by the presence of two lateral loops complemented by either diagonal or propeller loops. Bases of a given loop can establish interactions within the loop as well as with other loops and the flanking bases. This can lead to a formation of base alignments above and below the GQ stems. These base alignments are known to affect the loop structures and relative stabilities of different Htel GQ folds. We have carried out a total of 217 μs of classical (unbiased) molecular dynamics (MD) simulations starting from the available solution structures of Htel GQs to characterize structural dynamics of the lateral and diagonal loops, using several recent AMBER DNA force-field variants. As the loops are involved in diverse stacking and H-bonding interactions, their dynamics is slow, and extended sampling is required to capture different conformations. Nevertheless, although the simulations are far from being quantitatively converged, the data suggest that multiple 10 μs-scale simulations can provide a quite good assessment of the loop conformational space as described by the force field. The simulations indicate that the lateral loops may sample multiple coexisting conformations, which should be considered when comparing simulations with the NMR models as the latter include ensemble averaging. The adenine-thymine Watson-Crick arrangement was the most stable base pairing in the simulations. Adenine-adenine and thymine-thymine base pairs were also sampled but were less stable. The data suggest that the description of lateral and diagonal GQ loops in contemporary MD simulations is considerably more realistic than the description of propeller loops, though definitely not flawless.

Although atomistic explicit-solvent Molecular Dynamics (MD) is a popular tool to study protein-RNA recognition, satisfactory MD description of protein-RNA complexes is not always achieved. Unfortunately, it is often difficult to separate MD simulation instabilities primarily caused by the simple point-charge molecular mechanics (MM) force fields from problems related to the notorious uncertainties in the starting structures. Herein, we report a series of large-scale QM/MM calculations on the U1A protein-RNA complex. This experimentally well-characterized system has an intricate protein-RNA interface, which is very unstable in MD simulations. The QM/MM calculations identify several H-bonds poorly described by the MM method and thus indicate the sources of instabilities of the U1A interface in MD simulations. The results suggest that advanced QM/MM computations could be used to indirectly rationalize problems seen in MM-based MD simulations of protein-RNA complexes. As the most accurate QM method, we employ the computationally demanding meta-GGA density functional TPSS-D3(BJ)/def2-TZVP level of theory. Because considerably faster methods would be needed to extend sampling and to study even larger protein-RNA interfaces, a set of low-cost QM/MM methods is compared to the TPSS-D3(BJ)/def2-TZVP data. The PBEh-3c and B97-3c density functional composite methods appear to be suitable for protein-RNA interfaces. In contrast, HF-3c and the tight-binding Hamiltonians DFTB3-D3 and GFN-xTB perform unsatisfactorily and do not provide any advantage over the MM description. These conclusions are supported also by similar analysis of a simple HutP protein-RNA interface, which is well-described by MD with the exception of just one H-bond. Some other methodological aspects of QM/MM calculations on protein-RNA interfaces are discussed.

The neomycin sensing riboswitch is the smallest biologically functional RNA riboswitch, forming a hairpin capped with a U-turn loop-a well-known RNA motif containing a conserved uracil. It was shown previously that a U→C substitution of the eponymous conserved uracil does not alter the riboswitch structure due to C protonation at N3. Furthermore, cytosine is evolutionary permitted to replace uracil in other U-turns. Here, we use molecular dynamics simulations to study the molecular basis of this substitution in the neomycin sensing riboswitch and show that a structure-stabilizing monovalent cation-binding site in the wild-type RNA is the main reason for its negligible structural effect. We then use NMR spectroscopy to confirm the existence of this cation-binding site and to demonstrate its effects on RNA stability. Lastly, using quantum chemical calculations, we show that the cation-binding site is altering the electronic environment of the wild-type U-turn so that it is more similar to the cytosine mutant. The study reveals an amazingly complex and delicate interplay between various energy contributions shaping up the 3D structure and evolution of nucleic acids.

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.

HIV-1 reverse transcriptase (RT) possesses both DNA polymerase activity and RNase H activity that act in concert to convert single-stranded RNA of the viral genome to double-stranded DNA that is then integrated into the DNA of the infected cell. Reverse transcriptase-catalyzed reverse transcription critically relies on the proper generation of a polypurine tract (PPT) primer. However, the mechanism of PPT primer generation and the features of the PPT sequence that are critical for its recognition by HIV-1 RT remain unclear. Here, we used a chemical cross-linking method together with molecular dynamics simulations and single-molecule assays to study the mechanism of PPT primer generation. We found that the PPT was specifically and properly recognized within covalently tethered HIV-1 RT-nucleic acid complexes. These findings indicated that recognition of the PPT occurs within a stable catalytic complex after its formation. We found that this unique recognition is based on two complementary elements that rely on the PPT sequence: RNase H sequence preference and incompatibility of the poly(rA/dT) tract of the PPT with the nucleic acid conformation that is required for RNase H cleavage. The latter results from rigidity of the poly(rA/dT) tract and leads to base-pair slippage of this sequence upon deformation into a catalytically relevant geometry. In summary, our results reveal an unexpected mechanism of PPT primer generation based on specific dynamic properties of the poly(rA/dT) segment and help advance our understanding of the mechanisms in viral RNA reverse transcription., (© 2018 by The American Society for Biochemistry and Molecular Biology, Inc.)

The HutP protein from B. subtilis regulates histidine metabolism by interacting with an antiterminator mRNA hairpin in response to the binding of l-histidine and Mg 2+ . We studied the functional ligand-bound HutP hexamer complexed with two mRNAs using all-atom microsecond-scale explicit-solvent MD simulations performed with the Amber force fields. The experimentally observed protein-RNA interface exhibited good structural stability in the simulations with the exception of some fluctuations in an unusual adenine-threonine interaction involving two closely spaced H-bonds. We further investigated this interaction by comparing QM/MM and MM optimizations, using the QM region comprising almost 350 atoms described at the DFT-D3 level. The QM/MM method clearly improved the adenine-threonine interaction compared to MM, especially when the X-H bond lengths were frozen during the MM optimization to mimic the use of SHAKE in the MD simulations. Thus, both the MM approximation and the use of SHAKE can compromise the description of H-bonds at protein-RNA interfaces. The simulations also revealed a notable Mg 2+ -parameter dependence in the behavior of the ligand-binding pocket (LBP). With the SPC/E water model, the 12-6 Åqvist and Li&Merz parameters provided an entirely stable LBP structure, but the 12-6 Allnér and 12-6-4 Li&Merz parametrizations resulted in a progressive loss of direct nitrogen-Mg 2+ LBP coordination. The Allnér and Li&Merz 12-6 parametrizations were also tested with the TIP3P water model; the LBP was destabilized in both cases. This illustrates the difficulty of consistently describing different Mg 2+ interactions using nonpolarizable force fields. Overall, the simulations support the hypothesis that HutP protein becomes fully structured upon ligand binding. Subsequent RNA binding does not affect the protein structure, in keeping with the mechanism inferred from experimental structures.

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.

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., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (© 2016 Wiley Periodicals, Inc.)

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., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

We report an extensive set of explicit solvent molecular dynamics (MD) simulations (∼25 μs of accumulated simulation time) of the RNA kissing-loop complex of the HIV-1 virus initiation dimerization site. Despite many structural investigations by X-ray, NMR, and MD techniques, the position of the bulged purines of the kissing complex has not been unambiguously resolved. The X-ray structures consistently show bulged-out positions of the unpaired bases, while several NMR studies show bulged-in conformations. The NMR studies are, however, mutually inconsistent regarding the exact orientations of the bases. The earlier simulation studies predicted the bulged-out conformation; however, this finding could have been biased by the short simulation time scales. Our microsecond-long simulations reveal that all unpaired bases of the kissing-loop complex stay preferably in the interior of the kissing-loop complex. The MD results are discussed in the context of the available experimental data and we suggest that both conformations are biochemically relevant. We also show that MD provides a quite satisfactory description of this RNA system, contrasting recent reports of unsatisfactory performance of the RNA force fields for smaller systems such as tetranucleotides and tetraloops. We explain this by the fact that the kissing complex is primarily stabilized by an extensive network of Watson-Crick interactions which are rather well described by the force fields. We tested several different sets of water/ion parameters but they all lead to consistent results. However, we demonstrate that a recently suggested modification of van der Waals interactions of the Cornell et al. force field deteriorates the description of the kissing complex by the loss of key stacking interactions stabilizing the interhelical junction and excessive hydrogen-bonding interactions.

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., (© The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (Copyright © 2014 Elsevier B.V. All rights reserved.)

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.