Phase behavior of columnar DNA assemblies

Experimental detection of conformational transitions between forms of DNA: problems and prospects

Under different conditions, the DNA double helix can take different geometric forms. Of the large number of its conformations, in addition to the "canonical" B form, the A, C, and Z forms are widely known, and the D, Hoogsteen, and X forms are less known. DNA locally takes the A, C, and Z forms in the cell, in complexes with proteins. We compare different methods for detecting non-canonical DNA conformations: X-ray, IR, and Raman spectroscopy, linear and circular dichroism in both the infrared and ultraviolet regions, as well as NMR (measurement of chemical shifts and their anisotropy, scalar and residual dipolar couplings and inter-proton distances from NOESY (nuclear Overhauser effect spectroscopy) data). We discuss the difficulties in applying these methods, the problems of theoretical interpretation of the experimental results, and the prospects for reliable identification of non-canonical DNA conformations.

Order and interactions in DNA arrays: Multiscale molecular dynamics simulation

While densely packed DNA arrays are known to exhibit hexagonal and orthorhombic local packings, the detailed mechanism governing the associated phase transition remains rather elusive. Furthermore, at high densities the atomistic resolution is paramount to properly account for fine details, encompassing the DNA molecular order, the contingent ordering of counterions and the induced molecular ordering of the bathing solvent, bringing together electrostatic, steric, thermal and direct hydrogen-bonding interactions, resulting in the observed osmotic equation of state. We perform a multiscale simulation of dense DNA arrays by enclosing a set of 16 atomistically resolved DNA molecules within a semi-permeable membrane, allowing the passage of water and salt ions, and thus mimicking the behavior of DNA arrays subjected to external osmotic stress in a bathing solution of monovalent salt and multivalent counterions. By varying the DNA density, local packing symmetry, and counterion type, we obtain osmotic equation of state together with the hexagonal-orthorhombic phase transition, and full structural characterization of the DNA subphase in terms of its positional and angular orientational fluctuations, counterion distributions, and the solvent local dielectric response profile with its order parameters that allow us to identify the hydration force as the primary interaction mechanism at high DNA densities.

Modeling DNA condensation on freestanding cationic lipid membranes

Motivated by recent experimental observations of a rapid spontaneous DNA coil-globule transition on freestanding cationic lipid bilayers, we propose simple theoretical models for DNA condensation on cationic lipid membranes. First, for a single DNA rod, we examine the conditions of full wrapping of a cylindrical DNA-like semi-flexible polyelectrolyte by an oppositely charged membrane. Then, for two parallel DNA rods, we self-consistently analyze the shape and the extent of the membrane enveloping them, focusing on membrane elastic deformations and the membrane-DNA embracing angle, which enables us to compute the membrane-mediated DNA-DNA interactions. We examine the effects of the membrane composition and its charge density, which are the experimentally tunable parameters. We show that membrane-driven rod-rod attraction is more pronounced for higher charge densities and for smaller surface tensions of the membrane. Thus, we demonstrate that for a long DNA chain adhered to a cationic lipid membrane, such membrane-induced DNA-DNA attraction can trigger compaction of DNA.

Modulating the Arrangement of Charged Nanotubes by Ionic Strength in Salty Water

Despite the important role and potential application of charged cylindrical polyelectrolytes, biomacromolecules, and self-assembles, salt-modulated organization of those 1D charged nanostructures remains a topic relatively unexplored with an obscure underlying mechanism. In this Letter, the aggregation of oriented nanotubes self-assembled by ionic aromatic oligoamide in aqueous solution of NaCl over a wide concentration range is probed via small-angle X-ray scattering and a transmission electron microscope. The arrangement of nanotubes undergoes order-disorder transition sequences from an ordered rectangular phase to hexagonal packing and then to a lamellar gel. The observed transitions are understood by ionic effects on the electrostatic interaction between charged nanotubes and osmotic pressure due to ion partitioning. Above the physiological condition, electrostatic interactions are largely screened by the salts, while osmotic effects start to regulate the aggregation behavior and concomitantly deform the nanotubes. The study demonstrates rich phase behaviors of ordered, charged 1D nanostructures by tuning the ionic strength and underlying key physical principles.

Torsional deformation of double helix in interaction and aggregation of DNA

We incorporate sequence-dependent twisting between adjacent base pairs and torsional elasticity of double helix into the theory of DNA-DNA interaction. The results show that pairing and counterion-induced-aggregation of nonhomologous DNA are accompanied by considerable torsional deformation. The deformation tunes negatively charged phosphate strands and positively charged grooves on opposing molecules to stay "in register", substantially reducing nonideality of the helical structure of DNA. Its cost, however, makes interaction between nonhomologous DNA less energetically favorable. In particular, interaction between double helical DNA may result in sequence homology recognition and selective pairing of homologous fragments containing more than 100-200 base pairs. We also find a weak, but potentially measurable, increase in the expected counterion concentration required for aggregation of nonhomologous DNA and slightly higher solubility of such DNA above the critical concentration.

Understanding thermodynamic competitivity between biopolymer folding and misfolding under large-scale intermolecular interactions

Cooperativity is a hallmark of spontaneous biopolymer folding. The presence of intermolecular interactions could create off-pathway misfolding structures and suppress folding cooperativity. This raises the hypothesis that thermodynamic competitivity between off-pathway misfolding and on-pathway folding may intervene with cooperativity and govern biopolymer folding dynamics under conditions permitting large-scale intermolecular interactions. Here we report direct imaging and theoretical modeling of thermodynamic competitivity between biopolymer folding and misfolding under such conditions, using a two-dimensional array of proton-fueled DNA molecular motors packed at the maximal density as a model system. Time-resolved liquid-phase atomic force microscopy with enhanced phase contrast revealed that the misfolding and folding intermediates transiently self-organize into spatiotemporal patterns on the nanoscale in thermodynamic states far away from equilibrium as a result of thermodynamic competitivity. Computer simulations using a novel cellular-automaton network model provide quantitative insights into how large-scale intermolecular interactions correlate the structural dynamics of individual biomolecules together at the systems level.

Micro-rheology on (polymer-grafted) colloids using optical tweezers

Optical tweezers are experimental tools with extraordinary resolution in positioning (± 1 nm) a micron-sized colloid and in the measurement of forces (± 50 fN) acting on it-without any mechanical contact. This enables one to carry out a multitude of novel experiments in nano- and microfluidics, of which the following will be presented in this review: (i) forces within single pairs of colloids in media of varying concentration and valency of the surrounding ionic solution, (ii) measurements of the electrophoretic mobility of single colloids in different solvents (concentration, valency of the ionic solution and pH), (iii) similar experiments as in (i) with DNA-grafted colloids, (iv) the nonlinear response of single DNA-grafted colloids in shear flow and (v) the drag force on single colloids pulled through a polymer solution. The experiments will be described in detail and their analysis discussed.

Physics of DNA: unravelling hidden abilities encoded in the structure of ‘the most important molecule’

A comprehensive article “Structure and Interactions of Biological Helices”, published in 2007 in Reviews of Modern Physics, overviewed various aspects of the effect of DNA structure on DNA–DNA interactions in solution and related phenomena, with a thorough analysis of the theory of these effects. Here, an updated qualitative account of this area is presented without any sophisticated ‘algebra’. It overviews the basic principles of the structure-specific interactions between double-stranded DNA and focuses on the physics behind several related properties encoded in the structure of DNA. Among them are (i) DNA condensation and aptitude to pack into small compartments of cells or viral capcids, (ii) the structure of DNA mesophases, and (iii) the ability of homologous genes to recognize each other prior to recombination from a distance. Highlighted are some of latest developments of the theory, including the shape of the ‘recognition well’. The article ends with a brief discussion of the first experimental evidence of the protein-free homology recognition in a ‘test tube’.

Structural reorganization of parallel actin bundles by crosslinking proteins: incommensurate states of twist

We construct a coarse-grained model of parallel actin bundles crosslinked by compact globular bundling proteins, such as fascin and espin, necessary components of filopodial and mechanosensory bundles. Consistent with structural observations of bundles, we find that the optimal geometry for crosslinking is overtwisted, requiring a coherent structural change of the helical geometry of the filaments. We study the linker-dependent thermodynamic transition of bundled actin filaments from their native state to the overtwisted state and map out the "twist-state" phase diagram in terms of the availability as well as the flexibility of crosslinker proteins. We predict that the transition from the uncrosslinked to fully crosslinked state is highly sensitive to linker flexibility: flexible crosslinking smoothly distorts the twist state of bundled filaments, while rigidly crosslinked bundles undergo a phase transition, rapidly overtwisting filaments over a narrow range of free crosslinker concentrations. Additionally, we predict a rich spectrum of intermediate structures, composed of alternating domains of sparsely bound (untwisted) and strongly bound (overtwisted) filaments. This model reveals that subtle differences in crosslinking agents themselves modify not only the detailed structure of parallel actin bundles, but also the thermodynamic pathway by which they form.

Undulations enhance the effect of helical structure on DNA interactions

During the past decade, theory and experiments have provided clear evidence that specific helical patterns of charged groups and adsorbed (condensed) counterions on the DNA surface are responsible for many important features of DNA-DNA interactions in hydrated aggregates. The effects of helical structure on DNA-DNA interactions result from a preferential juxtaposition of the negatively charged sugar phosphate backbone with counterions bound within the grooves of the opposing molecule. Analysis of X-ray diffraction experiments confirmed the mutual alignment of parallel molecules in hydrated aggregates required for such juxtaposition. However, it remained unclear how this alignment and molecular interactions might be affected by intrinsic and thermal fluctuations, which cause structural deviations away from an ideal double helical conformation. We previously argued that the torsional flexibility of DNA allows the molecules to adapt their structure to accommodate a more electrostatically favorable alignment between molecules, partially compensating disruptive fluctuation effects. In the present work, we develop a more comprehensive theory, incorporating also stretching and bending fluctuations of DNA. We found the effects of stretching to be qualitatively and quantitatively similar to those of twisting fluctuations. However, this theory predicts more dramatic and surprising effects of bending. Undulations of DNA in hydrated aggregates strongly amplify rather than weaken the helical structure effects. They enhance the structural adaptation, leading to better alignment of neighboring molecules and pushing the geometry of the DNA backbone closer to that of an ideal helix. These predictions are supported by a quantitative comparison of the calculated and measured osmotic pressures in DNA.

Symmetries of interacting helices of charge

We analytically examine the pair interaction for parallel discrete helices of charge. Symmetry arguments allow for the free energy to be decomposed into a sum of terms, each of which has an intuitive geometric interpretation. Truncated Fourier expansions of these terms are shown to provide effective free energy expressions that are valid under very general circumstances. These expressions are used to briefly examine and characterize the azimuthal interactions within F-actin and A-DNA aggregates.

DNA interaction studies of new nano metal based anticancer agent: validation by spectroscopic methods

A new nano dimensional heterobimetallic Cu-Sn containing complex as a potential drug candidate was designed, synthesized and characterized by analytical and spectral methods. The electronic absorption and electron paramagnetic resonance parameters of the complex revealed that the Cu(II) ion exhibits a square pyramidal geometry with the two pyrazole nitrogen atoms, the amine nitrogen atom and the carboxylate oxygen of the phenyl glycine chloride ligand located at the equatorial sites and the coordinated chloride ion occupying an apical position. (119)Sn NMR spectral data showed a hexa-coordinated environment around the Sn(IV) metal ion. TEM, AFM and XRD measurements illustrate that the complex could induce the condensation of CT-DNA to a particulate nanostructure. The interaction of the Cu-Sn complex with CT-DNA was investigated by UV-vis absorption and emission spectroscopy, as well as cyclic voltammetric measurements. The results indicated that the complex interacts with DNA through an electrostatic mode of binding with an intrinsic binding constant K(b) = 8.42 x 10(4) M( - 1). The Cu-Sn complex exhibits effective cleavage of pBR322 plasmid DNA by an oxidative cleavage mechanism, monitored at different concentrations both in the absence and in the presence of reducing agents.

Self-assembled structures of Gaussian nematic particles

We investigate the stable crystalline configurations of a nematic liquid crystal made of soft parallel ellipsoidal particles interacting via a repulsive, anisotropic Gaussian potential. For this purpose, we use genetic algorithms (GA) in order to predict all relevant and possible solid phase candidates into which this fluid can freeze. Subsequently we present and discuss the emerging novel structures and the resulting zero-temperature phase diagram of this system. The latter features a variety of crystalline arrangements, in which the elongated Gaussian particles in general do not align with any one of the high-symmetry crystallographic directions, a compromise arising from the interplay and competition between anisotropic repulsions and crystal ordering. Only at very strong degrees of elongation does a tendency of the Gaussian nematics to align with the longest axis of the elementary unit cell emerge.

Salt-induced aggregation of stiff polyelectrolytes

Molecular dynamics simulation techniques are used to study the process of aggregation of highly charged stiff polyelectrolytes due to the presence of multivalent salt. The dominant kinetic mode of aggregation is found to be the case of one end of one polyelectrolyte meeting others at right angles, and the kinetic pathway to bundle formation is found to be similar to that of flocculation dynamics of colloids as described by Smoluchowski. The aggregation process is found to favor the formation of finite bundles of 10-11 filaments at long times. Comparing the distribution of the cluster sizes with the Smoluchowski formula suggests that the energy barrier for the aggregation process is negligible. Also, the formation of long-lived metastable structures with similarities to the raft-like structures of actin filaments is observed within a range of salt concentration.

Structure of toroidal DNA collapsed inside the phage capsid

The structure of DNA toroids made of individual DNA molecules of various lengths (3,000 to 55,000 bp) was studied, by using partially filled bacteriophage capsids in conjunction with cryoelectron microscopy. The tetravalent cation spermine was diffused through the capsid to condense the DNA under conditions that were chosen to produce a hexagonal packing. Our results demonstrate that the frustration arising between chirality and hexagonal packing leads to the formation of twist walls; the correlation between helices combined with their strong curvature impose variations of the DNA helical pitch.

DNA cholesteric phases: the role of DNA molecular chirality and DNA-DNA electrostatic interactions

DNA molecules form dense liquid-crystalline twisted phases both in vivo and in vitro. How the microscopic DNA chirality is transferred into intermolecular twist in these mesophases and what is the role of chiral DNA-DNA electrostatic interactions is still not completely clear. In this paper, we first give an extended overview of experimental observations on DNA cholesteric phases and discuss the factors affecting their stability. Then, we consider the effects of steric and electrostatic interactions of grooved helical molecules on the sign of cholesteric twist. We present some theoretical results on the strength of DNA-DNA chiral electrostatic interactions, on DNA-DNA azimuthal correlations in cholesteric phases, on the value of DNA cholesteric pitch, and on the regions of existence of DNA chiral phases stabilized by electrostatic interactions. We suggest for instance that 146 bp long DNA fragments with stronger affinities for the nucleosome formation can form less chiral cholesteric phases, with a larger left-handed cholesteric pitch. Also, the value of left-handed pitch formed in assemblies of homologous DNA fragments is predicted to be smaller than that of randomly sequenced DNAs. We expect also the cholesteric assemblies of several-kbp-long DNAs to require higher external osmotic pressures for their stability than twisted phases of short nucleosomal DNA fragments at the same DNA lattice density.

Ground states of colloidal molecular crystals on periodic substrates

Two-dimensional suspensions of spherical colloids subject to periodic external fields exhibit a rich variety of molecular crystalline phases. We study in simulations the ground state configurations of dimeric and trimeric systems, that are realized on square and triangular lattices, when either two or three macroions are trapped in each external potential minimum. Bipartite orders of the checkerboard or stripe types are reported together with more complex quadripartite orderings, and the shortcomings of envisioning the colloids gathered in a single potential minimum as a composite rigid object are discussed. This work also sheds light on simplifying assumptions underlying previous theoretical treatments and that made possible the mapping onto spin models.

Hexagonal order in crystalline and columnar phases of hard rods

We report a study of colloidal suspensions of highly monodisperse semiflexible chiral rodlike viruses, denoted fd, in the range of high concentrations. Small angle x-ray scattering experiments reveal the existence of two hexagonal phases: the first one is crystalline and the second one is hexatic columnar, as shown by its short-range positional order. The suspension of rodlike viruses is the first experimental system showing the whole phase sequence with increasing particle concentration theoretically predicted for systems of hard rods, ranging from the chiral nematic via the smectic to columnar and crystalline phases.

Polyelectrolyte-compression forces between spherical DNA brushes

Optical tweezers are employed to measure the forces of interaction within a single pair of DNA-grafted colloids, dependent on the molecular weight of the DNA chains, and the concentration and valence of the surrounding ionic medium. The resulting forces are short range and set in as the surface-to-surface distance between the colloidal cores reaches the value of the brush height. The measured force-distance relation is analyzed by means of a theoretical treatment that quantitatively describes the effects of compression of the chains on the surface of the opposite-lying colloid. Quantitative agreement with the experiment is obtained for all parameter combinations.

The role of electrostatics in the B to A transition of DNA: from solution to assembly

On the basis of a wealth of published experimental data and computer simulations, we build a simple physical model that allows us to rationalize the A to B transition of DNA in solution and in aggregates. In both cases we find that the electrostatic interactions are strong enough, alone, to induce the transition independently of other energetic contributions, e.g. those related to hydration. On the basis of this analysis we conclude that in ethanol/water mixtures, the effect responsible for the transition is the reduction of dielectric constant in the mixture. This is manifested in electrostatic self-energy terms that include the interaction of phosphate charges with condensed counterions. But in dense aggregates, electrostatics plays a dual role, giving rise to two competing effects. In the absence of groove localized counterions the electrostatic self-energy favours the B form, and the electrostatic interaction energy between neighbouring DNA favours the A form. However, the addition of enough counterions localized in the narrow groove reverses this. In dry aggregates of DNA both terms, in most cases, conspire to keep DNA in the A form. The analysis gives a broad picture of the B to A transition and sets a number of new research goals, particularly concerning simulations that may test our simple model for aggregates.

DNA melting in aggregates: impeded or facilitated?

How does DNA melt in columnar aggregate relative to its melting in diluted solution? Is the melting temperature increased or decreased with the aggregate density? Have DNA-DNA interactions, predominantly of electrostatic nature, an effect on the character of the melting transition? In attempt to answer these questions, we have incorporated the theory of electrostatic interactions between DNA duplexes into the simplest model of DNA melting. The analysis shows that the effect of aggregate density is very different for aggregates built of homologous (or identical) DNA fragments relative to the case of DNA with random base pair sequences. The putative attraction between homologous DNA helices hampers their melting and increases the melting temperature and can even dramatically change the character of the transition. In the aggregate of nonhomologous DNAs, the pattern of electrostatic interactions is more complicated, and their effect could be opposite; in some cases we may even expect electrostatically induced melting. These findings define new directions for melting experiments in dense DNA assemblies.

Biological Soft Materials

With one or two exceptions, biological materials are "soft", meaning that they combine viscous and elastic elements. This mechanical behavior results from self-assembled supramolecular structures that are stabilized by noncovalent interactions. It is an ongoing and profound challenge to understand the self-organization of biological materials. In many cases, concepts can be imported from soft-matter physics and chemistry, which have traditionally focused on materials such as colloids, polymers, surfactants, and liquid crystals. Using these ideas, it is possible to gain a new perspective on phenomena as diverse as DNA condensation, protein and peptide fibrillization, lipid partitioning in rafts, vesicle fusion and budding, and others, as discussed in this selective review of recent highlights from the literature.

Screwlike order, macroscopic chirality, and elastic distortions in high-density DNA mesophases

We investigate a new screwlike liquid-crystalline ordering in solutions of helical biopolymers and its influence on the state of individual molecules. In the resulting mesophase translational and rotational motions of molecules are coupled in screw fluctuations. We show that in contrast to the case of conventional chiral liquid crystals the elastic distortion does not twist the screw order but leads to overwinding of individual helical molecules. This explains the peculiarities of high-density DNA mesophases.

Lattice-gas model of DNA charge inversion by a positively charged polyelectrolyte

The model of DNA charge inversion by Nguyen and Shklovskii [T. T. Nguyen and B. I. Shklovskii, Phys. Rev. Lett. 89, 018101 (2002)] is extended. A single double-helix strand of DNA is represented by a lattice of negative charges at the positions of the protruding oxygens of the phosphates along the DNA backbone, and the adsorbed polyelectrolyte molecules are represented by charged dimers. A lattice-gas model is used in which dimers adsorbing either parallel or perpendicular to the lattice are treated as separate species, and the model allows for vacancies between adsorbed species. The mean field theory used is formulated as a saddle-point approximation of the exact functional integral representation of the grand canonical partition function, opening the way for the inclusion of the effects of charge fluctuation corrections.

Elongation/compaction of giant DNA caused by depletion interaction with a flexible polymer

Structural changes in giant DNA induced by the addition of the flexible polymer Polyethylene Glycol (PEG) were examined by the method of single-DNA observation. In dilute DNA conditions, individual DNA assumes a compact state via a discrete coil-globule transition, whereas in concentrated solution, DNA molecules exhibit an extended conformation via macroscopic phase segregation. The long-axis length of the stretched state in DNA is about 10(3) times larger than that of the compact state. Phase segregation at high DNA concentrations occurs at lower PEG concentrations than the compaction at low DNA concentrations. These opposite changes in the conformation of DNA molecule are interpreted in terms of the free energy, including depletion interaction.

Direct observation of azimuthal correlations between DNA in hydrated aggregates

This study revisits the classical x-ray diffraction patterns from hydrated, noncrystalline fibers originally used to establish the helical structure of DNA. We argue that changes in these diffraction patterns with DNA packing density reveal strong azimuthally dependent interactions between adjacent molecules up to approximately 40 A interaxial or approximately 20 A surface-to-surface separations. These interactions appear to force significant torsional "straightening" of DNA and strong azimuthal alignment of nearest neighbor molecules. The results are in good agreement with the predictions of recent theoretical models relating DNA-DNA interactions to the helical symmetry of their surface charge patterns.

Statistical mechanics of columnar DNA assemblies

Many physical systems can be mapped onto solved or "solvable" models of magnetism. In this work, we have mapped the statistical mechanics of columnar phases of ideally helical rigid DNA--subject to the earlier found unusual, frustrated pair potential (A.A. Kornyshev, S. Leikin, J. Chem. Phys. 107, 3656 (1997))--onto an exotic, unknown variant of the XY model on a fixed or restructurable lattice. Here, the role of the "spin" is played by the azimuthal orientation of the molecules. We have solved this model using a Hartree-Fock approximation, ground-state calculations, and finite-temperature Monte Carlo simulations. We have found peculiar spin order transitions, which may also be accompanied by positional restructuring, from hexagonal to rhombohedric lattices. Some of these have been experimentally observed in dense columnar aggregates. Note that DNA columnar phases are of great interest in biophysical research, not only because they are a useful in vitro tool for the study of DNA condensation, but also since these structures have been detected in living matter. Within the approximations made, our study provides insight into the statistical mechanics of these systems.

DNA-DNA interaction beyond the ground state

The electrostatic interaction potential between DNA duplexes in solution is a basis for the statistical mechanics of columnar DNA assemblies. It may also play an important role in recombination of homologous genes. We develop a theory of this interaction that includes thermal torsional fluctuations of DNA using field-theoretical methods and Monte Carlo simulations. The theory extends and rationalizes the earlier suggested variational approach which was developed in the context of a ground state theory of interaction of nonhomologous duplexes. It shows that the heuristic variational theory is equivalent to the Hartree self-consistent field approximation. By comparison of the Hartree approximation with an exact solution based on the QM analogy of path integrals, as well as Monte Carlo simulations, we show that this easily analytically-tractable approximation works very well in most cases. Thermal fluctuations do not remove the ability of DNA molecules to attract each other at favorable azimuthal conformations, neither do they wash out the possibility of electrostatic "snap-shot" recognition of homologous sequences, considered earlier on the basis of ground state calculations. At short distances DNA molecules undergo a "torsional alignment transition," which is first order for nonhomologous DNA and weaker order for homologous sequences.

Nonlinear effects in the torsional adjustment of interacting DNA

DNA molecules in solution, having negatively charged phosphates and countercations readsorbed on its surface, possess a distinct charge separation motif to interact electrostatically. If their double-helical structure were ideal, duplexes in parallel juxtaposition could choose azimuthal alignment providing attraction, or at least a reduction of repulsion, between them. But duplexes are not perfect staircases and the distortions of their helical structure correlate with their base pair texts. If the patterns of distortions on the opposing molecules are uncorrelated, the mismatch will accumulate as a random walk and attraction vanishes. Based on this idea, a model of recognition of homologous sequences has been proposed [A. A. Kornyshev and S. Leikin, Phys. Rev. Lett. 86, 3666 (2001)]. But DNA has torsional elasticity. How will this help to relax a mismatch between the charge distributions on two nonhomologous DNA's? In the same work, the solution of this problem has been mapped onto a frustrated sine Gordon equation in a nonlocal random field (where the latter represents a pattern of twist angle distortions on the opposing molecules), but the results had been obtained in the limit of torsionally rigid molecules. In the present paper, by solving this equation numerically, we find a strongly nonlinear relaxation mechanism which utilizes static kink-soliton modes triggered by the "random field." In the range of parameters where the solitons do not emerge, we find good agreement with the results of a variational study [A. G. Cherstvy, A. A. Kornyshev, and S. Leikin, J. Phys. Chem. B (to be published)]. We reproduce the first-order transitions in the interaxial separation dependence, but detect also second-order or weak first-order transitions for shorter duplexes. The recognition energy between two nonhomologous DNA sequences is calculated as a function of interaxial separation and the length of juxtaposition. The soliton-caused kinky length dependence is discussed in connection with plots of recombination frequency as a function of the length of homology.

Competition between compaction of single chains and bundling of multiple chains in giant DNA molecules

It has been established that in a dilute solution individual giant DNA molecules undergo a large discrete transition between an elongated coil state and a folded compact state. On the other hand, in concentrated solutions, DNA molecules assemble into various characteristic states, including multichain aggregate, liquid crystalline, ionic crystal, etc. In this study, we compared single-chain and multiple-chain events by observing individual chains using fluorescence microscopy. We used spermidine, SPD(3+), as a condensing agent for giant DNA. When the concentration of DNA is below 1 microM in base-pair units, individual DNA molecules exhibit a transition from an elongated state to a compact state. When the concentration of DNA is increased to 10 microM, a thick fiberlike assembly of multiple chains appears. AFM measurements of this thick fiber revealed that more than tens of DNA molecules form a bundle structure with parallel ordering of the chains. The transition between single-chain compaction and bundle formation with multiple-chain assemblies was reproduced by a theoretical calculation.

Water-dispersed lamellar phases of symmetric poly(styrene)-block-poly(acrylic acid) diblock copolymers: model systems for flat dense polyelectrolyte brushes

We investigate the static properties of a water-dispersed lamellar ( L) phase formed in the melt state with a nearly symmetric poly(styrene)-block-poly(acrylic acid) (PS- b-PAA) diblock copolymer. The PAA brush is considered as a model flat polyelectrolyte ( PE) brush of controlled surface density. Thanks to small-angle X-ray scattering, its behavior in water is studied as a function of (i) its ionization, through the pH of the dispersions which is increased by an addition of a known amount of a base, i.e. sodium hydroxyde NaOH, and (ii) in the presence of a monovalent salt, i.e. sodium chloride NaCl, of concentration C(S). At low pH, we find that the brush effectively behaves as a neutral brush. At high pH, the brush is in the so-called "osmotic regime", in which all sodium counterions are trapped within the brush volume and stretch the chains via an osmotic effect. The properties of such a brush in the presence of a monovalent salt, confirm this result, showing a C(S)(-1/3) dependence in the brush height L(O), in agreement with mean-field predictions. The L(O)- C(S) profiles at different ionizations give access to the actual brush internal charge fraction f. The results are found to be in very good quantitative agreement with experimental measures found in the literature, and can be completely and quantitatively described by Oosawa's approach to counterion condensation in a semi-dilute to concentrated solution of charged, rod-like chains.

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.

Azimuthal frustration and bundling in columnar DNA aggregates

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

X-ray diffraction characterization of the dense phases formed by nucleosome core particles

Multiple dense phases of nucleosome core particles (NCPs) were formed in controlled ionic conditions (15-160 mM monovalent salt, no divalent ions), under osmotic pressures ranging from 4.7 x 10(5) to 2.35 x 10(6) Pa. We present here the x-ray diffraction analysis of these phases. In the lamello-columnar phase obtained at low salt concentration (<25 mM), NCPs stack into columns that align to form bilayers, kept separated from one another by a layer of solvent. NCPs form a monoclinic lattice in the plane of the bilayer. For high salt concentration (>50 mM), NCPs order into either a two-dimensional columnar hexagonal phase or into three-dimensional orthorhombic (quasi-hexagonal) crystals. The lamellar and hexagonal (or quasi-hexagonal) organizations coexist in the intermediate salt range; their demixing requires a long time. For an applied pressure P = 4.7 10(5) Pa, the calculated NCPs concentration ranges from approximately 280 to 320 mg/ml in the lamello-columnar phase to 495 to 585 mg/ml in the three-dimensional orthorhombic phase. These concentrations cover the concentration of the living cell.

Counterion penetration and effective electrostatic interactions in solutions of polyelectrolyte stars and microgels

Counterion distributions and effective electrostatic interactions between spherical macroions in polyelectrolyte solutions are calculated via second-order perturbation (linear response) theory. By modeling the macroions as continuous charge distributions that are permeable to counterions, analytical expressions are obtained for counterion profiles and effective pair interactions in solutions of star-branched and microgel macroions. The counterions are found to penetrate stars more easily than microgels, with important implications for screening of bare macroion interactions. The effective pair interactions are Yukawa in form for separated macroions, but are softly repulsive and bounded for overlapping macroions. A one-body volume energy, which depends on the average macroion concentration, emerges naturally in the theory and contributes to the total free energy.