Effects of Na+ on the persistence length and excluded volume of T7 bacteriophage DNA

Evaporation-induced self-assembly of liquid crystal biopolymers

Evaporation-induced self-assembly (EISA) is a process that has gained significant attention in recent years due to its fundamental science and potential applications in materials science and nanotechnology. This technique involves controlled drying of a solution or dispersion of materials, forming structures with specific shapes and sizes. In particular, liquid crystal (LC) biopolymers have emerged as promising candidates for EISA due to their highly ordered structures and biocompatible properties after deposition. This review provides an overview of recent progress in the EISA of LC biopolymers, including DNA, nanocellulose, viruses, and other biopolymers. The underlying self-assembly mechanisms, the effects of different processing conditions, and the potential applications of the resulting structures are discussed.

Dilute polyelectrolyte solutions: recent progress and open questions

Polyelectrolytes are a class of polymers possessing ionic groups on their repeating units. Since counterions can dissociate from the polymer backbone, polyelectrolyte chains are strongly influenced by electrostatic interactions. As a result, the physical properties of polyelectrolyte solutions are significantly different from those of electrically neutral polymers. The aim of this article is to highlight key results and some outstanding questions in the polyelectrolyte research from recent literature. We focus on the influence of electrostatics on conformational and hydrodynamic properties of polyelectrolyte chains. A compilation of experimental results from the literature reveals significant disparities with theoretical predictions. We also discuss a new class of polyelectrolytes called poly(ionic liquid)s that exhibit unique physical properties in comparison to ordinary polyelectrolytes. We conclude this review by listing some key research challenges in order to fully understand the conformation and dynamics of polyelectrolytes in solutions.

Mechanical Constraint Effect on DNA Persistence Length

Persistence length is a significant criterion to characterize the semi-flexibility of DNA molecules. The mechanical constraints applied on DNA chains in new single-molecule experiments play a complex role in measuring DNA persistence length; however, there is a difficulty in quantitatively characterizing the mechanical constraint effects due to their complex interactions with electrostatic repulsions and thermal fluctuations. In this work, the classical buckling theory of Euler beam and Manning's statistical theories of electrostatic force and thermal fluctuation force are combined for an isolated DNA fragment to formulate a quantitative model, which interprets the relationship between DNA persistence length and critical buckling length. Moreover, this relationship is further applied to identify the mechanical constraints in different DNA experiments by fitting the effective length factors of buckled fragments. Then, the mechanical constraint effects on DNA persistence lengths are explored. A good agreement among the results by theoretical models, previous experiments, and present molecular dynamics simulations demonstrates that the new superposition relationship including three constraint-dependent terms can effectively characterize changes in DNA persistence lengths with environmental conditions, and the strong constraint-environment coupling term dominates the significant changes of persistence lengths; via fitting effective length factors, the weakest mechanical constraints on DNAs in bulk experiments and stronger constraints on DNAs in single-molecule experiments are identified, respectively. Moreover, the consideration of DNA buckling provides a new perspective to examine the bendability of short-length DNA.

The Effects of Flexibility on dsDNA–dsDNA Interactions

A detailed understanding of the physical mechanism of ion-mediated dsDNA interactions is important in biological functions such as DNA packaging and homologous pairing. We report the potential of mean force (PMF) or the effective solvent mediated interactions between two parallel identical dsDNAs as a function of interhelical separation in 0.15 M NaCl solution. Here, we study the influence of flexibility of dsDNAs on the effective interactions by comparing PMFs between rigid models and flexible ones. The role of flexibility of dsDNA pairs in their association is elucidated by studying the energetic properties of Na⁺ ions as well as the fluctuations of ions around dsDNAs. The introduction of flexibility of dsDNAs softens the vdW contact wall and induces more counterion fluctuations around dsDNAs. In addition, flexibility facilitates the Na⁺ ions dynamics affecting their distribution. The results quantify the extent of attraction influenced by dsDNA flexibility and further emphasize the importance of non-continuum solvation approaches.

Interconnecting solvent quality, transcription, and chromosome folding in Escherichia coli

All cells fold their genomes, including bacterial cells, where the chromosome is compacted into a domain-organized meshwork called the nucleoid. How compaction and domain organization arise is not fully understood. Here, we describe a method to estimate the average mesh size of the nucleoid in Escherichia coli. Using nucleoid mesh size and DNA concentration estimates, we find that the cytoplasm behaves as a poor solvent for the chromosome when the cell is considered as a simple semidilute polymer solution. Monte Carlo simulations suggest that a poor solvent leads to chromosome compaction and DNA density heterogeneity (i.e., domain formation) at physiological DNA concentration. Fluorescence microscopy reveals that the heterogeneous DNA density negatively correlates with ribosome density within the nucleoid, consistent with cryoelectron tomography data. Drug experiments, together with past observations, suggest the hypothesis that RNAs contribute to the poor solvent effects, connecting chromosome compaction and domain formation to transcription and intracellular organization.

Nucleosome plasticity is a critical element of chromatin liquid-liquid phase separation and multivalent nucleosome interactions

Liquid-liquid phase separation (LLPS) is an important mechanism that helps explain the membraneless compartmentalization of the nucleus. Because chromatin compaction and LLPS are collective phenomena, linking their modulation to the physicochemical features of nucleosomes is challenging. Here, we develop an advanced multiscale chromatin model-integrating atomistic representations, a chemically-specific coarse-grained model, and a minimal model-to resolve individual nucleosomes within sub-Mb chromatin domains and phase-separated systems. To overcome the difficulty of sampling chromatin at high resolution, we devise a transferable enhanced-sampling Debye-length replica-exchange molecular dynamics approach. We find that nucleosome thermal fluctuations become significant at physiological salt concentrations and destabilize the 30-nm fiber. Our simulations show that nucleosome breathing favors stochastic folding of chromatin and promotes LLPS by simultaneously boosting the transient nature and heterogeneity of nucleosome-nucleosome contacts, and the effective nucleosome valency. Our work puts forward the intrinsic plasticity of nucleosomes as a key element in the liquid-like behavior of nucleosomes within chromatin, and the regulation of chromatin LLPS.

Sequence-Dependent Three Interaction Site Model for Single- and Double-Stranded DNA

We develop a robust coarse-grained model for single- and double-stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inversion procedure using a large representative set of DNA structures extracted from the Protein Data Bank. Some of the parameters in the stacking interactions were calculated using a learning procedure, which ensured that the experimentally measured melting temperatures of dimers are faithfully reproduced. Without any further adjustments, the calculations based on the TIS model reproduce the experimentally measured salt and sequence-dependence of the size of single-stranded DNA (ssDNA), as well as the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon application of mechanical force, the extension of poly(dA) exhibits a plateau, which we trace to the formation of stacked helical domains. In contrast, the force-extension curve (FEC) of poly(dT) is entropic in origin and could be described by a standard polymer model. We also show that the persistence length of double-stranded DNA, formed from two complementary ssDNAs, is consistent with the prediction based on the worm-like chain. The persistence length, which decreases with increasing salt concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for stiff polyelectrolyte chains near the rod limit. Our model predicts the melting temperatures of DNA hairpins with excellent accuracy, and we are able to recover the experimentally known sequence-specific trends. The range of applications, which did not require adjusting any parameter after the initial construction based solely on PDB structures and melting profiles of dimers, attests to the transferability and robustness of the TIS model for ssDNA and dsDNA.

A new approach to precise thermodynamic characterization of hybridization properties of modified oligonucleotides: Comparative studies of deoxyribo- and glycine morpholine pentaadenines

The development of new derivatives and analogues of nucleic acids for the purposes of molecular biology, biotechnology, gene diagnostics, and medicine has been a hotspot for the last two decades. Methylenecarboxamide (glycine) morpholine oligomer analogues (gM) seem to be promising therapeutic candidates because of the ability to form sequence specific complexes with DNA and RNA. In this paper we describe new approaches to the determination of thermodynamic parameters for hybridization of tandem oligonucleotide complexes with the complementary template. It makes possible to determine changes in enthalpy and entropy corresponding to the binding of an individual oligomer with the template, and to the formation of cooperative contact at the helix-helix interface of two neighboring duplex fragments (in the nick). We have experimentally analyzed the series of model tandem complexes of different length at various oligomer concentrations, ionic strength, and pH. The analysis of thermodynamic parameters of complex formation for native and modified oligomers revealed higher Gibbs free energy values of hybridization and cooperative interaction of morpholine-containing complexes at the helix-helix interface under standard conditions (1M NaCl, pH7.2). Further comparative analysis of the hybridization properties of modified oligomers at ionic strength and pH allows us to determine the charge state of the morpholine backbone and the thermodynamic origin of the effects observed. It was found that the decrease in pH to 5.5 led to the protonation of internal morpholine nitrogens. The obtained results prove the veracity of the proposed model and the possibility to evaluate thermodynamic parameters of short native and modified oligomers with high accuracy.

Sensing conformational changes in DNA upon ligand binding using QCM-D. Polyamine condensation and Rad51 extension of DNA layers

Biosensors, in which binding of ligands is detected through changes in the optical or electrochemical properties of a DNA layer confined to the sensor surface, are important tools for investigating DNA interactions. Here, we investigate if conformational changes induced in surface-attached DNA molecules upon ligand binding can be monitored by the quartz crystal microbalance with dissipation (QCM-D) technique. DNA duplexes containing 59-184 base pairs were formed on QCM-D crystals by stepwise assembly of synthetic oligonucleotides of designed base sequences. The DNA films were exposed to the cationic polyamines spermidine and spermine, known to condense DNA molecules in bulk experiments, or to the recombination protein Rad51, known to extend the DNA helix. The binding and dissociation of the ligands to the DNA films were monitored in real time by measurements of the shifts in resonance frequency (Δf) and in dissipation (ΔD). The QCM-D data were analyzed using a Voigt-based model for the viscoelastic properties of polymer films in order to evaluate how the ligands affect thickness and shear viscosity of the DNA layer. Binding of spermine shrinks all DNA layers and increases their viscosity in a reversible fashion, and so does spermidine, but to a smaller extent, in agreement with its lower positive charge. SPR was used to measure the amount of bound polyamines, and when combined with QCM-D, the data indicate that the layer condensation leads to a small release of water from the highly hydrated DNA films. The binding of Rad51 increases the effective layer thickness of a 59 bp film, more than expected from the know 50% DNA helix extension. The combined results provide guidelines for a QCM-D biosensor based on ligand-induced structural changes in DNA films. The QCM-D approach provides high discrimination between ligands affecting the thickness and the structural properties of the DNA layer differently. The reversibility of the film deformation allows comparative studies of two or more analytes using the same DNA layer as demonstrated here by spermine and spermidine.

Ionic-content dependence of viscoelasticity of the lyotropic chromonic liquid crystal sunset yellow

A lyotropic chromonic liquid crystal (LCLC) is an orientationally ordered system made by self-assembled aggregates of charged organic molecules in water, bound by weak noncovalent attractive forces and stabilized by electrostatic repulsions. We determine how the ionic content of the LCLC, namely, the presence of mono- and divalent salts and pH enhancing agent, alter the viscoelastic properties of the LCLC. Aqueous solutions of the dye sunset yellow with a uniaxial nematic order are used as an example. By applying a magnetic field to impose orientational deformations, we measure the splay K1, twist K2, and bend K3 elastic constants and rotation viscosity γ1 as a function of concentration of additives. The data indicate that the viscoelastic parameters are influenced by ionic content in dramatic and versatile ways. For example, the monovalent salt NaCl decreases K3 and K2 and increases γ1, while an elevated pH decreases all the parameters. We attribute these features to the ion-induced changes in length and flexibility of building units of LCLC, the chromonic aggregates, a property not found in conventional thermotropic and lyotropic liquid crystals formed by covalently bound units of fixed length.

Construction and modeling of concatemeric DNA multilayers on a planar surface as monitored by QCM-D and SPR

The sequential hybridization of a 534 base pair DNA concatemer layer was monitored by QCM-D and SPR, and the QCM-D data were analyzed by Voigt viscoelastic models. The results show that Voigt-based modeling gives a good description of the experimental data but only if shear viscosity and elasticity are allowed to depend on the shear frequency. The derived layer thickness, shear viscosity and elasticity of the growing film give a representation of the DNA film in agreement with known bulk properties of DNA, and reveal a maximum in film viscosity when the molecules in the layer contain 75 base pairs. The experimental data during construction of a 3084 bp DNA concatemer layer were compared to predictions of the QCM-D response of a 1 μm thick film of rod-like polymers. A predicted nonmonotonous variation of dissipation with frequency (added mass) is in qualitative agreement with the experiments, but with a quantitative disagreement which likely reflects that the flexibility of such long DNA molecules is not included in the model.

Global analysis of the ground-state wrapping conformation of a charged polymer on an oppositely charged nano-sphere

We investigate the wrapping conformations of a single, strongly adsorbed polymer chain on an oppositely charged nano-sphere by employing a reduced (dimensionless) representation of a primitive chain-sphere model. This enables us to determine the global behavior of the chain conformation in a wide range of values for the system parameters including the chain contour length, its linear charge density and persistence length as well as the nano-sphere charge and radius, and also the salt concentration in the bathing solution. The structural behavior of a charged chain-sphere complex can be described in terms of a few distinct conformational symmetry classes separated by continuous or discontinuous transition lines which are determined by means of appropriately defined (order) parameters. Our results can be applied to a wide class of strongly coupled polymer-sphere complexes including, for instance, complexes that comprise a mechanically flexible or semiflexible polymer chain or an extremely short or long chain and, as a special case, include the biologically relevant example of DNA-histone complexes.

Do monovalent mobile ions affect DNA's flexibility at high salt content?

Numerous theoretical and experimental studies disagree on the impact of surrounding mobile ions on DNA conformational flexibility at high salt content. Specifically, it is not clear how the DNA persistence length varies when concentration of monovalent mobile ions is increased beyond the physiological value of ∼0.1 M. In the present Communication we address this biologically important issue computationally by means of molecular dynamics simulations. We utilize our recently developed chemically accurate coarse-grained model for the double-stranded DNA with explicit mobile ions. We find that in a range of moderate-to-high ionic concentrations, ∼0.1-1 M, DNA persistence length drops noticeably by ∼25%. Our results contradict some experimental works and the celebrated theory of Odijk, Skolnick and Fixman (Skolnick et al., Macromolecules, 1977, 10, 944), suggesting a negligible variation of DNA persistence length at these concentrations. On the other hand, our findings are in near quantitative agreement with a number of other theoretical and experimental studies. Combined with our recent work on elucidating the role of elastic and electrostatic effects in maintaining DNA shape, the results reported here may indicate that conceptually new understanding of DNA rigidity needs to be developed.

Effect of sodium and acetate ions on 8-hydroxyguanine formation in irradiated aqueous solutions of DNA and 2'-deoxyguanosine 5'-monophosphate

PURPOSE

The aim of this work was to study the combined effect of sodium and acetate ions on the radiation yield of 8-hydroxyguanine (8-OHG), one of the major DNA base lesions induced by free radicals.

MATERIALS AND METHODS

Aqueous solutions of DNA and 2'-deoxyguanosine 5'-monophosphate (dGMP) with various concentrations of sodium acetate and sodium perchlorate were γ-irradiated, enzymatically digested and analyzed by high-performance liquid chromatography (HPLC) methods.

RESULTS

It was found that both salts decrease the 8-OHG radiation yield in the concentration range studied for both DNA and dGMP, except in the case of dGMP wherein an increase in yield occurs in the concentration range from 0.1-1 mM. The dependence of the 8-hydroxy-2'-deoxyguanosine radiation yield on the concentration of both sodium acetate and sodium perchlorate have different shapes and have steeper slopes for the DNA compared with the dGMP solutions.

CONCLUSIONS

The observed decrease in the radiation yield of 8-OHG with increasing concentrations of sodium acetate is consistent with the hypothesis that sodium acetate produces two concentration-dependent effects in the DNA solutions: (1) A conformational change in the DNA caused by Na(+) counterions; and (2) free radical reactions related to the radiolysis of acetate ion.

Adsorption and desorption behaviors of DNA with magnetic mesoporous silica nanoparticles

The interaction between DNA and mesopores is one of the basic concerns when mesoporous silica nanoparticle (MSN) is used as a DNA carrier. In this work, we have synthesized a type of mesoporous silica nanoparticle that has a Fe(3)O(4) inner core and mesoporous silica shell. This magnetic mesoporous silica nanoparticle (denoted as M-MSN) offers us a convenient platform to manipulate the DNA adsorption and desorption processes as it can be easily separated from solution by applying a magnetic field. The DNA adsorption behavior is studied as a function of time in chaotropic salt solution. The maximum amount of adsorbed DNA is determined as high as 121.6 mg/g. We have also developed a method to separate the DNA adsorbed onto the external surface and into the mesopores by simply changing temperature windows. The desorption results suggest that, within the whole adsorbed DNA molecules, about 89.5% has been taken up by M-MSN mesopores. Through the dynamic light scattering experiment, we have found that the hydrodynamic size for M-MSN with DNA in its mesopores is higher than the naked M-MSN. Finally, the preliminary result of the adsorption mechanism study suggests that the DNA adsorption into mesopores may generate more intermolecular hydrogen bonds than those formed on the external surface.

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.

The DNA sequence-dependence of nucleosome positioning in vivo and in vitro

The contribution of histone-DNA interactions to nucleosome positioning in vivo is currently a matter of debate. We argue here that certain nucleosome positions, often in promoter regions, in yeast may be, at least in part, specified by the DNA sequence. In contrast other positions may be poorly specified. Positioning thus has both statistical and DNA-determined components. We further argue that the relative affinity of the octamer for different DNA sequences can vary and therefore the interaction of histones with the DNA is a 'tunable' property.

Swelling of biological and semiflexible polyelectrolytes

We have developed a theoretical model of swelling of semiflexible (biological) polyelectrolytes in salt solutions. Our approach is based on separation of length scales which allowed us to split a chain's electrostatic energy into two parts that describe local and remote electrostatic interactions along the polymer backbone. The local part takes into account interactions between charged monomers that are separated by distances along the polymer backbone shorter than the chain's persistence length. These electrostatic interactions renormalize chain persistence length. The second part includes electrostatic interactions between remote charged pairs along the polymer backbone located at distances larger than the chain persistence length. These interactions are responsible for chain swelling. In the framework of this approach we calculated effective chain persistence length and chain size as a function of the Debye screening length, chain degree of ionization, bare persistence length and chain degree of polymerization. Our crossover expression for the effective chain's persistence length is in good quantitative agreement with the experimental data on DNA. We have been able to fit experimental datasets by using two adjustable parameters: DNA ionization degree (α = 0.15-0.17) and a bare persistence length (l(p) = 40-44 nm).

Persistence lengths of DNA obtained from Brownian dynamics simulations

The persistence length of DNA has been studied for decades; however, experimentally obtained values of this quantity have not been entirely consistent. We report results from Brownian dynamics simulations that address this issue, validating and demonstrating the utility of an explicitly double-stranded model for mesoscale DNA dynamics. We find that persistence lengths calculated from rotational relaxation increase with decreasing ionic strength, corroborating experimental evidence, but contradicting results obtained from wormlike coil assumptions. Further, we find that natural curvature does not significantly affect the persistence length, corroborating cyclization efficiency measurements, but contradicting results from cryo-EM.

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.

The persistence length of DNA is reached from the persistence length of its null isomer through an internal electrostatic stretching force

To understand better the effect of electrostatics on the rigidity of the DNA double helix, we define DNA*, the null isomer of DNA, as the hypothetical structure that would result from DNA if its phosphate groups were not ionized. For the purposes of theoretical analysis, we model DNA* as identical to ordinary DNA but supplemented by a longitudinal compression force equal in magnitude but oppositely directed to the stretching (tension) force on DNA caused by phosphate-phosphate repulsions. The null isomer DNA* then becomes an elastically buckled form of fully ionized DNA. On this basis, we derive a nonadditive relationship between the persistence length P of DNA and the persistence length P* of its null isomer. From the formula obtained we can predict the value of P* if P is known, and we can predict the ionic strength dependence of P under the assumption that P* does not depend on ionic strength. We predict a value of P* for null DNA drastically lower than the value of P for DNA in its ordinary state of fully ionized phosphates. The predicted dependence of P on salt concentration is log-c over most of the concentration range, with no tendency toward a salt-independent value in the range of validity of the theory. The predictions are consistent with much of the persistence-length data available for DNA. Alternate theories of the Odijk-Skolnik-Fixman type, including one by the author, are considered skeptically on the grounds that the underlying model may not be realistic. Specifically, we doubt the accuracy for real polyelectrolytes of the Odijk-Skolnik-Fixman assumption that the polymer structure is invariant to changes in electrostatic forces.

Nanoscale mechanical and dynamical properties of DNA single molecules

Experimental evidence suggests DNA mechanical properties, in particular intrinsic curvature and flexibility, have a role in many relevant biological processes. Systematic investigations about the origin of DNA curvature and flexibility have been carried out; however, most of the applied experimental techniques need simplifying models to interpret the data, which can affect the results. Progress in the direct visualization of macromolecules allows the analysis of morphological properties and structural changes of DNAs directly from the digitised micrographs of single molecules. In addition, the statistical analysis of a large number of molecules gives information both on the local intrinsic curvature and the flexibility of DNA tracts at nanometric scale in relatively long sequences. However, it is necessary to extend the classical worm-like chain model (WLC) for describing conformations of intrinsically straight homogeneous polymers to DNA. This review describes the various methodologies proposed by different authors.

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.

The influence of DNA stiffness upon nucleosome formation

The rotational and translational positioning of nucleosomes on DNA is dependent to a significant extent on the physicochemical properties of the double helix. We have investigated the influence of the axial flexibility of the molecule on the affinity for the histone octamer by substituting selected DNA sequences with either inosine for guanosine or diaminopurine for adenine. These substitutions, respectively, remove or add a purine 2-amino group exposed in the minor groove and, respectively, decrease and increase the apparent persistence length. We observe that for all sequences tested inosine substitution, with one exception, increases the affinity for histone binding. Conversely diaminopurine substitution decreases the affinity. In the sole example where replacement of guanosine with inosine decreases the persistence length as well as the affinity for histones, the substitution concomitantly removes an intrinsic curvature of the DNA molecule. We show that, to a first approximation, the binding energy of DNA to histones at 1M NaCl is directly proportional to the persistence length. The data also indicate that a high local flexibility of DNA can favour strong rotational positioning.

DNA Codes for Nanoscience

The nanometer scale is a special place where all sciences meet and develop a particularly strong interdisciplinarity. While biology is a source of inspiration for nanoscientists, chemistry has a central role in turning inspirations and methods from biological systems to nanotechnological use. DNA is the biological molecule by which nanoscience and nanotechnology is mostly fascinated. Nature uses DNA not only as a repository of the genetic information, but also as a controller of the expression of the genes it contains. Thus, there are codes embedded in the DNA sequence that serve to control recognition processes on the atomic scale, such as the base pairing, and others that control processes taking place on the nanoscale. From the chemical point of view, DNA is the supramolecular building block with the highest informational content. Nanoscience has therefore the opportunity of using DNA molecules to increase the level of complexity and efficiency in self-assembling and self-directing processes.

The structural basis of DNA flexibility

Although the average physico-chemical properties of a long DNA molecule may approximate to those of a thin isotropic homogeneous rod, DNA behaves more locally as an anisotropic heterogeneous rod. This bending anisotropy is sequence dependent and to a first approximation reflects both the geometry and stability of individual base steps. The biological manipulation and packaging of the molecule often depend crucially on local variations in both bending and torsional flexibility. However, whereas the probability of DNA untwisting can be strongly correlated with a high bending flexibility, DNA bending, especially when the molecule is tightly wrapped on a protein surface, may be energetically favoured by a less flexible sequence whose preferred configuration conforms more closely to that of the complementary protein surface. In the latter situation the lower bending flexibility may be more than compensated for on binding by a reduced required deformation energy relative to a fully isotropic DNA molecule.

DNA persistence length revisited

DNA restriction fragments ranging from 79 to 789 base pairs in length have been characterized by transient electric birefringence (TEB) measurements at various temperatures between 4 and 43 degrees C. The DNA fragments do not contain runs of four or more adenine residues in a row and migrate with normal electrophoretic mobilities in polyacrylamide gels, indicating that they are not intrinsically curved or bent. The low ionic strength buffers used for the measurements contained 1 mM Tris Cl, pH 8.0, EDTA, and variable concentrations of Na(+) or Mg(2+) ions. The rotational relaxation times were obtained by fitting the TEB field-free decay signals with a nonlinear least-squared fitting program; the decay of the birefringence was monoexponential for fragments < or = 241 base pair (bp) in length and multiexponential for larger fragments. The terminal relaxation times, characteristic of the end-over-end rotation of the DNA molecules, were then used to determine the persistence length (p) and hydrodynamic radius (r) of DNA as a function of temperature and ionic strength, using several different hydrodynamic models. The specific values obtained for p and r are model dependent. The wormlike chain model of P. J. Hagerman and B. H. Zimm (Biopolymers 1981, Vol. 20, pp. 1481-1502) combined with the revised Broersma equation (J. Newman et al., Journal of Mol Biol 1997, Vol. 116, pp. 593-606) appears to be the most suitable for describing the flexibility of DNA in low ionic strength solutions. The values of p and r obtained from the global least squares fitting of this equation are independent of DNA length, and the deviations of the individual values from the average are reasonably small. The consensus r value calculated for DNA in various low ionic strength solutions containing 1 mM Tris buffer is 14.7 +/- 0.4 A at 20 degrees C. The consensus p values decrease from 814 approximately 564 A in solutions containing 1 mM Tris buffer plus 0.2-1 mM NaCl and decrease still further to 440 A in solutions containing 0.2 mM Mg(2+) ions. The persistence length exhibits a shallow maximum at 20 degrees C and decreases slowly upon either increasing or decreasing the temperature, regardless of the model used to fit the data. By contrast, the consensus values of the hydrodynamic radius are independent of temperature. The calculated persistence lengths and hydrodynamic radii are compared with other data in the literature.

Complexes of semiflexible polyelectrolytes and charged spheres as models for salt-modulated nucleosomal structures

We investigate the complexation behavior between a semiflexible charged polymer and an oppositely charged sphere with parameters appropriate for the DNA-histone system. We determine the ground state of a simple free energy expression (which includes electrostatic interactions on a linear level) numerically and use symmetry arguments to divide the obtained DNA configuration into broad classes, thereby obtaining global phase diagrams. We pay specific attention to the effects of salt concentration, DNA length variation, DNA charge renormalization, and externally applied force on the obtained complex structures.

Salt-induced DNA-histone complexation

We study numerically the binding of one semiflexible charged polymer onto an oppositely charged sphere. Using parameters appropriate for DNA-histone complexes, we find complete wrapping for intermediate salt concentrations only, in agreement with experiments. For high salt concentrations, a strongly discontinuous dewrapping occurs. For low salt concentrations, we find multiple conformational transitions, leading to an extended DNA configuration. The wrapped states are characterized by spontaneously broken rotational and mirror symmetries, giving rise to four distinct structures.

Shear-induced assembly of lambda-phage DNA

Recombinant DNA technology, which is based on the assembly of DNA fragments, forms the backbone of biological and biomedical research. Here we demonstrate that a uniform shear flow can induce and control the assembly of lambda-phage DNA molecules: increasing shear rates form integral DNA multimers of increasing molecular weight. Spontaneous assembly and grouping of end-blunted lambda-phage DNA molecules are negligible. It is suggested that shear-induced DNA assembly is caused by increasing the probability of contact between molecules and by stretching the molecules, which exposes the cohesive ends of the otherwise undeformed lambda-phage DNA molecules. We apply this principle to enhance the kinetics and extent of DNA concatenation in the presence of ligase. This novel approach to controlled DNA assembly could form the basis for improved approaches to gene-chip and recombinant DNA technologies and provide new insight into the rheology of associating polymers.

Advances in the separation of bacteriophages and related particles

Nondenaturing gel electrophoresis is used to both characterize multimolecular particles and determine the assembly pathways of these particles. Characterization of bacteriophage-related particles has yielded strategies for characterizing multimolecular particles in general. Previous studies have revealed means for using nondenaturing gel electrophoresis to determine both the effective radius and the average electrical surface charge density of any particle. The response of electrophoretic mobility to increasing the magnitude of the electrical field is used to detect rod-shaped particles. To increase the capacity of nondenaturing gel electrophoresis to characterize comparatively large particles, some current research is directed towards either determining the structure of gels used for electrophoresis or inducing steric trapping of particles in dead-end regions within the fibrous network that forms a gel. A trapping-dependent technique of pulsed-field gel electrophoresis is presented with which a DNA-protein complex can be made to electrophoretically migrate in a direction opposite to the direction of migration of protein-free DNA.

Static curvature and flexibility measurements of DNA with microscopy. A simple renormalization method, its assessment by experiment and simulation

We present the derivation of equations based on statistical polymer chain analysis and a method to quantify the average angle value of intrinsic bends and the local flexibility at a given locus on DNA fragments imaged by electron microscopy. DNA fragments of n base-pairs are considered as stiff chains of n jointed unit rigid rods. If the DNA fragments are composed of two branches A0Am and A0Bn, with, respectively, m and n base-pairs, where the standard deviations of the angle formed by two consecutive base-pairs are uniform over each branch, respectively, sigmathetaA and sigmathetaB, we show that the standard deviation of the angle AmA0Bn is: [formula: see text] where sigmatheta0 is the standard deviation of the angle at locus A0. This equation is established for small angular deviations by analysis of DNA at different scales and the validity of the methodology is controlled with the computation of the reduced chi2 statistical test. The length of the DNA fragments must be of the order of, or below, the persistence length, as determined by sets of statistics from computer simulations of DNA fragments. This is verified experimentally by a detailed analysis of the digitized contours of homogeneous linear 139 base-pair DNA fragments observed by electron microscopy. The images are compared to the reconstruction of DNA fragments from the measurements. The value found, sigma0=4.6 degrees/bp, is consistent with the well-accepted value for DNA in a plane. We discuss the relationship between the standard deviation of the measured angles and the flexibility at the base-pair level. This method is useful to quantify directly from microscopy techniques, such as electron or scanning force microscopy, the true bending angle, either intrinsic or induced by a ligand, and its associated flexibility at a given locus in any small DNA fragment.

Ionic effects on the elasticity of single DNA molecules

We used a force-measuring laser tweezers apparatus to determine the elastic properties of lambda-bacteriophage DNA as a function of ionic strength and in the presence of multivalent cations. The electrostatic contribution to the persistence length P varied as the inverse of the ionic strength in monovalent salt, as predicted by the standard worm-like polyelectrolyte model. However, ionic strength is not always the dominant variable in determining the elastic properties of DNA. Monovalent and multivalent ions have quite different effects even when present at the same ionic strength. Multivalent ions lead to P values as low as 250-300 A, well below the high-salt "fully neutralized" value of 450-500 A characteristic of DNA in monovalent salt. The ions Mg2+ and Co(NH3)63+, in which the charge is centrally concentrated, yield lower P values than the polyamines putrescine2+ and spermidine3+, in which the charge is linearly distributed. The elastic stretch modulus, S, and P display opposite trends with ionic strength, in contradiction to predictions of macroscopic elasticity theory. DNA is well described as a worm-like chain at concentrations of trivalent cations capable of inducing condensation, if condensation is prevented by keeping the molecule stretched. A retractile force appears in the presence of multivalent cations at molecular extensions that allow intramolecular contacts, suggesting condensation in stretched DNA occurs by a "thermal ratchet" mechanism.

Stretching DNA with optical tweezers

Force-extension (F-x) relationships were measured for single molecules of DNA under a variety of buffer conditions, using an optical trapping interferometer modified to incorporate feedback control. One end of a single DNA molecule was fixed to a coverglass surface by means of a stalled RNA polymerase complex. The other end was linked to a microscopic bead, which was captured and held in an optical trap. The DNA was subsequently stretched by moving the coverglass with respect to the trap using a piezo-driven stage, while the position of the bead was recorded at nanometer-scale resolution. An electronic feedback circuit was activated to prevent bead movement beyond a preset clamping point by modulating the light intensity, altering the trap stiffness dynamically. This arrangement permits rapid determination of the F-x relationship for individual DNA molecules as short as -1 micron with unprecedented accuracy, subjected to both low (approximately 0.1 pN) and high (approximately 50 pN) loads: complete data sets are acquired in under a minute. Experimental F-x relationships were fit over much of their range by entropic elasticity theories based on worm-like chain models. Fits yielded a persistence length, Lp, of approximately 47 nm in a buffer containing 10 mM Na1. Multivalent cations, such as Mg2+ or spermidine 3+, reduced Lp to approximately 40 nm. Although multivalent ions shield most of the negative charges on the DNA backbone, they did not further reduce Lp significantly, suggesting that the intrinsic persistence length remains close to 40 nm. An elasticity theory incorporating both enthalpic and entropic contributions to stiffness fit the experimental results extremely well throughout the full range of extensions and returned an elastic modulus of approximately 1100 pN.

Influence of optical probing with YOYO on the electrophoretic behavior of the DNA molecule

The influence of the fluorescent dye YOYO (1,1'-(4,4,8,8,-tetramethyl- 4,8-diazaundecamethylene)bis[4-[[3-methyl-benzo-1,3-oxazol-2 -yl] methylidene]-1,4-dihydroquinolinium] tetraiodide) on the electrophoretic behavior of the DNA molecule was investigated. This is important when using YOYO as a probe in capillary electrophoresis or in fluorescence microscopy studies of DNA with the purpose of studying the migration mechanism of DNA on the molecular level. We have measured the mobility and orientation dynamics (using the linear dichroism technique) for both pure DNA and the YOYO-DNA complex in agarose gel in order to compare their electrophoretic properties. Mobility decreases, the degree of orientation becomes lower, and the orientational dynamics slower, when YOYO binds to DNA. However, the dependence on field strength of the mobility, orientation and orientational dynamics, are similar for DNA and YOYO-DNA, indicating that the mode of migration does not change significantly upon binding YOYO to DNA. Furthermore, since our results show that the effect of YOYO on both the degree of orientation and orientational dynamics of the DNA can be measured and therefore be compensated for, it can be concluded that YOYO is an excellent optical probe for the study of the migrational behavior of DNA.

Overstretching B-DNA: the elastic response of individual double-stranded and single-stranded DNA molecules

Single molecules of double-stranded DNA (dsDNA) were stretched with force-measuring laser tweezers. Under a longitudinal stress of approximately 65 piconewtons (pN), dsDNA molecules in aqueous buffer undergo a highly cooperative transition into a stable form with 5.8 angstroms rise per base pair, that is, 70% longer than B form dsDNA. When the stress was relaxed below 65 pN, the molecules rapidly and reversibly contracted to their normal contour lengths. This transition was affected by changes in the ionic strength of the medium and the water activity or by cross-linking of the two strands of dsDNA. Individual molecules of single-stranded DNA were also stretched giving a persistence length of 7.5 angstroms and a stretch modulus of 800 pN. The overstretched form may play a significant role in the energetics of DNA recombination.