G4-iM Grinder: when size and frequency matter. G-Quadruplex, i-Motif and higher order structure search and analysis tool

The Rise of Bacterial G-Quadruplexes in Current Antimicrobial Discovery

Antimicrobial resistance (AMR) is a silent critical issue that poses several challenges to health systems. While the discovery of novel antibiotics is currently stalled and prevalently focused on chemical variations of the scaffolds of available drugs, novel targets and innovative strategies are urgently needed to face this global threat. In this context, bacterial G-quadruplexes (G4s) are emerging as timely and profitable targets for the design and development of antimicrobial agents. Indeed, they are expressed in regulatory regions of bacterial genomes, and their modulation has been observed to provide antimicrobial effects with translational perspectives in the context of AMR. In this work, we review the current knowledge of bacterial G4s as well as their modulation by small molecules, including tools and techniques suitable for these investigations. Finally, we critically analyze the needs and future directions in the field, with a focus on the development of small molecules as bacterial G4s modulators endowed with remarkable drug-likeness.

iM-Seeker: a webserver for DNA i-motifs prediction and scoring via automated machine learning

DNA, beyond its canonical B-form double helix, adopts various alternative conformations, among which the i-motif, emerging in cytosine-rich sequences under acidic conditions, holds significant biological implications in transcription modulation and telomere biology. Despite recognizing the crucial role of i-motifs, predictive software for i-motif forming sequences has been limited. Addressing this gap, we introduce 'iM-Seeker', an innovative computational platform designed for the prediction and evaluation of i-motifs. iM-Seeker exhibits the capability to identify potential i-motifs within DNA segments or entire genomes, calculating stability scores for each predicted i-motif based on parameters such as the cytosine tracts number, loop lengths, and sequence composition. Furthermore, the webserver leverages automated machine learning (AutoML) to effortlessly fine-tune the optimal i-motif scoring model, incorporating user-supplied experimental data and customised features. As an advanced, versatile approach, 'iM-Seeker' promises to advance genomic research, highlighting the potential of i-motifs in cell biology and therapeutic applications. The webserver is freely available at https://im-seeker.org.

Site-specific incorporation of a fluorescent nucleobase analog enhances i-motif stability and allows monitoring of i-motif folding inside cells

The i-motif is an intriguing non-canonical DNA structure, whose role in the cell is still controversial. Development of methods to study i-motif formation under physiological conditions in living cells is necessary to study its potential biological functions. The cytosine analog 1,3-diaza-2-oxophenoxazine (tCO) is a fluorescent nucleobase able to form either hemiprotonated base pairs with cytosine residues, or neutral base pairs with guanines. We show here that when tCO is incorporated in the proximity of a G:C:G:C minor groove tetrad, it induces a strong thermal and pH stabilization, resulting in i-motifs with Tm of 39ºC at neutral pH. The structural determination by NMR methods reveals that the enhanced stability is due to a large stacking interaction between the guanines of the tetrad with the tCO nucleobase, which forms a tCO:C+ in the folded structure at unusually-high pHs, leading to an increased quenching in its fluorescence at neutral conditions. This quenching is much lower when tCO is base-paired to guanines and totally disappears when the oligonucleotide is unfolded. By taking profit of this property, we have been able to monitor i-motif folding in cells.

Basic protein- and peptide-induced stabilization of long-loop DNA G-guadruplexes

The human genome contains many G-quadruplex-forming sequences, including sequences containing long single-stranded loops that are believed to be unfavorable for G-quadruplex formation. The intracellular environment of biological cells is crowded with proteins with charged surfaces. Understanding the effects of protein-rich environments is important for understanding the formation of G-quadruplexes in an intracellular environment. In this study, we investigated the structural stability of DNA G-quadruplexes in the presence of several types of globular proteins (lysozyme, cytochrome c, bovine serum albumin, myoglobin, histone proteins, and serum proteins), unstructured polypeptides (protamine and poly-l-lysine), and oligopeptides (RGG/RG-domain peptides and short repeated peptides). Thermal melting studies of G-quadruplex-forming oligonucleotides derived from the human telomeric repeat sequence revealed that environments containing high concentrations of proteins and peptides differently affected the G-quadruplex stability according to their loop lengths. We found that weak electrostatic interactions of G-quadruplex loops with basic proteins and peptides improved the stability of long-loop G-quadruplexes and the interactions were strengthened under crowded conditions simulated by dextran. The comparison of the effects of different types of proteins and peptides indicated that excluded volume interactions and structural flexibility of both DNA and polypeptide chains influenced the efficiency of their interactions. This study provides insights into long-loop G-quadruplex stability in a crowded intracellular environment and the recognition of G-quadruplexes by arginine-rich domains of G-quadruplex-binding proteins.

Prediction of DNA i-motifs via machine learning

i-Motifs (iMs), are secondary structures formed in cytosine-rich DNA sequences and are involved in multiple functions in the genome. Although putative iM forming sequences are widely distributed in the human genome, the folding status and strength of putative iMs vary dramatically. Much previous research on iM has focused on assessing the iM folding properties using biophysical experiments. However, there are no dedicated computational tools for predicting the folding status and strength of iM structures. Here, we introduce a machine learning pipeline, iM-Seeker, to predict both folding status and structural stability of DNA iMs. The programme iM-Seeker incorporates a Balanced Random Forest classifier trained on genome-wide iMab antibody-based CUT&Tag sequencing data to predict the folding status and an Extreme Gradient Boosting regressor to estimate the folding strength according to both literature biophysical data and our in-house biophysical experiments. iM-Seeker predicts DNA iM folding status with a classification accuracy of 81% and estimates the folding strength with coefficient of determination (R2) of 0.642 on the test set. Model interpretation confirms that the nucleotide composition of the C-rich sequence significantly affects iM stability, with a positive correlation with sequences containing cytosine and thymine and a negative correlation with guanine and adenine.

Preferential formation of Z-RNA over intercalated motifs in long noncoding RNA

Secondary structure is a principal determinant of lncRNA function, predominantly regarding scaffold formation and interfaces with target molecules. Noncanonical secondary structures that form in nucleic acids have known roles in regulating gene expression and include G-quadruplexes (G4s), intercalated motifs (iMs), and R-loops (RLs). In this paper, we used the computational tools G4-iM Grinder and QmRLFS-finder to predict the formation of each of these structures throughout the lncRNA transcriptome in comparison to protein-coding transcripts. The importance of the predicted structures in lncRNAs in biological contexts was assessed by combining our results with publicly available lncRNA tissue expression data followed by pathway analysis. The formation of predicted G4 (pG4) and iM (piM) structures in select lncRNA sequences was confirmed in vitro using biophysical experiments under near-physiological conditions. We find that the majority of the tested pG4s form highly stable G4 structures, and identify many previously unreported G4s in biologically important lncRNAs. In contrast, none of the piM sequences are able to form iM structures, consistent with the idea that RNA is unable to form stable iMs. Unexpectedly, these C-rich sequences instead form Z-RNA structures, which have not been previously observed in regions containing cytosine repeats and represent an interesting and underexplored target for protein-RNA interactions. Our results highlight the prevalence and potential structure-associated functions of noncanonical secondary structures in lncRNAs, and show G4 and Z-RNA structure formation in many lncRNA sequences for the first time, furthering the understanding of the structure-function relationship in lncRNAs.

A beginner's handbook to identify and characterize i-motif DNA

Genomic DNA exhibits an innate ability to manifest diverse sequence-dependent secondary structures, serving crucial functions in gene regulation and cellular equilibrium. While extensive research has confirmed the formation of G-quadruplex structures by guanine-rich sequences in vitro and in cells, recent investigations have turned the quadruplex community's attention to the cytosine (C)-rich complementary strands that can adopt unique tetra-stranded conformation, termed as intercalated motif or i-motif. I-motifs are stabilized by hemi-protonated C:CH+ base pairs under acidic conditions. Initially, the in vivo occurrence of i-motifs was underestimated because their formation is favored at non-physiological pH. However, groundbreaking research utilizing the structure-specific iMab antibody and high-throughput sequencing have recently detected their conserved dispersion throughout the genome, challenging previous assumptions. Given the evolving nature of this research field, it becomes imperative to conduct independent in vitro experiments aimed at identifying potential i-motif formation in C-rich sequences and consolidating the findings to address the properties of i-motifs. This chapter serves as an introductory guide for the swift identification of novel i-motifs, where we present an experimental framework for investigating and characterizing i-motif sequences in vitro. In this chapter, we selected a synthetic oligonucleotide (C7T3) sequence and outlined appropriate methodologies for annealing the i-motif structure into suitable buffers. Then, we validated its formation by CD (Circular Dichroism) and NMR (Nuclear Magnetic Resonance) spectroscopy. Finally, we provided a thorough account of the step-by-step procedures to investigate the effect of i-motif formation on the stalling or retardation of DNA replication using high resolution primer extension assays.

G-Quadruplexes in the Viral Genome: Unlocking Targets for Therapeutic Interventions and Antiviral Strategies

G-quadruplexes (G4s) are unique non-canonical four-stranded nucleic acid secondary structures formed by guanine-rich DNA or RNA sequences. Sequences with the potential to form quadruplex motifs (pG4s) are prevalent throughout the genomes of all organisms, spanning from prokaryotes to eukaryotes, and are enriched within regions of biological significance. In the past few years, the identification of pG4s within most of the Baltimore group viruses has attracted increasing attention due to their occurrence in regulatory regions of the genome and the subsequent implications for regulating critical stages of viral life cycles. In this context, the employment of specific G4 ligands has aided in comprehending the intricate G4-mediated regulatory mechanisms in the viral life cycle, showcasing the potential of targeting viral G4s as a novel antiviral strategy. This review offers a thorough update on the literature concerning G4s in viruses, including their identification and functional significance across most of the human-infecting viruses. Furthermore, it delves into potential therapeutic avenues targeting G4s, encompassing various G4-binding ligands, G4-interacting proteins, and oligonucleotide-based strategies. Finally, the article highlights both progress and challenges in the field, providing valuable insights into leveraging this unusual nucleic acid structure for therapeutic purposes.

QUFIND: tool for comparative prediction and mining of G4 quadruplexes overlapping with CpG islands

G-quadruplexes (G4s) are secondary structures in DNA that have been shown to be involved in gene regulation. They play a vital role in the cellular processes and several pathogens including bacteria, fungi, and viruses have also been shown to possess G4s that help them in their pathogenesis. Additionally, cross-talk among the CpG islands and G4s has been shown to influence biological processes. The virus-encoded G4s are affected by the mutational landscape leading to the formation/deletion of these G4s. Therefore, understanding and predicting these multivariate effects on traditional and non-traditional quadruplexes forms an important area of research, that is, yet to be investigated. We have designed a user-friendly webserver QUFIND (http://soodlab.com/qufinder/) that can predict traditional as well as non-traditional quadruplexes in a given sequence. QUFIND is connected with ENSEMBL and NCBI so that the sequences can be fetched in a real-time manner. The algorithm is designed in such a way that the user is provided with multiple options to customize the base (A, T, G, or C), size of the stem (2-5), loop length (1-30), number of bulges (1-5) as well as the number of mismatches (0-2) enabling the identification of any of the secondary structure as per their interest. QUFIND is designed to predict both CpG islands as well as G4s in a given sequence. Since G4s are very short as compared to the CpG islands, hence, QUFIND can also predict the overlapping G4s within CpG islands. Therefore, the user has the flexibility to identify either overlapping or non-overlapping G4s along with the CpG islands. Additionally, one section of QUFIND is dedicated to comparing the G4s in two viral sequences. The visualization is designed in such a manner that the user is able to see the unique quadruplexes in both the input sequences. The efficiency of QUFIND is calculated on G4s obtained from G4 high throughput sequencing data (n = 1000) or experimentally validated G4s (n = 329). Our results revealed that QUFIND is able to predict G4-quadruplexes obtained from G4-sequencing data with 90.06% prediction accuracy whereas experimentally validated quadruplexes were predicted with 97.26% prediction accuracy.

A Screening Protocol for Exploring Loop Length Requirements for the Formation of a Three Cytosine-Cytosine+ Base-Paired i-Motif

DNA sequences containing at least four runs of repetitive cytosines can fold into tetra-helical structures called i-Motifs (iMs). The interest in these DNA secondary structures is increasing due to their therapeutical and technological applications. Still, limited knowledge of their folding requirements is currently available. We developed a novel step-by-step pipeline for the systematic screening of putative iM-forming model sequences. Focusing on structures comprising only three cytosine-cytosine+ base pairs, we investigated what the minimal lengths of the loops required for formation of an intra-molecular iM are. Our data indicate that two and three nucleotides are required to connect the strands through the minor and majorgrooves of the iM, respectively. Additionally, they highlight an asymmetric behavior according to the distribution of the cytosines. Specifically, no sequence containing a single cytosine in the first and third run was able to fold into intra-molecular iMs with the same stability of those formed when the first and the third run comprise two cytosines. This knowledge represents a step forward toward the development of prediction tools for the proper identification of biologically functional iMs, as well as for the rational design of these secondary structures as technological devices.

Theranostic approach to specifically targeting the interloop region of BCL2 i-motif DNA by crystal violet

Ligands that recognise specific i-motif DNAs are helpful in cancer diagnostics and therapeutics, as i-motif formation can cause cancer. Although the loop regions of i-motifs are promising targets for ligands, the interaction between a ligand and the loop regions based on sequence information remains unexplored. Herein, we investigated the loop regions of various i-motif DNAs to determine whether these regions specifically interact with fluorescent ligands. Crystal violet (CV), a triphenylmethane dye, exhibited strong fluorescence with the i-motif derived from the promoter region of the human BCL2 gene in a sequence- and structure-specific manner. Our systematic sequence analysis indicated that CV was bound to the site formed by the first and third loops through inter-loop interactions between the guanine bases present in these loops. As the structural stability of the BCL2 i-motif was unaffected by CV, the local stabilisation of the loops by CV could inhibit the interaction of transcription factors with these loops, repressing the BCL2 expression of MCF-7 cells. Our finding suggests that the loops of the i-motif can act as a novel platform for the specific binding of small molecules; thus, they could be utilised for the theranostics of diseases associated with i-motif DNAs.

G-Quadruplexes as Key Transcriptional Regulators in Neglected Trypanosomatid Parasites

G-quadruplexes (G4s) are nucleic acid secondary structures that have been linked to the functional regulation of eukaryotic organisms. G4s have been extensively characterised in humans and emerging evidence suggests that they might also be biologically relevant for human pathogens. This indicates that G4s might represent a novel class of therapeutic targets for tackling infectious diseases. Bioinformatic studies revealed a high prevalence of putative quadruplex-forming sequences (PQSs) in the genome of protozoans, which highlights their potential roles in regulating vital processes of these parasites, including DNA transcription and replication. In this work, we focus on the neglected trypanosomatid parasites, Trypanosoma and Leishmania spp., which cause debilitating and deadly diseases across the poorest populations worldwide. We review three examples where G4-formation might be key to modulate transcriptional activity in trypanosomatids, providing an overview of experimental approaches that can be used to exploit the regulatory roles and relevance of these structures to fight parasitic infections.

Emerging roles of i-motif in gene expression and disease treatment

As non-canonical nucleic acid secondary structures consisting of cytosine-rich nucleic acids, i-motifs can form under certain conditions. Several i-motif sequences have been identified in the human genome and play important roles in biological regulatory functions. Due to their physicochemical properties, these i-motif structures have attracted attention and are new targets for drug development. Herein, we reviewed the characteristics and mechanisms of i-motifs located in gene promoters (including c-myc, Bcl-2, VEGF, and telomeres), summarized various small molecule ligands that interact with them, and the possible binding modes between ligands and i-motifs, and described their effects on gene expression. Furthermore, we discussed diseases closely associated with i-motifs. Among these, cancer is closely associated with i-motifs since i-motifs can form in some regions of most oncogenes. Finally, we introduced recent advances in the applications of i-motifs in multiple areas.

pH-Dependent Capping Interactions Induce Large-Scale Structural Transitions in i-Motifs

We study here a DNA oligonucleotide having the ability to form two different i-motif structures whose relative stability depends on pH and temperature. The major species at neutral pH is stabilized by two C:C+ base pairs capped by two minor groove G:C:G:C tetrads. The high pH and thermal stability of this structure are mainly due to the favorable effect of the minor groove tetrads on their adjacent positively charged C:C+ base pairs. At pH 5, we observe a more elongated i-motif structure consisting of four C:C+ base pairs capped by two G:T:G:T tetrads. Molecular dynamics calculations show that the conformational transition between the two structures is driven by the protonation state of key cytosines. In spite of large conformational differences, the transition between the acidic and neutral structures can occur without unfolding of the i-motif. These results represent the first case of a conformational switch between two different i-motif structures and illustrate the dramatic pH-dependent plasticity of this fascinating DNA motif.

Millimeter waves alter DNA secondary structures and modulate the transcriptome in human fibroblasts

As millimetre wave (MMW) frequencies of the electromagnetic spectrum are increasingly adopted in modern technologies such as mobile communications and networking, characterising the biological effects is critical in determining safe exposure levels. We study the exposure of primary human dermal fibroblasts to MMWs, finding MMWs trigger genomic and transcriptomic alterations. In particular, repeated 60 GHz, 2.6 mW cm^-2, 46.8 J cm^-2 d^-1 MMW doses induce a unique physiological response after 2 and 4 days exposure. We show that high dose MMWs induce simultaneous non-thermal alterations to the transcriptome and DNA structural dynamics, including formation of G-quadruplex and i-motif secondary structures, but not DNA damage.

[G-quadruplex structures in bacteria: functional properties and prospects for use as biotargets]

G-quadruplexes (G4), non-canonical secondary DNA structures, are intensively investigated for a long time. In eukaryotic organisms they play an important role in the regulation of gene expression and DNA repair. G4 have also been found in the genomes of numerous bacteria and archaea, but their functional role has not yet been fully explored. Nevertheless, their participation in the formation of antigenic variability, pathogenicity, antibiotic resistance and survival in extreme conditions has been established. Currently, many tools have been developed to detect potential G4 sequences and confirm their formation ability. Since the controlled formation and resolution of the quadruplex are significant means for the regulation of genes critical for survival, a promising direction is the search for ligands - compounds that can have a stabilizing effect on the quadruplex structure and thereby alter gene expression. Currently, a number of ligands are already known, their use stops the growth of pathogenic microorganisms. G4 ligands are of interest as potential antibiotics, which are extremely relevant due to the wide spread of drug resistant pathogens.

Genome-Wide Transcriptional Response of Mycobacterium smegmatis MC2155 to G-Quadruplex Ligands BRACO-19 and TMPyP4

G-quadruplexes (G4s) are non-canonical DNA structures that could be considered as potential therapeutic targets for antimicrobial compounds, also known as G4-stabilizing ligands. While some of these ligands are shown in vitro to have a stabilizing effect, the precise mechanism of antibacterial action has not been fully investigated. Here, we employed genome-wide RNA-sequencing to analyze the response of Mycobacterium smegmatis to inhibitory concentrations of BRACO-19 and TMPyP4 G4 ligands. The expression profile changed (FDR < 0.05, log₂FC > |1|) for 822 (515↑; 307↓) genes in M. smegmatis in response to BRACO-19 and for 680 (339↑; 341↓) genes in response to TMPyP4. However, the analysis revealed no significant ligand-induced changes in the expression levels of G4-harboring genes, genes under G4-harboring promoters, or intergenic regions located on mRNA-like or template strands. Meanwhile, for the BRACO-19 ligand, we found significant changes in the replication and repair system genes, as well as in iron metabolism genes which is, undoubtedly, evidence of the induced stress. For the TMPyP4 compound, substantial changes were found in transcription factors and the arginine biosynthesis system, which may indicate multiple biological targets for this compound.

G-quadruplex occurrence and conservation: more than just a question of guanine–cytosine content

G-quadruplexes are motifs found in DNA and RNA that can fold into tertiary structures. Until now, they have been studied experimentally mainly in humans and a few other species. Recently, predictions have been made with bacterial and archaeal genomes. Nevertheless, a global comparison of predicted G4s (pG4s) across and within the three living kingdoms has not been addressed. In this study, we aimed to predict G4s in genes and transcripts of all kingdoms of living organisms and investigated the differences in their distributions. The relation of the predictions with GC content was studied. It appears that GC content is not the only parameter impacting G4 predictions and abundance. The distribution of pG4 densities varies depending on the class of transcripts and the group of species. Indeed, we have observed that, in coding transcripts, there are more predicted G4s than expected for eukaryotes but not for archaea and bacteria, while in noncoding transcripts, there are as many or fewer predicted G4s in all species groups. We even noticed that some species with the same GC content presented different pG4 profiles. For instance, Leishmania major and Chlamydomonas reinhardtii both have 60% of GC content, but the former has a pG4 density of 0.07 and the latter 1.16.

Genome-wide characterization of i-motifs and their potential roles in the stability and evolution of transposable elements in rice

I-motifs (iMs) are non-canonical DNA secondary structures that fold from cytosine (C)-rich genomic DNA regions termed putative i-motif forming sequences (PiMFSs). The structure of iMs is stabilized by hemiprotonated C-C base pairs, and their functions are now suspected in key cellular processes in human cells such as genome stability and regulation of gene transcription. In plants, their biological relevance is still largely unknown. Here, we characterized PiMFSs with high potential for i-motif formation in the rice genome by developing and applying a protocol hinging on an iMab antibody-based immunoprecipitation (IP) coupled with high-throughput sequencing (seq), consequently termed iM-IP-seq. We found that PiMFSs had intrinsic subgenomic distributions, cis-regulatory functions and an intricate relationship with DNA methylation. We indeed found that the coordination of PiMFSs with DNA methylation may affect dynamics of transposable elements (TEs) among different cultivated Oryza subpopulations or during evolution of wild rice species. Collectively, our study provides first and unique insights into the biology of iMs in plants, with potential applications in plant biotechnology for improving important agronomic rice traits.

Modulating gene expression in breast cancer via DNA secondary structure and the CRISPR toolbox

Breast cancer is the most commonly diagnosed malignancy in women, and while the survival prognosis of patients with early-stage, non-metastatic disease is ∼75%, recurrence poses a significant risk and advanced and/or metastatic breast cancer is incurable. A distinctive feature of advanced breast cancer is an unstable genome and altered gene expression patterns that result in disease heterogeneity. Transcription factors represent a unique therapeutic opportunity in breast cancer, since they are known regulators of gene expression, including gene expression involved in differentiation and cell death, which are themselves often mutated or dysregulated in cancer. While transcription factors have traditionally been viewed as 'undruggable', progress has been made in the development of small-molecule therapeutics to target relevant protein-protein, protein-DNA and enzymatic active sites, with varying levels of success. However, non-traditional approaches such as epigenetic editing, transcriptional control via CRISPR/dCas9 systems, and gene regulation through non-canonical nucleic acid secondary structures represent new directions yet to be fully explored. Here, we discuss these new approaches and current limitations in light of new therapeutic opportunities for breast cancers.

G-Quadruplex DNA as a Target in Pathogenic Bacteria: Efficacy of an Extended Naphthalene Diimide Ligand and Its Mode of Action

Guanidine DNA quadruplex (G4-DNA) structures convey a distinctive layer of epigenetic information that is critical for regulating key biological activities and processes as transcription, replication, and repair in living cells. The information regarding their role and use as therapeutic drug targets in bacteria is still scarce. Here, we tested the biological activity of a G4-DNA ligand library, based on the naphthalene diimide (NDI) pharmacophore, against both Gram-positive and Gram-negative bacteria. For the best compound identified, NDI-10, a different action mechanism was described for Gram-positive or negative bacteria. This asymmetric activity profile could be related to the different prevalence of putative G4-DNA structures in each group, the influence that they can exert on gene expression, and the different roles of the G4 structures in these bacteria, which seem to promote transcription in Gram-positive bacteria and repress transcription in Gram-negatives.

DNA G-Quadruplexes Contribute to CTCF Recruitment

G-quadruplex (G4) sites in the human genome frequently colocalize with CCCTC-binding factor (CTCF)-bound sites in CpG islands (CGIs). We aimed to clarify the role of G4s in CTCF positioning. Molecular modeling data suggested direct interactions, so we performed in vitro binding assays with quadruplex-forming sequences from CGIs in the human genome. G4s bound CTCF with Kd values similar to that of the control duplex, while respective i-motifs exhibited no affinity for CTCF. Using ChIP-qPCR assays, we showed that G4-stabilizing ligands enhance CTCF occupancy at a G4-prone site in STAT3 gene. In view of the reportedly increased CTCF affinity for hypomethylated DNA, we next questioned whether G4s also facilitate CTCF recruitment to CGIs via protecting CpG sites from methylation. Bioinformatics analysis of previously published data argued against such a possibility. Finally, we questioned whether G4s facilitate CTCF recruitment by affecting chromatin structure. We showed that three architectural chromatin proteins of the high mobility group colocalize with G4s in the genome and recognize parallel-stranded or mixed-topology G4s in vitro. One of such proteins, HMGN3, contributes to the association between G4s and CTCF according to our bioinformatics analysis. These findings support both direct and indirect roles of G4s in CTCF recruitment.

Potential G-quadruplexes and i-Motifs in the SARS-CoV-2

Quadruplex structures have been identified in a plethora of organisms where they play important functions in the regulation of molecular processes, and hence have been proposed as therapeutic targets for many diseases. In this paper we report the extensive bioinformatic analysis of the SARS-CoV-2 genome and related viruses using an upgraded version of the open-source algorithm G4-iM Grinder. This version improves the functionality of the software, including an easy way to determine the potential biological features affected by the candidates found. The quadruplex definitions of the algorithm were optimized for SARS-CoV-2. Using a lax quadruplex definition ruleset, which accepts amongst other parameters two residue G- and C-tracks, 512 potential quadruplex candidates were discovered. These sequences were evaluated by their in vitro formation probability, their position in the viral RNA, their uniqueness and their conservation rates (calculated in over seventeen thousand different COVID-19 clinical cases and sequenced at different times and locations during the ongoing pandemic). These results were then compared subsequently to other Coronaviridae members, other Group IV (+)ssRNA viruses and the entire viral realm. Sequences found in common with other viral species were further analyzed and characterized. Sequences with high scores unique to the SARS-CoV-2 were studied to investigate the variations amongst similar species. Quadruplex formation of the best candidates were then confirmed experimentally. Using NMR and CD spectroscopy, we found several highly stable RNA quadruplexes that may be suitable therapeutic targets for the SARS-CoV-2.

Thermal and pH Stabilities of i-DNA: Confronting in vitro Experiments with Models and In-Cell NMR Data

Recent studies indicate that i-DNA, a four-stranded cytosine-rich DNA also known as the i-motif, is actually formed in vivo; however, a systematic study on sequence effects on stability has been missing. Herein, an unprecedented number of different sequences (271) bearing four runs of 3-6 cytosines with different spacer lengths has been tested. While i-DNA stability is nearly independent on total spacer length, the central spacer plays a special role on stability. Stability also depends on the length of the C-tracts at both acidic and neutral pHs. This study provides a global picture on i-DNA stability thanks to the large size of the introduced data set; it reveals unexpected features and allows to conclude that determinants of i-DNA stability do not mirror those of G-quadruplexes. Our results illustrate the structural roles of loops and C-tracts on i-DNA stability, confirm its formation in cells, and allow establishing rules to predict its stability.

Graphene oxide assisted light-up aptamer selection against Thioflavin T for label-free detection of microRNA

We selected an aptamer against a fluorogenic dye called Thioflavin T (ThT). Aptamers are single-stranded DNA that can bind a specific target. We selected the ThT aptamer using graphene oxide assisted SELEX and a low-cost Open qPCR instrument. We optimized, minimized, and characterized the best aptamer candidate against ThT. The aptamer, ThT dye, and the enzymatic strand displacement amplification (SDA) were used in a label-free approach to detect the micro RNA miR-215 in saliva and serum. The aptamer confers higher specificity than intercalating dyes but without expensive covalently modified DNA probes. This isothermal, low-cost, simple method can detect both DNA and RNA. The target, miR-215, was detected with a limit of detection of 2.6 nM.

Learning to Model G-Quadruplexes: Current Methods and Perspectives

G-quadruplexes have raised considerable interest during the past years for the development of therapies against cancer. These noncanonical structures of DNA may be found in telomeres and/or oncogene promoters, and it has been observed that the stabilization of such G-quadruplexes may disturb tumor cell growth. Nevertheless, the mechanisms leading to folding and stabilization of these G-quadruplexes are still not well established, and they are the focus of much current work in this field. In seminal works, stabilization was observed to be produced by cations. However, subsequent studies showed that different kinds of small molecules, from planar and nonplanar organic molecules to square-planar and octahedral metal complexes, may also lead to the stabilization of G-quadruplexes. Thus, the comprehension and rationalization of the interaction of these small molecules with G-quadruplexes are also important topics of current interest in medical applications. To shed light on the questions arising from the literature on the formation of G-quadruplexes, their stabilization, and their interaction with small molecules, synergies between experimental studies and computational works are needed. In this review, we mainly focus on in silico approaches and provide a broad compilation of different leading studies carried out to date by different computational methods. We divide these methods into twomain categories: (a) classical methods, which allow for long-timescale molecular dynamics simulations and the corresponding analysis of dynamical information, and (b) quantum methods (semiempirical, quantum mechanics/molecular mechanics, and density functional theory methods), which allow for the explicit simulation of the electronic structure of the system but, in general, are not capable of being used in long-timescale molecular dynamics simulations and, therefore, give a more static picture of the relevant processes.

Aureolic Acid Group of Agents as Potential Antituberculosis Drugs

Mycobacterium tuberculosis is one of the most dangerous pathogens. Bacterial resistance to antituberculosis drugs grows each year, but searching for new drugs is a long process. Testing for available drugs to find active against mycobacteria may be a good alternative. In this work, antibiotics of the aureolic acid group were tested on a model organism Mycobacterium smegmatis. We presumed that antibiotics of this group may be potential G4 ligands. However, this was not confirmed in our analyses. We determined the antimicrobial activity of these drugs and revealed morphological changes in the cell structure upon treatment. Transcriptomic analysis documented increased expression of MSMEG_3743/soj and MSMEG_4228/ftsW, involved in cell division. Therefore, drugs may affect cell division, possibly disrupting the function of the Z-ring and the formation of a septum. Additionally, a decrease in the transcription level of several indispensable genes, such as nitrate reductase subunits (MSMEG_5137/narI and MSMEG_5139/narX) and MSMEG_3205/hisD was shown. We concluded that the mechanism of action of aureolic acid and its related compounds may be similar to that bedaquiline and disturb the NAD+/NADH balance in the cell. All of this allowed us to conclude that aureolic acid derivatives can be considered as potential antituberculosis drugs.

pqsfinder web: G-quadruplex prediction using optimized pqsfinder algorithm

MOTIVATION

G-quadruplex is a DNA or RNA form in which four guanine-rich regions are held together by base pairing between guanine nucleotides in coordination with potassium ions. G-quadruplexes are increasingly seen as a biologically important component of genomes. Their detection in vivo is problematic; however, sequencing and spectrometric techniques exist for their in vitro detection. We previously devised the pqsfinder algorithm for PQS identification, implemented it in C++ and published as an R/Bioconductor package. We looked for ways to optimize pqsfinder for faster and user-friendly sequence analysis.

RESULTS

We identified two weak points where pqsfinder could be optimized. We modified the internals of the recursive algorithm to avoid matching and scoring many sub-optimal PQS conformations that are later discarded. To accommodate the needs of a broader range of users, we created a website for submission of sequence analysis jobs that does not require knowledge of R to use pqsfinder.

AVAILABILITY AND IMPLEMENTATION

https://pqsfinder.fi.muni.cz, https://bioconductor.org/packages/pqsfinder.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

How bioinformatics resources work with G4 RNAs

Quadruplexes (G4s) are of interest, which increases with the number of identified G4 structures and knowledge about their biomedical potential. These unique motifs form in many organisms, including humans, where their appearance correlates with various diseases. Scientists store and analyze quadruplexes using recently developed bioinformatic tools—many of them focused on DNA structures. With an expanding collection of G4 RNAs, we check how existing tools deal with them. We review all available bioinformatics resources dedicated to quadruplexes and examine their usefulness in G4 RNA analysis. We distinguish the following subsets of resources: databases, tools to predict putative quadruplex sequences, tools to predict secondary structure with quadruplexes and tools to analyze and visualize quadruplex structures. We share the results obtained from processing specially created RNA datasets with these tools.

Contact: mszachniuk@cs.put.poznan.pl

Supplementary information: Supplementary data are available at Briefings in Bioinformatics online.

Where are G-quadruplexes located in the human transcriptome?

It has been demonstrated that RNA G-quadruplexes (G4) are structural motifs present in transcriptomes and play important regulatory roles in several post-transcriptional mechanisms. However, the full picture of RNA G4 locations and the extent of their implication remain elusive. Solely computational prediction analysis of the whole transcriptome may reveal all potential G4, since experimental identifications are always limited to specific conditions or specific cell lines. The present study reports the first in-depth computational prediction of potential G4 region across the complete human transcriptome. Although using a relatively stringent approach based on three prediction scores that accounts for the composition of G4 sequences, the composition of their neighboring sequences, and the various forms of G4, over 1.1 million of potential G4 (pG4) were predicted. The abundance of G4 was computationally confirmed in both 5' and 3'UTR as well as splicing junction of mRNA, appreciate for the first time in the long ncRNA, while almost absent of most of the small ncRNA families. The present results constitute an important step toward a full understanding of the roles of G4 in post-transcriptional mechanisms.

A guide to computational methods for G-quadruplex prediction

Guanine-rich nucleic acids can fold into the non-B DNA or RNA structures called G-quadruplexes (G4). Recent methodological developments have allowed the characterization of specific G-quadruplex structures in vitro as well as in vivo, and at a much higher throughput, in silico, which has greatly expanded our understanding of G4-associated functions. Typically, the consensus motif G3+N1-7G3+N1-7G3+N1-7G3+ has been used to identify potential G-quadruplexes from primary sequence. Since, various algorithms have been developed to predict the potential formation of quadruplexes directly from DNA or RNA sequences and the number of studies reporting genome-wide G4 exploration across species has rapidly increased. More recently, new methodologies have also appeared, proposing other estimates which consider non-canonical sequences and/or structure propensity and stability. The present review aims at providing an updated overview of the current open-source G-quadruplex prediction algorithms and straightforward examples of their implementation.

Stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells

This review provides the biophysicochemical background and recent advances in stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells.