The contribution of transient counterion imbalances to DNA bending fluctuations

A two-sided model for DNA is employed to analyze fluctuations of the spatial distribution of condensed counterions and the effect of these fluctuations on transient bending. We analyze two classes of fluctuations. In the first, the number of condensed counterions on one side of the DNA remains at its average value, while on the other side, counterions are lost to bulk solution or gained from it. The second class of fluctuations is characterized by movement of some counterions from one side of the DNA to the other. The root-mean-square fluctuation for each class is calculated from counterion condensation theory. The amplitude of the root-mean-square fluctuation depends on the ionic strength as well as the length of the segment considered and is of the order 5-10%. Both classes of fluctuation result in transient bends toward the side of greater counterion density. The bending amplitudes are approximately 15% of the total root-mean-square bends associated with the persistence length of DNA. We are thus led to suggest that asymmetric fluctuations of counterion density contribute modestly but significantly toward the aggregate of thermalized solvent fluctuations that cause bending deformations of DNA free in solution. The calculations support the idea that counterions may exert some modulating influence on the fine structure of DNA.     

Electrostatic contribution to the energy of binding of aromatic ligands with DNA

The method of solution of the nonlinear Poisson-Boltzmann equation was used to calculate electrostatic energy of binding of various aromatic ligands with DNA oligomers of different length. Analysis of the electrostatic contribution was made in terms of a two-step DNA binding process: formation of the intercalation cavity and insertion of the ligand. The total electrostatic energy was also partitioned into components: the energy of atom-atom coulombic interactions and the energy of interaction with surrounding water. The results indicate that electrostatic interactions are, as a whole, unfavorable to the intercalation process and that a correct analysis of structure-energy interrelation for Ligand-DNA interactions should only be accomplished at the level of the components rather than at the level of total electrostatic energy.     

Electrostatic free energy landscapes for DNA helix bending

Nucleic acids are highly charged polyanionic molecules; thus, the ionic conditions are crucial for nucleic acid structural changes such as bending. We use the tightly bound ion theory, which explicitly accounts for the correlation and ensemble effects for counterions, to calculate the electrostatic free energy landscapes for DNA helix bending. The electrostatic free energy landscapes show that DNA bending energy is strongly dependent on ion concentration, valency, and size. In a Na⁺ solution, DNA bending is electrostatically unfavorable because of the strong charge repulsion on backbone. With the increase of the Na⁺ concentration, the electrostatic bending repulsion is reduced and thus the bending becomes less unfavorable. In contrast, in an Mg²⁺ solution, ion correlation induces a possible attractive force between the different parts of the helical strands, resulting in bending. The electrostatically most favorable and unfavorable bending directions are toward the major and minor grooves, respectively. Decreasing the size of the divalent ions enhances the electrostatic bending attraction, causing an increased bending angle, and shifts the most favorable bending to the direction toward the minor groove. The microscopic analysis on ion-binding distribution reveals that the divalent ion-induced helix bending attraction may come from the correlated distribution of the ions across the grooves in the bending direction.     

The contribution of phosphate-phosphate repulsions to the free energy of DNA bending

DNA bending is important for the packaging of genetic material, regulation of gene expression and interaction of nucleic acids with proteins. Consequently, it is of considerable interest to quantify the energetic factors that must be overcome to induce bending of DNA, such as base stacking and phosphate–phosphate repulsions. In the present work, the electrostatic contribution of phosphate–phosphate repulsions to the free energy of bending DNA is examined for 71 bp linear and bent-form model structures. The bent DNA model was based on the crystallographic structure of a full turn of DNA in a nucleosome core particle. A Green's function approach based on a linear-scaling smooth conductor-like screening model was applied to ascertain the contribution of individual phosphate–phosphate repulsions and overall electrostatic stabilization in aqueous solution. The effect of charge neutralization by site-bound ions was considered using Monte Carlo simulation to characterize the distribution of ion occupations and contribution of phosphate repulsions to the free energy of bending as a function of counterion load. The calculations predict that the phosphate–phosphate repulsions account for ~30% of the total free energy required to bend DNA from canonical linear B-form into the conformation found in the nucleosome core particle.     

Electrostatic contribution to the binding stability of protein-protein complexes

To investigate roles of electrostatic interactions in protein binding stability, electrostatic calculations were carried out on a set of 64 mutations over six protein-protein complexes. These mutations alter polar interactions across the interface and were selected for putative dominance of electrostatic contributions to the binding stability. Three protocols of implementing the Poisson-Boltzmann model were tested. In vdW4 the dielectric boundary between the protein low dielectric and the solvent high dielectric is defined as the protein van der Waals surface and the protein dielectric constant is set to 4. In SE4 and SE20, the dielectric boundary is defined as the surface of the protein interior inaccessible to a 1.4-A solvent probe, and the protein dielectric constant is set to 4 and 20, respectively. In line with earlier studies on the barnase-barstar complex, the vdW4 results on the large set of mutations showed the closest agreement with experimental data. The agreement between vdW4 and experiment supports the contention of dominant electrostatic contributions for the mutations, but their differences also suggest van der Waals and hydrophobic contributions. The results presented here will serve as a guide for future refinement in electrostatic calculation and inclusion of nonelectrostatic effects.     

The energetic contribution of induced electrostatic asymmetry to DNA bending by a site-specific protein

DNA bending can be promoted by reducing the net negative electrostatic potential around phosphates on one face of the DNA, such that electrostatic repulsion among phosphates on the opposite face drives bending toward the less negative surface. To provide the first assessment of energetic contribution to DNA bending when electrostatic asymmetry is induced by a site-specific DNA binding protein, we manipulated the electrostatics in the EcoRV endonuclease-DNA complex by mutation of cationic side chains that contact DNA phosphates and/or by replacement of a selected phosphate in each strand with uncharged methylphosphonate. Reducing the net negative charge at two symmetrically located phosphates on the concave DNA face contributes − 2.3 kcal mol^− 1 to − 0.9 kcal mol^− 1 (depending on position) to complex formation. In contrast, reducing negative charge on the opposing convex face produces a penalty of + 1.3 kcal mol^− 1. Förster resonance energy transfer experiments show that the extent of axial DNA bending (about 50°) is little affected in modified complexes, implying that modification affects the energetic cost but not the extent of DNA bending. Kinetic studies show that the favorable effects of induced electrostatic asymmetry on equilibrium binding derive primarily from a reduced rate of complex dissociation, suggesting stabilization of the specific complex between protein and markedly bent DNA. A smaller increase in the association rate may suggest that the DNA in the initial encounter complex is mildly bent. The data imply that protein-induced electrostatic asymmetry makes a significant contribution to DNA bending but is not itself sufficient to drive full bending in the specific EcoRV-DNA complex.     

Electrostatic contribution of serine phosphorylation to the Drosophila SLBP--histone mRNA complex

Unlike all other metazoan mRNAs, mRNAs encoding the replication-dependent histones are not polyadenylated but end in a unique 26 nucleotide stem-loop structure. The protein that binds the 3' end of histone mRNA, the stem-loop binding protein (SLBP), is essential for histone pre-mRNA processing, mRNA translation, and mRNA degradation. Using biochemical, biophysical, and nuclear magnetic resonance (NMR) experiments, we report the first structural insight into the mechanism of SLBP-RNA recognition. In the absence of RNA, phosphorylated and unphosphorylated forms of the RNA binding and processing domain (RPD) of Drosophila SLBP (dSLBP) possess helical secondary structure but no well-defined tertiary fold. Drosophila SLBP is phosphorylated at four out of five potential serine or threonine sites in the sequence DTAKDSNSDSDSD at the extreme C-terminus, and phosphorylation at these sites is necessary for histone pre-mRNA processing. Here, we provide NMR evidence for serine phosphorylation of the C-terminus using (31)P direct-detect experiments and show that both serine phosphorylation and RNA binding are necessary for proper folding of the RPD. The electrostatic effect of protein phosphorylation can be partially mimicked by a mutant form of SLBP wherein four C-terminal serines are replaced with glutamic acids. Hence, both RNA binding and protein phosphorylation are necessary for stabilization of the SLBP RPD.     

Polyelectrolyte contribution to the persistence length of DNA

The difference between the theories of Manning, on the one hand, and of Odijk and Skolnick and Fixman, on the other, for the polyelectrolyte contribution to the persistence length of DNA is shown to arise entirely from a subtle geometrical error in the theory of Manning. The corrected theory of Manning predicts a negligible polyelectrolyte contribution in 1.0M NaCl and only 33 Å in 0.01M NaCl, thus giving a change in total persistence length by a factor of only 1.07 over that range, in agreement with Odijk. Pertinent data in the literature indicate that the persistence length must change by a factor of ≤ 1.6 between 1.0 and 0.01M NaCl, and very likely by less than a factor of 1.4. Evidently, the intrinsic rigidity of the uncharged double-strand filament dominates the bending rigidity at NaCl concentrations above 0.01M.     

Supercoiling-induced DNA bending

Local DNA bending is a critical factor for numerous DNA functions including recognition of DNA by sequence-specific regulatory binding proteins. Negative DNA supercoiling increases both local and global DNA dynamics, and this dynamic flexibility can facilitate the formation of DNA-protein complexes. We have recently shown that apexes of supercoiled DNA molecules are sites that can promote the formation of an alternative DNA structure, a cruciform, suggesting that these positions in supercoiled DNA are under additional stress and perhaps have a distorted DNA geometry. To test this hypothesis, we used atomic force microscopy to directly measure the curvature of apical positions in supercoiled DNA. The measurements were performed for an inherently curved sequence formed by phased A tracts and a region of mixed sequence DNA. For this, we used plasmids in which an inverted repeat and A tract were placed at precise locations relative to each other. Under specific conditions, the inverted repeat formed a cruciform that was used as a marker for the unambiguous identification of the A tract location. When the A tract and cruciform were placed diametrically opposite, this yielded predominantly nonbranched plectonemic molecules with an extruded cruciform and A tract localized in the terminal loops. For both the curved A tract and mixed sequence nonbent DNA, their localization to an apex increased the angle of bending compared to that expected for DNA unconstrained in solution. This is consistent with increased helical distortion at an apical bend.     

Mechanisms of DNA bending

Genome packaging and gene regulation require DNA bending. Recent developments in the elucidation of the mechanisms involved in DNA bending include new X-ray structures (most notably that of the mammalian nucleosome) wherein DNA is bent, controversy surrounding interpretation of DNA-bending experiments with basic-leucine zipper proteins, studies of electrostatic effects in DNA bending, and the design of artificial DNA-bending ligands.     

DNA binding and bending to initiate packaging of phage lambda DNA

Physical and genetic studies verify that the DNA binding domain of protein gpNu1 (which initiates packaging of phage lambda DNA) is a winged helix-turn-helix (w HTH) and that gpNu1 dimers bind sites that are brought close through DNA bending.     

DNA bending and its relation to nucleosome positioning

X-ray and solution studies have shown that the conformation of a DNA double helix depends strongly on its base sequence. Here we show that certain sequence-dependent modulations in structure appear to determine the rotational positioning of DNA about the nucleosome. Three different experiments are described. First, a piece of DNA of defined sequence (169 base-pairs long) is closed into a circle, and its structure examined by digestion with DNAase I: the helix adopts a highly preferred configuration, with short runs of (A, T) facing in and runs of (G, C) facing out. Secondly, the same sequence is reconstituted with a histone octamer: the angular orientation around the histone core remains conserved, apart from a small uniform increase in helix twist. Finally, it is shown that the average sequence content of DNA molecules isolated from chicken nucleosome cores is non-random, as in a reconstituted nucleosome: short runs of (A, T) are preferentially positioned with minor grooves facing in, while runs of (G, C) tend to have their minor grooves facing out. The periodicity of this modulation in sequence content (10.17 base-pairs) corresponds to the helix twist in a local frame of reference (a result that bears on the change in linking number upon nucleosome formation). The determinants of translational positioning have not been identified, but one possibility is that long runs of homopolymer (dA) X (dT) or (dG) X (dC) will be excluded from the central region of the supercoil on account of their resistance to curvature.     

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.     

DNA bending by the bulge defect

Comparative gel electrophoresis measurements were used to characterize DNA bending in molecules containing an extra adenosine on one strand, the so-called bulge defect. We used oligomers containing A6 tracts separated from the bulged base by varying numbers of nucleotides to determine the direction and magnitude of the bulge bend. Helix unwinding by the bulge was determined from the electrophoretic anomaly as a function of the size of the repeated monomers. We conclude that the bulge bend is 21 degrees +/- 3 degrees, primarily in the direction of tilt away from the bulged base. The total helical advance of the DNA at the bulge site is smaller than would be the case if the complementary T were present, corresponding to an unwinding by 25 degrees +/- 6 degrees. These values are in good agreement with the results of NMR and energy minimization studies of the bulged base in double-helical deoxyoligonucleotides [Woodson, S. A., & Crothers, D.M. (1988) Biochemistry 27, 3130-3141]