LNA-induced dynamic stability in a therapeutic aptamer: insights from molecular dynamics simulations.

Modulation of structural and thermodynamic properties of nucleic acids with synthetic modifications is a promising area of research with possible applications in nanotechnology and nanotherapeutics. Locked nucleic acid (LNA) is one such modification in which the C4' and O2' atoms of the sugar moiety are connected through a methylene bridge. The LNA modified DNA aptamer RNV66, and its unmodified counterpart V7t1, both of which target the vascular endothelial growth factor (VEGF) implicated in oncogenic angiogenesis, have a G-rich tract that can fold into G-quadruplex structures. However, it is not understood why V7t1 has a polymorphic structure while its LNA modified counterpart RNV66 has a unique quadruplex fold with higher nuclease resistance, thermal stability and greater binding affinity for VEGF. In this work, we have performed extensive molecular dynamics simulations of RNV66 and V7t1 to study and compare the structural and dynamic consequences of the insertion of LNAs. It was observed that the increase in dynamic stability was significant in the presence of LNA residues and our protocol for combining different torsional parameters using OL15 for the DNA aptamer and parm99_LNA along with parmbsc0 and ßOL15 for the LNAs nicely reproduced the experimentally observed conformational features of RNV66. Our observations would help in further theoretical studies in understanding the lack of frustration in the folding of the LNA modified aptamer and its higher affinity for VEGF.Communicated by Ramaswamy H. Sarma.

Inosine and its methyl derivatives: Occurrence, biogenesis, and function in RNA.

Inosine is one of the most common post-transcriptional modifications. Since its discovery, it has been noted for its ability to contribute to non-Watson-Crick interactions within RNA. Rapidly accumulating evidence points to the widespread generation of inosine through hydrolytic deamination of adenosine to inosine by different classes of adenosine deaminases. Three naturally occurring methyl derivatives of inosine, i.e., 1-methylinosine, 2'-O-methylinosine and 1,2'-O-dimethylinosine are currently reported in RNA modification databases. These modifications are expected to lead to changes in the structure, folding, dynamics, stability and functions of RNA. The importance of the modifications is indicated by the strong conservation of the modifying enzymes across organisms. The structure, binding and catalytic mechanism of the adenosine deaminases have been well-studied, but the underlying mechanism of the catalytic reaction is not very clear yet. Here we extensively review the existing data on the occurrence, biogenesis and functions of inosine and its methyl derivatives in RNA. We also included the structural and thermodynamic aspects of these modifications in our review to provide a detailed and integrated discussion on the consequences of A-to-I editing in RNA and the contribution of different structural and thermodynamic studies in understanding its role in RNA. We also highlight the importance of further studies for a better understanding of the mechanisms of the different classes of deamination reactions. Further investigation of the structural and thermodynamic consequences and functions of these modifications in RNA should provide more useful information about their role in different diseases.

Topology-based classification of tetrads and quadruplex structures.

MOTIVATION: Quadruplexes attract the attention of researchers from many fields of bio-science. Due to a specific structure, these tertiary motifs are involved in various biological processes. They are also promising therapeutic targets in many strategies of drug development, including anticancer and neurological disease treatment. The uniqueness and diversity of their forms cause that quadruplexes show great potential in novel biological applications. The existing approaches for quadruplex analysis are based on sequence or 3D structure features and address canonical motifs only. RESULTS: In our study, we analyzed tetrads and quadruplexes contained in nucleic acid molecules deposited in Protein Data Bank. Focusing on their secondary structure topology, we adjusted its graphical diagram and proposed new dot-bracket and arc representations. We defined the novel classification of these motifs. It can handle both canonical and non-canonical cases. Based on this new taxonomy, we implemented a method that automatically recognizes the types of tetrads and quadruplexes occurring as unimolecular structures. Finally, we conducted a statistical analysis of these motifs found in experimentally determined nucleic acid structures in relation to the new classification. AVAILABILITY AND IMPLEMENTATION: https://github.com/tzok/eltetrado/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Computational and NMR studies of RNA duplexes with an internal pseudouridine-adenosine base pair.

Pseudouridine (Ψ) is the most common chemical modification present in RNA. In general, Ψ increases the thermodynamic stability of RNA. However, the degree of stabilization depends on the sequence and structural context. To explain experimentally observed sequence dependence of the effect of Ψ on the thermodynamic stability of RNA duplexes, we investigated the structure, dynamics and hydration of RNA duplexes with an internal Ψ-A base pair in different nearest-neighbor sequence contexts. The structures of two RNA duplexes containing 5'-GΨC/3'-CAG and 5'-CΨG/3'-GAC motifs were determined using NMR spectroscopy. To gain insight into the effect of Ψ on duplex dynamics and hydration, we performed molecular dynamics (MD) simulations of RNA duplexes with 5'-GΨC/3'-CAG, 5'-CΨG/3'-GAC, 5'-AΨU/3'-UAA and 5'-UΨA/3'-AAU motifs and their unmodified counterparts. Our results showed a subtle impact from Ψ modification on the structure and dynamics of the RNA duplexes studied. The MD simulations confirmed the change in hydration pattern when U is replaced with Ψ. Quantum chemical calculations showed that the replacement of U with Ψ affected the intrinsic stacking energies at the base pair steps depending on the sequence context. The calculated intrinsic stacking energies help to explain the experimentally observed sequence dependent changes in the duplex stability from Ψ modification.

RNA-Puzzles Round III: 3D RNA structure prediction of five riboswitches and one ribozyme.

RNA-Puzzles is a collective experiment in blind 3D RNA structure prediction. We report here a third round of RNA-Puzzles. Five puzzles, 4, 8, 12, 13, 14, all structures of riboswitch aptamers and puzzle 7, a ribozyme structure, are included in this round of the experiment. The riboswitch structures include biological binding sites for small molecules (S-adenosyl methionine, cyclic diadenosine monophosphate, 5-amino 4-imidazole carboxamide riboside 5'-triphosphate, glutamine) and proteins (YbxF), and one set describes large conformational changes between ligand-free and ligand-bound states. The Varkud satellite ribozyme is the most recently solved structure of a known large ribozyme. All puzzles have established biological functions and require structural understanding to appreciate their molecular mechanisms. Through the use of fast-track experimental data, including multidimensional chemical mapping, and accurate prediction of RNA secondary structure, a large portion of the contacts in 3D have been predicted correctly leading to similar topologies for the top ranking predictions. Template-based and homology-derived predictions could predict structures to particularly high accuracies. However, achieving biological insights from de novo prediction of RNA 3D structures still depends on the size and complexity of the RNA. Blind computational predictions of RNA structures already appear to provide useful structural information in many cases. Similar to the previous RNA-Puzzles Round II experiment, the prediction of non-Watson-Crick interactions and the observed high atomic clash scores reveal a notable need for an algorithm of improvement. All prediction models and assessment results are available at http://ahsoka.u-strasbg.fr/rnapuzzles/.

Dynamics and stability of GCAA tetraloops with 2-aminopurine and purine substitutions.

The contributions of various interactions in the GGCGCAAGCC hairpin containing a GCAA tetraloop were studied by computer simulations using the substitutions of functional groups. The guanosine (G) in the first tetraloop position or in the C-G closing base pair was replaced by 2-aminopurine (AP), and the individual tetraloop's adenosines (A) were replaced by purine (PUR). These substitutions eliminated particular hydrogen bonds thought to stabilize the GCAA tetraloop. For each substitution, molecular dynamics (MD) simulations were carried out in an aqueous solution with sodium counterions, using the CHARMM27 force field. The MD simulations showed that the substitutions in the first (G-->AP) and the third (A-->PUR) position of the GCAA tetraloop did not significantly influence the conformation of the hairpin. A long-lived bridging water molecule observed in the GCAA loop was present in both modified loops. The substitutions made in the last loop position (A-->PUR) or in the C-G base pair closing the tetraloop (G-->AP) to some extent influenced the loop structure and dynamics. These loops did not display the long-lived bridging water molecules. When the second A in the GCAA loop was replaced by PUR, the first A in the loop was observed in the anti or in the syn orientation about the glycosyl bond. The G to AP substitution in C-G base pair led to a change of their arrangement from the Watson-Crick to wobble. The MD simulations of the hairpin with C-AP wobble closing base pair showed increased conformational dynamics of the hairpin. The changes of hairpin formation free energy associated with the substitutions of individual bases were calculated by the free energy perturbation method. Our theoretical estimates suggest a larger destabilization for the G to AP substitutions in GCAA loop than for the substitutions of individual A's by PUR, which is in accordance with experimental tendency. The calculations predicted a similar free energy change for G to AP substitutions in the GCAA tetraloop and in the C-G closing base pair.

Effects of base substitutions in an RNA hairpin from molecular dynamics and free energy simulations.

Contributions of individual interactions in the GGCGCAAGCC hairpin containing a GCAA tetraloop were studied by computer simulations using base substitutions. The G in the first tetraloop position was replaced by inosine (I) or adenosine (A), and the G in the C-G basepair closing the tetraloop was replaced by I. These substitutions eliminate particular hydrogen bonds proposed in the nuclear magnetic resonance model of the GCAA tetraloop. Molecular dynamics simulations of the GCAA tetraloop in aqueous solvent displayed a well-defined hydrogen pattern between the first and last loop nucleotides (G and A) stabilized by a bridging water molecule. Substitution of G-->I in the basepair closing the tetraloop did not significantly influence the loop structure and dynamics. The ICAA loop maintained the overall structure, but displayed variation in the hydrogen-bond network within the tetraloop itself. Molecular dynamics simulations of the ACAA loop led to conformational heterogeneity of the resulting structures. Changes of hairpin formation free energy associated with substitutions of individual bases were calculated by the free energy perturbation method. The calculated decrease of the hairpin stability upon G-->I substitution in the C-G basepair closing the tetraloop was in good agreement with experimental thermodynamic data. Our theoretical estimates for G-->I and G-->A mutations located in the tetraloop suggest larger loop destabilization than corresponding experimental results. The extent of conformational sampling of the structures resulting from base substitutions and its impact on the calculated free energy was discussed.