Composite 5-methylations of cytosines modulate i-motif stability in a sequence-specific manner: Implications for DNA nanotechnology and epigenetic regulation of plant telomeric DNA

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.

Regulation of VEGF gene expression by bisacridine derivative through promoter i-motif for cancer treatment

BACKGROUND

Vascular endothelial growth factor (VEGF) is overexpressed in most malignant tumors, which has important impact on tumor angiogenesis and development. Its gene promoter i-motif structure formed by C-rich sequence can regulate gene expression, which is a promising new target for anti-tumor therapy.

METHODS

We screened various compounds and studied their effects on VEGF through extensive experiments, including SPR, MST, TO displacement, FRET, CD, ESI-MS, NMR, MTT, clone formation, qPCR, Western blot, dual-luciferase reporter assay, immunofluorescence, cell scrape, apoptosis, transwell assay, and animal model.

RESULTS

After extensive screening, bisacridine derivative B09 was found to have selective binding and stabilization to VEGF promoter i-motif, which could down-regulate VEGF gene expression. B09 showed potent inhibition on MCF-7 and HGC-27 cell proliferation and metastasis. B09 significantly inhibited tumor growth in xenograft mice model with HGC-27 cells, showing decreased VEGF expression analyzed through immunohistochemistry.

CONCLUSION

B09 could specifically regulate VEGF gene expression, possibly through interacting with promoter i-motif structure. As a lead compound, B09 could be further developed for innovative anti-cancer agent targeting VEGF.

i-Motif DNA: identification, formation, and cellular functions

An i-motif (iM) is a four-stranded (quadruplex) DNA structure that folds from cytosine (C)-rich sequences. iMs can fold under many different conditions in vitro, which paves the way for their formation in living cells. iMs are thought to play key roles in various DNA transactions, notably in the regulation of genome stability, gene transcription, mRNA translation, DNA replication, telomere and centromere functions, and human diseases. We summarize the different techniques used to assess the folding of iMs in vitro and provide an overview of the internal and external factors that affect their formation and stability in vivo. We describe the possible biological relevance of iMs and propose directions towards their use as target in biology.

In-cell NMR suggests that DNA i-motif levels are strongly depleted in living human cells

I-Motifs (iM) are non-canonical DNA structures potentially forming in the accessible, single-stranded, cytosine-rich genomic regions with regulatory roles. Chromatin, protein interactions, and intracellular properties seem to govern iM formation at sites with i-motif formation propensity (iMFPS) in human cells, yet their specific contributions remain unclear. Using in-cell NMR with oligonucleotide iMFPS models, we monitor iM-associated structural equilibria in asynchronous and cell cycle-synchronized HeLa cells at 37 °C. Our findings show that iMFPS displaying pHT < 7 under reference in vitro conditions occur predominantly in unfolded states in cells, while those with pHT > 7 appear as a mix of folded and unfolded states depending on the cell cycle phase. Comparing these results with previous data obtained using an iM-specific antibody (iMab) reveals that cell cycle-dependent iM formation has a dual origin, and iM formation concerns only a tiny fraction (possibly 1%) of genomic sites with iM formation propensity. We propose a comprehensive model aligning observations from iMab and in-cell NMR and enabling the identification of iMFPS capable of adopting iM structures under physiological conditions in living human cells. Our results suggest that many iMFPS may have biological roles linked to their unfolded states.

The roles of DNA methylation on pH dependent i-motif (iM) formation in rice

I-motifs (iMs) are four-stranded non-B DNA structures containing C-rich DNA sequences. The formation of iMs is sensitive to pH conditions and DNA methylation, although the extent of which is still unknown in both humans and plants. To investigate this, we here conducted iMab antibody-based immunoprecipitation and sequencing (iM-IP-seq) along with bisulfite sequencing using CK (original genomic DNA without methylation-related treatments) and hypermethylated or demethylated DNA at both pH 5.5 and 7.0 in rice, establishing a link between pH, DNA methylation and iM formation on a genome-wide scale. We found that iMs folded at pH 7.0 displayed higher methylation levels than those formed at pH 5.5. DNA demethylation and hypermethylation differently influenced iM formation at pH 7.0 and 5.5. Importantly, CG hypo-DMRs (differentially methylated regions) and CHH (H = A, C and T) hyper-DMRs alone or coordinated with CG/CHG hyper-DMRs may play determinant roles in the regulation of pH dependent iM formation. Thus, our study shows that the nature of DNA sequences alone or combined with their methylation status plays critical roles in determining pH-dependent formation of iMs. It therefore deepens the understanding of the pH and methylation dependent modulation of iM formation, which has important biological implications and practical applications.

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.

Regulation of c-Kit gene transcription selectively by bisacridine derivative through promoter dual i-motif structures

BACKGROUND

c-Kit protein is a signal transduction protein involved in multiple signal pathways, which play an important role in a variety of cellular events such as cell proliferation, apoptosis and differentiation. Special DNA secondary structures on the promoter of c-Kit gene, including G-quadruplex and i-motif structures, could act as "molecular switch" for gene transcriptional regulation, which are potentially important target for development of new anti-cancer drugs.

METHODS

We screened and evaluated the effect of compounds on c-Kit through several experiments, including SPR, FRET, CD, MST, NMR, dual-luciferase reporter assay, Western blot, qPCR, immunofluorescence, MTT assay, colony formation, cell scrape, cell apoptosis, cell cycle analysis, and transwell assay.

RESULTS

After extensive screening, we found that bisacridine derivative B05 had selective binding and stabilization to dual i-motif structures on c-Kit gene promoter, which could down-regulate c-Kit gene transcription and translation, resulting in inhibition of cell proliferation and metastasis. B05 exhibited potent anti-tumor activity on HGC-27 cells, and strongly suppressed tumor growth in HGC-27 xenograft mice model.

CONCLUSIONS

B05 could interact with c-Kit promoter dual i-motif structures with excellent selectivity, which make it possible for selective regulation of gene transcription and translation. B05 could be further developed for selective anti-cancer agent targeting c-Kit promoter i-motifs.

GENERAL SIGNIFICANCE

i-Motifs on different proto-oncogene promoters are diversified, and especially binding of dual i-motifs on the same promoter simultaneously could significantly down-regulate gene transcription with decreased dosage, and therefore increasing the selectivity. This new strategy shed bight light on development of selective DNA-targeting ligands.

How Can CpG Methylations, and Pair-to-Pair Correlations between the Main (Gene) and the Opposite Strands, Suggest a Bending DNA Loop: Insights into the 5'-UTR of DAT1

A working hypothesis issues from patterns of methylation in the 5'-UTR of the DAT1 gene. We considered relationships between pairs of CpGs, of which one on the main-gene strand and another on the complementary opposite strand (COS). We elaborated on data from ADHD children: we calculated all possible combinations of probabilities (estimated by multiplying two raw values of methylation) in pairs of CpGs from either strand. We analyzed all correlations between any given pair and all other pairs. For pairs correlating with M6-M6COS, some pairs had cytosines positioning to the reciprocal right (e.g., M3-M2COS and M6-M5COS), other pairs had cytosines positioning to the reciprocal left (e.g., M2-M3COS; M5-M6COS). Significant pair-to-pair correlations emerged between main-strand and COS CpG pairs. Through graphic representations, we hypothesized that DNA folded to looping conformations: the C1GG C2GG C3GG and C5G C6G motifs would become close enough to allow cytosines 1-2-3 to interact with cytosines 5-6 (on both strands). Data further suggest a sliding, with left- and right-ward oscillations of DNA strands. While thorough empirical verification is needed, we hypothesize simultaneous methylation of main-strand and COS DNA ("methylation dynamics") to serve as a promising biomarker.

Multiscale Modeling of Phosphate···π Contacts in RNA U-Turns Exposes Differences between Quantum-Chemical and AMBER Force Field Descriptions

Phosphate···π, also called anion···π, contacts occur between nucleobases and anionic phosphate oxygens (OP2) in r(GNRA) and r(UNNN) U-turn motifs (N = A,G,C,U; R = A,G). These contacts were investigated using state-of-the-art quantum-chemical methods (QM) to characterize their physicochemical properties and to serve as a reference to evaluate AMBER force field (AFF) performance. We found that phosphate···π interaction energies calculated with the AFF for dimethyl phosphate···nucleobase model systems are less stabilizing in comparison with double-hybrid DFT and that minimum contact distances are larger for all nucleobases. These distance stretches are also observed in large-scale AFF vs QM/MM computations and classical molecular dynamics (MD) simulations on several r(gcGNRAgc) tetraloop hairpins when compared to experimental data extracted from X-ray/cryo-EM structures (res. ≤ 2.5 Å) using the WebFR3D bioinformatic tool. MD simulations further revealed shifted OP2/nucleobase positions. We propose that discrepancies between the QM and AFF result from a combination of missing polarization in the AFF combined with too large AFF Lennard-Jones (LJ) radii of nucleobase carbon atoms in addition to an exaggerated short-range repulsion of the r^-12 LJ repulsive term. We compared these results with earlier data gathered on lone pair···π contacts in CpG Z-steps occurring in r(UNCG) tetraloops. In both instances, charge transfer calculations do not support any significant n → π* donation effects. We also investigated thiophosphate···π contacts that showed reduced stabilizing interaction energies when compared to phosphate···π contacts. Thus, we challenge suggestions that the experimentally observed enhanced thermodynamic stability of phosphorothioated r(GNRA) tetraloops can be explained by larger London dispersion.

Influence of 5-Halogenation on the Base-Pairing Energies of Protonated Cytidine Nucleoside Analogue Base Pairs: Implications for the Stabilities of Synthetic i-Motif Structures for DNA Nanotechnology Applications

DNA nanotechnology has been employed to develop devices based on i-motif structures. The protonated cytosine-cytosine base pairs that stabilize i-motif conformations are favored under slightly acidic conditions. This unique property has enabled development of the first DNA molecular motor driven by pH changes. The ability to alter the stability and pH transition range of such DNA molecular motors is desirable. Understanding how i-motif structures are influenced by modifications, and which modifications enhance stability and/or affect the pH characteristics, are therefore of great interest. Here, the influence of 5-halogenation of the cytosine nucleobases on the base pairing of protonated cytidine nucleoside analogue base pairs is examined using complementary threshold collision-induced dissociation techniques and computational methods. The nucleoside analogues examined here include the 5-halogenated forms of the canonical DNA and RNA cytidine nucleosides. Comparisons among these systems and to the analogous canonical base pairs previously examined enable the influence of 5-halogenation and the 2'-hydroxy substituent on the base pairing to be elucidated. 5-Halogenation of the cytosine nucleobases is found to enhance the strength of base pairing of DNA base pairs and generally weakens the base pairing for RNA base pairs. Trends in the strength of base pairing indicate that both inductive and polarizability effects influence the strength of base pairing. Overall, the present results suggest that 5-halogenation, and in particular, 5-fluorination and 5-iodination, provide effective means of stabilizing DNA i-motif conformations for applications in nanotechnology, whereas only 5-iodination is effective for stabilizing RNA i-motif conformations but the enhancement in stability is less significant.

Stability and context of intercalated motifs (i-motifs) for biological applications

DNA is naturally dynamic and can self-assemble into alternative secondary structures including the intercalated motif (i-motif), a four-stranded structure formed in cytosine-rich DNA sequences. Until recently, i-motifs were thought to be unstable in physiological cellular environments. Studies demonstrating their existence in the human genome and role in gene regulation are now shining light on their biological relevance. Herein, we review the effects of epigenetic modifications on i-motif structure and stability, and biological factors that affect i-motif formation within cells. Furthermore, we highlight recent progress in targeting i-motifs with structure-specific ligands for biotechnology and therapeutic purposes.

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.

Revealing structural peculiarities of homopurine GA repetition stuck by i-motif clip

Non-canonical forms of nucleic acids represent challenging objects for both structure-determination and investigation of their potential role in living systems. In this work, we uncover a structure adopted by GA repetition locked in a parallel homoduplex by an i-motif. A series of DNA oligonucleotides comprising GAGA segment and C3 clip is analyzed by NMR and CD spectroscopies to understand the sequence-structure-stability relationships. We demonstrate how the relative position of the homopurine GAGA segment and the C3 clip as well as single-base mutations (guanine deamination and cytosine methylation) affect base pairing arrangement of purines, i-motif topology and overall stability. We focus on oligonucleotides C3GAGA and methylated GAGAC3 exhibiting the highest stability and structural uniformity which allowed determination of high-resolution structures further analyzed by unbiased molecular dynamics simulation. We describe sequence-specific supramolecular interactions on the junction between homoduplex and i-motif blocks that contribute to the overall stability of the structures. The results show that the distinct structural motifs can not only coexist in the tight neighborhood within the same molecule but even mutually support their formation. Our findings are expected to have general validity and could serve as guides in future structure and stability investigations of nucleic acids.

Effect of DNA modifications on the transition between canonical and non-canonical DNA structures in CpG islands during senescence

Patterns and levels of DNA modifications play important roles in senescence. Two major epigenetic modifications of DNA, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), target CpG sites. Importantly, CpG concentrated regions, known as CpG islands, contain GC-rich sequences, which have the potential to fold into non-canonical DNA structures such as i-motifs and G-quadruplexes. In this study, we investigated the effect of 5mC and 5hmC modifications on the transition between a duplex, and i-motif and G-quadruplexes. To examine the transition, we firstly investigated the stability and structure of the i-motif and G-quadruplexes, considering the molecular environment in senescent cells. Analyses of their stability showed that the modifications did not drastically affect the stability. However, noteworthily, the modification can weaken the (de)stabilisation effect on G-quadruplexes caused by cosolute(s) and cations. Circular dichroism analyses indicated that the surrounding environments, including the molecular crowding and the type of cations such as K⁺ and Na⁺, regulate the topology of G-quadruplexes, while neither 5mC nor 5hmC had a drastic effect. On the other hand, the modifications changed the transition between duplexes and quadruplexes. Unmodified DNA preferred to fold into quadruplexes, whereas DNA with 5mC and 5hmC preferred to fold into duplexes in the absence of PEG200; on the other hand, DNA with or without modifications tended to fold into i-motifs under crowded conditions. Furthermore, an investigation of quadruplexes forming sequences in CpG islands, which are hyper- or hypomethylated during senescence, followed by gene ontology enrichment analysis for each gene group classified by the presence of quadruplexes, showed a difference in function between genes with and without quadruplexes in the CpG region. These results indicate that it is important to consider the effects of patterns and levels of DNA modifications on the transition between canonical and non-canonical DNA structures to understand gene regulation by epigenetic modification during senescence.

The modification of DNA can regulate the transition between a duplex and quadruplexes during senescence responding to surrounding environments.

VEGF aptamer/i-motif-grafted multi-functional SPION nanocarrier for chemotherapeutic/phototherapeutic synergistic research

Chemotherapeutic agents and photosensitizers often suffer from poor tumor selectivity, high side toxicity, or low water solubility. To address these problems, various drug delivery systems (DDS) have been explored but most of them are toxic, difficult to synthesize, or of single function. In order to design a highly biocompatible, conveniently prepared, multi-functional drug delivery system, herein, an aptamer of vascular endothelial growth factor (VEGF) and a cytosine (C)-DNA fragment were grafted on the surface of superparamagnetic iron oxide nanoparticles (SPION), and then a chemotherapeutic agent daunomycin (DNM) and a photosensitizer 5, 10, 15, 20-tetra (phenyl-4-N-methyl-4-pyridyl) porphyrin (TMPyP) were self-assembled with the hybridized VEGF-based DNA structure. By loading DNM and TMPyP, the DDS displayed strong chemotherapeutic/phototherapeutic capability against cancer cells via mechanisms such as mitochondrial dysfunction and ROS elevation, which triggered the apoptosis of the tumor cells. The dual delivery of chemotherapeutical agents and photosensitizers with aptamer/C-rich DNA successfully integrated the functions of pH stimuli-responsive drug release and chemotherapeutic/phototherapeutic modalities into one single system and thus could be considered as an ideal drug delivery vehicle with great potential in clinic.

VEGF aptamer/i-motif-based drug co-delivery system for combined chemotherapy and photodynamic therapy

BACKGROUND

Nucleic acids used as drug delivery systems (DDS) have gained attention because of their biosafety and effortless synthesis. G-quadruplex (G4) structured aptamer such as AS1411 was frequently employed to deliver photosensitizers or chemotherapeutic agents while other aptamers were seldomly reported in this field.

METHODS

Herein, a chemical anticancer drug daunomycin (DNM), and a photosensitizer 5, 10, 15, 20-tetra (phenyl-4-N-methyl-4-pyridyl) porphyrin (TMPyP) were physically assembled with a novel DNA structure composed of an aptamer of vascular endothelial growth factor (VEGF) and a cytosine (C)-rich DNA fragment (gc-34). Spectral and molecular mimicking methods were employed to research the drug loading/releasing process. The in vitro cytotoxicity was studied by MTT, ROS, cell cycle, and cell apoptotic assays and the in vivo anticancer efficiency was evaluated by the inhibitive effect on the cancerous growth of MCF-7 tumor-bearing nude mice.

RESULTS

The G4-structured VEGF aptamer delivered TMPyP successfully for the first time. The designed DDS displayed sensitive VEGF/pH controlled drug release. The co-delivery of DNM and TMPyP exhibited high ROS production, significant cell cycle arresting and evident cell apoptosis, and displayed superior cytotoxicity against tumor cells compared with individual agents in vitro. In vivo studies showed that the dual-drug loaded system can greatly inhibit tumor growth with chemotherapeutic/photodynamic synergistic effects.

CONCLUSION

The co-delivery of DNM and TMPyP with aptamer/C-rich DNA successfully integrates the functions of VEGF/pH stimuli-responsive drug release and chemotherapeutic/phototherapeutic modalities into one single system, and may have great potential in cancer treatment.

Influence of 5-Methylation and the 2'- and 3'-Hydroxy Substituents on the Base Pairing Energies of Protonated Cytidine Nucleoside Analogue Base Pairs: Implications for the Stabilities of i-Motif Structures

Repetitive nucleic acid sequences, which occur in abundance throughout the mammalian genome, are of enormous research interest due to their potential to adopt fascinating and unusual molecular structures such as the i-motif. In remarkable contrast to the DNA double helix, i-motif conformations are stabilized by protonated cytosine base pairs, (Cyt)H⁺(Cyt), that are centrally located in the core of the i-motif and intercalated vertically in an antiparallel fashion. An in-depth understanding of how modifications influence the stability of i-motif conformations is a prerequisite to understanding their biological functions and the development of effective means of tuning their stability for specific medical and technological applications. Here, the influence of the 2'- and 3'-hydroxy substituents of the sugar moieties and 5-methylation of the cytosine nucleobases on the base-pairing interactions of protonated cytidine nucleoside analogue base pairs, (xCyd)H⁺(xCyd), are examined by complementary threshold collision-induced dissociation techniques and computational methods. The xCyd nucleosides examined include the canonical DNA and RNA cytidine nucleosides, 2'-deoxycytidine (dCyd) and cytidine (Cyd), as well as several modified cytidine nucleoside analogues, 2',3'-dideoxycytidine (ddCyd), 5-methyl-2'-deoxycytidine (m⁵dCyd), and 5-methylcytidine (m⁵Cyd). Comparisons among these model base pairs indicate that the 2'- and 3'-hydroxy substituents of the sugar moieties have very little influence on the strength of the base-pairing interactions, whereas 5-methylation of the cytosine nucleobases is found to enhance the strength of the base-pairing interactions. The increase in stability resulting from 5-methylation is only modest but is more than twice as large for the DNA than RNA protonated cytidine base pair. Overall, present results suggest that canonical DNA i-motif conformations should be more stable than analogous RNA i-motif conformations and that 5-methylation of cytosine residues, a significant epigenetic marker, provides greater stabilization to DNA than RNA i-motif conformations.

Evolution of Diverse Strategies for Promoter Regulation

DNA is fundamentally important for all cellular organisms due to its role as a store of hereditary genetic information. The precise and accurate regulation of gene transcription depends primarily on promoters, which vary significantly within and between genomes. Some promoters are rich in specific types of bases, while others have more varied, complex sequence characteristics. However, it is not only base sequence but also epigenetic modifications and altered DNA structure that regulate promoter activity. Significantly, many promoters across all organisms contain sequences that can form intrastrand hairpins (cruciforms) or four-stranded structures (G-quadruplex or i-motif). In this review we integrate recent studies on promoter regulation that highlight the importance of DNA structure in the evolutionary adaptation of promoter sequences.

The i-Motif as a Molecular Target: More Than a Complementary DNA Secondary Structure

Stretches of cytosine-rich DNA are capable of adopting a dynamic secondary structure, the i-motif. When within promoter regions, the i-motif has the potential to act as a molecular switch for controlling gene expression. However, i-motif structures in genomic areas of repetitive nucleotide sequences may play a role in facilitating or hindering expansion of these DNA elements. Despite research on the i-motif trailing behind the complementary G-quadruplex structure, recent discoveries including the identification of a specific i-motif antibody are pushing this field forward. This perspective reviews initial and current work characterizing the i-motif and providing insight into the biological function of this DNA structure, with a focus on how the i-motif can serve as a molecular target for developing new therapeutic approaches to modulate gene expression and extension of repetitive DNA.