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.

Utility of intercalator displacement assays for screening of ligands for i-motif DNA structures.

Cytosine rich sequences can form intercalated, i-motif DNA structures stabilized by hemi-protonated cytosine:cytosine base pairing. These sequences are often located in regulatory regions of genes such as promoters. Ligands targeting i-motif structures may provide potential leads for treatments for genetic disease. The focus on ligands interacting with i-motif DNA has been increasing in recent years. Here, we describe the fluorescent intercalator displacement (FID) assay using thiazole orange binding i-motif DNA and assess the binding affinity of a ligand to the i-motif DNA by displacing thiazole orange. This provides a time and cost-effective high throughput screening of ligands against secondary DNA structures for hit identification.

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.

QN-302 demonstrates opposing effects between i-motif and G-quadruplex DNA structures in the promoter of the S100P gene.

GC-rich sequences can fold into G-quadruplexes and i-motifs and are known to control gene expression in many organisms. The potent G-quadruplex experimental anticancer drug QN-302 down-regulates a number of cancer-related genes, in particular S100P. Here we show this ligand has strong opposing effects with i-motif DNA structures and is one of the most potent i-motif destabilising agents reported to date. QN-302 down-regulates the expression of numerous cancer-related genes by pan-quadruplex targeting. QN-302 exhibits exceptional combined synergistic effects compared to many other G-quadruplex and i-motif interacting compounds. This work further emphasises the importance of considering G-quadruplex and i-motif DNA structures as one dynamic system.

Replication-induced DNA secondary structures drive fork uncoupling and breakage.

Sequences that form DNA secondary structures, such as G-quadruplexes (G4s) and intercalated-Motifs (iMs), are abundant in the human genome and play various physiological roles. However, they can also interfere with replication and threaten genome stability. Multiple lines of evidence suggest G4s inhibit replication, but the underlying mechanism remains unclear. Moreover, evidence of how iMs affect the replisome is lacking. Here, we reconstitute replication of physiologically derived structure-forming sequences to find that a single G4 or iM arrest DNA replication. Direct single-molecule structure detection within solid-state nanopores reveals structures form as a consequence of replication. Combined genetic and biophysical characterisation establishes that structure stability and probability of structure formation are key determinants of replisome arrest. Mechanistically, replication arrest is caused by impaired synthesis, resulting in helicase-polymerase uncoupling. Significantly, iMs also induce breakage of nascent DNA. Finally, stalled forks are only rescued by a specialised helicase, Pif1, but not Rrm3, Sgs1, Chl1 or Hrq1. Altogether, we provide a mechanism for quadruplex structure formation and resolution during replication and highlight G4s and iMs as endogenous sources of replication stress.

The Potent G-Quadruplex-Binding Compound QN-302 Downregulates S100P Gene Expression in Cells and in an In Vivo Model of Pancreatic Cancer.

The naphthalene diimide compound QN-302, designed to bind to G-quadruplex DNA sequences within the promoter regions of cancer-related genes, has high anti-proliferative activity in pancreatic cancer cell lines and anti-tumor activity in several experimental models for the disease. We show here that QN-302 also causes downregulation of the expression of the S100P gene and the S100P protein in cells and in vivo. This protein is well established as being involved in key proliferation and motility pathways in several human cancers and has been identified as a potential biomarker in pancreatic cancer. The S100P gene contains 60 putative quadruplex-forming sequences, one of which is in the promoter region, 48 nucleotides upstream from the transcription start site. We report biophysical and molecular modeling studies showing that this sequence forms a highly stable G-quadruplex in vitro, which is further stabilized by QN-302. We also report transcriptome analyses showing that S100P expression is highly upregulated in tissues from human pancreatic cancer tumors, compared to normal pancreas material. The extent of upregulation is dependent on the degree of differentiation of tumor cells, with the most poorly differentiated, from more advanced disease, having the highest level of S100P expression. The experimental drug QN-302 is currently in pre-IND development (as of Q1 2023), and its ability to downregulate S100P protein expression supports a role for this protein as a marker of therapeutic response in pancreatic cancer. These results are also consistent with the hypothesis that the S100P promoter G-quadruplex is a potential therapeutic target in pancreatic cancer at the transcriptional level for QN-302.

Identification of sugar-containing natural products that interact with i-motif DNA.

There are thousands of compounds shown to interact with G-quadruplex DNA, yet very few which target i-motif (iM) DNA. Previous work showed that tobramycin can interact with iM- DNA, indicating the potential for sugar-molecules to target these structures. Computational approaches indicated that the sugar-containing natural products baicalin and geniposidic acid had potential to target iM-DNA. We assessed the DNA interacting properties of these compounds using FRET-based DNA melting and a fluorescence-based displacement assay using iM-DNA structures from the human telomere and the insulin linked polymorphic region (ILPR), as well as complementary G-quadruplex and double stranded DNA. Both baicalin and geniposidic acid show promise as iM-interacting compounds with potential for use in experiments into the structure and function of i-motif forming DNA sequences and present starting points for further synthetic development of these as probes for iM-DNA.

Activation of the P2Y2 receptor regulates bone cell function by enhancing ATP release.

Bone cells constitutively release ATP into the extracellular environment where it acts locally via P2 receptors to regulate bone cell function. Whilst P2Y2 receptor stimulation regulates bone mineralisation, the functional effects of this receptor in osteoclasts remain unknown. This investigation used the P2Y2 receptor knockout (P2Y2R-/- ) mouse model to investigate the role of this receptor in bone. MicroCT analysis of P2Y2R-/- mice demonstrated age-related increases in trabecular bone volume (≤48%), number (≤30%) and thickness (≤17%). In vitro P2Y2R-/- osteoblasts displayed a 3-fold increase in bone formation and alkaline phosphatase activity, whilst P2Y2R-/- osteoclasts exhibited a 65% reduction in resorptive activity. Serum cross-linked C-telopeptide levels (CTX, resorption marker) were also decreased (≤35%). The resorption defect in P2Y2R-/- osteoclasts was rescued by the addition of exogenous ATP, suggesting that an ATP deficit could be a key factor in the reduced function of these cells. In agreement, we found that basal ATP release was reduced up to 53% in P2Y2R-/- osteoclasts. The P2Y2 receptor agonists, UTP and 2-thioUTP, increased osteoclast activity and ATP release in wild-type but not in P2Y2R-/- cells. This indicates that the P2Y2 receptor may regulate osteoclast function indirectly by promoting ATP release. UTP and 2-thioUTP also stimulate ATP release from osteoblasts suggesting that the P2Y2 receptor exerts a similar function in these cells. Taken together, our findings are consistent with the notion that the primary action of P2Y2 receptor signalling in bone is to regulate extracellular ATP levels.