Noncovalent DNA binding of bis(1,10-phenanthroline)copper(I) and related compounds

The noncovalent DNA binding of the bis(1,10-phenanthroline)copper(I) complex [(Phen)2CuI] was examined under anaerobic conditions by absorption and circular dichroism spectroscopy, and viscometry, as a function of phenanthroline concentration. Analyses according to the McGhee-von Hippel method indicated that binding exhibited both neighbor-exclusion and positive cooperativity effects, with a neighbor-exclusion parameter n approximately 2 and a cooperativity parameter omega approximately 4. The association constant for (Phen)2CuI binding decreased with increasing concentration of phenanthroline in excess over that required to stoichiometrically generate (Phen)2CuI, indicating that free phenanthroline was a weak competitive inhibitor of (Phen)2CuI binding. The maximal association constant for DNA binding of (Phen)2CuI in 0.2 M NaCl and 9.8% ethanol, extrapolated to zero concentration of excess phenanthroline, was 4.7 x 10(4) M-1 (DNA base pairs). The magnitude of the neighbor-exclusion parameter, the changes in spectral properties of (Phen)2CuI induced by DNA binding, and the increase in DNA solution viscosity upon (Phen)2CuI addition are consistent with a model for DNA binding by (Phen)2CuI involving partial intercalation of one phenanthroline ring of the complex between DNA base pairs in the minor groove as suggested previously [Veal & Rill (1989) Biochemistry 28, 3243-3250]. Viscosity measurements indicated that the mono(phenanthroline)copper(I) complex also binds to DNA by intercalation; however, no spectroscopic or viscometric evidence was found for DNA binding of free phenanthroline or the bis(2,9-dimethyl-1,10-phenanthroline)copper(I) complex. DNA binding of free phenanthroline may be cooperative and induced by prior binding of (Phen)2CuI.     

DNA binding of tilorone: 1H NMR and calorimetric studies of the intercalation

The fluorene derivative tilorone has received great attention as a DNA intercalator and has been widely recognized as an inducer of interferon. The biological activity of tilorone is known to be related to its binding mode with DNA; however, few structural and thermodynamic studies have elaborated on this issue. This paper presents two-dimensional (2-D) NMR and isothermal titration calorimetry (ITC) for the tilorone/DNA complex, coupled with circular dichroism (CD) spectroscopy and viscosity measurements. NMR investigation suggests that tilorone binds to DNA through intercalation, showing greater affinity for insertion between AT base pairs than between CG pairs. CD spectral changes were observed for T/B (tilorone/DNA base pair molar ratio) ratios greater than the stoichiometric ratio generally expected for intercalators (i.e., T/B = 0.5, according to the neighbor-exclusion principle). However, there was a clear plateau in the CD intensity between T/B < 0.35 and T/B > 0.45. From comparison with NMR and other measurements, we postulate that CD changes below the plateau should be related to the intercalation and the latter to electrostatic interactions and nonspecific bindings. ITC data showed that DeltaH < -TDeltaS < 0, which indicated that tilorone/DNA binding is enthalpy controlled. The magnitude of Kb (the binding constant) was of the same order as that of ethidium bromide. The stoichiometric number, obtained from ITC, CD, and UV data, implied a relatively smaller value (0.28-0.35) than that of the neighbor-exclusion principle. This is because side chains located in the groove disrupt further intercalation to the adjacent sites.     

Molecular dynamics of structural transitions and intercalation in DNA

The large-scale flexibility of DNA and the intercalation of actinomycin D have been studied by computer simulation using molecular dynamics. The stretching and unwinding of B and Z forms of DNA and intercalation in B-DNA were examined through molecular dynamics simulations, and the energetics of transitions were calculated by the conformational energy minimization method. The principal results of this research are as follows: (1) A dynamic conformational pathway is presented for longitudinal stretching and unwinding of the double helix to open an intercalation site. (2) Large-scale transitions are possible in both B and Z forms of DNA through a conformationally allowed kinetic pathway. (3) The stretching and untwisting of a 5′(CG)3′ step is energetically more favorable than for a GC step in B-DNA. (4) The formation of an adjacent second cavity in B-DNA requires larger energy than the formation of the first cavity, affirming the neighbor-exclusion principle of intercalation. (5) Docking an intercalated actinomycin D in the stretched structure is shown to be geometrically and energetically feasible.     

Anticooperative Binding Governs the Mechanics of Ethidium-Complexed DNA

DNA intercalators bind nucleic acids by stacking between adjacent basepairs. This causes a considerable elongation of the DNA backbone as well as untwisting of the double helix. In the past few years, single-molecule mechanical experiments have become a common tool to characterize these deformations and to quantify important parameters of the intercalation process. Parameter extraction typically relies on the neighbor-exclusion model, in which a bound intercalator prevents intercalation into adjacent sites. Here, we challenge the neighbor-exclusion model by carefully quantifying and modeling the force-extension and twisting behavior of single ethidium-complexed DNA molecules. We show that only an anticooperative ethidium binding that allows for a disfavored but nonetheless possible intercalation into nearest-neighbor sites can consistently describe the mechanical behavior of intercalator-bound DNA. At high ethidium concentrations and elevated mechanical stress, this causes an almost complete occupation of nearest-neighbor sites and almost a doubling of the DNA contour length. We furthermore show that intercalation into nearest-neighbor sites needs to be considered when estimating intercalator parameters from zero-stress elongation and twisting data. We think that the proposed anticooperative binding mechanism may also be applicable to other intercalating molecules.     

Interaction of a macrocyclic bisacridine with DNA

The binding of the macrocycle SDM to DNA was investigated by visible spectroscopy, stopped-flow kinetics, and NMR spectroscopy. SDM is composed of two 9-aminoacridines linked via the amino groups by a spermine side chain and via the 4-positions by a N,N'-[(methylthio)ethyl]succinamide side chain [Zimmerman, S. C., Lamberson, C. R., Cory, M., & Fairley, T. A. (1989) J. Am. Chem. Soc. 111, 6805-6809]. The visible spectrum of SDM bound to poly[d(A-T)]2 or poly[d(G-C)]2 is red-shifted relative to the spectrum of SDM alone and displays considerable hypochromicity. Results from titrations of SDM with polymer indicate a binding site size of three base pairs per macrocycle. The dissociation constant for SDM bound to either poly[d(A-T)]2 or poly[d(G-C)]2 is an order of magnitude lower than that for a similar bisacridine linked only by a spermine side chain. In addition, the dependence of the dissociation constant on ionic strength is significantly reduced. NMR studies of SDM complexes with poly[d(A-T)]2 or a tetramer, d(CGCG)2, show that intercalation is the mode of binding. The magnitudes of the chemical shift differences for SDM aromatic protons in the free and bound states support intercalation with the acridine ring systems essentially parallel to the long axis of the base pairs. Cross peaks from NOESY spectra of the SDM complex with d(CGCG)2 further support this mode of binding and provide information on the structure of the complex. The results are analyzed for consistency with each of three binding models: (i) bisintercalation with the two side chains in the same groove; (ii) bisintercalation according to the neighbor-exclusion principle with the two side chains in opposite grooves; and (iii) bisintercalation with two side chains in opposite grooves but with violation of the neighbor-exclusion principle. Model i is found to be unlikely on the basis of all evidence obtained, including preliminary modeling studies. Both models ii and iii can be reconciled with the experimental evidence and from a modeling standpoint are energetically feasible.     

A microscopic model of complex formation of intercalators with DNA: analysis of NMR spectroscopy data

A model of intercalation complex formation of ligands with an arbitrary oligonucleotide sequence has been proposed, which takes in the explicit form the microscopic constants and the parameters of the cooperativity of binding to particular sites of the oligomer into account. The model has been adapted for the analysis of NMR spectroscopy data.     

DNA bis-intercalation: Theoretical calculation of binding curves

A statistical mechanical calculation of the binding properties of DNA bis-intercalators is presented, based on the sequence-generating function method of Lifson. The effects of binding by intercalation of one or both chromophores of a bifunctional intercalating agent are examined. The secular equation for a general model that includes the effects of neighbor (nearest and non-nearest) exclusion and/or cooperativity in the binding of both singly and doubly intercalated ligands is derived. Numerical results for binding curves are presented for a more restricted model in which each type of bound ligand rigorously excludes its nearest neighbor and the total number of sites covered by a doubly intercalated ligand is variable. At low values of free ligand concentration bis-intercalation dominates the binding process, while at high value of free ligand concentration, intercalation of only one chromophore per ligand becomes significant due to the unavailability of contiguous free sites required for bis-intercalation. Also, depending on the binding parameters, the free energy of the system can be lowered by a loss of doubly intercalated ligands in favor of singly intercalated ones. Corresponding to this transition in binding mode, the average number of sites occupied by a bound ligand decreases from that characteristic of bis-intercalation to that characteristic of mono-intercalation as free ligand concentration increases. An analysis of Scatchard plots describing bis-intercalation is presented.     

Cooperative DNA binding by lambda integration protein--a key component of specificity

Quantitative analysis of nitrocellulose filter binding data by the method of Clore, Gronenborn and Davies [(1982) J. Mol. Biol. 155, 447-466] has been used to show that lambda integration protein (Int) exhibits cooperativity in binding to specific recognition sites within the attachment site region (lambda attP) of bacteriophage lambda DNA. Optimal values of the equilibrium constant obtained were 3.0(+/- 1.0) X 10(10) M-1 for the P' site using a model of three sites with equal affinity and 1.9(+/- 0.4) X 10(10) M-1 for the P1 site on a two-site model. The value of the cooperativity parameter alpha is 172(+106)(-66) in all cases. The occurrence of a consensus recognition sequence is necessary but not sufficient for strong binding; cooperative interaction between Int molecules binding to adjacent members of an array of binding sites is also essential. The occurrence of binding site arrays distinguishes lambda attP very clearly from other DNA sequences containing single recognition sites by chance.     

Imino 1H- and 31P-NMR analysis of the interaction of propidium and ethidium with DNA

At low temperature and low salt concentration, both imino proton and ³¹p-nmr spectra of DNA complexes with the intercalators ethidium and propidium are in the slow-exchange region. Increasing temperature and/or increasing salt concentration results in an increase in the site exchange rate. Ring-current effects from the intercalated phenanthridinium ring of ethidium and propidium cause upfield shifts of the imino protons of A · T and G · C base pairs, which are quite similar for the two intercalators. The limiting induced chemical shifts for propidium and ethidium at saturation of DNA binding sites are approximately 0.9 ppm for A · T and 1.1 ppm for G · C base pairs. The similarity of the shifts for ethidium and propidium, in both the slow- and fast-exchange regions over the entire titration of DNA, shows that a binding model for propidium with neighbor-exclusion binding and negative ligand cooperativity is correct. The fact that a unique chemical shift is obtained for imino protons at intercalated sites over the entire titration and that no unshifted imino proton peaks remain at saturation binding of ethidium and propidium supports a neighbor-exclusion binding model with intercalators bound at alternating sites rather than in clusters on the double helix. Addition of ethidium and propidium to DNA results in downfield shifts in ³¹P-nmr spectra. At saturation ratios of intercalator to DNA base pairs in the titration, a downfield shoulder (approximately ?2.7 ppm) is apparent, which accounts for approximately 15% of the spectral area. The main peak is at ?3.9 to ?4.0 ppm relative to ?4.35 in uncomplexed DNA. The simplest neighbor-binding model predicts a downfield peak with approximately 50% of the spectral area and an upfield peak, near the chemical shift for uncomplexed DNA, with 50% of the area. This is definitely not the case with these intercalators. The observed chemical shifts and areas for the DNA complexes can be explained by models, for example, that involve spreading the intercalation-induced unwinding of the double helix over several base pairs and/or a DNA sequence- and conformation-dependent heterogeneity in intercalation-induced chemical shifts and resulting exchange rates.     

DNA sequence dependent binding modes of 4',6-diamidino-2-phenylindole (DAPI)

The interactions of DAPI with natural DNA and synthetic polymers have been investigated by hydrodynamic, DNase I footprinting, spectroscopic, binding, and kinetic methods. Footprinting results at low ratios (compound to base pair) are similar for DAPI and distamycin. At high ratios, however, GC regions are blocked from enzyme cleavage by DAPI but not by distamycin. Both poly[d(G-C)]2 and poly[d(A-T)]2 induce hypochromism and shifts of the DAPI absorption band to longer wavelengths, but the effects are larger with the GC polymer. NMR shifts of DAPI protons in the presence of excess AT and GC polymers are significantly different, upfield for GC and mixed small shifts for AT. The dissociation rate constants and effects of salt concentration on the rate constants are also quite different for the AT and the GC polymer complexes. The DAPI dissociation rate constant is larger with the GC polymer but is less sensitive to changes in salt concentration than with the AT complex. Binding of DAPI to the GC polymer and to poly[d(A-C)].poly[d(G-T)] exhibits slight negative cooperativity, characteristic of a neighbor-exclusion binding mode. DAPI binding to the AT polymer is unusually strong and exhibits significant positive cooperativity. DAPI has very different effects on the bleomycin-catalyzed cleavage of the AT and GC polymers, a strong inhibition with the AT polymer but enhanced cleavage with the GC polymer. All of these results are consistent with two totally different DNA binding modes for DAPI in regions containing consecutive AT base pairs versus regions containing GC or mixed GC and AT base pair sequences. The binding mode at AT sites has characteristics which are similar to those of the distamycin-AT complex, and all results are consistent with a cooperative, very strong minor groove binding mode. In GC and mixed-sequence regions the results are very similar to those observed with classical intercalators such as ethidium and indicate that DAPI intercalates in DNA sequences which do not contain at least three consecutive AT base pairs.     

Analysis of cooperativity and ion effects in the interaction of quinacrine with DNA

The interaction of quinacrine with calf thymus DNA was monitored at several different ionic strengths using spectrophotometric and equilibrium dialysis techniques. The binding results can be explained, assuming each base pair is a potential binding site, using a model containing two negative cooperative effects: (1) ligand exclusion at binding sites adjacent to a filled binding site and (2) ligand–ligand negative cooperativity at adjacent filled binding sites. The logarithm of the observed equilibrium constant (K_obs) determined by this model varies linearily with log[Na⁺], as predicted by the ion condensation theory for polyelectrolytes. When the log K_obs plot is correlated for sodium release by DNA in the intercalation conformational change, the predicted number of ion pairs between the ligand and DNA is approximately two, as expected for the quinacrine dication. Even though K_obs depends strongly on ionic strength, the ligand negative cooperativity parameter ω was found to be indpendent of ionic strength within experimental error. This finding is also in agreement with the ion condensation theory, which predicts a relatively constant amount of condensed counterion on the DNA double helix over this ionic strength range. Drugs would, therefore, experience a relatively constant ionic environment when complexed to DNA even though the ionic conditions of the solvent could change considerably.     

Electrostatic factors in DNA intercalation

The factors that determine the binding of a chromophore between the base pairs in DNA intercalation complexes are dissected. The electrostatic potential in the intercalation plane is calculated using an accurate ab initio based distributed multipole electrostatic model for a range of intercalation sites, involving different sequences of base pairs and relative twist angles. There will be a significant electrostatic contribution to the binding energy for chromophores with a predominantly positive electrostatic potential, but this varies significantly with sequence, and somewhat with twist angle. The usefulness of these potential maps for understanding the binding of intercalators is explored by calculating the electrostatic binding energy for 9-aminoacridine, ethidium, and daunomycin in a variety of model binding sites. The electrostatic forces play a major role in the positioning of an intercalating 9-aminoacridine and a significant stabilizing role in the binding of ethidium in its sterically constrained position, but the intercalation of daunomycin is determined by the side-chain binding. Sequence preferences are likely to be determined by a complex and subtle mixture of effects, with electrostatics being just one component. The electrostatic binding energy is also unlikely to be a major determinant of the twist angle, as its variation with angle is modest for most intercalation sites. Overall, the electrostatic potential maps give guidance on how positively charged chromophores can be chemically adapted by heteroatomic substitution to optimise their binding.     

Porphyrin intercalation and non-specific 'edge on' outside binding to natural DNA

A theoretical two-mode binding model for porphyrin binding to natural DNA is presented. One of the binding modes is assumed to be base sequence specific with binding sites n base-pairs long. The other binding mode has binding sites which consist of only one base-pair and can involve cooperativity. The model fits satisfactorily to data for H2TMPyP-4, Cu(II)TMPyP-3 and Cu(II)TMPyP-4 binding to calf thymus DNA in both a high (mu congruent to 1.0 M) and a low (mu congruent to 0.2 M) ionic strength buffer. The results show that the fraction of porphyrin bound in the non-specific mode reaches a maximum at certain input DNA to porphyrin concentrations ratios. The value of this maximum decreased, and its position shifted to higher DNA to porphyrin concentration ratios for binding in the high ionic strength buffer. The value of the cooperativity parameter obtained through the fitting process suggests that the non-specific binding is positively cooperative. The results are compared with the data analysed using other techniques.     

The antitumor agent mitoxantrone binds cooperatively to DNA: evidence for heterogeneity in DNA conformation

The equilibrium binding of the antitumor compound DHAQ, or mitoxantrone [1,4-dihydroxy-5,8-bis[[2-[(2-hydroxyethyl)amino]ethyl]amino]-9,10- anthracenedione], to various DNAs has been examined by optical titration and equilibrium dialysis methods. At low r (bound drug/DNA base pair) values, r less than 0.03, DHAQ binds, in a highly cooperative manner, to calf thymus and Micrococcus lysodeikticus DNAs. The binding isotherms for the interaction of DHAQ with Clostridium perfringens DNA and poly(dA-dT).poly(dA-dT) exhibit a small positive slope at low r values, suggestive of cooperative binding. In contrast, the binding of DHAQ to poly(dG-dC).poly(dG-dC) shows no evidence of cooperative binding even at very low r values. At higher r values (r greater than 0.05), the binding of DHAQ to all the DNAs studied is characterized by a neighbor-exclusion process. A model is proposed to account for the two modes of binding exhibited in the cooperative binding isotherms. The main feature of the proposed model is that local sequence and structural heterogeneity of the DNA give rise to sets of binding sites to which DHAQ binds in a highly cooperative manner, while the majority of the DNA sites bind DHAQ via a neighbor-exclusion process. This two-site model reproduces the observed binding isotherms and leads to the conclusion that DHAQ binds in clusters to selected regions of DNA. It is suggested that clustering may play a role in the physiological activity of drugs.     

The two modes of binding of Ru(phen)(2)dppz(2+) to DNA: thermodynamic evidence and kinetic studies

The binding of Ru(phen)(2)dppz(2+) (dppz=dipyrido[3,2-a:2',3'-c]phenazine) to DNA was investigated at pH 7.0 and 25 degrees C using stopped-flow and spectrophotometric methods. Equilibrium measurements show that two modes of binding, whose characteristics depend on the polymer to dye ratio (C(P)/C(D)), are operative. The binding mode occurring for values of C(P)/C(D) higher than 3 exhibits positive cooperativity, which is confirmed by kinetic experiments. The reaction parameters are K=2 x 10(3)M(-1), omega=550, n=1, k(r)=(1.9+/-0.5) x 10(7)M(-1)s(-1) and k(d)=(9.5+/-2.5)x10(3)s(-1) at I=0.012 M. The results are discussed in terms of prevailing surface interaction with DNA grooves accompanied by partial intercalation of the dppz residue. The other binding mode becomes operative for C(P)/C(D)<3 and the equilibria analysis shows this is an ordinary intercalation mode (K=1.3 x 10(6) M(-1), n=1.5 at I=0.012 M and K=2 x 10(5) M(-1), n=1.2 at I=0.21 M). Similar behaviour is displayed by double-stranded poly(A).     

Analysis of ligand binding process using binding capacity concept

Binding capacity is the homotropic second derivative of the binding potential with respect to the chemical potential of the ligand. It provides a measure of steepness of the binding isotherm and represents the extent of cooperativity. In the present study, the shape of the binding capacity curve for various systems was investigated and the relation between binding capacity and the extent of cooperativity examined. In this regard, a novel linear graphical method was introduced for binding data analysis. The stoichiometry of binding and the extent of cooperativity can be determined by this method. This method has been successfully applied to various systems such as binding of oxygen to hemoglobin, warfarin to human serum albumin and dodecyltrimethylammonium bromide to alpha-amylase.     

Methylene Blue Binding to DNA with Alternating AT Base Sequence: Minor Groove Binding is Favored over Intercalation

Abstract

The results presented in this paper on methylene blue (MB) binding to DNA with AT alternating base sequence complement the data obtained in two former modeling studies of MB binding to GC alternating DNA. In the light of the large amount of experimental data for both systems, this theoretical study is focused on a detailed energetic analysis and comparison in order to understand their different behavior. Since experimental high-resolution structures of the complexes are not available, the analysis is based on energy minimized structural models of the complexes in different binding modes. For both sequences, four different intercalation structures and two models for MB binding in the minor and major groove have been proposed. Solvent electrostatic effects were included in the energetic analysis by using electrostatic continuum theory, and the dependence of MB binding on salt concentration was investigated by solving the non-linear Poisson-Boltzmann equation. We find that the relative stability of the different complexes is similar for the two sequences, in agreement with the interpretation of spectroscopic data. Subtle differences, however, are seen in energy decompositions and can be attributed to the change from symmetric 5′-YpR-3′ intercalation to minor groove binding with increasing salt concentration, which is experimentally observed for the AT sequence at lower salt concentration than for the GC sequence. According to our results, this difference is due to the significantly lower non-electrostatic energy for the minor groove complex with AT alternating DNA, whereas the slightly lower binding energy to this sequence is caused by a higher deformation energy of DNA. The energetic data are in agreement with the conclusions derived from different spectroscopic studies and can also be structurally interpreted on the basis of the modeled complexes. The simple static modeling technique and the neglect of entropy terms and of non-electrostatic solute-solvent interactions, which are assumed to be nearly constant for the compared complexes of MB with DNA, seem to be justified by the results.     

Methylene blue binding to DNA with alternating AT base sequence: minor groove binding is favored over intercalation

The results presented in this paper on methylene blue (MB) binding to DNA with AT alternating base sequence complement the data obtained in two former modeling studies of MB binding to GC alternating DNA. In the light of the large amount of experimental data for both systems, this theoretical study is focused on a detailed energetic analysis and comparison in order to understand their different behavior. Since experimental high-resolution structures of the complexes are not available, the analysis is based on energy minimized structural models of the complexes in different binding modes. For both sequences, four different intercalation structures and two models for MB binding in the minor and major groove have been proposed. Solvent electrostatic effects were included in the energetic analysis by using electrostatic continuum theory, and the dependence of MB binding on salt concentration was investigated by solving the non-linear Poisson-Boltzmann equation. We find that the relative stability of the different complexes is similar for the two sequences, in agreement with the interpretation of spectroscopic data. Subtle differences, however, are seen in energy decompositions and can be attributed to the change from symmetric 5'-YpR-3' intercalation to minor groove binding with increasing salt concentration, which is experimentally observed for the AT sequence at lower salt concentration than for the GC sequence. According to our results, this difference is due to the significantly lower non-electrostatic energy for the minor groove complex with AT alternating DNA, whereas the slightly lower binding energy to this sequence is caused by a higher deformation energy of DNA. The energetic data are in agreement with the conclusions derived from different spectroscopic studies and can also be structurally interpreted on the basis of the modeled complexes. The simple static modeling technique and the neglect of entropy terms and of non-electrostatic solute-solvent interactions, which are assumed to be nearly constant for the compared complexes of MB with DNA, seem to be justified by the results.     

Binding isotherms and the interaction between proflavine and a DNA of high G · C content

Binding isotherms corresponding to several situations of ligand binding to a linear polymer are calculated, including cases of cooperativity or anticooperativity between the bound ligand states, multiple binding modes that are competitive or non competitive, and possible exclusion of an arbitrary number of adjacent sites upon occupancy of a site by a single ligand. The sequence generating function method of Lifson and Bradley is used, requiring the assumption that no end effects are involved. The case of strong binding of the dye proflavine to a DNA of high G. C content, that of M. lysodeikticus, is considered in detail, and a single model capable of reconciling the available kinetic and equilibrium data on this system, involving two competing binding modes, is discussed.