Isoalloxazine derivatives promote photocleavage of natural RNAs at G.U base pairs embedded within helices

RNA Secondary Structure Analysis Using RNAstructure

RNAstructure is a user-friendly program for the prediction and analysis of RNA secondary structure. It is available as a web server, a program with a graphical user interface, or a set of command line tools. The programs are available for Microsoft Windows, macOS, or Linux. This article provides protocols for prediction of RNA secondary structure (using the web server, the graphical user interface, or the command line) and high-affinity oligonucleotide binding sites to a structured RNA target (using the graphical user interface). © 2023 Wiley Periodicals LLC. Basic Protocol 1: Predicting RNA secondary structure using the RNAstructure web server Alternate Protocol 1: Predicting secondary structure and base pair probabilities using the RNAstructure graphical user interface Alternate Protocol 2: Predicting secondary structure and base pair probabilities using the RNAstructure command line interface Basic Protocol 2: Predicting binding affinities of oligonucleotides complementary to an RNA target using OligoWalk.

Structure and thermodynamics of a UGG motif interacting with Ba2+ and other metal ions: accommodating changes in the RNA structure and the presence of a G(syn)-G(syn) pair

The self-complementary triplet 5'UGG3'/5'UGG3' is a particular structural motif containing noncanonical G-G pair and two U·G wobble pairs. It constitutes a specific structural and electrostatic environment attracting metal ions, particularly Ba2+ ions. Crystallographic research has shown that two Ba2+ cations are located in the major groove of the helix and interact directly with the UGG triplet. A comparison with the unliganded structure has revealed global changes in the RNA structure in the presence of metal ions, whereas thermodynamic measurements have shown increased stability. Moreover, in the structure with Ba2+, an unusual noncanonical G(syn)-G(syn) pair is observed instead of the common G(anti)-G(syn). We further elucidate the metal binding properties of the UGG/UGG triplet by performing crystallographic and thermodynamic studies using DSC and UV melting with other metal ions. The results explain the preferences of the UGG sequence for Ba2+ cations and point to possible applications of this metal-binding propensity.

The multiple flavors of GoU pairs in RNA

Wobble GU pairs (or GoU) occur frequently within double‐stranded RNA helices interspersed within the standard G═C and A─U Watson‐Crick pairs. However, other types of GoU pairs interacting on their Watson‐Crick edges have been observed. The structural and functional roles of such alternative GoU pairs are surprisingly diverse and reflect the various pairings G and U can form by exploiting all the subtleties of their electronic configurations. Here, the structural characteristics of the GoU pairs are updated following the recent crystallographic structures of functional ribosomal complexes and the development in our understanding of ribosomal translation.

Selective tumor cell death induced by irradiated riboflavin through recognizing DNA G-T mismatch

Riboflavin (vitamin B2) has been thought to be a promising antitumoral agent in photodynamic therapy, though the further application of the method was limited by the unclear molecular mechanism. Our work reveals that riboflavin was able to recognize G–T mismatch specifically and induce single-strand breaks in duplex DNA targets efficiently under irradiation. In the presence of riboflavin, the photo-irradiation could induce the death of tumor cells that are defective in mismatch repair system selectively, highlighting the G–T mismatch as potential drug target for tumor cells. Moreover, riboflavin is a promising leading compound for further drug design due to its inherent specific recognition of the G–T mismatch.

Sequence-specific RNA Photocleavage by Single-stranded DNA in Presence of Riboflavin

Constant efforts have been made to develop new method to realize sequence-specific RNA degradation, which could cause inhibition of the expression of targeted gene. Herein, by using an unmodified short DNA oligonucleotide for sequence recognition and endogenic small molecue, vitamin B2 (riboflavin) as photosensitizer, we report a simple strategy to realize the sequence-specific photocleavage of targeted RNA. The DNA strand is complimentary to the target sequence to form DNA/RNA duplex containing a G•U wobble in the middle. The cleavage reaction goes through oxidative elimination mechanism at the nucleoside downstream of U of the G•U wobble in duplex to obtain unnatural RNA terminal, and the whole process is under tight control by using light as switch, which means the cleavage could be carried out according to specific spatial and temporal requirements. The biocompatibility of this method makes the DNA strand in combination with riboflavin a promising molecular tool for RNA manipulation.

5-Deazaflavin derivatives as inhibitors of p53 ubiquitination by HDM2

Based on previous reports of certain 5-deazaflavin derivatives being capable of activating the tumour suppressor p53 in cancer cells through inhibition of the p53-specific ubiquitin E3 ligase HDM2, we have conducted an structure–activity relationship (SAR) analysis through systematic modification of the 5-deazaflavin template. This analysis shows that HDM2-inhibitory activity depends on a combination of factors. The most active compounds (e.g., 15) contain a trifluoromethyl or chloro substituent at the deazaflavin C9 position and this activity depends to a large extent on the presence of at least one additional halogen or methyl substituent of the phenyl group at N10. Our SAR results, in combination with the HDM2 RING domain receptor recognition model we present, form the basis for the design of drug-like and potent activators of p53 for potential cancer therapy.

RNA self-cleavage activated by ultraviolet light-induced oxidation

A novel UV-C-light-induced ribozyme activity was discovered within the highly structured 5'-genomic regions of both Hepatitis C Virus (HCV) and the related Classic Swine Fever Virus (CSFV). Cleavage is mediated by exposure to UV-C light but not by exogenous oxygen radicals. It is also very selective, occurring at base positions HCV C(79) and CSFV A(45) in some molecules and at the immediately adjacent 5'-positions HCV U(78) and CSFV U(44) in others. Among other reaction products, the majority of biochemically active products detected contained 3'-phosphate and 5'-phosphate-end groups at the newly generated termini, along with a much lower amount of 3'-hydroxyl end group. While preservation of an E-loop RNA structure in the vicinity of the cleavage site was a requisite for HCV RNA self-cleavage, this was not the case for CSFV RNA. The short size of the reactive domains (~33 nt), which are compatible with primitive RNA motifs, and the lack of sequence homology, indicate that as-yet unidentified UV-activated ribozymes are likely to be found throughout structured RNAs, thereby providing clues to whether early RNA self-cleavage events were mediated by photosensitive RNA structures.

A multifunctional bioconjugate module for versatile photoaffinity labeling and click chemistry of RNA

A multifunctional reagent based on a coumarin scaffold was developed for derivatization of naive RNA. The alkylating agent N3BC [7-azido-4-(bromomethyl)coumarin], obtained by Pechmann condensation, is selective for uridine. N3BC and its RNA conjugates are pre-fluorophores which permits controlled modular and stepwise RNA derivatization. The success of RNA alkylation by N3BC can be monitored by photolysis of the azido moiety, which generates a coumarin fluorophore that can be excited with UV light of 320 nm. The azidocoumarin-modified RNA can be flexibly employed in structure-function studies. Versatile applications include direct use in photo-crosslinking studies to cognate proteins, as demonstrated with tRNA and RNA fragments from the MS2 phage and the HIV genome. Alternatively, the azide function can be used for further derivatization by click-chemistry. This allows e.g. the introduction of an additional fluorophore for excitation with visible light.

Small-molecule inhibitors of MDM2 as new anticancer therapeutics

It has long been known that traditional anticancer radio- and chemotherapies in part work through direct or indirect activation of the p53 tumour suppressor pathway. However, many of these strategies are nonselective and genotoxic. The emerging understanding of the pathways that regulate p53 has led to the notion that it should be possible to activate the p53 pathway in ways that are inherently nongenotoxic. Important targets for pharmacological interference in this respect are MDM2 and MDMX, key negative regulators of p53. Genetic and pharmacologic studies suggest that blocking the physical interaction of these proteins with p53, or inhibiting the catalytic role of MDM2 in tagging p53 for proteasomal degradation, both of which lead to an increase in the transcriptional activity of p53, may indeed be an efficient and safe way to eradicate tumour cells that retain wild-type p53. Here we review the rationale for such strategies, as well as the current state in the discovery and development of drugs that reactivate p53 by inhibiting its inhibitors MDM2 and MDMX. The first compounds that have been shown in model systems to be able selectively to kill cancer cells in this way are now entering clinical trials and the promise of MDM2 inhibitors as a new therapeutic anticancer modality should therefore become clear in the not-too-distant future.

RNA secondary structure analysis using RNAstructure

RNAstructure is a user-friendly program for the prediction and analysis of RNA secondary structure under Microsoft Windows. This unit provides protocols for RNA secondary structure prediction and prediction of high-affinity oligonucleotide binding sites to a structured RNA target.

Communication between RNA folding domains revealed by folding of circularly permuted ribozymes

To study the role of sequence and topology in RNA folding, we determined the kinetic folding pathways of two circularly permuted variants of the Tetrahymena group I ribozyme, using time-resolved hydroxyl radical footprinting. Circular permutation changes the distance between interacting residues in the primary sequence, without changing the native structure of the RNA. In the natural ribozyme, tertiary interactions in the P4-P6 domain form in 1 s, while interactions in the P3-P9 form in 1-3 min at 42 degrees C. Permutation of the 5' end to G111 in the P4 helix allowed the stable P4-P6 domain to fold in 200 ms at 30 degrees C, five times faster than in the wild-type RNA, while the other domains folded five times more slowly (5-8 min). By contrast, circular permutation of the 5' end to G303 in J8/7 decreased the folding rate of the P4-P6 domain. In this permuted RNA, regions joining P2, P3 and P4 were protected in 500 ms, while the P3-P9 domain was 60-80% folded within 30 s. RNase T(1) digestion and FMN photocleavage showed that circular permutation of the RNA sequence alters the initial ensemble of secondary structures, thereby changing the tertiary folding pathways. Our results show that the natural 5'-to-3' order of the structural domains in group I ribozymes optimizes structural communication between tertiary domains and promotes self-assembly of the catalytic center.

Reaction mechanisms of riboflavin triplet state with nucleic acid bases

ESR and laser flash photolysis studies have determined a reasonable order of reactivity of nucleotides with triplet riboflavin (3Rb*) for the first time. ESR detection of triplet state reactivity of Rb with nucleoside, polynucleotide and DNA has been obtained simultaneously. In addition, ESR spin elimination measurement of the reactivity of 3Rb* with nucleotides in good accord with laser flash photolysis determination of the corresponding rate constants offers a simple and reliable method to detect the reactivities of nucleic acids and its components with photoexcited flavins. Kinetic, ESR and thermodynamic studies have demonstrated that Rb should be a strong endogenous photosensitizer capable of oxidizing all nucleic acid bases, and preferentially two purine nucleotides with high rate constants.

Using an RNA secondary structure partition function to determine confidence in base pairs predicted by free energy minimization

A partition function calculation for RNA secondary structure is presented that uses a current set of nearest neighbor parameters for conformational free energy at 37 degrees C, including coaxial stacking. For a diverse database of RNA sequences, base pairs in the predicted minimum free energy structure that are predicted by the partition function to have high base pairing probability have a significantly higher positive predictive value for known base pairs. For example, the average positive predictive value, 65.8%, is increased to 91.0% when only base pairs with probability of 0.99 or above are considered. The quality of base pair predictions can also be increased by the addition of experimentally determined constraints, including enzymatic cleavage, flavin mono-nucleotide cleavage, and chemical modification. Predicted secondary structures can be color annotated to demonstrate pairs with high probability that are therefore well determined as compared to base pairs with lower probability of pairing.

Chemical and functional diversity of small molecule ligands for RNA

Functional RNAs such as ribosomal RNA and structured domains of mRNA are targets for small molecule ligands that can act as modulators of the RNA biological activity. Natural ligands for RNA display a bewildering structural and chemical complexity that has yet to be matched by synthetic RNA binders. Comparison of natural and artificial ligands for RNA may help to direct future approaches to design and synthesize potent novel scaffolds for specific recognition of RNA targets.

The coenzyme thiamine pyrophosphate inhibits the self-splicing of the group I intron

Effects of the coenzyme thiamine pyrophosphate and its analogs on the inhibition of self-splicing of primary transcripts of the phage T4 thymidylate synthase gene (td) were investigated. Of all compounds tested, the coenzyme thiamine pyrophosphate was the most potent inhibitor and the order of inhibitory efficiency for compounds tested was as follows: thiamine pyrophosphate>thiamine monophosphate>thiamine>thiochrome. Increasing guanosine concentration overcame the suppression of self-splicing by thiamine pyrophosphate close to the level of normal splicing. Kinetic analysis demonstrated that thiamine pyrophosphate acts as a competitive inhibitor for the td intron RNA with a Ki of 2.2mM. The splicing specificity inhibition by thiamine pyrophosphate is predominantly due to changes in Km.

Photoinduced cleavage by a rhodium complex at G.U mismatches and exposed guanines in large and small RNAs

Photoinduced cleavage reactions by the rhodium complex tris(4,7-diphenyl-1,10-phenanthroline)rhodium(III) [Rh(DIP)(3)(3+)] with three RNA hairpins, r(GGGGU UCGCUC CACCA) (16 nucleotide, tetraloop(Ala2)), r(GGGGCUAUAGCUCUAGCUC CACCA) (24 nucleotide, microhelix(Ala)), and r(GGCGGUUAGAUAUCGCC) (17 nucleotide, 790 loop), and full-length (1542 nucleotide) 16S rRNA from Escherichia coli were investigated. The cleavage reactions were monitored by gel electrophoresis and the sites of cleavage by Rh(DIP)(3)(3+) were determined by comparisons with chemical or enzymatic sequencing reactions. In general, RNA backbone scission by the metal complex was induced at G.U mismatches and at exposed G residues. The cleavage activity was observed on the three small RNA hairpins as well as on the isolated 1542-nucleotide ribosomal RNA.

Specific cleavage of hyper-edited dsRNAs

Extended double-stranded DNA (dsRNA) duplexes can be hyper-edited by adenosine deaminases that act on RNA (ADARs). Long uninterrupted dsRNA is relatively uncommon in cells, and is frequently associated with infection by DNA or RNA viruses. Moreover, extensive adenosine to inosine editing has been reported for various viruses. A number of cellular antiviral defence strategies are stimulated by dsRNA. An additional mechanism to remove dsRNA from cells may involve hyper-editing of dsRNA by ADARs, followed by targeted cleavage. We describe here a cytoplasmic endonuclease activity that specifically cleaves hyper-edited dsRNA. Cleavage occurs at specific sites consisting of alternating IU and UI base pairs. In contrast, unmodified dsRNA and even deaminated dsRNAs that contain four consecutive IU base pairs are not cleaved. Moreover, dsRNAs in which alternating IU and UI base pairs are replaced by isomorphic GU and UG base pairs are not cleaved. Thus, the cleavage of deaminated dsRNA appears to require an RNA structure that is unique to hyper-edited RNA, providing a molecular target for the disposal of hyper-edited viral RNA.

NAD+ inhibits the self-splicing of the group I intron

We investigated the effects of the coenzyme NAD+ (nicotinamide adenine dinucleotide) and its analogs on the self-splicing of primary transcripts of the phage T4 thymidylate synthase gene (td). Of all the nicotinamide coenzymes and analogs tested, NADP+ was the strongest inhibitor, with a potency approximately threefold that of NAD+. Kinetic analysis demonstrated that NAD+ acts as a mixed type noncompetitive inhibitor for the td intron RNA with a K(i) of 4.1 mM. The splicing specificity inhibition by NAD+ is predominantly due to changes in Km and kcat, and was Mg2+ concentration dependent. The results suggest that both the ADP and nicotinamide moieties are the key structural features in NAD+ responsible for the inhibition of splicing.

The flavin coenzymes: a new class of group I intron inhibitors

Effects of the coenzyme flavin mononucleotide (FMN) and its analogs on the self-splicing of primary transcripts of the phage T4 thymidylate synthase gene (td) have been investigated. Among all flavins and analogs tested, the lumichrome was the most inhibitory. The kinetic analysis demonstrated that FMN acts as a competitive inhibitor for the td intron RNA with a Ki of 1.86 mM although it does not possess a guanidino group in its structure. FMN is able to inhibit the first step of the self-splicing, thus identifying FMN as a novel class of group I intron splicing inhibitors. The specificity of the splicing inhibition by FMN is predominantly due to changes in Km but not k(cat). The splicing inhibition is believed to be due to the interference with the affinity of GTP for the intron RNA. The analysis of the inhibitory concentration and structural examination suggests that the key structural features in FMN responsible for the inhibition of splicing may be an alloxazine group.

Spectinomycin inhibits the self-splicing of the group 1 intron RNA

Effects of the aminoglycoside spectinomycin on the self-splicing of primary transcripts of the phage T4 thymidylate synthase gene (td) have been investigated. The kinetic analysis demonstrated that spectinomycin acts as a mixed noncompetitive inhibitor for the td intron RNA with a K(i) of 7.2 mM. Increasing the spectinomycin concentration raised the K(m) values with the corresponding decrease of V(max) and k(cat) values. The specificity of the splicing inhibition by spectinomycin is due to changes in both K(m) and k(cat). The splicing inhibition by spectinomycin is dependent on pH changes and Mg(2+) concentration, indicating electrostatic interactions with the intron RNA. It has been proposed that the key structural features in spectinomycin responsible for the inhibition of splicing may be the hydroxyl groups on the antibiotic.

On the wobble GoU and related pairs

The wobble GoU pairs have been implicated in several biological processes where RNA molecules play a key role. We review the geometrical and conformational properties of wobble GoU pairs on the basis of available crystal structures of RNAs at high resolution. The similarities with the wobble A+oC pairs and UoU pairs are illustrated, while the differences with the recently discovered bifurcated G x U pairs are contrasted.

A molecular dynamics simulation of the flavin mononucleotide-RNA aptamer complex

We report on an unrestrained molecular dynamics simulation of the flavin mononucleotide (FMN)-RNA aptamer. The simulated average structure maintains both cross-strand and intermolecular FMN-RNA nuclear Overhauser effects from the nmr experiments and has all qualitative features of the nmr structure including the G10-U12-A25 base triple and the A13-G24, A8-G28, and G9-G27 mismatches. However, the relative orientation of the hairpin loop to the remaining part of the molecule differs from the nmr structure. The simulation predicts that the flexible phosphoglycerol part of FMN moves toward G27 and forms hydrogen bonds. There are structurally long-lived water molecules in the FMN binding pocket forming hydrogen bonds within FMN and between FMN and RNA. In addition, long-lived water is found bridging primarily RNA backbone atoms. A general feature of the environment of long-lived "structural" water is at least two and in most cases three or four potential acceptor atoms. The 2'-OH group of RNA usually acts as an acceptor in interactions with the solvent. There are almost no intrastrand O2'H(n) vertical ...O4'(n + 1) hydrogen bonds within the RNA backbone. In the standard case the preferred orientation of the 2'-OH hydrogen atoms is approximately toward O3' of the same nucleotide. However, a relatively large number of conformations with the backbone torsional angle gamma in the trans orientation is found. A survey of all experimental RNA x-ray structures shows that this backbone conformation occurs but is less frequent than found in the simulation. Experimental nmr RNA aptamer structures have a higher fraction of this conformation as compared to the x-ray structures. The backbone conformation of nucleotide n + 1 with the torsional angle gamma in the trans orientation leads to a relatively short distance between 2'-OH(n) and O5'(n + 1), enabling hydrogen-bond formation. In this case the preferred orientation of the 2'-OH hydrogen atom is approximately toward O5'(n + 1). We find two relatively short and dynamically stable types of backbone-backbone next-neighbor contacts, namely C2'(H)(n) vertical ...O4'(n + 1) and C5'(H)(n + 1) vertical ...O2'(n). These interactions may affect both backbone rigidity and thermodynamic stability of RNA helical structures.

Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure

An improved dynamic programming algorithm is reported for RNA secondary structure prediction by free energy minimization. Thermodynamic parameters for the stabilities of secondary structure motifs are revised to include expanded sequence dependence as revealed by recent experiments. Additional algorithmic improvements include reduced search time and storage for multibranch loop free energies and improved imposition of folding constraints. An extended database of 151,503 nt in 955 structures? determined by comparative sequence analysis was assembled to allow optimization of parameters not based on experiments and to test the accuracy of the algorithm. On average, the predicted lowest free energy structure contains 73 % of known base-pairs when domains of fewer than 700 nt are folded; this compares with 64 % accuracy for previous versions of the algorithm and parameters. For a given sequence, a set of 750 generated structures contains one structure that, on average, has 86 % of known base-pairs. Experimental constraints, derived from enzymatic and flavin mononucleotide cleavage, improve the accuracy of structure predictions.

Docking of cationic antibiotics to negatively charged pockets in RNA folds

The binding of aminoglycosides to RNA provides a paradigm system for the analysis of RNA-drug interactions. The electrostatic field around three-dimensional RNA folds creates localized and defined negatively charged regions which are potential docking sites for the cationic ammonium groups of aminoglycosides. To explore in RNA folds the electronegative pockets suitable for aminoglycoside binding, we used calculations of the electrostatic field and Brownian dynamics simulations of cation diffusion. We applied the technique on those RNA molecules experimentally known to bind aminoglycosides, namely, two tobramycin aptamers (Wang, Y.; Rando, R. R. Chem. Biol. 1995, 2, 281-290): the aminoglycoside-binding region in 16S ribosomal RNA (Moazed, S.; Noller, H. F. Nature 1987, 327, 389-394) and the TAR RNA from human immunodeficiency virus (Mei, H.-Y.; et al. Bioorg. Med. Chem. Lett. 1995, 5, 2755-2760). For the aptamers and ribosomal RNA, for which the binding sites of the aminoglycosides are known, a good agreement between negatively charged pockets and the binding positions of the drugs was found. On the basis of variations between neomycin-like and kanamycin-like aminoglycosides in the interaction with the electrostatic field of ribosomal RNA, we propose a model for the different binding specificities of these two classes of drugs. The spatial congruence between the electronegative pockets in RNA folds and binding positions of aminoglycosides was used to dock aminoglycosides to ribosomal and TAR RNAs. Molecular dynamics simulations were used to analyze possible RNA-drug interactions. Aminoglycosides inhibit the binding of the viral Tat protein to TAR RNA; however, the drug-binding sites are still unknown. Thus, our docking approach provides first structural models for TAR-aminoglycoside complexes. The RNA-drug interactions observed in the modeled complexes support the view that the antibiotics might lock TAR in a conformation with low affinity for the Tat protein, explaining the experimentally found aminoglycoside inhibition of the Tat-TAR interaction (Mei, H.-Y.; et al. Bioorg. Med. Chem. Lett. 1995, 5, 2755-2760).

Engineering precision RNA molecular switches

Ligand-specific molecular switches composed of RNA were created by coupling preexisting catalytic and receptor domains via structural bridges. Binding of ligand to the receptor triggers a conformational change within the bridge, and this structural reorganization dictates the activity of the adjoining ribozyme. The modular nature of these tripartite constructs makes possible the rapid construction of precision RNA molecular switches that trigger only in the presence of their corresponding ligand. By using similar enzyme engineering strategies, new RNA switches can be made to operate as designer molecular sensors or as a new class of genetic control elements.

A nickel complex cleaves uridine in folded RNA structures: application to E. coli tmRNA and related engineered molecules

To gain more insight about Escherichia coli tmRNA structure, NiCR, a square planar macrocyclic nickel (II) complex, was used to probe guanine N7 exposure. On the basis of this additional structural information, a refined secondary structure of the molecule is proposed. In addition to its known specificity for guanine N7, we show here that the chemical probe can also cleave at specific uridine residues. In contrast to the alkaline-labile modification of guanine, the reactivity of NiCR at these uridine residues results in direct strand scission. To better characterize the uridine cleavage sites and assess the importance of the RNA structure for the reaction to occur, smaller RNA molecules derived from one pseudoknot (PK4) of E. coli tmRNA containing two uridine cleavage sites were engineered and probed. It is shown that this pseudoknot can fold by itself in solution and that the expected uridine residues are also cleaved by the nickel complex, suggesting that only a local sequence and/or structural context is required for cleavage. In E. coli tmRNA, the five uridine cleavage sites are located in double-stranded regions. These sites contain a G-U wobble base-pair and a downstream uridine which is cleaved. Using smaller RNAs derived from one stem of PK4, systematic changes in the proposed recognition motif indicate that the G-U pair is required for cleavage. Furthermore, there is no cleavage if the G-U pair is reversed. If the recognition motif is moved within the stem, the cleavage site moves accordingly. Additionally, if the recognition motif is changed such that the G-U pair is flanked by two uridine residues, the reactivity occurs only at the 3' uridine. Radical quenching studies have indicated that sulfate radical, as in the case of guanine oxidation, is involved in uridine oxidation. Although additional studies are required to better characterize the reaction, this paper reports a novel specificity for a chemical probe which may be useful for investigating structural motifs involving G-U pairs in folded RNAs.

Simulations of the molecular dynamics of nucleic acids

The growing amount of high quality molecular dynamics simulations generated using the latest methodological developments and force fields has led to a sharper understanding of the forces underlying the dynamics of biomolecular systems, as well as to stimulating insights into the structure and catalysis of nucleic acids. It is now clear that inclusion of long-range electrostatic interactions and of the aqueous and ionic environment is necessary for producing realistic and accurate simulations. Yet, many papers hint at a force field and protocol dependence of the results and thus contain the seeds for the future improvements that will be necessary for deepening our understanding of recognition phenomena and folding of nucleic acids.

RNA as a drug target: chemical, modelling, and evolutionary tools

Dramatic technical progress in RNA synthesis and structure determination has allowed several difficulties inherent to the preparation, handling and structural analysis of RNA to be overcome, and this has led to a wealth of information about RNA structure and its relationship with biological function. It is now fully recognized that RNA molecules intervene at all stages of cell life, not only because of key sequence motifs but also because of intricate three-dimensional folds. This realization has promoted RNA to a potential therapeutic target. As in protein motifs recognizing nucleic acids, groups of the molecule interacting with RNA contribute to specific binding through defined hydrogen bonds and van der Waals docking, while other parts contribute to the driving force of binding via less specific electrostatic interactions accompanied by water and ion displacement.