G-quadruplex virtual drug screening: A review

Bioactive nutraceuticals as G4 stabilizers: potential cancer prevention and therapy-a critical review

G-quadruplexes (G4) are non-canonical, four-stranded, nucleic acid secondary structures formed in the guanine-rich sequences, where guanine nucleotides associate with each other via Hoogsteen hydrogen bonding. These structures are widely found near the functional regions of the mammalian genome, such as telomeres, oncogenic promoters, and replication origins, and play crucial regulatory roles in replication and transcription. Destabilization of G4 by various carcinogenic agents allows oncogene overexpression and extension of telomeric ends resulting in dysregulation of cellular growth-promoting oncogenesis. Therefore, targeting and stabilizing these G4 structures with potential ligands could aid cancer prevention and therapy. The field of G-quadruplex targeting is relatively nascent, although many articles have demonstrated the effect of G4 stabilization on oncogenic expressions; however, no previous study has provided a comprehensive analysis about the potency of a wide variety of nutraceuticals and some of their derivatives in targeting G4 and the lattice of oncogenic cell signaling cascade affected by them. In this review, we have discussed bioactive G4-stabilizing nutraceuticals, their sources, mode of action, and their influence on cellular signaling, and we believe our insight would bring new light to the current status of the field and motivate researchers to explore this relatively poorly studied arena.

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.

Probing the Competition between Duplex, G-Quadruplex and i-Motif Structures of the Oncogenic c-Myc DNA Promoter Region

Development of probe systems that provide unique spectral signatures for duplex, G-quadruplex (GQ) and i-motif (iM) structures is very important to understand the relative propensity of a G-rich-C-rich promoter region to form these structures. Here, we devise a platform using a combination of two environment-sensitive nucleoside analogs namely, 5-fluorobenzofuran-modified 2'-deoxyuridine (FBF-dU) and 5-fluoro-2'-deoxyuridine (F-dU) to study the structures adopted by a promoter region of the c-Myc oncogene. FBF-dU serves as a dual-purpose probe containing a fluorescent and 19 F NMR label. When incorporated into the C-rich sequence, it reports the formation of different iMs via changes in its fluorescence properties and 19 F signal. F-dU incorporated into the G-rich ON reports the formation of a GQ structure whose 19 F signal is clearly different from the signals obtained for iMs. Rewardingly, the labeled ONs when mixed with respective complementary strands allows us to determine the relative population of different structures formed by the c-Myc promoter by the virtue of the probe's ability to produce distinct and resolved 19 F signatures for different structures. Our results indicate that at physiological pH and temperature the c-Myc promoter forms duplex, random coil and GQ structures, and does not form an iM. Whereas at acidic pH, the mixture largely forms iM and GQ structures. Taken together, our system will complement existing tools and provide unprecedented insights on the population equilibrium and dynamics of nucleic acid structures under different conditions.

G4-QuadScreen: A Computational Tool for Identifying Multi-Target-Directed Anticancer Leads against G-Quadruplex DNA

The study presents 'G4-QuadScreen', a user-friendly computational tool for identifying MTDLs against G4s. Also, it offers a few hit MTDLs based on in silico and in vitro approaches. Multi-tasking QSAR models were developed using linear discriminant analysis and random forest machine learning techniques for predicting the responses of interest (G4 interaction, G4 stabilization, G4 selectivity, and cytotoxicity) considering the variations in the experimental conditions (e.g., G4 sequences, endpoints, cell lines, buffers, and assays). A virtual screening with G4-QuadScreen and molecular docking using YASARA (AutoDock-Vina) was performed. G4 activities were confirmed via FRET melting, FID, and cell viability assays. Validation metrics demonstrated the high discriminatory power and robustness of the models (the accuracy of all models is ~>90% for the training sets and ~>80% for the external sets). The experimental evaluations showed that ten screened MTDLs have the capacity to selectively stabilize multiple G4s. Three screened MTDLs induced a strong inhibitory effect on various human cancer cell lines. This pioneering computational study serves a tool to accelerate the search for new leads against G4s, reducing false positive outcomes in the early stages of drug discovery. The G4-QuadScreen tool is accessible on the ChemoPredictionSuite website.

Higher-order G-quadruplexes in promoters are untapped drug targets

G-quadruplexes (G4s) are four-stranded nucleic acid secondary structures that form within guanine-rich regions of chromatin. G4 motifs are abundant in the genome, with a sizable proportion (∼40%) existing within gene promoter regions. G4s are proven epigenetic features that decorate the promoter landscape as binding centers for transcription factors. Stabilizing or disrupting promoter G4s can directly influence adjacent gene transcription, making G4s attractive as indirect drug targets for hard-to-target proteins, particularly in cancer. However, no G4 ligands have progressed through clinical trials, mostly owing to off targeting effects. A major hurdle in G4 drug discovery is the lack of distinctiveness of the small monomeric G4 structures currently used as receptors. This mini review describes and contrasts monomeric and higher-order G-quadruplex structure and function and provides a rationale for switching focus to the higher-order forms as selective molecular targets. The human telomerase reverse transcriptase (hTERT) core promoter G-quadruplex is then used as a case study that highlights the potential for higher-order G4s as selective indirect inhibitors of hard-to-target proteins in cancer.

Inhibited complete folding of consecutive human telomeric G-quadruplexes

Noncanonical DNA structures, termed G-quadruplexes, are present in human genomic DNA and are important elements in many DNA metabolic processes. Multiple sites in the human genome have G-rich DNA stretches able to support formation of several consecutive G-quadruplexes. One of those sites is the telomeric overhang region that has multiple repeats of TTAGGG and is tightly associated with both cancer and aging. We investigated the folding of consecutive G-quadruplexes in both potassium- and sodium-containing solutions using single-molecule FRET spectroscopy, circular dichroism, thermal melting and molecular dynamics simulations. Our observations show coexistence of partially and fully folded DNA, the latter consisting of consecutive G-quadruplexes. Following the folding process over hours in sodium-containing buffers revealed fast G-quadruplex folding but slow establishment of thermodynamic equilibrium. We find that full consecutive G-quadruplex formation is inhibited by the many DNA structures randomly nucleating on the DNA, some of which are off-path conformations that need to unfold to allow full folding. Our study allows describing consecutive G-quadruplex formation in both nonequilibrium and equilibrium conditions by a unified picture, where, due to the many possible DNA conformations, full folding with consecutive G-quadruplexes as beads on a string is not necessarily achieved.

Probing the Efficiency of 13-Pyridylalkyl Berberine Derivatives to Human Telomeric G-Quadruplexes Binding: Spectroscopic, Solid State and In Silico Analysis

The interaction between the series of berberine derivatives 1-5 (NAX071, NAX120, NAX075, NAX077 and NAX079) and human telomeric G-quadruplexes (G4), which are able to inhibit the Telomerase enzyme's activity in malignant cells, was investigated. The derivatives bear a pyridine moiety connected by a hydrocarbon linker of varying length (n = 1-5, with n number of aliphatic carbon atoms) to the C13 position of the parent berberine. As for the G4s, both bimolecular 5'-TAGGGTTAGGGT-3' (Tel12) and monomolecular 5'-TAGGGTTAGGGTTAGGGTTAGGG-3' (Tel23) DNA oligonucleotides were considered. Spectrophotometric titrations, melting tests, X-ray diffraction solid state analysis and in silico molecular dynamics (MD) simulations were used to describe the different systems. The results were compared in search of structure-activity relationships. The analysis pointed out the formation of 1:1 complexes between Tel12 and all ligands, whereas both 1:1 and 2:1 ligand/G4 stoichiometries were found for the adduct formed by NAX071 (n = 1). Tel12, with tetrads free from the hindrance by the loop, showed a higher affinity. The details of the different binding geometries were discussed, highlighting the importance of H-bonds given by the berberine benzodioxole group and a correlation between the strength of binding and the hydrocarbon linker length. Theoretical (MD) and experimental (X-ray) structural studies evidence the possibility for the berberine core to interact with one or both G4 strands, depending on the constraints given by the linker length, thus affecting the G4 stabilization effect.

G-Quadruplex Linked DNA Guides Selective Transfection into Nucleolin-Overexpressing Cancer Cells

Gene therapy is a promising approach for treating tumors. Conventional approaches of DNA delivery depending on non-viral or viral vectors are unsatisfactory due to the concerns of biosafety and cell-targeting efficiency. The question how to deliver DNA into tumor cells efficiently and selectively is a major technological problem in tumor gene therapy. Here, we develop a vector-free gene transfer strategy to deliver genes effectively and selectively by taking advantage of targeting nucleolin. Nucleolin, a shuttle protein moving between cell membrane, cytoplasm and nuclei, is overexpressed in tumor cells. It has a natural ligand G-quadruplex (Gq). Gq-linked DNA (Gq-DNA) is likely to be internalized by ligand dependent uptake mechanisms independently of vectors after neutralizing negative charges of cell membrane by targeting nucleolin. This strategy is referred to as Gq-DNA transfection. Benefiting from its high affinity to nucleolin, Gq-DNA can be effectively delivered into nucleolin-positive tumor cells even nuclei. Gq-DNA transfection is characterized by low cytotoxicity, high efficiency, ease of synthesis, high stability in serum, direct access into nuclei, and specific nucleolin-positive tumor cell targeting.

Drug discovery of small molecules targeting the higher-order hTERT promoter G-quadruplex

DNA G-quadruplexes (G4s) are now widely accepted as viable targets in the pursuit of anticancer therapeutics. To date, few small molecules have been identified that exhibit selectivity for G4s over alternative forms of DNA, such as the ubiquitous duplex. We posit that the lack of current ligand specificity arises for multiple reasons: G4 atomic models are often small, monomeric, single quadruplex structures with few or no druggable pockets; targeting G-tetrad faces frequently results in the enrichment of extended electron-deficient polyaromatic end-pasting scaffolds; and virtual drug discovery efforts often under-sample chemical search space. We show that by addressing these issues we can enrich for non-standard molecular templates that exhibit high selectivity towards G4s over other forms of DNA. We performed an extensive virtual screen against the higher-order hTERT core promoter G4 that we have previously characterized, targeting 12 of its unique loop and groove pockets using libraries containing 40 million drug-like compounds for each screen. Using our drug discovery funnel approach, which utilizes high-throughput fluorescence thermal shift assay (FTSA) screens, microscale thermophoresis (MST), and orthogonal biophysical methods, we have identified multiple unique G4 binding scaffolds. We subsequently used two rounds of catalogue-based SAR to increase the affinity of a disubstituted 2-aminoethyl-quinazoline that stabilizes the higher-order hTERT G-quadruplex by binding across its G4 junctional sites. We show selectivity of its binding affinity towards hTERT is virtually unaffected in the presence of near-physiological levels of duplex DNA, and that this molecule downregulates hTERT transcription in breast cancer cells.

Long promoter sequences form higher-order G-quadruplexes: an integrative structural biology study of c-Myc, k-Ras and c-Kit promoter sequences

We report on higher-order G-quadruplex structures adopted by long promoter sequences obtained by an iterative integrated structural biology approach. Our approach uses quantitative biophysical tools (analytical ultracentrifugation, small-angle X-ray scattering, and circular dichroism spectroscopy) combined with modeling and molecular dynamics simulations, to derive self-consistent structural models. The formal resolution of our approach is 18 angstroms, but in some cases structural features of only a few nucleotides can be discerned. We report here five structures of long (34-70 nt) wild-type sequences selected from three cancer-related promoters: c-Myc, c-Kit and k-Ras. Each sequence studied has a unique structure. Three sequences form structures with two contiguous, stacked, G-quadruplex units. One longer sequence from c-Myc forms a structure with three contiguous stacked quadruplexes. A longer c-Kit sequence forms a quadruplex-hairpin structure. Each structure exhibits interfacial regions between stacked quadruplexes or novel loop geometries that are possible druggable targets. We also report methodological advances in our integrated structural biology approach, which now includes quantitative CD for counting stacked G-tetrads, DNaseI cleavage for hairpin detection and SAXS model refinement. Our results suggest that higher-order quadruplex assemblies may be a common feature within the genome, rather than simple single quadruplex structures.

Major Achievements in the Design of Quadruplex-Interactive Small Molecules

Organic small molecules that can recognize and bind to G-quadruplex and i-Motif nucleic acids have great potential as selective drugs or as tools in drug target discovery programs, or even in the development of nanodevices for medical diagnosis. Hundreds of quadruplex-interactive small molecules have been reported, and the challenges in their design vary with the intended application. Herein, we survey the major achievements on the therapeutic potential of such quadruplex ligands, their mode of binding, effects upon interaction with quadruplexes, and consider the opportunities and challenges for their exploitation in drug discovery.

Structured Waters Mediate Small Molecule Binding to G-Quadruplex Nucleic Acids

The role of G-quadruplexes in human cancers is increasingly well-defined. Accordingly, G-quadruplexes can be suitable drug targets and many small molecules have been identified to date as G-quadruplex binders, some using computer-based design methods and co-crystal structures. The role of bound water molecules in the crystal structures of G-quadruplex-small molecule complexes has been analyzed in this study, focusing on the water arrangements in several G-quadruplex ligand complexes. One is the complex between the tetrasubstituted naphthalene diimide compound MM41 and a human intramolecular telomeric DNA G-quadruplex, and the others are in substituted acridine bimolecular G-quadruplex complexes. Bridging water molecules form most of the hydrogen-bond contacts between ligands and DNA in the parallel G-quadruplex structures examined here. Clusters of structured water molecules play essential roles in mediating between ligand side chain groups/chromophore core and G-quadruplex. These clusters tend to be conserved between complex and native G-quadruplex structures, suggesting that they more generally serve as platforms for ligand binding, and should be taken into account in docking and in silico studies.

Evaluating Molecular Docking Software for Small Molecule Binding to G-Quadruplex DNA

G-quadruplexes are four-stranded nucleic acid secondary structures of biological significance and have emerged as an attractive drug target. The G4 formed in the MYC promoter (MycG4) is one of the most studied small-molecule targets, and a model system for parallel structures that are prevalent in promoter DNA G4s and RNA G4s. Molecular docking has become an essential tool in structure-based drug discovery for protein targets, and is also increasingly applied to G4 DNA. However, DNA, and in particular G4, binding sites differ significantly from protein targets. Here we perform the first systematic evaluation of four commonly used docking programs (AutoDock Vina, DOCK 6, Glide, and RxDock) for G4 DNA-ligand binding pose prediction using four small molecules whose complex structures with the MycG4 have been experimentally determined in solution. The results indicate that there are considerable differences in the performance of the docking programs and that DOCK 6 with GB/SA rescoring performs better than the other programs. We found that docking accuracy is mainly limited by the scoring functions. The study shows that current docking programs should be used with caution to predict G4 DNA-small molecule binding modes.

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

Combining Electrospray Mass Spectrometry (ESI-MS) and Computational Techniques in the Assessment of G-Quadruplex Ligands: A Hybrid Approach to Optimize Hit Discovery

Guanine-rich sequences forming G-quadruplexes (GQs) are present in several genomes, ranging from viral to human. Given their peculiar localization, the induction of GQ formation or GQ stabilization with small molecules represents a strategy for interfering with crucial biological functions. Investigating the recognition event at the molecular level, with the aim of fully understanding the triggered pharmacological effects, is challenging. Native electrospray ionization mass spectrometry (ESI-MS) is being optimized to study these noncovalent assemblies. Quantitative parameters retrieved from ESI-MS studies, such as binding affinity, the equilibrium binding constant, and sequence selectivity, will be overviewed. Computational experiments supporting the ESI-MS investigation and boosting its efficiency in the search for GQ ligands will also be discussed with practical examples. The combination of ESI-MS and in silico techniques in a hybrid high-throughput-screening workflow represents a valuable tool for the medicinal chemist, providing data on the quantitative and structural aspects of ligand–GQ interactions.

Folding and persistence times of intramolecular G-quadruplexes transiently embedded in a DNA duplex

G-quadruplex (G4) DNA structures have emerged as important regulatory elements during DNA metabolic transactions. While many in vitro studies have focused on the kinetics of G4 formation within DNA single-strands, G4 are found in vivo in double-stranded DNA regions, where their formation is challenged by the complementary strand. Since the energy of hybridization of Watson-Crick structures dominates the energy of G4 folding, this competition should play a critical role on G4 persistence. To address this, we designed a single-molecule assay allowing to measure G4 folding and persistence times in the presence of the complementary strand. We quantified both folding and unfolding rates of biologically relevant G4 sequences, such as the cMYC and cKIT oncogene promoters, human telomeres and an avian replication origin. We confirmed that G4s are found much more stable in tested replication origin and promoters than in human telomere repeats. In addition, we characterized how G4 dynamics was affected by G4 ligands and showed that both folding rate and persistence time increased. Our assay opens new perspectives for the measurement of G4 dynamics in double-stranded DNA mimicking a replication fork, which is important to understand their role in DNA replication and gene regulation at a mechanistic level.

Structural recognition of the MYC promoter G-quadruplex by a quinoline derivative: insights into molecular targeting of parallel G-quadruplexes

DNA G-Quadruplexes (G4s) formed in oncogene promoters regulate transcription. The oncogene MYC promoter G4 (MycG4) is the most prevalent G4 in human cancers. However, the most studied MycG4 sequence bears a mutated 3′-residue crucial for ligand recognition. Here, we report a new drug-like small molecule PEQ without a large aromatic moiety that specifically binds MycG4. We determined the NMR solution structures of the wild-type MycG4 and its 2:1 PEQ complex, as well as the structure of the 2:1 PEQ complex of the widely used mutant MycG4. Comparison of the two complex structures demonstrates specific molecular recognition of MycG4 and shows the clear effect of the critical 3′-mutation on the drug binding interface. We performed a systematic analysis of the four available complex structures involving the same mutant MycG4, which can be considered a model system for parallel G4s, and revealed for the first time that the flexible flanking residues are recruited in a conserved and sequence-specific way, as well as unused potential for selective ligand-G4 hydrogen-bond interactions. Our results provide the true molecular basis for MycG4-targeting drugs and new critical insights into future rational design of drugs targeting MycG4 and parallel G4s that are prevalent in promoter and RNA G4s.

A joint method for the screening of pharmacological chaperones for phenylalanine hydroxylase

Phenylalanine hydroxylase (PAH) deficiency (PAHD) is an autosomal recessive disorder that causes severe injury to the nervous system, the treatment of which mainly depends on dietary therapy. The limited treatment options for PAHD are an incentive to develop new methods to identify more efficient therapeutic drugs, such as agonists which could improve PAH activity. In this study, we aimed to establish a rapid and convenient method for the screening and verification of PAH agonists. We compared fluorospectrophotometry and tandem mass spectrometry for detection of enzymatic formation of tyrosine, finding that the latter was a more sensitive method. We optimized immunoprecipitation purification conditions and measurement conditions of PAH activity. The optimal ratio between PAH protein and magnetic beads was 500 μg protein per 20 μL beads, and the optimized conditions for the detection of PAH enzymatic activity included the presence of 75 μM coenzyme ((6R)-l-erythro-5,6,7,8-tetrahydrobiopterin) and 30 min reaction time. Based on virtual screening, we screened ten candidate agonists from the FDA drug library. Three of these (nefopam, fluocinonide, and risperidone) were found to activate the enzyme in a dose-dependent manner (0.1-10 μM) by the joint method. We tested the efficacy of the three agonists on three PAH mutations (p.I65T, p.H107R, and p.D101N) that influence enzyme activity, and found that risperidone could specifically activate D101N-mutated enzyme. In conclusion, we established a joint method that is highly reliable, cost-effective, labor-saving, and time-saving. And we also found a specific agonist for D101N-mutated PAH by this joint method which may assist the development of clinical treatment for PAHD patients with different enzyme deficiencies.

Unlocking G-Quadruplexes as Antiviral Targets

Guanine-rich DNA and RNA sequences can fold into noncanonical nucleic acid structures called G-quadruplexes (G4s). Since the discovery that these structures may act as scaffolds for the binding of specific ligands, G4s aroused the attention of a growing number of scientists. The versatile roles of G4 structures in viral replication, transcription, and translation suggest direct applications in therapy or diagnostics. G4-interacting molecules (proteins or small molecules) may also affect the balance between latent and lytic phases, and increasing evidence reveals that G4s are implicated in generally suppressing viral processes, such as replication, transcription, translation, or reverse transcription. In this review, we focus on the discovery of G4s in viruses and the role of G4 ligands in the antiviral drug discovery process. After assessing the role of viral G4s, we argue that host G4s participate in immune modulation, viral tumorigenesis, cellular pathways involved in virus maturation, and DNA integration of viral genomes, which can be potentially employed for antiviral therapeutics. Furthermore, we scrutinize the impediments and shortcomings in the process of studying G4 ligands and drug discovery. Finally, some unanswered questions regarding viral G4s are highlighted for prospective future projects. SIGNIFICANCE STATEMENT: G-quadruplexes (G4s) are noncanonical nucleic acid structures that have gained increasing recognition during the last few decades. First identified as relevant targets in oncology, their importance in virology is now increasingly clear. A number of G-quadruplex ligands are known: viral transcription and replication are the main targets of these ligands. Both viral and cellular G4s may be targeted; this review embraces the different aspects of G-quadruplexes in both host and viral contexts.

Targeting MYC: From understanding its biology to drug discovery

The MYC oncogene is considered to be a high priority target for clinical intervention in cancer patients due to its aberrant activation in more than 50% of human cancers. Direct small molecule inhibition of MYC has traditionally been hampered by its intrinsically disordered nature and lack of both binding site and enzymatic activity. In recent years, however, a number of strategies for indirectly targeting MYC have emerged, guided by the advent of protein structural information and the growing set of computational tools that can be used to accelerate the hit to lead process in medicinal chemistry. In this review, we provide an overview of small molecules developed for clinical applications of these strategies, which include stabilization of the MYC guanine quadruplex, inhibition of BET factor BRD4, and disruption of the MYC:MAX heterodimer. The recent identification of novel targets for indirect MYC inhibition at the protein level is also discussed.

Properties and biological impact of RNA G-quadruplexes: from order to turmoil and back

Guanine-quadruplexes (G4s) are non-canonical four-stranded structures that can be formed in guanine (G) rich nucleic acid sequences. A great number of G-rich sequences capable of forming G4 structures have been described based on in vitro analysis, and evidence supporting their formation in live cells continues to accumulate. While formation of DNA G4s (dG4s) within chromatin in vivo has been supported by different chemical, imaging and genomic approaches, formation of RNA G4s (rG4s) in vivo remains a matter of discussion. Recent data support the dynamic nature of G4 formation in the transcriptome. Such dynamic fluctuation of rG4 folding-unfolding underpins the biological significance of these structures in the regulation of RNA metabolism. Moreover, rG4-mediated functions may ultimately be connected to mechanisms underlying disease pathologies and, potentially, provide novel options for therapeutics. In this framework, we will review the landscape of rG4s within the transcriptome, focus on their potential impact on biological processes, and consider an emerging connection of these functions in human health and disease.

Enhanced sampling molecular dynamics simulations correctly predict the diverse activities of a series of stiff-stilbene G-quadruplex DNA ligands

Ligands with the capability to bind G-quadruplexes (G4s) specifically, and to control G4 structure and behaviour, offer great potential in the development of novel therapies, technologies and functional materials. Most known ligands bind to a pre-formed topology, but G4s are highly dynamic and a small number of ligands have been discovered that influence these folding equilibria. Such ligands may be useful as probes to understand the dynamic nature of G4 in vivo, or to exploit the polymorphism of G4 in the development of molecular devices. To date, these fascinating molecules have been discovered serendipitously. There is a need for tools to predict such effects to drive ligand design and development, and for molecular-level understanding of ligand binding mechanisms and associated topological perturbation of G4 structures. Here we study the G4 binding mechanisms of a family of stiff-stilbene G4 ligands to human telomeric DNA using molecular dynamics (MD) and enhanced sampling (metadynamics) MD simulations. The simulations predict a variety of binding mechanisms and effects on G4 structure for the different ligands in the series. In parallel, we characterize the binding of the ligands to the G4 target experimentally using NMR and CD spectroscopy. The results show good agreement between the simulated and experimentally observed binding modes, binding affinities and ligand-induced perturbation of the G4 structure. The simulations correctly predict ligands that perturb G4 topology. Metadynamics simulations are shown to be a powerful tool to aid development of molecules to influence G4 structure, both in interpreting experiments and to help in the design of these chemotypes.

Synthesis and investigation of quadruplex-DNA-binding, 9-O-substituted berberine derivatives

A small series of five novel berberine derivatives was synthesized by the Cu-catalyzed click reaction of 9-propargyladenine with 9-O-(azidoalkyl)berberine derivatives. The association of the resulting berberine-adenine conjugates with representative quadruplex-forming oligonucleotides 22AG dA(G3TTA)3G3 and a2 d(ACAG4TGTG4)2 was examined with photometric and fluorimetric titrations, thermal DNA denaturation analysis, and CD spectroscopy. The results from the spectrometric titrations indicated the formation of 2:1 or 1:1 complexes (ligand:G4-DNA) with log K b values of 10-11 (2:1) and 5-6 (1:1), which are typical for berberine derivatives. Notably, a clear relationship between the binding affinity of the ligands with the length of the alkyl linker chain, n, was not observed. However, depending on the structure, the ligands exhibited different effects when bound to the G4-DNA, such as fluorescent light-up effects and formation of ICD bands, which are mostly pronounced with a linker length of n = 4 (with a2) and n = 5 (with 22AG), thus indicating that each ligand-G4-DNA complex has a specific structure with respect to relative alignment and conformational flexibility of the ligand in the binding site. It was shown exemplarily with one representative ligand from the series that such berberine-adenine conjugates exhibit a selective binding, specifically a selectivity to quadruplex DNA in competition with duplex DNA, and a preferential thermal stabilization of the G4-DNA forms 22AG and KRAS. Notably, the experimental data do not provide evidence for a significant effect of the adenine unit on the binding affinity of the ligands, for example, by additional association with the loops, presumably because the adenine residue is sterically shielded by the neighboring triazole unit.

No G-Quadruplex Structures in the DNA of Parvovirus B19: Experimental Evidence versus Bioinformatic Predictions

Parvovirus B19 (B19V), an ssDNA virus in the family Parvoviridae, is a human pathogenic virus, responsible for a wide range of clinical manifestations, still in need of effective and specific antivirals. DNA structures, including G-quadruplex (G4), have been recognised as relevant functional features in viral genomes, and small-molecule ligands binding to these structures are promising antiviral compounds. Bioinformatic tools predict the presence of potential G4 forming sequences (PQSs) in the genome of B19V, raising interest as targets for antiviral strategies. Predictions locate PQSs in the genomic terminal regions, in proximity to replicative origins. The actual propensity of these PQSs to form G4 structures was investigated by circular dichroism spectroscopic analysis on synthetic oligonucleotides of corresponding sequences. No signature of G4 structures was detected, and the interaction with the G4 ligand BRACO-19 (N,N'-(9-{[4-(dimethylamino)phenyl]amino}acridine-3,6-diyl)bis(3-pyrrolidin-1-ylpropanamide) did not appear consistent with the stabilisation of G4 structures. Any potential role of PQSs in the viral lifecycle was then assessed in an in vitro infection model system, by evaluating any variation in replication or expression of B19V in the presence of the G4 ligands BRACO-19 and pyridostatin. Neither showed a significant inhibitory activity on B19V replication or expression. Experimental challenge did not support bioinformatic predictions. The terminal regions of B19V are characterised by relevant sequence and symmetry constraints, which are functional to viral replication. Our experiments suggest that these impose a stringent requirement prevailing over the propensity of forming actual G4 structures.

Probing the Importance of the G-Quadruplex Grooves for the Activity of the Anti-HIV-Integrase Aptamer T30923

In this paper, we report studies concerning four variants of the G-quadruplex forming anti-HIV-integrase aptamer T30923, in which specific 2′-deoxyguanosines have been singly replaced by 8-methyl-2′-deoxyguanosine residues, with the aim to exploit the methyl group positioned in the G-quadruplex grooves as a steric probe to investigate the interaction aptamer/target. Although, the various modified aptamers differ in the localization of the methyl group, NMR, circular dichroism (CD), electrophoretic and molecular modeling data suggest that all of them preserve the ability to fold in a stable dimeric parallel G-quadruplex complex resembling that of their natural counterpart T30923. However, the biological data have shown that the T30923 variants are characterized by different efficiencies in inhibiting the HIV-integrase, thus suggesting the involvement of the G-quadruplex grooves in the aptamer/target interaction.

An amphiphilic BODIPY-based selective probe for parallel G4 DNA targeting via disaggregation-induced emission

An amphiphilic BODIPY-based probe, AB-1, was established for parallel G4 DNA targeting based on the concept of triggered disaggregation-induced emission.

Modulation of DNA structure formation using small molecules

Genome integrity is essential for proper cell function such that genetic instability can result in cellular dysfunction and disease. Mutations in the human genome are not random, and occur more frequently at "hotspot" regions that often co-localize with sequences that have the capacity to adopt alternative (i.e. non-B) DNA structures. Non-B DNA-forming sequences are mutagenic, can stimulate the formation of DNA double-strand breaks, and are highly enriched at mutation hotspots in human cancer genomes. Thus, small molecules that can modulate the conformations of these structure-forming sequences may prove beneficial in the prevention and/or treatment of genetic diseases. Further, the development of molecular probes to interrogate the roles of non-B DNA structures in modulating DNA function, such as genetic instability in cancer etiology are warranted. Here, we discuss reported non-B DNA stabilizers, destabilizers, and probes, recent assays to identify ligands, and the potential biological applications of these DNA structure-modulating molecules.

Structural flexibility versus rigidity of the aromatic unit of DNA ligands: binding of aza- and azoniastilbene derivatives to duplex and quadruplex DNA

The known azastilbene (E)-1,2-di(quinolin-3-yl)ethane (2a) and the novel azoniastilbene derivatives (E)-2-(2-(naphthalen-2-yl)vinyl)quinolizinium (2b) and (E)-3,3'-(ethane-1,2-diyl)bis(1-methylquinolinin-1-ium) (2c) were synthesized. Their interactions with duplex and quadruplex DNA (G4-DNA) were studied by photometric, fluorimetric, polarimetric and flow-LD analysis, and by thermal DNA denaturation studies, as well as by ¹H-NMR spectroscopy. The main goal of this study was a comparison of these conformationally flexible compounds with the known G4-DNA-binding diazoniadibenzo[b,k]chrysenes, that have a comparable π-system extent, but a rigid structure. We have observed that the aza- and azoniastilbene derivatives 2a-c, i.e. compounds with almost the same spatial dimensions and steric demand, bind to DNA with an affinity and selectivity that depends significantly on the number of positive charges. Whereas the charge neutral derivative 2a binds unspecifically to the DNA backbone of duplex DNA, the ionic compounds 2b and 2c are typical DNA intercalators. Notably, the bis-quinolinium derivative 2c binds to G4-DNA with moderate affinity (K_b = 4.8 × 10⁵ M^-1) and also stabilizes the G4-DNA towards thermal denaturation (ΔT_m = 11 °C at ligand-DNA ratio = 5.0). Strikingly, the corresponding rigid counterpart, 4a,12a-diazonia-8,16-dimethyldibenzo[b,k]chrysene, stabilizes the G4-DNA to an even greater extent under identical conditions (ΔT_m = 27 °C). These results indicate that the increased flexibility of a G4-DNA ligand does not necessarily lead to stronger interactions with the G4-DNA as compared with rigid ligands that have essentially the same size and π system extent.

Exploration of the Structure and Recognition of a G-quadruplex in the her2 Proto-oncogene Promoter and Its Transcriptional Regulation

G-quadruplexes in oncogene promoters provide putative targets for transcriptional regulation. The structure of a putative G-quadruplex sequence (S1: GGAGAAGGAGGAGGTGGAGGAGGAGGG) in potassium solution in the her2 promoter has been resolved mainly through nuclear magnetic resonance (NMR) spectroscopy. By application of various NMR spectra, we proved the formation of a four-layer G-quadruplex composing of two G-tetrads and two G/A-mixed planes with a four-residues loop (A3-G4-A5-A6). Further evidence from a luciferase reporter assay, Q-RT-PCR and Western blotting indicates that S1 G-quadruplex formation can repress her2 promoter activity, and a selected G-quadruplex ligand cβ can enhance the repression by down regulating her2 transcription and expression. These findings provide a G-quadruplex target and perspective implications in her2 transcriptional regulation.

Therapeutic Inhibition of Myc in Cancer. Structural Bases and Computer-Aided Drug Discovery Approaches

Myc (avian myelocytomatosis viral oncogene homolog) represents one of the most sought after drug targets in cancer. Myc transcription factor is an essential regulator of cell growth, but in most cancers it is overexpressed and associated with treatment-resistance and lethal outcomes. Over 40 years of research and drug development efforts did not yield a clinically useful Myc inhibitor. Drugging the "undruggable" is problematic, as Myc inactivation may negatively impact its physiological functions. Moreover, Myc is a disordered protein that lacks effective binding pockets on its surface. It is well established that the Myc function is dependent on dimerization with its obligate partner, Max (Myc associated factor X), which together form a functional DNA-binding domain to activate genomic targets. Herein, we provide an overview of the knowledge accumulated to date on Myc regulation and function, its critical role in cancer, and summarize various strategies that are employed to tackle Myc-driven malignant transformation. We focus on important structure-function relationships of Myc with its interactome, elaborating structural determinants of Myc-Max dimer formation and DNA recognition exploited for therapeutic inhibition. Chronological development of small-molecule Myc-Max prototype inhibitors and corresponding binding sites are comprehensively reviewed and particular emphasis is placed on modern computational drug design methods. On the outlook, technological advancements may soon provide the so long-awaited Myc-Max clinical candidate.

Binding Study of the Fluorescent Carbazole Derivative with Human Telomeric G-Quadruplexes

The carbazole ligand 3 was synthesized, characterized and its binding interactions with human telomeric (22HT) G-quadruplex DNA in Na⁺ and K⁺-containing buffer were investigated by ultraviolet-visible (UV-Vis) spectrophotometry, fluorescence, circular dichroism (CD) spectroscopy, and DNA melting. The results showed that the studied carbazole ligand interacted and stabilized the intramolecular G-quadruplexes formed by the telomeric sequence in the presence of sodium and potassium ions. In the UV-Vis titration experiments a two-step complex formation between ligand and G-quadruplex was observed. Very low fluorescence intensity of the carbazole derivative in Tris HCl buffer in the presence of the NaCl or KCl increased significantly after addition of the 22HT G4 DNA. Binding stoichiometry of the ligand/G-quadruplex was investigated with absorbance-based Job plots. Carbazole ligand binds 22HT with about 2:1 stoichiometry in the presence of sodium and potassium ions. The binding mode appeared to be end-stacking with comparable binding constants of ~10⁵ M⁻¹ as determined from UV-Vis and fluorescence titrations data. The carbazole ligand is able to induce formation of G4 structure of 22HT in the absence of salt, which was proved by CD spectroscopy and melting studies. The derivative of carbazole 3 shows significantly higher cytotoxicity against breast cancer cells then for non-tumorigenic breast epithelial cells. The cytotoxic activity of ligand seems to be not associated with telomerase inhibition.