- Open Access
ElTetrado: a tool for identification and classification of tetrads and quadruplexes
BMC Bioinformatics volume 21, Article number: 40 (2020)
Quadruplexes are specific structure motifs occurring, e.g., in telomeres and transcriptional regulatory regions. Recent discoveries confirmed their importance in biomedicine and led to an intensified examination of their properties. So far, the study of these motifs has focused mainly on the sequence and the tertiary structure, and concerned canonical structures only. Whereas, more and more non-canonical quadruplex motifs are being discovered.
Here, we present ElTetrado, a software that identifies quadruplexes (composed of guanine- and other nucleobase-containing tetrads) in nucleic acid structures and classifies them according to the recently introduced ONZ taxonomy. The categorization is based on the secondary structure topology of quadruplexes and their component tetrads. It supports the analysis of canonical and non-canonical motifs. Besides the class recognition, ElTetrado prepares a dot-bracket and graphical representations of the secondary structure, which reflect the specificity of the quadruplex’s structure topology. It is implemented as a freely available, standalone application, available at https://github.com/tzok/eltetrado.
The proposed software tool allows to identify and classify tetrads and quadruplexes based on the topology of their secondary structures. It complements existing approaches focusing on the sequence and 3D structure.
One of the key problems that can be solved by structural bioinformatics is the discovery and analysis of motifs that form in molecular structures. Recurring patterns are searched for at all levels of the architecture: in sequence, secondary and tertiary structure [1–8]. Motif’s discovery entails the research into its influence on the molecule functions. Among the motifs that are currently being researched in many areas of life sciences are quadruplexes. These structures were first found in DNA, in the late 1960s. However, only after their association with targeted therapies against cancer and neurodegenerative diseases, the study on these motifs started to expand . Today, computationally-supported analysis of quadruplex characteristics constitutes a relatively new challenge in structural bioinformatics and computational biology [10, 11].
The building unit of the quadruplex is a tetrad (called also a quartet). Single tetrad is made up of four nucleotide residues arranged in a plane. Each of them forms two non-canonical bonds (one along the Watson-Crick edge, and one along the Hoogsteen edge) with the others. If at least two tetrads are stacked upon one another at a distance of about 3.3 Å, the quadruplex is created . Its structure is considered at the sequence level, examining the N-tracts which are part of the motif. The 3D structure analysis encompasses the properties of tetrads, stems and loops that connect external quartets. The structural diversity of quadruplexes is large and new motifs are constantly being discovered [13–18]. Meanwhile, most studies have been limited to their specific subsets. In particular, many researchers focused on canonical quadruplexes built of Guanines (so-called G4s). The latter ones were the subject of analysis and formalism proposed in .
Recently, we have proposed a new approach to the analysis and classification of tetrads and quadruplexes . Its foundation is the secondary structure topology of these motifs. Our concept has resulted from the observation of patterns (O-shaped, N-shaped, and Z-shaped tracks) occurring in the secondary structure diagrams of quadruplexes generated and annotated by RNApdbee 2.0, a web server from the RNApolis toolset [21, 22]. We have defined six classes of tetrad topologies (so called ONZ taxonomy) that correspond to pairings between quartet-forming nucleotides. Next, we have introduced ONZ-based formalism for quadruplexes . It addresses both canonical and non-canonical motifs. Finally, we have designed two-line dot-bracket that allows for unambiguous encoding of the secondary structure of tetrads and quadruplexes, and we have adjusted arc diagrams to represent a variety of their topologies.
Here, we present ElTetrado, a software tool that accompanies our novel classification approach, and the results of its application in processing nucleic acids structures from Protein Data Bank. ElTetrado identifies tetrads and quadruplexes in nucleic acid structures and assigns them to categories according to ONZ taxonomy. It also provides dot-bracket encoded secondary structures of these motifs and the corresponding arc diagrams. The application is freely available at https://github.com/tzok/eltetrado.
In the preliminary computational experiments with ElTetrado, we aimed to find the distribution of tetrads and quadruplexes in different ONZ categories . We focused on unimolecular motifs, as the classification is based on the order of nucleotides in the chain, which is unambiguously determined only in the case of single-stranded structures. The input dataset was constructed from PDB-deposited nucleic acid structures on 18 April, 2019. The research showed that in unimolecular structures O-type tetrads constituted the majority (75%). We have assumed that a similar situation should occur in bi- and tetramolecular motifs.
This assumption has become the starting point for the development of a procedure for the classification of tetrads and quadruplexes made up of more than one strand. The procedure optimizes the processing order of the motif’s strands so as to maximize the allocation of tetrads to the O category. The result of its application is the assignment of the ONZ categories to bi- and tetramolecular tetrads and quadruplexes as well. Strand reordering can be switched off by calling ElTetrado with --no-reorder parameter. In such a program call, bi- and tetramolecular motifs are not classified according to the ONZ taxonomy.
An example outcome of ElTetrado running with and without the reordering procedure is shown in Fig. 1. The analysed structure of DNA/RNA hybrid (PDB ID: 1N7A)  is composed of eight strands which form ten tetrads. ElTetrado called with --no-reorder parameter is not able to classify the tetrads and generates the arc diagram as shown in Fig. 1a. If the parameter is not passed, the reordering procedure is executed. In this case, ElTetrado identifies two regular 5-tetrad quadruplexes assigned to Op category (Fig. 1b). The first one includes four O+ type tetrads (dark blue) and one O- type (light blue) quartet. The second quadruplex has three O+ tetrads and two tetrads in O- group.
Using the reordering procedure, we executed ElTetrado for a set of tetrads and quadruplexes identified in nucleic acid structures deposited in the Protein Data Bank as of 18 December, 2019 (Tables 1, 2, and 3). In this experiment, we were particularly interested in the coverage of various categories in the ONZ taxonomy by bi- and tetramolecular motifs. The analysed dataset contained 328 structures, 148 of which were determined by X-ray crystallography (107 DNAs, 38 RNAs and 3 modified nucleic acids) and 180 by NMR (169 DNAs, 9 RNAs, 2 modified). All X-ray structures in this collection have good resolution. Its minimum value is 0.61 Å, maximum 3.793 Å, and average resolution is 2.148 Å. In the considered structures, ElTetrado identified 1445 tetrads and 452 quadruplexes in total. The analysis showed that most tetrads (namely 1035) in the collection were in O-type classes (O+ type tetrads accounted for 86% of this group, while O- tetrads for 14%). Moreover, ElTetrado found 197 tetrads from N categories (65% of them were in N+, 35% in N-), and 213 from Z (half in half in classes Z+ and Z-). Among quadruplexes, we found 424 regular motifs and 28 irregular ones. The subset of regular quadruplexes contained 312 cases from class O (72% of them were in Op, 16% in Oa, and 12% in Oh), 58 from class N (88% of them were in Na, 12% in Nh), and 54 from class Z (98% of them in Za, 2% in Zh). Among irregular quadruplexes, 46% of them were in Mp and 54% in Mh category. Let us note that among bi- and tetramolecular quadruplexes, none was in negative (-) subcategory. These quadruplexes belonged to positive (+) or hybrid (*) subtypes (Table 3).
Recent study of the topology of quadruplex secondary structure bore fruit with new approach to analyse and categorize these motifs . ElTetrado is the first software tool that - apart from identification of tetrads and quadruplexes - can assign them to classes defined in ONZ taxonomy. It can also prepare their unequivocal representation in dot-bracket notation and arc diagram. We believe that by providing an insight into the secondary structure features, ElTetrado complements formalisms proposed for canonical G4s . Thus, our method will contribute to a better understanding of quadruplex architecture and improve its full description.
ElTetrado processes PDB and mmCIF files to identify quadruplexes and their component tetrads in nucleic acid structures (Fig. 2). It applies DSSR  to collect the preliminary information about base pairs and stacking. Next, it scans the properties of intra-tetrad interactions to determine tetrad categories within ONZ taxonomy. Once the tetrads are classified, the algorithm identifies quadruplexes. These motifs are next categorized based on the classes of their component quartets and the order of nucleotides in the N-tract. Finally, ElTetrado prepares the dot-bracket representation of the secondary structure of identified quadruplexes. It also applies R4RNA package  to create the output arc diagram.
Tetrad classification procedure
There are six classes of tetrad topologies in ONZ taxonomy: O+, O-, N+, N-, Z+ and Z- . They differ in hydrogen-bonding interactions, i.e., in each class, the other nucleotide residues form base pairs, along the other edges - Watson-Crick, Hoogsteen or Sugar edge.
Let us have four residues building tetrad T. We will mark them with letters A, B, C, D. Suppose that in the sequence, A is the first from 5’-end of these four residues. Assume also that the following base pairs are formed in T: (A,B), (B,C), (C,D), (D,A), where the first residue of the pair binds with the second one along its Watson-Crick edge, while the second with the first along its Hoogsteen or Sugar edge (for example, in pair (A,B), A binds along its Watson-Crick edge, B binds along its Hoogsteen edge). Let us now check in what order the residues are arranged in the nucleotide sequence from 5’ to 3’ end. Depending on this arrangement, tetrad T is assigned to a given category:
if 5’ - A - B - C - D - 3’ then T ∈ O+
if 5’ - A - D - C - B - 3’ then T ∈ O-
if 5’ - A - B - D - C - 3’ then T ∈ N+
if 5’ - A - C - D - B - 3’ then T ∈ N-
if 5’ - A - C - B - D - 3’ then T ∈ Z+
if 5’ - A - D - B - C - 3’ then T ∈ Z-.
A pseudocode of tetrad classification procedure applied in ElTetrado is given in Algorithm 1.
Quadruplex classification procedure
In the case of quadruplexes, we have distinguished between regular and irregular structures. A regular motif contains tetrads of one general type, i.e. O, N or Z. We pay attention to the order of nucleotides in N-tracts, thus having parallel (p), antiparallel (a) or hybrid (h) motifs (cf. Algorithm 2). This way we define nine classes of regular quadruplexes: Op, Oa, Oh, Np, Na, Nh, Zp, Za, Zh.
Let us have a quadruplex Q composed of three tetrads, T i(A i-B i-C i-D i), T j(A j-B j-C j-D j) and T k(A k-B k-C k-D k), all of them in the O+ category. Quadruplex Q is assigned to a given category if the following order of residues is preserved in the nucleic acid strand:
if 5’ - A iAjAk - B iBjBk - C iCjCk - D iDjDk - 3’ then Q ∈ Op
if two tracts have different order than the other two (for example, 5’ - A iAjAk - B kBjBi - C iCjCk - D kDjDi - 3’) then Q ∈ Oa
if one tract has a different order than the other three (for example, 5’ - A iAjAk - B kBjBi - C kCjCi - D kDjDi - 3’) then Q ∈ Oh.
If a quadruplex is built from tetrads of different types, it is assigned to category M (mixed) that collects irregular motifs. Regular and irregular categories of quadruplexes can be further divided into positive (+), negative (-) and hybrid (*) subcategories, depending on the + and - signs in the classification of quadruplex-building tetrads. For example, if a parallel quadruplex contains only O+ tetrads it is assigned to Op+ group, if it is built of O- tetrads only we assign it to class Op-, and if it contains a mixture of O+ and O- tetrads then it belongs to Op* subcategory.
For the purposes of quadruplexes with unclassified tetrads we have also introduced the R class (remaining ones) . However, recent computational experiments have not shown the existence of motifs that could belong to this class.
Algorithm 3 presents the quadruplex classification procedure executed in ElTetrado.
ElTetrado is a standalone application implemented in Python 3, using its standard libraries. Therefore it works under all common operating systems. In order to start working with the program, the users need to download it onto the local drive from https://github.com/tzok/eltetrado.
We recommend that, apart from ElTetrado, the users should download the DSSR binary  and place it in the same local directory. DSSR is utilized for the preliminary analysis of base pairs in the input 3D structure. Its local installation allows the users to control DSSR execution. For example, one can decide to pass --symmetry parameter to x3dna-dssr binary when dealing with X-ray structures, which is necessary for some quadruplexes. Alternatively, if the users do not want to have a local version of DSSR binary, they can obtain the DSSR output in JSON format from any place and use them as input data for ElTetrado (with --dssr-json parameter).
ElTetrado is started from the command line. The users enter the program name and either --pdb followed by an input file name (the file should be in PDB or mmCIF format), or --dssr-json followed by a path to JSON file generated by DSSR, or both switches at once. Optionally, users can pass the following input parameter(s) that allow controlling ElTetrado operation mode: --stacking-mismatchThis parameter can take one of three values: 0, 1, 2 (default). It controls the algorithm’s tolerance of mismatches when analyzing tetrad stacking. If stacking-mismatch equals 0, then the ideal geometry is enforced. It means that exactly all 4 pairs of nucleotides should be stacked for two tetrads to be considered the quadruplex or N4-helix. If stacking-mismatch equals 2 (default value), then the room is left for some flexibility and imperfect stacking. --relaxed-stem-definitionIf this parameter is passed, ElTetrado treats sequentially neighbouring tetrads as part of the same quadruplex, regardless of their mutual stacking interactions. --strictIf this parameter is set, ElTetrado looks for cis Watson-Crick–Hoogsteen base pairs only when looking for tetrads. By default, this mode is switched off what allows ElTetrado to find tetrads linked by base pairs of different types. --no-reorderThe ONZ classification is based on strand direction from 5’ to 3’. When multiple strands contribute to build the quadruplex, their order influences the assigned class. Therefore, ElTetrado checks all permutations of chain orders in bi- and tetramolecular quadruplexes before classifying them. If this parameter is set, chain order remains the same as in the input file and bi- and tetramolecular quadruplexes are not classified. --complete-2dRunning ElTetrado with this parameter causes that in the arc diagram, apart from the quadruplex, canonical base pairs occurring in the remaining part of the structure will be represented. This option allows seeing the quadruplex in the context of the entire structure.
The arcs in the arc diagram are assigned colors depending on the category of represented tetrads. By default, class O is blue, class N is green, and class Z is orange. Positive (+) and negative (-) subcategories are distinguished by light and dark shades of the same color, respectively. For example, O+ is dark blue, O- is light blue in the diagram. Unclassified tetrads are colored grey. Canonical base pairs that do not belong to the quadruplex are represented by black arcs. However, the users can define in the configuration file any set of colors for annotating ONZ classes in the diagram.
Output data includes the following information about quadruplexes identified in the input structure: (i) quadruplex category in ONZ-based taxonomy, (ii) the list of tetrads that compose the motif with their ONZ classification, (iii) quadruplex secondary structure in dot-bracket notation, (iv) arc diagram representing quadruplex secondary structure. Example output data for the input PDB structure of dimeric RNA quadruplex (PDB id: 1MY9)  is presented on the bottom panel of Fig. 2.
Availability and requirements
Project name: ELTetrado Project home page:https://github.com/tzok/eltetradoOperating system(s): Platform independent Programming language: Python 3 Other requirements: DSSR binary License: MIT Any restrictions to use by non-academics: no restrictions
Availability of data and materials
All the data and materials are available at the project home page.
Macromolecular crystallographic information file
Nuclear magnetic resonance
O-shaped, N-shaped, Z-shaped
Protein data bank
Leontis NB, Westhof E. Geometric nomenclature and classification of RNA base pairs. RNA. 2001; 7(4):499–512. https://doi.org/10.1093/nar/gks339.
Jonassen I, Eidhammer I, Conklin D, Taylor WR. Structure motif discovery and mining the PDB. Bioinformatics. 2002; 18(2):362–7. https://doi.org/10.1093/bioinformatics/18.2.362.
Popenda M, Blazewicz M, Szachniuk M, Adamiak RW. RNA FRABASE version 1.0: an engine with a database to search for the three-dimensional fragments within RNA structures. Nucleic Acids Res. 2008; 36(D1):386–91. https://doi.org/10.1093/nar/gkm786.
Djelloul M, Denise A. Automated motif extraction and classification in RNA tertiary structures. RNA. 2008; 14(12):2489–97. https://doi.org/10.1261/rna.1061108.
Petrov AI, Zirbel CL, Leontis NB. Automated classification of RNA 3D motifs and the RNA 3D Motif Atlas. RNA. 2013; 19(10):1327–40. https://doi.org/10.1261/rna.039438.113.
Chojnowski G, Waleń T, Bujnicki JM. RNA bricks - a database of RNA 3D motifs and their interactions. Nucleic Acids Res. 2013; 42(D1):123–31. https://doi.org/10.1093/nar/gkt1084.
Miskiewicz J, Szachniuk M. Discovering structural motifs in miRNA precursors from the Viridiplantae kingdom. Molecules. 2018; 23(6):1367. https://doi.org/10.3390/molecules23061367.
Reinharz V, Soulé A, Westhof E, Waldispühl J, Denise A. Mining for recurrent long-range interactions in RNA structures reveals embedded hierarchies in network families. Nucleic Acids Res. 2018; 46(8):3841–51. https://doi.org/10.1093/nar/gky197.
Cammas A, Millevoi S. RNA G-quadruplexes: emerging mechanisms in disease. Nucleic Acids Res. 2016; 45(4):1584–95. https://doi.org/10.1093/nar/gkw1280.
Gudanis D, Popenda L, Szpotkowski K, Kierzek R, Gdaniec Z. Structural characterization of a dimer of RNA duplexes composed of 8-bromoguanosine modified CGG trinucleotide repeats: a novel architecture of RNA quadruplexes. Nucleic Acids Res. 2016; 44(5):2409–16. https://doi.org/10.1093/nar/gkv1534.
O’Hagan MP, Mergny J-L, Waller ZAE. G-quadruplexes in prague: A bohemian rhapsody. Biochimie. 2018; 147:170–80. https://doi.org/10.1016/j.biochi.2018.02.004.
Fay MM, Lyons SM, Ivanov P. RNA G-quadruplexes in biology: principles and molecular mechanisms. J Mol Biol. 2017; 429(14):2127–47. https://doi.org/10.1016/j.jmb.2017.05.017.
Huppert JL. Prevalence of quadruplexes in the human genome. Nucleic Acids Res. 2005; 33(9):2908–16. https://doi.org/10.1093/nar/gki609.
Todd AK. Highly prevalent putative quadruplex sequence motifs in human DNA. Nucleic Acids Res. 2005; 33(9):2901–7. https://doi.org/10.1093/nar/gki553.
Scaria V, Hariharan M, Arora A, Maiti S. Quadfinder: server for identification and analysis of quadruplex-forming motifs in nucleotide sequences. Nucleic Acids Res. 2006; 34(Web Server):683–5. https://doi.org/10.1093/nar/gkl299.
Xu Y, Suzuki Y, Komiyama M. Click chemistry for the identification of G-quadruplex structures: Discovery of a DNA-RNA G-quadruplex. Angewandte Chemie International Edition. 2009; 48(18):3281–4. https://doi.org/10.1002/anie.200806306.
González V, Guo K, Hurley L, Sun D. Identification and characterization of nucleolin as a c-mycG-quadruplex-binding protein. J Biol Chem. 2009; 284(35):23622–35. https://doi.org/10.1074/jbc.m109.018028.
Lyons SM, Gudanis D, Coyne SM, Gdaniec Z, Ivanov P. Identification of functional tetramolecular RNA G-quadruplexes derived from transfer RNAs. Nat Commun. 2017; 8(1). https://doi.org/10.1038/s41467-017-01278-w.
Webba da Silva M. Geometric formalism for DNA quadruplex folding. Chem Eur J. 2007; 13(35):9738–45. https://doi.org/10.1002/chem.200701255.
Popenda M, Miskiewicz J, Sarzynska J, Zok T, Szachniuk M. Topology-based classification of tetrads and quadruplex structures. Bioinformatics. 2019. https://doi.org/10.1093/bioinformatics/btz738.
Zok T, Antczak M, Zurkowski M, Popenda M, Blazewicz J, Adamiak RW, Szachniuk M. RNApdbee 2.0: multifunctional tool for RNA structure annotation. Nucleic Acids Res. 2018; 46(W1):30–5. https://doi.org/10.1093/nar/gky314.
Szachniuk M. RNApolis: computational platform for RNA structure analysis. Found Comput Dec Sci. 2019; 44(2):241–57. https://doi.org/10.2478/fcds-2019-0012.
Ravelli RBG, Leiros H. -K. S., Pan B, Caffrey M, McSweeney S. Specific radiation damage can be used to solve macromolecular crystal structures. Structure. 2003; 11(2):217–24. https://doi.org/10.1016/S0969-2126(03)00006-6.
Lu X. -J., Bussemaker HJ, Olson WK. DSSR: an integrated software tool for dissecting the spatial structure of RNA. Nucleic Acids Res. 2015; 43(21):142. https://doi.org/10.1093/nar/gkv716.
Lai D, Proctor JR, Zhu JYA, Meyer IM. R-CHIE: a web server and R package for visualizing RNA secondary structures. Nucleic Acids Res. 2012; 40(12):95. https://doi.org/10.1093/nar/gks241.
Liu H, Matsugami A, Katahira M, Uesugi S. A dimeric RNA quadruplex architecture comprised of two G:G(:A):G:G(:A) hexads, G:G:G:G tetrads and UUUU loops. J Mol Biol. 2002; 322(5):955–70. https://doi.org/10.1016/s0022-2836(02)00876-8.
This research was carried in the European Centre for Bioinformatics and Genomics (Poznan, Poland).
This research was supported by the National Science Centre, Poland [2016/23/B/ST6/03931]. The authors also acknowledge partial support by the statutory funds of Poznan University of Technology, Polish Ministry of Science and Higher Education, and the Institute of Bioorganic Chemistry, PAS within intramural financing program. The funding bodies did not play any role in the design of the study, collection, analysis and interpretation of data, and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zok, T., Popenda, M. & Szachniuk, M. ElTetrado: a tool for identification and classification of tetrads and quadruplexes. BMC Bioinformatics 21, 40 (2020). https://doi.org/10.1186/s12859-020-3385-1
- Secondary structure topology
- Dot-bracket notation
- Arc diagram