Next Article in Journal
Differential Regulation of Gonadotropins as Revealed by Transcriptomes of Distinct LH and FSH Cells of Fish Pituitary
Next Article in Special Issue
Non-Canonical Helical Structure of Nucleic Acids Containing Base-Modified Nucleotides
Previous Article in Journal
Granulomatous Inflammation in ANCA-Associated Vasculitis
Previous Article in Special Issue
Towards Profiling of the G-Quadruplex Targeting Drugs in the Living Human Cells Using NMR Spectroscopy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Interaction of Curaxin CBL0137 with G-Quadruplex DNA Oligomers

1
Department of Food, Environmental and Nutritional Sciences (DEFENS), University of Milan (Università degli Studi di Milano), 20133 Milan, Italy
2
National Institute of Fundamental Studies, Kandy 20000, Sri Lanka
3
Scientia Advice di Roberto Artali, 20832 Desio, Italy
4
Institute for Advanced Chemistry of Catalonia (IQAC), CSIC, Networking Center on Bioengineering, Biomaterials and Nanomedicine (CIBER-BBN), 08034 Barcelona, Spain
5
Department of Chemical Engineering and Analytical Chemistry, University of Barcelona, 08028 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(12), 6476; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22126476
Submission received: 29 April 2021 / Revised: 8 June 2021 / Accepted: 15 June 2021 / Published: 17 June 2021

Abstract

:
Curaxins and especially the second-generation derivative curaxin CBL0137 have important antitumor activities in multiple cancers such as glioblastoma, melanoma and others. Although most of the authors suggest that their mechanism of action comes from the activation of p53 and inactivation of NF-kB by targeting FACT, there is evidence supporting the involvement of DNA binding in their antitumor activity. In this work, the DNA binding properties of curaxin CBL0137 with model quadruplex DNA oligomers were studied by 1H NMR, CD, fluorescence and molecular modeling. We provided molecular details of the interaction of curaxin with two G-quadruplex structures, the single repeat of human telomere d(TTAGGGT)4 and the c-myc promoter Pu22 sequence. We also performed 1H and 31P NMR experiments were also performed in order to investigate the interaction with duplex DNA models. Our data support the hypothesis that the interaction of curaxin with G-quadruplex may provide a novel insight into the DNA-binding properties of CBL0137, and it will be helpful for the design of novel selective DNA-targeting curaxin analogues.

1. Introduction

Curaxins are a small group of substances endowed with anticancer activity [1]. They were identified in a search for non-genotoxic antiproliferative compounds, simultaneously acting on two tumor targets. The first hits found were quinacrine and other antimalarials, but SAR studies led rapidly to more active compounds, represented by curaxins, a class of substituted carbazoles with electron-withdrawing groups at positions 3 and 6, and an aminoalkyl chain at the nitrogen 9. Curaxins are able to activate the tumor suppressing protein p53 and to suppress the anti-apoptotic nuclear factor NF-κB [2,3]. Importantly, curaxins were found to be more toxic to tumor than to normal cells. Among curaxins, CBL0137 (Figure 1), appeared the most suitable for further development, on the basis of its metabolic stability, solubility and activity in vivo. CBL0137 suppressed tumor growth in xenograft models of colon (DLD-1), renal cell carcinoma (Caki-1) [3], medulloblastoma [4], small-cell lung cancer [5,6], melanoma (Mel-7) and transplanted surgical samples from patients with pancreatic ductal adenocarcinoma [7].
CBL0137 was also found to be active in mixed-lineage leukemia [8] and against cancer stem cells and to potentiate efficacy of gemcitabine in pancreatic cancer [9].
At present, CBL0137 is under phase I clinical trials in patients with hematological malignancies and solid tumors [10].
The cytotoxicity of curaxins, which are associated with the absence of DNA damage, stimulated studies directed to unveil their mechanism of action. Eukaryotic DNA is packed into chromatin, which is a highly ordered complex of DNA and histone proteins. Its basic unit, the nucleosome, consists of an octamer of histones, known as the histone or nucleosome core, which is bound to a 147 base pairs of DNA fragment.
Recent studies suggested that CBL0137 interacts with DNA by modifying the shape of the DNA helix, thus increasing the interbase-pair distance. As a consequence, DNA unwinding and detachment from the histone octamer occur, eventually leading to nucleosome disassembly both in vitro and in cells. The destabilization of the nucleosomes induces the intervention of a histone chaperone, FACT (Facilitates Chromatin Transcription), which binds tightly to chromatin (c-trapping) [2,11,12]. This results in activating phosphorylation of p53 by FACT-associated CK2 and reducing NF-kB–dependent transcription, because of depletion of soluble FACT. Recently, it was reported that such effects are associated with the induction of the interferon response to epigenetic derepression of the cellular “repeatome” [13].
Kantidze et al. [14] found that CBL0137 alters DNA topology leading to the inability of the transcriptional repressor CTCF to bind efficiently to its cognate DNA sites. This effect on CTCF binding results in partial disruption of chromatin loops and in large-scale perturbations in the 3D genome organization.
More recently, Lu et al. [15] tried to clarify how curaxins alter the genomic DNA structure and affect the DNA binding property of key proteins, such as CTCF and FACT. They found that CBL0137 strongly and persistently binds to dsDNA, inducing a huge barrier for DNA unzipping during replication and transcription, thus causing the distinct binding response of CTCF and FACT on DNA.
Further investigations revealed that several pathways, such as inhibition of the self-renewal of cancer stem cells/tumor-initiating cells through NOTCH1 activation and downregulation of heat shock factor 1 (HSF1), thereby increasing tumor cell apoptosis, are involved in the process [5,14,16]. Moreover, CBL037 highly suppresses the expression of c-MYC family genes. Sergeev et al. [17] reported that this curaxin significantly inhibits in vitro DNA methylation by eukaryotic DNA methyltransferase Dnmt3a at low micromolar concentrations. These effects are attributed to the intercalation of CBL0137 into DNA [3,11], although there is no direct evidence for this type of interaction. The change of topology of DNA by binding with CBL0137 was deduced from CD experiments [11], whereas molecular modeling studies showed a possible protruding of the side chains of the carbazole nucleus into the major groove of DNA, with the carbazole N-side chain filling the minor groove [11].
The above reported effects, including inhibition of c-MYC expression, DNA methyltransferase inhibition, and chromatin remodelling, would be consistent with DNA G-quadruplex binding [18,19].
G-quadruplexes are non-canonical nucleic acids secondary structures that may form in G-rich sequences under physiological conditions. Their structural building block is the G-quartet, a planar array of four guanines paired through Hoogsteen bonds. G-quadruplexes play a role in several key cellular processes, including gene transcription, chromatin epigenetics and DNA recombination. G-quadruplex DNA is found in key regulatory regions of the cell such as promoters of proto-oncogenes (c-myc, bcl-2 and c-Kit) [20]. Stabilization of the folded G-quadruplexes due to ligand interactions is proposed to inhibit the binding of transcription factors, leading to downstream silencing of oncogene expression [21,22,23]. Notably, also human telomeric sequences are able to form G-quadruplex structures that are not recognized by telomerase, an enzyme involved in telomere elongation [24].
The planar structure with an aminoalkyl side chain of CBL0137, similar to other non-DNA-damaging G-quadruplex ligands [23], supports the hypothesis that G-quadruplex interaction can have a role in the curaxin activity; however, so far, no evidence has been given for such an interaction.
To confirm our hypothesis, we undertook an investigation of the binding of CBL0137 to G-quadruplex DNA structures of telomeres and promoter oncogenes by exploiting fluorescence spectroscopy, CD, NMR and molecular modelling. We used as models the single repeat sequence of human telomere, d(TTAGGGT)4, and the G to T mutated Pu22 sequence (Pu22T14T23) of the c-myc oncogene, which is overexpressed in a wide range of human tumors.
To further investigate the possible intercalation mode of curaxin into the double helix DNA, as previously reported by some authors (References [3,11]), the study was extended to the self-complementary double helix oligomers d(CGTACG)2 and d(AAGAATTCTT)2 as models for CG and AT-rich sequences, respectively.
Our work provides a novel insight into the DNA-binding properties of CBL0137 that may be relevant in the important anticancer activity of curaxines.

2. Results and Discussion

To increase the solubility of CBL0137 into the aqueous medium, its hydrochloric salt 1 was prepared by treatment of the free amine with HCl in dioxane.

2.1. Interaction of Curaxin with Telomere d(TTAGGGT)4 Quadruplex

The imino proton region of the 1H NMR spectrum of d(TTAGGGT)4 showed three signals between 10 and 12 ppm that are indicative of the formation of a single G-quadruplex species with three G-quartet planes. As curaxin was added to the d(TTAGGGT)4 solution, the NH imino protons moved upfield and G4 and G6 signals became broad even at low ratio R = [curaxin]/[DNA] = 0.25/0.75. At R = 2.0 the imino proton of G4 remained very broad, while the G6 became sharp again. This behavior can be explained by the binding of curaxin to G4 and G6 tetrads, with an intermediate exchange between free and bound state at the level of G4. The greatest variation of chemical shift is observed for G6 signal (Δδ = −0.55 ppm) (Figure 2).
In order to better define the geometry of the complex, a series of 2D NMR experiments were performed. NOESY and TOCSY experiments allowed us to identify the curaxin (Table S1) and d(TTAGGGT)4 protons in the complex (Table S2). Despite of the overlapping of some of curaxin and oligonucleotide signals, several intermolecular NOE interactions were detected. The contacts involved both aromatic and side chain protons of the curaxin with aromatic and ribose protons of d(TTAGGGT)4 at A3G4 and G6T7 sites (Table 1 and Figure 3).
The NMR studies were complemented by a molecular docking simulation, followed by molecular dynamics (MD) optimization. Curaxin was docked at both A3G4 and G6T7 sites (Figure 4A). At the A3G4 binding site the ligand adopts a quite centre-symmetrical location, allowing the formation of π-π stacking interactions with all the bases of the upper A3 and lower G4 tetrads. In this site the complex is stabilized by a hydrogen bond between the charged quaternary nitrogen group of curaxin and N7A3, at a distance of 2.44 Å. One of the two benzene rings of curaxin lies above the K+ ion, resulting in a strong cation-π interaction (4.90 Å). The interaction pattern is completed by two other cation-π interactions formed between the quaternary nitrogen of the ligand and the aromatic component of the A3 and G4 bases (Figure 4B).
Otherwise, at the G6T7 site the ligand does not adopt a center-symmetrical stacking interaction but it is rather shifted towards one half of the G6 tetrad (Figure 4C). In this orientation, the ligand gives π-π stacking interactions with G6 and T7 and is stabilized by an attractive charge interaction between OP2G6 and the charged quaternary nitrogen group of the ligand. The quaternary nitrogen is also involved in a cation-π interaction with the aromatic component of T7.
The best docked conformations of the complexes at the ApG and GpT intercalation sites are in good agreement with the reported NOE contacts (Table 1).

2.2. Interaction of Curaxin with Pu22T14T23 G-Quadruplex

An important and generalized line broadening of the guanine NH imino protons was observed upon the titration of Pu22T14T23 with curaxin, even at low ratio R = 0.25/1.0. This can be due to a strong interaction of the curaxin with the nucleotide, producing an intermediate exchange process between free and bound state on the NMR timescale. In particular, the NH signals belonging to the tetrad G9-G13-G18-G22 disappeared almost completely. At a ratio R > 1.0, the NH signals sharpened and for R = 2.0 only one set of resonance was present (Figure 5). By increasing the R value until a ratio of 4, the spectrum did not change significantly. All the signals moved upfield, but the most relevant chemical shift variation was observed for G7, G11 and G16 residues at 5′-end, and for G22, G18 and G13 at 3′-end. These findings indicate that a single conformation of the complex occurs in solution and suggest that the binding sites are at the level of the external tetrads.
The proton assignment and the inter-residue NOE connectivities characterizing the three tetrads in the complex are described in the Experimental section and the values are reported in Tables S3 and S4.
Many NOE interactions between curaxin and the nucleotide were revealed in the NOESY spectra (Table 2, Figure 6 and Figure S1). A large portion of the curaxin molecule, going from H1 and H2 to H7 and H8, including the side chain at N9, has contacts with the aromatic protons of the guanines G7, G11 and G20 at 5′-end. At the same terminal, the CH3CO groups of curaxin show strong interactions with the imino H1 protons of the guanines of the above tetrad.
The 3′-end terminal appears more compact. H4,5 and CH3CO groups present strong NOE interactions with the imino protons H1 of G18 and H1 of G22 and/or G13, respectively. In addition, both aromatic and methyl protons show significant contacts with the tail of the flanking chain A24 and A25.
These results indicate that curaxin binds the Pu22T14T23 quadruplex over the two external tetrads. The location at 3′-end appears more stable, while the interaction at 5′-end terminal seems to be characterized by a higher mobility of the ligand. This relative mobility is also suggested by the finding of additional weak NOE interactions involving the ribose H1′ proton of G7 and T4 units with H1,8 and the side chain of the ligand.
The three-dimensional models for the curaxin-Pu22T14T23 complexes were obtained by performing molecular docking experiments, followed by a Molecular Dynamics (MD) optimization of the resulting complexes (Figure 7). In both 3′-end and 5′-end positions, the curaxin molecule is arranged along the main groove of Pu22.
At 3’-end, the ligand is positioned towards the center of the tetrad and is stabilized by a π–π interaction involving G13 and one of the two benzene rings. The side chain faces A25, forming two hydrogen bonds between the quaternary nitrogen protons and N1A25, with distances of 2.69 Å and 2.57 Å. This explains the downfield shift of A25 H2 (Δδ + 0.21 ppm), which has lost the shielding by the aromatic system of the guanine G22. No particular interactions involving the two CH3CO residues are observed. (Figure 7A,B).
At 5’-end, curaxin is stabilized by an extensive network of π–π interactions involving the underlying 5’-end G-tetrad, with the tricyclic moiety located near the center of the tetrad (Figure 7C,D). The aromatic system interacts with the π systems of G5, G7, G11 and G16. The complex is held in place by two cation-π interaction. The first between the potassium ion and the pyrrole moiety (4.93 Å) and the second between the quaternary nitrogen of the side chain and the aromatic system of G7. The system is further stabilized by a hydrogen bond between one carbonyl group (CH3CO) and G7 H1 (3.06 Å), and by a bidentate hydrogen bond between the quaternary nitrogen of the curaxin and G5 N7 (2.82 Å) and A6 N1 (2.70 Å).
The best conformations of the complexes at 5’-end and 3’-end are in agreement with the reported NOE contacts (Table 2).

2.3. CD and Fluorescence Studies of the Complexes between Curaxin and PuT14T23 and d(TTAGGGT)4

The CD spectra of Pu22T14T23 showed a negative band at 245 nm and a positive band at 262 nm, which are characteristic of a parallel G-quadruplex topology. Upon heating, the G-quadruplex structure unfolded with a Tm equal to 89 ± 1 °C, which indicates a high thermal stability of this structure at the experimental conditions (Figure S2a). The addition of curaxin did not affect the intensity and the shape of the CD spectra of the G-quadruplex (Figure S2b). The ellipticity trace at 265 nm (Figure S2c) indicates clearly that the ligand did not compromise the stability of the structure. The determined Tm value (91 °C) suggests that curaxin produces a certain stabilization, preserving the overall G-quadruplex structure.
Considering the high stability of the tested G-quadruplex structure, CD analysis was performed also without added KCl (5 mM potassium phosphate buffer, pH 7.1), so to better investigate the compound effect. The Tm values in absence and presence of curaxine (1:3 molar ratio) were 70.0 ± 0.3 and 75.1 ± 0.4 °C, respectively (Figure S3). This confirmed that curaxine effectively stabilizes the folded G-quadruplex structure.
The fluorescence-monitored titration involved both the titration of curaxin with Pu22T14T23 and the titration of Pu22T14T23 with curaxin (Figure 8a and Figure S4a, respectively). In Figure 8a, a decreasing of the fluorescence signal intensity is observed, whereas the opposite occurs in Figure S4a, where an increasing of the fluorescence signal intensity is detected. In both cases, the estimation of the stoichiometry and the binding constants (Kb) relative to the interaction with curaxin were performed with the EQUISPEC program, which is based on the multivariate analysis of the whole spectra measured along the titration. When DNA is the titrating agent, the 1:1 complex is favored over the formation of complexes with higher stoichiometries (Figure 8b). On the contrary, complexes with a higher number of curaxin molecules occur when the DNA is titrated with the ligand [25]. From the titration curves, Kb equal to 5.2 ± 1.3 × 106 M−1 and 7.8 ± 2.4 × 1012 M−2 were calculated for the 1:1 and 1:2 (DNA:curaxin) complexes, respectively.
Similar results were obtained when the interaction of curaxin with d(TTAGGGT)4 was studied.
The titrations involving curaxin and d(TTAGGGT)4 showed similar trends to the interaction with Pu22T14T23. Hence, the addition of DNA to curaxin produced a decrease of fluorescence, whereas the opposite effect was observed along the titration of curaxin with DNA. The data analysis was carried out in a similar way to Pu22T14T23, showing the formation of 1:1 and 1:2 (DNA:curaxin) complexes with Kb equal to 1.0 ± 1.3 × 106 M−1 and 10.2 ± 1.0 × 1010 M−2, respectively. These values are slightly lower than in the case of those calculated for the interaction with Pu22T14T23.
The titration studies were also performed on non-G-quadruplex DNA oligonucleotides, both single and double stranded. The results of fluorescence-monitored titration of curaxine with the self-complementary sequence 5′-CGTACG-3′ showed the formation of a 1:1 complex with a binding constant of 2.0·106 × M−1 (Figure S5). When an unfolded strand (5′-CTCTCTACTACCCTTCTGCTC-3′) was titrated, the fluorescence decrease of the ligand upon addition of DNA was clearly lower (Figure S6). This fact could reflect a weaker interaction of the ligand with the unfolded strand than with the G-quadruplex or duplex structure.

2.4. 1H and 31P NMR of Curaxin with Double Helix B-DNA d(CGTACG)2 and d(AAGAATTCTT)2

To further investigate the intercalation mode of curaxin into the double-helix DNA, the study was extended to models of CG and AT-rich sequences, the self-complementary double helix oligomers d(CGTACG)2 and d(AAGAATTCTT)2, respectively.
In the 1H NMR spectra, both B-DNA oligonucleotides displayed, in the 1H NMR spectra, signals from 12 to 14 ppm, assigned to the NH imino protons of the GC and AT base pairs, confirming that both oligomers adopt a double-helix conformation.
The titration of d(CGTACG)2 with curaxin showed, at first, the appearance in the region of 12–14 ppm, of three sets of NH imino signals (Figure S7). Increasing the R value until 3.0, the three sets of signals collapsed in one pattern of three signals assigned to a bound species. This is explained by the formation of two different complexes in slow chemical exchange, on the NMR time scale, together with the presence of the free species (Figure S7). The final situation at R = 3 shows a single complex, with up field shift of the aromatic protons of G2 and C5 units (Table S8). This suggests a binding of curaxin at the level of these residues validated by the NOE interactions between aromatic protons of curaxin with aromatic protons G2, T3 and A4 (Table S6).
The 31P NMR spectra of d(CGTACG)2 with curaxin showed a significant low field shift of all the signals. It is known that 31P resonance is a sensitive probe to detect changes in the phosphorus chain of the oligonucleotide [26,27]. Thus, a low field shift is an indication of intercalation processes. In our case, the significant shift found for all the signals showed that the intercalation of curaxin occurs at different sites, with an exchange among them, suggested by the broadening of the signals (Figure S8).
As concerns the interaction of curaxin with AT-rich sequence, the titration experiment led to a generalized line broadening of 1H and 31P signals (Figure S8). At R ≥ 1.0, a precipitate did not allow us to complete the experiment. The broadening of the signals should suggest some kind of interaction, but no other information can be deduced. The 31P spectra present insignificant shift variation, excluding an intercalation process.

3. Materials and Methods

3.1. Ligand

Curaxin CBL0137 was purchased from Biosynth Carbosynth. The corresponding hydrochloride 1 was prepared by treatment with HCl 4M in dioxane and was used for the titration experiments.
To a solution of curaxin CBL0137 (5 mg, 0.015 mmoli) in dioxane (0.4 mL), 4M HCl in dioxane (5 µL) was added, and the mixture was stirred 30 min, at RT. The solvent was removed to obtain 5.5 mg of desired compound.

3.2. Nuclear Magnetic Resonance

The NMR sample of G-quadruplex Pu22T14T23 d(TGAGGGTGGGTAGGGTGGGTAA) was prepared at 0.34 mM concentration in 25 mM KH2PO4, 70 mM KCl, pH 6.9, 10% D2O. d(TTAGGGT)4 was prepared at 0.45 mM concentration in G-quadruplex, in 25 mM KH2PO4, 150 mM KCl and 1 mM EDTA, pH 6.7, 10% D2O. The NMR samples of double helixes d(CGTACG)2 and d(AAGAATTCTT)2 were prepared at 0.3 mM concentration in 10 mM NaH2PO4, 100 mM NaCl, pH 7.0, 10% D2O. The oligonucleotide samples were heated to 85 °C for 1 min and then cooled at room temperature overnight. Curaxin was dissolved in DMSO-d6 at concentration of 27 mM.
We acquired 1H NMR spectra at 15 and 25 °C with a Bruker AV600, 600MHz spectrometer, equipped with a TXI probe with z-gradient, and processed with TOPSPIN 2.1 software. We recorded 31P NMR spectra at 242.94 MHz with a broad-band probe and referenced at MDA (external reference). Two-dimensional NOESY spectra were acquired with mixing time between 150 and 400 ms and 2D TOCSY spectra with mixing time of 60 ms. Heteronuclear one-bond 1H/13C (HSQC) was carried out in 1H detection mode with broad band decoupling in the 13C domain.
Proton resonance assignments of the free d(TTAGGGT)4 [28], Pu22T14T23 [29], d(CGTACG)2 [30] and d(AAGAATTCTT)2 [31] sequences were performed on the basis of previous assignments by TOCSY and NOESY experiments. The assignments in the complexes were carried out by following the same procedure (Tables S2, S3 and S5). In particular, for G-quadruplexes, the guanine protons through the sequential NOE interactions and the inter-residue NOE connectivities between H1 and H8 resonances are characteristic of the three tetrads, thus confirming the conserved quadruplex structure. The assignment of the thymine H6 protons was followed by the NOE interactions with the methyl signals, which were easily recognized at high fields. The adenine aromatic protons were assigned through their 1H/13C coupling by HSQC experiments. Curaxin protons were assigned, in K-buffer solution pH 6.7, by an integrated series of 2D experiments, such as ROESY, TOCSY and COSY (Table S1). Peak assignment was carried out with Sparky programs.

3.3. Molecular Modeling Studies

The coordinates for the starting models of Pu22T14T23 and d(TTAGGGT)4 were obtained from the NMR structure deposited in the Protein Data Bank (accession code: 2L7V for Pu22T14T23 and 1NZM for d(TTAGGGT)4) [32,33].
The molecular docking calculations were performed by using the AutoDock 4.2 software [34]. The Lamarckian Genetic Algorithm [35], and the AutoDock Toolkit (ADT) [36] were used to further process the ligand and the Pu22T14T23 and d(TTAGGGT)4 models. The Gasteiger–Marsili charges [37] were added to the ligand by using ADT, and the phosphorus atoms in the DNA were parameterized by using the Cornell parameters. The solvation parameters were added to the system by means of the Addsol utility of AutoDock. The initial population consisted of 100 randomly placed individuals for each docking run, with a maximum number of 250 energy evaluations and an elitism value of 1, a mutation rate of 0.02 and a crossover rate of 0.80. The local search for the ligand was conducted by using 250 independent docking runs by applying the so-called pseudo-Solis and Wets algorithm with a maximum of 250 iterations per local search. The grid maps used in the actual docking process were calculated with Autogrid and cantered between the two K+ ions, with a grid dimensions of 80 × 80 × 80 Å (spacing of 0.01 Å). The docking results were scored by using an in-house version of the simpler intermolecular energy function based on the Weiner force field, and the results differing by less than 1.0 Å in positional root-mean-square deviation (rmsd) were clustered together and represented by the most favorable free energy of binding.
The best poses obtained from the previous phase were optimized through 5.0 ns of molecular dynamics (MD) using an OpenCL version of the GROMACS package with a modified version of the 53A6 GROMOS force field [38,39] and running on a dual-Xeon workstation (8 core) equipped with an NVIDIA GPU containing about 5000 CUDA® cores. The systems were placed in the centre of a box with boundaries at 2.0 nm apart from all atoms. Then 3′-end and 5′-end terminal nucleotide topologies were modified according to Ricci et al. [40], the counterions (K+ ions) were random placed and SPC water molecules were added to the systems. The fully solvated systems were optimized through 100 ps of position restrained MD, followed by a heating ramp of short (100 ps) consecutive simulations at 50, 100, 150, 200, 250, and 300 K. The production simulations consisted of 5 ns of partially restrained MD at 310 K (time step of 0.002 ps).
The Lennard–Jones interactions were calculated by using a two-range switch interaction (cut-off radius of 0.9 and 1.1 nm), while the constraints were obtained by using the Lincs [41] and SETTLE [42] algorithms.
A Berendsen thermostat (coupling time of 0.1 ps) was applied to the systems [43], and the electrostatic interactions were calculated by using the Particle Mesh Ewald (PME) [44,45] method (Coulomb cut-off radius of 1.2 nm).

3.4. CD and Fluorescence

CD spectra were recorded on a Jasco J-810 spectropolarimeter equipped with a Peltier temperature control unit (Seelbach, Germany). The DNA solution (Pu22T14T23 or d(TTAGGGT)4 was transferred to a covered cell and ellipticity was recorded with a heating rate of approximately 0.4 °C·min−1. Simultaneously, CD spectra were recorded every 5 °C from 220 to 310 nm. The spectrum of the buffer was subtracted. Each sample was allowed to equilibrate at the initial temperature for 30 min before the melting experiment began. In all experiments, the concentration of DNA was kept constant (3 µM) whereas the concentration of the considered ligands was increased. The medium consisted of 25 mM KH2PO4 and 70 mM (Pu22T14T23) or 150 mM (d(TTAGGGT)4) KCl [46].
Molecular fluorescence spectra were measured with a JASCO FP-6200 spectrofluorimeter. The temperature was controlled at 20 °C using a water bath. The fluorescence spectra were monitored by using a quartz cuvette with a 10-mm path length, with the excitation and emission slits set at 10 nm, and the scan speed at 250 nm/min. Measurements were taken at 334 nm excitation wavelength. The buffer consisted of 25 mM phosphate buffer (pH 6.9) and 70 mM KCl. In all experiments, the concentration of curaxin was 3 µM, whereas the concentration of the considered DNA sequence was increased. The determination of the ratio ligand:DNA and the calculation of the binding constants were conducted from the fluorescence data recorded along titrations of ligands with DNAs by using the EQUISPEC program [47]. This program is based on the multivariate analysis of the whole spectra measured along the titration.

4. Conclusions

Curaxins are a class of substituted carbazoles that exert anticancer activity by complex and diverse mechanisms. Among curaxins, CBL0137 is under phase I clinical trials in patients with hematological malignancies and solid tumors.
Recent findings have claimed the involvement of DNA binding in curaxin antitumor activity [14,15].
In particular, numerous in vivo and in vitro experimental evidences, including inhibition of c-MYC expression, DNA methyltransferase inhibition, and chromatin remodelling, could be consistent with curaxin DNA G-quadruplex binding [4,14,15,17,18,19].
To gain a deeper insight into the potential DNA-binding properties of curaxin CBL0137 and to support the cellular evidences, we investigated its ability to interact with G-quadruplex DNA.
An NMR study was performed with the single repeat sequence of human telomers (TTAGGGT)4 and with Pu22T14T23, a model of the c-myc promoter Pu22 sequence. The results evidenced that the binding of curaxin to these G-quadruplex structures is significant. In the first case, curaxin was located over the G4 and G6 tetrads, with some exchange process between free and bound state at the level of G4. With the Pu22 sequence, the binding occurs at 3′-end and 5′-end, over the two external tetrads of G-quadruplex. The location at 3′-end appears more stable, while the interaction at 5′-end terminal is characterized by a higher mobility of the ligand. This behavior is in line with that observed for the telomere quadruplex (TTAGGGT)4. The results obtained from the detection of numerous NOE contacts between curaxin and both the G-quadruplex structures, were validated by molecular modeling studies.
The Kb values obtained for the 1:1 interaction of ligand 1 with Pu22T14T23 and d(TTAGGGT)4 are in the order of 106 M−1, which indicating a significant interaction between this ligand and both structures. At high concentrations of curaxin, the formation of the 1:2 DNA:ligand complex is favored, especially in the case of Pu22T14T23. Curaxin is able to intercalate between the base pairs of B-DNA with a preference for CG-rich sequence. The relative small dimension of the flat aromatic portion of the molecule favors the entrance between the base pairs promoting the exchange among different sites. Overall, the results are in agreement with the literature’s available cellular data [4,14,15,17] and suggest that the interaction with G-quadruplex DNA may play a role in the anticancer activity of curaxin CBL0137. The molecular models here built for the complexes with DNA G-quadruplex structures could be a precious source of inspiration for the design of curaxin-related more active and selective DNA-binding ligands.

Supplementary Materials

Author Contributions

S.M. and S.D. conceived the study and designed the experiments. L.M. prepared compound 1. R.A. performed molecular modeling studies. R.G., A.A. and C.F. performed CD and fluorescence experiments under the supervision of R.E., L.M.M. and S.M. performed the NMR experiments. L.M. and S.M. organized data. S.M. and S.D. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by PIANO DI SOSTEGNO ALLA RICERCA 2020—Linea 2 azione B (DEFENS) and by the Spanish Ministerio de Ciencia e Innnvacion (PID2019-107158GB-I00) and recognition from the Autonomous Catalan government (2017SGR114).

Acknowledgments

We would like to thank Rosanna Mondelli for her help during manuscript preparation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Di Bussolo, V.; Minutolo, F. Curaxins: A New Family of Non-genotoxic Multitargeted Anticancer Agents. ChemMedChem 2011, 6, 2133–2136. [Google Scholar] [CrossRef] [PubMed]
  2. Gurova, K.V.; Hill, J.E.; Guo, C.; Prokvolit, A.; Burdelya, L.G.; Samoylova, E.; Khodyakova, A.V.; Ganapathi, R.; Ganapathi, M.; Tararova, N.D.; et al. Small molecules that reactivate p53 in renal cell carcinoma reveal a NF-κB-dependent mechanism of p53 suppression in tumors. Proc. Natl. Acad. Sci. USA 2005, 102, 17448–17453. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Gasparian, A.V.; Burkhart, C.A.; Purmal, A.A.; Brodsky, L.; Pal, M.; Saranadasa, M.; Bosykh, D.A.; Commane, M.; Guryanova, O.A.; Pal, S.; et al. Curaxins: Anticancer Compounds That Simultaneously Suppress NF-κB and Activate p53 by Targeting FACT. Sci. Transl. Med. 2011, 3, 95ra74. [Google Scholar] [CrossRef] [PubMed]
  4. Wang, J.; Sui, Y.; Li, Q.; Zhao, Y.; Dong, X.; Yang, J.; Liang, Z.; Han, Y.; Tang, Y. Effective inhibition of MYC-amplified group 3 medulloblastoma by FACT-targeted curaxin drug CBL0137. Cell Death Dis. 2020, 11, 1029. [Google Scholar] [CrossRef]
  5. De, S.; Lindner, D.J.; Coleman, C.J.; Wildey, G.; Dowlati, A.; Stark, G.R. The FACT inhibitor CBL0137 Synergizes with Cisplatin in Small-Cell Lung Cancer by Increasing NOTCH1 Expression and Targeting Tumor-Initiating Cells. Cancer Res. 2018, 78, 2396–2406. [Google Scholar] [CrossRef] [Green Version]
  6. Lindner, D.J.; Wildey, G.; Parker, Y.; Dowlati, A.; Stark, G.R.; De, S. CBL0137 increases the targeting efficacy of Rovalpituzumab tesirine against tumour-initiating cells in small cell lung cancer. Br. J. Cancer 2021, 124, 893–895. [Google Scholar] [CrossRef] [PubMed]
  7. Jin, M.Z.; Xia, B.R.; Xu, Y.; Jin, W.L. Curaxin CBL0137 Exerts Anticancer Activity via Diverse Mechanisms. Front. Oncol. 2018, 8, 598. [Google Scholar] [CrossRef]
  8. Somers, K.; Kosciolek, A.; Bongers, A.; El-Ayoubi, A.; Karsa, M.; Mayoh, C.; Henderson, M.J. Potent antileukemic activity of curaxin CBL0137 against MLL-rearranged leukemia. Int. J. Cancer 2020, 146, 1902–1916. [Google Scholar] [CrossRef]
  9. Burkhart, C.; Fleyshman, D.; Kohrn, R.; Commane, M.; Garrigan, J.; Kurbatov, V.; Gurova, K.V. Curaxin CBL0137 eradicates drug resistant cancer stem cells and potentiates efficacy of gemcitabine in preclinical models of pancreatic cancer. Oncotarget 2014, 5, 11038–11053. [Google Scholar] [CrossRef] [Green Version]
  10. Available online: https://www.clinicaltrials.gov/ct2/results?cond=&term=CBL0137 (accessed on 29 April 2021).
  11. Safina, A.; Cheney, P.; Pal, M.; Brodsky, L.; Ivanov, A.; Kirsanov, K.; Lesovaya, E.; Naberezhnov, D.; Nesher, E.; Koman, I.; et al. FACT is a sensor of DNA torsional stress in eukaryotic cells. Nucleic Acids Res. 2017, 45, 1925–1945. [Google Scholar] [CrossRef] [Green Version]
  12. Nesher, E.; Safina, A.; Aljahdali, I.; Portwood, S.; Wang, E.S.; Koman, I.; Wang, J.; Gurova, K.V. Role of Chromatin Damage and Chromatin Trapping of FACT in Mediating the Anticancer Cytotoxicity of DNA-Binding Small-Molecule Drugs. Cancer Res. 2018, 78, 1431–1443. [Google Scholar] [CrossRef] [Green Version]
  13. Leonova, K.; Safina, A.; Nesher, E.; Sandlesh, P.; Pratt, R.; Burkhart, C.; Lipchick, B.; Gitlin, I.; Frangou, C.; Koman, I.; et al. TRAIN (Transcription of Repeats Activates INterferon) in response to chromatin destabilization induced by small molecules in mammalian cells. eLife 2018, 7, e30842. [Google Scholar] [CrossRef]
  14. Kantidze, O.L.; Luzhin, A.V.; Nizovtseva, E.V.; Safina, A.; Valieva, M.E.; Golov, A.K.; Velichko, A.K.; Lyubitelev, A.V.; Feofanov, A.V.; Gurova, K.V.; et al. The anti-cancer drugs curaxins target spatial genome organization. Nature Commun. 2019, 10, 1441–1452. [Google Scholar] [CrossRef] [Green Version]
  15. Lu, K.; Liu, C.; Liu, Y.; Luo, A.; Chen, J.; Lei, Z.; Kong, J.; Xiao, X.; Zhang, S.; Wang, Y.-Z.; et al. Curaxin-Induced DNA Topology Alterations Trigger the Distinct Binding Response of CTCF and FACT at the Single-Molecule Level. Biochemistry 2021, 60, 494–499. [Google Scholar] [CrossRef] [PubMed]
  16. Kim, M.; Neznanov, N.; Wilfong, C.D.; Fleyshman, D.I.; Purmal, A.A.; Haderski, G.; Stanhope-Baker, P.; Burkhart, C.A.; Gurova, K.V.; Gudkov, A.V.; et al. Preclinical Validation of a Single-Treatment Infusion Modality That Can Eradicate Extremity Melanomas. Cancer Res. 2016, 76, 6620–6630. [Google Scholar] [CrossRef] [Green Version]
  17. Sergeev, A.; Vorobyova, A.; Yakubovskaya, M.; Kirsanova, O.; Gromova, E. Novel anticancer drug curaxin CBL0137 impairs DNA methylation by eukaryotic DNA methyltransferase Dnmt3a. Bioorg. Med. Chem. Lett. 2020, 30, 127296–127300. [Google Scholar] [CrossRef] [PubMed]
  18. Wang, W.; Hu, S.; Gu, Y.; Yan, Y.; Stovall, D.B.; Li, D.; Sui, G. Human MYC G-quadruplex: From discovery to a cancer therapeutic target. BBA Rev. Cancer 2020, 1874, 188410. [Google Scholar] [CrossRef] [PubMed]
  19. Jara-Espejo, M.; Line, S.R. DNA G-quadruplex stability, position and chromatin accessibility are associated with CpG island methylation. FEBS J. 2020, 287, 483–495. [Google Scholar] [CrossRef]
  20. Varizhuk, A.; Isaakova, E.I.; Pozmogova, G. DNA G-Quadruplexes (G4s) Modulate Epigenetic (Re) Programming and Chromatin Remodeling. BioEssays 2019, 41, 1900091. [Google Scholar] [CrossRef] [Green Version]
  21. Siddiqui-Jain, A.; Grand, C.L.; Bearss, D.J.; Hurley, L.H. Direct evidence for a G-quadruplex in a promoter region and its targeting with a small molecule to repress c-MYC transcription. Proc. Natl. Acad. Sci. USA 2002, 99, 11593–11598. [Google Scholar] [CrossRef] [Green Version]
  22. Agrawal, P.; Lin, C.; Mathad, R.I.; Carver, M.; Yang, D. The Major G-Quadruplex Formed in the Human BCL-2 Proximal Promoter Adopts a Parallel Structure with a 13-nt Loop in K+ Solution. J. Am. Chem. Soc. 2014, 136, 1750–1753. [Google Scholar] [CrossRef]
  23. Musso, L.; Mazzini, S.; Rossini, A.; Castagnoli, L.; Scaglioni, L.; Artali, R.; Di Nicola, M.; Zunino, F.; Dallavalle, S. c-MYC G-quadruplex binding by the RNA polymerase I inhibitor BMH-21 and analogues revealed by a combined NMR and biochemical approach. BBA Gen. Subj. 2018, 1862, 615–629. [Google Scholar] [CrossRef]
  24. Carvalho, J.; Mergny, J.L.; Salgado, G.F.; Queiroz, J.A.; Cruz, C. G-quadruplex, Friend or Foe: The Role of the G-quartet in Anticancer Strategies. Trends. Mol. Med. 2020, 26, 848–861. [Google Scholar] [CrossRef] [PubMed]
  25. Gabelica, V.; Maeda, R.; Fujimoto, T.; Yaku, H.; Murashima, T.; Sugimoto, N.; Miyoshi, D. Multiple and Cooperative Binding of Fluorescence Light-up Probe Thioflavin T with Human Telomere DNA G-Quadruplex. Biochemistry 2013, 52, 5620–5628. [Google Scholar] [CrossRef]
  26. Misiak, M.; Heldt, M.; Szeligowska, M.; Mazzini, S.; Scaglioni, L.; Grabe, G.J.; Serocki, M.; Lica, J.; Switalska, M.; Wietrzyk, J.; et al. Molecular basis for the DNA damage induction and anticancer activity of asymmetrically substituted anthrapyridazone PDZ-7. Oncotarget 2017, 8, 105137–105154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Gorenstein, D.G. Phosphorus-31 NMR. In Principles and Applications; Academic Press: New York, NY, USA, 1984. [Google Scholar]
  28. Gavathiotis, E.; Searle, M.S. Structure of the parallel-stranded DNA quadruplex d (TTAGGGT) 4 containing the human telomeric repeat evidence for A-tetrad formation from NMR and molecular dynamics simulations. Org. Biomol. Chem. 2003, 1, 1650–1656. [Google Scholar] [CrossRef]
  29. Ambrus, A.; Chen, D.; Dai, J.; Jones, R.A.; Yang, D. Solution Structure of the Biologically Relevant G Quadruplex Element in the Human c-MYC Promoter. Implications for G-Quadruplex Stabilization. Biochemistry 2004, 44, 2048–2058. [Google Scholar] [CrossRef] [PubMed]
  30. Mazzini, M.; Mondelli, R.; Ragg, E. Structure and dynamics of intercalation complexes of anthracyclines with d (CGATCG) 2 and d(CGTACG)2. 2D-1H and 31P NMR investigations. J. Chem. Soc. Perkin Trans. 1998, 2, 1983–1991. [Google Scholar] [CrossRef]
  31. Mazzini, M.; Bellucci, M.C.; Mondelli, R. Mode of Binding of the Cytotoxic Alkaloid Berberine with theDouble Helix Oligonucleotide d (AAGAATTCTT) 2. Bioorg. Med. Chem. 2003, 11, 505–514. [Google Scholar] [CrossRef]
  32. Dai, J.; Carver, M.; Hurley, L.H.; Yang, D. Solution structure of a 2:1 quindoline/c-MYC G-quadruplex: Insights into G-quadruplex-interactive small molecule drug design. J. Am. Chem. Soc. 2011, 133, 17673–17680. [Google Scholar] [CrossRef] [Green Version]
  33. Gavathiotis, E.; Heald, R.A.; Stevens, M.F.; Searle, M.S. Drug recognition and stabilisation of the parallel-stranded DNA quadruplex d (TTAGGGT) 4 containing the human telomeric repeat. J. Mol. Biol. 2003, 334, 25–36. [Google Scholar] [CrossRef]
  34. Huey, R.; Morris, G.M.; Olson, A.J.; Goodsell, D.S. A semiempirical free energy force field with charge-based desolvation. J. Comput. Chem. 2007, 28, 1145–1152. [Google Scholar] [CrossRef]
  35. Solis, F.J.; Wets, R.J.B. Minimization by Random Search Techniques. Math. Oper. Res. 1981, 6, 19–30. [Google Scholar] [CrossRef]
  36. Sanner, M.F. Python: A programming language for software integration and development. J. Mol. Graph. Model. 1999, 17, 57–61. [Google Scholar] [PubMed]
  37. Gasteiger, J.; Marsili, M.M. Iterative partial equalization of orbital electronegativity-a rapid access to atomic charges. Tetrahedron 2008, 36, 3219–3228. [Google Scholar] [CrossRef]
  38. Lindahl, E.; Hess, B.; Van Der Spoel, D. GROMACS 3.0: A package for molecular simulation and trajectory analysis. J. Mol. Model. 2001, 7, 306–317. [Google Scholar] [CrossRef]
  39. Oostenbrink, C.; Soares, T.A.; Van Der Vegt, N.F.A.; Van Gunsteren, W.F. Validation of the 53A6 GROMOS force field. Eur. Biophys. J. 2005, 34, 273–284. [Google Scholar] [CrossRef] [Green Version]
  40. Ricci, G.C.; De Andrade, A.S.C.; Mottin, M.; Netz, P.A. Molecular dynamics of DNA: Comparison of force fields and terminal nucleotide definitions. J. Phys. Chem. B 2010, 114, 9882–9893. [Google Scholar] [CrossRef]
  41. Hess, B.; Bekker, H.; Berendsen, J.J.C.; Fraaije, J.G.E.M. LINCS: A linear constraint solver for molecular simulations. J. Comput. Chem. 1997, 18, 1463–1472. [Google Scholar] [CrossRef]
  42. Miyamoto, S.; Kollman, P.A. Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models. J. Comput. Chem. 1992, 13, 952–962. [Google Scholar] [CrossRef]
  43. Berendsen, H.J.C.; Postma, J.P.M.; van Gunsteren, W.F.; Di Nola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef] [Green Version]
  44. Darden, T.; York, D.; Pedersen, L.G. Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089–10092. [Google Scholar] [CrossRef] [Green Version]
  45. Essmann, U.; Perera, L.; Berkowitz, M.L.; Darden, T.; Lee, H.; Pedersen, L.G. A smooth particle mesh Ewald method. J. Chem. Phys. 1995, 103, 8577–8592. [Google Scholar] [CrossRef] [Green Version]
  46. Mazzini, S.; Gargallo, R.; Musso, L.; De Santis, F.; Aviñó, A.; Scaglioni, L.; Eritja, R.; Di Nicola, M.; Zunino, F.; Amatulli, A.; et al. Stabilization of c-KIT G-Quadruplex DNA Structures by the RNA Polymerase I Inhibitors BMH-21 and BA-41. Int. J. Mol. Sci. 2019, 20, 4927. [Google Scholar] [CrossRef] [Green Version]
  47. Dyson, R.M.; Kaderli, S.; Lawrance, G.A.; Maeder, M. Second order global analysis: The evaluation of series of spectrophotometric titrations for improved determination of equilibrium constants. Anal. Chim. Acta 1997, 353, 381–393. [Google Scholar] [CrossRef]
Figure 1. Structure of curaxin CBL0137.
Figure 1. Structure of curaxin CBL0137.
Ijms 22 06476 g001
Figure 2. (a) Imino proton region of the 1D NMR titration spectra of d(TTAGGGT)4 with curaxin at 25 °C in H2O/D2O (9:1), 25 mM KH2PO4, 150 mM KCl and 1 mM EDTA, pH 6.7, at different R = [drug]/[DNA] ratios; (b) schematic representation of d(TTAGGGT)4/curaxin (in red) complex.
Figure 2. (a) Imino proton region of the 1D NMR titration spectra of d(TTAGGGT)4 with curaxin at 25 °C in H2O/D2O (9:1), 25 mM KH2PO4, 150 mM KCl and 1 mM EDTA, pH 6.7, at different R = [drug]/[DNA] ratios; (b) schematic representation of d(TTAGGGT)4/curaxin (in red) complex.
Ijms 22 06476 g002
Figure 3. Selected region of the 2D NOESY spectrum of d(TTAGGGT)4/curaxin complex. (a) G6 imino proton displays intermolecular NOEs between CH3T7 and CH3(isopropyl) of curaxin; (b) A3, G6 and T7 aromatic protons and methyl groups of T2 and T7 of T2AG3T display intermolecular NOEs with CH3 (isopropyl) and aromatic protons of curaxin; (c) schematic representation of d(TTAGGGT)4/curaxin (in red) complex.
Figure 3. Selected region of the 2D NOESY spectrum of d(TTAGGGT)4/curaxin complex. (a) G6 imino proton displays intermolecular NOEs between CH3T7 and CH3(isopropyl) of curaxin; (b) A3, G6 and T7 aromatic protons and methyl groups of T2 and T7 of T2AG3T display intermolecular NOEs with CH3 (isopropyl) and aromatic protons of curaxin; (c) schematic representation of d(TTAGGGT)4/curaxin (in red) complex.
Ijms 22 06476 g003
Figure 4. Graphical representations of the curaxin-d(TTAGGGT)4 complexes at the ApG and GpT intercalation sites, obtained by Molecular Docking and optimized by Molecular Dynamics (MD). (A) Side view of the ghostly white solvent accessible surface (SAS) of the d(TTAGGGT)4 quadruplex. The nucleotides are represented in stick and filled rings: adenine in red, guanine in green and thymine in blue. The ligand is represented as van der Waals (vdW) spheres and colored following the CPK code. The optimized conformations of the ligand are represented in (B) for the complex at ApG and in (C) for the complex at GpT. Potassium ions are represented by their vdW spheres (K+ in purple), while the ligand is depicted in stick and colored following the CPK code. The nucleotides are represented as filled plates: adenine in red, guanine in green and thymine in blue.
Figure 4. Graphical representations of the curaxin-d(TTAGGGT)4 complexes at the ApG and GpT intercalation sites, obtained by Molecular Docking and optimized by Molecular Dynamics (MD). (A) Side view of the ghostly white solvent accessible surface (SAS) of the d(TTAGGGT)4 quadruplex. The nucleotides are represented in stick and filled rings: adenine in red, guanine in green and thymine in blue. The ligand is represented as van der Waals (vdW) spheres and colored following the CPK code. The optimized conformations of the ligand are represented in (B) for the complex at ApG and in (C) for the complex at GpT. Potassium ions are represented by their vdW spheres (K+ in purple), while the ligand is depicted in stick and colored following the CPK code. The nucleotides are represented as filled plates: adenine in red, guanine in green and thymine in blue.
Ijms 22 06476 g004
Figure 5. Imino proton region of the 1D NMR titration spectra of Pu22T14T23 with curaxin at 25 °C in H2O/D2O (9:1), 25 mM KH2PO4 and 70 mM KCl at pH 6.9, at different R = [drug]/[DNA] ratios.
Figure 5. Imino proton region of the 1D NMR titration spectra of Pu22T14T23 with curaxin at 25 °C in H2O/D2O (9:1), 25 mM KH2PO4 and 70 mM KCl at pH 6.9, at different R = [drug]/[DNA] ratios.
Ijms 22 06476 g005
Figure 6. Selected region of the 2D NOESY spectrum of Pu22T14T23/curaxin complex. (a) Some aromatic protons of Pu22T14T23 display intermolecular NOEs with CH3(isopropyl) of curaxin; (b) some aromatic protons of Pu22T14T23 display intermolecular NOEs with aromatic protons of curaxin and (c) schematic representation of Pu22-T14 T23 oligomer G-quadruplex.
Figure 6. Selected region of the 2D NOESY spectrum of Pu22T14T23/curaxin complex. (a) Some aromatic protons of Pu22T14T23 display intermolecular NOEs with CH3(isopropyl) of curaxin; (b) some aromatic protons of Pu22T14T23 display intermolecular NOEs with aromatic protons of curaxin and (c) schematic representation of Pu22-T14 T23 oligomer G-quadruplex.
Ijms 22 06476 g006
Figure 7. Representation of the Pu22T14T23 complexes with curaxin at 5’-end and 3’-end, obtained by Molecular Docking and optimized by Molecular Dynamics (MD). At the center of the figure Pu22T14T23, represented by its solvent accessible surface (SAS), color-coded by the underlying nucleotide (adenine in red, guanine in green and thymine in blue) and complexed with curaxin (represented in stick) at both 3’-end and 5’-end. On the left, lateral (A) and top (B) representation of the ligand conformation at the 3’-end, while on the right we can see the lateral (C) and bottom (D) representation of curaxin at the 5’-end. Ligand and potassium ions are represented by their van der Waals spheres (ligand colored in CPK, K+ in purple), while the nucleotide units of Pu22T14T23 are represented as filled rings: adenine in red, guanine in green and thymine in blue.
Figure 7. Representation of the Pu22T14T23 complexes with curaxin at 5’-end and 3’-end, obtained by Molecular Docking and optimized by Molecular Dynamics (MD). At the center of the figure Pu22T14T23, represented by its solvent accessible surface (SAS), color-coded by the underlying nucleotide (adenine in red, guanine in green and thymine in blue) and complexed with curaxin (represented in stick) at both 3’-end and 5’-end. On the left, lateral (A) and top (B) representation of the ligand conformation at the 3’-end, while on the right we can see the lateral (C) and bottom (D) representation of curaxin at the 5’-end. Ligand and potassium ions are represented by their van der Waals spheres (ligand colored in CPK, K+ in purple), while the nucleotide units of Pu22T14T23 are represented as filled rings: adenine in red, guanine in green and thymine in blue.
Ijms 22 06476 g007
Figure 8. Molecular fluorescence-monitored titration of curaxin with Pu22T14T23. (a) Experimental spectra measured along the titration. Numbers in inset indicate the DNA:curaxin ratio; (b) experimental (symbols) and fitted (line) fluorescence at 450 nm. A 1:1 (DNA:curaxin) stoichiometry was used to model the data.
Figure 8. Molecular fluorescence-monitored titration of curaxin with Pu22T14T23. (a) Experimental spectra measured along the titration. Numbers in inset indicate the DNA:curaxin ratio; (b) experimental (symbols) and fitted (line) fluorescence at 450 nm. A 1:1 (DNA:curaxin) stoichiometry was used to model the data.
Ijms 22 06476 g008
Table 1. Intermolecular NOE and distances from modelling in the curaxin-d(TTAGGGT)4 a complexes.
Table 1. Intermolecular NOE and distances from modelling in the curaxin-d(TTAGGGT)4 a complexes.
ApG Binding Site
NOEd (Å) b
Curaxind(TTAGGGT)4
1,8-HA3H1′4.60
2,7-HA3H84.26
4,5-HG4H1′5.88
CH3 isoA3H83.20
GpT Binding Site
1,8-HG6H85.45
1,8-HG6H14.36
1,8-HT7Me4.74
1,8-HT7H13.13
4,5-HG6H13.55
CH3 isoG6H82.30
CH3 isoG6H14.91
CH3 isoT7H64.84
CH3 isoG5H1′5.91
a Acquired at 25 °C in H2O-D2O (90:10 v/v), 25 mM K-phosphate buffer, 150 mM KCl and 1 mM EDTA, at pH 6.7. b Distances obtained by molecular modelling of the complex.
Table 2. Intermolecular NOEs and distances from modelling in the curaxin-Pu22T14T23 a complex.
Table 2. Intermolecular NOEs and distances from modelling in the curaxin-Pu22T14T23 a complex.
3′-Binding Site NOEd (Å) b
CuraxinPu22T14T23
4,5-HG18H15.69
4,5-HA24H26.60
CH3COG13H15.38
CH3COG22H14.68
CH3 isoA24H26.21 c
CH3 isoA25H24.25
CH3COG18H12.26
5′-Binding Site NOE
1,8-HG7H1′7.39 c
2,7-H G11H83.97
1,8-HG11H85.34
CH3 isoG20H82.75
CH3 isoT4H1′6.22 c
CH3COG7H14.90
CH3COG11H13.45
a Acquired at 25 °C in H2O-D2O (90:10 v/v), 25 mM KH2PO4, 70 mM KCl, pH 6.9. b Distances obtained by molecular modelling of the complex. c The long distance is explained with the mobility of the ligand and the tails of the nucleotide.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dallavalle, S.; Mattio, L.M.; Artali, R.; Musso, L.; Aviñó, A.; Fàbrega, C.; Eritja, R.; Gargallo, R.; Mazzini, S. Exploring the Interaction of Curaxin CBL0137 with G-Quadruplex DNA Oligomers. Int. J. Mol. Sci. 2021, 22, 6476. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22126476

AMA Style

Dallavalle S, Mattio LM, Artali R, Musso L, Aviñó A, Fàbrega C, Eritja R, Gargallo R, Mazzini S. Exploring the Interaction of Curaxin CBL0137 with G-Quadruplex DNA Oligomers. International Journal of Molecular Sciences. 2021; 22(12):6476. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22126476

Chicago/Turabian Style

Dallavalle, Sabrina, Luce M. Mattio, Roberto Artali, Loana Musso, Anna Aviñó, Carme Fàbrega, Ramon Eritja, Raimundo Gargallo, and Stefania Mazzini. 2021. "Exploring the Interaction of Curaxin CBL0137 with G-Quadruplex DNA Oligomers" International Journal of Molecular Sciences 22, no. 12: 6476. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22126476

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop