Binding modes and pathway of RHPS4 to human telomeric G-quadruplex and duplex DNA probed by all-atom molecular dynamics simulations with explicit solvent

Determination of ligand selectivity to G-tetrad through an AMCA fluorescence quenching approach

The selective binding of ligand molecules towards the 5' and 3' ends of G-quadruplex (G4) may differentially affect the physiological function of G4s. However, there is still a lack of sensitive and low-cost approaches to accurately measure the binding preference of ligands on G4s, although multiple ways have been developed to evaluate the interaction between ligands and G4s. Here, we propose a new protocol named G4-AFQ to test the selectivity of ligands towards the two terminal G-tetrads of G4s. In this protocol, the fluorophore AMCA is respectively modified at the 5' or 3' end of G4, and which end of AMCA fluorescence is quenched means that the ligand binds to the G-tetrad at that end. Through G4-AFQ, the affinity constant of ligands towards the binding site can also be obtained. Compared with the commonly used nuclear magnetic resonance (NMR) method, G4-AFQ is more convenient, sensitive, cost-effective, and suitable for the measurement of the vast majority of G4 ligands, with a great potential for widespread application.

From Crypts to Cancer: A Holistic Perspective on Colorectal Carcinogenesis and Therapeutic Strategies

Colorectal cancer (CRC) represents a significant global health burden, with high incidence and mortality rates worldwide. Recent progress in research highlights the distinct clinical and molecular characteristics of colon versus rectal cancers, underscoring tumor location's importance in treatment approaches. This article provides a comprehensive review of our current understanding of CRC epidemiology, risk factors, molecular pathogenesis, and management strategies. We also present the intricate cellular architecture of colonic crypts and their roles in intestinal homeostasis. Colorectal carcinogenesis multistep processes are also described, covering the conventional adenoma-carcinoma sequence, alternative serrated pathways, and the influential Vogelstein model, which proposes sequential APC, KRAS, and TP53 alterations as drivers. The consensus molecular CRC subtypes (CMS1-CMS4) are examined, shedding light on disease heterogeneity and personalized therapy implications.

G-quadruplex ligands as therapeutic agents against cancer, neurological disorders and viral infections

G-quadruplexes (G4s) within the human genome have undergone extensive molecular investigation, with a strong focus on telomeres, gene promoters and repetitive regulatory sequences. G4s play central roles in regulating essential biological processes, including telomere maintenance, replication, transcription and translation. Targeting these molecular processes with G4-binding ligands holds substantial therapeutic potential in anticancer treatments and has also shown promise in treating neurological, skeletal and muscular disorders. The presence of G4s in bacterial and viral genomes also suggests that G4-binding ligands could be a critical tool in fighting infections. This review provides an overview of the progress and applications of G4-binding ligands, their proposed mechanisms of action, challenges faced and prospects for their utilization in anticancer treatments, neurological disorders and antiviral activities.

Targeting telomeres: advances in telomere maintenance mechanism-specific cancer therapies

Cancer cells establish replicative immortality by activating a telomere-maintenance mechanism (TMM), be it telomerase or the alternative lengthening of telomeres (ALT) pathway. Targeting telomere maintenance represents an intriguing opportunity to treat the vast majority of all cancer types. Whilst telomerase inhibitors have historically been heralded as promising anticancer agents, the reality has been more challenging, and there are currently no therapeutic options for cancer types that use ALT despite their aggressive nature and poor prognosis. In this Review, we discuss the mechanistic differences between telomere maintenance by telomerase and ALT, the current methods used to detect each mechanism, the utility of these tests for clinical diagnosis, and recent developments in the therapeutic strategies being employed to target both telomerase and ALT. We present notable developments in repurposing established therapeutic agents and new avenues that are emerging to target cancer types according to which TMM they employ. These opportunities extend beyond inhibition of telomere maintenance, by finding and exploiting inherent weaknesses in the telomeres themselves to trigger rapid cellular effects that lead to cell death.

Exploring the Interaction of G-quadruplex Binders with a (3 + 1) Hybrid G-quadruplex Forming Sequence within the PARP1 Gene Promoter Region

The enzyme PARP1 is an attractive target for cancer therapy, as it is involved in DNA repair processes. Several PARP1 inhibitors have been approved for clinical treatments. However, the rapid outbreak of resistance is seriously threatening the efficacy of these compounds, and alternative strategies are required to selectively regulate PARP1 activity. A noncanonical G-quadruplex-forming sequence within the PARP1 promoter was recently identified. In this study, we explore the interaction of known G-quadruplex binders with the G-quadruplex structure found in the PARP gene promoter region. The results obtained by NMR, CD, and fluorescence titration, also confirmed by molecular modeling studies, demonstrate a variety of different binding modes with small stabilization of the G-quadruplex sequence located at the PARP1 promoter. Surprisingly, only pyridostatin produces a strong stabilization of the G-quadruplex-forming sequence. This evidence makes the identification of a proper (3+1) stabilizing ligand a challenging goal for further investigation.

Studies on the Interactions of 3,11-Difluoro-6,8,13-trimethyl-8H-quino[4,3,2-kl]acridinium and Insulin with the Quadruplex-Forming Oligonucleotide Sequence a2 from the Insulin-Linked Polymorphic Region

Recently, several quadruplex-DNA-forming sequences have been identified in the insulin-linked polymorphic region (ILPR), which is a guanine-rich oligonucleotide sequence in the promoter region of insulin. The formation of this non-canonical quadruplex DNA (G4-DNA) has been shown to be involved in the biological activity of the ILPR, specifically with regard to its interplay with insulin. In this context, this contribution reports on the investigation of the association of the quadruplex-forming ILPR sequence a2 with insulin as well as with the well-known G4-DNA ligand 3,11-difluoro-6,8,13-trimethyl-8H-quino[4,3,2-kl]acridinium (1), also named RHPS4, by optical and NMR spectroscopy. CD- and NMR-spectroscopic measurements confirmed the preferential formation of an antiparallel quadruplex structure of a2 with four stacked guanine quartets. Furthermore, ligand 1 has high affinity toward a2 and binds by terminal π stacking to the G1-G11-G15-G25 quartet. In addition, the spectroscopic studies pointed to an association of insulin to the deoxyribose backbone of the loops of a2.

G-Quadruplexes and Their Ligands: Biophysical Methods to Unravel G-Quadruplex/Ligand Interactions

Progress in the design of G-quadruplex (G4) binding ligands relies on the availability of approaches that assess the binding mode and nature of the interactions between G4 forming sequences and their putative ligands. The experimental approaches used to characterize G4/ligand interactions can be categorized into structure-based methods (circular dichroism (CD), nuclear magnetic resonance (NMR) spectroscopy and X-ray crystallography), affinity and apparent affinity-based methods (surface plasmon resonance (SPR), isothermal titration calorimetry (ITC) and mass spectrometry (MS)), and high-throughput methods (fluorescence resonance energy transfer (FRET)-melting, G4-fluorescent intercalator displacement assay (G4-FID), affinity chromatography and microarrays. Each method has unique advantages and drawbacks, which makes it essential to select the ideal strategies for the biological question being addressed. The structural- and affinity and apparent affinity-based methods are in several cases complex and/or time-consuming and can be combined with fast and cheap high-throughput approaches to improve the design and development of new potential G4 ligands. In recent years, the joint use of these techniques permitted the discovery of a huge number of G4 ligands investigated for diagnostic and therapeutic purposes. Overall, this review article highlights in detail the most commonly used approaches to characterize the G4/ligand interactions, as well as the applications and types of information that can be obtained from the use of each technique.

Spectroscopic and In Silico Studies on the Interaction of Substituted Pyrazolo[1,2-a]benzo[1,2,3,4]tetrazine-3-one Derivatives with c-Myc G4-DNA

Herein we describe a combined experimental and in silico study of the interaction of a series of pyrazolo[1,2-a]benzo[1,2,3,4]tetrazin-3-one derivatives (PBTs) with parallel G-quadruplex (GQ) DNA aimed at correlating their previously reported anticancer activities and the stabilizing effects observed by us on c-myc oncogene promoter GQ structure. Circular dichroism (CD) melting experiments were performed to characterize the effect of the studied PBTs on the GQ thermal stability. CD measurements indicate that two out of the eight compounds under investigation induced a slight stabilizing effect (2-4 °C) on GQ depending on the nature and position of the substituents. Molecular docking results allowed us to verify the modes of interaction of the ligands with the GQ and estimate the binding affinities. The highest binding affinity was observed for ligands with the experimental melting temperatures (Tms). However, both stabilizing and destabilizing ligands showed similar scores, whilst Molecular Dynamics (MD) simulations, performed across a wide range of temperatures on the GQ in water solution, either unliganded or complexed with two model PBT ligands with the opposite effect on the Tms, consistently confirmed their stabilizing or destabilizing ability ascertained by CD. Clues about a relation between the reported anticancer activity of some PBTs and their ability to stabilize the GQ structure of c-myc emerged from our study. Furthermore, Molecular Dynamics simulations at high temperatures are herein proposed for the first time as a means to verify the stabilizing or destabilizing effect of ligands on the GQ, also disclosing predictive potential in GQ-targeting drug discovery.

An in silico investigation of the binding modes and pathway of APTO-253 on c-KIT G-quadruplex DNA

The stability of c-KIT G-quadruplex DNA via ligands has been a significant concern in the growing field of cancer therapy. Thus, it is very important to understand the mechanism behind the high binding affinity of the small drug molecules on the c-KIT G-quadruplex DNA. In this study, we have investigated the binding mode and pathway of the APTO-253 ligand on the c-KIT G-quadruplex DNA employing a total of 10 μs all atom molecular dynamics simulations and further 8.82 μs simulations via the umbrella sampling method using both OL15 and BSC1 latest force fields for DNA structures. From the cluster structure analysis, mainly three binding pathways i.e., top, bottom and side loop stacking modes are identified. Moreover, RMSD, RMSF and 2D-RMSD values indicate that the c-KIT G-quadruplex DNA and APTO-253 molecules are stable throughout the simulation run. Furthermore, the number of hydrogen bonds in each tetrad and the distance between the two central K⁺ cations confirm that the c-KIT G-quadruplex DNA maintains its conformation in the process of complex formation with the APTO-253 ligand. The binding free energies and the minimum values in the potential of mean forces suggest that the binding processes are energetically favorable. Furthermore, we have found that the bottom stacking mode is the most favorable binding mode among all the three modes for the OL15 force field. However, for the BSC1 force field, both the top and bottom binding modes of the APTO-253 ligand in c-KIT G-quadruplex DNA are comparable to each other. To investigate the driving force for the complex formation, we have noticed that the van der Waals (vdW) and π-π stacking interactions are mainly responsible. Our detailed studies provide useful information for the discovery of novel drugs in the field of stabilization of G-quadruplex DNAs.

Insights into G-Quadruplex-Hemin Dynamics Using Atomistic Simulations: Implications for Reactivity and Folding

Guanine quadruplex nucleic acids (G4s) are involved in key biological processes such as replication or transcription. Beyond their biological relevance, G4s find applications as biotechnological tools since they readily bind hemin and enhance its peroxidase activity, creating a G4-DNAzyme. The biocatalytic properties of G4-DNAzymes have been thoroughly studied and used for biosensing purposes. Despite hundreds of applications and massive experimental efforts, the atomistic details of the reaction mechanism remain unclear. To help select between the different hypotheses currently under investigation, we use extended explicit-solvent molecular dynamics (MD) simulations to scrutinize the G4/hemin interaction. We find that besides the dominant conformation in which hemin is stacked atop the external G-quartets, hemin can also transiently bind to the loops and be brought to the external G-quartets through diverse delivery mechanisms. The simulations do not support the catalytic mechanism relying on a wobbling guanine. Similarly, the catalytic role of the iron-bound water molecule is not in line with our results; however, given the simulation limitations, this observation should be considered with some caution. The simulations rather suggest tentative mechanisms in which the external G-quartet itself could be responsible for the unique H₂O₂-promoted biocatalytic properties of the G4/hemin complexes. Once stacked atop a terminal G-quartet, hemin rotates about its vertical axis while readily sampling shifted geometries where the iron transiently contacts oxygen atoms of the adjacent G-quartet. This dynamics is not apparent from the ensemble-averaged structure. We also visualize transient interactions between the stacked hemin and the G4 loops. Finally, we investigated interactions between hemin and on-pathway folding intermediates of the parallel-stranded G4 fold. The simulations suggest that hemin drives the folding of parallel-stranded G4s from slip-stranded intermediates, acting as a G4 chaperone. Limitations of the MD technique are briefly discussed.

Enhanced Rigidity Changes Ultraviolet Absorption: Effect of a Merocyanine Binder on G-Quadruplex Photophysics

The urge to discover selective fluorescent binders to G-quadruplexes (G4s) for rapid diagnosis must be linked to understand the effect that those have on the DNA photophysics. Herein, we report on the electronic excited states of a bound merocyanine dye to c-Myc G4 using extensive multiscale quantum mechanics/molecular mechanics calculations. We find that the absorption spectra of c-Myc G4, both without and with the intercalated dye, are mainly composed of exciton states and mixed local/charge-transfer states. The presence of merocyanine hardly affects the energy range of the guanine absorption or the number of guanines excited. However, it triggers a substantial amount (16%) of detrimental pure charge-transfer states involving oxidized guanines. We identify the rigidity introduced by the probe in G4, reducing the overlap among guanines, as the one responsible for the changes in the exciton and charge-transfer states, ultimately leading to a redshift of the absorption maximum.

Molecular Dynamics Study on the Binding of an Anticancer DNA G-Quadruplex Stabilizer, CX-5461, to Human Telomeric, c-KIT1, and c-Myc G-Quadruplexes and a DNA Duplex

DNA G-quadruplex (G4) stabilizer, CX-5461, is in phase I/II clinical trials for advanced cancers with BRCA1/2 deficiencies. A FRET-melting temperature increase assay measured the stabilizing effects of CX-5461 to a DNA duplex (∼10 K), and three G4 forming sequences negatively implicated in the cancers upon its binding: human telomeric (∼30 K), c-KIT1 (∼27 K), and c-Myc (∼25 K). Without experimentally solved structures of these CX-5461-G4 complexes, CX-5461's interactions remain elusive. In this study, we performed a total of 73.5 μs free ligand molecular dynamics binding simulations of CX-5461 to the DNA duplex and three G4s. Three binding modes (top, bottom, and side) were identified for each system and their thermodynamic, kinetic, and structural nature were deciphered. The molecular mechanics/Poisson Boltzmann surface area binding energies of CX-5461 were calculated for the human telomeric (-28.6 kcal/mol), c-KIT1 (-23.9 kcal/mol), c-Myc (-22.0 kcal/mol) G4s, and DNA duplex (-15.0 kcal/mol) systems. These energetic differences coupled with structural differences at the 3' site explained the different melting temperatures between the G4s, while CX-5461's lack of intercalation to the duplex explained the difference between the G4s and duplex. Based on the interaction insight, CX-5461 derivatives were designed and docked, showing higher selectivity to the G4s over the duplex.

Synthesis of naphthalimide derivatives with potential anticancer activity, their comparative ds- and G-quadruplex-DNA binding studies and related biological activities

Two new cytotoxic 1,8-naphthalimide derivatives have been synthesized and characterized. Their biological activities as cytotoxicity and antimicrobial activities and inhibitory activities against DNA-polymerase were evaluated. The interactions of compounds with double-stranded- and quadruple-DNA have been studied by UV-Vis, fluorescent intercalator displacement, competition dialysis, circular dichroism and the findings were compared with the parent naphthalimide and the other compounds. The results show that both compounds (1 and 2) and the parent compound NI have strong cytotoxic activities against Beas-2B, MCF-7, HepG2 and MDA-MB-231 cancer cell lines, antimicrobial activities against Staphylococcus aureus ATCC 29213, Enterococcus faecalis ATCC 29212 and inhibitory activities towards Taq-polymerase and transcriptase. These novel cationic compounds 1 and 2 can stabilize G-quadruplexes DNA according to thermal denaturation experiments, they change the 3D structure of the DNA (see details in CD experiments) and they exhibit different binding affinities for q-DNA and ds-DNA revealed by spectrophotometric titrations and competitive dialysis studies.

Three-Dimensional Structure of RNA Monomeric G-Quadruplex Containing ALS and FTD Related G4C2 Repeat and Its Binding with TMPyP4 Probed by Homology Modeling based on Experimental Constraints and Molecular Dynamics Simulations

The G-quadruplex-forming hexanucleotide repeat expansion (HRE), d(G4C2)_n, within the human C9orf72 gene is the root cause for familial amyotrophic lateral sclerosis-frontotemporal dementia (ALS-FTD). A recent study has shown that TMPyP4 has good potential to work as a RNA G-quadruplex binder in treating ALS and FTD. Although the high-resolution structure of the monomeric DNA antiparallel G-quadruplex form of the monomeric hexanucleotide repeat was recently solved, the RNA parallel G-quadruplex structure and its complex with TMPyP4 are not available yet. In this study, we first constructed the homology model for the parallel monomeric RNA G-quadruplex of r(G4C2)₃G4 based on experimental constraints and the parallel monomeric G-quadruplex DNA crystal structure. Although the G-tetra core of the homology model was stable observed in 15 μs molecular dynamics (MD) simulations, we observed that the loops adopt additional conformations besides the initial crystal conformation, where TMPyP4 binding was found to reduce the loop fluctuation of the RNA monomeric G-quadruplex. Next, we probed the elusive binding behavior of TMPyP4 to the RNA monomeric G-quadruplex. Encouragingly, the binding modes observed are similar to the modes observed in two experimental complexes of a parallel DNA G-quadruplex with TMPyP4. We also constructed a Markov state model to provide insights into the binding pathways. Together, the findings from our study may assist future development of G-quadruplex-specific ligands in the treatment of neurodegenerative diseases like ALS and FTD.

To probe the binding pathway of a selective compound (D089-0563) to c-MYC Pu24 G-quadruplex using free ligand binding simulations and Markov state model analysis

D089-0563 is a highly promising anti-cancer compound that selectively binds the transcription-silencing G-quadruplex element (Pu27) at the promoter region of the human c-MYC oncogene; however, its binding mechanism remains elusive.

The investigation of the G-quadruplex aptamer selectivity to Pb2+ ion: a joint molecular dynamics simulation and density functional theory study

The aptamers with the ability to form a G-quadruplex structure can be stable in the presence of some ions. Hence, study of the interactions between such aptamers and ions can be beneficial to determine the highest selective aptamer toward an ion. In this article, molecular dynamics (MD) simulations and quantum mechanics (QM) calculations have been applied to investigate the selectivity of the T30695 aptamer toward Pb²⁺ in comparison with some ions. The Free Energy Landscape (FEL) analysis indicates that Pb²⁺ has remained inside the aptamer during the MD simulation, while the other ions have left it. The Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) binding energies prove that the conformational stability of the aptamer is the highest in the presence of Pb²⁺. According to the compaction parameters, the greatest compressed ion-aptamer complex, and hence, the highest ion-aptamer interaction have been induced in the presence of Pb²⁺. The contact maps clarify the closer contacts between the nucleotides of the aptamer in the presence of Pb²⁺. The density functional theory (DFT) results show that Pb²⁺ forms the most stable complex with the aptamer, which is consistent with the MD results. The QM calculations reveal that the N-H bonds and the O…H distances are the longest and the shortest, respectively, in the presence of Pb²⁺. The obtained results verify that the strongest hydrogen bonds (HBs), and hence, the most compressed aptamer structure are induced by Pb²⁺. Besides, atoms in molecules (AIM) and natural bond orbital (NBO) analyses confirm the results.Communicated by Ramaswamy H. Sarma.

Evaluation of the selectivity of G-quadruplex ligands in living cells with a small molecule fluorescent probe

G-quadruplex has been an emerging target for drug design due to its physiologically important roles in oncology. A number of quadruplex-interactive ligands have been developed by synthetic and medicinal chemists over the past decades. However, the great challenge still remains that the method for detecting the specific targeting of these ligands to the G-quadruplex structures in cells is still lacking. Herein, a detection system for directly identifying the specific targeting of a ligand to DNA G-quadruplexes in cells was constructed by using a small-molecular fluorescent probe (IMT) as a fluorescent indicator. Four typical ligands have been successfully evaluated, demonstrating the promising application of this detection system in the screening and evaluation of quadruplex-specific therapeutic agents.

Binding of BRACO19 to a Telomeric G-Quadruplex DNA Probed by All-Atom Molecular Dynamics Simulations with Explicit Solvent

Although BRACO19 is a potent G-quadruplex binder, its potential for clinical usage is hindered by its low selectivity towards DNA G-quadruplex over duplex. High-resolution structures of BRACO19 in complex with neither single-stranded telomeric DNA G-quadruplexes nor B-DNA duplex are available. In this study, the binding pathway of BRACO19 was probed by 27.5 µs molecular dynamics binding simulations with a free ligand (BRACO19) to a DNA duplex and three different topological folds of the human telomeric DNA G-quadruplex (parallel, anti-parallel and hybrid). The most stable binding modes were identified as end stacking and groove binding for the DNA G-quadruplexes and duplex, respectively. Among the three G-quadruplex topologies, the MM-GBSA binding energy analysis suggested that BRACO19′s binding to the parallel scaffold was most energetically favorable. The two lines of conflicting evidence plus our binding energy data suggest conformation-selection mechanism: the relative population shift of three scaffolds upon BRACO19 binding (i.e., an increase of population of parallel scaffold, a decrease of populations of antiparallel and/or hybrid scaffold). This hypothesis appears to be consistent with the fact that BRACO19 was specifically designed based on the structural requirements of the parallel scaffold and has since proven effective against a variety of cancer cell lines as well as toward a number of scaffolds. In addition, this binding mode is only slightly more favorable than BRACO19s binding to the duplex, explaining the low binding selectivity of BRACO19 to G-quadruplexes over duplex DNA. Our detailed analysis suggests that BRACO19′s groove binding mode may not be stable enough to maintain a prolonged binding event and that the groove binding mode may function as an intermediate state preceding a more energetically favorable end stacking pose; base flipping played an important role in enhancing binding interactions, an integral feature of an induced fit binding mechanism.

Tuning the ring strain effect in acridine derivatives on binding affinity with G-quadruplex-DNA: A computational and experimental study

Search for inhibitors to stabilize the telomeric G-quadruplex in order to deter telomerase activity is an active area of research. Inhibitors play an important role to initiate the tumor cell mortalization process. This work reports for the first time of acridine derivative with four membered ammonium rings at the side chain to surpass the binding ability against BRACO-19 with G-quadruplex-DNA. It is known in the literature that acridine based molecule BRACO-19 can effectively bind with G-quadruplex-DNA. The computational study performed in this study revealed that the binding ability of acridine based molecule can be augmented with subtle variation in the molecular structure of the drug like candidates. Steered molecular dynamics (SMD) performed with the acridine derivatives and G-quadruplex DNA showed the importance of ring strain to the side chain of those ligand molecules. The rupture force analysis, hydrogen bonding interactions and the calculated free energies in MM-PBSA method suggest that ligand 3 is superior than that of BRACO-19. The synthesized ligand 3 and BRACO-19 showed the binding constants obtained from ITC measurements are 4 × 10⁶ mol^-1 and 2.6 × 10⁶, which corroborates the computational findings.

Binding of Telomestatin, TMPyP4, BSU6037, and BRACO19 to a Telomeric G-Quadruplex-Duplex Hybrid Probed by All-Atom Molecular Dynamics Simulations with Explicit Solvent

A promising anticancer therapeutic strategy is the stabilization of telomeric G-quadruplexes using G-quadruplex-binding small molecules. Although many G-quadruplex-specific ligands have been developed, their low potency and selectivity to G-quadruplexes over duplex remains unsolved. Recently, a crystal structure of a telomeric 3' quadruplex-duplex hybrid was reported and the quadruplex-duplex interface was suggested to a good target to address the issues. However, there are no high-resolution complex structures reported for G-quadruplex ligands except for a docked BSU6037. In this study, molecular dynamic (MD) binding simulations with a free ligand were used to study binding poses and dynamics of four representative ligands: telomestatin, TMPyP4, BSU6037, and BRACO19. The MD data showed that BSU6037 was able to fully intercalate into the interface whereas TMPyP4 and BRACO19 could only maintain partial intercalation into the interface and telomestatin only binds at the quadruplex and duplex ends. Both linear ligands, BSU6037 and BRACO19, were able to interact with the interface, yet they were not selective over duplex DNA. The DNA geometry, binding modes, and binding pathways were systematically characterized, and the binding energy was calculated and compared for each system. The interaction of the ligands to the interface was by the means of an induced-fit binding mechanism rather than a lock-key mechanism, consisting of the DNA unfolding at the interface to allow entrance of the drug and then the refolding and repacking of the DNA and the ligand to further stabilize the G-quadruplex. On the basis of the findings in this study, modifications were suggested to optimize the interface binding for TMPyp4 and telomestatin.

Probing the Binding Pathway of BRACO19 to a Parallel-Stranded Human Telomeric G-Quadruplex Using Molecular Dynamics Binding Simulation with AMBER DNA OL15 and Ligand GAFF2 Force Fields

Human telomeric DNA G-quadruplex has been identified as a good therapeutic target in cancer treatment. G-quadruplex-specific ligands that stabilize the G-quadruplex have great potential to be developed as anticancer agents. Two crystal structures (an apo form of parallel stranded human telomeric G-quadruplex and its holo form in complex with BRACO19, a potent G-quadruplex ligand) have been solved, yet the binding mechanism and pathway remain elusive. In this study, we simulated the binding of a free BRACO19 molecule to the apo form of the G-quadruplex using the latest AMBER DNA (OL15) and ligand (GAFF2) force fields. Three binding modes have been identified: top stacking, bottom intercalation, and groove binding. Bottom intercalation (51% of the population) resembles the bottom binding pose in the complex crystal structure very well. The groove binding mode is less stable than the bottom binding mode and is likely to be an intermediate state leading to the bottom binding mode. A flip-insertion mechanism was observed in the bottom intercalation mode, during which flipping of the bases outward makes space for ligand insertion, after which the bases flip back to increase the stability of the complex. In addition to reproducing the base-flipping behavior for some loop residues upon ligand binding, the direct alignment type of the ATAT-tetrad was observed in our simulations for the first time. These successes provide initial support for using this combination of the OL15 and GAFF2 force fields to study quadruplex-ligand interactions.