QGRS Mapper: a web-based server for predicting G-quadruplexes in nucleotide sequences

The quadruplex structures formed by guanine-rich nucleic acid sequences have received significant attention recently because of growing evidence for their role in important biological processes and as therapeutic targets. G-quadruplex DNA has been suggested to regulate DNA replication and may control cellular proliferation. Sequences capable of forming G-quadruplexes in the RNA have been shown to play significant roles in regulation of polyadenylation and splicing events in mammalian transcripts. Whether quadruplex structure directly plays a role in regulating RNA processing requires investigation. Computational approaches to study G-quadruplexes allow detailed analysis of mammalian genomes. There are no known easily accessible user-friendly tools that can compute G-quadruplexes in the nucleotide sequences. We have developed a web-based server, QGRS Mapper, that predicts quadruplex forming G-rich sequences (QGRS) in nucleotide sequences. It is a user-friendly application that provides many options for defining and studying G-quadruplexes. It performs analysis of the user provided genomic sequences, e.g. promoter and telomeric regions, as well as RNA sequences. It is also useful for predicting G-quadruplex structures in oligonucleotides. The program provides options to search and retrieve desired gene/nucleotide sequence entries from NCBI databases for mapping G-quadruplexes in the context of RNA processing sites. This feature is very useful for investigating the functional relevance of G-quadruplex structure, in particular its role in regulating the gene expression by alternative processing. In addition to providing data on composition and locations of QGRS relative to the processing sites in the pre-mRNA sequence, QGRS Mapper features interactive graphic representation of the data. The user can also use the graphics module to visualize QGRS distribution patterns among all the alternative RNA products of a gene simultaneously on a single screen. QGRS Mapper can be accessed at

Analysis of G-Quadruplex Formation in mRNA Transcripts of Phospholemman/FXYD1

G-quadruplexes are higher-order nucleic acid structures formed by tetrads of guanine bases (G-tetrads) through non-canonical base interactions. Two G-tetrads are stabilised by a potassium-ion sandwiched between the tetrads. It has emerged from recent studies that G-quadruplexes occur widely throughout the human genome and have significant biological roles. In this study the FXYD1 pre-mRNA encoding the protein Phospholemman (PLM) is investigated. PLM is highly expressed in cardiomyocytes and forms a third subunit of the Na+/K+ pump (NKA). PLM is a major phosphorylation target and thus regulates NKA activity. FXYD1 pre-mRNA was investigated for its ability to form G-quadruplexes. By computational analysis, it was found that FXYD1 can fold into G-quadruplex and multiple sequence alignment of ortholog FXYD1 sequences indicated that G-quadruplex-forming potential is conserved in evolution, hinting at a potential regulatory mechanism of FXYD1 expression. Comparative analysis confirmed that FXYD1-009, a variant of FXYD1, is a product of alternative splicing of FXYD1’s pre-mRNA. G-quadruplex formation in human and bovine FXYD1-derived oligonucleotides was detected experimentally by non-denaturing poly acrylamide gel electrophoresis that showed an increased mobility rate of G-quadruplexes in contrast to controls. Further analysis by fluorescence emission spectroscopy confirmed G-quadruplex formation in the human and bovine FXYD1-oligonucleotides that was triggered by the presence of K+ ions. The results provided clear evidence of G-quadruplex formation in vitro and together with evolutionary conservation point to potential role in regulating expression of FXYD1 possibly through alternative splicing and thus regulate indirectly the activity of Na+/K+-ATPase. Further in-vivo works should address whether alternative splicing of FXYD1 to FXYD1-009 is associated with G-quadruplex formation

Formation of a Unique Cluster of G-Quadruplex Structures in the HIV-1 nef Coding Region: Implications for Antiviral Activity

G-quadruplexes are tetraplex structures of nucleic acids that can form in G-rich sequences. Their presence and functional role have been established in telomeres, oncogene promoters and coding regions of the human chromosome. In particular, they have been proposed to be directly involved in gene regulation at the level of transcription. Because the HIV-1 Nef protein is a fundamental factor for efficient viral replication, infectivity and pathogenesis in vitro and in vivo, we investigated G-quadruplex formation in the HIV-1 nef gene to assess the potential for viral inhibition through G-quadruplex stabilization. A comprehensive computational analysis of the nef coding region of available strains showed the presence of three conserved sequences that were uniquely clustered. Biophysical testing proved that G-quadruplex conformations were efficiently stabilized or induced by G-quadruplex ligands in all three sequences. Upon incubation with a G-quadruplex ligand, Nef expression was reduced in a reporter gene assay and Nef-dependent enhancement of HIV-1 infectivity was significantly repressed in an antiviral assay. These data constitute the first evidence of the possibility to regulate HIV-1 gene expression and infectivity through G-quadruplex targeting and therefore open a new avenue for viral treatment. © 2013 Perrone et al

Ionic modulation of QPX stability as a nano - switch regulating gene expression in neurons

G-quadruplexes (G-QPX) have been the subject of intense research due to their unique structural configuration and potential applications, particularly their functionality in biological process as a novel type of nano–switch. They have been found in critical regions of the human genome such as telomeres, promoter regions, and untranslated regions of RNA. About 50% of human DNA in promoters has G-rich regions with the potential to form G-QPX structures. A G-QPX might act mechanistically as an ON/OFF switch, regulating gene expression, meaning that the formation of G-QPX in a single strand of DNA disrupts double stranded DNA, prevents the binding of transcription factors (TF) to their recognition sites, resulting in gene down-regulation. Although there are numerous studies on biological roles of G-QPXs in oncogenes, their potential formation in neuronal cells, in particular upstream of transcription start sites, is poorly investigated. The main focus of this research is to identify stable G-QPXs in the 97bp active promoter region of the choline acetyltransferase (ChAT) gene, the terminal enzyme involved in synthesis of the neurotransmitter acetylcholine, and to clarify ionic modulation of G-QPX nanostructures through the mechanism of neural action potentials. Different bioinformatics analyses (in silico), including the QGRS, quadparser and G4-Calculator programs, have been used to predict stable G-QPX in the active promoter region of the human ChAT gene, located 1000bp upstream from the TATA box. The results of computational studies (using those three different algorithms) led to the identification of three consecutive intramolecular G-QPX structures in the negative strand (ChAT G17-2, ChAT G17, and ChAT G29) and one intramolecular G-QPX structure in the positive strand (ChAT G30). Also, the results suggest the possibility that nearby G-runs in opposed DNA strands with a short distance of each other may be able to form a stable intermolecular G-QPX involving two DNA complementary strands (ds ChAT G21). Formation of G-QPX structures, by blocking the availability of the transcription factor binding site (TFBS) on double stranded DNA, can interfere with transcriptional activation. This suggests that there is competition between TFBS binding to dsDNA and the conversion to high order non-B form secondary structures (G-QPXs) in the active promoter region. TFBS mapping analysis of the active promoter region of the human ChAT gene revealed that it contains multiple consensus AP-2a and Sp1 binding sites and consensus sites for other TF, including multiple sites for GR-alpha, Pax-5, p53 and GC box proteins. To get a better understanding of how modulation of G-QPX structures might affect the ChAT promoter activity, an artificial GFP reporter vector (modified GFP) was constructed, synthesized and used for reporter gene measurement. As known human ChAT promoter activators, nerve growth factors (HNGF and TGB) and cytokines (IL-ß and TNF-a) were used for activation of the artificial promoter driving GFP. Also, the G-QPX stabilizing drug TMPYP4 and aconitine, a Na+ channel opening drug, were used as G-QPX stability modulating factors. It was observed that aconitine potentiated the action of the transcriptional activator NGF, suggesting that the effect of sodium is contrary to that of TMPYP4, i.e., that an increase in promoter activity may be due to instability of G-QPX structures in a high Na+ environment, which results in melting these structures, enabling dsDNA formation required for the binding of TF to their recognition sites for initiation of transcription. The results were confirmed in several independent sets of experiments, using GFP reporter gene measurement by plate reader, by flow cytometry and using fluorescent microscopy. Moreover, quantitative RT-PCR was conducted to evaluate the effect of the same factors under similar conditions on the actual ChAT mRNA expression. It was observed that TMPY4 knocked down the ChAT mRNA expression by 87%, suggesting that G-QPX stabilization inhibits promoter activity as expected and that aconitine along with HNGF increases ChAT mRNA expression up to 2.8 fold. Aconitine-mediated influx of Na+ ions, possibly by inhibiting the formation of stable G-QPX structures, resulted in an Unique G-QPX structures can be stabilized with certain metal cations or small cationic molecule ligands such as TMPYP4, through occupying the space between the layers of G-tetrads. Although G-QPX are reported to have high stability in potassium solution, the diversity of G-QPX structures (due to diversity in sequence and size of G-runs, sequence and size of loops) will lead to diversity in physical behavior of G-QPX structures. Therefore, to get a clear image of folding topology and stability of identified G-QPX structures, physical studies including CD spectroscopy and AFM imaging were conducted. CD results showed that the identified ChAT G-QPX structures formed a hybrid, stable configuration in potassium environment (10mM) while being instable in sodium solution (100mM). AFM imaging demonstrated star-shaped structures (involving clusters of DNA strands) due to incubation with TMPYP4, where a greater number of these G-rich sequences have converted to G-QPX structures. The results of both an artificial engineered reporter gene system and actual ChAT mRNA expression (in vitro), plus physical characterization studies, strongly support the novel hypothesis that a neural action potential ionic mechanism regulates G-QPX formation/ deformation in the promoter region, due to movement of monovalent cations across the membrane, which is consistent with gene silencing and expression during neuronal resting and firing

Predicting and understanding the stability of G-quadruplexes

Motivation: G-quadruplexes are stable four-stranded guanine-rich structures that can form in DNA and RNA. They are an important component of human telomeres and play a role in the regulation of transcription and translation. The biological significance of a G-quadruplex is crucially linked with its thermodynamic stability. Hence the prediction of G-quadruplex stability is of vital interest

QuadBase2: web server for multiplexed guanine quadruplex mining and visualization

DNA guanine quadruplexes or G4s are non-canonical DNA secondary structures which affect genomic processes like replication, transcription and recombination. G4s are computationally identified by specific nucleotide motifs which are also called putative G4 (PG4) motifs. Despite the general relevance of these structures, there is currently no tool available that can allow batch queries and genome-wide analysis of these motifs in a user-friendly interface. QuadBase2 (quadbase.igib.res.in) presents a completely reinvented web server version of previously published QuadBase database. QuadBase2 enables users to mine PG4 motifs in up to 178 eukaryotes through the EuQuad module. This module interfaces with Ensembl Compara database, to allow users mine PG4 motifs in the orthologues of genes of interest across eukaryotes. PG4 motifs can be mined across genes and their promoter sequences in 1719 prokaryotes through ProQuad module. This module includes a feature that allows genome-wide mining of PG4 motifs and their visualization as circular histograms. TetraplexFinder, the module for mining PG4 motifs in user-provided sequences is now capable of handling up to 20 MB of data. QuadBase2 is a comprehensive PG4 motif mining tool that further expands the configurations and algorithms for mining PG4 motifs in a user-friendly way

The disruptive positions in human G-quadruplex motifs are less polymorphic and more conserved than their neutral counterparts

Specific guanine-rich sequence motifs in the human genome have considerable potential to form four-stranded structures known as G-quadruplexes or G4 DNA. The enrichment of these motifs in key chromosomal regions has suggested a functional role for the G-quadruplex structure in genomic regulation. In this work, we have examined the spectrum of nucleotide substitutions in G4 motifs, and related this spectrum to G4 prevalence. Data collected from the large repository of human SNPs indicates that the core feature of G-quadruplex motifs, 5′-GGG-3′, exhibits specific mutational patterns that preserve the potential for G4 formation. In particular, we find a genome-wide pattern in which sites that disrupt the guanine triplets are more conserved and less polymorphic than their neutral counterparts. This also holds when considering non-CpG sites only. However, the low level of polymorphisms in guanine tracts is not only confined to G4 motifs. A complete mapping of DNA three-mers at guanine polymorphisms indicated that short guanine tracts are the most under-represented sequence context at polymorphic sites. Furthermore, we provide evidence for a strand bias upstream of human genes. Here, a significantly lower rate of G4-disruptive SNPs on the non-template strand supports a higher relative influence of G4 formation on this strand during transcription

The disruptive positions in human G-quadruplex motifs are less polymorphic and more conserved than their neutral counterparts

Specific guanine-rich sequence motifs in the human genome have considerable potential to form four-stranded structures known as G-quadruplexes or G4 DNA. The enrichment of these motifs in key chromosomal regions has suggested a functional role for the G-quadruplex structure in genomic regulation. In this work, we have examined the spectrum of nucleotide substitutions in G4 motifs, and related this spectrum to G4 prevalence. Data collected from the large repository of human SNPs indicates that the core feature of G-quadruplex motifs, 5′-GGG-3′, exhibits specific mutational patterns that preserve the potential for G4 formation. In particular, we find a genome-wide pattern in which sites that disrupt the guanine triplets are more conserved and less polymorphic than their neutral counterparts. This also holds when considering non-CpG sites only. However, the low level of polymorphisms in guanine tracts is not only confined to G4 motifs. A complete mapping of DNA three-mers at guanine polymorphisms indicated that short guanine tracts are the most under-represented sequence context at polymorphic sites. Furthermore, we provide evidence for a strand bias upstream of human genes. Here, a significantly lower rate of G4-disruptive SNPs on the non-template strand supports a higher relative influence of G4 formation on this strand during transcription

Mapping the sequences of potential guanine quadruplex motifs.

The knowledge that potential guanine quadruplex sequences (PQs) are non-randomly distributed in relation to genomic features is now well established. However, this is for a general potential quadruplex motif which is characterized by short runs of guanine separated by loop regions, regardless of the nature of the loop sequence. There have been no studies to date which map the distribution of PQs in terms of primary sequence or which categorize PQs. To this end, we have generated clusters of PQ sequence groups of various sizes and various degrees of similarity for the non-template strand of introns in the human genome. We started with 86 697 sequences, and successively merged them into groups based on sequence similarity, carrying out 66 clustering cycles before convergence. We have demonstrated here that by using complete linkage hierarchical agglomerative clustering such PQ sequence categorization can be achieved. Our results give an insight into sequence diversity and categories of PQ sequences which occur in human intronic regions. We also highlight a number of clusters for which interesting relationships among their members were immediately evident and other clusters whose members seem unrelated, illustrating, we believe, a distinct role for different sequence types