Contact networks in RNA: a structural bioinformatics study with a new tool

When will RNA get its AlphaFold moment?

The protein structure prediction problem has been solved for many types of proteins by AlphaFold. Recently, there has been considerable excitement to build off the success of AlphaFold and predict the 3D structures of RNAs. RNA prediction methods use a variety of techniques, from physics-based to machine learning approaches. We believe that there are challenges preventing the successful development of deep learning-based methods like AlphaFold for RNA in the short term. Broadly speaking, the challenges are the limited number of structures and alignments making data-hungry deep learning methods unlikely to succeed. Additionally, there are several issues with the existing structure and sequence data, as they are often of insufficient quality, highly biased and missing key information. Here, we discuss these challenges in detail and suggest some steps to remedy the situation. We believe that it is possible to create an accurate RNA structure prediction method, but it will require solving several data quality and volume issues, usage of data beyond simple sequence alignments, or the development of new less data-hungry machine learning methods.

Lone Pair…π Contacts and Structure Signatures of r(UNCG) Tetraloops, Z-Turns, and Z-Steps: A WebFR3D Survey

Z-DNA and Z-RNA have long appeared as oddities to nucleic acid scientists. However, their Z-step constituents are recurrently observed in all types of nucleic acid systems including ribosomes. Z-steps are NpN steps that are isostructural to Z-DNA CpG steps. Among their structural features, Z-steps are characterized by the presence of a lone pair…π contact that involves the stacking of the ribose O4′ atom of the first nucleotide with the 3′-face of the second nucleotide. Recently, it has been documented that the CpG step of the ubiquitous r(UNCG) tetraloops is a Z-step. Accordingly, such r(UNCG) conformations were called Z-turns. It has also been recognized that an r(GAAA) tetraloop in appropriate conditions can shapeshift to an unusual Z-turn conformation embedding an ApA Z-step. In this report, we explore the multiplicity of RNA motifs based on Z-steps by using the WebFR3D tool to which we added functionalities to be able to retrieve motifs containing lone pair…π contacts. Many examples that underscore the diversity and universality of these motifs are provided as well as tutorial guidance on using WebFR3D. In addition, this study provides an extensive survey of crystallographic, cryo-EM, NMR, and molecular dynamics studies on r(UNCG) tetraloops with a critical view on how to conduct database searches and exploit their results.

DrawTetrado to create layer diagrams of G4 structures

Motivation

Quadruplexes are specific 3D structures found in nucleic acids. Due to the exceptional properties of these motifs, their exploration with the general-purpose bioinformatics methods can be problematic or insufficient. The same applies to visualizing their structure. A hand-drawn layer diagram is the most common way to represent the quadruplex anatomy. No molecular visualization software generates such a structural model based on atomic coordinates.

Results

DrawTetrado is an open-source Python program for automated visualization targeting the structures of quadruplexes and G4-helices. It generates static layer diagrams that represent structural data in a pseudo-3D perspective. The possibility to set color schemes, nucleotide labels, inter-element distances or angle of view allows for easy customization of the output drawing.

Availability and implementation

The program is available under the MIT license at https://github.com/RNApolis/drawtetrado.

Understanding the role of non-Watson-Crick base pairs in DNA-protein recognition: Structural and energetic aspects using crystallographic database analysis and quantum chemical calculation

Specific recognition of DNA base sequences by proteins is vital for life-cycles of all organisms. In a large number of crystal structures of protein-DNA complexes, DNA conformation significantly deviates from the canonical B-DNA structure. A key question is whether such alternate conformations exist prior to protein binding and one is selected for complexation or the structure observed is induced by protein binding. Non-canonical base pairs, such as Hoogsteen base pairs, are often observed in crystal structures of protein-DNA complexes. We decided to explore whether the occurrence of such non-canonical base pairs in protein-DNA complexes is induced by the protein or is selected from pre-existing conformations. Detailed quantum chemical calculations with dispersion-corrected density functional theory (DFT-D) indicated that most of the non-canonical base pairs with DNA bases are stable even in the absence of the interacting amino acids. However, the G:G Hoogsteen base pair, which also appears in the telomere structure, appears to be unstable in the absence of other stabilizing agents, such as positively charged amino acids. Thus, the stability of many of the non-canonical base pair containing duplexes may be close to the canonical B-DNA structure and hence energetically accessible in the ground state; suggesting that the selection from pre-existing conformations may be an important mechanism for observed non-canonical base pairs in protein-DNA complexes.