ONQUADRO: a database of experimentally determined quadruplex structures

RNA-Puzzles Round V: blind predictions of 23 RNA structures

RNA-Puzzles is a collective endeavor dedicated to the advancement and improvement of RNA three-dimensional structure prediction. With agreement from structural biologists, RNA structures are predicted by modeling groups before publication of the experimental structures. We report a large-scale set of predictions by 18 groups for 23 RNA-Puzzles: 4 RNA elements, 2 Aptamers, 4 Viral elements, 5 Ribozymes and 8 Riboswitches. We describe automatic assessment protocols for comparisons between prediction and experiment. Our analyses reveal some critical steps to be overcome to achieve good accuracy in modeling RNA structures: identification of helix-forming pairs and of non-Watson-Crick modules, correct coaxial stacking between helices and avoidance of entanglements. Three of the top four modeling groups in this round also ranked among the top four in the CASP15 contest.

Evaluating the Performance of Peak Calling Algorithms Available for Intracellular G-Quadruplex Sequencing

DNA G-quadruplexes (G4) are non-canonical DNA structures that play key roles in various biological processes. Antibody-dependent sequencing is an important tool for identifying intracellularly formed DNA G4s, and peak calling is a crucial step in processing the sequencing data. As the applicability of existing peak calling algorithms to intracellular G4 data has not been previously assessed, we systematically compared and evaluated these algorithms to determine those best suited for G4 detection. We selected seven representative candidates from 43 published peak calling algorithms for detailed evaluation. The performance of each candidate on six published intracellular G4 sequencing datasets (GSE107690, GSE145090, GSE133379, GSE178668ChIP-seq, GSE178668CUT&Tag, GSE221437) were assessed by precision and recall against customized benchmarks integrating results from multiple algorithms, as well as consistency with known G4 information (pG4 predicted by pqsfinder, oG4 from GSE63874, and multi-cell-line conserved G4s) and epigenetic signals. We identified MACS2, PeakRanger, and GoPeaks as the most effective algorithms for analyzing intracellular G4 sequencing data, and attributed their superior performance partially to the distribution model of sequencing reads/fragments used in the hypothesis testing step of the peak calling procedures. These findings provide guidance and rationale for selecting peak callers appropriate for intracellular G4 data.

Resolving the Structure of a Guanine Quadruplex in TMPRSS2 Messenger RNA by Circular Dichroism and Molecular Modeling

The presence of a guanine quadruplex in the opening reading frame of the messenger RNA coding for the transmembrane serine protease 2 (TMPRSS2) may pave the way to original anticancer and host-oriented antiviral strategy. Indeed, TMPRSS2 in addition to being overexpressed in different cancer types, is also related to the infection of respiratory viruses, including SARS-CoV-2, by promoting the cellular and viral membrane fusion through its proteolytic activity. The design of selective ligands targeting TMPRSS2 messenger RNA requires a detailed knowledge, at atomic level, of its structure. Therefore, we have used an original experimental-computational protocol to predict the first resolved structure of the parallel guanine quadruplex secondary structure in the RNA of TMPRSS2, which shows a rigid core flanked by a flexible loop. This represents the first atomic scale structure of the guanine quadruplex structure present in TMPRSS2 messenger RNA.

Strategy for modeling higher-order G-quadruplex structures recalcitrant to NMR determination

Guanine-rich nucleic acids can form intramolecularly folded four-stranded structures known as G-quadruplexes (G4s). Traditionally, G4 research has focused on short, highly modified DNA or RNA sequences that form well-defined homogeneous compact structures. However, the existence of longer sequences with multiple G4 repeats, from proto-oncogene promoters to telomeres, suggests the potential for more complex higher-order structures with multiple G4 units that might offer selective drug-targeting sites for therapeutic development. These larger structures present significant challenges for structural characterization by traditional high-resolution methods like multi-dimensional NMR and X-ray crystallography due to their molecular complexity. To address this current challenge, we have developed an integrated structural biology (ISB) platform, combining experimental and computational methods to determine self-consistent molecular models of higher-order G4s (xG4s). Here we outline our ISB method using two recent examples from our lab, an extended c-Myc promoter and long human telomere G4 repeats, that highlights the utility and generality of our approach to characterizing biologically relevant xG4s.

Entangled World of DNA Quadruplex Folds

DNA quadruplexes participate in many biological functions. It takes up a variety of folds based on the sequence and environment. Here, a meticulous analysis of experimentally determined 437 quadruplex structures (433 PDBs) deposited in the PDB is carried out. The analysis reveals the modular representation of the quadruplex folds. Forty-eight unique quadruplex motifs (whose diversity arises out of the propeller, bulge, diagonal, and lateral loops that connect the quartets) are identified, leading to simple to complex inter/intramolecular quadruplex folds. The two-layered structural motifs are further classified into 33 continuous and 15 discontinuous motifs. While the continuous motifs can directly be extended to a quadruplex fold, the discontinuous motif requires an additional loop(s) to complete a fold, as illustrated here with examples. Similarly, higher-order quadruplex folds can also be represented by continuous or discontinuous motifs or their combinations. Such a modular representation of the quadruplex folds may assist in custom engineering of quadruplexes, designing motif-based drugs, and the prediction of the quadruplex structure. Furthermore, it could facilitate understanding of the role of quadruplexes in biological functions and diseases.

Crystal structure of a tetrameric RNA G-quadruplex formed by hexanucleotide repeat expansions of C9orf72 in ALS/FTD

The abnormal GGGGCC hexanucleotide repeat expansions (HREs) in C9orf72 cause the fatal neurodegenerative diseases including amyotrophic lateral sclerosis and frontotemporal dementia. The transcribed RNA HREs, short for r(G4C2)n, can form toxic RNA foci which sequestrate RNA binding proteins and impair RNA processing, ultimately leading to neurodegeneration. Here, we determined the crystal structure of r(G4C2)2, which folds into a parallel tetrameric G-quadruplex composed of two four-layer dimeric G-quadruplex via 5'-to-5' stacking in coordination with a K+ ion. Notably, the two C bases locate at 3'- end stack on the outer G-tetrad with the assistance of two additional K+ ions. The high-resolution structure reported here lays a foundation in understanding the mechanism of neurological toxicity of RNA HREs. Furthermore, the atomic details provide a structural basis for the development of potential therapeutic agents against the fatal neurodegenerative diseases ALS/FTD.

Comprehensive analysis of intramolecular G-quadruplex structures: furthering the understanding of their formalism

G-quadruplexes (G4) are helical structures found in guanine-rich DNA or RNA sequences. Generally, their formalism is based on a few dozen structures, which can produce some inconsistencies or incompleteness. Using the website ASC-G4, we analyzed the structures of 333 intramolecular G4s, of all types, which allowed us to clarify some key concepts and present new information. To each of the eight distinguishable topologies corresponds a groove-width signature and a predominant glycosidic configuration (gc) pattern governed by the directions of the strands. The relative orientations of the stacking guanines within the strands, which we quantified and related to their vertical gc successions, determine the twist and tilt of the helices. The latter impact the minimum groove widths, which represent the space available for lateral ligand binding. The G4 four helices have similar twists, even when these twists are irregular, meaning that they have various angles along the strands. Despite its importance, the vertical gc succession has no strict one-to-one relationship with the topology, which explains the discrepancy between some topologies and their corresponding circular dichroism spectra. This study allowed us to introduce the new concept of platypus G4s, which are structures with properties corresponding to several topologies.

EndoQuad: a comprehensive genome-wide experimentally validated endogenous G-quadruplex database

G-quadruplexes (G4s) are non-canonical four-stranded structures and are emerging as novel genetic regulatory elements. However, a comprehensive genomic annotation of endogenous G4s (eG4s) and systematic characterization of their regulatory network are still lacking, posing major challenges for eG4 research. Here, we present EndoQuad (https://EndoQuad.chenzxlab.cn/) to address these pressing issues by integrating high-throughput experimental data. First, based on high-quality genome-wide eG4s mapping datasets (human: 1181; mouse: 24; chicken: 2) generated by G4 ChIP-seq/CUT&Tag, we generate a reference set of genome-wide eG4s. Our multi-omics analyses show that most eG4s are identified in one or a few cell types. The eG4s with higher occurrences across samples are more structurally stable, evolutionarily conserved, enriched in promoter regions, mark highly expressed genes and associate with complex regulatory programs, demonstrating higher confidence level for further experiments. Finally, we integrate millions of functional genomic variants and prioritize eG4s with regulatory functions in disease and cancer contexts. These efforts have culminated in the comprehensive and interactive database of experimentally validated DNA eG4s. As such, EndoQuad enables users to easily access, download and repurpose these data for their own research. EndoQuad will become a one-stop resource for eG4 research and lay the foundation for future functional studies.

The Nucleic Acid Knowledgebase: a new portal for 3D structural information about nucleic acids

The Nucleic Acid Knowledgebase (nakb.org) is a new data resource, updated weekly, for experimentally determined 3D structures containing DNA and/or RNA nucleic acid polymers and their biological assemblies. NAKB indexes nucleic acid-containing structures derived from all major structure determination methods (X-ray, NMR and EM), including all held by the Protein Data Bank (PDB). As the planned successor to the Nucleic Acid Database (NDB), NAKB's design preserves all functionality of the NDB and provides novel nucleic acid-centric content, including structural and functional annotations, as well as annotations from and links to external resources. A variety of custom interactive tools have been developed to enable rapid exploration and drill-down of NAKB's content.

EGFR suppression contributes to growth inhibitory activity of G-quadruplex ligands in non-small cell lung cancers

Non-small cell lung carcinomas (NSCLCs) commonly harbor activating mutations in the epidermal growth factor receptor (EGFR). Drugs targeting the tyrosine kinase activity of EGFR have shown effectiveness in inhibiting the growth of cancer cells with EGFR mutations. However, the development of additional mutations in cancer cells often leads to the persistence of the disease, necessitating alternative strategies to overcome this challenge. We explored the efficacy of stabilizing the G-quadruplex structure formed in the promoter region of EGFR as a means to suppress its expression and impede the growth of cancer cells with EGFR mutations. We revealed that the carbazole derivative BMVC-8C3O effectively suppressed EGFR expression and demonstrated significant growth inhibition in EGFR-mutated NSCLC cells, both in cell culture and mouse xenograft models. Importantly, the observed repression of EGFR expression and growth inhibition were not exclusive to carbazole derivatives, as several other G-quadruplex ligands exhibited similar effects. The growth-inhibitory activity of BMVC-8C3O is attributed, at least in part, to the repression of EGFR, although it is possible that additional cellular targets are also affected. Remarkably, the growth-inhibitory effect was observed even in osimertinib-resistant cells, indicating that BMVC-8C3O holds promise for treating drug-resistant NSCLC. Our findings present a promising and innovative approach for inhibiting the growth of NSCLC cells with EGFR mutations by effectively suppressing EGFR expression. The demonstrated efficacy of G-quadruplex ligands in this study highlights their potential as candidates for further development in NSCLC therapy.

G4Bank: A database of experimentally identified DNA G-quadruplex sequences

G-quadruplex (G4), a non-canonical nucleic acid structure, has been suggested to play a key role in important cellular processes including transcription, replication and cancer development. Recently, high-throughput sequencing approaches for G4 detection have provided a large amount of experimentally identified G4 data that reveal genome-wide G4 landscapes and enable the development of new methods for predicting potential G4s from sequences. Although several existing databases provide G4 experimental data and relevant biological information from different perspectives, there is no dedicated database to collect and analyze DNA G4 experimental data genome-widely. Here, we constructed G4Bank, a database of experimentally identified DNA G-quadruplex sequences. A total of 6,915,983 DNA G4s were collected from 13 organisms, and state-of-the-art prediction methods were performed to filter and analyze the G4 data. Therefore, G4Bank will facilitate users to access comprehensive G4 experimental data and enable sequence feature analysis of G4 for further investigation. The database of the experimentally identified DNA G-quadruplex sequences can be accessed at http://tubic.tju.edu.cn/g4bank/ .

Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding

DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.

G-quadruplex DNA: A Longer Story

G-quadruplexes (G4s) are distinctive four-stranded DNA or RNA structures found within cells that are thought to play functional roles in gene regulation and transcription, translation, recombination, and DNA damage/repair. While G4 structures can be uni-, bi-, or tetramolecular with respect to strands, folded unimolecular conformations are most significant in vivo. Unimolecular G4 can potentially form in sequences with runs of guanines interspersed with what will become loops in the folded structure: 5'G_xL_yG_xL_yG_xL_yG_x, where x is typically 2-4 and y is highly variable. Such sequences are highly conserved and specifically located in genomes. In the folded structure, guanines from each run combine to form planar tetrads with four hydrogen-bonded guanine bases; these tetrads stack on one another to produce four strand segments aligned in specific parallel or antiparallel orientations, connected by the loop sequences. Three types of loops (lateral, diagonal, or "propeller") have been identified. The stacked tetrads form a central cavity that features strong coordination sites for monovalent cations that stabilize the G4 structure, with potassium or sodium preferred. A single monomeric G4 typically forms from a sequence containing roughly 20-30 nucleotides. Such short sequences have been the primary focus of X-ray crystallographic or NMR studies that have produced high-resolution structures of a variety of monomeric G4 conformations. These structures are often used as the basis for drug design efforts to modulate G4 function.We believe that the focus on monomeric G4 structures formed by such short sequences is perhaps myopic. Such short sequences for structural studies are often arbitrarily selected and removed from their native genomic sequence context, and then are often changed from their native sequences by base substitutions or deletions intended to optimize the formation of a homogeneous G4 conformation. We believe instead that G-quadruplexes prefer company and that in a longer natural sequence context multiple adjacent G4 units can form to combine into more complex multimeric G4 structures with richer topographies than simple monomeric forms. Bioinformatic searches of the human genome show that longer sequences with the potential for forming multiple G4 units are common. Telomeric DNA, for example, has a single-stranded overhang of hundreds of nucleotides with the requisite repetitive sequence with the potential for formation of multiple G4s. Numerous extended promoter sequences have similar potentials for multimeric G4 formation. X-ray crystallography and NMR methods are challenged by these longer sequences (>30 nt), so other tools are needed to explore the possible multimeric G4 landscape. We have implemented an integrated structural biology approach to address this challenge. This approach integrates experimental biophysical results with atomic-level molecular modeling and molecular dynamics simulations that provide quantitatively testable model structures. In every long sequence we have studied so far, we found that multimeric G4 structures readily form, with a surprising diversity of structures dependent on the exact native sequence used. In some cases, stable hairpin duplexes form along with G4 units to provide an even richer landscape. This Account provides an overview of our approach and recent progress and provides a new perspective on the G-quadruplex folding landscape.

Insights into the Small Molecule Targeting of Biologically Relevant G-Quadruplexes: An Overview of NMR and Crystal Structures

G-quadruplexes turned out to be important targets for the development of novel targeted anticancer/antiviral therapies. More than 3000 G-quadruplex small-molecule ligands have been described, with most of them exerting anticancer/antiviral activity by inducing telomeric damage and/or altering oncogene or viral gene expression in cancer cells and viruses, respectively. For some ligands, in-depth NMR and/or crystallographic studies were performed, providing detailed knowledge on their interactions with diverse G-quadruplex targets. Here, the PDB-deposited NMR and crystal structures of the complexes between telomeric, oncogenic or viral G-quadruplexes and small-molecule ligands, of both organic and metal-organic nature, have been summarized and described based on the G-quadruplex target, from telomeric DNA and RNA G-quadruplexes to DNA oncogenic G-quadruplexes, and finally to RNA viral G-quadruplexes. An overview of the structural details of these complexes is here provided to guide the design of novel ligands targeting more efficiently and selectively cancer- and virus-related G-quadruplex structures.

QUADRatlas: the RNA G-quadruplex and RG4-binding proteins database

RNA G-quadruplexes (RG4s) are non-canonical, disease-associated post-transcriptional regulators of gene expression whose functions are driven by RNA-binding proteins (RBPs). Being able to explore transcriptome-wide RG4 formation and interaction with RBPs is thus paramount to understanding how they are regulated and exploiting them as potential therapeutic targets. Towards this goal, we present QUADRatlas (https://rg4db.cibio.unitn.it), a database of experimentally-derived and computationally predicted RG4s in the human transcriptome, enriched with biological function and disease associations. As RBPs are key to their function, we mined known interactions of RG4s with such proteins, complemented with an extensive RBP binding sites dataset. Users can thus intersect RG4s with their potential regulators and effectors, enabling the formulation of novel hypotheses on RG4 regulation, function and pathogenicity. To support this capability, we provide analysis tools for predicting whether an RBP can bind RG4s, RG4 enrichment in a gene set, and de novo RG4 prediction. Genome-browser and table views allow exploring, filtering, and downloading the data quickly for individual genes and in batch. QUADRatlas is a significant step forward in our ability to understand the biology of RG4s, offering unmatched data content and enabling the integrated analysis of RG4s and their interactions with RBPs.

RNAloops: a database of RNA multiloops

MOTIVATION

Knowledge of the 3D structure of RNA supports discovering its functions and is crucial for designing drugs and modern therapeutic solutions. Thus, much attention is devoted to experimental determination and computational prediction targeting the global fold of RNA and its local substructures. The latter include multi-branched loops-functionally significant elements that highly affect the spatial shape of the entire molecule. Unfortunately, their computational modeling constitutes a weak point of structural bioinformatics. A remedy for this is in collecting these motifs and analyzing their features.

RESULTS

RNAloops is a self-updating database that stores multi-branched loops identified in the PDB-deposited RNA structures. A description of each loop includes angular data-planar and Euler angles computed between pairs of adjacent helices to allow studying their mutual arrangement in space. The system enables search and analysis of multiloops, presents their structure details numerically and visually, and computes data statistics.

AVAILABILITY AND IMPLEMENTATION

RNAloops is freely accessible at https://rnaloops.cs.put.poznan.pl.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.