Effect of G-Quadruplex Polymorphism on the Recognition of Telomeric DNA by a Metal Complex

The physiological role(s) played by G-quadruplexes renders these ‘non-canonical’ DNA secondary structures interesting new targets for therapeutic intervention. In particular, the search for ligands for selective recognition and stabilization of G-quadruplex arrangements has led to a number of novel targeted agents. An interesting approach is represented by the use of metal-complexes, their binding to DNA being modulated by ligand and metal ion nature, and by complex stoichiometry. In this work we characterized thermodynamically and stereochemically the interactions of a Ni(II) bis-phenanthroline derivative with telomeric G-quadruplex sequences using calorimetric, chiroptical and NMR techniques. We employed three strictly related sequences based on the human telomeric repeat, namely Tel22, Tel26 and wtTel26, which assume distinct conformations in potassium containing solutions. We were able to monitor specific enthalpy/entropy changes according to the structural features of the target telomeric sequence and to dissect the binding process into distinct events. Interestingly, temperature effects turned out to be prominent both in terms of binding stoichiometry and ΔH/ΔS contributions, while the final G-quadruplex-metal complex architecture tended to merge for the examined sequences. These results underline the critical choice of experimental conditions and DNA sequence for practical use of thermodynamic data in the rational development of effective G-quadruplex binders.


Introduction
The potential of nucleic acids to fold into non canonical secondary structures and the assessment of their role in regulating physiological processes is increasingly becoming object of interest [1][2][3]. Conformation is generally strictly related to DNA sequence but DNA-protein interactions can largely affect it [4,5]. Additionally, changes in ionic strength, salt composition, pH as well as interactions with small molecules can drive nucleic acid structural transitions and trap a given sequence in a specific energetically favored conformation.
A relevant example of highly ordered DNA structures is provided by G-quadruplexes, peculiar structural arrangements which can be assumed by G-rich sequences [6]. They derive from the overlapping of planar arrays in which four guanines are paired together through a network of Hoogsteen bonds and are further stabilized by the presence of a monocationic atom (mainly Na + or K + ) [7].
Genome-wide analysis reveals that G-rich sequences are not randomly distributed along the genome but are clustered in specific regions which mainly correspond to gene promoters [8,9]. This finding suggested the structural equilibrium duplex-quadruplex as an additional level of control of protein expression. Since often they encode for oncogenes, selective induction of Gquadruplex structures by small ligands has been investigated as a novel chemotherapeutic approach [10]. An additional G-enriched site is represented by telomeres which form the termini of the chromosomes and, in humans, are composed of repetitive TTAGGG sequences oriented 59 to 39 [11]. They are involved in maintaining chromosomal stability and genome integrity [12,13]. Interestingly, in over 85% of human cancer cells the cell ability to indefinitely replicate and become immortal is achieved through the activation of telomerase, a RNA-dependent DNA polymerase responsible for telomere elongation. Telomerase is inactive in somatic cells, thus it has become an attractive target for anticancer therapy. Among different possible approaches, we are particularly interested in indirect inhibition of telomerase activity [14,15]. Indeed, the peculiar G-quadruplex conformational arrangement that the G-rich telomeric sequence can assume is not recognized by the enzyme which processes only a single stranded template. For this reason, small molecules able to induce and stabilize the G-quadruplex form of the telomeric sequence can prevent the hybridization of the telomerase RNA template onto the primer and thus inhibit the enzymatic activity [16]. Accordingly, G-quadruplex induction has been confirmed to stimulate cellular senescence, apoptosis or autophagy [17][18][19]. Additionally, modification of the telomere structural equilibria can further alter the recognition by telomere-directed proteins leading to short term cytotoxic effects [15].
Up-to date, most of the reported G-quadruplex binders are ligands structurally based on extended planar aromatic arrays in which the p-delocalized system allows stacking interactions with the external guanine tetrads. Generally, ligand scaffolds are further functionalised with positively charged side chains that enhance the ligand affinity by interacting with the negative phosphates of the DNA backbone [20,21].
Although the majority of these compounds are organic systems, recently metal ion complexes were elegantly used as potential Gquadruplex binders and stabilizers [22][23][24][25]. In these complexes, the metal center binds to specific ligands according to the electronic configuration of the species involved, which ultimately produces specific geometries around the coordination sphere [26,27]. Additionally, the coordination of the metal to aromatic ligands (such as canonical G-quadruplex ligands) can withdraw the electron density yielding an electron-deficient system with increased p-p interacting capabilities towards the G-quartet planes. As a result, the G-quadruplex stabilization properties of the metal complex are distinct from those of the ligand itself due to the charged metal center and the stoichiometry/stereochemistry of the metal complex arrangement.
Phenanthroline is a known efficient ligand for several metal ions able to form metal complexes of different geometries. In particular, extensive literature data describes the ability of the phenanthroline moiety to interact with duplex DNA with characteristic affinity and binding mode when bound to metal ions [26].
In previous works, we exported this model toward Gquadruplex DNA with the idea that only selected metal ions could coordinate phenanthroline units according to geometries suitable to enhance the G-quadruplex recognition and disfavor the binding to the double helix [28,29]. In agreement with other authors [30], our investigation showed Cu(II) and Ni(II) as the most effective transition metals in providing phenanthroline complexes selective toward G-quadruplex [28,29]. Moreover, the evaluation of several phenanthroline-based free ligands and of their Ni(II) complexes indicated the assembly of two phenanthroline moieties around the metal center as the key structural feature required to enhance such a selectivity. Based on these results, we proposed a model according to which the extended planar surface achieved through metal coordination promotes stacking onto a guanine quartet. In particular, we suggested the terminal G-tetrads as the potential DNA ligand binding sites.
To investigate in more details this binding mode, we decided to perform a thermodynamic characterization of the interactions of these metal complexes with the telomeric G-quadruplex. In particular, here we focus on the Ni(II) complex formed by the phenathroline derivative K34, (K34) 2 Ni(II) (Figure 1). Indeed this metal complex has physical and chemical properties, such as water solubility and thermal stability, which make it a good candidate for biologically relevant studies.
Our previous target sequence was the telomeric sequence Tel22. However, in potassium containing solution, its G-quadruplex form is highly polymorphic [31]. This property can affect the quantitative description of the ligand binding process. Distinctly, NMR studies reveal that two related sequences wtTel26, and Tel26 in the same experimental conditions assume prevalent well-defined different hybrid forms [32,33]. Noteworthy, these sequences comprise the reference sequence Tel22 but they contain additional 59 and 39-flanking dinucleotides which are involved in pairing and stacking interactions thus providing capping structures which lock DNA preferentially in one conformation [34]. For this reason, we performed calorimetric and circular dichroism titrations using all these three sequences based on the telomeric one: the reference Tel22, wtTel26 and Tel26.
Using this approach we were able not only to monitor different ligand-DNA affinities according to the structural features of the target telomeric sequence but also to dissect the binding process into distinct events.
Before each experiment, G-quadruplex forming sequences were dissolved in 10 mM Tris, 20 mM KCl at pH 7.5, extensively washed with the same buffer using Millipore Filter Units (MWCO: 3000), diluted to the desired concentration, heated at 95uC for 5 minutes and slowly cooled to room temperature to insure folding of the sequence. Oligonucleotides concentration was determined by UV-Vis from absorbance at 260 nm using their reported molar extinction coefficient.
Ligand, metal ion and metal complex solutions K34 was synthesized as previously reported [29]. Stock solutions (4 mM) were prepared in 10 mM Tris, 20 mM KCl at pH 7.5 and diluted to the required concentrations with the same buffer. NiCl 2 (0.8 M) was dissolved in deionized water and metal ion concentration was determined by ICP (Inductively Coupled Plasma, Optima 3000 DV Perkin Elmer). This solution was further diluted in 10 mM Tris, 20 mM KCl at pH 7.5. Solutions of (K34) 2 Ni(II) were prepared by mixing the required volumes of ligand and metal ion solutions.

Circular Dichroism Studies
Circular dichroism spectra were recorded on a Jasco J-810 spectropolarimeter equipped with a Peltier temperature controller in 10 mM Tris, 20 mM KCl at pH 7.5 using a 10 mm pathlength cell. Before data acquisition, DNA solutions (4 mM strand concentration) were heated at 95uC for 5 min and left to cool down at room temperature over night (o.n.). Each reported spectrum represents the average of 2 scans recorded with 1-nm step resolution. Observed CD signals were converted to mean residue ellipticity [H] = deg x cm 2 6 dmol 21 (Mol. Ellip.).
For scT22, the saturation fraction f (DCD/DCD max ) was analyzed according to equation (1) where L o and R o are the ligand and DNA concentration expressed in residues, respectively, and n is the complex stoichiometry. Thermal denaturation experiments were performed by recording the optical signal while increasing the temperature at 0.8uC/ min. The melted solution was then cooled down at the same temperature change rate to check for hysteresis. Melting temperatures (Tm) were calculated from the first derivatives of the melting profiles. Each curve was repeated at least three times and errors were60.4uC. DTm were calculated by subtracting the Tm value recorded in the presence of the ligand from the corresponding value in the absence of ligand.

Isothermal Titration Calorimetry (ITC)
ITC titrations were performed on a MicroCal VP ITC instrument in 10 mM Tris, 20 mM KCl at pH 7.5 at 25uC and 37uC. Working solutions were degassed for 5 minutes prior the use. Volumes of 10 ml of K34 (1 mM) or (K34) 2 Ni(II) (0.5 mM) were injected into a solution of previously folded DNA (25 mM). The ITC titration settings were: injection volume 10 ml, spacing between ligand injection 360 s, injection time 10 s, stirring speed 345 rpm, equilibration time 60 s.
Before data analysis, raw data were corrected for the heat of dilution. Heats were integrated and binding parameters were calculated according to one or two binding site model using Origin Software.
Data analysis provides DH (reaction enthalpy change, kcal. mol 21 ), K a (binding constant, M 21 ), and n (number of bound ligands) whereas the Gibbs energy and the entropic contribution were calculated using the relationships DG = 2RT ln K a and DG = DH 2 TDS, respectively.
For the interaction of our metal complex with Tel26 at 25uC the best fitting results were obtained by a two sequential binding sites model always. According to it the number of sequential sites must be exactly integral and thus it is held constant during the fitting procedure.

Nuclear Magnetic Resonance (NMR)
NMR experiments were performed at 25 uC on a Bruker DMX 600 spectrometer, equipped with a 5 mm TXI probe with gradients, and the data were processed using the TOPSPIN 2.0 software.
1H one-dimensional spectra were acquired with 32 k scans, 1.52 s acquisition time, 1 s relaxation delay, 18 ppm spectral width.
Suppression of the water signal was achieved using the WATERGATE sequence before acquisition. For proton assignments see [32,33].

ITC Based Thermodynamic Studies
To obtain information on the binding of (K34) 2 Ni(II) to the Gquadruplex telomeric sequence we performed Isothermal Titration Calorimetry (ITC) analysis which provides a direct evaluation of thermodynamic parameters relevant to describe biomolecular interactions [35,36]. As previously reported, this metal complex was shown to be stable in our working conditions [29] and did not tend to dissociate as evidenced by spectroscopic titrations ( Figure  S1). To evaluate the contribution of DNA structural arrangements in the binding process the analysis was performed with Tel22 and, in addition, with Tel26 and wtTel26 which preferentially assume a Hybrid 1 and Hybrid 2 folding, respectively [32][33][34]. These two G-quadruplex forms share common structural elements such as the overlapping of three G-tetrads and a mixed parallel/ antiparallel orientation of the four strands. However, they differ for the loops arrangement and the relative strand orientation. Moreover, specific capping structures are formed. In particular, a T:A:T triple capping is present in the Hybrid 2 structure while an A:T base pair and an adenine triple capping are found in the Hybrid 1 arrangement [32,33].
In our experimental conditions, addition of (K34) 2 Ni(II) to any tested G-quadruplex folded sequences resulted in heat release which indicates the occurrence of an exothermic process (Figure 2).
At 25uC, data collected in the presence of Tel22 were poorly reproducible, thus, unfortunately, they could not be used for a safe comparison. Nevertheless data acquired with Tel26 and wtTel26 underlined interesting features of the (K34) 2 Ni(II) binding process which share several analogies but also quite peculiar differences as a function of the sequence of the nucleic acid target.
Recognition of the wtTel26 sequence (largely folded in a Hybrid 2 arrangement), showed a sigmoidal profile consistent with a single set of binding sites (Figure 2A). Data were analyzed accordingly, thus assuming any potential site of interaction as equivalent and independent. Such an analysis fits well the experimental data and indicates the binding of two metal complexes per G-quadruplex molecule. The resulting best-fit thermodynamic parameters are reported in Table 1.
The calorimetric titration of Tel26 (folded in a Hybrid 1 arrangement) with the same ligand (K34) 2 Ni(II) provided a binding isotherm, derived from the integrated heat data, which corresponds to two distinct binding events ( Figure 2B). Data analysis indicated that the best fitting was obtained using a two sequential binding model. In particular the experimentally recorded profile is well described by two sequential processes, each resultant from the interaction of one metal complex to one DNA target. Thus, a final stoichiometry of two (K34) 2 Ni(II) molecules associated to one G-quadruplex template was confirmed. Thermodynamic parameters (Table 1) describing the first binding event showed a Ka of 8.66*10 5 M 21 which is higher, albeit of the same order of magnitude, in comparison to the one observed with wtTel26. The binding of the second metal complex is characterized by a lower binding constant associated to a larger favourable enthalpy change accompanied by a negative DS value.
By increasing the working temperature up to 37uC, the binding model of (K34) 2 Ni(II) towards wtTel26 is essentially conserved ( Figure S2A). Although two temperature data points are not sufficient to properly analyze the d(DH)/dT relationship, the negative temperature dependence of DH points to a negative DCp, which can be related to a hydrophobic effect or to the occurrence of temperature dependent conformational equilibria [37,38]. Interestingly, the binding affinity and, consequently, the resulting negative free energy are not largely affected by changes in the working temperature (Table 2). If we analyze the binding process in terms of enthalpic and entropic contributions we note that both of them are favourable for the binding at 25uC. However, the entropic term goes from positive (17.2 cal/mol at 25uC) to close to zero (0.4 cal/mol) at 37uC. This results in a shift from a binding process driven preferentially by entropic contributions to an enthalpically promoted one. (1) Remarkably, this binding profile is generally shared also by Tel26 when titrated with (K34) 2 Ni(II) at 37uC. Indeed, the two binding events described at 25uC are no more clearly dissectible when the titration is performed at higher temperatures (Table 2, Figure S2B). A final stoichiometry close to 2:1 (the small increase over 2 possibly reflects tendency to ligand aggregation) was preserved but it corresponds to the recognition of two equivalent binding sites. In particular, a comparison of the binding constants at 25uC and 37uC suggests that the lower affinity event recorded at lower temperature is favoured by increasing the working temperature. Indeed, the binding constant for the two equivalent binding events at 37uC, is comparable to the one associated to the high affinity event at 25uC. This can be explained by an incremented accessibility to the second binding site.
As a result, a more effective recognition of (K34) 2 Ni(II) for the Hybrid 1 arrangement was evidenced although with a conserved progressive shift to an entropically less favorable reaction by increasing the working temperature (Figure 3).
At temperatures above 25uC we were able to examine also the Tel22 sequence ( Figure S2C). As already seen with Tel26 and wtTel26, the binding process monitored at 37uC was found to be exothermic and could be described by the binding of the metal complex to two equivalent sites. In this case, a stoichiometry lower than 2 might reflect ligand-mediated G-quadruplex-G-quadruplex  stacking. Interestingly, Ka value is almost one order of magnitude higher than those observed with the other two tested Gquadruplex templates ( Table 2). It is useful to underline that this reflects a relevant entropic term which, associated to the negative DH, favourably contributes to the binding process ( Figure 3).
To assess the peculiarities of the binding of our metal complex towards G-quadruplex structures, we performed also ITC titrations using the ligand K34 in the absence of the Ni(II) ion ( Figure S3). They showed a remarkable suppression of the heat release thus confirming the metal complex as the preferential binder for G-quadruplex structures in agreement with our previously reported data [29].

Circular Dichroism conformational studies
The calorimetric data clearly indicated that the folding of the analyzed targets highly influences the G-quadruplex binding process of (K34) 2 Ni(II). The capping residues, which are peculiar of the longer tested sequences, might be responsible for such a variation. Indeed, they can represent distinct binding sites. Alternatively, the (K34) 2 Ni(II)-G-quadruplex binding sites can be conserved among the three tested structures and in this case a common structural domain such as the tetrads would represent a reasonable target. As a result, the observed modulation in thermodynamic parameters can at least partly be connected to a different propensity of DNA to undergo structural modifications upon binding. In this connection CD titrations were performed to monitor conformational changes occurring upon ligand-macromolecule binding.
The recorded dichroic spectra of all tested oligonucleotides are characterized by two positive bands deriving from the 3+1 arrangement, one centered at 290 nm and the other at 265-268 nm. However, they mainly differ in terms of relative intensity and resolution that reflect their distinct folding (Figure 4).
Upon addition of (K34) 2 Ni(II), modifications in the CD spectra of all tested oligonucleotides are detectable even at low ligand versus G-quadruplex ratios (Figures 4 and 5). In agreement with the calorimetric data the signal intensity increased according to one hyperbolic curve and these changes reached saturation at a stoichiometric ratio close to 2:1.
Due to the spectral properties of tested DNA and of (K34) 2 Ni(II), it is not possible to attribute the observed spectral changes exclusively to the G-quadruplex structural rearrangement. Indeed, overlapping of DNA CD transitions and induced CD of bound metal complex can occur. This prevents a safe quantitative analysis of the binding isotherms derived from the optical signal variations. However, our data support a less efficient process of binding and structural rearrangement for the wtTel26 sequence ( Figure 5).
Interestingly, at ligand saturation, CD spectra appear to converge toward a similar signature, characterized by a positive band at 280 nm with a shoulder at 295 nm. This signature is peculiar of the G-quadruplex bound form. Indeed, likely due to its positive charge, (K34) 2 Ni(II) can also interact with different DNA structures. However, the resulting CD spectra of the metal complex bound to non G-quadruplex forming DNAs are quite distinct. As an example, data referring to the interaction of (K34) 2 Ni(II) with an oligonucleotide sequence of the same length and base composition as Tel22 but unable to fold into Gquadruplex is reported in Figure 4D. In this case, a binding constant Ka 2.7 * 10 5 M 21 , which involves two bases per metal complex, was observed. A direct comparison of this value with those obtained by ITC in the presence of G-quadruplex folded substrates is not safe. However, as clearly evidenced in Figure 5, the linearly arranged DNA is recognized less efficiently by the tested metal complex.

Induction and stabilization of G-quadruplex structure
The CD spectral changes associated to the addition of (K34) 2 Ni(II) to G-quadruplex folded telomeric sequences were conserved also when titrations were performed at temperatures above 25uC (37uC and 45uC) ( Figure 6). Data referring to the highest temperature are rather interesting. Indeed, the melting profiles obtained while recording the CD signal at 290 nm, indicate that the two longer oligonucleotides, wtTel26 and Tel26, irrespectively of their preferential folding into a unique defined structure, are thermally less stable than Tel22. In particular, their melting temperatures are 44.2uC and 47.2uC, respectively, thus indicating that in solution at 45uC they are partially unfolded. Nevertheless, (K34) 2 Ni(II) was able to produce the same Gquadruplex-bound CD spectra as those recorded at lower temperatures. This suggests that the metal complex is able to shift the DNA structural equilibrium towards a final folded form, which corresponds to the one observed at 25uC. Additionally, this form is stable at physiological temperature. Indeed, melting profiles of (K34) 2 Ni(II):G-quadruplex (molar ratio 2:1) recorded at 290 nm confirmed a single melting transition shifted to 61.0uC (DTm = 13.8uC) using Tel26 and to 59.4uC (DTm = 15.2uC) using wtTel26. Distinctly, in the presence of the more thermally stable Tel22 (Tm = 56.2uC), a thermal shift to 66.6uC (DTm = 10.4uC) was observed upon metal complex binding. This thermal stabilization is not unexpected since (K34) 2 Ni(II) interacts with G-quadruplex. However, the herein reported increments of DNA melting temperature further sustained the preferential binding of our metal complex to the G-quadruplex folded forms in comparison to the unfolded one.

Discussion
The thermodynamic signature for the binding of several small molecules to DNA is extensively investigated since it can help understanding the features of the recognition process [39,40]. However, distinctly from results reported for double stranded DNA, literature data concerning G-quadruplex binders evidenced unpredictable behaviours which do not allow correlating enthalpic/entropic contributions to binding mode [41,42]. A rational for such variability rests on the polymorphic nature of G-quadruplexes, which is extremely sensitive to buffer composition and DNA sequence. Additionally, different complexes can be formed depending upon ligand nature. As an example, even for a short model sequence like Tel22, although a hybrid-type folding appears to be predominant in potassium containing solutions, the coexistence of different conformations in mutual equilibrium has to be considered [43]. Moreover, different arrangements can be separated by small energy barriers (as in the case of basket-hybrid transition) thus allowing an easy shuffling from one form to the other [44]. To overcome these difficulties, in addition to the reference telomeric sequence Tel22 we used two related sequences, wtTel26 and Tel26, which prevalently fold into two distinct hybrid-type structures. Although only about 75% of the wtTel26 sequence assumes an Hybrid 2 fold whereas essentially Tel26 is fully present in solution in a Hybrid 1 conformation, these two structures do not easily interconvert [32,33]. Thus, due to their distinct three-dimensional folding they actually represent two structurally different targets for a ligand. In any event it is important to remind that the herein reported thermodynamic analysis represents the sum of multiple contributions, including both binding and structural rearrangement.
The analysis of the recognition process showed that (K34) 2 Ni(II) interacts efficiently with all the tested sequences according to a binding path characterized by a negative DH. Additionally, only modest changes in the binding free energy were observed within all the tested templates.
These common hallmarks do not exclude the interaction of (K34) 2 Ni(II) with different portions of the tested G-quadruplex structures. However, the binding stoichiometry provided by ITC and CD analysis indicates two (K34) 2 Ni(II) molecules per Gquadruplex structure. Additionally, dichroic studies showed that, irrespectively of the starting DNA folding, the overall final Gquadruplex-ligand complexes seem to share significant analogies. As above pointed out no DNA structural information can be safely derived by CD spectroscopy. To better dissect the issue of structural similarities among the two complexes, 600 MHz NMR studies were performed ( Figure S4). Unfortunately, this powerful technique could be only partly useful due to the large quenching and broadening of the signals due to the presence of the paramagnetic metal center. Nevertheless, comparison of the NMR spectra in the imino protons region clearly suggests that the Tel26 and wtTel26 bound forms share several common peaks ( Figure S4B). Interestingly, differences are related to peaks (located at 12.00, 11.05, 10.9 and 11.45, 11,4, 10.55 ppm for H1 and H2, respectively) which result conserved in the free and bound forms.
These peaks generally correspond to the imino proton signals relative to the guanines (G11, G17, G6 and G17, G23, G18 for H1 and H2, respectively) located in the central tetrad facing the parallel moiety of the oligonucleotide structure [32,33].
These results suggest that the binding modes might be conserved among the three tested DNA sequences. In particular, a recognition process which takes place at the terminal G-quartets where p-p stacking interactions can occur between the assembled phenathroline moieties and the exposed aromatic portion of the bases can be inferred.
This shared binding mode can actually help justifying also the differences observed in the binding process. These are: N a binding pathway which is template dependent. Indeed, at 25 uC, using the Hybrid 1 folded substrate, the binding process occurs through two sequential steps (on two independent sites of the G-quadruplex structure) whereas two equivalent binding sites are identified on the Hybrid 2 structure. N a higher affinity showed by (K34) 2 Ni(II) for Tel22 in comparison to the longer Tel26 sequences.
It is known from NMR studies that the two tested 26-mer telomeric sequences are selectively stabilized in a defined hybrid folding by the presence of capping structures at the top and bottom of the G-quadruplex core which stacks on the terminal tetrads. Thus it is feasible that these domains impair the accessibility of the metal complex to the target tetrads. Additionally, since these cappings are distinct in Tel26 as compared to wtTel26, it is not odd to assume that the formally conserved binding sites (the G-tetrads) can be seen as similar (as in the case of wtTel26) or different (as in the case of Tel26) by (K34) 2 Ni(II).
Interestingly, in this instance, it appears that a modest increment in the working temperature (thus a modest increment in DNA flexibility) is sufficient to level off such a difference. This points to the need for (K34) 2 Ni(II) to compete with the capping moieties above and below the G-tetrad core to gain access to the planar surface. On the opposite, Tel22, which lacks interfering terminal residues, does not require uncapping to grant exposure of the tetraplex surface, thus exhibiting higher affinity for the metal complex.
By keeping in mind that for Tel26 and wtTel26the capping moieties are actually the structural elements required to select  defined topologies among those assumed in solution by Tel22, it turns out reasonable that capping displacement by the ligand can allow the oligonucleotide to rearrange to a conformation best fit for optimal binding.
A final comment is deserved by the entropic contribution to the binding process. Only a limited number of studies deal with the interaction of metal complexes with G-quadruplexes [45][46][47]. In all instances a negative DH associated to a predominantly favorable entropic contribution was reported. Due to the charged nature of metal complexes, this likely reflects significant DNA counterions release produced upon binding in addition to solvation and hydrogen bonding effects. However, with our tested sequences, upon increasing the working temperature, an enthalpy-entropy compensation occurred (more negative DH, less positive DS), a well documented behaviour typical of ligandreceptor binding [48]. Interestingly, this effect is more pronounced when working close to the G-quadruplex melting temperature ( Figure S5 and Table S1). Since our metal complex is actually able to induce refolding of the melted oligonucleotides, the compensation effect may be reinforced by the introduction of structural constraints in the nucleic acid upon interaction with (K34) 2 Ni(II).
In conclusion, our study demonstrated how different Gquadruplex structures, produced by recurrence of a G-quadruplex forming sequence in closely related oligonucleotides, do not grant unique quantitative information on ligand binding. Consequently, an anticipation of possible pharmacological implications of the recognition process is very hard to adequately rationalize, even not considering metabolic effects.
Here, we show how the binding process can be thermodynamically affected not only in terms of enthalpy and entropy changes but also in terms of binding site equivalence. We have seen that the uncapped Tel22 can bind to a ligand with higher affinity than the structurally constrained 26-mer sequences, that increasing temperature leads to similar binding strength notwithstanding the exothermic nature of the process, that individual binding sites can be turned from similar to dissimilar on the same template for the same ligand just by slightly adjusting temperature conditions. Combination of these events allows distinct G-quadruplex template structures to adopt similar arrangements when in complex with an effective small molecule binder.
The more we learn, the less simple is foreseeing thermodynamics and stereochemistry of induced fit effects and the consequent biological/pharmacological implications of an apparently straightforward process such as ligand G-quadruplex recognition.