The averaged values of *ΔH*, *T0*, *ΔS* at different concentrations of Na^{+}.

## Abstract

Determination of free energy of double helix formation from two single-stranded polynucleotides and estimation of energetics of different low-molecular compounds binding to nucleic acids provide valuable tools for understanding of mechanisms that govern noncovalent binding of ligands to their receptor targets. In order to completely understand the molecular forces that drive and stabilize double helix formation and its complexes with ligands, thermodynamic studies are needed to complement the structural data. Structural characterization of a number of DNA-ligand complexes by X-ray and high-resolution NMR method provides key insight relating to the properties of complex formation, but structural data alone, even when coupled with the most sophisticated current computational methods, cannot fully define the driving forces for binding inte-ractions (or interac-tions) or even accurately predict their binding affinities. Thermodynamics provides quantitative data of use in elucidating these driving forces and for evaluating and understanding at a deeper level the effects of substituent changes on binding affinity.

### Keywords

- free energy
- double helix
- helix-coil transition
- transition thermodynamic parameters
- DNA-ligand complexes
- binding parameters

## 1. Introduction

The 3D structure of solids by the change of environmental conditions may convert to a phase with quite different physical parameters describing the resulting state of matter. Transitions from one phase to another are accompanied by absorption or release of heat and sharply defined changes of energetic characteristics of the matter. At the fifties of last century, the biologically important molecules, nucleic acids and proteins, have been discovered — the structures of which were like one-dimensional linear aperiodic crystals [1]. The phase transition in linear crystals was theoretically treated at the twenties of last century [2]. According to this theory, the thermodynamic equilibrium is impossible for two homogeneous phases sharing common frontiers. Proper demonstration of the theorem efficacy was given much later, when the linear crystal to coil (helix-coli) transition of proteins and nucleic acid was investigated [3].

The unique feature of nucleic acid chains is their folding manner that encloses functional groups, i.e., purine and pyrimidine bases, so as to protect them inside a rigid and monotonous double-helix structure. At present, it is well established that DNA, the “major” molecule in the living cells, is polymorphous, and while functioning, the biopolymer may be in several forms: B-, A-, Z-, coil, etc., of which only Z-form was found to be a left-handed helix [4, 5]. There are two different types of structural transitions in DNA one of which (helix-coil, A-coil, Z-coil) is accompanied by unwinding of double helix (translation and replication, etc.). The second type of transitions (B-B^{1}, B-Z, B-A, A-Z, etc.) is realized by certain structural changes in sugar-phosphate backbone and base-pairs (bp) of DNA without unwinding the helix of the biopolymer. To understand the biological role of the existence of various forms of DNA, it is important to know the thermodynamic parameters of the phase transitions, particularly the value of free energy changed (*ΔF*), which is very difficult to obtain directly from the experiment. To estimate the *ΔF* value, enthalpy (*ΔH*) and entropy (*ΔS*) of transitions as usual are experimentally determined that are the constituents of free energy. We shall discuss below the experimental ways of estimating the values of these major thermodynamic parameters.

## 2. Main body

### 2.1. Theory

Along with genetic information realization in vivo (replication, transcription, translation), the molecule of DNA is being subjected to different conformational transitions. Moreover, there are no conformational transitions in “pure” molecule: it is always surrounded and interacts with huge number of various low-molecular compounds, which in turn, interacting with DNA, can stabilize or destabilize different conformational states of polymer molecule. To judge if this or other ligands stabilize or not different conformations of DNA, DNA conformational transition in the complex under any external factor inducing this transition should be studied (temperature, pH, chemical effect, etc.) and compared with the pure molecule transition.

Nowadays, it may be strictly established that these transitions (B-coil, A-coil, B-A, B-Z, Z-A, etc.) carry a cooperative character. The transition cooperativity is a direct consequence of the fact that the transition occurs in quasi-one-dimensional aperiodic crystal: in this case, the real phase transition is excluded.

Analysis of numerous experimental data, as well as some general representations about helix-coil transition, condition the possibility to formulate DNA main model, which is applied for theoretical observation of its melting. The model is sufficiently simple – DNA is one dimensional system that forms pairs of bases and each of them may be only in two states: helical and coil-like. Lengthening of the helical region per pair is accompanied by free energy value change Δ*F*. The value of Δ*F* determines the constant of this process:

where *R* is gas constant, *T* is temperature, and *σI* and *σ*_{i+1} are concentrations of molecules containing helical regions from i and i + 1 pairs of bases, respectively. In the transition point T_{0} σ =1, consequently, Δ*F* turns to zero. In the vicinity of this point, Δ*F* linearly depends on the temperature:

where *ΔH* and *ΔS* are changes of enthalpy and entropy, respectively.

Formation of new melted region in helical part is connected to appearing of additional boundaries between helical and melted regions and requires additional changes of free energy.

Value of *ΔF* determines the cooperativity of the system and

is called a cooperativity factor. If *F0* = 0, the cooperativity is absent. When *F0*→∞, the system is exposed to phase transition. At 0 < *F0*<∞, the transition carries a cooperative character and the higher is *F0*, the more favorable are long helical and melted regions and correspondingly the melting interval decreases.

Observed model, known in statistical physics as Ising model, physically corresponds to a case of single-stranded homopolymer. Let us observe this model applying the method of more probable distribution [6].

It is known that equilibrium values of physical magnitudes are corresponded to their most probable values at the given energy of the system. They can be found from the condition of “nonequilibrium” free-energy minimum:

where *W* is the number of states corresponding to the given energy of *E*.

Linear homopolymer consisting of *N* rings is observed. Each of these rings may be in one of these two states: melted – coil-like and helical. Macroscopic state of such system at the certain T temperature is given by three parameters: N_{2} is number of helical rings (in the second state), N_{1} is the number of rings in coil-like state (1), and n is the number of regions consisting of rings 1 or 2. It is clear that *N*_{1} + *N*_{2} = *N*; moreover, the case of infinite homopolymer is observed N→∞. If *F1* and *F2* are free energies of rings being in melted and helical states, respectively, *F0/2* is the free energy of boundary between helical and coil-like rings, the whole energy of the system is:

Number of microstates corresponding to given values of N_{1}, N_{2}, and n will be equal to

where *W2* and *W1* are numbers of modes by which helical and coil-like (melted) rings at the given values of N_{1}, N_{2}, and n may be distributed:

In the observed case in Eqs. (7) and (8), the unit can be neglected (N→∞). In this case,

Replacing (5) and (9) in (4) and applying Stirling’s formula, we will obtain

Equilibrium values of N_{1}, N_{2}, and n are determined from conditions

and

If to mark *ξ* = exp(−*Fσ*/*RT*) and *σ* = exp(*ΔF*/*RT*), where *ΔF* = *F*_{1} − *F*_{2} is a free energy change at helix-coil transition, from Eqs. (11) and (12), we will obtain

at the condition of total ring number constancy (*N*_{1} + *N*_{2} = *N*):

The obtained equations have dependences of N_{1}, N_{2}, and n on σ. Jointly solving Eqs. (13) and (14), the equation of ring part being in helical state θ = N_{2}/N_{1} is

The Eq. (15) describes the helix-coil transition curve.

From definition of the transition interval width, we will obtain:

where *ΔH* is enthalpy and T_{0} is the transition temperature.

T_{0} is determined from the condition that in transition point, the free energy change is equal to zero

From the Eq. (2), we will obtain

where *ΔS* is the difference of entropy in melted and helical states.

It should be mentioned that in the case of *F0* = 0, i.e., at the absence of interaction between rings (ξ=1), the formula (16) transmits to Boltzmann’s distribution. This case responds to cooperativity absence. At boundary energy increasing (decreasing of ξ) the system “becomes” cooperative; the melting interval decreases Eq. (16). In the threshold case when ξ→0, the system is entirely cooperative, but the transition is sharp: *ΔT*→0.

One of the fundamental predictions is that in the transition interval, polynucleotide chain is divided into alternate helical and coil-like regions, the length of which depends on the value of ξ. The average length of the helical region is equal to

In the transition point, θ = 1/2 and the average length of helical (and coil-like) region is equal to:

One of the first attempts to estimate the cooperativity factor value was presented in [7] by comparison of experimentally obtained value of ξ for homopolymer to the theory. It was shown that the value is in interval 10^{−4}–10^{−5}. For heteropolymer, the estimation is less precise since the melting interval width dependence on ξ in this case is logarithmic [6]. Uncertainty in values of *ΔT* depending on the cooperativity factor in the cases of different models shows that it is necessary to calculate and compare to experiment such characteristics of the helix-coil transition, which do not depend on ξ in wide change interval of this parameter.

Such invariant values are changes of melting temperature and melting interval width invoked by DNA binding to low-molecular compounds (ligands) [6, 8, 9, 10, 11]. From the point of view of the effect on DNA double-helix stability, ligands that are able to form complexes with polymers may be divided into stabilizers and destabilizers. Comparison of the melting curves of “pure” and ligands bound to >DNA can give information about the character of ligand binding to DNA: if the complex melting temperature (T_{m}) is higher than T_{0} for pure DNA, stabilization occurs, and if T_{m} decreases, then, destabilization occurs. Independently on the chosen model, molecules, possessing high affinity to double-helical polynucleotide, will stabilize the native structure and molecules, well binding to coil-like DNA — destabilize polymer double helix. What concerns to the melting interval, in both cases it increases as compared to that of pure polymer.

One of the predictions of the theory is that the melting interval width dependence on ligand concentration should have bell-like shape. It is explained by the fact that at small concentrations of ligands, ΔT of complexes increases due to the redistribution of ligands between helical and coil-like regions, which takes place during denaturation process with ligand concentration enhancement in accordance to their affinity to those regions. This redistribution results in additional stabilization of remained helical (or formed denatured) regions, and the melting process is extended. Due to confinement of number of the binding sites on DNA, the further increasing of concentration of ligands leads to difficulties of redistribution process and the melting interval width again decreases. In the boundary case when all binding sites are occupied by ligands, the melting interval width increment tends to zero. In the observed case, it is assumed that each pair of bases in polymer may be a binding site for ligand [12].

It is followed from the above-mentioned case that maximum of bell-like curve of the melting interval width increment dependence on ligand concentration corresponds to concentration of the ligand on DNA equal to half of the binding sites. The treated theory was compared with the experiment of complex melting, where as a ligand acridine dyes and actinomycin [13], native (destabilizer) and denatured (stabilizer), RNAase, heavy metal ions were used.

The effect of ligands on the helix-coil transition in polynucleotide in the case of random number of the binding sites has been studied. The chosen model in [14] is the following. We will assume that in solution, there are polymer molecules with fixed values of N_{1} (number of rings in coil-like state) and N_{2} (number of helical regions), the total number of rings N remains constant:

Let add ligands into solution with polymer that can bind both with coil-like and with native regions of DNA and can be in solution in nonbound state as well. If K_{2} and K_{1} are numbers of ligands bound to helical and coil-like regions, respectively, and K_{0} is the number of nonbound ligands, it is obvious that total number of ligands K per molecule satisfies the condition:

Let us mark the number of pairs of bases per binding site for denatured and native parts of the molecules as r_{1} and r_{2}, respectively. In this case, the number of binding sites for the respective regions will be equal to N_{1}/r_{1} and N_{2}/r_{2}. Taking this fact into consideration for nonequilibrium free energy, we will have:

where Ψ_{1} and Ψ_{2} are free energies of ligand bond with coil-like and helical parts of polymer, Ψ_{0} is the free energy, N_{0} is the number of binding sites in solution for free, nonbound ligand to polymer, S_{0}(N_{1}, N_{2}, n) is the entropic member bound to pure polynucleotide, W(N, K) function is determined by:

Taking into account the Eq. (24) and neglecting the unit (when N→∞ case is observed) for additional entropic member in (23) responsible for redistributing entropy of ligands, we will obtain:

It is obvious that the equation obtained from the condition (∂*G*/∂*n*) = 0 remains as it was in the absence of ligand. It means that the average length of helical region *ν*_{2} = *N*_{2}/*n* at given denaturation degree does not change when the ligand is added. In its turn, it means that ligand does not change the boundary energy. On the other hand, the equation obtained from ∂*G*/∂*N*_{1} = 0 condition does not change:

where c_{1} = K_{1}/N_{1}, c_{2} = K_{2}/N_{2} are concentrations of ligands for denatured and coil-like parts of polymer, respectively.

The developed theory gives dependencies of the experimentally observed transition parameters (the melting interval width *ΔT* and melting temperature *Tm*) on the binding parameters of a ligand with DNA (the binding constant *K* and the binding site *rq*) and the concentration of ligands [15].

where

where *T0* and *Δ0T* are the melting temperature and melting interval width for DNA in the absence of the ligand, *Tm* and *ΔT* are the same parameters for DNA-ligand complexes; *ΔH* is the enthalpy of the transition; *pq* = *Kq*/*K*_{1}, where *K1* is the binding constant for the first type (arbitrary chosen) of interaction of the ligand with one of the DNA forms; *Kq* is the binding constant for the q-th binding type (q = 2,…,m), which is expressed by the following equation:

*c*_{0} = *k*_{0}/*N*_{0} < < 1 is the concentration of unbound ligand, where *x0* is the equation solution at q = 0.5 (*q = Nf/N (f = 1, 2)* is the fraction of either forms of the polynucleotides (B, Z, A, coil, etc.) within the melting interval).

c = 2D/P, where D is the total concentration of ligand in solution and P is that of bases of DNA.

### 2.2. Experiment

DNA is a one-dimensional aperiodic crystal [1]. Therefore, as it was mentioned above, the true phase transition in such molecules could not occur. The two phases formed during the transition will tend to be mixed as continuously decreasing parts of the system. Such conversion is known as cooperative phase transition, two thermodynamic parameters of which are characterized by temperature of transition T_{0} and width of transition ΔT, on the contrary of real phase transition, which is realized at fixed temperature.

#### 2.2.1. Helix-coil transition

All nucleotides in the native state of DNA are in helix form, which has much lower free energy, i.e., high stability, than any other states that DNA assumes to be at room temperature and other ordinary physiological conditions. In the nonbound state, the nucleotide chain to which the nucleotide base pains are attached has freedom of motion. The bound or nonbound states may be classified in terms of “helix” and “coil” states, respectively, and the transition from one phase to another is called helix-coil transition or melting.

Unfolding of the double helix of DNA is produced as an effect of temperature (T), pH, ionic strength (μ), and denaturants [16, 17]. The process is accompanied by the transition of the bound state of bp to nonbound state, which propagates from more stable to less stable groups. It has been found experimentally that the transition of DNA occurs in a very sharp manner (the transition is highly co-operative), which is characterized by two physical parameters: the melting temperature, T_{0} and the width of transition, ΔT. The sharpness of the transition depends on the value of junction free energy, F_{j}. The true phase transition (transition of crystal structures) occurs only at the case of F_{j}→∞ (the junction energy is infinitely large). These parameters change for complexes of DNA with “low-weight” compounds (ligands). Integration of the structural, kinetic, and thermodynamic data of ligand-nucleic acids interaction is necessary to clearly understand the mechanisms of ligand-nucleic acid complex formation. Such investigations are much important to characterize the binding mode, sequence specificity, and understanding in detail designing new generation of drugs affecting the gene expression. Structural data obtained by X-ray crystallography and NMR for many drug-nucleic-acid complexes were successfully used for estimating ligands that attempt to correlate structure with binding affinity. It was established that upon binding, the ligands interact with substrate as a rigid compound, which is advantageous for revealing thermodynamic contribution from structural data [18]. Data of the structures of ligand-DNA complexes obtained by X-ray crystallography and NMR methods showed the more possible way to much ligand shape with the receptors of substrates and represent only one aspect of the complex formation. That is, the binding site will be occupied by ligand complemented it in terms of shape, charge, and other binding components [19], neglecting the energetic characteristics of binding process. So, the structural data alone cannot define the driving forces for binding and predicting the binding affinities. To understand the molecular mechanism and energetics of ligand-nucleic acid interaction, knowledge of thermodynamic parameters provide data elucidating the driving forces of complex formation process [20]. A complete thermodynamic profile for a system of interest requires determination of the free energy, enthalpy, and entropy.

The free energy ΔG is the key thermodynamics parameter, dictating the direction of biomolecular equilibria. If its sign is negative, the binding reaction or conformational transition will proceed spontaneously to an extent governed by the magnitude of ΔG. If its sign is positive, the magnitude of ΔG specifies the energy needed to drive the reaction to form product. The free energy is a balance between enthalpy and entropy. The enthalpy change reflects the amount of heat energy required for achievement a particular state, and the entropy measures how easily that energy might be distributed among various molecular energy levels. For binding reactions, negative enthalpy values are common (but not omnipresent), reflecting a tendency for the system to fall to lower energy levels by bond formation. Positive entropy values are common for binding reactions, reflecting a natural tendency for disruption of order. All binding reactions must overcome inescapable entropic penalties resulting from the loss of rotational and translational degrees of freedom.

The binding enthalpy (*ΔH*) can be detected using isothermal titration calorimetric (ITC) or differential scanning calorimetric (DSC) methods [18, 19, 20, 21]. The methods have several advantages for measuring binding energetic parameters at the same time having distinct difficulties, the dominant of which is high concentration of nucleic acids that require large quantities of expensive products, and besides, the possible aggregation makes very difficult to explain the experimental results [21, 22, 23]. DSC and ITC are laborious and time-consuming methods that often relegate calorimetric ones to be used as a secondary screening method. To overcome these limitations, several attempts have been made to improve the throughput of calorimetric and thermodynamic measurements. Mentioned difficulties for detecting the thermodynamic parameters of ligand-nucleic acid interaction may be overcome by applying methods, which are experimentally easy to perform, where very low concentrations of nucleic acids are used, which exclude the very unwanted process of aggregation [24].

The quantitative analyses of the effect of different substances (ligands) such as ions, antibiotics, dyes, proteins, etc. made it possible to suggest a simple method named “tie calorimetry” to estimate ΔH of conformational transitions [25, 26, 27, 28]. It has been shown that the enthalpy of helix-coil transition or melting (per base pair) could be determined from the experiments on DNA melting with ligands by the following general formula (32) (This formula is valid for all known types of conformational transition in one-dimensional crystals and for all types of ligands):

where *δTm* = *Tm* − *T*_{0} and *δ*Δ*T* = Δ*T* − Δ_{0}*T*, *Tm* and Δ*T* are the melting temperature and width of transition for DNA when the ligand is added to the solution, *T0* and Δ_{0}*T* are the same quantities for DNA without ligand, c = 2D/P is the total number of ligand molecules in solution (D) divided by the total number of DNA base pairs (2P), and R is the gas constant. Eq. (32) is absolutely general, and its validity does not depend on the values of the thermodynamics parameters of complexes such as binding constants of ligand with DNA, the number of binding sites on the biopolymer, etc. This was covered comprehensively in [25, 26, 27, 28]. On the other hand, obtained data showed that if the ligand complexes preferably bind with one of the conformations of DNA, the calculations become very simple [21, 22, 23], and for calculations, it is enough to compare the theoretical formula with the experiment either for *δTm* or *δ*Δ*T* [29]. Ethidium bromide, a very well-known ligand, binds preferably with the helix DNA [26, 27]. This enabled us at very low ligand concentration (c→0) with the combination of the area method [28, 29] to estimate ΔH for DNAs of two different GC contents at different Na^{+} concentrations with very high accuracy.

It was shown that at very small concentration of the ligand, the shift of the melting temperature (*δTm*) and widening of the melting curve (*δ*Δ*T*) are represented by the following equations:

For the coefficient A, the following formula is obtained:

where *r2* and *r1*, are the number of binding sites on the duplex and single-stranded DNA, respectively, *K2* and *K1* are the binding constants of ligand with helix and coil states of DNA, respectively, and P is the concentration of phosphate groups of DNA: p = K_{2}/K_{1}.

The only condition for validity of Eqs. (33) and (34) is c→0.

From Eq. (35), it follows that if

These formulas show that if *δΔT* is twice greater than *δTm*, the binding constant of ligand with one of the conformation of DNA is much greater than that of the other conformation. In this case, *ΔH* may be estimated on by Eqs. (34) or (35). The accuracy of *ΔH* value depends on the accuracy of experimental estimations of the *δΔT* and *δTm* values at different concentrations of ligand (different c). Therefore, the error is large (1.5–2 kcal/mol) when *ΔH* is calculated by Eq. (32) [28].

The accuracy of *ΔH* estimation is much higher if the “area” method is used for obtaining *δTm*. The method may be explained as following. The DNA melting temperature *T0* may be defined as the first moment of the differential melting curve (−dϑ/dT):

After integration, we have the following expression:

Here,

is numerically equal to the square limited by the melting curve *(1-ϑ),* the temperature axis, and the *T = TGC* vertical line. It follows from Eq. (39) that *T0* varies if the shape and place of melting curve change. Both the shape and the place of the melting curve change if ligand is added to the DNA solution. In this case, melting temperature of the complex *Tm* is found as:

where *Tm* is the temperature and (1-ϑ^{*}) is the melting curve of DNA-ligand complex.

It follows from Eqs. (39) and (40) that the variation in temperature can be expressed as:

where δs is the area limited by melting curves of DNA (left curve) and DNA-ligand complex (right curve) (Figure 1).

Substituting Eq. (43) to Eq. (36) for the enthalpy of helix-coil transition, one gets:

It should be noted that measuring of δs should be done at very small concentrations of ligand (c < 3⋅10^{−2}) [23], where *ΔH* is independent of the chosen concentrations of ligand. The dependence of *ΔH* (in kcal/mol) on Na^{+} is shown in Figure 2.

The values of *ΔH* obtained here agree excellently with calorimetric data [30]. The *ΔS* value may be calculated taking into account that at the transition mid-point (T_{0}), the free energies of the phases (for example, helix and coil) are equal to each other. Therefore,

Knowing the value of transition point *T0*, one can calculate the value of *ΔS* by equation

The averaged values of *ΔH*, *T0*, *ΔS* at different concentrations of Na^{+} are presented in Table 1.

Na^{+} M | −lgNa^{+} | T_{0} ^{°}C | Cl.perfr. ΔH kcal/mol | ΔS ent. unit | T_{0} ^{°}C | M. Lysod. ΔH kcal/mol | ΔS ent. unit | |
---|---|---|---|---|---|---|---|---|

1 | 10^{−3} | 3.0 | 45.1 | 7.2 ± 0.5 | 22.6 | 68.4 | 7.7 ± 0.4 | 22.5 |

5 | 10^{−3} | 2.3 | 56.3 | 8.2 ± 0.3 | 24.9 | 77.6 | 8.6 ± 0.3 | 24.5 |

1 | 10^{−2} | 2.0 | 61.0 | 8.7 ± 0.2 | 26.0 | 81.8 | 9.0 ± 0.2 | 25.4 |

5 | 10^{−1} | 1.0 | 77.3 | 9.8 ± 0.2 | 28.0 | 94.9 | 10.1 ± 0.2 | 27.4 |

#### 2.2.2. B-A transition

Investigations show that B-A conformational transition is cooperative, and it is realized in big amount of nucleotides conversion from B-form to A-form. Since the transition is independent of temperature and GC content of biopolymer, the Ising model has been employed for theoretical description of the process. The B-A transition initiated by any external factor a is represented by:

where *ϑ=f(a)* is the dependence of part of A-form on –a factor, s-form is the equilibrium constant of transition, *σ* = (exp − *ε*_{0}*T*), ε_{0} is the energy of junction, and T is the absolute temperature. Eq. (47) gives

where *Q* is a constant showing the steepness of free energy *ΔG* of A- and B-forms at the transition region. Knowing *Q* and Δa, the *ΔG* of A-B transition can be determined. It was shown that for pure water (100%) environment, *ΔGAB* = 1 *kcal*/*mol* [31, 32]. The obtained data coincide with the experimental results of [33], where the junctions of A- and B-forms are considered as a tie, which stabilizes the duplex.

#### 2.2.3. B-Z transition

Poly[d(G-C)] in a 55% ethanol solution exhibits the B-Z transition when the temperature increased [15, 33]. A polyamine, AEPDA, stabilizes Z-form and binds to it much stronger than to the B-form (Figure 3). Results show the temperature effect on the B-Z equilibrium without the polyamine (Figure 3a) and in its presence at a concentration of one molecule per 50 base pairs (Figure 3b). Obviously, the B-form of the polymer is stabilized by the rise of temperature in both cases. The pattern of CD spectra and the presence of a distinct isodichroic point at 301 nm show that only B- and Z-forms are involved in the equilibrium.

The transition profiles (Figure 4) show that, when the polyamine is added to the Z-form of poly[d(G-C)], the Z-form is stabilized and also the transition interval ΔT significantly widens. So, the polyamine is a “tie” for the Z-form. In this case, the transition enthalpy can be calculated by measuring the shift of the transition point (*δT0*) and the widening of the transition curve *δΔT* (see the Eqs. (36)–(38)).

Experiments showed that the ratio *δΔT/δT = 2,* which is independent on ionic strength [15, 34, 35, 36]. Thus, the *ΔH* value for the poly[d(G-C)] in 55% ethanol, *ΔHBZ* **= −**1.4 kcal/mol, is independent on the ionic strength.

Eqs. (37) and (38) are restricted to the only condition of c→0. K_{Z}> > 1, K_{B}> > 1, and at this case, we obtain Eqs. (6) and (7), which show that if the widening of the transition curve is twice as great as the shift of the transition point, *ΔHBZ* can be determined independently by each equation.

Figure 5 shows that for the polyamine, the ratio *δΔT/δT0* = 2; therefore, this ligand is perfect for thermodynamic investigations of a B-Z transition. Table 2 presents the data on the polyamine action calculated within the range of ionic strength form of 0.5–2 mM NaCl. These data show that *δT0* and *δΔT* values do not depend on NaCl concentration. Consequently, *ΔHbz* is independent of ionic strength under these conditions. Figure 5 shows the *δΔT* and *δT0* of a B-Z transition plotted as a function of c = 2D/P.

c = 2D/P | c = 0.01 | c = 0.02 | ||
---|---|---|---|---|

NaCl, M | δT_{0} | δΔT | δT_{0} | δΔT |

5 × 10^{−4} | 3.0 | 4.6 | 5.6 | 11.6 |

2.7 | 5.2 | |||

2.0 | 4.6 | 5.0 | 10.1 | |

2.6 | 4.8 | 5.3 | 10.8 | |

10^{−3} | 2.1 | 5.0 | ||

2.6 | 4.4 | |||

2.0 | 6.0 | |||

2.2 | 5.1 | |||

1.4 × 10^{−3} | 2.7 | 6.6 | ||

2 × 10^{−3} | 3.0 | 5.6 | ||

3.0 | 5.0 | |||

3.0 | 5.3 |

Since the slopes of the lines are related as 2:1, the enthalpy of the transition may be obtained using Eqs. (37) or (38):

We consider the values of *ΔHbz* obtained in this work reliable, which is in good agreement with the recently obtained same value of *ΔHBZ* using an independent method based on the polyelectrolyte theory of the B-Z transition [37].

#### 2.2.4. Thermodynamic parameters of binding: binding constants (K and binding site size r, the number of DNA base pairs corresponding to a binding site)

A traditional method of obtaining the interaction thermodynamic parameters is the Scatchard’s analysis of the ligands binding data, which consists of plotting the r/c_{f} value versus *r,* where *r* is the ratio of the bound ligand to DNA base pair concentration and c_{f} is the free ligand concentration [37]. This method has two major drawbacks. The first is the uncertainty in the *cf* value [37, 38], and the second is the existence of two different models of interaction of the ligands with DNA in the case of nonlinear Scatchard plots. One model assumes the presence of more than one type of independent binding sites, and the other model suggests interaction between bound ligands [34].

Our theory suggests another method for obtaining the binding parameters of the ligands interacting with DNA [15].

The binding parameters *(K* and r) are the parameters of the theory. They can be evaluated from comparison of the theory with experiment. The shape of the curves of dependencies of the inverse melting temperature *δ(l/Tm)* and of the melting range width *δ(ΔT/Tm2)* on the concentration of ligands is very different and sensitive to different values of r*:* the binding site size and p_{q} = K_{q}/K_{1}, where K_{1} is the binding constant for the first type (arbitrarily chosen) of interaction of the ligand with one of DNA forms, K_{q} is the binding constant for the q-th binding type (q = 2 …….m, m types of bending are considered). Figure 6 shows that the effect of the value of p_{q} is very significant on the shape of the dependence of δ(ΔT/T_{m}^{2}) on c. Therefore, the parameters may be determined, so as to provide the best fit between the theory and experiment. Major criteria for fitting are the position, the shape, and the size of maximum of the experimentally obtained *δ(ΔT/Tm2)* on **c** or *δ(l/Tm).*

We applied the conjugated gradient method for the theoretical analysis of the obtained experimental data of helix-coil transition of the complexes EtBr and AMD with DNA. The binding parameters were determined to provide the best fit between the calculated dependence of *δ(ΔT/T*_{m}^{2}) on c and the observed one (Figure 7). The analysis of the obtained data shows that EtBr and AMD may form at least five types of complexes of which three types with helix DNA and two types with coiled DNA at 10^{−2} M Na^{+} (Table 3). Another theoretical parameter is the ratio of the binding constants *pq = Kj*/*Kt* (q = 2,…,m). *Kq* values are readily calculated if one of the binding constant is known.

Experimental conditions | Ligand | Binding site size n | Binding constant Kq×10^{−3} | ||||||||
---|---|---|---|---|---|---|---|---|---|---|---|

Helix | Coil | Helix | Coil | ||||||||

ns | nw1 | nw2 | ns | nw1 | Ks | Kw1 | Kw2 | Ks | Kw1 | ||

10^{−2} M Na^{+} | EtBr and AMD | 6 | 3.3 | 2.5 | 5 | 1.5 | 2.7 | 0.07 | 0.06 | 0.07 | 0.013 |

2.2 × 10^{−3} M Na^{+} | EtBr | 6 | 3.3 | 2.5 | 5 | 1.5 | 5.0 | 0.6 | 0.4 | 0.05 | 0.04 |

The calculated values of K_{q} for EtBr and AMD are presented in Table 3. The values of n and K_{q} obtained agree with the values determined from independent experiments [15].

Index “s” corresponds to the “strong” binding mode and “w” to the “weak” binding mode.

## 3. Conclusion

In this work, it has been shown that “tie calorimetry” possesses a number of advantages. The measurements can be carried out in such concentrations that the intermolecular interactions and denaturation effect on medium pH are neglected. From the above mentioned, it is concluded that the helix-coil transition enthalpy can be calculated by the “tie calorimetry” according to the formula (32) and the only condition is that c << 1. On the other hand, determining the value of *ΔH*, the respective value of *ΔS* may be calculated by the formula (18). Calorimetry is a direct technique especially suitable when *ΔHbz* is high. By contrast, the “tie” calorimetry is most suitable at low *ΔHBZ* values. It follows from Eqs. (33), (34) or (36), (37) that the low *ΔHBZ* value results in a great change in the position of the transition point or in the transition width. Therefore, these methods are complementary.

Besides, it may be said that the “tie” calorimetric method is simple and very easy to be performed. It is absolute and no graduation is required for it. The method is based on the measuring of differential experimental values, which excluded systematic errors. Therefore, combination of both “area” and “tie” calorimetric methods makes it possible to establish the energetic parameters of phase transitions with very high accuracy.