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.

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.

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.

The chromodomain (CD) of HP1 proteins is an established H3K9me3 reader that also binds H1, EHMT2 and H3K23 lysine-methylated targets. Structural experiments have provided atomistic pictures of its recognition of the conserved ARKme3S/T motif, but structural dynamics' contribution to the recognition may have been masked by ensemble averaging.

Examination of thermal decomposition of street samples of cocaine and methamphetamine shows that typical products detected in previous studies are accompanied by a wide palette of simple volatile compounds easily detectable by spectral techniques. These molecules increase smoke toxicity and their spectral detection can be potentially used for identification of drug samples by well-controlled laboratory thermolysis in temperature progression. In our study, street samples of cocaine and methamphetamine have been thermolyzed under vacuum over the temperature range of 350-650 °C. The volatile products (CO, HCN, CH4, C2H4, etc.) have been monitored by high-resolution Fourier-transform infrared (FTIR) spectrometry in this temperature range. The decomposition mechanism has been additionally examined theoretically by quantum-chemical calculations for the highest temperature achieved experimentally in our study and beyond. Prior to analysis, the street samples have also been characterized by FTIR, Raman spectroscopy, energy-dispersive X-ray spectroscopy, and melting point determination.

Accumulation and selection of nucleotides is one of the most challenging problems surrounding the origin of the first RNA molecules on our planet. In the current work we propose that guanosine 3',5' cyclic monophosphate could selectively crystallize upon evaporation of an acidic prebiotic pool containing various other nucleotides. The conditions of the evaporative crystallization are fully compatible with the subsequent acid catalyzed polymerization of this cyclic nucleotide reported in earlier studies and may be relevant in a broad range of possible prebiotic environments. Albeit cytidine 3',5' cyclic monophosphate has the ability to selectively accumulate under the same conditions, its crystal structure is not likely to support polymer formation.

Guanine-rich DNA sequences tend to form four-stranded G-quadruplex structures. Characteristic glycosidic conformational patterns along the G-strands, such as the 5'-syn-anti-syn-anti pattern observed with the Oxytricha nova telomeric G-quadruplexes, have been well documented. However, an explanation for these featured glycosidic patterns has not emerged. This work presents MD simulation and free energetic analyses for simplified two-quartet [d(GG)](4) models and suggests that the four base pair step patterns show quite different relative stabilities: syn-anti > anti-anti > anti-syn > syn-syn. This suggests the following rule: when folding, anti-parallel G-quadruplexes tend to maximize the number of syn-anti steps and avoid the unfavorable anti-syn and syn-syn steps. This rule is consistent with most of the anti-parallel G-quadruplex structures in the Protein Databank (PDB). Structural polymorphisms of G-quadruplexes relate to these glycosidic conformational patterns and the lengths of the G-tracts. The folding topologies of G2- and G4-tracts are not very polymorphic because each strand tends to populate the stable syn-anti repeat. G3-tracts, on the other hand, cannot present this repeating pattern on each G-tract. This leads to smaller energy differences between different geometries and helps explain the extreme structural polymorphism of the human telomeric G-quadruplexes.

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 gHBfixopt. 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 widespread Mn2+-sensing yybP-ykoY riboswitch controls the expression of bacterial Mn2+ homeostasis genes. Here, we first determine the crystal structure of the ligand-bound yybP-ykoY riboswitch aptamer from Xanthomonas oryzae at 2.96 Å resolution, revealing two conformations with docked four-way junction (4WJ) and incompletely coordinated metal ions. In >100 µs of MD simulations, we observe that loss of divalents from the core triggers local structural perturbations in the adjacent docking interface, laying the foundation for signal transduction to the regulatory switch helix. Using single-molecule FRET, we unveil a previously unobserved extended 4WJ conformation that samples transient docked states in the presence of Mg2+. Only upon adding sub-millimolar Mn2+, however, can the 4WJ dock stably, a feature lost upon mutation of an adenosine contacting Mn2+ in the core. These observations illuminate how subtly differing ligand preferences of competing metal ions become amplified by the coupling of local with global RNA dynamics.

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.

RNA molecules play a key role in countless biochemical processes. RNA interactions, which are of highly diverse nature, are determined by the fact that RNA is a highly negatively charged polyelectrolyte, which leads to intimate interactions with an ion atmosphere. Although RNA molecules are formally single-stranded, canonical (Watson-Crick) duplexes are key components of folded RNAs. A double-stranded (ds) RNA is also important for the design of RNA-based nanostructures and assemblies. Despite the fact that the description of canonical dsRNA is considered the least problematic part of RNA modeling, the imperfect shape and flexibility of dsRNA can lead to imbalances in the simulations of larger RNAs and RNA-containing assemblies. We present a comprehensive set of molecular dynamics (MD) simulations of four canonical A-RNA duplexes. Our focus was directed toward the characterization of the influence of varying ion concentrations and of the size of the solvation box. We compared several water models and four RNA force fields. The simulations showed that the A-RNA shape was most sensitive to the RNA force field, with some force fields leading to a reduced inclination of the A-RNA duplexes. The ions and water models played a minor role. The effect of the box size was negligible, and even boxes with a small fraction of the bulk solvent outside the RNA hydration sphere were sufficient for the simulation of the dsRNA.

A-tracts are functionally important DNA sequences which induce helix bending and have peculiar structural properties. While A-tract structure has been qualitatively well characterized, their mechanical properties remain controversial. A-tracts appear structurally rigid and resist nucleosome formation, but seem flexible in DNA looping. In this work, we investigate mechanical properties of symmetric AnTn and asymmetric A2n tracts for n = 3, 4, 5 using two types of coarse-grained models. The first model represents DNA as an ensemble of interacting rigid bases with non-local quadratic deformation energy, the second one treats DNA as an anisotropically bendable and twistable elastic rod. Parameters for both models are inferred from microsecond long, atomic-resolution molecular dynamics simulations. We find that asymmetric A-tracts are more rigid than the control G/C-rich sequence in localized distortions relevant for nucleosome formation, but are more flexible in global bending and twisting relevant for looping. The symmetric tracts, in contrast, are more rigid than asymmetric tracts and the control, both locally and globally. Our results can reconcile the contradictory stiffness data on A-tracts and suggest symmetric A-tracts to be more efficient in nucleosome exclusion than the asymmetric ones. This would open a new possibility of gene expression manipulation using A-tracts.

DNA G-hairpins are potential key structures participating in folding of human telomeric guanine quadruplexes (GQ). We examined their properties by standard MD simulations starting from the folded state and long T-REMD starting from the unfolded state, accumulating ∼130 μs of atomistic simulations. Antiparallel G-hairpins should spontaneously form in all stages of the folding to support lateral and diagonal loops, with sub-μs scale rearrangements between them. We found no clear predisposition for direct folding into specific GQ topologies with specific syn/anti patterns. Our key prediction stemming from the T-REMD is that an ideal unfolded ensemble of the full GQ sequence populates all 4096 syn/anti combinations of its four G-stretches. The simulations can propose idealized folding pathways but we explain that such few-state pathways may be misleading. In the context of the available experimental data, the simulations strongly suggest that the GQ folding could be best understood by the kinetic partitioning mechanism with a set of deep competing minima on the folding landscape, with only a small fraction of molecules directly folding to the native fold. The landscape should further include non-specific collapse processes where the molecules move via diffusion and consecutive random rare transitions, which could, e.g. structure the propeller loops.

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.

Cyclobutane pyrimidine dimers (CpDs) are among the most common DNA lesions occurring due to the interaction with ultraviolet light. While photolyases have been well known as external factors repairing CpDs, the intrinsic self-repairing capabilities of the GAT[double bond, length as m-dash]T DNA sequence were discovered only recently and are still largely obscure. Here, we elucidate the mechanistic details of this self-repair process by means of MD simulations and QM/MM computations involving the algebraic diagrammatic construction to the second order [ADC(2)] method. We show that local UV-excitation of guanine may be followed by up to three subsequent electron transfers, which may eventually enable efficient CpD ring opening when the negative charge resides on the T[double bond, length as m-dash]T dimer. Consequently, the molecular mechanism of GAT[double bond, length as m-dash]T self-repair can be envisaged as sequential electron transfer (SET) occurring downhill along the slope of the S1 potential energy surface. Even though the general features of the SET mechanism are retained in both of the studied stacked conformers, our optimizations of different S1/S0 state crossings revealed minor differences which could influence their self-repair efficiencies. We expect that such assessment of the availability and efficiency of the SET process in other DNA oligomers could hint towards other sequences exhibiting similar photochemical properties. Such explorations will be particularly fascinating in the context of the origins of biomolecules on Earth, owing to the lack of external repairing factors in the Archean age.

Genomic sequences susceptible to form G-quadruplexes (G4s) are always flanked by other nucleotides, but G4 formation in vitro is generally studied with short synthetic DNA or RNA oligonucleotides, for which bases adjacent to the G4 core are often omitted. Herein, we systematically studied the effects of flanking nucleotides on structural polymorphism of 371 different oligodeoxynucleotides that adopt intramolecular G4 structures. We found out that the addition of nucleotides favors the formation of a parallel fold, defined as the 'flanking effect' in this work. This 'flanking effect' was more pronounced when nucleotides were added at the 5'-end, and depended on loop arrangement. NMR experiments and molecular dynamics simulations revealed that flanking sequences at the 5'-end abolish a strong syn-specific hydrogen bond commonly found in non-parallel conformations, thus favoring a parallel topology. These analyses pave a new way for more accurate prediction of DNA G4 folding in a physiological context.

Guanine radical cation (G•+) is a key intermediate in many oxidative processes occurring in nucleic acids. Here, by combining mixed Quantum Mechanical/Molecular Mechanics calculations and Molecular Dynamics (MD) simulations, we study how the structural behaviour of a tract GGG(TTAGGG)3 (hereafter Tel21) of the human telomeric sequence, folded in an antiparallel quadruple helix, changes when one of the G bases is ionized to G•+ (Tel21+). Once assessed that the electron-hole is localized on a single G, we perform MD simulations of twelve Tel21+ systems, differing in the position of G•+ in the sequence. When G•+ is located in the tetrad adjacent to the diagonal loop, we observe substantial structural rearrangements, which can decrease the electrostatic repulsion with the inner Na+ ions and increase the solvent exposed surface of G•+. Analysis of solvation patterns of G•+ provides new insights on the main reactions of G•+, i.e. the deprotonation at two different sites and hydration at the C8 atom, the first steps of the processes producing 8oxo-Guanine. We suggest the main structural determinants of the relative reactivity of each position and our conclusions, consistent with the available experimental trends, can help rationalizing the reactivity of other G-quadruplex topologies.

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.

The hepatitis delta virus (HDV) ribozyme is a member of the class of small, self-cleaving catalytic RNAs found in a wide range of genomes from HDV to human. Both pre- and post-catalysis (precursor and product) crystal structures of the cis-acting genomic HDV ribozyme have been determined. These structures, together with extensive solution probing, have suggested that a significant conformational change accompanies catalysis. A recent crystal structure of a trans-acting precursor, obtained at low pH and by molecular replacement from the previous product conformation, conforms to the product, raising the possibility that it represents an activated conformer past the conformational change. Here, using fluorescence resonance energy transfer (FRET), we discovered that cleavage of this ribozyme at physiological pH is accompanied by a structural lengthening in magnitude comparable to previous trans-acting HDV ribozymes. Conformational heterogeneity observed by FRET in solution appears to have been removed upon crystallization. Analysis of a total of 1.8 µsec of molecular dynamics (MD) simulations showed that the crystallographically unresolved cleavage site conformation is likely correctly modeled after the hammerhead ribozyme, but that crystal contacts and the removal of several 2'-oxygens near the scissile phosphate compromise catalytic in-line fitness. A cis-acting version of the ribozyme exhibits a more dynamic active site, while a G-1 residue upstream of the scissile phosphate favors poor fitness, allowing us to rationalize corresponding changes in catalytic activity. Based on these data, we propose that the available crystal structures of the HDV ribozyme represent intermediates on an overall rugged RNA folding free-energy landscape.

The utility of molecular dynamics (MD) simulations to model biomolecular structure, dynamics, and interactions has witnessed enormous advances in recent years due to the availability of optimized MD software and access to significant computational power, including GPU multicore computing engines and other specialized hardware. This has led researchers to routinely extend conformational sampling times to the microsecond level and beyond. The extended sampling time has allowed the community not only to converge conformational ensembles through complete sampling but also to discover deficiencies and overcome problems with the force fields. Accuracy of the force fields is a key component, along with sampling, toward being able to generate accurate and stable structures of biopolymers. The Amber force field for nucleic acids has been used extensively since the 1990s, and multiple artifacts have been discovered, corrected, and reassessed by different research groups. We present a direct comparison of two of the most recent and state-of-the-art Amber force field modifications, bsc1 and OL15, that focus on accurate modeling of double-stranded DNA. After extensive MD simulations with five test cases and two different water models, we conclude that both modifications are a remarkable improvement over the previous bsc0 force field. Both force field modifications show better agreement when compared to experimental structures. To ensure convergence, the Drew-Dickerson dodecamer (DDD) system was simulated using 100 independent MD simulations, each extended to at least 10 μs, and the independent MD simulations were concatenated into a single 1 ms long trajectory for each combination of force field and water model. This is significantly beyond the time scale needed to converge the conformational ensemble of the internal portions of a DNA helix absent internal base pair opening. Considering all of the simulations discussed in the current work, the MD simulations performed to assess and validate the current force fields and water models aggregate over 14 ms of simulation time. The results suggest that both the bsc1 and OL15 force fields render average structures that deviate significantly less than 1 Å from the average experimental structures. This can be compared to similar but less exhaustive simulations with the CHARMM 36 force field that aggregate to the ∼90 μs time scale and also perform well but do not produce structures as close to the DDD NMR average structures (with root-mean-square deviations of 1.3 Å) as the newer Amber force fields. On the basis of these analyses, any future research involving double-stranded DNA simulations using the Amber force fields should employ the bsc1 or OL15 modification.