Effect of a low-dielectric interior on DNA electrostatic response to twisting and bending

We study the effects of a low-dielectric core of rod-like macromolecules on their electrostatic persistence lengths. We use the exact solution of the linear Poisson-Boltzmann equation for the potential of a charge on the surface of a low-dielectric cylinder. We apply the results to the B-DNA molecule, modeled as a double helical array of discrete charges wound on the surface of a low-dielectric rod. For this charge geometry, we calculate the change in the electrostatic twist persistence as compared to DNA with a water-permeable core. We also discuss possible effects of the low-dielectric molecular core on DNA bending persistence.     

Density functional theory for planar electric double layers: closing the gap between simple and polyelectrolytes

We report a nonlocal density functional theory (NLDFT) for polyelectrolyte solutions within the primitive model; i.e., the solvent is represented by a continuous dielectric medium, and the small ions and polyions by single and tangentially connected charged hard spheres, respectively. The excess Helmholtz energy functional is derived from a modified fundamental measure theory for hard-sphere repulsion, an extended first-order thermodynamic perturbation theory for chain connectivity, and a quadratic functional Taylor expansion for electrostatic correlations. With the direct and cavity correlation functions of the corresponding monomeric systems as inputs, the NLDFT predicts the segment-level microscopic structures and adsorption isotherms of polyelectrolytes at oppositely charged surfaces in good agreement with molecular simulations. In particular, it faithfully reproduces the layering structures of polyions, charge inversion, and overcharging that cannot be captured by alternative methods including the polyelectrolyte Poisson-Boltzmann equation and an earlier version of DFT. The NLDFT has also been used to investigate the influences of the small ion valence, polyion chain length, and size disparity between polyion segments and counterions on the microscopic structure, mean electrostatic potential, and overcharging in planar electric double layers containing polyelectrolytes.     

Complexation of semiflexible chains with oppositely charged cylinder

We study the complexation of long thin semiflexible polymer chains with an oppositely charged cylinder. Starting from the linear Poisson-Boltzmann equation, we calculate the electrostatic potential and the energy of such a charge distribution. We find that sufficiently flexible chains prefer to wrap around the cylinder in a helical manner, when their charge density is smaller than that of the cylinder. The optimal value of the helical pitch is found by minimization of the sum of electrostatic and bending energies. The dependence of the pitch on the number of chains, their rigidity, and salt concentration in solution is analyzed. We discuss our results in the light of recent experiments on DNA complexation with cylindrical dendronized polymers.     

Encapsulation of a polyelectrolyte chain by an oppositely charged spherical surface

Using the ground state dominance approximation and a variational theory, we study the encapsulation of a polyelectrolyte chain by an oppositely charged spherical surface. The electrostatic attraction between the polyelectrolyte and the surface and the entropy loss of the encapsulated polyelectrolyte chain dictate the optimum conditions for encapsulation. Two scenarios of encapsulation are identified: entropy-dominated and adsorption-dominated encapsulation. In the entropy-dominated encapsulation regime, the polyelectrolyte chain is delocalized, and the optimum radius of the encapsulating sphere decreases with increasing the attraction. In the adsorption-dominated encapsulation regime, the polyelectrolyte chain is strongly localized near the surface, and the optimum radius increases with increasing the attraction. After identifying a universal encapsulation parameter, the dependencies of the optimum radius on the salt concentration, surface charge density, polymer charge density, and polymer length are explored.     

Nanoparticle adsorption on a weak polyelectrolyte. Stiffness, pH, charge mobility, and ionic concentration effects investigated by Monte Carlo simulations

Monte Carlo simulations have been used to study two different models for a weak linear polyelectrolyte in the presence of nanoparticles: (i) a rodlike and (ii) a flexible polyelectrolytes. The use of simulated annealing has made it possible to simulate a polyelectrolyte chain in the presence of several nanoparticles by improving conformation sampling and avoiding multiple minima problems when dense conformations are produced. Nanoparticle distributions along the polymer backbone were analyzed versus the ionic concentration, polyelectrolyte stiffness, and nanoparticle surface charge. Titration curves were calculated and the influences of the ionic concentration, solution pH, and number of adsorbed nanoparticles on the acid/base polyelectrolyte properties have been systematically investigated. The subtle balance of attractive and repulsive interactions has been discussed, and some characteristic conformations are presented. The comparison of the two limit models provides a good representation of the stiffness influence on the complex formation. In some conditions, overcharging was obtained and presented with respect to both the polyelectrolyte and nanoparticle as the central element. Finally, the charge mobility influence along the polyelectrolyte backbone was investigated by considering annealed and quenched polyelectrolyte chains.     

Overcharging, charge inversion, and reentrant condensation: using highly charged polyelectrolytes in tetravalent salt solutions as an example of study

We study salt-induced charge overcompensation and charge inversion of flexible polyelectrolytes via computer simulations and demonstrate the importance of ion excluded volume. Reentrant condensation takes place when the ion size is comparable to monomer size, and happens in a middle region of salt concentration. In a high-salt region, ions can overcharge a chain near its surface and charge distribution around a chain displays an oscillatory behavior. Unambiguous evidence obtained by electrophoresis shows that charge inversion does not necessarily appear with overcharging and occurs when the ion size is not big. These findings suggest a disconnection of resolubilization of polyelectrolyte condensates at high salt concentration with charge inversion.     

Overcharging of nanoparticles in electrolyte solutions

Monte Carlo simulations are performed to investigate the effects of salt concentration, valence and size of small ions, surface charge density, and Bjerrum length on the overcharging of isolated spherical nanoparticles within the framework of a primitive model. It is found that charge inversion is most probable in solutions containing multivalent counterions at high salt concentrations. The maximum strength of overcharging occurs near the nanoparticle surface where counterions and coions have identical local concentrations. The simulation results also suggest that both counterion size and electrostatic correlations play major roles for the occurrence of overcharging.     

Strong and weak adsorptions of polyelectrolyte chains onto oppositely charged spheres

We investigate the complexation of long thin polyelectrolyte (PE) chains with oppositely charged spheres. In the limit of strong adsorption, when strongly charged PE chains adapt a definite wrapped conformation on the sphere surface, we analytically solve the linear Poisson-Boltzmann equation and calculate the electrostatic potential and the energy of the complex. We discuss some biological applications of the obtained results. For weak adsorption, when a flexible weakly charged PE chain is localized next to the sphere in solution, we solve the Edwards equation for PE conformations in the Hulthen potential, which is used as an approximation for the screened Debye-Huckel potential of the sphere. We predict the critical conditions for PE adsorption. We find that the critical sphere charge density exhibits a distinctively different dependence on the Debye screening length than for PE adsorption onto a flat surface. We compare our findings with experimental measurements on complexation of various PEs with oppositely charged colloidal particles. We also present some numerical results of the coupled Poisson-Boltzmann and self-consistent field equation for PE adsorption in an assembly of oppositely charged spheres.     

Theory of competitive counterion adsorption on flexible polyelectrolytes: divalent salts

The counterion distribution around an isolated flexible polyelectrolyte in the presence of a divalent salt is evaluated using the adsorption model [M. Muthukumar, J. Chem. Phys. 120, 9343 (2004)] that considers the Bjerrum length, salt concentration, and local dielectric heterogeneity as physical variables in the system. Self-consistent calculations of effective charge and size of the polymer show that divalent counterions replace condensed monovalent counterions in competitive adsorption. The theory further predicts that at modest physical conditions for a flexible polyelectrolytes such as sodium polystyrene sulfonate in aqueous solutions polymer charge is compensated and reversed with increasing divalent salt. Consequently, the polyelectrolyte shrinks and reswells. Lower temperatures and higher degrees of dielectric heterogeneity between chain backbone and solvent enhance condensation of all species of ions. Complete diagrams of states for the effective charge calculated as functions of the Coulomb strength and salt concentration suggest that (a) overcharging requires a minimum Coulomb strength and (b) progressively higher presence of salt recharges the polymer due to either electrostatic screening (for low Coulomb strengths) or coion condensation (for high Coulomb strengths). Consideration of ion-bridging by divalent counterions leads to a first-order collapse of polyelectrolytes in modest presence of divalent salts and at higher Coulomb strengths. The authors' theoretical predictions are in agreement with the generic results from experiments and simulations.     

Effective charges of polyelectrolytes in a salt-free solution based on counterion chemical potential

The phenomenon of counterion condensation around a flexible polyelectrolyte chain with N monomers is investigated by Monte Carlo simulations in terms of the degree of ionization alpha, which is proportional to the effective charge. It is operationally defined as the ratio of observed to intrinsic counterion concentration, alpha = co/ci. The observed counterion concentration in the dilute polyelectrolyte solution is equivalent to an electrolyte solution of concentration co with the same counterion chemical potential. It can be determined directly by thermodynamic experiments such as ion-selective electrode. With the polyelectrolyte fixed at the center of the spherical Wigner-Seitz cell, the polymer conformation, counterion distribution, and chemical potential can be obtained. Our simulation shows that the degree of ionization rises as the polymer concentration decreases. This behavior is opposite to that calculated from the infinitely long charged rod model, which is often used to study counterion condensation. Moreover, we find that, for a specified line charge density, alpha decreases with an increment in chain length and chain flexibility. In fact, the degree of ionization is found to decline with increasing polymer fractal dimension, which can be tuned by varying bending modulus and solvent quality. Those results can be qualitatively explained by a simple model of two-phase approximation.     

Single polyelectrolyte macromolecule in the salt solution: Effect of escaped counter ions

We developed simplest theory of swelling of polyelectrolyte macromolecules in the solution of low‐molecular salts. The novel feature taken into account is the electrostatic interaction between the macro‐ion and counter ions which escaped from the interior of the macromolecular coil and distributed in the whole volume of solution. The phase diagram in the variables solvent quality vs. salt concentrationis derived. One of the unexpected findings is that the total charge within the sphere surrounding the macro‐ion changes in a non‐monotonical manner with the increase of the salt concentration. In particular, our calculations predict that the electrophoretic mobility of a polyelectrolyte macromolecule can demonstrate non‐monotonic behavior with an increase of salt concentration. The results of the calculations are compared with those obtained under the assumption of electroneutrality of the macromolecule when all of the counter ions are kept within the effective volume of macromolecule.     

The Effect of Surface Charge on Adsorption of a Cationic Polyelectrolyte

The regularities of adsorption of a cationic polyelectrolyte, poly(diallyldimethylammonium chloride), on the surface of fused quartz are studied at different values of solution pH by capillary electrokinetics. It is shown that the polyelectrolyte adsorption on a negatively charged surface depends on the value of the surface charge and increases with its growth. At a low charge value (pH 3.8), the polyelectrolyte adsorption increases the quartz surface charge. The driving forces of the adsorption are both electrostatic interaction and forces of nonelectrostatic nature, probably hydrophobic interactions and a change in entropy due to the displacement of counterions from a double layer. The adsorption of poly(diallyldimethylammonium chloride) on quartz from alkaline and neutral solutions is irreversible, which indicates the key role of the electrostatic interaction. At low values of the surface charge, the nonelectrostatic interactions play the main role, thereby resulting in polyelectrolyte desorption.     

Chain stiffness, salt valency, and concentration influences on titration curves of polyelectrolytes: Monte Carlo simulations

Monte Carlo simulations have been used to study two different models of a weak linear polyelectrolyte surrounded by explicit counterions and salt particles: (i) a rigid rod and (ii) a flexible chain. We focused on the influence of the pH, chain stiffness, salt concentration, and valency on the polyelectrolyte titration process and conformational properties. It is shown that chain acid-base properties and conformational properties are strongly modified when multivalent salt concentration variation ranges below the charge equivalence. Increasing chain stiffness allows to minimize intramolecular electrostatic monomer interactions hence improving the deprotonation process. The presence of di and trivalent salt cations clearly promotes the chain degree of ionization but has only a limited effect at very low salt concentration ranges. Moreover, folded structures of fully charged chains are only observed when multivalent salt at a concentration equal or above charge equivalence is considered. Long-range electrostatic potential is found to influence the distribution of charges along and around the polyelectrolyte backbones hence resulting in a higher degree of ionization and a lower attraction of counterions and salt particles at the chain extremities.     

Brownian dynamics simulations of polyelectrolyte adsorption in shear flow: effects of solvent quality and charge patterning

We probe the effects of solvent quality and charge patterning on polyelectrolyte adsorption in shear flow using Brownian dynamics simulations with hydrodynamic interaction (HI). The polyelectrolyte is modeled as a freely jointed bead-rod chain, and electrostatic and non-electrostatic interactions are accounted for by using screened Coulombic and Lennard-Jones potentials, respectively. In the absence of flow, the conformation of a polyelectrolyte molecule adsorbed onto a uniformly charged surface changes from flat to globular with an increase in bead-bead attraction (hydrophobicity), consistent with prior experimental observations. In the presence of flow, migration due to bead-wall HI and, as a consequence, desorption decrease with an increase in bead-bead attraction, implying that flow-induced desorption is more difficult under poor-solvent conditions. When bead-bead non-electrostatic attraction is strong, desorption can be enhanced by increasing bead-bead electrostatic repulsion. Analogous to the effect of bead-surface electrostatic attraction, an increase in the strength of bead-surface non-electrostatic attraction reduces desorption. We also study the effect of shear flow on the adsorption of a polyelectrolyte molecule onto surfaces decorated with periodic arrays of charged patches. An increase in patch periodicity increases desorption even when the effective surface charge density is kept the same. The results of this work suggest mechanisms for controlling the desorption of polyelectrolyte molecules in shear flows.     

Release of lysozyme from the branched polyelectrolyte-lysozyme complexation

On the basis of the discretely charged sphere model of lysozyme, the release behavior of lysozyme from the branched polyelectrolyte-lysozyme complexation is investigated by adding salt and changing the pH values of the solution. It is found that, with the increase of the salt ionic strength of the solution, the lysozymes are gradually released from the oppositely charged polyelectrolyte as a result of the screening of electrostatic attraction between the two ionic species by adding the salt. Interestingly, there exists a critical salt ionic strength at which all proteins are released from the branched polyelectrolyte, and the polyelectrolyte-protein complexation is broken completely. Beyond the critical value, the increase of the salt ionic strength causes self-association of the proteins released from the branched polyelectrolyte-protein complexation. The self-association of the protein is detrimental in biological systems. By calculating the second virial coefficient, we found that the optimal salt content for the dispersion of proteins coincides with the critical ionic strength, because the second virial coefficient reaches its maximum at the critical ionic strength. Similarly, increasing the pH value of the solution can also release the lysozymes from the polyelectrolyte, because the increase of pH value of the solution changes the charge distribution and net charge of the lysozyme, weakens the attraction between lysozymes mediated by polyelectrolyte, and finally leads to the dissolution of the complexation of branched polyelectrolyte with lysozymes in strong alkaline solution. In addition, by exploring the effect of architecture of the polyelectrolyte on the release behavior of proteins, we found that it is more difficult to release proteins from the branched polyelectrolyte than from the linear polyelectrolyte.     

Is a small number of charge neutralizations sufficient to bend nucleosome core DNA onto its superhelical ramp?

X-ray diffraction structures of the nucleosome core particle along with a variety of experiments are consistent with the idea that an important source of the free energy holding DNA to the superhelical ramp on the histone octamer surface is obtained from a relatively small amount of electrostatic neutralization of the DNA phosphate charge by positively charged histone groups, especially arginine residues. Here we present a theoretical analysis of a simple model that emphasizes the competition between the high degree of bending of the stiff DNA molecule required for its tight curvature on the histone octamer and the neutralization of the DNA phosphate charge by basic histone residues. Our calculation accounts for the strong influence of condensed counterions on the electrostatic interactions. We find that the minimum amount of free energy required to bend DNA into axial conformity with the superhelical ramp at physiological salt concentration can be provided by a scant 6% neutralization of the phosphate charge, in close correspondence to the stoichiometric neutralization of phosphate charge by the arginine side chain that intrudes into the inward-facing minor groove of each DNA double helical turn.     

An Exact Solution to the Electrostatic Interaction between an Ion-Penetrable Sphere and an Ion-Penetrable Rod

No exact solution for the free energy of electrostatic interaction for a charged sphere and rod geometry in an electrolyte solution has yet been proposed. This geometry is interesting because it can be applied to describe macromolecules interacting with a random fiber-matrix for modeling of hindered transport in diffusional systems. Here we present an analytical approach that yields an exact solution to the problem for ion-penetrable-also called "soft"-sphere and infinitely long rod. This solution is compared to a published finite-element analysis of the same system with nonpenetrable-also called "hard"-sphere and infinitely long rod maintaining a constant surface charge density restriction. For any ionic strength or ratio of rod radius to sphere radius the ion-penetrable method yields an electrostatic free energy of interaction which is lower than that given by the analysis for hard bodies. This free energy is significantly lower for most parameter value combinations and therefore suggests that one should carefully examine the system being modeled to determine if it is approximated better by a hard body or ion-penetrable body approach. Copyright 2000 Academic Press.     

Molecular dynamics simulations of polyelectrolyte multilayering on a charged particle

Molecular dynamics simulations of polyelectrolyte multilayering on a charged spherical particle revealed that the sequential adsorption of oppositely charged flexible polyelectrolytes proceeds with surface charge reversal and highlighted electrostatic interactions as the major driving force of layer deposition. Far from being completely immobilized, multilayers feature a constant surge of chain intermixing during the deposition process, consistent with experimental observations of extensive interlayer mixing in these films. The formation of multilayers as well as the extent of layer intermixing depends on the degree of polymerization of the polyelectrolyte chains and the fraction of charge on its backbone. The presence of ionic pairs between oppositely charged macromolecules forming layers seems to play an important role in stabilizing the multilayer film.     

Experimental and theoretical evidence of overcharging of calcium silicate hydrate

Electrokinetic measurements such as electrophoresis may show an inversion of the effective surface charge of colloidal particle called overcharging. This phenomenon has been studied by various theoretical approaches but up to now very few attempts of confrontation between theory and experiment have been conducted. In this work we report electrophoretic measurements as well as Monte Carlo simulations of the electrokinetic potential for the surface of calcium silicate hydrate (CSH), which is the major constituent of hydrated cement. In the simulations, the surface charge of CSH nanoparticles in equilibrium with the ionic solution is determined by a single site characteristic and electrostatic interactions between all explicit charges at the surface and in the electric double layer. We will show that ordinary electrostatic interactions are enough to describe all experimental observations. Actually, an excellent agreement is found between experimental and simulated results without any fitting parameter, both with respect to surface titration and electrokinetic behaviour. The agreement extends over a wide range of electrostatic coupling, from a weakly charged surface with mainly monovalent counter-ions to a highly charged one with divalent counter-ions.