Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in four 40 base pair deoxyoligonucleotides

Effects of different end sequences on melting, circular dichroism spectra (CD), and enzyme binding properties were investigated for four 40 base pair, non-self-complementary duplex DNA oligomers. The center sequences of these oligoduplexes have either of two 22 base pair modules flanked on both sides by sequences differing in AT content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from a van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming that the melting transition is two-state. Melting free energies (20 degrees C) evaluated from DSC melting experiments on the four duplex DNAs ranged from -52.2 to -77.5 kcal/mol. Free energies based on the van't Hoff analysis were -37.9 to -58.8 kcal/mol. Although the values are different, trends in the melting free energies of the four duplex DNAs as a function of sequence were identical in both DSC and optical analyses. Subject to several assumptions, values for the initiation free energy were estimated for each duplex, defined as DeltaG(int) = DeltaG(cal) - DeltaG(pred), where DeltaG(cal) is the experimental free energy at 20 degrees C determined from the experimentially measured values of the transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal). The predicted free energy of the sequence, DeltaG(pred)(20 degrees C), is obtained using published nearest-neighbor sequence stability values. For three of the four duplexes, values of DeltaG(int) are essentially nil. In contrast, the duplex with 81.8% GC has a considerably higher estimate of DeltaG(int) = 7.1 kcal/mol. The CD spectra for the six duplexes collected over the wavelength range from 200 to 320 nm are also sequence-dependent. Factor analysis of the CD spectra by singular value decomposition revealed that the experimental CD spectra could be reconstructed from linear combinations of two minor and one major subspectra. Changes in the coefficients of the major subspectrum for different sequences reflect incremental sequence-dependent variations of the CD spectra. Equilibrium binding by BamHI restriction endonuclease to the 40 base pair DNAs whose central eight base pairs contain the recognition sequence for BamHI restriction enzyme bounded by A.T base pairs, 5'-A-GGATCC-A-3' was investigated. Binding assays were performed by titering BamHI against a constant concentration of each of the duplex DNA substrates, in the absence of Mg(2+), followed by analysis by gel retardation. Under the conditions employed, the enzyme binds but does not cleave the DNAs. Results of the assays revealed two binding modes with retarded gel mobilities. Binding isotherms for the fraction of bound DNA species versus enzyme concentration for each binding mode were constructed and analyzed with a simple two-step equilibrium binding model. This analysis provided semiquantitative estimates on the equilibrium binding constants for BamHI to the four DNAs. Values obtained for the binding constants varied only 7-fold and ranged from 6 x 10(-)(8) to 42 x 10(-)(8) M, with binding free energies from -8.6 to -9.7 (+/- 0.2) kcal/mol depending on the sequence that flanks the enzyme binding site. Unlike what was found earlier in binding studies of the 22 base pair duplexes that constitute the core modules of the present 40-mers [Riccelli, P. V., Vallone, P. M., Kashin, I., Faldasz, B. D., Lane, M. J., and Benight, A. S. (1999) Biochemistry 38, 11197-11208], no obvious relationship between binding and stability was found for these longer DNAs. Apparently, effects of flanking sequence stability on restriction enzyme binding may only be measurable in very short duplex deoxyoligonucl     

Gilvocarcin V exhibits both equilibrium DNA binding and UV light induced DNA adduct formation which is sequence context dependent.

The relative degree of both equilibrium binding and of ultraviolet light induced adduct formation for the antitumor antibiotic gilvocarin V with two hexaecamer DNA sequence isomers, d[ATATATAGCTATATAT]2 and d[AAAAAAAGCTTTTTTT]2, was assessed. The experiments reveal that gilvocarin V binds, under equilibrium conditions, and reacts, in the presence of exogenously applied UV light, more efficiently with the alternating purine:pyrimidine sequence hexadecamer than the homopurine:homopyrimidine duplex at identical gilvocarcin V to DNA duplex ratios. DNAse I digests of adduct containing duplexes derived from the d[AAAAAAAGCTTTTTTT]2 duplex, identified and isolated using gel shift assays employing denaturing polyacrylamide gels, confirm that gilvocarcin V adducts can be formed with thymine residues but suggest that adduct formation with either adenine or guanine residues is also possible.     

A 145-base pair DNA sequence that positions itself precisely and asymmetrically on the nucleosome core.

A 145-bp DNA sequence, cloned from Escherichia coli, was reconstituted into nucleosome core particles by a number of methods. The behaviour of the resulting complex upon sucrose gradient sedimentation and nucleoprotein gel electrophoresis closely resembled that of control bulk nucleosome core particles. DNase I digestion of the 32P-end-labelled complex revealed the 10-bp periodicity of cleavages expected for DNA bound on a histone surface. The narrow cleavage sites observed (1 bp wide) imply that the sequence occupies a single preferred position on the nucleosome core, accurate to the level of single base pairs. By relating the digestion pattern observed to the pattern of site protection found for random sequence nucleosomes, the DNA position was found to be offset by 17 bp from that in the normal core particle. A number of experiments argue against the involvement of length or end effects and suggest that it is some feature of the DNA sequence itself that determines this precise positioning of DNA on the nucleosome.     

The interaction of core histones with DNA: equilibrium binding studies.

The binding of core histone proteins to DNA, measured as a function of [NaCl[ is a reversible process. Dissociation and reassociation occurs in two stages. Between 0.7 and 1.2 M NaCl H2a H2b bind non-cooperatively as an equimolar complex with deltaGo = 1.6 Kcals/mole at 4 degree C and 1.0 M NaCl. Between 1.2 and 2.0 M NaCl H3 and H4 bind cooperatively as an equimolar complex with delta Go = 7.4 Kcal/mole at 4 degree C and 1.0 M NaCl. The proper binding of H2a and H2b requires the presence of bound H3 and H4. Nuclease digestion of the H3-H4 DNA produces a tetramer of H3-H4 bound to fragments of DNA 145, 125 and 104 base pairs long. Thus an H3-H4 tetramer can protect fragments of DNA as long as those found in complete core particles and must therefore span the nucleosome core particle.     

Porphyrin and metalloporphyrin binding to DNA polymers: rate and equilibrium binding studies

Interactions of meso-tetrakis(4-N-methylpyridiniumyl)porphyrin [TMpyP(4)] with poly[d(G-C)].poly[d(G-C)] [poly[d(G-C)2] and poly[d(A-T)].poly[d(A-T)] [poly[d(A-T)2] were studied by equilibrium dialysis and stopped-flow dissociation kinetics as a function of [Na+]. Metalloderivatives of TMpyP(4), NiTMpyP(4), and ZnTMpyP(4) were also investigated. The apparent equilibrium binding constants (Kobs) were approximately the same for TMpyP(4) binding to either poly[d(G-C)2] or poly[d(A-T)2] and decreased with increasing [Na+]. The slopes of the plots of log Kobs vs log [Na+] were similar, with values close to -2.7. Contrary to implications in previously reported studies, these data do not indicate that TMpyP(4) prefers to bind to GC sites at low ionic strength and to AT sites at high ionic strength. In contrast, binding of ZnTMpyP(4) to these two polymers is very different. Comparisons of Kobs values at 0.065 M [Na+] indicate that ZnTMpyP(4) binding to AT sites is approximately 200 times more favorable than binding to GC sites, a finding in agreement with previous qualitative observations. Although the binding of the Zn species to the GC polymer was too weak for us to assess the salt effect, the plot of log Kobs vs log [Na+] gave a slope of -2.0 for ZnTMpyP(4) binding to poly[d(A-T)2]. Application of condensation theory for polyelectrolytes suggests similar charge interactions for ZnTMpyP(4) and for TMpyP(4) binding to poly[d(A-T)2]. Likewise, the rates of dissociation from poly[d(A-T)2] were similar for TMpyP(4) and ZnTMpyP(4) [and also NiTMpyP(4)]. However, whereas TMpyP(4) [and NiTMpyP(4)] dissociation from poly[d(G-C)2] was measurable, that for ZnTMpyP(4) was too fast to measure.(ABSTRACT TRUNCATED AT 250 WORDS)     

Raman spectroscopic studies of the DNA cro binding site conformation, free and bound to cro protein.

Raman spectra of the DNA binding site for cro repressor protein were obtained in the presence and absence of bound cro protein. The 17 base pair fragment is a consensus sequence of the six cro binding sites in phage lambda, except that the second base to the right of the center of pseudosymmetry is altered. Analysis of the spectrum of the free DNA indicates that the molecule exists in a B-like conformation with deviations from the usual B form occurring mainly in the bands assigned to A-T vibrations. The spectrum of the bound DNA was obtained by subtracting the spectrum of free cro from the spectrum of the complex which was estimated to be 90% bound. The DNA undergoes significant structural changes upon binding to the protein; most notable of these changes is a destacking of the G-C bases reflected by increases in the 1240, 1262, and 1320 cm-1 bands. A decrease in the 1361 cm-1 band that occurs has also been assigned to a destacking in guanine bases. The appearance of a 705 cm-1 band and the decrease and downshift of the 670 cm-1 band are consistent with the appearance of A-like character in the A-T region of the binding site when the protein binds; however, the spectra indicate that the entire binding site remains in a distorted B-like conformation. We use the 705 cm-1 band to estimate A-like character because the 800-850 cm-1 region is obscured by interference from strong protein bands. Other shifts in both intensity and position cannot be assigned to characteristic changes in conformation and therefore must be attributed to the protein influencing the structure in a novel way.     

Calorimetric and spectroscopic studies of aminoglycoside binding to AT-rich DNA triple helices

Calorimetric and fluorescence techniques were used to characterize the binding of aminoglycosides-neomycin, paromomycin, and ribostamycin, with 5′-dA₁₂-x-dT₁₂-x-dT₁₂-3′ intramolecular DNA triplex (x = hexaethylene glycol) and poly(dA)·2poly(dT) triplex. Our results demonstrate the following features: (1) UV thermal analysis reveals that the T_m for triplex decreases with increasing pH value in the presence of neomycin, while the T_m for the duplex remains unchanged. (2) The binding affinity of neomycin decreases with increased pH, although there is an increase in observed binding enthalpy. (3) ITC studies conducted in two buffers (sodium cacodylate and MOPS) yield the number of protonated drug amino groups (Δn) as 0.29 and 0.40 for neomycin and paromomycin interaction with 5′-dA₁₂-x-dT₁₂-x-dT₁₂-3′, respectively. (4) The specific heat capacity change (ΔC_p) determined by ITC studies is negative, with more negative values at lower salt concentrations. From 100 mM to 250 mM KCl, the ΔC_p ranges from −402 to −60 cal/(mol K) for neomycin. At pH 5.5, a more positive ΔC_p is observed, with a value of −98 cal/(mol K) at 100 mM KCl. ΔC_p is not significantly affected by ionic strength. (5) Salt dependence studies reveal that there are at least three amino groups of neomycin participating in the electrostatic interactions with the triplex. (6) FID studies using thiazole orange were used to derive the AC₅₀ (aminoglycoside concentration needed to displace 50% of the dye from the triplex) values. Neomycin shows a seven fold higher affinity than paromomycin and eleven fold higher affinity than ribostamycin at pH 6.8. (7) Modeling studies, consistent with UV and ITC results, show the importance of an additional positive charge in triplex recognition by neomycin. The modeling and thermodynamic studies indicate that neomycin binding to the DNA triplex depends upon significant contributions from charge as well as shape complementarity of the drug to the DNA triplex Watson–Hoogsteen groove.     

Influence of supercoiling and sequence context on operator DNA binding with lac repressor

The dissociation of the repressor-operator complex from a series of negatively supercoiled plasmid DNAs was examined as a function of the sequence context, orientation, and spacing. The plasmids were grouped into four classes, each with common sequence context. The highest dissociation rate constants were observed for the plasmids containing only a single operator (or pseudooperator) sequence, while approximately 10-fold lower rate constants were measured for plasmids with the I gene pseudooperator in conjunction with either the Z gene pseudooperator or the primary operator. Comparison of the behavior of these two classes of plasmids demonstrated the importance of two operator sequences and supported a model of DNA loop formation to stabilize the repressor-operator complex (Whitson, P. A., and Matthews, K. S. (1986) Biochemistry 25, 3845-3852; Whitson, P. A., Olson, J. S., and Matthews, K. S. (1986) Biochemistry 25, 3852-3858; Whitson, P. A., Hsieh, W. T., Wells, R. D., and Matthews, K. S. (1987) J. Biol. Chem. 262, 4943-4946; Kr?mer, H., Niem?ller, M., Amouyal, M., Revet, B., von Wilcken-Bergmann, B., and Müller-Hill, B. (1987) EMBO J. 6, 1481-1491). The third class, with intermediate dissociation rate constants, was comprised of plasmids which contained the primary operator and the higher affinity pseudooperator normally located in the Z gene. Neither the additional presence of the I gene pseudooperator nor the orientation of the primary operator relative to the Z gene pseudooperator significantly affected the dissociation rate constants. The binding characteristics of this group of plasmids demonstrated the essential role of the Z gene pseudooperator in the formation of intramolecular ternary complex and suggested an in vivo function for this pseudooperator. Plasmids containing two primary operator sequences were the class with lowest dissociation rate constants from lac repressor, and minimal effects of salt or spacing on dissociation of this class were observed. These data are consistent with formation of an intramolecular complex with a looped DNA segment stabilized by the combination of increased local concentration of binding sites and torsional stresses on the DNA which favor binding in supercoiled DNA.     

Thermodynamic and structural properties of pentamer DNA.DNA, RNA.RNA, and DNA.RNA duplexes of identical sequence

Four pentamers with the general sequence 5'CU(T)GU(T)G/5'CACAG have been prepared by chemical synthesis in order to generate duplex structures with common sequences. The four duplexes studied include the DNA.DNA duplex (5'dCACAG/5'dCTGTG) and the RNA.RNA duplex (5'rCUGUG/5'rCACAG) as well as the two corresponding DNA.RNA heteroduplexes (5'rCUGUG/5'dCACAG and 5'CACAG/5'dCTGTG). The measured entropy, enthalpy, and free energy changes upon melting are reported for each pentamer and compared to the predicted values where possible. Results show that the two DNA.RNA heteroduplexes are destabilized (delta G degrees 25 = -4.2 +/- 0.4 kcal/mol) relative to either the DNA.DNA duplex (delta G degrees 25 = -4.8 +/- 0.5 kcal/mol) or the RNA.RNA duplex (delta G degrees 25 = -5.8 +/- 0.6 kcal/mol). Circular dichroism spectra indicate that the RNA and the two heteroduplexes adopt an A-form conformation, while the DNA conformation is B-form. Imino proton NMR spectra also show that the heteroduplex structures resemble the RNA.RNA duplex.     

DNA binding studies of antibiotic drug cephalexin using spectroscopic and molecular docking techniques

The purpose of this study was to explore an accurate characterization of the binding interaction of antibiotic drug cephalexin with calf thymus DNA (CT-DNA) as a relevant biological target by using UV absorption, fluorescence spectroscopy and circular dichroism (CD) in vitro under simulated physiological conditions (pH = 7.4) and also through a molecular modeling study. The results showed that the drug interacts with the DNA helix via a minor groove binding mode. The thermodynamic parameters were calculated and showed that the reaction between the drug and CT-DNA was exothermic. In addition, the drug enforced traceable changes in the viscosity of DNA. The molecular modeling results indicated that cephalexin forcefully binds to the minor groove of DNA with a relative binding energy of ?21.02?kJ mol^?1. The obtained theoretical results were in good agreement with those obtained from experimental studies.     

Thermodynamic analysis of saccharide binding to snake gourd (Trichosanthes anguina) seed lectin. Fluorescence and absorption spectroscopic studies.

The interaction of different saccharides with the snake gourd (Trichosanthes anguina) seed lectin (SGSL) was investigated by fluorescence spectroscopy. Binding of 4-methylumbelliferyl-beta-D-galactopyranoside (MeUmb beta Gal) to SGSL resulted in a significant increase in the fluorescence emission intensity of the sugar at 376 nm, and this change was used to estimate the association constants for the binding interaction. Interestingly, the increase in emission intensity changed with a change in temperature, increasing from 19.2% at 20 degrees C to 80.2% at 40 degrees C. At 20 degrees C the association constant, K(a), for the MeUmb beta Gal-SGSL interaction was found by fluorescence titration to be 5.8 x 10(4) M(-1). From the temperature dependence of the association constants, the changes in enthalpy (Delta H) and entropy (Delta S) associated with binding of MeUmb beta Gal to SGSL were estimated to be -80.85 kJ.mol(-1) and -184.0 J.mol(-1).K(-1), respectively. Binding of unlabeled sugars was investigated by monitoring the decrease in fluorescence intensity when they were added to a mixture of SGSL and MeUmb beta Gal. The Ka values for different sugars were determined at several temperatures, and Delta H and Delta S were determined from the van't Hoff plots. Enthalpy-entropy compensation was noticed in all cases. The results indicate that saccharide binding to SGSL is enthalpy-driven and the negative contribution from entropy is, in general, quite high.     

Triple helical DNA in a duplex context and base pair opening

It is fundamental to explore in atomic detail the behavior of DNA triple helices as a means to understand the role they might play in vivo and to better engineer their use in genetic technologies, such as antigene therapy. To this aim we have performed atomistic simulations of a purine-rich antiparallel triple helix stretch of 10 base triplets flanked by canonical Watson–Crick double helices. At the same time we have explored the thermodynamic behavior of a flipping Watson–Crick base pair in the context of the triple and double helix. The third strand can be accommodated in a B-like duplex conformation. Upon binding, the double helix changes shape, and becomes more rigid. The triple-helical region increases its major groove width mainly by oversliding in the negative direction. The resulting conformations are somewhere between the A and B conformations with base pairs remaining almost perpendicular to the helical axis. The neighboring duplex regions maintain a B DNA conformation. Base pair opening in the duplex regions is more probable than in the triplex and binding of the Hoogsteen strand does not influence base pair breathing in the neighboring duplex region.