Quantitative analysis and prediction of G-quadruplex forming sequences in double-stranded DNA.

G-quadruplex (GQ) is a four-stranded DNA structure that can be formed in guanine-rich sequences. GQ structures have been proposed to regulate diverse biological processes including transcription, replication, translation and telomere maintenance. Recent studies have demonstrated the existence of GQ DNA in live mammalian cells and a significant number of potential GQ forming sequences in the human genome. We present a systematic and quantitative analysis of GQ folding propensity on a large set of 438 GQ forming sequences in double-stranded DNA by integrating fluorescence measurement, single-molecule imaging and computational modeling. We find that short minimum loop length and the thymine base are two main factors that lead to high GQ folding propensity. Linear and Gaussian process regression models further validate that the GQ folding potential can be predicted with high accuracy based on the loop length distribution and the nucleotide content of the loop sequences. Our study provides important new parameters that can inform the evaluation and classification of putative GQ sequences in the human genome.

G-quadruplex formation in double strand DNA probed by NMM and CV fluorescence.

G-quadruplexes (GQs) are alternative DNA secondary structures that can form throughout the human genome and control the replication and transcription of important regulatory genes. Here, we established an ensemble fluorescence assay by employing two GQ-interacting compounds, N-methyl mesoporphyrin IX (NMM) and Crystal Violet (CV). This enables quantitative measurement of the GQ folding propensity and conformation specificity in both single strand (ss) and double strand (ds) DNA. Our GQ mapping indicates that the likelihood of GQ formation is substantially diminished in dsDNA, likely due to the competition from the Watson-Crick base pairing. Unlike GQ folding sequence in ssDNA which forms both parallel and antiparallel GQs, dsDNA displays only parallel folding. Additionally, we employed single molecule FRET to obtain a direct quantitation of stably formed-, weakly folded and unfolded GQ conformations. The findings of this study and the method developed here will enable identifying and classifying potential GQ-forming sequences in human genome.

Cancer. The transcription factor GABP selectively binds and activates the mutant TERT promoter in cancer.

Reactivation of telomerase reverse transcriptase (TERT) expression enables cells to overcome replicative senescence and escape apoptosis, which are fundamental steps in the initiation of human cancer. Multiple cancer types, including up to 83% of glioblastomas (GBMs), harbor highly recurrent TERT promoter mutations of unknown function but specific to two nucleotide positions. We identified the functional consequence of these mutations in GBMs to be recruitment of the multimeric GA-binding protein (GABP) transcription factor specifically to the mutant promoter. Allelic recruitment of GABP is consistently observed across four cancer types, highlighting a shared mechanism underlying TERT reactivation. Tandem flanking native E26 transformation-specific motifs critically cooperate with these mutations to activate TERT, probably by facilitating GABP heterotetramer binding. GABP thus directly links TERT promoter mutations to aberrant expression in multiple cancers.

Telomeric overhang length determines structural dynamics and accessibility to telomerase and ALT-associated proteins.

The G-rich single-stranded DNA at the 3' end of human telomeres can self-fold into G-quaduplex (GQ). However, telomere lengthening by telomerase or the recombination-based alternative lengthening of telomere (ALT) mechanism requires protein loading on the overhang. Using single-molecule fluorescence spectroscopy, we discovered that lengthening the telomeric overhang also increased the rate of dynamic exchanges between structural conformations. Overhangs with five to seven TTAGGG repeats, compared with four repeats, showed much greater dynamics and accessibility to telomerase binding and activity and loading of the ALT-associated proteins RAD51, WRN, and BLM. Although the eight repeats are highly dynamic, they can fold into two GQs, which limited protein accessibility. In contrast, the telomere-specific protein POT1 is unique in that it binds independently of repeat number. Our results suggest that the telomeric overhang length and dynamics may contribute to the regulation of telomere extension via telomerase action and the ALT mechanism.