Search

The second derivative techniques, which are used in spectroscopy to reveal poorly resolved bands, were adjusted for processing of calorimetric melting profiles of DNAs from higher organisms. The procedure strongly increased the sensitivity of the profiles to various distortions of the double helix and allowed clear resolution of the constituents arising from the melting of satellite DNA (minor narrow peaks) and the rest of genome (dominant wide peak). As an illustration of the advantages of the procedure, selective suppression of narrow satellite peaks was demonstrated as DNA was chemically modified by the antitumor drug cisplatin. As the cisplatin/nucleotide molar ratio increased from 0.001 to 0.05, the peaks decreased in height until they fully disappeared. The influence of cisplatin on the shape of the main peak, which corresponds to the rest of the genome, was much weaker. Transplatin, an ineffective analog of cisplatin, suppressed all constituents of the melting profile without selectivity.

Many factors that change the temperature position and interval of the DNA helix-coil transition often also alter the shape of multi-peak differential melting curves (DMCs). For DNAs with a multi-peak DMC, there is no agreement on the most useful definition for the melting temperature, Tm, and temperature melting width, ΔT, of the entire DNA transition. Changes in Tm and ΔT can reflect unstable variation of the shape of the DMC as well as alterations in DNA thermal stability and heterogeneity. Here, experiments and computer modeling for DNA multi-peak DMCs varying under different factors allowed testing of several methods of defining Tm and ΔT. Indeed, some of the methods give unreasonable "jagged" Tm and ΔT dependences on varying relative concentration of DNA chemical modifications (rb), [Na(+)], and GC content. At the same time, Tm determined as the helix-coil transition average temperature, and ΔT, which is proportional to the average absolute temperature deviation from this temperature, are suitable to characterize multi-peak DMCs. They give smoothly varying theoretical and experimental dependences of Tm and ΔT on rb, [Na(+)], and GC content. For multi-peak DMCs, Tm value determined in this way is the closest to the thermodynamic melting temperature (the helix-coil transition enthalpy/entropy ratio).

Antitumor activity of cisplatin is exerted by covalent binding to DNA. For comparison, studies of cisplatin-DNA complexes often employ the very similar but inactive transplatin. In this work, thermal and thermodynamic properties of DNA complexes with these compounds were studied using differential scanning calorimetry (DSC) and computer modeling. DSC demonstrates that cisplatin decreases thermal stability (melting temperature, Tm) of long DNA, and transplatin increases it. At the same time, both compounds decrease the enthalpy and entropy of the helix–coil transition, and the impact of transplatin is much higher. From Pt/nucleotide molar ratio rb = 0.001, both compounds destroy the fine structure of DSC profile and increase the temperature melting range (ΔT). For cisplatin and transplatin, the dependences δTm vs rb differ in sign, while δΔT vs rb are positive for both compounds. The change in the parameter δΔT vs rb demonstrates the GC specificity in the location of DNA distortions. Our experimental results and calculations show that 1) in contrast to [Pt(dien)Cl]Cl, monofunctional adducts formed by transplatin decrease the thermal stability of long DNA at [Na+] > 30 mM; 2) interstrand crosslinks of cisplatin and transplatin only slightly increase Tm; 3) the difference in thermal stability of DNA complexes with cisplatin vs DNA complexes with transplatin mainly arises from the different thermodynamic properties of their intrastrand crosslinks. This type of crosslink appears to be responsible for the antitumor activity of cisplatin. At any [Na+] from interval 10–210 mM, cisplatin and transplatin intrastrand crosslinks give rise to destabilization and stabilization, respectively.

Many experiments demonstrate that regions with higher GC-content in natural DNAs unwind at higher temperatures adsorbing more heat than equivalently sized regions with lower GC-content. This simple observation implies that normalized calorimetric melting profiles (calorimetric cDMCs) will not be equivalent differential melting curves (DMCs). We propose simple expressions for long natural and random DNA sequences to reciprocally convert DMCs and corresponding calorimetric cDMCs. The expressions are confirmed by the Poland–Fixman–Freire method and an approach based upon mixtures of homopolymeric duplexes. Using these expressions and experimental calorimetric data, we demonstrate that the average relative deviation between DMC and cDMC is proportional to the temperature melting range of the helix-coil transition Δ T . Corresponding difference between melting temperatures is proportional to Δ T 2 . In general, sequence and ionic conditions influence the deviation through their effect on Δ T . On the basis of the developed approach, we propose a method to determine the thermodynamic melting temperature (ratio of calorimetric enthalpy and entropy of the helix-coil transition) for natural DNAs from optical DMCs without calorimetric experiments.

The Poland-Fixman-Freire formalism was adapted for modeling of calorimetric DNA melting profiles, and applied to plasmid pBR 322 and long random sequences. We studied the influence of the difference (HGC -HAT ) between the helix-coil transition enthalpies of AT and GC base pairs on the calorimetric melting profile and on normalized calorimetric melting profile. A strong alteration of DNA calorimetrical profile with HGC -HAT was demonstrated. In contrast, there is a relatively slight change in the normalized profiles and in corresponding ordinary (optical) normalized differential melting curves (DMCs). For fixed HGC -HAT , the average relative deviation (S) between DMC and normalized calorimetric profile, and the difference between their melting temperatures (Tcal -Tm ) are weakly dependent on peculiarities of the multipeak fine structure of DMCs. At the same time, both the deviation S and difference (Tcal -Tm ) enlarge with the temperature melting range of the helix-coil transition. It is shown that the local deviation between DMC and normalized calorimetric profile increases in regions of narrow peaks distant from the melting temperature.

Although many anticancer drugs exert their biological activity by forming DNA interstrand crosslinks (ICLs), the thermodynamics of biologically relevant long crosslinked DNAs has not been intensively studied in contrast to short duplexes. Here, we carry out computer modeling of the shift of melting temperature of long DNAs caused by ICLs taking into account crosslinking effect in itself and concomitant local alterations in the free energy (δG) of the helix-coil transition at sites of ICLs. Depending on δG, DNA interstrand crosslinks at per nucleotide concentration r = 0.05 can change the melting temperature by value from -17 to +47°C, and the influence weakly depends on DNA sequence and GC content. A change in melting temperature caused by introduction of interstrand crosslinking in modified DNA at sites of modifications also depends on δG and varies from 0 to +12°C. Comparison with experiment for the three platinum crosslinking compounds demonstrates utility of the theoretical method for understanding how crosslinking compounds can influence the melting behavior. On the basis of the method, interdependence of local distortions and crosslinking in itself was studied for thermal effect of ICLs. A method for evaluating the nature of the structural alteration that produces a change in thermal stability for short crosslinked DNA is also proposed. The methods can be used for comparative thermodynamic characterization of various DNA crosslinking agents.

It is generally accepted that DNA conserves its secondary structure after a freeze-thaw cycle. A negligible amount of degradation occurs after this procedure. Degradation becomes appreciable only after multiple cycles of freezing and thawing. In this study, we have found that a single freeze-thaw cycle in alkaline medium (pH>or=10.8) gives rise to denaturation of calf thymus DNA, although the melting temperature of intact DNA in the solution used for the freeze-thaw experiments is higher than 60 degrees C. The degree of denaturation is almost independent of the regime of freezing. The melting curve obtained after DNA is frozen at -2 degrees C and then thawed is almost the same as after a freezing carried out in liquid nitrogen (-196 degrees C). However, incubation in the same solution at 0 degrees C for 24 hours without freezing does not give rise to any denaturation. The degree of denaturation caused by freezing increases with pH (if pH>or=10.8) and decreases with Na2CO3 concentration at fixed pH and [Na+], although Na2CO3 decreases the melting temperature of intact DNA. A preliminary treatment of DNA with cisplatin or transplatin (0.01 Pt atoms per nucleotide) gives rise to a full recovery of the DNA secondary structure after freezing and thawing similar to what occurs after heating DNA to 100 degrees C and quick cooling. Possible mechanisms that may cause DNA denaturation during a freeze-thaw cycle in alkaline medium are discussed.

A computer modeling of thermodynamic properties of a long DNA of N base pairs that includes omega interstrand crosslinks (ICLs), or omega chemical modifications involving one strand (monofunctional adducts, intrastrand crosslinks) has been carried out. It is supposed in our calculation that both types of chemical modifications change the free energy of the helix-coil transition at sites of their location by deltaF. The value deltaF0 corresponds to stabilization, i.e., to the increase in melting temperature. It is also taken into account that ICLs form additional loops in melted regions and prohibit strand dissociation after full DNA melting. It is shown that the main effect of interstrand crosslinks on the stability of long DNA's is caused by the formation of additional loops in melted regions. This formation increases DNA melting temperature (Tm) much stronger than replacing omega base pairs of AT type with GC. A prohibition of strand dissociation after crosslinking, which strongly elevates the melting temperature of oligonucleotide duplexes, does not influence melting behavior of long DNA's (Nor=1000 bp). As was demonstrated earlier for the modifications involving one or the other strand, the dependence of the shift of melting temperature deltaTm on the relative number of modifications r=omega/(2N) is a linear function for any deltaF, and deltaTm(r) identical with 0 for the ideal modifications (deltaF=0). We have shown that deltaTm(r) is the same for periodical and random distribution if the absolute value of deltaF is less 2 kcal. The absolute value of deltaTm(r) at deltaF2 kcal and deltaF-2 kcal is higher for periodical distribution. For interstrand crosslinks, the character of the dependence deltaTm(r) is quite different. It is nonlinear, and the shape of the corresponding curve is strongly dependent on deltaF. For "ideal" interstrand crosslinks (deltaF=0), the function deltaTm(r) is not zero. It is monotone positive nonlinear, and its slope decreases with r. If r0.004, then the entropy stabilizing effect of interstrand crosslinking itself exceeds the influence of a distortion of the double helix at sites of their location. The resulting deltaTm(r) is positive even in the case of the infinite destabilization at sites of the ICLs (deltaF---infinity). In general, stabilizing influence of interstrand crosslinks is almost fully compensated for by local structural distortions caused by them if 0r0.01 and -infinitydeltaFapproximately -5 kcal per mole of ICLs. The absolute value of deltaTm(r) in this case is approximately 10 times lower than for ideal ICLs as well as for local destabilization without crosslinking. Interstrand crosslinks formed by cisplatin do not change the melting temperature of long DNA's because of this compensation effect.

Unselective and reversible adsorption of ligands on DNA for a model of binding proposed by Zasedatelev, Gursky, and Volkenshtein is considered. In this model, the interaction between neighboring ligands located at the distance of i binding centers is characterized by the statistical weight ai. Each ligand covers L binding centers. For this model, expressions for binding averages are represented in a new simple form. This representation is convenient for the calculation of the fraction of inter-ligand distances of i binding centers fd(i) and the fraction of binding centers included in the distances of i binding centers fbc(i) for various types of interaction between bound ligands. It is shown that, for non-cooperative binding, contact cooperativity and long-range cooperativity, the fraction of the zero inter-ligand distance fd(0) is maximal at any relative concentration of bound ligands (r). Calculations demonstrate that, at low r, fd(0) approximately r.ao, and fd(i) approximately r at 1or=i0.3.r(-1) for non-cooperative binding (ai=1 for all i) and contact inter-ligand interaction (ao not equal to 1, ai=1 for i not equal 0). The function fd(i) decreases in this i interval, but the decrease is very slow. If i1/r-L, then fd(i) rapidly decreases with i at any r for all types of inter-ligand interaction. At high ligand concentration (r is close to rmax=L(-1)), fd(0) is close to unity and fd(i) rapidly decreases with i for any type of inter-ligand interaction. For strong contact cooperativity, fd(0) is close to unity in a much lager r interval ((0.5-1).rmax), and fd(1) approximately ao(-1) at r approximately 0.5.rmax. In the case of long-range interaction between bound ligands, the dependence fd(i) is more complex and has a maximum at i approximately (1/r-L)1/2 for anti-cooperative binding. fbc(i) is maximal at i approximately 1/r-L for all types of binding except the contact cooperativity. A strong asymmetry in the influence of contact cooperativity and anticooperativity on the ligand distribution along DNA is demonstrated.

We describe a new approach to calculate the binding of flexible peptides and unfolded proteins to multicomponent lipid membranes. The method is based on the transfer matrix formalism of statistical mechanics recently described as a systematic tool to study DNA-protein-drug binding in gene regulation. Using the energies of interaction of the individual polymer segments with different membrane lipid species and the scaling corrections due to polymer looping, we calculate polymer adsorption characteristics and the degree of sequestration of specific membrane lipids. The method is applied to the effector domain of the MARCKS (myristoylated alanine rich C kinase substrate) protein known to be involved in signal transduction through membrane binding. The calculated binding constants of the MARCKS(151-175) peptide and a series of related peptides to mixed PC/PS/PIP2 membranes are in satisfactory agreement with in vitro experiments. Copyright  2008 European Peptide Society and John Wiley &S ons, Ltd.

Addition of Na(2)CO(3) to almost salt-free DNA solution (5.10(-5)M EDTA, pH=5.7, T(m)=26.5 degrees C) elevates both pH and the DNA melting temperature (T(m)) if Na(2)CO(3) concentration is less than 0.004 M. For 0.004 M Na(2)CO(3), T(m)=58 degrees C is maximal and pH=10.56. Further increase in concentration gives rise to a monotonous decrease in T(m) to 37 degrees C for 1M Na(2)CO(3) (pH=10.57). Increase in pH is also not monotonous. The highest pH=10.87 is reached at 0.04 M Na(2)CO(3) (T(m)=48.3 degrees C). To reveal the cause of this DNA destabilization, which happens in a narrow pH interval (10.56/10.87) and a wide Na(2)CO(3) concentration interval (0.004/1M), a procedure has been developed for determining the separate influences on T(m) of Na(+), pH, and anions formed by Na(2)CO(3) (HCO(3)(-) and CO(3)(2-)). Comparison of influence of anions formed by Na(2)CO(3) on DNA stability with Cl(-) (anion inert to DNA stability), ClO(4)(-) (strong DNA destabilizing "chaotropic" anion) and OH(-) has been carried out. It has been shown that only Na(+) and pH influence T(m) in Na(2)CO(3) solution at concentrations lower than 0.001 M. However, the T(m) decrease with concentration for [Na(2)CO(3)]>/=0.004 M is only partly caused by high pH=10.7. Na(2)CO(3) anions also exert a strong destabilizing influence at these concentrations. For 0.1M Na(2)CO(3) (pH=10.84, [Na(+)]=0.2M, T(m)=42.7 degrees C), the anion destabilizing effect is higher 20 degrees C. For NaClO(4) (ClO(4)(-) is a strong "chaotropic" anion), an equal anion effect occurs at much higher concentrations approximately 3M. This means that Na(2)CO(3) gives rise to a much stronger anion effect than other salts. The effect is pH dependent. It decreases fivefold at neutral pH after addition of HCl to 0.1M Na(2)CO(3) as well as after addition of NaOH for pH greater than 11.2.

DNA chemical modifications caused by the binding of some antitumor drugs give rise to a very strong local stabilization of the double helix. These sites melt at a temperature that is well above the melting temperatures of ordinary AT and GC base pairs. In this work we have examined the melting behavior of DNA containing very stable sites. Analytical expressions were derived and used to evaluate the thermodynamic properties of homopolymer DNA with several different distributions of stable sites. The results were extended to DNA with a heterogeneous sequence of AT and GC base pairs. The results were compared to the melting properties of DNA with ordinary covalent interstrand cross-links. It was found that, as with an ordinary interstrand cross-link, a single strongly stabilized site makes a DNA's melting temperature (T(m)) independent of strand concentration. However in contrast to a DNA with an interstrand cross-link, a strongly stabilized site makes the DNA's T(m) independent of DNA length and equal to T(infinity), the melting temperature of an infinite length DNA with the same GC-content and without a stabilized site. Moreover, at a temperature where more than 80% of base pairs are melted, the number of ordinary (non-modified) helical base pairs (n) is independent of both the DNA length and the location of the stabilized sites. For this condition, n(T) = (2 omega-a)S/(1-S) and S = exp[DeltaS(T(infinity)-T)/(RT)] where omega is the number of strongly stabilized sites in the DNA chain, a is the number of DNA ends that contain a stabilized site, and DeltaS, T, and R are the base pair entropy change, the temperature, and the universal gas constant per mole. The above expression is valid for a temperature interval that corresponds to n0.2N for omega=1, and n0.1N for omega1, where N is the number of ordinary base pairs in the DNA chain.

Short melted regions less than 100 base pairs (bp) in length are rarely found in the differential melting curves (DMC) of natural DNAs. Therefore, it is supposed that their characteristics do not affect DNA melting behavior. However, in our previous study, a strong influence of the form of the entropy factor of small loops on melting of cross-linked DNAs was established (D. Y. Lando, A. S. Fridman et al., Journal of Biomolecular Structure and Dynamics, 1997, Vol. 15, pp. 141–150; Journal of Biomolecular Structure and Dynamics, 1998, Vol. 16, pp. 59–67). Quite different dependencies of the melting temperature on the relative concentration of interstrand cross-links were obtained for the loop entropy factors given by the Fixman–Freire (Jacobson–Stockmayer) and Wartell–Benight relations. In the present study, the influence of the entropy factor of small loops on the melting of natural DNAs, cross-linked DNAs and periodical double-stranded polynucleotides is compared using computer simulation. A fast combined computational method for calculating DNA melting curves was developed for this investigation. It allows us to assign an arbitrary dependence of the loop entropy factor on the length of melted regions for the terms corresponding to small loops (less than τ bp in length). These terms are calculated using Poland's approach. The Fixman–Freire approach is used for long loops. Our calculations have shown that the temperature dependence of the average length of interior melted regions (loops) has a maximum at T ≈ Tm (Tm is the DNA melting temperature) in contrast to the dependence of the total average length of melted regions, which increases almost monotonously. Computer modeling demonstrates that prohibition of formation of loops less than τ base pairs in length does not markedly change the DMC for τ < 150 bp. However, the same prohibition strongly affects the average length of internal melted regions for much smaller τ's. The effect is already noticeable for τ = 1 bp and increases with τ. A tenfold increase in the entropy factor of all loops with length less than τ bp causes a noticeable alteration of the DMC for τ ≥ 30 bp. It is shown that DMCs are identical for the Wartell–Benight and for the Fixman–Freire (Jacobson–Stockmayer) form of the loop entropy factor. However, for low degree of denaturation, the average length of internal melted regions is 40% lower for the Wartell–Benight form due to the fluctuational opening of short AT-rich regions less than 10 bp in length. The same calculations carried out for periodical polynucleotides demonstrate a much stronger difference in melting behavior for different forms of entropy factors of short loops. The strongest difference occurs if the length of stable GC-rich and unstable AT-rich stretches is equal to 30 bp. However, the comparison carried out in this work demonstrates that the entropy factor of short loops influences melting behavior of cross-linked DNA much stronger than of unmodified DNA with random or periodical sequences. © 2001 John Wiley & Sons, Inc. Biopolymers 58: 374–389, 2001

Influence of long-range interactions between ligands bound to DNA molecule on the character of their adsorption is studied using computer modeling. For this investigation, two calculation procedures are developed. They are based upon the method of the free energy minimum and on the partition function method. The both procedures demonstrate that in the case of a strong enough attraction between all the bound ligands their binding to DNA has the character of phase transition of the first kind. There is a break in the binding curve c(c0) where c - relative concentration of bound ligands, c0 - molar concentration of free ligands. The break occurs because there is an interval of central degrees of binding (approximately 50% of the maximum c value) that are prohibited for individual DNA molecules. Such a transition might be caused by some types of DNA condensation. Attraction between the neighboring ligands only, adjacent or/and separated by double helix regions, does not cause this effect.

The influence of cis-diamminedichloroplatinum(II) (cis-DDP) and trans-diamminedichloroplatinum(II) (trans-DDP) binding on the helix-coil transition of DNAs with different GC composition is investigated experimentally and theoretically (26% ≤ GC ≤ 72%; 0.003 ⩽ rb ⩽ 3; where rb is the number of moles of bound DDP per mole of nucleotide). In contrast to previous studies, it is demonstrated that cis-DDP increases the DNA melting temperature (Tm) if GC content and rb are low enough. Stabilization occurs for rb = 0.003 if GC = 26% or 42%. Destabilization takes place for higher GC and rb values. Trans-DDP causes only stabilization independently of GC content if rb ⩽ 0.3. For the highest GC content (72%), very low concentration rb = 0.003) of both compounds gives rise to a strong shift of the melting temperature which is equal to −4.5°C for cis-DDP and +3.5°C for trans-DDP. It is shown that well known bidentate and monodentate, kinetically inert platinum binding cannot bring about such a strong Tm alteration. Computer simulation of the helix-coil transition of DNA-ligand complexes is used to explain the observed effects.

Covalent and strong coordination binding to DNA of a large number of antitumour drugs and other compounds leads to interstrand cross-link formation. To investigate cross-link influence on double helix stability, two methods are developed for the calculation of melting curves. The first method is based on Poland's approach. It requires computer time proportional to u.N, where u is the average distance (in base pairs) between neighboring cross-links and N is the number of base pairs in the DNA chain. The method is more suitable when u is not large, and small loops formed by interstrand cross-links in melted regions strongly affect DNA melting. The computer time for the second method, based on the Fixman-Freire approach, does not depend on the number of cross-links and is proportional to I.N (I is the number of exponential functions used for a decomposition of the loop entropy factor). It is more appropriate when N and u are large, and therefore particular values of the entropy factors of small loops do not influence DNA melting behavior.

Penetrating into cell nuclei, antitumor drug cisplatin sequentially forms various intermediate and final adducts destroying local DNA structure. The demonstrated disappearance of the fine structure of melting curve of long DNAs along with a strong decrease in melting enthalpy conforms to the structural impact. However, the negative thermal effect (δ T m ) caused by cisplatin is relatively small if neutral medium is used in melting experiments. Cisplatin's inactive analogs transplatin and diethylenetriaminechloroplatinum {Pt[(dien)Cl]Cl} also distort DNA structure but their thermal effect is even positive. We have found that the use of alkaline medium in melting experiments strengthens the negative thermal effect for cisplatin. For transplatin and Pt[(dien)Cl]Cl, the thermal effect becomes negative that makes it qualitatively consistent with structural distortions. Those changes are explained by elimination of nonspecific electrostatic stabilization of DNA under platination. Additionally, alkaline medium fixes intermediate states of DNA platination and makes them stable against heating. These results allowed us to monitor δ T m under binding of platinum compounds to DNA and their further transformation. The kinetic and thermal characteristics of monofunctional and bifunctional adducts were evaluated. It has been demonstrated that monofunctional adducts of cisplatin, transplatin and Pt[(dien)Cl]Cl produce approximately the same thermal destabilization. Cisplatin intrastrand crosslinks cause a two-fold stronger thermal destabilization than its monofunctional adducts. The value of δ T m for cisplatin's final adducts is ten times larger than for transplatin. This difference mainly comes from the much stronger thermal destabilizing power of cisplatin's intrastrand crosslinks, which are responsible for antitumor activity of this compound.

Binding of the antitumor compound cisplatin to DNA locally distorts the double helix. These distortions correlate with a decrease in DNA melting temperature (Tm). However, the influence of cisplatin on DNA stability is more complex because it decreases the DNA charge density. In this way, cisplatin increases the melting temperature and partially compensates for the destabilizing influence of structural distortions. The stabilization is stronger at low Na+ ion concentration. Due to this compensation, the total decrease in the DNA melting temperature after cisplatin binding is much lower than the decrease caused by the distortions themselves, especially at low [Na+]. It is shown in this study that, besides Na+ concentration, pH also strongly influences the value of a change in the melting temperature caused by cisplatin. In alkaline medium (pH=10.5-10.8), a fall in the melting temperature caused by platination is enhanced several times with respect to neutral medium. Such a stronger drop in Tm is explained by a decrease in pK values of base pairs caused by lowering the charge density under platination that facilitates proton release. At neutral pH, the proton release is low for both control and platinated DNA and does not influence the melting behavior. Therefore, lowering in the charge density under platination, besides stabilization, gives additional destabilization just in alkaline medium. Destabilization caused by structural distortions due to this pH induced compensation of stabilizing effect is more pronounced. In the presence of carbonate ion, destabilization caused by high pH value is strengthened. As a decrease in DNA charge density, interstrand crosslinking caused by cisplatin also increases the DNA stability due to loss in the entropy of the melted state. However, computer modeling of DNA stability demonstrates that interstrand crosslinks formed by cisplatin do not stabilize long DNA. It is shown that the increase in Tm caused by interstrand crosslinking itself is compensated for by a local destabilization of the double helix at the sites of location of interstrand crosslinks formed by cisplatin.

DNA interstrand cross-links are usually formed due to bidentate covalent or coordination binding of a cross-linking agent to nucleotides of different strands. However interstrand linkages can be also caused by any type of chemical modification that gives rise to a strong local stabilization of the double helix. These stabilized sites conserve their helical structure and prevent local and total strand separation at temperatures above the melting of ordinary AT and GC base pairs. This local stabilization makes DNA melting fully reversible and independent of strand concentration like ordinary covalent interstrand cross-links. The stabilization can be caused by all the types of chemical modifications (interstrand cross-links, intrastrand cross-links or monofunctional adducts) if they give rise to a strong enough local stabilization of the double helix. Our calculation demonstrates that an increase in stability by 25 to 30 kcal in the free energy of a single base pair of the double helix is sufficient for this "cross-linking effect" (i.e. conserving the helicity of this base pair and preventing strand separation after melting of ordinary base pairs). For the situation where there is more then one stabilized site in a DNA duplex (e.g., 1 stabilized site per 1000 bp), a lower stabilization per site is sufficient for the "cross-linking effect" (18 - 20 kcal). A substantial increase in DNA stability was found in various experimental studies for some metal-based anti-tumor compounds. These compounds may give rise to the effect described above. If ligand induced stabilization is distributed among several neighboring base pairs, a much lower minimum increase per stabilized base pair is sufficient to produce the cross-linking effect (1 bp- 24.4 kcal; 5 bp- 5.3 kcal; 10 bp- 2.9 kcal, 25 bp- 1.4 kcal; 50 bp- 1.0 kcal). The relatively weak non-covalent binding of histones or protamines that cover long regions of DNA (20- 40 bp) can also cause this effect if the salt concentration of the solution is sufficiently low to cause strong local stabilization of the double helix. Stretches of GC pairs more than 25 bp in length inserted into poly(AT) DNA also exhibit properties of stabilizing interstrand cross-links.

A theoretical method for computer modeling of DNA condensation caused by ligand binding is developed. In the method, starting (s) and condensed (c) states are characterized by different free energies for ligand free DNA (F(s) and F(c) respectively), ligand binding constants (K(s) and K(c)) and stoichiometry dependent parameters (c(sm) and c(cm) - maximum relative concentration of bound ligands (per base pair) for starting and condensed state respectively). The method allows computation of the dependence of the degree of condensation (the fraction of condensed DNA molecules) on ligand concentration. Calculations demonstrate that condensation transition occurs under an increase in ligand concentration if F(s)F(c) (i.e. S(sc) = exp [- (F(c) - F(s)) / (RT)], the equilibrium constant of the s-c transition, is low (S(sc)1)) and K(s)K(c). It was also found that condensation is followed by decondensation at high ligand concentration if the condensed DNA state provides the number of sites for ligand binding less than the starting state (c(sm)c(cm)). A similar condensation-decondensation effect was found in recent experimental studies. We propose its simple explanation.

Long-range interaction between all the ligands bound to DNA molecule may give rise to adsorption with the character of phase transition of the first kind (D. Y. Lando, V. B. Teif, J. Biomol. StructDynam. 18, 903-911 (2000)). In this case, the binding curve, c(c(o)), is characterized by a sudden change of the relative concentration of bound ligands ((c)) at a critical concentration of free (unbound) ligands, c(o)=c(ocr), from a low c value to a high one where c(o) is molar concentration of free ligands. Such a transition might be caused by some types of DNA condensation or changes in DNA topology. For the study of the conditions necessary for adsorption with the character of phase transition, a calculation procedure based on the method of the free energy minimum is developed. The ligand size and two types of interactions between ligands adsorbed on DNA molecule are taken into consideration: long-range interaction between all the ligands bound to DNA and contact interactions between neighboring ligands. It was found that a) Stronger long-range interaction is required for longer ligands to induce phase transition that is occurred at greater c(ocr) values; b) Pure contact interaction between neighboring ligands can not itself initiate phase transition. However contact cooperativity strongly decreases the threshold value of energy of long-range interaction necessary to give rise to the transition.

Interactions of meso-tetra-(4-N-oxyethylpyridyl) porphyrin (TOEPyP(4)), its 3-N analog (TOEPyP(3)) and their Co, Cu, Ni, Zn metallocomplexes with duplex DNA have been investigated by uv/visible absorbance and circular dichrosim spectroscopies. Results reveal the interactions of these complexes with duplex DNA are of two types. (1) External binding of duplex DNA by metalloporphyrins containing Zn and Co, and (2) Binding of duplex DNA both externally and internally (by intercalation) by porphyrins not containing metals, and metalloporphyrins containing Cu and Ni. Results indicate that (4N-oxyethylpyridyl) porphyrins intercalate more preferably in the structure of duplex DNA and have weaker external binding than 3N-porphyrins.

DNA is a promising construction material for the engineering of artificial nanostractured devices [1]. One of possible DNA implications in bionanodevices is detecting of metal ions. This possibility is based on the fact that metal ions may preferentially bind to definite DNA conformation and thus metal ion binding may give rise to transition between A-, B-, or Z-DNA [2]. Another approach based on utilizing the specific DNA sequence required to detect specific metals was reported recently [3]. In this method single-stranded DNA forms “pocket” that accepts only lead ions. Here we consider a potential DNA-based sensor detecting metal ions which is based on the phenomenon of DNA condensation.

A theoretical method is developed for calculation of melting curves of covalent complexes of DNA with antitumor drugs. The method takes into account all the types of chemical modifications of the double helix caused by platinum compounds and DNA alkylating agents: 1) monofunctional adducts bound to one nucleotide; 2) intrastrand cross-links which appear due to bidentate binding of a drug molecule to two nucleotides that are included into the same DNA strand; 3) interstrand cross-links caused by bidentate binding of a molecule to two nucleotides of different strands. The developed calculation method takes into account the following double helix alterations at sites of chemical modifications: 1) a change in stability of chemically modified base pairs and neighboring ones, that is caused by all the types of chemical modifications; 2) a change in the energy of boundaries between helical and melted regions at sites of chemical modification (local alteration of the factor of cooperativity of DNA melting), that is caused by all the types of chemical modifications, too; 3) a change in the loop entropy factor of melted regions that include interstrand cross-links; 4) the prohibition of divergence of DNA strands in completely melted DNA molecules, which is caused by interstrand cross-links only. General equations are derived, and three calculation methods are proposed to calculate DNA melting curves and the parameters that characterize the helix-coil transition.

In our previous papers I and II (D. Y. Lando et al, J. Biomol. Struct. Dynam. (1997) v. 15, N1, p. 129-140, p. 141-150), two methods were developed for calculation of melting curves of cross-linked DNA. One of them is based on Poland's and another on the Fixman-Freire approach. In the present communication, III, a new theoretical method is developed for computation of differential melting curves of DNAs cross-linked by anticancer drugs and their inactive analogs. As Poland's approach, the method allows study of the influence of the loop entropy factor, delta(n), on melting behavior (n is the length of a loop in base pairs). However the method is much faster and requires computer time that inherent for the most rapid Fixman-Freire calculation approach. In contrast to the computation procedures described before in communications I and II, the method is suitable for computation of differential melting curves in the case of long DNA chains, arbitrary loop entropy factors of melted regions and arbitrary degree of cross-linking including very low values that occur in vivo after administration of antitumor drugs. The method is also appropriate for DNAs without cross-links. The results of calculation demonstrate that even very low degree of cross-linking alters the DNA differential melting curve. Cross-linking also markedly strengthens the influence of particular function delta(n) upon melting behavior.

In the previous paper (D.Y. Lando, J. Biomol. Struct. Dynam, 15, 129-140 (1997)) the melting of cross-linked DNA with N base pairs and omega interstrand cross-links has been considered theoretically. In the present study on the basis of these results, two simple schemes are developed for the computation of melting curves of cross-linked DNA. The investigation of influence of interstrand linking on DNA stability has been carried out by computer simulation. It is shown that the relative concentration of cross-links, CCT = omega/N, their distribution along a DNA molecule, and particular values of the entropy factors of small loops formed by cross-links in melted regions strongly affect the DNA melting temperature, Tm. On the contrary, for DNA without cross-links, a ten-fold increase or decrease in the entropy factors of small loops does not cause the Tm variation. The comparison of the results of calculation with experimental data suggests that the majority of types of cross-link neither maintain ordered parallel orientation of bases in melted regions nor increase considerably the thermostability of cross-linked base pairs. Four different ways of influence of interstrand cross-linking on the DNA double helix stability are considered. It is shown that cross-linking significantly enhances the influence of single strand stiffness in melted regions on DNA melting behavior.

Strong local stabilization of the double helix can be an alternative to the interstrand crosslinks formed in DNA by some antitumor drugs. Therefore, we have carried out a computer modeling of the thermodynamic properties of DNA that contains strongly stabilized sites (SSSs). Melting of DNA locally fastened by SSSs was compared with DNA that includes interstrand crosslinks. Using experimental data from the literature and the results of our calculations, we have shown that SSSs really exist: some irreversibly bound protein molecules and chemical modifications caused by some ruthenium and antitumor platinum compounds form such sites in DNA. The theoretical calculated results for the increase of melting temperatures for random distribution of SSSs are consistent with experimental data.