The evolution of protein-coding genes has both structural and regulatory components. The first can be assessed by measuring the ratio of non-synonymous to synonymous nucleotide substitutions. The second component can be measured as the normalized proportion of transposable elements that are used as regulatory elements. For the first time, we characterized in parallel the regulatory and structural evolutionary profiles for 10,890 human genes and 2972 molecular pathways. We observed a ~0.1 correlation between the structural and regulatory metrics at the gene level, which appeared much higher (~0.4) at the pathway level. We deposited the data in the publicly available database RetroSpect. We also analyzed the evolutionary dynamics of six cancer pathways of two major axes: Notch/WNT/Hedgehog and AKT/mTOR/EGFR. The Hedgehog pathway had both components slower, whereas the Akt pathway had clearly accelerated structural evolution. In particular, the major hub nodes Akt and beta-catenin showed both components strongly decreased, whereas two major regulators of Akt TCL1 and CTMP had outstandingly high evolutionary rates. We also noticed structural conservation of serine/threonine kinases and the genes related to guanosine metabolism in cancer signaling: GPCRs, G proteins, and small regulatory GTPases (Src, Rac, Ras); however, this was compensated by the accelerated regulatory evolution.

RNA polymerase III (Pol III) transcribes tRNA genes using type II promoters. The internal control regions contain a Box A and a Box B, which are recognized by TFIIIC. The 5′-flanking regions of tRNA genes clearly play a role in the regulation of transcription, but consensus sequences in it have been found only in some plants and S. pombe; although, the TATA binding protein (TBP) is a component of the TFIIIB complex in all eukaryotes. Archaea utilize an ortholog of the TBP. The goal of this work is the detection of the positions of intragenic and extragenic promoters of Pol III, which regulate the transcription of tRNA genes in eukaryotes and archaea. For this purpose, we analyzed textual and some structural, mechanical, and physicochemical properties of the DNA in the 5′-flanking regions of tRNA genes, as well as in 30 bp at the beginning of genes and 60 bp at the end of genes in organisms possessing the TBP or its analog (eukaryotes, archaea) and organisms not possessing the TBP (bacteria). Representative tRNA gene sets of 11 organisms were taken from the GtRNAdb database. We found that the consensuses of A- and B-boxes in organisms from all three domains are identical; although, they differ in the conservativism of some positions. Their location relative to the ends of tRNA genes is also identical. In contrast, the structural and mechanical properties of DNA in the 5′-flanking regions of tRNA genes differ not only between organisms from different domains, but also between organisms from the same domain. Well-expressed TBP binding positions are found only in S. pombe and A. thaliana. We discuss possible reasons for the variability of the 5′-flanking regions of tRNA genes.

RNA polymerase II (POL II) is responsible for the transcription of messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs). Previously, we have shown the evolutionary invariance of the structural features of DNA in the POL II core promoters of the precursors of mRNAs. In this work, we have analyzed the POL II core promoters of the precursors of lncRNAs in Homo sapiens and Mus musculus genomes. Structural analysis of nucleotide sequences in positions −50, +30 bp in relation to the TSS have shown the extremely heterogeneous 3D structure that includes two singular regions - hexanucleotide “INR” around the TSS and octanucleotide “TATA-box” at around ~−28 bp upstream. Thus, the 3D structure of core promoters of lncRNA resembles the architecture of the core promoters of mRNAs; however, textual analysis revealed differences between promoters of lncRNAs and promoters of mRNAs, which lies in their textual characteristics; namely, the informational entropy at each position of the nucleotide text of lncRNA core promoters (by the exception of singular regions) is significantly higher than that of the mRNA core promoters. Another distinguishing feature of lncRNA is the extremely rare occurrence in the TATA box of octanucleotides with the consensus sequence. These textual differences can significantly affect the efficiency of the transcription of lncRNAs.

Nucleotide motifs in the core promoters of eukaryotic protein-coding genes transcribed by RNA polymerase II (Pol II) play an important role in the transcription process. We analyzed the role of an octanucleotide located in the TATA box position. Depending on whether this octanucleotide can form a complex with the TATA-binding protein (TBP), the promoter is classified as either TATA-containing or TATA-less. We analyzed the differences in the primary and spatial structures, as well as their dynamics, in TATA-containing and TATA-less promoters of mammals and plants. We divided the complete promoter sets of six organisms (H. sapiens, M. musculus, C. familiaris, A. thaliana, Z. mays, and H. vulgare) from the EPDnew database into TATA-containing and TATA-less fractions. The sizes of the TATA-containing promoter fractions are significantly smaller than those of the TATA-less fractions in all studied organisms, except in A. thaliana, where the sizes of both fractions are approximately equal. We characterized promoter architecture using variation profiles of various base-pair step parameters, minor-groove width, and the conformational dynamics of native DNA. The architectures of TATA-containing and TATA-less promoters differ significantly. The possible mechanistic influence of DNA structural features on the formation of the pre-initiation complex (PIC) in both types of promoters is discussed.