The averaged values of * ΔH*,

*,*T

_{0}

*at different concentrations of Na*ΔS

^{+}.

Open access peer-reviewed chapter

Submitted: 21 June 2017 Reviewed: 27 September 2017 Published: 18 July 2018

DOI: 10.5772/intechopen.71313

From the Edited Volume

Edited by Juan Carlos Moreno-Piraján

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.

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

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

Advertisement
## 2. Main body

### 2.1. Theory

σ = σ i + 1 σ i = exp Δ F / RT E1Δ F = Δ H + T Δ S E2σ = exp − F 0 / 2 RT E3F = E − T ln W E4F = F 1 N 1 + F 2 N 2 + F 0 n E5W = W 1 − W 2 E6W 1 = N 1 − 1 ! n − 1 ! N 1 − n ! E7W 2 = N 2 − 1 ! n − 1 ! N 2 − n ! E8W N 1 N 2 n = N 1 ! N 2 ! n ! N 1 − n ! n ! N 2 − n ! E9F = F 1 N 1 + F 2 N 2 + F 0 n − T N 1 ln N 1 − N 1 − n ln N 1 − n + N 2 ln N 2 − N 2 − n ln N 2 − n − 2 n ln n E10∂ F ∂ n N 1 N 2 = 0 E11∂ F ∂ N 1 n = 0 E121 ξ = N 1 n − 1 N 2 n − 1 E13σ = 1 − n N 2 1 − n N 1 E141 − 2 θ 1 − θ θ = 1 ξ ⋅ 1 − σ σ E15Δ T = 4 ξ T 0 2 Δ H E16Δ F = 0 E17T 0 = Δ H Δ S E18ν 2 = N 2 n = 1 ξ ⋅ θ 1 − θ E19ν 0 = 1 ξ E20N 1 + N 2 = N E21K 1 + K 2 + K 0 = K E22F = F 1 N 1 + F 2 N 2 + F 0 n + ψ 1 K 1 + ψ 2 K 2 + ψ 0 K 0 − TS 0 + ln W 1 N 1 r 1 K 1 + ln W 2 N 2 r 2 K 2 + ln W 0 N 0 K 0 E23W N K = N − 1 ! K − 1 ! N − K ! E24S add . = N 1 r 1 ! N 2 r 2 ! N 0 ! K 1 ! N 1 r 1 − K 1 ! N 2 r 2 − K 2 ! K 2 ! N 0 − K 0 ! K 0 ! E251 − n N 2 1 − n N 1 = σ 1 − c 1 r 1 1 / r 1 1 − c 2 r 2 1 / r 2 E26δ 1 T m = 1 Δ H ln ∏ j = l + 1 m x 0 + p j 1 / r j ∏ i = 1 l x 0 + p i 1 / r i x 0 ∑ i = 1 l 1 r i − ∑ j = l + 1 m 1 r j E27δ Δ T T m 2 = 1 x 0 Δ H ∑ j = l + 1 m p j r j x 0 + p j − ∑ i = 1 l p i r i x 0 + p i ∂ x ∂ θ θ = 1 / 2 E28δ 1 T m = 1 T m − 1 T 0 δ Δ T T m 2 = Δ T T m 2 − Δ 0 T T 0 2 E29K q = c q ⋅ r q c 0 1 − c q ⋅ r q , q = 1 … m E301 − θ ∑ i = 1 l p i r i x + p i + θ ∑ j = l + 1 m p j r j x + p j = c E31### 2.2. Experiment

#### 2.2.1. Helix-coil transition

Δ H = R ⋅ lim c → 0 δ Δ T / δ T m 2 T 0 2 ⋅ c , E32δ T m = Λ RT 0 2 c Δ H E33δ Δ T = Λ 2 RT 0 2 c Δ H E34Λ = 2 r 1 / r 2 p − 1 r 1 / r 2 p + 1 ⋅ K 1 / r 1 P + K 2 / r 2 P 4 + K 1 / r 1 P + K 2 / r 2 P E35K 1 / r 1 P < < 4 K 2 / r 2 P > > 4 E36δ T m = 2 RT 0 2 c ΔH E37δΔT = 4 RT 0 2 c ΔH . E38T 0 = ∫ T AT T GC − dϑ dT TdT E39T 0 = T GC − ∫ T AT T GC 1 − ϑ dT E40s = ∫ T AT T GC 1 − ϑ dT E41T m = T GC − ∫ T AT T GC 1 − ϑ ∗ dT , E42δT = δs = ∫ T AT T GC 1 − ϑ dT − ∫ T AT T GC 1 − ϑ ∗ dT , E43Δ H = 2 RT 0 2 c δs E44Δ G = Δ H − T 0 Δ S = 0 . E45Δ S = Δ H T 0 . E46

#### 2.2.2. B-A transition

1 − 2 ϑ 1 − ϑ ϑ = 1 − s s σ AB E47Δ a = ∂ ϑ ∂ a a = a 0 − 1 = 4 Q ν 0 E48#### 2.2.3. B-Z transition

Δ H BZ = − 2 RT 0 2 δ T 0 ⋅ c = − 1.4 ± 0.2 kcal / mol Δ H BZ = − 4 RT 0 2 δ Δ H ⋅ c = − 1.4 ± 0.2 kcal / mol E49#### 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)

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 Δ

where * R* is gas constant,

where * ΔH* and

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 * F* = 0, the cooperativity is absent. When

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

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

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

where * W* and

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(−

at the condition of total ring number constancy (_{1} + _{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

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 * F* = 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:

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*/∂

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

where

where * T* and

_{0} = _{0}/_{0} < < 1 is the concentration of unbound ligand, where * x* is the equation solution at q = 0.5 (

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

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.

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 * δT* =

It was shown that at very small concentration of the ligand, the shift of the melting temperature (* δT*) and widening of the melting curve (

For the coefficient A, the following formula is obtained:

where * r* and

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

The accuracy of * ΔH* estimation is much higher if the “area” method is used for obtaining

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

where * T* is the temperature and (1-ϑ

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

The values of * ΔH* obtained here agree excellently with calorimetric data [30]. The

Knowing the value of transition point * T*, one can calculate the value of

The averaged values of * ΔH*,

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 |

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,

where * Q* is a constant showing the steepness of free energy

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 (* δT*) and the widening of the transition curve

Experiments showed that the ratio * δΔT/δT = 2,* which is independent on ionic strength [15, 34, 35, 36]. Thus, the

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, * ΔH* can be determined independently by each equation.

Figure 5 shows that for the polyamine, the ratio * δΔT/δT* = 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

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 * ΔH* obtained in this work reliable, which is in good agreement with the recently obtained same value of

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

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

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 _{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 * p*/

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

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

_{s} | _{w1} | _{w2} | _{s} | _{w1} | _{s} | _{w1} | _{w2} | _{s} | _{w1} | ||

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.

Advertisement
## 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

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.

- 1.
Schrodinger E. What Is Life? The Physical Aspect of the Living Cell. Dublin: Cambridge University Press; 1944. 194 p - 2.
Landau LD, Lifshits EM. Statistic Physics. Moscow: Nauka; 1964 - 3.
Vedenov AM, Dichne AM, Frank-Kameneteskii MD. Helix-Coli transition of DNA. Uspechi physicheskich nauk (In Russian). 1971; 105 :479-519 - 4.
Saenger W. Principles of Nucleic Acids Structure. New York, Berlin, Heidelberg, Tokyo: Springer-Verlag; 1984. 515 p - 5.
Dickerson RE. DNA structure from A to Z. Methods in Enzymology. 1992; 211 :67-111 - 6.
Vardevanyan PO. Structural Transitions in DNA and DNA-Protein Complexes at Different Functional States of Genome [Thesis]. Yerevan: Yerevan State University; 1990 - 7.
Crothers DM, Zimm BH. Theory of the melting transition of synthetic polynucleotides: Evaluation of the stacking free energy. Journal of Molecular Biology. 1964; 9 :1-9 - 8.
Vardevanyan PO, Antonyan AP, Parsadanyan MA, Davtyan HG, Karapetyan AT. The binding of ethidium bromide with DNA: Interaction with single- and double-stranded structures. Experimental and Molecular Medicine. 2003; 35 :527-533 - 9.
Vardevanyan PO, Antonyan AP, Hambardzumyan LA, Shahinyan MA, Karapetian AT. Thermodynamic analysis of DNA complexes with methylene blue, ethidium bromide and Hoechst 33258. Biopolymers and Cell. 2013; 29 :515-520 - 10.
Vardevanyan P, Antonyan A, Parsadanyan M, Shahinyan M, Melqonyan G. Behavior of ethidium bromide-Hoechst 33258-DNA and ethidium bromide-methylene blue-DNA triple systems by means of UV melting. International Journal of Spectroscopy. 2015; 2015 :1-5 - 11.
Chaires JB. Energetics of drug–DNA interactions. Biopolymers. 1997; 44 :201-215 - 12.
Permogorov VI, Frank-Kamenetskii MD, Serdyukova LA, Lazurkin YS. Determination of helix-coil transition heat from the melting curves of deoxyribonucleic acids containing additional bonds between chains. Highmol. Comp. (In Russian). 1965; 7 :362-369 - 13.
Wadkins RM, Jovin TM. Actinomycin D and 7-aminoactinomysin D binding to sigle-stranded DNA. The Biochemist. 1991; 30 :9469-9478 - 14.
Frank-Kamenetskii MD, Karapetyan AT. To the theory of melting od DNA complexes with low-molecular compounds. Journal of Molecular Biology (In Russian). 1972; 6 :621-627 - 15.
Karapetian AT, Mehrabian NM, Terzikian GA, Vardevanian PO, Antonian AP, Borisova OF, Frank-Kamenetskii MD. Theoretical treatment of melting of complexes of DNA with ligands having several types of binding sites on helical and single-stranded DNA. Journal of Biomolecular Structure & Dynamics. 1996; 14 :275-283 - 16.
Vardevanyan PO, Antonyan AP, Parsadanyan MA, Shahinyan MA, Hambardzumyan LA, Torosyan MA, Karapetian AT. The influence of GC/AT composition on intercalating and semi-intercalating binding of Ethidium bromide to DNA. Journal of the Brazilian Chemical Society. 2012; 23 :2016-2020 - 17.
Vardevanyan PO, Antonyan AP, Parsadanyan MA, Pirumyan KV, Muradyan AM, Karapetian AT. Influence of ionic strength on Hoechst 33258 binding with DNA. Journal of Biomolecular Structure and Dynamics. 2008; 25 :641-646 - 18.
Chaires JB. Calorimetry and thermodynamics in drug design. Annual Review of Biophysics. 2008; 37 :138-151 - 19.
Chaires JB. Drug-DNA interactions. Current Opinion in Structural Biology. 1998; 8 :314-320 - 20.
Herry CM. Structure-based drug design. Chemical & Engineering News. 2001; 79 :69-74 - 21.
Suh D, Chaires JB. Criteria for the mode of binding of DNA binding agents. Bioorganic & Medicinal Chemistry. 1995; 3 :723-728 - 22.
Delban F, Quadrifolio F, Giancotti V, Crescenzi V. Comparative microcalorimetric dilatometric analysis of the interactions of quinacrine, chloroquine and ethidium bromide with DNA. Biopolymers. 1982; 21 :331-341 - 23.
Babayan Y, Manzini G, Xodo LE, Quadrifolio F. Base specificity in the interaction of ethidium and synthetic polyribonucleotides. Nucleic Acids Research. 1987; 15 :5803-5812 - 24.
Vardevanyan PO, Antonyan AP, Parsadanyan MA, Davtyan HG, Boyajyan ZR, Karapetyan AT. Complex-formation of ethidium bromide with poly[D(A-T)]⋅poly[D(A-T)]. Journal of Biomolecular Structure and Dynamics 2005; 22 :465‐470 - 25.
Lazurkin YS, Frank-Kamenetskii MD, Trifonov EN. Melting of DNA: Its study and application as a research method. Biopolymers. 1970; 9 :1253-1306 - 26.
Lazurkin YS. Physical Methods of Investigation of Proteins and Nucleic Acids. (in Russian). Moscow: Nauka; 1967 - 27.
Karapetyan AT, Vardevanyan PO, Terzikyan GA, Frank-Kameneteskii MD. Theory of helix-coil transition on DNA-ligand complexes: The effect to two types of interaction of ligand on the parameters of transition. Journal of Biomolecular Structure and Dynamics. 1990; 8 :123-130 - 28.
Karapetyan AT, Permogorov BU, Frank-Kameneteskii MD. In: Andronicashvili EH, editor. Conformational Changes of Biopolymers in Solutions. Moscow: Nauka; 1973. p. 43-49 - 29.
Karapetyan AT, Vardevanyan PO, Frank-Kameneteskii MD. Enthalpy of helix-coil transition of DNA: Dependence on Na ^{+}concentration and GC-content. Journal of Biomolecular Structure and Dynamics. 1990;8 :131-138 - 30.
Klump H, Ackermann T. Influence of the base composition of DNA on the transition enthalpy. Biopolymers. 1971; 10 :513-522 - 31.
Ivanov VI, Minchenkova LE, Minyat EE, Frank-Kamenetskii MD, Schyolkina AK. The B to A transition of DNA in solution. Journal of Molecular Biology. 1974; 87 :817-833 - 32.
Ivanov VI, Krylov DU, Minyat EE, Minchenkova LE. B-A transition in DNA. Journal of Biomolecular Structure & Dynamics. 1983; 1 :453-460 - 33.
Ivanov VI, Karapetian AT, Minyat EE. Structure and expression. In: Sarma RH, Sarma MH, editors. DAN and its Drugs Complexes. Vol. 2. Guilderland, N.Y.: Adenine Press; 1987. p. 205-216 - 34.
Cantor CR, Schimmel PR. Biophysical Chemistry. Part II. The Behavior of Biological Macromolecules. N.Y: W.H. Freeman and company; 1980. 365 p - 35.
Nechipurenko YD. Analysis of Binding of Biologically Active Compounds to Nucleic Acids. (in Russian). Ijevsk: Moscow; 2015. 188 p - 36.
Nechipurenko YD. Analysis of binding of ligands to nucleic acids. Molecular Biophysics 2014; 59 :12‐36 - 37.
Vardevanyan PO, Antonyan AP, Manukyan GA, Karapetian AT, Shchelkina AK. Borisova OF. Molecular Biology (Russia). 2000; 34 :310-315 - 38.
McGhee JD, von Hippel PH. Theoretical aspects of DNA-protein interactions: Co-operative and non-cooperative binding of large ligands to a one-dimensional homogeneous lattice. Journal of Molecular Biology. 1974; 86 :469‐489 - 39.
Karapetian AT, Mehrabian NM, Terzikian GA, Antonian AP, Vardevanian PO, Frank-Kamenetskii MD. Journal of Biomolecular Structure & Dynamics. 1998; 14 :229-265

Submitted: 21 June 2017 Reviewed: 27 September 2017 Published: 18 July 2018

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.