Abstract
The binding mode of telomestatin to G-quadruplex DNA has been investigated using electrospray mass spectrometry, by detecting the intact complexes formed in ammonium acetate. The mass measurements show the incorporation of one extra ammonium ion in the telomestatin complexes. Experiments on telomestatin alone also show that the telomestatin alone is able to coordinate cations in a similar way as a crown ether. Finally, density functional theory calculations suggest that in the G-quadruplex-telomestatin complex, potassium or ammonium cations are located between the telomestatin and a G-quartet. This study underlines that monovalent cation coordination capabilities should be integrated in the rational design of G-quadruplex binding ligands.
1. Introduction
The formation of G-quadruplex folds by telomeric DNA is thought to play a role in telomere regulation. It has been shown that G-quadruplex ligands binding specifically to the telomeric G-quadruplex structure effectively alter telomere capping and cause the senescence or apoptosis of cancer cells [1–5]. A variety of ligands have now been described as G-quadruplex binders, but a key issue in ligand design is often the specificity for G-quadruplexes over duplex sequences [4, 6–8]. Identifying binding modes that make a ligand a specific and highly active G-quadruplex binder is crucial for the rational design of novel molecules.
Telomestatin (Figure 1) is one of the most emblematic G-quadruplex ligands. The molecule was first extracted from Streptomyces anulatus 3533-SV4 [9]. It is highly specific for G-quadruplexes, with no significant binding to duplexes [10–12]. Telomestatin was found to effectively inhibit the DNA binding of telomere-associated proteins such as telomerase [13], POT1 and TRF2 [14], and even Topo III in ALT cell lines [15]. It therefore induces telomere shortening and apoptosis [16–19] not only via telomerase inhibition but also via telomere uncapping, and therefore has potential activity against many cancer cell types.
(a)
(b)
(c)
(d)
Telomestatin binds to G-quadruplexes, among which is the human telomeric G-quadruplex, by external stacking [12]. One G-quadruplex unit can therefore accommodate two telomestatin ligands, one on each end. A recent modeling study showed that telomestatin has a tendency to capture a potassium ion, either from the G-quadruplex itself or from the solution [20]. Here we show that mass spectrometry provides experimental evidence for the accommodation of one extra cation when a telomestatin molecule is bound to a G-quadruplex. This will be illustrated for three typical G-quadruplexes: the tetramolecular quadruplex, the 4-repeat telomeric sequence , and the Pu22myc promoter sequence . The typical folds adopted by each of these three G-quadruplexes in ammonium acetate were studied previously [21–23] and are summarized in Figure 1.
2. Experimental Section
2.1. Materials
All oligonucleotides were purchased from Eurogentec (Seraing, Belgium), solubilized in water doubly distilled in house, and the 400 M stock solutions were stored at −20°C. The solvents used include methanol (absolute, HPLC grade, Biosolve, Valkenswaard, The Netherlands), bi-distilled water, and aqueous ammonium acetate (5 M stock solution from Fluka, diluted with bi-distilled water). KCl for the evaluation of the complexation of telomestatin alone was puriss, p.a., 99.5% (T) (Fluka). G-quadruplexes were formed by annealing (heating the oligonucleotides for 5 minutes at 85°C, followed by slow cooling to room temperature) in 150 mM ammonium acetate. The G-quadruplex-forming oligonucleotides were (annealed at 200 M single strand to form 50 M tetramolecular G-quadruplex); the telomeric sequence and the Pu22myc promoter sequence (annealed at 50 M single strand to form intramolecular G-quadruplexes). Telomestatin was isolated and purified as described elsewhere [9, 24] to obtain a 1 mM stock solution in DMSO. For the electrospray mass spectrometry analysis of the complexes, the final injected solutions were 5 M in G-quadruplex and 5 to 10 M in telomestatin (only 10 M results are shown), in 80/20 (v/v) aqueous ammonium acetate (150 mM)/methanol.
2.2. Mass Spectrometry
Electrospray mass spectrometry experiments were performed on a Q-TOF Ultima Global (Waters, Manchester, UK). The spectra of the intact G-quadruplexes and their noncovalent complexes with telomestatin were recorded in the negative ion mode (capillary voltage = −2.2 V, source and desolvation temperatures = 70°C, cone = 100 V, RF Lens1 Energy = 45 V, source pirani pressure = 3.94 mbar, collision energy = 10 V), smoothed (mean function, 3 × 20 channels) and subtracted (polynomial order 99, 0.1% below curve). The spectra of telomestatin in the absence of G-quadruplex were recorded in the positive ion mode (capillary voltage = +2.8 V, source and desolvation temperatures = 80 and 100°C, resp., cone = 100 V, RF Lens1 Energy = 50 V, source pirani pressure = 3.33 mbar and collision energy = 7 V) and were not subjected to smoothing or background subtraction.
2.3. Calculations
For the [telomestatin + cation] binary complexes, the ammonium, potassium, and sodium cations were manually docked within the telomestatin ring, and the resulting structures were optimized using density functional theory (DFT) with the hybrid functional B3LYP and the 6-31G(d,p) basis set. For the ternary complexes between [telomestatin + cation + one G-quartet], the telomestatin was manually docked on top of an optimized structure of a G-quartet coordinated with the cation (ammonium, potassium, sodium). The ternary complexes were then optimized using DFT B3LYP at the 6-31G(d,p) level. All electronic structure calculations were performed using the Gaussian 03 rev. D.02 software suite [25]. Comparison with a larger basis set (6-311 + G(d,p)) was performed for one of the complexes (G-tetrad + K + telomestatin), and the results were similar in terms of both energies (0.8 kcal/mol) and geometries (RMSD 0.16 Ǻ). 6-31G(d,p) basis set was therefore used for all calculations. Other hybrid functional BHandHLYP and new meta-GGA hybrid MPWB1K [26] have also been tested for comparison with B3LYP.
3. Results and Discussion
When operated in soft source conditions, electrospray mass spectrometry allows detecting intact noncovalent complexes [27–30]. In the analysis of quadruplex-ligand complexes, it therefore allows determining the number of strands, the number of ligands, and the number of cations in each complex. Electrospray mass spectrometry of nucleic acid noncovalent complexes is typically performed in ammonium acetate solution in order to obtain clean spectra [31]. Ammonium ions present in the counter-ion shell around phosphates are lost during the final stages of desolvation in the electrospray source, even in soft conditions (low acceleration voltages). In contrast, ammonium ions bound sufficiently tightly to the G-quadruplex, such as the ammonium ions trapped between the G-quartets, will persist at higher acceleration voltages in the source than the nonspecifically bound ones [21].
Figure 2 shows the electrospray mass spectra of three typical G-quadruplexes: (a) , (b) the telomeric sequence , and (c) the Pu22myc promoter sequence and on their 1 : 1 complexes with telomestatin (d–f, resp.). The injected mixtures are 5 M in each G-quadruplex 10 M in telomestatin, in 80/20 aqueous ammonium acetate (150 mM)/methanol, and the spectra were recorded using soft source conditions to preserve the specifically bound ammonium ions. The free G-quadruplex (Figure 2()) contains three ammonium ions: the 5− charge state is found at = 1500.22, corresponding to . The major peaks of the telomeric (Figure 2()) and Pu22myc (Figure 2()) G-quadruplexes at charge state 5− are corresponding to the intramolecular G-quadruplex with two ammonium ions, at = 1521.22 and 1419.47, respectively. For the charge state = 5, with the average mass of telomestatin being = 582.5 Da, the between the G-quadruplex and its complex with one telomestatin, with no change in the amount of ammonium ions incorporated, is expected to be equal to 582.5/5 = 116.5. In contrast, the observed between the major peak of the quadruplex and the major peak of its complex with one ligand is equal to 119.9 (compare Figure 2() with Figure 2(), and Figure 2() with Figure 2(), Figure 2() with Figure 2()). This corresponds to the addition of one telomestatin molecule, one extra ammonium ion, and the subtraction of one proton for the charge balance (119.9 = (585.5 + 18−1)/5). The complex with one telomestatin ligand therefore systematically contains one more ammonium ion than the corresponding unbound G-quadruplex. This extra ammonium ion is lost when the acceleration voltage is increased.
In our previous report on the MS detection of telomestatin binding to telomeric DNA [11], we have missed this “extra ammonium” for two reasons. Firstly, for the 3.5-repeat telomeric sequence studied previously it is difficult to preserve two inner ammoniums even in soft conditions. Secondly, soft conditions could not be used because a long duplex had to be detected simultaneously with the G-strand, and the electrospray source conditions were chosen as a compromise.
Electrospray mass spectrometry is also powerful to analyze caged supramolecular complexes such as crown ethers bound to cations [32–35]. To probe whether telomestatin is able to coordinate a cation already in the absence of G-quadruplex, we used ESI-MS in the positive ion mode on telomestatin solutions. Figure 3(a) shows the spectrum obtained with telomestatin dilution in bi-distilled water. The signal-to-noise ratio of all peaks is weak, indicating low protonation and cationization efficiencies. Surprisingly, we found that the major peak was a doubly charged ion at = 686.1 (base peak), whose isotopic distribution matches with that of a complex between two telomestatin ligands and one lead ion adduct. The fragment ion spectrum shows the loss of one telomestatin, and the isotopic distribution of the resulting [Telomestatin + Pb] complex confirms the presence of lead. Traces of lead come from the purification of telomestatin from Streptomyces anulatus 3533-SV4 [24]. Another weak doubly charged peak is tentatively assigned to a complex between two telomestatin and two ions, and the fragment ion spectrum shows only the singly protonated telomestatin.
(a)
(b)
(c)
We then doped the solution with either ammonium acetate or with potassium chloride, in substoichiometric amounts, to probe if telomestatin was able to coordinate monovalent cations such as those typically coordinated to G-quartets. In the presence of ammonium ions, the lead complex totally disappears and the complex [telomestatin + ] is detected at = 600. In the presence of potassium ions, the complex [telomestatin + ] is detected at = 621, but the lead complex has not disappeared completely.
All experimental results suggest that telomestatin has significant affinity for monovalent cations like the ammonium ion, and this could influence its binding mode to the G-quadruplex DNA. We performed DFT calculations in order to ascertain the possible coordination geometries of the monovalent cations to telomestatin. The geometries of the [telomestatin + cation] binary complexes are shown in Figures 4(a)–4(c), and the geometries of the ternary complexes between [telomestatin + cation + one G-quartet] are shown in Figures 4(d)–4(f). In the isolated [telomestatin + cation] complexes, all cations are coordinated in the plane of the telomestatin. In the ternary complexes with the G-quartet, the cation clearly moves towards the G-quartet. The structures of the complexes with and are similar, with the cation coordinated midway between the telomestatin and the G-quartet. Sodium, on the other hand, sits closer to the G-quartet than to the telomestatin.
(a)
(b)
(c)
(d)
(e)
(f)
The mode of cation coordination to the telomestatin—G-quartet system is similar to the coordination mode already described for successive G-quartets [36]. In terms of coordination geometries, potassium tends to sit between G-quartets while sodium tends to fit in the middle of one G-quartet because it is smaller. and have similar ionic radii [37, 38] and therefore behave similarly, and the same trend is observed for our telomestatin-G-quartet complex. The optimized geometries are similar for all functionals tested (B3LYP, BHandHLYP, MPBW1K) (see Figure in Supplementary Material available online at doi: 10.4061/2010/121259), with the cation in the plane of telomestatin in the absence of G-quartet, and between telomestatin and the G-quartet in the ternary complex. However, the functionals have a greater influence on the computed interaction energies. The root mean square distances (RMSDs) for two-by-two comparisons of hybrid functional, and the interaction energies of , and with telomestatin alone, telomestatin + one G-quartet are given in supplementary Tables and , respectively.
4. Conclusion
We have therefore shown that telomestatin readily coordinates monovalent cations such as and , and that telomestatin retains this cation while binding to G-quadruplexes. The observed stoichiometry and the calculations are consistent with a cation trapped midway between the telomestatin and the G-quartet. Telomestatin therefore acts like an analog of a G-quartet. This study underlines that monovalent cation coordination capabilities should be integrated in the rational design of G-quadruplex binding ligands.
Acknowledgments
The authors wish it to be known that, in their opinion, the first two authors should be regarded as joint first authors. This work was supported by the Fonds de la Recherche Scientifique (FNRS) (FRFC 2.4.623.05 to Edwin De Pauw, CC.1.5.286.09.F to Valérie Gabelica, research associate position to Valérie Gabelica, postdoctoral fellowship to Frédéric Rosu, and logistician fellowship to Nicolas Smargiasso), and by the University of Liége (starting Grant FRSD-08/10 to V.G.).
Supplementary Materials
Comparison between three density functionals: overlay of structures (Figure S1), RMSD (Table S1) and interaction energies (Table S2).