Anomalous Laterally Stressed Kinetically Trapped DNA Surface Conformations

DNA kinking is inevitable for the highly anisotropic 1D-1D electrostatic interaction with the one-dimensionally periodically charged surface. The double helical structure of the DNA kinetically trapped on positively charged monomolecular films comprising the lamellar templates is strongly laterally stressed and extremely perturbed at the nanometer scale. The DNA kinetic trapping is not a smooth 3D-> 2D conformational flattening but is a complex nonlinear in-plane mechanical response (bending, tensile and unzipping) driven by the physics beyond the scope of the applicability of the linear worm-like chain approximation. Up to now, the DNA molecule adsorbed on a surface was believed to always preserve its native structure. This belief implies a negligible contribution of lateral surface forces during and after DNA adsorption although their impact has never been elucidated. High-resolution atomic force microscopy was used to observe that stiff DNA molecules kinetically trapped on monomolecular films comprising one-dimensional periodically charged lamellar templates as a single layer or as a sublayer are oversaturated by sharp discontinuous kinks and can also be locally melted and supercoiled. We argue that kink/anti-kink pairs are induced by an overcritical lateral bending stress (> 30 pNnm) inevitable for the highly anisotropic 1D-1D electrostatic interaction of DNA and underlying rows of positive surface charges. In addition, the unexpected kink-inducing mechanical instability in the shape of the template-directed DNA confined between the positively charged lamellar sides is observed indicating the strong impact of helicity. The previously reported anomalously low values of the persistence length of the surface-adsorbed DNA are explained by the impact of the surface-induced low-scale bending. The sites of the local melting and supercoiling are convincingly introduced as other lateral stress-induced structural DNA anomalies by establishing a link with DNA high-force mechanics. The results open up the study in the completely unexplored area of the principally anomalous kinetically trapped DNA surface conformations in which the DNA local mechanical response to the surface-induced spatially modulated lateral electrostatic stress is essentially nonlinear. The underlying rich and complex in-plane nonlinear physics acts at the nanoscale beyond the scope of applicability of the worm-like chain approximation.

Multi-shell model of ion-induced nucleic acid condensation

We present a semi-quantitative model of condensation of short nucleic acid (NA) duplexes induced by trivalent cobalt(iii) hexammine (CoHex) ions. The model is based on partitioning of bound counterion distribution around single NA duplex into "external" and "internal" ion binding shells distinguished by the proximity to duplex helical axis. In the aggregated phase the shells overlap, which leads to significantly increased attraction of CoHex ions in these overlaps with the neighboring duplexes. The duplex aggregationfree energy is decomposed into attractive and repulsive components in such a way that they can be represented by simple analytical expressions with parameters derived from molecular dynamic simulations and numerical solutions of Poisson equation. The attractive term depends on the fractions of bound ions in the overlapping shells and affinity of CoHex to the "external" shell of nearly neutralized duplex. The repulsive components of the free energy are duplex configurational entropy loss upon the aggregation and the electrostatic repulsion of the duplexes that remains after neutralization by bound CoHex ions. The estimates of the aggregationfree energy are consistent with the experimental range of NA duplex condensation propensities, including the unusually poor condensation of RNA structures and subtle sequence effects upon DNAcondensation. The model predicts that, in contrast to DNA, RNA duplexes may condense into tighter packed aggregates with a higher degree of duplex neutralization. An appreciable CoHex mediated RNA-RNA attraction requires closer inter-duplex separation to engage CoHex ions (bound mostly in the "internal" shell of RNA) into short-range attractive interactions. The model also predicts that longer NA fragments will condense more readily than shorter ones. The ability of this model to explain experimentally observed trends in NAcondensation lends support to proposed NAcondensation picture based on the multivalent "ion binding shells."

Direct measurement of the dielectric polarization properties of DNA

The electric polarizability of DNA, represented by the dielectric constant, is a key intrinsic property that modulates DNA interaction with effector proteins. Surprisingly, it has so far remained unknown owing to the lack of experimental tools able to access it. Here, we experimentally resolved it by detecting the ultraweak polarization forces of DNA inside single T7 bacteriophages particles using electrostatic force microscopy. In contrast to the common assumption of low-polarizable behavior like proteins (εr ∼ 2-4), we found that the DNA dielectric constant is ∼ 8, considerably higher than the value of ∼ 3 found for capsid proteins. State-of-the-art molecular dynamic simulations confirm the experimental findings, which result in sensibly decreased DNA interaction free energy than normally predicted by Poisson-Boltzmann methods. Our findings reveal a property at the basis of DNA structure and functions that is needed for realistic theoretical descriptions, and illustrate the synergetic power of scanning probe microscopy and theoretical computation techniques.

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.

Atomic force microscopy study of the interaction of DNA and nanostructured beta-Gallia rutile

The ability to attach DNA molecules to solid planar substrates is desired for imaging the molecule and for building DNA-mediated nanostructures. The deposition of DNA on [001] rutile and beta-gallia rutile (BGR) substrates from buffer solutions containing various divalent cations was studied using tapping mode atomic force microscopy (AFM). beta-Gallia rutile intergrowths were prepared by spin-coating gallium isopropoxide onto [001]-oriented TiO2 single-crystal slabs and heating above 1350 degrees C for >24 h, resulting in the formation of intergrowth lines along the {210} planes in the parent rutile structure. Rutile and BGR intergrowth substrates were exposed to various buffered solutions containing DNA and the following divalent cations: Ca(II), Co(II), Cu(II), Fe(II), Mg(II), Mn(II), Ni(II), and Zn(II). Among all the cations examined, only Ni(II) resulted in the attachment of DNA on the rutile surfaces. DNA attachment to BGR surfaces was strong enough to allow AFM imaging when the deposition buffer contained one of the following cations: Co(II), Mg(II), Mn(II), Ni(II), and Zn(II). For all of these cations, DNA attachment occurred preferentially, but not exclusively, along BGR intergrowth lines. When buffers without cation additions and those containing Ca(II), Cu(II), and Fe(II) were used, DNA failed to bind the BGR surfaces strongly enough to allow AFM imaging. The mechanism(s) by which DNA attaches to the BGR surface is (are) not well understood but may involve the incorporation of divalent cations at the tunnel sites of the BGR intergrowths.

Electrostatic correlations and fluctuations for ion binding to a finite length polyelectrolyte

A statistical mechanical model is presented which explicitly accounts for the fluctuations, the electrostatic, and the excluded volume correlations for ions bound to a polyelectrolyte such as DNA. The method can be employed to treat a wide range of ionic conditions including multivalent ions. The microscopic framework of the theory permits the use of realistic finite length and grooved structural model for the polyelectrolyte and modeling of the finite size of the bound ions. Test against Monte Carlo simulations suggests that the theory can give accurate predictions for the ion distribution and the thermodynamic properties. For multivalent ions, the theory makes improved predictions as compared with the mean-field approach. Moreover, for long polyelectrolyte and dilute salt concentration, the theory predicts ion binding properties that agree with the counterion condensation theory.

Ion-mediated nucleic acid helix-helix interactions

Salt ions are essential for the folding of nucleic acids. We use the tightly bound ion (TBI) model, which can account for the correlations and fluctuations for the ions bound to the nucleic acids, to investigate the electrostatic free-energy landscape for two parallel nucleic acid helices in the solution of added salt. The theory is based on realistic atomic structures of the helices. In monovalent salt, the helices are predicted to repel each other. For divalent salt, while the mean-field Poisson-Boltzmann theory predicts only the repulsion, the TBI theory predicts an effective attraction between the helices. The helices are predicted to be stabilized at an interhelix distance approximately 26-36 A, and the strength of the attractive force can reach -0.37 k(B)T/bp for helix length in the range of 9-12 bp. Both the stable helix-helix distance and the strength of the attraction are strongly dependent on the salt concentration and ion size. With the increase of the salt concentration, the helix-helix attraction becomes stronger and the most stable helix-helix separation distance becomes smaller. For divalent ions, at very high ion concentration, further addition of ions leads to the weakening of the attraction. Smaller ion size causes stronger helix-helix attraction and stabilizes the helices at a shorter distance. In addition, the TBI model shows that a decrease in the solvent dielectric constant would enhance the ion-mediated attraction. The theoretical findings from the TBI theory agree with the experimental measurements on the osmotic pressure of DNA array as well as the results from the computer simulations.

Adsorption of monovalent and multivalent cations and anions on DNA molecules

Adsorption of monovalent and multivalent cations and anions on a deoxyribose nucleic acid (DNA) molecule from a salt solution is investigated by computer simulation. The ions are modeled as charged hard spheres, the DNA molecule as a point charge pattern following the double-helical phosphate strands. The geometrical shape of the DNA molecules is modeled on different levels ranging from a simple cylindrical shape to structured models which include the major and minor grooves between the phosphate strands. The densities of the ions adsorbed on the phosphate strands in the major and in the minor grooves are calculated. First, we find that the adsorption pattern on the DNA surface depends strongly on its geometrical shape: counterions adsorb preferentially along the phosphate strands for a cylindrical model shape, but in the minor groove for a geometrically structured model. Second, we find that an addition of monovalent salt ions results in an increase of the charge density in the minor groove while the total charge density of ions adsorbed in the major groove stays unchanged. The adsorbed ion densities are highly structured along the minor groove while they are almost smeared along the major groove. Furthermore, for a fixed amount of added salt, the major-groove cationic charge is independent of the counterion valency. For increasing salt concentration the major groove is neutralized while the total charge adsorbed in the minor groove is constant. DNA overcharging is detected for multivalent salts. Simulations for larger ion radii, which mimic the effect of ion hydration, indicate an increased adsorbtion of cations in the major groove.

Nonmonotonic variation with salt concentration of the second virial coefficient in protein solutions

The osmotic virial coefficient B2 of globular protein solutions is calculated as a function of added salt concentration at fixed pH by computer simulations of the "primitive model." The salt and counterions as well as a discrete charge pattern on the protein surface are explicitly incorporated. For parameters roughly corresponding to lysozyme, we find that B2 first decreases with added salt concentration up to a threshold concentration, then increases to a maximum, and then decreases again upon further raising the ionic strength. Our studies demonstrate that the existence of a discrete charge pattern on the protein surface profoundly influences the effective interactions and that linear and nonlinear Poisson Boltzmann theories fail for large ionic strength. The observed nonmonotonicity of B2 is compared with experiments. Implications for protein crystallization are discussed.

DNA properties investigated by dynamic force microscopy

In this work, we show that by varying the experimental conditions, the driving amplitude, a dynamic force microscope allows DNA properties to be selectively imaged. The substrate on which the DNA is fixed is a silica surface grafted with silane molecules terminated with amine groups. Use of small oscillation amplitudes favors the attractive interaction between the tip and the sample, while the use of large amplitudes renders the contribution of the attractive interaction negligible. In particular, at small amplitudes, the images show that the attractive interaction is strongly enhanced along the DNA. This enhancement is found to be amenable with a model considering a narrow strip of randomly oriented dipoles on each side of the molecule. This work should provide new insights on the DNA interaction and conformational changes with localized charges.

Electrostatic mechanisms of DNA deformation

The genomes of higher cells consist of double-helical DNA, a densely charged polyelectrolyte of immense length. The intrinsic physical properties of DNA, as well as the properties of its complexes with proteins and ions, are therefore of fundamental interest in understanding the functions of DNA as an informational macromolecule. Because individual DNA molecules often exceed 1 cm in length, it is clear that DNA bending, folding, and interaction with nuclear proteins are necessary for packaging genomes in small volumes and for integrating the nucleotide sequence information that guides genetic readout. This review first focuses on recent experiments exploring how the shape of the densely charged DNA polymer and asymmetries in its surrounding counterion distribution mutually influence one another. Attention is then turned to experiments seeking to discover the degree to which asymmetric phosphate neutralization can lead to DNA bending in protein-DNA complexes. It is argued that electrostatic effects play crucial roles in the intrinsic, sequence-dependent shape of DNA and in DNA shapes induced by protein binding.

Effective interaction between helical biomolecules

The effective interaction between two parallel strands of helical biomolecules, such as deoxyribose nucleic acids (DNA), is calculated using computer simulations of the "primitive" model of electrolytes. In particular we study a simple model for B-DNA incorporating explicitly its charge pattern as a double-helix structure. The effective force and the effective torque exerted onto the molecules depend on the central distance and on the relative orientation. The contributions of nonlinear screening by monovalent counterions to these forces and torques are analyzed and calculated for different salt concentrations. As a result, we find that the sign of the force depends sensitively on the relative orientation. For intermolecular distances smaller than 6 A it can be both attractive and repulsive. Furthermore, we report a nonmonotonic behavior of the effective force for increasing salt concentration. Both features cannot be described within linear screening theories. For large distances, on the other hand, the results agree with linear screening theories provided the charge of the biomolecules is suitably renormalized.

Stretching of single collapsed DNA molecules

The elastic response of single plasmid and lambda phage DNA molecules was probed using optical tweezers at concentrations of trivalent cations that provoked DNA condensation in bulk. For uncondensed plasmids, the persistence length, P, decreased with increasing spermidine concentration before reaching a limiting value 40 nm. When condensed plasmids were stretched, two types of behavior were observed: a stick-release pattern and a plateau at approximately 20 pN. These behaviors are attributed to unpacking from a condensed structure, such as coiled DNA. Similarly, condensing concentrations of hexaammine cobalt(III) (CoHex) and spermidine induced extensive changes in the low and high force elasticity of lambda DNA. The high force (5-15 pN) entropic elasticity showed worm-like chain (WLC) behavior, with P two- to fivefold lower than in low monovalent salt. At lower forces, a 14-pN plateau abruptly appeared. This corresponds to an intramolecular attraction of 0.083-0.33 kT/bp, consistent with osmotic stress measurements in bulk condensed DNA. The intramolecular attractive force with CoHex is larger than with spermidine, consistent with the greater efficiency with which CoHex condenses DNA in bulk. The transition from WLC behavior to condensation occurs at an extension about 85% of the contour length, permitting looping and nucleation of condensation. Approximately half as many base pairs are required to nucleate collapse in a stretched chain when CoHex is the condensing agent.

Thermodynamics of the conformational transition of biopolyelectrolytes: the case of specific affinity of counterions

A formal development of the Counterion Condensation theory (CC) of linear polyelectrolytes has been performed to include specific (chemical) affinity of condensed counterions, for polyelectrolyte charge density values larger than the critical value of condensation. It has been conventionally assumed that each condensed counterion exhibits an affinity free-energy difference for the polymer, (DeltaG(aff)). Moreover, the model assumes that the enthalpic and entropic contributions to DeltaG(aff), i.e., DeltaH(aff) and DeltaS(aff), are both independent of temperature, ionic strength and polymer concentration. Equations have been derived relative to the case of the thermally induced, ionic strength dependent, conformational transition of a biopolyelectrolyte between two conformations for which chemical affinity is supposed to take place. The experimental data of the intramolecular conformational transition of the ionic polysaccharide kappa-carrageenan in dimethylsulfoxide (DMSO) have been successfully compared with the theoretical predictions. This novel approach provides the enthalpic and entropic affinity values for both conformations, together with the corresponding thermodynamic functions of nonpolyelectrolytic origin pertaining to the biopolymer backbone change per se, i.e., DeltaH(n.pol) and DeltaS(n.pol), according to a treatment previously shown to be successful for lower values of the biopolyelectrolyte linear charge density. The ratio of DeltaH(n.pol) to DeltaS(n.pol) was found to be remarkably constant independent of the value of the dielectric constant of the solvent, from formamide to water to DMSO, pointing to the identity of the underlying conformational process.

DNA structure: cations in charge?

Recent X-ray diffraction, NMR spectroscopy and molecular mechanics results suggest that monovalent cations selectively partition into the minor groove of AT-tracts in DNA. These observations are consistent with DNA deformation by electrostatic collapse around areas of uneven cation density. This model predicts the occurrence of known DNA deformations, such as AT-tract bending and changes in the minor-groove width.