Calculations of the spatial distribution of charge density in the environment of DNA

A technique for modeling the structured environmental charge distribution about isolated polyions of arbitrary geometry is presented and applied to B-DNA. It describes the three-dimensional variation of the continuous space charge and allows estimation of local electrostatic potentials and fields that the electrolytic environment induces at nuclei of the polyion. Calculations involve an iterative solution to the set of equations coupling electrostatic potential and average charge density in space. By dividing the region around a DNA segment into finite volume elements, sets of numerically stable atmospheric charge densities have been obtained over a range of concentrations of added monovalent salt. Results are in good agreement with those of Poisson-Boltzmann calculations on comparable systems and are consistent with findings from Monte Carlo simulations of DNA.     

Local conformational variations observed in B-DNA crystals do not improve base stacking: computational analysis of base stacking in a d(CATGGGCCCATG)(2) B<-->A intermediate crystal structure

The crystal structure of d(CATGGGCCCATG)₂ shows unique stacking patterns of a stable B↔A-DNA intermediate. We evaluated intrinsic base stacking energies in this crystal structure using an ab initio quantum mechanical method. We found that all crystal base pair steps have stacking energies close to their values in the standard and crystal B-DNA geometries. Thus, naturally occurring stacking geometries were essentially isoenergetic while individual base pair steps differed substantially in the balance of intra-strand and inter-strand stacking terms. Also, relative dispersion, electrostatic and polarization contributions to the stability of different base pair steps were very sensitive to base composition and sequence context. A large stacking flexibility is most apparent for the CpA step, while the GpG step is characterized by weak intra-strand stacking. Hydration effects were estimated using the Langevin dipoles solvation model. These calculations showed that an aqueous environment efficiently compensates for electrostatic stacking contributions. Finally, we have carried out explicit solvent molecular dynamics simulation of the d(CATGGGCCCATG)₂ duplex in water. Here the DNA conformation did not retain the initial crystal geometry, but moved from the B↔A intermediate towards the B-DNA structure. The base stacking energy improved in the course of this simulation. Our findings indicate that intrinsic base stacking interactions are not sufficient to stabilize the local conformational variations in crystals.     

Electrostatic effects on the stability of condensed DNA in the presence of divalent cations.

Cylindrical cell model Poisson-Boltzmann (P-B) calculations are used to evaluate the electrostatic contributions to the relative stability of various DNA conformations (A, B, C, Z, and single-stranded (ss) with charge spacings of 3.38 and 4.2 A) as a function of interhelix distance in a concentrated solution of divalent cations. The divalent ion concentration was set at 100 mM, to compare with our earlier reports of spectroscopic and calorimetric experiments, which demonstrate substantial disruption of B-DNA geometry. Monovalent cations neutralize the DNA phosphates in two ways, corresponding to different experimental situations: 1) There is no significant contribution to the ionic strength from the neutralizing cations, corresponding to DNA condensation from dilute solution and to osmotic stress experiments in which DNA segments are brought into close proximity to each other in the presence of a large excess of buffer. 2) The solution is uniformly concentrated in DNA, so that the neutralizing cations add significantly to those in the buffer at close DNA packing. In case 1), conformations with lower charge density (Z and ssDNA) have markedly lower electrostatic free energies than B-DNA as the DNA molecules approach closely, due largely to ionic entropy. If the divalent cations bind preferentially to single-stranded DNA or a distorted form of B-DNA, as is the case with transition metals, the base pairing and stacking free energies that stabilize the double helix against electrostatic denaturation may be overcome. Strong binding to the bases is favored by the high concentration of divalent cations at the DNA surface arising from the large negative surface potential; the surface concentration increases sharply as the interhelical distance decreases. In case 2), the concentration of neutralizing monovalent cations becomes very large and the electrostatic free energy difference between secondary structures becomes small as the interhelical spacing decreases. Such high ionic concentrations will be expected to modify the stability of DNA by changing water activity as well as by screening electrostatic interactions. This may be the root of the decreased thermal stability of DNA in the presence of high concentrations of magnesium ions.     

B-DNA structure is intrinsically polymorphic: even at the level of base pair positions

Increasingly exact measurement of single crystal X-ray diffraction data offers detailed characterization of DNA conformation, hydration and electrostatics. However, instead of providing a more clear and unambiguous image of DNA, highly accurate diffraction data reveal polymorphism of the DNA atomic positions and conformation and hydration. Here we describe an accurate X-ray structure of B-DNA, painstakingly fit to a multistate model that contains multiple competing positions of most of the backbone and of entire base pairs. Two of ten base-pairs of CCAGGCCTGG are in multiple states distinguished primarily by differences in slide. Similarly, all the surrounding ions are seen to fractionally occupy discrete competing and overlapping sites. And finally, the vast majority of water molecules show strong evidence of multiple competing sites. Conventional resolution appears to give a false sense of homogeneity in conformation and interactions of DNA. In addition, conventional resolution yields an average structure that is not accurate, in that it is different from any of the multiple discrete structures observed at high resolution. Because base pair positional heterogeneity has not always been incorporated into model-building, even some high and ultrahigh-resolution structures of DNA do not indicate the full extent of conformational polymorphism.     

Sequence and temperature effect on hydrogen bond disruption in DNA determined by a statistical analysis

We use the modified self-consistent phonon approximation theory to calculate temperature dependent interbase hydrogen bond disruption profiles for a number of six base pair repeating sequence infinite B-DNA polymers with various guanine-cytosine/adenine-thymine ratios. For comparison we also include results we have obtained in our earlier work on several B-DNA homopolymers, copolymers and a four-base-pair repeating sequence polymer. Our theory gives a statistical estimate of thermal fluctuational disruption probability of individual hydrogen bonds in individual base pairs in DNA as a function of temperature. The calculated probabilities show no sequence dependence at premelting temperatures, in agreement with proton exchange measurements. These probabilities however become very sensitive to base sequence at temperatures close to the observed melting temperatures. Multi-phasic critical transitions are found in which a portion of base pairs are disrupted at temperatures below the final disruption temperature. These transitions include localized as well as non-localized base pair opening. The localized transitions involve disruption of a few base-pairs at every other location without large scale base unstacking, and they may not appear in the observed UV curves with current resolution. On the other hand the overall disruption behavior is consistent with observations. The midpoint transition temperatures are close to the observed melting temperatures and these temperatures show the observed linear dependence on guanine-cytosine content. Our calculations indicate that our theory can be used effectively to calculate H-bond disruption behavior of different DNA sequences. Received: 20 February 1996 / Accepted: 2 May 1996     

Nonplanar DNA Base Pairs

Abstract

Three-dimensional structures of a representative set of more than 30 hydrogen-bonded nucleic acids base pairs have been studied by reliable ab initio quantum mechanical methods. We show that many hydrogen-bonded nucleic acid base pairs are intrinsically nonplanar, mainly due to the partial sp³ hybridization of nitrogen atoms of their amino groups and secondary electrostatic interactions. This finding extends the variability of intermolecular interactions of DNA bases in that i) flexibility of the base pairs is larger than has been assumed before, and ii) attractive proton-proton acceptor interactions oriented out of the base pair plane are allowed. For example, all four G…A mismatch base pairs are propeller twisted, and the energy preferences for the nonplanar structures range from less than 0.1 kcal/mol to 1.8 kcal/mol. We predict that nonplanarity of the amino group of guanine in the G(anti)…A(anti) pair of the ApG step of the d(CCAAGATTGG)₂ crystal structure is an important stabilizing factor that improves the energy of this structure by almost 3 kcal/mol. Currently used empirical potentials are not accurate enough to properly cover the interactions associated with amino-group and base-pair nonplanarity.     

The electrostatic field of DNA: the role of the nucleic acid conformation.

Calculations of the electrostatic field of DNA in two very different double helical conformations, A and Z, are reported and compared with the results previously obtained for B-DNA. Striking contrasts between these fields and the associated electrostatic potentials are brought into evidence. One of the major differences is that while the deepest potentials are generally located in the grooves of DNA, the strongest fields are associated with the phosphate groups. The results of screening the nucleic acids by counterions are also presented.     

The dependence of the surface electrostatic potential of B-DNA on environmental factors

The electrostatic potential of B-DNA is calculated on its surface envelope for two homopolymeric base pair sequences using models representing the effects of both counterion binding and of aqueous solution. The influence of these two factors on the resulting potentials is established and the significance of calculations which omit such effects is discussed.     

Self-consistent reaction field computation of the reactive characteristics of DNA bases in water

The effect of the solvent in the molecular properties of DNA bases has been explored by using a self-consistent reaction field (SCRF) method based on the AM1 (Austin Model 1) Hamiltonian and a modified version of the high level Miertus–Scrocco–Tomasi (MST) algorithm. MST/AM1 estimates of free energies of hydration compare qualitatively well with the available experimental data, as well as with the results obtained from molecular dynamic simulations. Furthermore, the changes in the dipole predicted by the MST/AM1 method are in good agreement with Monte Carlo/quantum mechanical data, as well as with AM1-SM2 (Soluation Model 2) estimates. AM1/MST calculations of Mulliken, and electrostatic charges, dipoles, molecular electrostatic potentials and molecular interaction potentials in both vacuum and solution allowed us to quantify the effect of the water on the reactive characteristics of the DNA bases. This effect is large and complex, and cannot be neglected in theoretical calculations where an accurate representation of the DNA bases is needed. The possibility of including the polarization effect of the water into force-field simulations of DNA structures is discussed. © 1993 John Wiley & Sons, Inc.     

Characterization of the base stacking interactions in DNA by means of Lennard-Jones empirical potentials

Three empirical potentials of the Lennard-Jones type taken from literature were used to calculate van der Waals contributions to the base-pair couples stacking energies in B-DNA and A-DNA type double helical conformations. The information obtained can be summarized as follows: (1) Purine-pyrimidine and purine-purine (pyrimidine-pyrimidine in the complementary strand) sequences preferred right-handed helical arrangement, whereas pyrimidine-purine sequences favoured left-handed (C-G) or unwound (T-A) stacking geometry; in the latter case this only held for B- but not A-DNA (the C-G sequence was not studied in A-DNA owing to difficulties (see below) with the G amino group in B-DNA); (2) Positive propeller twist of base-pairs was stable in both B- and A-DNA; the thymine methyl group promoted the propeller and this effect was strongest in the A-T step; (3) Tilt of base pairs occurred around zero in B-DNA and between 15-20 degrees C in A-DNA, in agreement with the experimental observations; (4) Vertical separation of base pairs was optimal within 0.33-0.34 nm for B-DNA and around 0.29 nm for A-DNA using the 9-6 potential. The 12-6 potential gave similar results with B-DNA as the 9-6 potential if, however, base pairs were separated by 0.35-0.36 nm; (5) The calculated effect of the guanine amino group was substantially stronger than expected on the basis of data derived from X-ray diffraction studies of oligonucleotide single crystals; (6) In comparison with the 9-6 potential, the 12-6 potential provided more strict energy minima. In summary, the empirical potentials reproduce, at least semiquantitatively, many but not all DNA properties; this should be taken into account whenever the potentials are used for prediction purposes.     

Influence of divalent magnesium ion on DNA: molecular dynamics simulation studies

A large amount of experimental evidence is available on the effect of magnesium ions on the structure and stability of DNA double helix. Less is known, however, on how these ions affect the stability and dynamics of the molecule. The static time average pictures from X-ray structures or the quantum chemical energy minimized structures lack understanding of the dynamic DNA–ion interaction. The present work addresses these questions by molecular dynamics simulation studies on two DNA duplexes and their interaction with magnesium ions. Results show typical B-DNA character with occasional excursions to deviated states. We detected expected stability of the duplexes in terms of backbone conformations and base pair parameter by the CHARMM-27 force field. Ion environment analysis shows that Mg²⁺ retains the coordination sphere throughout the simulation with a preference for major groove over minor. An extensive analysis of the influence of the Mg²⁺ ion shows no evidence of the popular predictions of groove width narrowing by dipositive metal ion. The major groove atoms show higher occupancy and residence time compared to minor groove for magnesium, where no such distinction is found for the charge neutralizing Na⁺ ions. The determining factor of Mg²⁺ ion’s choice in DNA binding site evolves as the steric hindrance faced by the bulky hexahydrated cation where wider major groove gets the preference. We have shown that in case of binding of Mg²⁺ to DNA non electrostatic contributions play a major role.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:5     

Electrostatic deformation of DNA by a DNA-binding protein

Complementary electrostatic interactions between negatively charged B-DNA and a positively charged array on the lambda Cro repressor protein are shown to substantially contribute to the formation energy of sequence-specific and nonspecific Cro-DNA complexes. The electrostatic interactions favor Cro binding to a bent form of DNA, a geometry which optimizes hydrogen-bonding contacts between Cro and exposed base pair groups in the DNA major groove.     

Electrostatic interactions across a charged lipid bilayer

We present theoretical work in which the degree of electrostatic coupling across a charged lipid bilayer in aqueous solution is analyzed on the basis of nonlinear Poisson–Boltzmann theory. In particular, we consider the electrostatic interaction of a single, large macroion with the two apposed leaflets of an oppositely charged lipid bilayer where the macroion is allowed to optimize its distance to the membrane. Three regimes are identified: a weak and a high macroion charge regime, separated by a regime of close macroion–membrane contact for intermediate charge densities. The corresponding free energies are used to estimate the degree of electrostatic coupling in a lamellar cationic lipid–DNA complex. That is, we calculate to what extent the one-dimensional DNA arrays in a sandwich-like lipoplex interact across the cationic membranes. We find that, in spite of the low dielectric constant inside a lipid membranes, there can be a significant electrostatic contribution to the experimentally observed cross-bilayer orientational ordering of the DNA arrays. Our approximate analytical model is complemented and supported by numerical calculations of the electrostatic potentials and free energies of the lamellar lipoplex geometry. To this end, we solve the nonlinear Poisson–Boltzmann equation within a unit cell of the lamellar lipoplex using a new lattice Boltzmann method. Dedicated to Prof. K. Arnold on the occasion of his 65th birthday.     

Mediation of the A/B-DNA helix transition by G-tracts in the crystal structure of duplex CATGGGCCCATG

The crystal structure of the DNA dodecamer duplex CATGGGCCCATG lies on a structural continuum along the transition between A- and B-DNA. The dodecamer possesses the normal vector plot and inclination values typical of B-DNA, but has the crystal packing, helical twist, groove width, sugar pucker, slide and x-displacement values typical of A-DNA. The structure shows highly ordered water structures, such as a double spine of water molecules against each side of the major groove, stabilizing the GC base pairs in an A-like conformation. The different hydration of GC and AT base pairs provides a physical basis for solvent-dependent facilitation of the A↔B helix transition by GC base pairs. Crystal structures of CATGGGCCCATG and other A/B-DNA intermediates support a ‘slide first, roll later’ mechanism for the B→A helix transition. In the distribution of helical parameters in protein–DNA crystal structures, GpG base steps show A-like properties, reflecting their innate predisposition for the A conformation.     

The Dependence of the Surface Electrostatic Potential of B-DNA on Environmental Factors

Abstract

The electrostatic potential of B-DNA is calculated on its surface envelope for two homopolymeric base pair sequences using models representing the effects of both counterion binding and of aqueous solution. The influence of these two factors on the resulting potentials is established and the significance of calculations which omit such effects is discussed.     

Structural and energetic consequences of oxidation of d(ApGpGpGpTpT) telomere repeat unit in complex with TRF1 protein

The configuration hyperspace of canonical and oxidized 14-mers of B-DNA comprising telomere repeat units d(ApGpGpGpTpT) was sampled over 40 ns via molecular dynamic (MD) simulations. The energetic and structural consequences of TRF1 binding to telomere B-DNA were compared with non-complexed systems. Energetic properties of analyzed pairs, di- and tri-nucleotide steps occurring in central telomere repeat unit were estimated by means of advanced quantum chemistry computations including not only BSSE corrections, electron correlation contributions but also non-negligible many-body terms. These data along with bases pair and base step parameters distributions allow for quantization of consequences of oxidation and/or TRF1 binding to telomere repeat units. Occurrence of 8-oxoguanine in central telomeric triad (CTT) is the source of high stiffness if compared to non-modified oligomer. The origin of this property comes from significantly alteration of intermolecular interactions introduced by 8-oxoguanine. The increased stability observed for base–base interactions are accumulated and characterizes also di- and tri-nucleotides. The observed changes in the intermolecular interactions originate from structural alterations imposed by TRF1 binding to canonical and oxidized telomere B-DNA. First and most direct consequence of TRF1 binding to oxidized telomere repeat unit is alteration of shift-slide correlations if compared to canonical system. This in turn leads to large differences in purine-purine overlapping in oxidized structures. Thus, oxidized telomere B-DNA double strands are sensitive to interactions with protein ligands and numerous structural and energetic changes are imposed on base pairs forming CTT.     

Effect of bulky lesions on DNA: solution structure of a DNA duplex containing a cholesterol adduct

The three-dimensional solution structure of two DNA decamers of sequence d(CCACXGGAAC)-(GTTCCGGTGG) with a modified nucleotide containing a cholesterol derivative (X) in its C1 '(chol)alpha or C1 '(chol)beta diastereoisomer form has been determined by using NMR and restrained molecular dynamics. This DNA derivative is recognized with high efficiency by the UvrB protein, which is part of the bacterial nucleotide excision repair, and the alpha anomer is repaired more efficiently than the beta one. The structures of the two decamers have been determined from accurate distance constraints obtained from a complete relaxation matrix analysis of the NOE intensities and torsion angle constraints derived from J-coupling constants. The structures have been refined with molecular dynamics methods, including explicit solvent and applying the particle mesh Ewald method to properly evaluate the long range electrostatic interactions. These calculations converge to well defined structures whose conformation is intermediate between the A- and B-DNA families as judged by the root mean square deviation but with sugar puckerings and groove shapes corresponding to a distorted B-conformation. Both duplex adducts exhibit intercalation of the cholesterol group from the major groove of the helix and displacement of the guanine base opposite the modified nucleotide. Based on these structures and molecular dynamics calculations, we propose a tentative model for the recognition of damaged DNA substrates by the UvrB protein.     

Free energy calculation on base specificity of drug--DNA interactions: application to daunomycin and acridine intercalation into DNA

We present the results of free energy perturbation/molecular dynamics studies on B-DNA.daunomycin and B-DNA.9-aminoacridine complexes as well as on B-DNA itself in order to calculate the free energy differences between complexes having different base pair sequences. The results generally reproduce the trends observed experimentally, i.e., preferences of acridine and daunomycin to bind to a specific base sequence in the DNA. This is encouraging, given the simplicity of the molecular mechanical/dynamical model in which solvent is not explicitly included.     

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.