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

Force-induced unzipping of DNA in the presence of solvent molecules

The melting of double-stranded DNA (dsDNA) in the presence of solvent molecules is a fundamental process with significant implications for understanding the thermal and mechanical behavior of DNA and its interactions with the surrounding environment. The solvents play an essential role in the structural transformation of DNA subjected to a pulling force. In this study, we simulate the thermal and force induced denaturation of dsDNA and elucidate the solvent dependent melting behavior, identifying key factors that influence the stability of DNA melting in presence of solvent molecules. Using a statistical model, we first find the melting profile of short heterogeneous DNA molecules in the presence of solvent molecules in Force ensemble. We also investigate the effect of solvent's strengths on the melting profile of DNA. In the force ensemble, we consider two homogeneous DNA chains and apply the force on different locations along the chain in the presence of solvent molecules. Different pathways manifest the melting of the molecule in both ensembles, and we found several interesting features of melting DNA in a constant force ensemble, such as lower critical force when the chain is pulled from the base pair close to a solvent molecule. The results provide new insights into the force-induced unzipping of DNA and could be used to develop new methods for controlling the unzipping process. By providing a better understanding of melting and unzipping of dsDNA in the presence of solvent molecules, this study provides valuable guidelines for predicting DNA thermodynamic quantities and for designing DNA nanostructures.

Ceftriaxone induces glial EAAT-2 promotor region via NF-kB conformational changes: An interaction analysis using HADDOCK

Excitotoxicity, depletion of energy metabolites, and ionic imbalance are the major factors involved in neurodegeneration mediated through excitatory amino acid transporter-2 (EAAT-2) dysfunction in ischemic insult. Recent studies have revealed that ceftriaxone expresses EAAT-2 via nuclear transcription factor kappa-B (NF-kB) signaling pathway, stimulation of EAAT-2 expression in the ischemic, and excitotoxic conditions that could provide potential benefits to control neurodegeneration. In this study, we have predicted the in silico model for interaction between NF-kB and EAAT-2 promoter region to rule out the conformational changes for the expression of EAAT-2 protein. Using homology-built model of NF-kB, we identified ceftriaxone-induced conformational changes in gene locus -272 of DNA where NF-kB binding with EAAT-2 promoter region through protein-DNA docking calculation. The interaction profile and conformational dynamics occurred between ceftriaxone predocked and postdocked conformations of NF-kB with DNA employing HADDOCK 2.2 web server followed by 250 ns long all atom explicit solvent molecular dynamics simulations. Both the protein and DNA exhibited modest conformational changes with respect to HADDOCK score, energy terms (desolvation energy [E_desolv ]), van der waal energy (E_vdw ), electrostatic energy (E_elec ), restraints energy (E_air ), buried surface area, root mean square deviation, RMSF, radius of gyration, total hydrogen bonds when ceftriaxone pre- and postdocked NF-kB conformations were bound to DNA. Hence, the conformational changes in the C-terminal domain could be the reason for EAAT-2 expression through ceftriaxone specific binding pocket of -272 of DNA.

Semi-empirical and Linear-Scaling DFT Methods to Characterize duplex DNA and G-quadruplexes in Presence of Interacting Small Molecules

The computational study of DNA and its interaction with ligands is a highly relevant area of research, with significant consequences for developing new therapeutic strategies. However, the computational description of...

Mechanical Flexibility of DNA: A Quintessential Tool for DNA Nanotechnology

The mechanical properties of DNA have enabled it to be a structural and sensory element in many nanotechnology applications. While specific base-pairing interactions and secondary structure formation have been the most widely utilized mechanism in designing DNA nanodevices and biosensors, the intrinsic mechanical rigidity and flexibility are often overlooked. In this article, we will discuss the biochemical and biophysical origin of double-stranded DNA rigidity and how environmental and intrinsic factors such as salt, temperature, sequence, and small molecules influence it. We will then take a critical look at three areas of applications of DNA bending rigidity. First, we will discuss how DNA's bending rigidity has been utilized to create molecular springs that regulate the activities of biomolecules and cellular processes. Second, we will discuss how the nanomechanical response induced by DNA rigidity has been used to create conformational changes as sensors for molecular force, pH, metal ions, small molecules, and protein interactions. Lastly, we will discuss how DNA's rigidity enabled its application in creating DNA-based nanostructures from DNA origami to nanomachines.

DNA size in confined environments

For short DNA molecules in crowded environments, we evaluate macroscopic parameters such as the average end-to-end distance and the twist conformation by tuning the strength of the site specific confinement driven by the crowders.

Understanding the Elementary Steps in DNA Tile-Based Self-Assembly

Although many models have been developed to guide the design and implementation of DNA tile-based self-assembly systems with increasing complexity, the fundamental assumptions of the models have not been thoroughly tested. To expand the quantitative understanding of DNA tile-based self-assembly and to test the fundamental assumptions of self-assembly models, we investigated DNA tile attachment to preformed "multi-tile" arrays in real time and obtained the thermodynamic and kinetic parameters of single tile attachment in various sticky end association scenarios. With more sticky ends, tile attachment becomes more thermostable with an approximately linear decrease in the free energy change (more negative). The total binding free energy of sticky ends is partially compromised by a sequence-independent energy penalty when tile attachment forms a constrained configuration: "loop". The minimal loop is a 2 × 2 tetramer (Loop4). The energy penalty of loops of 4, 6, and 8 tiles was analyzed with the independent loop model assuming no interloop tension, which is generalizable to arbitrary tile configurations. More sticky ends also contribute to a faster on-rate under isothermal conditions when nucleation is the rate-limiting step. Incorrect sticky end contributes to neither the thermostability nor the kinetics. The thermodynamic and kinetic parameters of DNA tile attachment elucidated here will contribute to the future improvement and optimization of tile assembly modeling, precise control of experimental conditions, and structural design for error-free self-assembly.

Investigating the effect of charged amino acids on DNA conformation in EcoRI-DNA complex: a molecular dynamics simulation study

Sequence-specific binding of proteins to DNA is essential for almost all the cellular processes like transcription, translation, replication, etc. One among the various mechanisms that has been identified so far that contributes to the specificity in protein-DNA interaction is the DNA conformational change. Electrostatic neutralization of the phosphate groups by the positively charged amino acid residues in proteins is thought to bring about such conformational changes in DNA. Here, we employ molecular dynamics simulations to examine the effect of charge on amino acids Lys113, Arg145, and Asp91 which are attached to the scissile phosphate on the conformation of DNA in EcoRI-DNA complex. The results indicate that the charge of these amino acids is essential for maintaining the local conformation of DNA in the EcoRI-bound form. Interestingly, we observe that the positively charged amino acids Lys113 and Arg145 have a long-range influence on the DNA conformation, whereas the negatively charged amino acid Asp91 has only a localized effect on the DNA conformation. The charge on the amino acids also alters the collective dynamics of EcoRI. Collectively, the results shed light on the diversity of the effect of charges on DNA conformation as well as on protein dynamics.

Exploring LacI-DNA dynamics by multiscale simulations using the SIRAH force field

The lac repressor protein (LacI) together with its target regulatory sequence are a common model for studying DNA looping and its implications on transcriptional control in bacteria. Owing to the molecular size of this system, standard all-atom (AA) simulations are prohibitive for achieving relevant biological time scales. As an alternative, multiscale models, which combine AA descriptions at particular regions with coarse-grained (CG) representations of the remaining components, were used to address this computational challenge while preserving the relevant details of the system. In this work, we implement a new multiscale approach based on the SIRAH force field to gain deeper insights into the dynamics of the LacI-DNA system. Our methodology allows for a dual resolution treatment of the solute and solvent, explicitly representing the protein, DNA, and solvent environment without compromising the AA region. Starting from the P1 loop configuration in an undertwisted conformation, we were able to observe the transition to the more stable overtwisted state. Additionally, a detailed characterization of the conformational space sampled by the DNA loop was done. In agreement with experimental and theoretical evidence, we observed the transient formation of kinks at the loop, which were stabilized by the presence of counterions at the minor groove. We also show that the loop's intrinsic flexibility can account for reported FRET measurements and bent conformations required to bind the CAP transcription factor.

Ultrafast dynamics of solvation and charge transfer in a DNA-based biomaterial

Charge migration along DNA molecules is a key factor for DNA-based devices in optoelectronics and biotechnology. The association of a significant amount of water molecules in DNA-based materials for the intactness of the DNA structure and their dynamic role in the charge-transfer (CT) dynamics is less documented in contemporary literature. In the present study, we have used a genomic DNA-cetyltrimethyl ammonium chloride (CTMA) complex, a technological important biomaterial, and Hoechest 33258 (H258), a well-known DNA minor groove binder, as fluorogenic probe for the dynamic solvation studies. The CT dynamics of CdSe/ZnS quantum dots (QDs; 5.2 nm) embedded in the as-prepared and swollen biomaterial have also been studied and correlated with that of the timescale of solvation. We have extended our studies on the temperature-dependent CT dynamics of QDs in a nanoenvironment of an anionic, sodium bis(2-ethylhexyl)sulfosuccinate reverse micelle (AOT RMs), whereby the number of water molecules and their dynamics can be tuned in a controlled manner. A direct correlation of the dynamics of solvation and that of the CT in the nanoenvironments clearly suggests that the hydration barrier within the Arrhenius framework essentially dictates the charge-transfer dynamics.

Mechanical properties of DNA-like polymers

The molecular structure of the DNA double helix has been known for 60 years, but we remain surprisingly ignorant of the balance of forces that determine its mechanical properties. The DNA double helix is among the stiffest of all biopolymers, but neither theory nor experiment has provided a coherent understanding of the relative roles of attractive base stacking forces and repulsive electrostatic forces creating this stiffness. To gain insight, we have created a family of double-helical DNA-like polymers where one of the four normal bases is replaced with various cationic, anionic or neutral analogs. We apply DNA ligase-catalyzed cyclization kinetics experiments to measure the bending and twisting flexibilities of these polymers under low salt conditions. Interestingly, we show that these modifications alter DNA bending stiffness by only 20%, but have much stronger (5-fold) effects on twist flexibility. We suggest that rather than modifying DNA stiffness through a mechanism easily interpretable as electrostatic, the more dominant effect of neutral and charged base modifications is their ability to drive transitions to helical conformations different from canonical B-form DNA.

Binding and cleavage of DNA with the restriction enzyme EcoR1 using time-resolved second harmonic generation

The binding of EcoR1 to a 90-bp DNA duplex attached to colloidal microparticles and the subsequent cleavage by the enzyme was observed in real time and label-free with time-resolved second harmonic (SH) spectroscopy. This method provides a unique way to investigate biomolecular interactions based on its sensitivity to changes in structure and electrical charge on formation of a complex and subsequent dynamics. The binding of EcoR1 to the recognition sequence in DNA appears as a rapid increase in the SH signal, which is attributed to the enzyme-induced change in the DNA conformation, going from a rod-like to a bent shape. In the presence of the cofactor Mg(2+), the subsequent decay in the SH signal was monitored in real time as the following processes occurred: cleavage of DNA, dissociation of the enzyme from the DNA, and diffusion of the 74-bp fragment into the bulk solution leaving the 16-bp fragment attached to the microparticle. The observed decay was dependent on the concentration of Mg(2+), which functions as a cofactor and as an electrolyte. With SH spectroscopy the rehybridization dynamics between the rehybridized microparticle bound and free cleaved DNA fragments was observed in real time and label-free following the cleavage of DNA. Collectively, the experiments reported here establish SH spectroscopy as a powerful method to investigate equilibrium and time-dependent biological processes in a noninvasive and label-free way.

Contribution of nonlocal interactions to DNA elasticity

A nonlocal harmonic elastic rod model is proposed to describe the elastic behavior of short DNA molecules. We show that the nonlocal interactions contribute to effective bending energy of the molecule and affect its apparent persistence length. It is also shown that the anomalous behavior which has been observed in all-atom molecular dynamic simulations [A. K. Mazur, Biophys. J. 134, 4507 (2006)] can be a consequence of both nonlocal interactions between DNA base pairs and the intrinsic curvature of DNA.

DNA curvature and flexibility in vitro and in vivo

It has been more than 50 years since the elucidation of the structure of double-helical DNA. Despite active research and progress in DNA biology and biochemistry, much remains to be learned in the field of DNA biophysics. Predicting the sequence-dependent curvature and flexibility of DNA is difficult. Applicability of the conventional worm-like chain polymer model of DNA has been challenged. The fundamental forces responsible for the remarkable resistance of DNA to bending and twisting remain controversial. The apparent 'softening' of DNA measured in vivo in the presence of kinking proteins and superhelical strain is incompletely understood. New methods and insights are being applied to these problems. This review places current work on DNA biophysics in historical context and illustrates the ongoing interplay between theory and experiment in this exciting field.

DNA on a tube: electrostatic contribution to stiffness

Two simple models are used to estimate the electrostatic contributions to the stiffness of short DNA fragments. The first model views DNA as two strands that are appropriately parametrized and are wrapped helically around a straight cylinder radius equal to the radius of the DNA molecule. The potential energy of the DNA due to phosphate-phosphate electrostatic interactions is evaluated assuming that the charges interact through Debye-Hückel potentials. This potential energy is compared with the potential energy as computed using our second model in which DNA is viewed as two helical strands wrapping around a curved tube whose cross-section is a disk of radius equal to the radius of the DNA. We find that the electrostatic persistence length for B-DNA molecules in the range of 105-130 bp is 125.64 angstroms (37 bp) and 76.05 angstroms (23 bp) at 5 and 10 mM monovalent salt concentration, respectively. If the condensed fraction theta is taken to be 0.715 at 10 mM, then the electrostatic persistence length is 108.28 angstroms (32 bp), while that based on taking into account end effects is 72.87 angstroms (21 bp). At 5 mM monovalent salt, the total persistence length for DNA fragments in this length range is approximately 575.64 angstroms (171 bp), using the best estimate for nonelectrostatic contribution to persistence length. Electrostatic effects thus contribute 21.8% to DNA stiffness at 5 mM for fragments between 105- to 130-bp. In contrast, electrostatics are calculated to make a negligible contribution to the DNA persistence length at physiological monovalent cation concentration. The results are compared with counterion condensation models and experimental data.

p53-Induced DNA bending: the interplay between p53-DNA and p53-p53 interactions

Specific p53 binding-induced DNA bending and its underlying responsible forces are crucial for the understanding of selective transcription activation. Diverse p53-response elements exist in the genome; however, it is not known what determines the DNA bending and to what extent. In order to gain knowledge of the forces that govern the DNA bending, molecular dynamics simulations were performed on a series of p53 core domain tetramer-DNA complexes in which each p53 core domain was bound to a DNA quarter site specifically. By varying the sequence of the central 4-base pairs of each half-site, different DNA bending extents were observed. The analysis showed that the dimer-dimer interactions in p53 were similar for the complexes; on the other hand, the specific interactions between the p53 and DNA, including the interactions of Arg280, Lys120, and Arg248 with the DNA, varied more significantly. In particular, the Arg280 interactions were better maintained in the complex with the CATG-containing DNA sequence and were mostly lost in the complex with the CTAG-containing DNA sequence. Structural analysis shows that the base pairings for the CATG sequence were stable throughout the simulation trajectory, whereas those for the CTAG sequence were partially dissociated in part of the trajectory, which affected the stability of the nearby Arg280-Gua base interactions. Thus, DNA bending depends on the balance between the p53 dimer-dimer interactions and p53-DNA interactions, which is in turn related to the DNA sequence and DNA flexibility.

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(2+) 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.

Interplay of ion binding and attraction in DNA condensed by multivalent cations

We have measured forces generated by multivalent cation-induced DNA condensation using single-molecule magnetic tweezers. In the presence of cobalt hexammine, spermidine, or spermine, stretched DNA exhibits an abrupt configurational change from extended to condensed. This occurs at a well-defined condensation force that is nearly equal to the condensation free energy per unit length. The multivalent cation concentration dependence for this condensation force gives the apparent number of multivalent cations that bind DNA upon condensation. The measurements show that the lower critical concentration for cobalt hexammine as compared to spermidine is due to a difference in ion binding, not a difference in the electrostatic energy of the condensed state as previously thought. We also show that the resolubilization of condensed DNA can be described using a traditional Manning–Oosawa cation adsorption model, provided that cation–anion pairing at high electrolyte concentrations is taken into account. Neither overcharging nor significant alterations in the condensed state are required to describe the resolubilization of condensed DNA. The same model also describes the spermidine³⁺/Na⁺ phase diagram measured previously.

Monovalent cation binding by curved DNA molecules containing variable numbers of a-tracts

Monovalent cation binding by DNA A-tracts, runs of four or more contiguous adenine or thymine residues, has been determined for two curved approximately 200 basepair (bp) restriction fragments, one taken from the M13 origin of replication and the other from the VP1 gene of SV40. These two fragments have previously been shown to contain stable, centrally located bends of 44 degrees and 46 degrees , respectively, located within approximately 60 bp "curvature modules" containing four or five irregularly spaced A-tracts. Transient electric birefringence measurements of these two fragments, sequence variants containing reduced numbers of A-tracts in the SV40 curvature module or changes in the residues flanking the A-tracts in the M13 curvature module, have been combined with the free solution electrophoretic mobilities of the same fragments using known equations to estimate the effective charge of each fragment. The effective charge is reduced, on average, by one-third charge for each A-tract in the curvature module, suggesting that each A-tract binds a monovalent cation approximately one-third of the time. Monovalent cation binding to two or more A-tracts is required to observe significant curvature of the DNA helix axis.

DNA bending induced by carbocyclic sugar analogs constrained to the north conformation

DNA bending caused by introduction of carbocyclic sugars constrained to the north conformation was studied, using explicit solvent molecular dynamic (MD) simulations. The native Drew-Dickerson (DD) dodecamer and its three modifications containing north carbocyclic sugars in the 7th (T7*), 8th (T8*) or both 7th and 8th (T7T8*) nucleotide positions were examined. Introduction of the carbocyclic sugar results in A-form conformations for the alpha, beta, chi, zeta, and sugar pucker backbone parameters in the modified nucleotides. Increased steric repulsion between the sugar and its parent base in the modified oligonucleotides impacts the roll and cup dinucleotide step parameters, increasing the bending of the oligomer axis. Increased buckling of the substituted nucleotides disrupts the usual stabilizing base stacking interactions. The level of overall bending depends on the number and position of carbocyclic sugars introduced in the DNA sequence. Single sugar substitutions are unable to induce substantial bending due to the neighboring unmodified nucleotides counterbalancing the distortion. Significant bending can, however, be induced by two consecutive north sugars (T7T8*), which is in agreement with experimental results. The modified oligomers populate a wide range of bend angles, indicating that they maintain flexibility in the bent state. The present results suggest that insertion of carbocyclic sugars into DNA or RNA duplexes can be used to engineer bending of the duplexes without impacting the electrostatic or chemical properties of the phosphodiester backbone, thereby serving as excellent tools for experimental elucidation of nucleic acid structure-function relationships.

Quantitative analysis of monovalent counterion binding to random-sequence, double-stranded DNA using the replacement ion method

A variation of affinity capillary electrophoresis, called the replacement ion (RI) method, has been developed to measure the binding of monovalent cations to random sequence, double-stranded (ds) DNA. In this method, the ionic strength is kept constant by gradually replacing a non-binding ion in the solution with a binding ion and measuring the mobility of binding and non-binding analytes as a function of binding ion concentration. The method was validated by measuring the binding of Li+ ions to adenosine nucleotides; the apparent dissociation constants obtained by the RI method are comparable to literature values obtained by other methods. The binding of Tris+, NH4+, Li+, Na+, and K+ to dsDNA was then investigated. The apparent dissociation constants observed for counterion binding to a random-sequence 26-base pair (bp) oligomer ranged from 71 mM for Tris+ to 173 mM for Na+ and K+. Hence, positively charged Tris buffer ions will compete with other monovalent cations in Tris-buffered solutions. The bound cations identified in this study may correspond to the strongly correlated, tightly bound ions recently postulated to exist as a class of ions near the surface of dsDNA (Tan, Z.-J., and Chen, S.-J. (2006) Biophys. J. 91, 518-536). Monovalent cation binding to random-sequence dsDNA would be expected to occur in addition to any site-specific binding of cations to A-tracts or other DNA sequence motifs. Single-stranded DNA oligomers do not bind the five tested cations under the conditions investigated here.

Computer Simulation of the Ionic Atmosphere around Z-DNA

We describe a coarse-grained model for Z-DNA that mimics the DNA shape with a relatively small number of repulsive interaction sites. In addition, negative charges are placed at the phosphate positions. The ionic atmosphere around this grooved Z-DNA model is then investigated with Monte Carlo simulation. Cylindrically averaged concentration profiles as well as the spatial distribution of ions have been calculated. The results are compared to those for other DNA models differing in the repulsive core. This allows the examination of the effect of the DNA shape in the ionic distribution. It is seen that the penetrability of the ions to the DNA groove plays an important role in the ionic distribution. The results are also compared with those reported for B-DNA. In both conformers the ions are structured in alternating layers of positive and negative charge. In Z-DNA the layers are more or less concentric to the molecular axis. Besides, no coions enter into the single groove of this conformer. On the contrary, the alternating layers of B-DNA are also structured along the axial coordinate with some coions penetrating into the major groove. In both cases we have found five preferred locations of the counterions and two for the coions. The concentration of counterions reaches its absolute maximum at the narrow Z-DNA groove and at the minor groove of B-DNA, the value of the maximum being higher in the Z conformer.

DNA bending by bHLH charge variants

We wish to understand the role of electrostatics in DNA stiffness and bending. The DNA charge collapse model suggests that mutual electrostatic repulsions between neighboring phosphates significantly contribute to DNA stiffness. According to this model, placement of fixed charges near the negatively charged DNA surface should induce bending through asymmetric reduction or enhancement of these inter-phosphate repulsive forces. We have reported previously that charged variants of the elongated basic-leucine zipper (bZIP) domain of Gcn4p bend DNA in a manner consistent with this charge collapse model. To extend this result to a more globular protein, we present an investigation of the dimeric basic-helix–loop–helix (bHLH) domain of Pho4p. The 62 amino acid bHLH domain has been modified to position charged amino acid residues near one face of the DNA double helix. As observed for bZIP charge variants, DNA bending toward appended cations (away from the protein:DNA interface) is observed. However, unlike bZIP proteins, DNA is not bent away from bHLH anionic charges. This finding can be explained by the structure of the more globular bHLH domain which, in contrast to bZIP proteins, makes extensive DNA contacts along the binding face.

Stability of gold nanoparticle-bound DNA toward biological, physical, and chemical agents

Positively charged trimethylammonium-modified mixed monolayer protected clusters (MMPCs) interact with DNA by complementary electrostatic binding, serving as efficient DNA delivery systems. The stability of gold nanoparticle-bound DNA toward biological, physical, and chemical agents is investigated. The MMPC-bound DNA is efficiently protected from DNAse I digestion and experiences nicking/cleavage-induced morphology changes with higher concentrations of DNAse I. Significant protection of MMPC-bound DNA was also observed in a physical sonication assay. However, the MMPC-bound DNA was found to show enhanced cleavage upon exposure to chemically induced radicals. The latter may indicate that bound DNA is bent and wrapped on the surface of the cationic MMPC.

Effect of organic cosolvents on the free solution mobility of curved and normal DNA molecules

The free solution mobilities of curved and normal 199-bp DNA fragments have been measured in buffer solutions containing various quantities of the organic cosolvents methanol, ethanol, 2-propanol, 2-methyl-2,4-pentanediol (MPD), ethylene glycol, and ACN, using CE. The curved fragment, taken from the VP1 gene of SV40, contains five unevenly spaced A- and T-tracts in a centrally located "curvature module"; the A- and T-tracts have been mutated to other sequences in the normal 199-bp fragment. The free solution mobility of the curved 199-bp fragment is significantly lower than that of its normal counterpart in aqueous solutions [Stellwagen, E., Lu, Y. J., Stellwagen, N. C., Nucleic Acids Res. 2005, 33, 4425-4432]. The mobilities of both the curved and normal fragments decrease with increasing cosolvent concentration, due to the effect of the cosolvent on the viscosity and dielectric constant of the solution. The mobility differences between the curved and normal 199-bp fragments and the mobility ratios decrease approximately linearly with the increasing mole fraction of cosolvent in the solution. Hence, MPD and other organic cosolvents affect DNA electrophoretic mobility by a common mechanism, most likely the preferential hydration of the DNA surface that occurs in aqueous cosolvents. The gradual loss of the anomalously slow mobility of the curved 199-bp fragment with increasing cosolvent concentration, combined with other data in the literature, suggests that preferential hydration gradually widens the narrow A-tract minor groove, releasing site-bound counterions in the minor groove and shifting the conformation toward that of normal DNA.

Curved DNA molecules migrate anomalously slowly in free solution

The electrophoretic mobility of a curved DNA restriction fragment taken from the VP1 gene in the SV40 minichromosome has been measured in polyacrylamide gels and free solution, using capillary electrophoresis. The 199 bp restriction fragment has an apparent bend angle of 46 ± 2° located at SV40 sequence position 1922 ± 2 bp [Lu Y.J., Weers B.D. and Stellwagen N. C. (2005) Biophys. J., 88, 1191–1206]. The ‘curvature module’ surrounding the apparent bend center contains five unevenly spaced A- and T-tracts, which are responsible for the observed curvature. The parent 199 bp fragment and sequence mutants containing at least one A-tract in the curvature module migrate anomalously slowly in free solution, as well as in polyacrylamide gels. Hence, the anomalously slow mobilities observed for curved DNA molecules in polyacrylamide gels are due in part to their anomalously slow mobilities in free solution. Analysis of the gel and free solution mobility decrements indicates that each A- or T-tract contributes independently, but not equally, to the curvature of the 199 bp fragment and its A-tract mutants. The relative contribution of each A- or T-tract to the observed curvature depends on its spacing with respect to the first A-tract in the curvature module.