Left-handed DNA crossovers. Implications for DNA-DNA recognition and structural alterations

The close approach of DNA segments participates in many biological functions including DNA condensation and DNA processing. Previous crystallographic studies have shown that B-DNA self-fitting by mutual groove-backbone interaction produces right-handed DNA crossovers. These structures have opened new perspectives on the role of close DNA-DNA interactions in the architecture and activity the DNA molecule. In the present study, the analysis of the crystal packing of two B-DNA decamer duplexes d(CCIIICCCGG) and d(CCGCCGGCGG) reveals the existence of new modes of DNA crossing. Symmetric left-handed crossovers are produced by mutual fitting of DNA grooves at the crossing point. New sequence patterns contribute to stabilize longitudinal fitting of the sugar-phosphate backbone into the major groove. In addition, the close approach of DNA segments greatly influences the DNA conformation in a sequence dependent manner. This study provides new insights into the role of DNA sequence and structure in DNA-DNA recognition. In providing detailed molecular views of DNA crossovers of opposite chirality, this study can also help to elucidate the role of symmetry and chirality in the recognition of complex DNA structures by protein dimers or tetramers, such as topoisomerase II and recombinase enzymes. These results are discussed in the context of the possible relationships between DNA condensation and DNA processing.     

Helical Chirality: a Link between Local Interactions and Global Topology in DNA

DNA supercoiling plays a major role in many cellular functions. The global DNA conformation is however intimately linked to local DNA-DNA interactions influencing both the physical properties and the biological functions of the supercoiled molecule. Juxtaposition of DNA double helices in ubiquitous crossover arrangements participates in multiple functions such as recombination, gene regulation and DNA packaging. However, little is currently known about how the structure and stability of direct DNA-DNA interactions influence the topological state of DNA. Here, a crystallographic analysis shows that due to the intrinsic helical chirality of DNA, crossovers of opposite handedness exhibit markedly different geometries. While right-handed crossovers are self-fitted by sequence-specific groove-backbone interaction and bridging Mg²⁺ sites, left-handed crossovers are juxtaposed by groove-groove interaction. Our previous calculations have shown that the different geometries result in differential stabilisation in solution, in the presence of divalent cations. The present study reveals that the various topological states of the cell are associated with different inter-segmental interactions. While the unstable left-handed crossovers are exclusively formed in negatively supercoiled DNA, stable right-handed crossovers constitute the local signature of an unusual topological state in the cell, such as the positively supercoiled or relaxed DNA. These findings not only provide a simple mechanism for locally sensing the DNA topology but also lead to the prediction that, due to their different tertiary intra-molecular interactions, supercoiled molecules of opposite signs must display markedly different physical properties. Sticky inter-segmental interactions in positively supercoiled or relaxed DNA are expected to greatly slow down the slithering dynamics of DNA. We therefore suggest that the intrinsic helical chirality of DNA may have oriented the early evolutionary choices for DNA topology.     

Differential stability of DNA crossovers in solution mediated by divalent cations

The assembly of DNA duplexes into higher-order structures plays a major role in many vital cellular functions such as recombination, chromatin packaging and gene regulation. However, little is currently known about the molecular structure and stability of direct DNA–DNA interactions that are required for such functions. In nature, DNA helices minimize electrostatic repulsion between double helices in several ways. Within crystals, B-DNA forms either right-handed crossovers by groove–backbone interaction or left-handed crossovers by groove–groove juxtaposition. We evaluated the stability of such crossovers at various ionic concentrations using large-scale atomistic molecular dynamics simulations. Our results show that right-handed DNA crossovers are thermodynamically stable in solution in the presence of divalent cations. Attractive forces at short-range stabilize such crossover structures with inter-axial separation of helices less than 20 Å. Right-handed crossovers, however, dissociate swiftly in the presence of monovalent ions only. Surprisingly, left-handed crossovers, assembled by sequence-independent juxtaposition of the helices, appear unstable even at the highest concentration of Mg²⁺studied here. Our study provides new molecular insights into chiral association of DNA duplexes and highlights the unique role divalent cations play in differential stabilization of crossover structures. These results may serve as a rational basis to understand the role DNA crossovers play in biological processes.     

Does topoisomerase II specifically recognize and cleave hairpins,cruciforms and crossovers of DNA?

DNA topoisomerase II is an enzyme that specializes in DNA disentanglement. It catalyzes the interconversion of DNA between different topological states. This event requires the passage of one duplex through another one via a transient double-strand break. Topoisomerase II is able to process any type of DNA, including structures such as DNA juxtapositions (crossovers), DNA hairpins or cruciforms, which are recognized with high specificity. In this review, we focused our attention on topoisomerase II recognizing DNA substrates that possess particular geometries. A strong cleavage site, as we identified in pBR322 DNA in the presence of ellipticine (site 22), appears to be characterized by a cruciform structure formed from two stable hairpins. The same sequence could also constitute a four-way junction structure stabilized by interactions involving ATC sequences. The latter have been shown to be able to promote Holliday junctions. We reviewed the recent literature that deals with the preferential recognition of crossovers by various topoisomerases. The single molecule relaxation experiments have demonstrated the differential abilities of the topoisomerases to recognize crossovers. It appears that enzymes, which distinguish the chirality of the crossovers, possess specialized domains dedicated to this function. We also stress that the formation of crossovers is dependent on the presence of adequate stabilizing sequences. Investigation of the impact of such structures on enzyme activity is important in order to both improve our knowledge of the mechanism of action of the topoisomerase II and to develop new inhibitors of this enzyme.     

Secondary conformational polymorphism of nucleic acids as a possible functional link between cellular parameters and DNA packaging processes

Circular dichroism and electron microscopy studies of various in vitro DNA packaging systems indicate that all the factors which induce and modulate the secondary conformation of DNA molecules are capable of eliciting nucleic acids condensation processes into tight, highly ordered tertiary structures as well as altering the extent of order and compactness within the resulting species. Specifically, such factors include the ionic strength, the presence of particular dehydrating agents and polyamines, as well as the pH values. It is proposed that slight alterations of these parameters induce the formation of short non-B-DNA segments that propagate as a perturbation along the B-DNA double helix. The structural fluctuations of the dsDNA molecules that result from the conformational discontinuities formed at the junction sites between the B motif and the conformationally altered segments alter the elastic response of the nucleic acids and facilitate cooperative condensation processes. Moreover, the type and frequency of the structurally modified clusters interspersed within the B conformation and determined by the environmental parameters are shown to provide a means for continuous regulation of the extent and mode of DNA packaging. The ionic strength and hydrophobic environment in the close vicinity of the DNA molecules are controlled and modulated in vivo by DNA-binding proteins such as histones and protamines; similarly, pH values and polyamine concentrations are constantly regulated in living systems. It is suggested, therefore, that the secondary structural polymorphism which characterizes the DNA molecules might display a regulatory role by acting as a functional link between cellular parameters and the extent, mode, and timing of nucleic acid packaging processes.     

Interactions Between Amino Groups in DNA. An Ab Initio Study and a Comparison with Empirical Potentials

Abstract

An ab initio quantum chemical analysis of the close amino group contacts, existing in many DNA crystal structures, is presented. The calculations are made at the Hartree-Fock (HF) level with medium 6–31G* and 6–31 G(NH₂*) basis sets as well as with inclusion of correlation energy using the second order Møller-Plesset theory (MP2) with the 6–31G* basis set. We demonstrate that the model system (methylamine dimer, cytosine dimer) amino groups are forced to adopt significantly non-planar geometry to stabilize their mutual interaction. Comparison is made with a representative set of empirical potentials including AMBER, CHARMM and GROMOS. The empirical potentials are not reliable enough to analyze the amino group contacts occurring in the DNA double helices. We propose that the mutual amino group interactions contribute to the conformational variability of the CpG and ApT B-DNA steps.     

Local sequential minimization of double stranded B-DNA using Monte Carlo annealing

A software algorithm has been developed to investigate the folding process in B-DNA structures in vacuum under a simple and accurate force field. This algorithm models linear double stranded B-DNA sequences based on a local, sequential minimization procedure. The original B-DNA structures were modeled using initial nucleotide structures taken from the Brookhaven database. The models contain information at the atomic level allowing one to investigate as accurately as possible the structure and characteristics of the resulting DNA structures. A variety of DNA sequences and sizes were investigated containing coding and non-coding, random and real, homogeneous or heterogeneous sequences in the range of 2 to 40 base pairs. The force field contains terms such as angle bend, Lennard-Jones, electrostatic interactions and hydrogen bonding which are set up using the Dreiding II force field and defined to account for the helical parameters such as twist, tilt and rise. A close comparison was made between this local minimization algorithm and a global one (previously published) in order to find out advantages and disadvantages of the different methods. From the comparison, this algorithm gives better and faster results than the previous method, allowing one to minimize larger DNA segments. DNA segments with a length of 40 bases need approximately 4 h, while 2.5 weeks are needed with the previous method. After each minimization the angles between phosphate–oxygen-carbon A₁, the oxygen–phosphate–oxygen A₂ and the average helical twists were calculated. From the generated fragments it was found that the bond angles are A₁=150°±2°and A₂=130°±10°, while the helical twist is 36.6°±2° in the A strand and A₁=150°±6° and A₂=130±6° with helical twist 39.6°±2° in the B strand for the DNA segment with the same sequence as the Dickerson dodecamer.Figure The final minimized DNA segment of the Dickerson dodecamer sequence represented by ball drawings and viewed (left) perpendicular and (right) down the helical axis     

Molecular Dynamics Studies of Sequence-directed Curvature in Bending Locus of Trypanosome Kinetoplast DNA

Abstract

The macroscopic curvature induced in the double helical B-DNA by regularly repeated adenine tracts (A-tracts) plays an exceptional role in structural studies of DNA because this effect presents the most well documented example of sequence specific conformational modulations. Recently, a new hypothesis of its physical origin has been put forward. According to it, the intrinsic bends in B-DNA may represent one of the consequences of the compressed frustrated state of its backbone. The compressed backbone hypothesis agrees with many data and explains some controversial experimental observations. The original arguments of this theory came out from MD simulations of a DNA fragment with a strong bending propensity. Its sequence, however, was not experimental. It was constructed empirically so as to maximize the magnitude of bending in calculations. To make sure that our computations reproduce the experimental effect we carried out similar simulations with an A-tract repeat of a natural base pair sequence found in a bent locus of a minicircle DNA. We demonstrate spontaneous development of static curvature in the course of MD simulations excluding any initial bias except the base pair sequence. Its direction and magnitude agree with experimental estimates. The results confirm earlier qualitative conclusions and agree with the hypothesis of a compressed backbone as the origin of static bending in B-DNA.     

Selective Binding of Synthetic Polypeptides to DNA of Varying Composition and Sequence: Effect of Minor Groove Binding Drugs

Abstract

Repetitive basic polypeptides containing lysine or arginine as every third amino acid were shown to cause DNA condensation at physiological salt concentration connected with selective DNA binding with respect to DNA composition and sequence. This selectivity is very similar to that existing in the case of histone H1 and other basic proteins and does not depend on polypeptide chain conformation. The effect of the minor groove binding drugs netropsin and distamycin was tested to elucidate the origin of the binding selectivity. The results suggest that the binding preferences are due to the variations in the conformation in various types of B-DNA that depend on DNA composition and sequence. The most important factor affecting the selectivity is probably the value of the negative electrostatic potential in the minor groove.     

Some characteristics of crossing over in induced recombination between chromosomes of wheat and rye

Allopolyploid wheat (Triticum aestivum L.) carries three pairs of homoeologous genomes but its meiotic pairing is diploid-like. This is the effect of the Ph (pairing homoeologous) system which restricts chromosome pairing to strictly homologous. Ph1 is the locus with the strongest effect. Disabling Ph1 permits pairing between homoeologues and is routinely used in chromosome engineering to introgress alien variation into breeding stocks. Whereas the efficiency of Ph1 and the general pattern of homoeologous crossovers in its absence are quite well known from numerous studies, other characteristics of such crossovers remain unknown. This study analyzed the crossover points in four sets of the ph1b-induced recombinants between wheat homologues as well as between three wheat and rye (Secale cereale) homoeologous chromosome arms, and compared them to crossovers between homologues in a reference wheat population. The results show the Ph1 locus also controls crossing over of homologues, and the general patterns of homologous (with Ph1) and homoeologous (with ph1b) crossing over are the same. In all intervals analyzed, homoeologous crossovers fell within the range of frequency distribution of homologous crossovers among individual families of the reference population. No specific DNA sequence characteristics were identified that could be recognized by the Ph1 locus; the only difference between homologous and homoeologous crossing over appears to be in frequency. It is concluded that the Ph1 locus likely recognizes DNA sequence similarity; crossing over is permitted between very similar sequences. In the absence of Ph1 dissimilarities are ignored, in proportion to the level of the sequence divergence.     

Unusual DNA conformation at low pH revealed by NMR: parallel-stranded DNA duplex with homo base pairs.

We have investigated the conformational potentials of several DNA oligonucleotides containing sequences related to 5'-CGA in neutral pH and low pH (< 5.0) conditions. One-dimensional proton NMR spectra show that d(CGATCG), d(TCGATCGA), and d(CGATCGATCG) exhibit new sets of resonances at low pH (approximately 3.8-4.4), when compared to those from the neutral pH samples. The low pH form and the neutral pH form are in slow equilibrium. Analyses of the data suggest that these sequences under low pH conditions adopt structures distinct from B-DNA. Two-dimensional nuclear Overhauser effect spectroscopy (2D NOESY) data from the DNA hexamer d(CGATCG) of the neutral and low pH samples were used to analyze their respective structures in solution. An iterative NOE spectral-driven refinement procedure, SPEDREF [Robinson, H., & Wang, A. H.-J. (1992) Biochemistry 31, 3524-3533], was used to show that the neutral pH structure is close to canonical B-DNA. In contrast, analysis of the low pH form using the 2D NOESY data suggests that its structure is consistent with a right-handed parallel-stranded (PS) double helix with symmetrical non-Watson-Crick (C+:C, G:G, A:A, T:T) homo base pairs. Supporting evidence for the PS helix includes the asymmetric inversion-recovery relaxation times associated with the two ends of the helix. The structure is favored by the 5'-CGA sequence in which the cytosines provide the C+:C pairing for the nucleation step and the GpA step is significantly stabilized by the interstrand G-A stacking interactions.(ABSTRACT TRUNCATED AT 250 WORDS)     

Molecular dynamic simulations of environment and sequence dependent DNA conformations: the development of the BMS nucleic acid force field and comparison with experimental results

Molecular dynamic (MD) simulations using the BMS nucleic acid force field produce environment and sequence dependent DNA conformations that closely mimic experimentally derived structures. The parameters were initially developed to reproduce the potential energy surface, as defined by quantum mechanics, for a set of small molecules that can be used as the building blocks for nucleic acid macromolecules (dimethyl phosphate, cyclopentane, tetrahydrofuran, etc.). Then the dihedral parameters were fine tuned using a series of condensed phase MD simulations of DNA and RNA (in zero added salt, 4M NaCl, and 75% ethanol solutions). In the tuning process the free energy surface for each dihedral was derived from the MD ensemble and fitted to the conformational distributions and populations observed in 87 A- and B-DNA x-ray and 17 B-DNA NMR structures. Over 41 nanoseconds of MD simulations are presented which demonstrate that the force field is capable of producing stable trajectories, in the correct environments, of A-DNA, double stranded A-form RNA, B-DNA, Z-DNA, and a netropsin-DNA complex that closely reproduce the experimentally determined and/or canonical DNA conformations. Frequently the MD averaged structure is closer to the experimentally determined structure than to the canonical DNA conformation. MD simulations of A- to B- and B- to A-DNA transitions are also shown. A-DNA simulations in a low salt environment cleanly convert into the B-DNA conformation and converge into the RMS space sampled by a low salt simulation of the same sequence starting from B-DNA. In MD simulations using the BMS force field the B-form of d(GGGCCC)2 in a 75% ethanol solution converts into the A-form. Using the same methodology, parameters, and conditions the A-form of d(AAATTT)2 correctly converts into the B-DNA conformation. These studies demonstrate that the force field is capable of reproducing both environment and sequence dependent DNA structures. The 41 nanoseconds (nsec) of MD simulations presented in this paper paint a global picture which suggests that the DNA structures observed in low salt solutions are largely due to the favorable internal energy brought about by the nearly uniform screening of the DNA electrostatics. While the conformations sampled in high salt or mixed solvent environments occur from selective and asymmetric screening of the phosphate groups and DNA grooves, respectively, brought about by sequence induced ion and solvent packing.     

First evidence on induced topological changes in supercoiled DNA by an aluminium <Emphasis Type="SmallCaps">d</Emphasis>-aspartate complex

Aluminium is a debatable and suspected etiological factor in neurodegenerative disorders. Aluminium–amino acid complexes also play an important role in the complex biology of the metal. Recent reports indicate the presence of d-aspartate and d-glutamate in aging brain, human breast tumors, core amyloid plaques and neurofibrillary tangles of Alzheimer's brain. This stereoinversion from the l- to the d-enantiomer is enhanced by Al. Further, the observation that Al is localized in the chromatin region encouraged the present study of the interaction of Al–amino acid complexes with DNA. This study used circular dichroism of supercoiled DNA and showed that Al–d-Asp caused a native B-DNA to C-DNA conformational change, while Al–l-Asp, Al–l-Glu and Al–d-Glu did not alter the native B-DNA conformation. This differential DNA binding property of Al–amino acid complexes is assigned to the stereoisomerism and chirality of the complexes. Interestingly, polyamines like spermine further induced an asymmetric condensation of the "limit C-motif" induced by Al–d-Asp to a -DNA. The results were supported by computer modeling, gel studies and ethidium bromide binding. We also propose a mechanism of Al–d-Asp binding and its ability to modulate DNA topology.     

Azimuthal frustration and bundling in columnar DNA aggregates

The interaction between two stiff parallel DNA molecules is discussed using linear Debye-Hückel screening theory with and without inclusion of the dielectric discontinuity at the DNA surface, taking into account the helical symmetry of DNA. The pair potential furthermore includes the amount and distribution of counterions adsorbed on the DNA surface. The interaction does not only depend on the interaxial separation of two DNA molecules, but also on their azimuthal orientation. The optimal mutual azimuthal angle is a function of the DNA-DNA interaxial separation, which leads to azimuthal frustrations in an aggregate. On the basis of the pair potential, the positional and orientational order in columnar B-DNA assemblies in solution is investigated. Phase diagrams are calculated using lattice sums supplemented with the entropic contributions of the counterions in solution. A variety of positionally and azimuthally ordered phases and bundling transitions is predicted, which strongly depend on the counterion adsorption patterns.     

DNA structure and polymerase fidelity.

The accuracy of DNA replication results from both the intrinsic DNA polymerase fidelity and the DNA sequence. Although the recent structural studies on polymerases have brought new insights on polymerase fidelity, the role of DNA sequence and structure is less well understood. Here, the analysis of the crystal structures of hotspots for polymerase slippage including (CA)n and (A)n tracts in different intermolecular contexts reveals that, in the B-form, these sequences share common structural alterations which may explain the high rate of replication errors. In particular, a two-faced "Janus-like" structure with shifted base-pairs in the major groove but an apparent normal geometry in the minor groove constitutes a molecular decoy specifically suitable to mislead the polymerases. A model of the rat polymerase beta bound to this structure suggests that an altered conformation of the nascent template-primer duplex can interfere with correct nucleotide incorporation by affecting the geometry of the active site and breaking the rules of base-pairing, while at the same time escaping enzymatic mechanisms of error discrimination which scan for the correct geometry of the minor groove.In contrast, by showing that the A-form greatly attenuates the sequence-dependent structural alterations in hotspots, this study suggests that the A-conformation of the nascent template-primer duplex at the vicinity of the polymerase active site will contribute to fidelity. The A-form may play the role of a structural buffer which preserves the correct geometry of the active site for all sequences. The detailed comparison of the conformation of the nascent template-primer duplex in the available crystal structures of DNA polymerase-DNA complexes shows that polymerase beta, the least accurate enzyme, is unique in binding to a B-DNA duplex even close to its active site. This model leads to several predictions which are discussed in the light of published experimental data.     

Geometry of the DNA strands within the RecA nucleofilament: role in homologous recombination

Homologous recombination consists of exchanging DNA strands of identical or almost identical sequence. This process is important for both DNA repair and DNA segregation. In prokaryotes, it involves the formation of long helical filaments of the RecA protein on DNA. These filaments incorporate double-stranded DNA from the cell's genetic material, recognize sequence homology and promote strand exchange between the two DNA segments. DNA processing by these nucleofilaments is characterized by large amplitude deformations of the double helix, which is stretched by 50% and unwound by 40% with respect to B-DNA. In this article, information concerning the structure and interactions of the RecA, DNA and ATP molecules involved in DNA strand exchange is gathered and analyzed to present a view of their possible arrangement within the filament, their behavior during strand exchange and during ATP hydrolysis, the mechanism of RecA-promoted DNA deformation and the role of DNA deformation in the process of homologous recombination. In particular, the unusual characteristics of DNA within the RecA filament are compared to the DNA deformations locally induced by architectural proteins which bind in the DNA minor groove. The possible role and location of two flexible loops of RecA are discussed.     

The zab domain of the human RNA editing enzyme ADAR1 recognizes Z-DNA when surrounded by B-DNA

The Zab domain of the editing enzyme ADAR1 binds tightly and specifically to Z-DNA stabilized by bromination or supercoiling. A stoichiometric amount of protein has been shown to convert a substrate of suitable sequence to the Z form, as demonstrated by a characteristic change in the CD spectrum of the DNA. Now we show that Zab can bind not only to isolated Z-forming d(CG)(n) sequences but also to d(CG)(n) embedded in B-DNA. The binding of Zab to such sequences results in a complex including Z-DNA, B-DNA, and two B-Z junctions. In this complex, the d(CG)(n) sequence, but not the flanking region, is in the Z conformation. The presence of Z-DNA was detected by cleavage with a Z-DNA specific nuclease, by undermethylation using Z-DNA sensitive SssI methylase, and by circular dichroism. It is possible that Zab binds to B-DNA with low affinity and flips any favorable sequence into Z-DNA, resulting in a high affinity complex. Alternatively, Zab may capture Z-DNA that exists transiently in solution. The binding of Zab to potential as well as established Z-DNA segments suggests that the range of biological substrates might be wider than previously thought.     

Sequence specific hydrogen bond of DNA with denaturants affects its stability: Spectroscopic and simulation studies

BackgroundSeveral different small molecules have been used to target the DNA helix in order to treat the diseases caused by its mutation. Guanidinium(Gdm⁺) and urea based drugs have been used for the diseases related to central nervous system, also as the anti-inflammatory and chemotherapeutic agent. However, the role of Gdm⁺ and urea in the stabilization/destabilization of DNA is not well understood.MethodsSpectroscopic techniques along with molecular dynamics (MD) simulation have been performed on different sequences of DNA in the presence of guanidinium chloride (GdmCl) and urea to decode the binding of denaturants with DNA and the role of hydrogen bond with the different regions of DNA in its stability/destability.Results and conclusionOur study reveals that, Gdm⁺ of GdmCl and urea both intrudes into the groove region of DNA along with the interaction with its phosphate backbone. However, interaction of Gdm⁺ and urea with the nucleobases in the groove region is different. Gdm⁺ forms the intra-strand hydrogen bond with the central region of the both sequences of DNA whereas inter-strand hydrogen bond along with water assisted hydrogen bond takes place in the case of urea. The intra-strand hydrogen bond formation capability of Gdm⁺ with the nucleobases in the minor groove of DNA decreases its groove width which probably causes the stabilization of B-DNA in GdmCl. In contrast, the propensity of the formation of inter-strand hydrogen bond of urea with the nucleobases in the groove region of DNA without affecting the groove width destabilizes B-DNA as compared to GdmCl. This study depicts that the opposite effect of GdmCl and urea on the stability is a general property of B-DNA. However, the extent of stabilization/destabilization of DNA in Gdm⁺ and urea depend on its sequence probably due to the difference in the intra/inter-strand hydrogen bonding with different bases present in both the sequences of DNA.General significanceThe information obtained from this study will be useful for the designing of Gdm⁺ based drug molecule which can target the DNA more specifically and selectively.     

Base-pair opening and spermine binding--B-DNA features displayed in the crystal structure of a gal operon fragment: implications for protein-DNA recognition.

A sequence that is represented frequently in functionally important sites involving protein-DNA interactions is GTG/CAC, suggesting that the trimer may play a role in regulatory processes. The 2.5 A resolution structure of d(CGGTGG)/d(CCACCG), a part of the interior operator (OI, nucleotides +44 to +49) of the gal operon, co-crystallized with spermine, is described herein. The crystal packing arrangement in this structure is unprecedented in a crystal of B-DNA, revealing a close packing of columns of stacked DNA resembling a 5-stranded twisted wire cable. The final structure contains one hexamer duplex, 17 water molecules and 1.5 spermine molecules per crystallographic asymmetric unit. The hexamer exhibits base-pair opening and shearing at T.A resulting in a novel non-Watson-Crick hydrogen-bonding scheme between adenine and thymine in the GTG region. The ability of this sequence to adopt unusual conformations in its GTG region may be a critical factor conferring sequence selectivity on the binding of Gal repressor. In addition, this is the first conclusive example of a crystal structure of spermine with native B-DNA, providing insight into the mechanics of polyamine-DNA binding, as well as possible explanations for the biological action of spermine.     

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.