An Iterative Approach to Placing Counterions Around DNA

Abstract

An iterative approach, in which the effect of placing counter ions around DNA influences the electrostatic potential that the other subsequently approaching ions feel, has been used to place sodium ions around polynucleotides. The main focus of this report is to study the sequence and structure dependence on the distribution of ions around DNA, particularly that of tightly bound ions. The interesting results of the calculations are that there is significant sequence dependence on the electrostatic potentials in the B form of DNA, whereas relatively less difference in A form. In the case of Z form, the cations bridge the inter-stand phosphates along the minor groove.     

A negative electrostatic determinant mediates the association between the Escherichia coli trp repressor and its operator DNA.

The electrostatic potential surfaces were characterized for trp repressor models that bind to DNA with sequence specificity, without specificity, and not at all. Comparisons among the surfaces were used to isolate protein surface features likely to be important in DNA binding. Models that differ in protein conformation and tryptophan-analogue binding consistently showed positive potential associated with the protein surfaces that interact with the DNA major groove. However, negative potential is associated with the trp repressor surface that contacts the DNA minor groove. This negative potential is significantly neutralized in the protein conformation that is bound to DNA. Positive potential is also associated with the tryptophan binding-site surface, a consequence of the tryptophan- or tryptophan analogue-induced allosteric change. This protein region is complementary to the strongest negative potential associated with the DNA phosphate backbone and is also present in the isolated protein structure from the protein-DNA complex. The effects of charge-change mutation, pH dependence, and salt dependence on the electrostatic potential surfaces were also examined with regard to their effects on protein-DNA binding constants. A consistent model is formed that defines a role for long-range electrostatics early in the protein-DNA association process and complements previous structural, molecular association, and mutagenesis studies.     

Cationic DMPC/DMTAP lipid bilayers: molecular dynamics study

Cationic lipid membranes are known to form compact complexes with DNA and to be effective as gene delivery agents both in vitro and in vivo. Here we employ molecular dynamics simulations for a detailed atomistic study of lipid bilayers consisting of a mixture of cationic dimyristoyltrimethylammonium propane (DMTAP) and zwitterionic dimyristoylphosphatidylcholine (DMPC). Our main objective is to examine how the composition of the DMPC/DMTAP bilayers affects their structural and electrostatic properties in the liquid-crystalline phase. By varying the mole fraction of DMTAP, we have found that the area per lipid has a pronounced nonmonotonic dependence on the DMTAP concentration, with a minimum around the point of equimolar DMPC/DMTAP mixture. We show that this behavior has an electrostatic origin and is driven by the interplay between positively charged TAP headgroups and the zwitterionic phosphatidylcholine (PC) heads. This interplay leads to considerable reorientation of PC headgroups for an increasing DMTAP concentration, and gives rise to major changes in the electrostatic properties of the lipid bilayer, including a significant increase of total dipole potential across the bilayer and prominent changes in the ordering of water in the vicinity of the membrane. Moreover, chloride counterions are bound mostly to PC nitrogens implying stronger screening of PC heads by Cl ions compared to TAP headgroups. The implications of these findings are briefly discussed.     

Electrostatic Potential of B-DNA: Effect of Interionic Correlations

Modified Poisson-Boltzmann (MPB) equations have been numerically solved to study ionic distributions and mean electrostatic potentials around a macromolecule of arbitrarily complex shape and charge distribution. Results for DNA are compared with those obtained by classical Poisson-Boltzmann (PB) calculations. The comparisons were made for 1:1 and 2:1 electrolytes at ionic strengths up to 1 M. It is found that ion-image charge interactions and interionic correlations, which are neglected by the PB equation, have relatively weak effects on the electrostatic potential at charged groups of the DNA. The PB equation predicts errors in the long-range electrostatic part of the free energy that are only ∼1.5 kJ/mol per nucleotide even in the case of an asymmetrical electrolyte. In contrast, the spatial correlations between ions drastically affect the electrostatic potential at significant separations from the macromolecule leading to a clearly predicted effect of charge overneutralization.     

Interpreting protein/DNA interactions: distinguishing specific from non-specific and electrostatic from non-electrostatic components

We discuss the effectiveness of existing methods for understanding the forces driving the formation of specific protein-DNA complexes. Theoretical approaches using the Poisson-Boltzmann (PB) equation to analyse interactions between these highly charged macromolecules to form known structures are contrasted with an empirical approach that analyses the effects of salt on the stability of these complexes and assumes that release of counter-ions associated with the free DNA plays the dominant role in their formation. According to this counter-ion condensation (CC) concept, the salt-dependent part of the Gibbs energy of binding, which is defined as the electrostatic component, is fully entropic and its dependence on the salt concentration represents the number of ionic contacts present in the complex. It is shown that although this electrostatic component provides the majority of the Gibbs energy of complex formation and does not depend on the DNA sequence, the salt-independent part of the Gibbs energy--usually regarded as non-electrostatic--is sequence specific. The CC approach thus has considerable practical value for studying protein/DNA complexes, while practical applications of PB analysis have yet to demonstrate their merit.     

Origin of DNA helical structure and its sequence dependence

Conformational analysis of DNA shows that the origin of the B-form double helix can be attributed in large part to the atomic charge pattern in the base pairs. The charge patterns favor specific helical stacking of the base pairs. Base pairs alone--without backbones--have a strong tendency to form helix, indicating that the backbones play a rather passive role in determining the basic helical structure of DNA. It is mainly the electrostatic interactions determined by the charge pattern on base pairs that stabilize a particular helical conformation. The charge pattern in the base pairs appears to be responsible for much of the sequence dependence of DNA conformation, rather than steric clashes.     

Electrostatic field of the large fragment of Escherichia coli DNA polymerase I

The electrostatic field of the large fragment of Escherichia coli DNA polymerase I (Klenow fragment) has been calculated by the finite difference procedure on a 2 A grid. The potential field is substantially negative at physiological pH (reflecting the net negative charge at this pH). The largest regions of positive potential are in the deep crevice of the C-terminal domain, which is the proposed binding site for the DNA substrate. Within the crevice, the electrostatic potential has a partly helical form. If the DNA is positioned to fulfil stereochemical requirements, then the positive potential generally follows the major groove and (to a lesser extent) the negative potential is in the minor groove. Such an arrangement could stabilize DNA configurations related by screw symmetry. The histidine residues of the Klenow fragment give the positive field of the groove a sensitivity to relatively small pH changes around neutrality. We suggest that the histidine residues could change their ionization states in response to DNA binding, and that this effect could contribute to the protein-DNA binding energy.     

Association kinetics with coupled diffusion: III. Ionic-strength dependence of the lac repressor-operator association

The repressor- operator association is treated in a model where the represser molecule can find its specific binding site the operator, on a large DNA chain by performing a one-dimensional diffusion along the chain. The ionic-strength depen deuce is calculated by introducing a screened electrostatic potential around the DNA chain and coupling the free diffusion of the represser in this potential to the proposed one-dimensional diffusion along the chain. The main influence on the association rate comes from the competitive binding of ions to the unspecific DNA sites. It is also demonstrated that during the time that the represser is bound in a global sense, the diffusion along the chain will be made up of a strictly one-dimensional motion over fairly short distances, interspersed with many local dissociations during which the repressor in essence is free in solution.     

Sequence-dependent binding of metal ion to DNA oligomeres. A comparison of molecular electrostatic potentials with NMR data

Experimentally observed sequence-selective binding of metal ion to DNA oligonucleotides have been compared with variations of electrostatic potential (EP) along the helix. Calculations of EP have been performed for three atomic models of the oligonucleotide duplex [d(CGCGAATTCGCG)2] using several variants of EP calculations, including a solution of non-linear Poisson-Boltzmann equation (NPBE). N7 atom of guanine adjacent to adenine base was identified as a region with the most negative electrostatic potential in the major groove. The EP value for the Me ion binding site surpasses the value for N7 of other guanines by 10-26% depending on particular duplex conformation. Qualitatively, the sequence dependent variations of EP near guanine N7 atoms are in agreement with the sequence-selective behavior of Mn(II) and Zn(II) ions as revealed by NMR experiments. But the difference in EP between the two most negative regions near guanine N7 atoms does not exceed 1.25 kT/e. Simple model suggests that metal ions are capable to form ion-hydrate complexes with G-Pu steps of DNA duplex. These complexes are formed via one Me...G and five Me...water coordination bonds with water molecules hydrogen bonded to two adjacent purine bases in the same chain. We suppose that such a stereospecific structural possibility is the main factor which control the sequence-selectivity in the metal ion binding. A combination of both mechanisms allows to explain sequence specific Mn(II) and Zn(II) binding to a set of oligonucleotides.     

Dynamics of water and ions around DNA: What is so special about them?

Water around biomolecules is special for behaving strangely – both in terms of structure and dynamics, while ions are found to control various interactions in biomolecules such as DNA, proteins and lipids. The questions that how water and ions around these biomolecules behave in terms of their structure and dynamics, and how they affect the biomolecular functions have triggered tremendous research activities worldwide. Such activities not only unfolded important static and dynamic properties of water and ions around these biomolecules, but also provoked heated debate regarding their explanation and role in biological functions. DNA, being negatively charged, interacts strongly with the surrounding dipolar water and positively charged counterions, leading to complex electrostatic coupling of water and ions with the DNA. Recent time-resolved fluorescence Stokes shift experiments and related computer simulation studies from our and other laboratories have unfolded some unique dynamic characteristics of water and ions near different structures of DNA. These results are discussed here to showcase the specialty of water and ion dynamics around DNA.     

Onset of DNA aggregation in presence of monovalent and multivalent counterions

We address theoretically aggregation of DNA segments by multivalent polyamines such as spermine and spermidine. In experiments, the aggregation occurs above a certain threshold concentration of multivalent ions. We demonstrate that the dependence of this threshold on the concentration of DNA has a simple form. When the DNA concentration c(DNA) is smaller than the monovalent salt concentration, the threshold multivalent ion concentration depends linearly on c(DNA), having the form alphac(DNA) + beta. The coefficients alpha and beta are related to the density profile of multivalent counterions around isolated DNA chains, at the onset of their aggregation. This analysis agrees extremely well with recent detailed measurements on DNA aggregation in the presence of spermine. From the fit to the experimental data, the number of condensed multivalent counterions per DNA chain can be deduced. A few other conclusions can then be reached: 1), the number of condensed spermine ions at the onset of aggregation decreases with the addition of monovalent salt; 2), the Poisson-Boltzmann theory overestimates the number of condensed multivalent ions at high monovalent salt concentrations; and 3), our analysis of the data indicates that the DNA charge is not overcompensated by spermine at the onset of aggregation.     

Biological Macromolecular Dynamics

Abstract

Experimentally observed sequence-selective binding of metal ion to DNA oligonucleotides have been compared with variations of electrostatic potential (EP) along the helix. Calculations of EP have been performed for three atomic models of the oligonucleotide duplex [d(CGCGAATTCGCG)2] using several variants of EP calculations, including a solution of non-linear Poisson-Boltzmann equation (NPBE). N7 atom of guanine adjacent to adenine base was identified as a region with the most negative electrostatic potential in the major groove. The EP value for the Me ion binding site surpasses the value for N7 of other guanines by 10–26% depending on particular duplex conformation. Qualitatively, the sequence dependent variations of EP near guanine N7 atoms are in agreement with the sequence-selective behavior of Mn(II) and Zn(II) ions as revealed by NMR experiments. But the difference in EP between the two most negative regions near guanine N7 atoms does not exceed 1.25 kT/e. Simple model suggests that metal ions are capable to form ion-hydrate complexes with G-Pu steps of DNA duplex. These complexes are formed via one Me…G and five Me…water coordination bonds with water molecules hydrogen bonded to two adjacent purine bases in the same chain. We suppose that such a stereospecific structural possibility is the main factor which control the sequence-selectivity in the metal ion binding. A combination of both mechanisms allows to explain sequence specific Mn(II) and Zn(II) binding to a set of oligonucleotides.     

Brownian dynamics simulations of ion atmospheres around polyalanine and B-DNA: effects of biomolecular dielectric

We have extended an earlier Brownian dynamics simulation algorithm for simulating the structural dynamics of ions around biomolecules to accommodate dielectric inhomogeneity. The electrostatic environment of a biomolecule immersed in water was obtained by numerically solving the Poisson equation with the biomolecule treated as a low dielectric region and the solvent treated as a high dielectric region. Instead of using the mean-field type approximations of ion interactions as in the Poisson-Boltzmann model, the ions were treated explicitly by allowing them to evolve dynamically under the electrostatic field of the biomolecule. This model thus accounts for ion-ion correlations and the finite-size effects of the ions. For a 13-residue alpha-helical polyalanine and a 12-base-pair bp B-form DNA, we found that the choice of the dielectric constant of the biomolecule has much larger effects on the mean ionic structure around the biomolecule than on the fluctuational and dynamical properties of the ions surrounding the biomolecule.     

DNA minor groove electrostatic potential: influence of sequence-specific transitions of the torsion angle gamma and deoxyribose conformations

The structural adjustments of the sugar-phosphate DNA backbone (switching of the γ angle (O5′–C5′–C4′–C3′) from canonical to alternative conformations and/or C2′-endo → C3′-endo transition of deoxyribose) lead to the sequence-specific changes in accessible surface area of both polar and non-polar atoms of the grooves and the polar/hydrophobic profile of the latter ones. The distribution of the minor groove electrostatic potential is likely to be changing as a result of such conformational rearrangements in sugar-phosphate DNA backbone. Our analysis of the crystal structures of the short free DNA fragments and calculation of their electrostatic potentials allowed us to determine: (1) the number of classical and alternative γ angle conformations in the free B-DNA; (2) changes in the minor groove electrostatic potential, depending on the conformation of the sugar-phosphate DNA backbone; (3) the effect of the DNA sequence on the minor groove electrostatic potential. We have demonstrated that the structural adjustments of the DNA double helix (the conformations of the sugar-phosphate backbone and the minor groove dimensions) induce changes in the distribution of the minor groove electrostatic potential and are sequence-specific. Therefore, these features of the minor groove sizes and distribution of minor groove electrostatic potential can be used as a signal for recognition of the target DNA sequence by protein in the implementation of the indirect readout mechanism.     

Interaction of a G-DNA quadruplex with mono- and divalent cations. A force field calculation

The formation and stability of four-stranded DNA in solution is specifically dependent on the type of cations present. The interaction potential of a model quadruplex structure with different mono- and divalent ions was determined by force field calculations. Though the electrostatic contribution to the total energy is mainly responsible for the stabilisation of the cations within the quadruplex channel, it is the van der Waals interaction at short distances that determines the specific characteristics of the different cations. An explicit consideration of the solvent indicates that the position of water molecules in close proximity to the DNA channel have a strong influence on the form of the potential, and hence on the capability of the cations for leaving and re-entering the cavity. The effect of cation size, as expressed through their Lennard-Jones parameters, is discussed.     

Hydration effects on the electrostatic potential around tuftsin

The electrostatic potential and component dielectric constants from molecular dynamics (MD) trajectories of tuftsin, a tetrapeptide with the amino acid sequence Thr–Lys–Pro–Arg in water and in saline solution are presented. The results obtained from the analysis of the MD trajectories for the total electrostatic potential at points on a grid using the Ewald technique are compared with the solution to the Poisson–Boltzmann (PB) equation. The latter was solved using several sets of dielectric constant parameters. The effects of structural averaging on the PB results were also considered. Solute conformational mobility in simulations gives rise to an electrostatic potential map around the solute dominated by the solute monopole (or lowest order multipole). The detailed spatial variation of the electrostatic potential on the molecular surface brought about by the compounded effects of the distribution of water and ions close to the peptide, solvent mobility, and solute conformational mobility are not qualitatively reproducible from a reparametrization of the input solute and solvent dielectric constants to the PB equation for a single structure or for structurally averaged PB calculations. Nevertheless, by fitting the PB to the MD electrostatic potential surfaces with the dielectric constants as fitting parameters, we found that the values that give the best fit are the values calculated from the MD trajectories. Implications of using such field calculations on the design of tuftsin peptide analogues are discussed. © 1999 John Wiley & Sons, Inc. Biopoly 50: 133–143, 1999     

Ions around DNA: Monte Carlo estimates of distribution with improved electrostatic potentials

Coulomb's law for the electrostatic interactions between ions is modified when discontinuities in dielectric constant (relative permittivity) occur. In a DNA solution such a discontinuity occurs at the interface between the DNA molecular helix and the surrounding water. We take the modified interaction potentials from a previous report [Macromolecules (1986) 19 , 1186–1194] and use them with the Monte Carlo method to find the distribution of univalent and bivalent counterions around the DNA helix in the absence of coions (i.e., no added salt). In comparing the ion distribution with the modified potential to the distribution without, we find that the effects of the modifications to the potentials are considerable. The modifications tend to drive the ions out of the grooves of the helix, especially out of the major groove. This result comes partly from the repulsion exerted on the ions by the low-permittivity helix and partly from the concentration of the field of the phosphates at the surface of the helix, a concentration that is also caused by the discontinuity in permittivity.     

The effect of a variable dielectric coefficient and finite ion size on Poisson-Boltzmann calculations of DNA-electrolyte systems.

The results of variable dielectric coefficient Poisson-Boltzmann calculations of the counter-ion concentration in the vicinity of an all-atom model of the B-form of DNA are presented with an emphasis on the importance of spatial variations in the dielectric properties of the solvent, particularly at the macro-ion-solvent interface. Calculations of the distribution of hard-sphere electrolyte ions of various dimensions are reported. The presence of a dielectric boundary significantly increases the magnitude of the electrostatic potential with a concomitant increase in the accumulation of small counter-ions in the groove regions of DNA. Because ions with radii greater than 2 A have restricted access to the minor groove, the effect there is less significant than it is within the major groove. Changes in the dielectric coefficient for the electrolyte solution, allowing variation from 10 to 25, 40, 60, and 78.5 within the first 7.4 A of the surface of DNA, substantially increases the calculated surface concentration of counter-ions of all sizes. A lower dielectric coefficient near the macro-ion surface also tends to increase the counter-ion density in regions where the electrostatic potential is more negative than -kT. Regardless of the choice of dielectric coefficient, the number of ions in regions where the electrostatic potential is less than -kT remains the same for 0.153 M added 1-1 monovalent electrolyte as for the case without added salt.(ABSTRACT TRUNCATED AT 250 WORDS)     

Binding of DNA to large fragment of DNA polymerase I: identification of strong and weak electrostatic forces and their biological implications.

Examination of the electrostatic potential of a modeled complex, consisting of the Klenow fragment of E. coli DNA polymerase I and DNA template-primer, suggested the presence of two distinct interacting regions. The one displaying a strong electropositive potential field is generated by side chains of basic amino acid pairs and is directed towards the major groove site in DNA. The second electrostatic potential field around DNA is somewhat weaker and appears to be exerted by a pair of vicinal side chains of acidic and basic amino acids. The distribution of charges in this manner appears well suited for the binding of enzyme to the template-primer required in the enzymatic synthesis of DNA.