Nonspecific interactions in dye binding to DNA. Influence of alcohols and amides

The binding of a few drugs (ethidium bromide, propidium diiodide, proflavine and actinomycin D) to DNA has been investigated in aqueous solutions to which cosolvents of different polarity have been added. It is found that both alcohols (less polar than water) and amides (more polar) lower the binding constant according to a linear relationship between the intercalation free energy and cosolvent concentration. The main action of cosolvents cannot be described in terms of electrostatic effects, since they predict much smaller changes in the binding constant than those observed. It appears instead that relevant solvation effects are responsible for the binding strength of the different dyes to DNA. As a general result, it is found that solvation effects largely contribute to the intercalation free energy, thereby weakening the influence of nonspecific interactions at the intercalation site.     

Optimization of electrostatic interactions in protein-protein complexes

In this article, we present a statistical analysis of the electrostatic properties of 298 protein-protein complexes and 356 domain-domain structures extracted from the previously developed database of protein complexes (ProtCom, http://www.ces.clemson.edu/compbio/protcom). For each structure in the dataset we calculated the total electrostatic energy of the binding and its two components, Coulombic and reaction field energy. It was found that in a vast majority of the cases (>90%), the total electrostatic component of the binding energy was unfavorable. At the same time, the Coulombic component of the binding energy was found to favor the complex formation while the reaction field component of the binding energy opposed the binding. It was also demonstrated that the components in a wild-type (WT) structure are optimized/anti-optimized with respect to the corresponding distributions, arising from random shuffling of the charged side chains. The degree of this optimization was assessed through the Z-score of WT energy in respect to the random distribution. It was found that the Z-scores of Coulombic interactions peak at a considerably negative value for all 654 cases considered while the Z-score of the reaction field energy varied among different types of complexes. All these findings indicate that the Coulombic interactions within WT protein-protein complexes are optimized to favor the complex formation while the total electrostatic energy predominantly opposes the binding. This observation was used to discriminate WT structures among sets of structural decoys and showed that the electrostatic component of the binding energy is not a good discriminator of the WT; while, Coulombic or reaction field energies perform better depending upon the decoy set used.     

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.     

N,N-diethylaminoalkoxy-4,7-phenanthrolines and 1,8-diazaflourene derivatives, a novel class of potential interferon inducers and antiviral agents: interactions with nucleic acids in vitro and cellular activities

Tilorone aza-analogues, derivatives of 4,7-phenanthroline and 1,8-diazafluorene, were examined as DNA-complexing agents by spectral and electrophoretic methods. The binding process includes at least two types of interactions: electrostatic and, possibly, intercalation. Complex formation with the denatured DNA was also observed, but its nature remained unsolved. Binding and thermodynamic parameters were determined. All ligands studied showed weak antiviral activity and essentially no interferon induction when assayed in vitro and in vivo. It was concluded that interferon induction by tilorone may involve specific cell receptors or intermediaries.     

Anatomy of specific interactions between lambda repressor and operator DNA

Recognition of specific DNA sequences by proteins is essential for regulation of gene expression. To fully understand the recognition mechanism, it is necessary to understand not only the structure of the specific protein-DNA interactions but also the energetics. We therefore performed a computer analysis in which a phage DNA-binding protein, lambda repressor, was used to examine the changes in binding free energy (DeltaDeltaG) and its energy components caused by single base mutations. We then determined which of the calculated energy components best correlated with the experimental data. The experimental DeltaDeltaG values were well reproduced by the calculations. Component analysis revealed that the electrostatic and hydrogen bond energies were most strongly correlated with the experimental data. Among the 51 single base-substitution mutants examined, positive DeltaDeltaG values, corresponding to weakened binding, were caused by the loss of favorable electrostatic interactions and hydrogen bonds, the introduction of steric collisions and electrostatic repulsion, the loss of favorable interactions with a thymine methyl group, and the increase of unfavorable hydration energy from isolated DNA. This analysis also showed distinct patterns of recognition at A-T and G-C positions, as different combinations of energy components were involved in DeltaDeltaG caused by the two substitution types. We have thus been able to identify the energy components that most strongly correlate with sequence-dependent DeltaDeltaG and determine their contribution to the specificity of DNA sequence recognition by the lambda repressor. Application of this method to other systems should provide additional insight into the molecular mechanism of protein-DNA recognition.     

Anthracycline antibiotic arugomycin binds in both grooves of the DNA helix simultaneously: an NMR and molecular modelling study.

Perturbations to the 1H and 31P chemical shifts of DNA resonances together with twenty-four intermolecular nuclear Overhauser effects show that the anthracycline antibiotic arugomycin intercalates between the basepairs of the hexamer duplex d(5'-GCATGC)2 at the 5'-CpA and 5'-TpG binding sites. In the complex two drug molecules are bound per duplex with full retention of the dyad symmetry. Arugomycin adopts a threaded binding orientation with chains of sugars positioned in both the major and minor groove of the helix simultaneously. The complex is stabilized by hydrogen bonding, electrostatic and van der Waals interactions principally in the major groove and involving substituents on the rigidly oriented bicycloamino-glucose sugar of the antibiotic. A specific hydrogen bond is identified between the C2'-hydroxyl and the guanine N7 at the intercalation site. Together, interactions in the major groove appear to account for the intercalation specificity of arugomycin that requires both a guanine and thymine at the intercalation site. We are unable to identify any sequence specific interactions between the minor groove and the arugarose sugar (S1) which binds only weakly, through van der Walls contacts, over the d(GCA).d(TGC) trinucleotide sequence. The data indicate that the sugar chains of arugomycin are flexible and play little part in the interaction of the antibiotic with DNA. The intensity of sequential internucleotide NOEs identifies the intercalation site as being assymmetric. A family of conformers computed using restrained energy minimisation and molecular dynamics indicate that basepair buckling is a feature of the anthracycline intercalation site that may serve to maximise intermolecular van der Waals interactions by wrapping the basepairs around the antibiotic chromophore.     

Energetics of complex formation of the DNA hairpin structure d(GCGAAGC) with aromatic ligands

The energy contributions of various physical interactions to the total Gibbs energy of complex formation of the biologically important DNA hairpin d(GCGAAGC) with aromatic antitumor antibiotics daunomycin and novantron and the mutagens ethidium and proflavine have been calculated. It has been shown that the relatively small value of the total energy of binding of the ligands to the hairpin is the sum of components great in absolute value and different in sign. The contributions of van der Waals interactions and both intra- and intermolecular hydrogen bonds and bonds with aqueous environment have been studied. According to the calculations, the hydrophobic and van der Waals components are energetically favorable in complex formation of the ligands with the DNA pairpin d(GCGAAGC), whereas the electrostatic (with consideration of hydrogen bonds) and entropic components are unfavorable.     

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.     

Electrostatic contributions to the binding free energy of the lambdacI repressor to DNA.

A model based on the nonlinear Poisson-Boltzmann (NLPB) equation is used to study the electrostatic contribution to the binding free energy of the lambdacI repressor to its operator DNA. In particular, we use the Poisson-Boltzmann model to calculate the pKa shift of individual ionizable amino acids upon binding. We find that three residues on each monomer, Glu34, Glu83, and the amino terminus, have significant changes in their pKa and titrate between pH 4 and 9. This information is then used to calculate the pH dependence of the binding free energy. We find that the calculated pH dependence of binding accurately reproduces the available experimental data over a range of physiological pH values. The NLPB equation is then used to develop an overall picture of the electrostatics of the lambdacI repressor-operator interaction. We find that long-range Coulombic forces associated with the highly charged nucleic acid provide a strong driving force for the interaction of the protein with the DNA. These favorable electrostatic interactions are opposed, however, by unfavorable changes in the solvation of both the protein and the DNA upon binding. Specifically, the formation of a protein-DNA complex removes both charged and polar groups at the binding interface from solvent while it displaces salt from around the nucleic acid. As a result, the electrostatic contribution to the lambdacI repressor-operator interaction opposes binding by approximately 73 kcal/mol at physiological salt concentrations and neutral pH. A variety of entropic terms also oppose binding. The major force driving the binding process appears to be release of interfacial water from the protein and DNA surfaces upon complexation and, possibly, enhanced packing interactions between the protein and DNA in the interface. When the various nonelectrostatic terms are described with simple models that have been applied previously to other binding processes, a general picture of protein/DNA association emerges in which binding is driven by the nonpolar interactions, whereas specificity results from electrostatic interactions that weaken binding but are necessary components of any protein/DNA complex.     

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.     

Thermodynamic analysis of complex formation of ethidium bromide with DNA in water solutions

A thermodynamic analysis of two types of binding of ethidium bromide with DNA: intercalation and binding to the outer surface of a biopolymer has been performed by spectrophotometry. It has been shown that the dominant contribution to the energy of external binding of the ligand to DNA is made by hydrophobic interactions, which lead to less negative values of enthalpy and entropy and more severe negative changes in the heat capacity of complex formation as compared with the intercalation type of binding.     

Structural Analysis of HMGD-DNA Complexes Reveals Influence of Intercalation on Sequence Selectivity and DNA Bending

The ubiquitous, eukaryotic, high-mobility group box (HMGB) chromosomal proteins promote many chromatin-mediated cellular activities through their non-sequence-specific binding and bending of DNA. Minor-groove DNA binding by the HMG box results in substantial DNA bending toward the major groove owing to electrostatic interactions, shape complementarity, and DNA intercalation that occurs at two sites. Here, the structures of the complexes formed with DNA by a partially DNA intercalation-deficient mutant of Drosophila melanogaster HMGD have been determined by X-ray crystallography at a resolution of 2.85 Å. The six proteins and 50 bp of DNA in the crystal structure revealed a variety of bound conformations. All of the proteins bound in the minor groove, bridging DNA molecules, presumably because these DNA regions are easily deformed. The loss of the primary site of DNA intercalation decreased overall DNA bending and shape complementarity. However, DNA bending at the secondary site of intercalation was retained and most protein-DNA contacts were preserved. The mode of binding resembles the HMGB1 box A-cisplatin-DNA complex, which also lacks a primary intercalating residue. This study provides new insights into the binding mechanisms used by HMG boxes to recognize varied DNA structures and sequences as well as modulate DNA structure and DNA bending.     

Coordinate ion pair formation between EcoRI endonuclease and DNA

The free energy of the binding reaction between EcoRI restriction endonuclease and a specific cognate dodecadeoxynucleotide (d(CGCGAATTCGCG)) has contributions from both electrostatic and nonelectrostatic components. These contributions were dissected by measuring the effects of varying salt concentration on the equilibrium binding constant and applying the thermodynamic analyses of Record et al. (Record, M. T., Jr., Lohman, T. M., and deHaseth, P. L. (1976) J. Mol. Biol. 107, 145-158). Endonuclease mutation S187 (Arg 187 to Ser) (Greene, P. J., Gupta, M., Boyer, H. W., Brown, W. E., and Rosenberg, J. M. (1981) J. Biol. Chem. 256, 2143-2153) did not significantly affect the nonelectrostatic component but did perturb the electrostatic contribution to the binding energy (we are numbering the amino acid residues according to the DNA sequence). The former was determined by extrapolating the linear portion of the salt dependence curve (0.125 to 0.25 M KCl) to 1 M ionic strength, with the same result for both wild type and S187 endonucleases at both pH 6.0 and 7.4 (-8.5 +/- 1.5 kcal/mol or greater than 50% of the total binding free energy). The slopes of these same curves yield estimates of eight ionic interactions between wild type endonuclease and the DNA at both pH values. By contrast, binding of EcoRI-S187 to dodecanucleotide involves six charge-charge interactions at pH 6.0. Only two ionic interactions are observed at pH 7.4. This was unexpected since gel permeation chromatography demonstrated that the recognition complex for both wild type and S187 proteins contains an enzyme dimer and a DNA duplex. EcoRI-S187 endonuclease retains wild type DNA sequence specificity, and the rate of the phosphodiester hydrolysis step is also unchanged. Thus, electrostatic interactions are functionally separable from sequence recognition and strand cleavage. Our results also establish that arginine 187 plays a key role in the electrostatic function and suggest that it might be located at the DNA-protein interface. The disproportionate loss of ion pairs at pH 7.4 can be rationalized by a model which suggests that six conformationally mobile ionic groups on the protein act in a coordinated manner during the interaction with DNA.     

Energetics of complex formation of the dna hairpin structure d(GCGAAGC) with aromatic ligands

The energy contributions of various physical interactions to the total Gibbs energy of complex formation of the biologically important DNA hairpin d(GCGAAGC) with aromatic antitumor antibiotics daunomycin and novantrone and mutagens ethidium bromide and proflavin have been calculated. It has been shown that the relatively small value of the total energy of binding of the ligands to the hairpin is the sum of components great in absolute value and different in sign. The contributions of van der Waals interactions and both intra- and intermolecular hydrogen bonds and bonds with aqueous environment have been studied. According to the calculations, the hydrophobic and van der Waals components are energetically favorable in complex formation of the ligands with the hairpin d(GCGAAGC), whereas the electrostatic (with consideration of hydrogen bonds) and entropic components are unfavorable.     

Contribution of enthalpy to the energetics of complex formation of aromatic ligands with DNA

The energy contributions of electrostatic, van der Waals interactions, hydrogen bonds, and interactions of charge transfer type to the enthalpy of complex formation of the double-stand DNA with the antitumor antibiotics daunomycin, nogalamycin, and novantron, as well as the mutagens ethidium bromide and proflavine have been calculated. According to the calculations, the van der Waals component (except for nogalamycin) is energetically favorable during complex formation of the antibiotics with DNA, and the contributions of H bonds and electrostatic interactions are unfavorable, with the probability of charge transfer in the complexes being low. It has been shown that the relatively low value of the experimental enthalpy of binding is the sum of components greater in absolute value and different in the sign, which is the cause of large errors in estimating the total enthalpy of complex formation of aromatic ligands with DNA.     

Studies on the interaction of extended terpyridyl and triazine metal complexes with DNA

The interaction with duplex DNA of a series of 12 previously synthesised ruthenium (II) complexes incorporating extended terpyridines or analogous ligands is reported. The interaction was investigated using absorption spectroscopy and also, when possible, luminescence spectroscopy. Binding parameters for the interactions were estimated using established models. It was found that binding affinities are relatively large and their magnitude is somewhat dependent on the metal centre attached to the ligand. Viscosity studies reveal that all the complexes interact with DNA through groove binding and not intercalation. The variation in binding affinities of these complexes is rationalised by a consideration of electrostatic and steric factors.     

Energetics of sequence-specific protein-DNA association: computational analysis of integrase Tn916 binding to its target DNA

The N-terminal domain of the bacterial integrase Tn916 specifically recognizes the 11 bp DNA target site by positioning the face of a three-stranded beta-sheet into the major groove. Binding is linked to structural adaptation. We have characterized INT-DBD binding to DNA in detail by calorimetry [Milev, S., Gorfe, A., Karshikoff, A., Clubb, R. T., Bosshard, H. R., and Jelesarov, I. (2003) Biochemistry 42, 3481-3491]. Our thermodynamic analysis has indicated that the major driving force of association is the hydrophobic effect while polar interactions contribute less. To gain more comprehensive information about the binding process, we performed a computational analysis of the binding free energy and report here the results. A hybrid molecular mechanics/continuum approach was followed. The total binding free energy is predicted with reasonable accuracy. The calculations confirm that nonpolar effects stabilize the protein-DNA complex while electrostatics opposes binding. Structural changes optimizing surface complementarity are costly in terms of energy. The energetic consequences from the replacement of nine DNA-contacting residues by alanine were investigated. The calculations correctly predict the binding affinity decrease of eight mutations and the destabilizing effect of one wild-type residue. Bulky side chains stabilize the wild-type complex through packing interactions and favorable nonpolar dehydration, but the net nonpolar energy changes do not correlate with the relative affinity loss upon mutation. Discrete protein-DNA electrostatic interactions may be net stabilizing or net destabilizing depending on the local environment. In contrast to nonpolar energy changes, the magnitude of the electrostatic free energy ranks the mutations according to the experimentally measured DeltaDeltaG. Free energy decomposition analysis from a structural perspective leads to detailed information about the thermodynamic strategy used by INT-DBD for sequence-specific DNA binding.