Specific interactions versus counterion condensation. 2. Theoretical treatment within the counterion condensation theory

Polyuronates such as pectate and alginate are very well-known examples of biological polyelectrolytes undergoing, upon addition of divalent cations, an interchain association that acts as the junction of an eventually formed stable hydrogel. In the present paper, a thermodynamic model based on the counterion condensation theory has been developed to account for this cation-induced chain pairing of negatively charged polyelectrolytes. The strong interactions between cross-linking ions and uronate moieties in the specific binding site have been described in terms of chemical bonding, with complete charge annihilation between the two species. The chain-pairing process is depicted as progressively increasing with the concentration of cross-linking counterions and is thermodynamically defined by the fraction of each species. On these bases, the total Gibbs energy of the system has been expressed as the sum of the contributions of the Gibbs energy of the (single) chain stretches and of the (associated) dimers, weighted by their respective fractions 1 - theta and theta. In addition, the model assumes that the condensed divalent counterions exhibit an affinity free-energy for the chain, G(C)(aff,0), and the junction, G(D)(aff,0), respectively. Moreover, a specific Gibbs energy of chemical bonding, G(bond,0), has been introduced as the driving force for the formation of dimers. The model provides the mathematical formalism for calculating the fraction, theta, of chain dimers formed and the amount of ions condensed and bound onto the polyelectrolyte when two different types of counterions (of equal or different valence) are present. The effect of the parameter G(bond,0) has been investigated and, in particular, its difference from G(C,D)(aff,0) was found to be crucial in determining the distribution of the ions into territorial condensation and chemical bonding, respectively. Finally, the effect of the variation of the molar ratio between cross-linking ions and uronic groups in the specific binding sites, sigma0, was evaluated. In particular, a remarkable decrease in the amount of condensed counterions has been pointed out in the case of sigma0 = 1/3, with respect to the value of sigma0 = 1/4, characterizing the traditional "egg-box" structure, as a result of the drop of the charge density of the polyelectrolyte induced by complete charge annihilation.     

Negative surface charge near sodium channels of nerve: divalent ions, monovalent ions, and pH.

Evidence is given for a high density of negative surface charge near the sodium channel of myelinated nerve fibres. The voltage dependence of peak sodium permeability is measured in a voltage clamp. The object is to measure voltage shifts in sodium activation as the following external variables are varied: divalent cation concentration and type, monovalent concentration, and pH. With equimolar substitution of divalent ions the order of effectiveness for giving a positive shift is: Ba equals Sr less than Mg less than Ca less than Co approximately equal to Mn less than Ni less than Zn. A tenfold increase of concentration of any of these ions gives a shift of +20 to +25 mV. At low pH, the shift with a tenfold increase in Ca-2+ is much less than at normal pH, and conversely for high pH. Soulutions with no added divalent ions give a shift of minus 18 mV relative to 2 mM Ca-2+. Removal of 7/8 of the cations from the calcium-free solution gives a further shift of minue 35 mV. All shifts are explained quantitatively by assuming that changes in an external surface potential set up by fixed charges near the sodium channel produce the shifts. The model involves a diffuse double layer of counterions at the nerve surface and some binding of H+ions and divalent ions to the fixed charges. Three types of surface groups are postulated: (1) an acid pKa equals 2.88 charge density minus 0.9 nm- minus 2; (i) an acid pKa equals 4.58, charge density minus 0.58 nm- minus 2; (3) a base pKa equals 6.28, charge density +0.33 nm- minus 2. The two acid groups also bind Ca-2+ ions with a dissociation constant K equals 28 M. Reasonable agreement can also be obtained with a lower net surface charge density and stronger binding of divalent ions and H+ ions.     

Effect of finite ionic size on the solution of the Poisson-Boltzmann equation: Application to the binding of divalent metal ions to DNA

The nonlinear Poisson-Boltzmann equation is solved for a cylindrical polyelectrolyte solution containing mono- and divalent counterions and monovalent coions. The finite size of the ions is taken into account by the introduction of the distances of closest approach between the ionic charges and the surface of the polyelectrolyte. The choice of these distances is based on the physicochemical properties of the polyelectrolyte and ions in solution. The effects of the finite ionic size on the distribution of the counterions around the polyelectrolyte and on the local ion concentration and the integrated charge fraction of the divalent cations in the vicinity of the polyelectrolyte are discussed. Theoretical predictions regarding the overall extent of binding and the extent of inner-sphere binding of divalent counterions to rodlike polyions are compared with the results of nmr studies of the binding of divalent metal ions to DNA.     

Interactions of mono- and divalent counterions with alkali- and enzyme-deesterified pectins in salt-free solutions

The free fractions of monovalent and divalent counterions were determined on salt-free solutions of pectins. The effects of charge density, distribution of the carboxyl groups, polymer concentration, and the nature of the counterion were investigated by determinating the calcium and sodium activity coefficients (with specific electrodes) and by measuring the transport parameters (by conductimetry). Poor agreement for calcium ions was found with the Manning theory. The strong binding of these ions to highly charged polymers, which is ascribed to a dimerization process was demonstrated in very dilute solutions.     

Adsorption of DNA to mica mediated by divalent counterions: a theoretical and experimental study

The adsorption of DNA molecules onto a flat mica surface is a necessary step to perform atomic force microscopy studies of DNA conformation and observe DNA-protein interactions in physiological environment. However, the phenomenon that pulls DNA molecules onto the surface is still not understood. This is a crucial issue because the DNA/surface interactions could affect the DNA biological functions. In this paper we develop a model that can explain the mechanism of the DNA adsorption onto mica. This model suggests that DNA attraction is due to the sharing of the DNA and mica counterions. The correlations between divalent counterions on both the negatively charged DNA and the mica surface can generate a net attraction force whereas the correlations between monovalent counterions are ineffective in the DNA attraction. DNA binding is then dependent on the fractional surface densities of the divalent and monovalent cations, which can compete for the mica surface and DNA neutralizations. In addition, the attraction can be enhanced when the mica has been pretreated by transition metal cations (Ni(2+), Zn(2+)). Mica pretreatment simultaneously enhances the DNA attraction and reduces the repulsive contribution due to the electrical double-layer force. We also perform end-to-end distance measurement of DNA chains to study the binding strength. The DNA binding strength appears to be constant for a fixed fractional surface density of the divalent cations at low ionic strength (I < 0.1 M) as predicted by the model. However, at higher ionic strength, the binding is weakened by the screening effect of the ions. Then, some equations were derived to describe the binding of a polyelectrolyte onto a charged surface. The electrostatic attraction due to the sharing of counterions is particularly effective if the polyelectrolyte and the surface have nearly the same surface charge density. This characteristic of the attraction force can explain the success of mica for performing single DNA molecule observation by AFM. In addition, we explain how a reversible binding of the DNA molecules can be obtained with a pretreated mica surface.     

Evaluation of the counterion condensation theory of polyelectrolytes.

We compare free energies of counterion distributions in polyelectrolyte solutions predicted from the cylindrical Poisson-Boltzmann (PB) model and from the counterion condensation theories of Manning: CC1 (Manning, 1969a, b), which assumes an infinitely thin region of condensed counterions, and CC2 (Manning, 1977), which assumes a region of finite thickness. We consider rods of finite radius with the linear charge density of B-DNA in 1-1 valent and 2-2 valent salt solutions. We find that under all conditions considered here the free energy of the CC1 and the CC2 models is higher than that of the PB model. We argue that counterion condensation theory imposes nonphysical constraints and is, therefore, a poorer approximation to the underlying physics based on continuum dielectrics, point-charge small ions, Poisson electrostatics, and Boltzmann distributions. The errors in counterion condensation theory diminish with increasing distance from, or radius of, the polyion.     

Density functional theory for the nonspecific binding of salt to polyelectrolytes: thermodynamic properties

The thermodynamics of the nonspecific binding of salt to a polyelectrolyte molecule is studied using a density functional approach. The polyelectrolyte molecule is modeled as an infinite, inflexible, and impenetrable charged cylinder and the counterions and co-ions are modeled as charged hard spheres of equal diameter. The density functional theory is based on a hybrid approach where the hard-sphere contribution to the one-particle correlation function is evaluated nonperturbatively and the ionic contribution to the one-particle correlation function is evaluated perturbatively. The advantage of the approach is that analytical expressions are available for all the correlation functions. The calculated single ion preferential interaction coefficients, excess free energy, and activity coefficients show a nonmonotonic variation as a function of polyion charge in the presence of divalent ions. These properties display considerable departure from the predictions of the nonlinear Poisson-Boltzmann (NLPB) equation, with qualitative differences in some cases, which may be attributed to correlation effects neglected in the NLPB theory.     

Site binding as a local equilibrium. Mn2+ binding to poly(acrylic acid) as a function of the degree of ionization

Mn2+ binding to poly(acrylic acid) at different degrees of ionization, alpha, has been studied from the frequency dependence of the water protons' relaxation rates T1(-1) and T2(-1). Site binding is treated as an equilibrium with the concentration of free ions at the immediate vicinity (CIV) of the polyion. The CIV is calculated as the solution of the Poisson-Boltzmann equation at the surface of the cylindrical polyion. A single value of K is shown to fit the results at all values of alpha. The amount of site binding is higher than the total amount of condensed divalent counterions predicted for a finite polyion concentration in the presence of monovalent counterions by Manning's theory.     

Dielectric properties of polyelectrolytes. II. A theory of dielectric increment due to ion fluctuation by a matrix method

The mean square of dipole moment of a linear macromolecule which is responsible for dielectric increment of aqueous polyelectrolyte solutions is calculated by means of a matrix method in which ion binding at discrete sites and the nearest-neighbor interaction are taken into account. On the basis of the relationship between polarization of poly-ion and fluctuation of bound counterions the present theory indicates that the loosely bound ions result in larger increment and otherwise smaller increment. Also, the theory predicts that the dielectric increment has a maximum at an intermediate monovalent–divalent ion ratio when both species coexist. These results are consistent with experiments on polyacrylic acid neutralized with NaOH and Ca(OH)₂. At large contents of divalent ions the effect of chelation is also discussed.     

Specific interactions versus counterion condensation. 1. Nongelling ions/polyuronate systems

The characteristics of the interaction between nongelling divalent cations (typically Mg(2+)) and polyuronates have been explored by means of isothermal calorimetry. In particular, three polyuronates mimicking separately guluronan (polyguluronate, polyG), mannuronan (polymannuronate, polyM), and polyalternating (polyMG), the three block-components of natural alginate samples, have been treated with divalent ions, and the enthalpy of mixing was determined for different values of the [M(2+)]/[Polym](rep.unit) ratio. Despite the absence of a site-specific chemical bonding between the two, as confirmed by circular dichroism spectroscopy, a substantial deviation of the experimental enthalpy of mixing from the theoretical behavior, as predicted by the classical counterion condensation (CC) theory, was observed. Such deviation has been interpreted in terms of a "generic" nonbonding affinity of the condensed divalent counterion for the polyelectrolytes. The mathematical formalism of the CC theory was extended to include a contribution to the (reduced) free energy and enthalpy arising from the counterion affinity, g(aff,0) and h(aff,0), and allowed the parametrical calculation of the fraction of divalent counterions condensed as function of the reduced thermodynamic quantity g(aff,0). A best fit procedure of the experimental enthalpy of mixing allowed the g(aff,0) and h(aff,0) pair to be estimated for each of the different polyuronates considered, revealing differences in the three samples. In qualitative terms, the results obtained seem to suggest a notable contribution of the desolvation process (i.e., release of structured water as a consequence of the interaction between the divalent counterion and the uronate group) to the enthalpy of affinity for polyM which is counterbalanced and overcome by an ion pairing term (i.e., partial formation of ion-ion and/or ion-dipole bonds) for polyG and polyMG, respectively.     

Influence of metal ions on the ribonuclease P reaction. Distinguishing substrate binding from catalysis.

A high yield, photoactivated cross-linking reaction between a modified tRNA and RNase P RNA was used as a quantitative assay of substrate binding affinity. The cross-linking assay allows the effects of metal ions on substrate binding to be measured independently and in the absence of the pre-tRNA cleavage reaction. The results of this assay, in conjunction with the conventional cleavage assay, support the following conclusions about the nature of the RNase P RNA-tRNA binding interaction. (i) Monovalent cations act primarily to enhance enzyme-substrate binding, presumably by functioning as counterions. This enhancement can be attributed to a reduction in the tRNA off-rate. (ii) Although divalent cation is required for cleavage, the enzyme-substrate complex can form in the absence of divalent cation; the essential role of divalent cation in the reaction is thus catalytic. (iii) Ca2+ is as efficient as Mg2+ in promoting binding but supports catalysis only at a low rate.     

A commentary on the screened-Oseen, counterion-condensation formalism of polyion electrophoresis

The use of linear theory, in particular, counterion condensation (CC) theory, in describing electrophoresis of polyelectrolyte chains, is criticized on several grounds. First, there are problems with CC theory in describing the equilibrium distribution of ions around polyelectrolytes. Second, CC theory is used to treat ion relaxation in a linear theory with respect to the polyion charge despite the fact that ion relaxation arises as a consequence of nonlinear charge effects. This nonlinearity has been well established by several investigators over the last 70 years for spherical, cylindrical, and arbitrarily shaped model polyions. Third, current use of CC theory ignores the electrophoretic hindrance as well as the ion relaxation for condensed counterions and only includes such interactions for uncondensed counterions. Because most of the condensed counterions lie outside the shear surface of the polyion (in the example of DNA), the assumption of ion condensation is artificial and unphysical. Fourth, the singular solution, based on a screened Oseen tensor, currently used in the above mentioned theories is simply wrong and fails to account for the incompressibility of the solvent. The actual singular solution, which has long been available, is discussed. In conclusion, it is pointed out that numerical alternatives based on classic electrophoresis theory (J.T.G. Overbeek, Kolloid-Beih, 1943, 54:287-364) are now available.     

Cation binding to membranes: Competition between mono-, di- and trivalent cations

The binding of mono-, di- and trivalent cations to negatively charged surfaces is studied within the framework of a modified Gouy-Chapman equation. For any given combination of ions of the above valences, the existence and uniqueness of the solution for the surface potential is shown. The treatment provides the surface potential and charge density. For a system containing only monovalent and divalent ions, analytical solutions are given. When trivalent ions are also present, a procedure based on numerical integration is described. The distance dependence of the electrostatic potential for planar surfaces is given. The calculations provide the amount of cations tightly bound and the amount trapped in the double layer region. The competition between cations for binding to surfaces is elucidated.     

Concentration dependence of equivalent conductivity of polyelectrolyte

Equivalent conductivity of aqueous solutions of alternating copolymer of iso-butyl vinyl ether and maleic acid, [poly(iso BVE-co-MA)] was studied, especially its polymer-concentration dependence. Various species of counterions such as quaternary ammoniumions (NMe4+, NEt4+, NPr4+, NBu4+) and divalent ions (Ca2+, Sr2+, Ba2+) were employed besides alkali metal ions. The applicability of Manning's conductivity theory was examined for the case of univalent counterions at various degrees of neutralization (beta). A major discrepancy against the theory was observed at beta = 1.0, while a comparatively good agreement was found at beta around 0.5. This suggests that the rod-like polyion model, which is the basis of the theory, is applicable near beta = 0.5, where polyions are most expanded. The low conductivities in the case of quaternary ammonium counterions suggested the ion-binding due to hydrophobic interaction with alkyl side chains. Molecular weight dependence was not appreciably observed near beta = 0.5 similarly to usual polyelectrolytes, but it appeared slightly at beta = 1.0.     

DNA like-charge attraction and overcharging by divalent counterions in the presence of divalent co-ions

Strongly correlated electrostatics of DNA systems has drawn the interest of many groups, especially the condensation and overcharging of DNA by multivalent counterions. By adding counterions of different valencies and shapes, one can enhance or reduce DNA overcharging. In this paper, we focus on the effect of multivalent co-ions, specifically divalent co-ions such as SO\(_{4}^{2-}\). A computational experiment of DNA condensation using Monte Carlo simulation in grand canonical ensemble is carried out where the DNA system is in equilibrium with a bulk solution containing a mixture of salt of different valency of co-ions. Compared to systems with purely monovalent co-ions, the influence of divalent co-ions shows up in multiple aspects. Divalent co-ions lead to an increase of monovalent salt in the DNA condensate. Because monovalent salts mostly participate in linear screening of electrostatic interactions in the system, more monovalent salt molecules enter the condensate leads to screening out of short-range DNA–DNA like charge attraction and weaker DNA condensation free energy. The overcharging of DNA by multivalent counterions is also reduced in the presence of divalent co-ions. Strong repulsions between DNA and divalent co-ions and among divalent co-ions themselves lead to a depletion of negative ions near the DNA surface as compared to the case without divalent co-ions. At large distances, the DNA–DNA repulsive interaction is stronger in the presence of divalent co-ions, suggesting that divalent co-ions’ role is not only that of simple stronger linear screening.     

Intrinsic binding effects in mixed counterionic polyelectrolyte systems: extension of condensation theory and comparison with voltammetric data

Voltammetric speciation data for the potassium/zinc/polymethacrylate system, recently obtained for various charge densities of the polyelectrolyte (Díaz-Cruz et al., Anal. Chim. Acta, 264 (1992) 163) and for different concentrations of monovalent counterions (van den Hoop and van Leeuwen, Anal. Chim. Acta, 273 (1993) 275), are compared with theoretical predictions computed according to a new thermodynamic model developed by Paoletti et al. (Biophys. Chem., 41 (1991) 73) and recently extended by Benegas and Paoletti (in preparation). The model allows: (i) the simultaneous condensation of both monovalent and divalent counterions and (ii) can account for a certain specific affinity of the polyelectrolyte for one type of the counterion over the other. For various charge densities of the polyelectrolyte, experimentally obtained speciation data for the K/Zn/PMA system agree well with theoretical predictions by considering an extra reduced molar affinity energy of -4RT for the Zn(2+) polyelectrolyte binding. The agreement between experimental and theoretical values for the distribution of Zn(2+) ions over the free and bound state becomes less perfect for relatively high concentrations of monovalent counterions.     

Polyelectrolyte theory. I. Counterion accumulation,site-binding,and their insensitivity to polyelectrolyte shape in solutions containing finite salt concentrations

We have computed the Poisson-Boltzmann distribution of counterions around polyelectrolytes in solutions containing finite salt concentrations. The polyelectrolytes considered here are highly charged in the sense that ξ > 1, ξ being the linear charge density parameter for cylinders, which is generalized by us to other shapes. Contrary to the situation at zero salt concentration, the counterion distribution is not strongly shape dependent, being similar for cylinders or spheres which have the same superficial charge density and radius of curvature R_c. The distribution resembles that in the neighborhood of a plane with the same charge density. Three regions are distinguished. (1) In the “inner region” which extends up to a distance R_c/2ξ from the surface, the counterion distribution is essentially salt independent. The counterion concentration in the immediate vicinity of the polyelectrolyte surface (CIV) is quite high, typically 1–10M, and proportional to the square of the surface charge density, which is its main determinant. (2) An intermediate region extends out to a distance where the electrostatic potential is equal to κT/e. This distance is comparable to λ for plane and cylinder, and smaller for the sphere. (3) In the outer region, the distribution is hardly influenced by the details of the inner region, on which it cannot, therefore, give much information. Colligative properties are dependent on the distribution in the outer region and are fairly well predicted even by a rudimentary theory. The large value of the CIV implies that site binding must often be significant. It can be computed by applying the mass-action law to site-bound counterions in equilibrium with the counterions in the neighborhood, whose concentration is the CIV, the relevant equilibrium constant being that for the binding of counterions to isolated monomer sites. Because the CIV is insensitive to salt concentration, this will also be the case for site binding. With the graphs provided, one can compute the extent of sitebinding within the Poisson-Boltzmann framework. The “condensation radius,” i.e., the radius encompassing a counterionic charge 1 ? ξ^?1 around a cylinder, is found to be large. It varies with salt concentration and tends to infinity as the salt is diluted. Neither this radius nor the charge fraction 1 ? ξ^?1 of condensation theory plays any special role in the counterion distribution. The “finite-salt” results apply to salt concentrations, typically as low as 1–10 mM. This encompasses, among others, all experiments on biological polyelectrolytes.     

Dielectric properties of polyelectrolytes. III. Effect of divalent cations on dielectric increment of polyacids

Dielectric dispersion of polyacrylic acid (PAA) and polystyrene sulfonic acid (PSS) was measured in the presence of divalent cations. Effects of divalent ions were studied by neutralization with varying ratios of sodium hydroxide and divalent base concentration, addition of salts of divalent cations, and neutralization with divalent bases only. Two dispersion regions were observed in all cases, i.e., low-frequency dispersion (10²–10⁴ Hz) and high-frequency dispersion (10⁵–10⁶ Hz). The dielectric increment increases in the presence of sodium and alkaline earth metal ions together, but not with sodium and transition metal ions. This is due to the increment of low-frequency dispersion and is attributable to the fluctuation of bound counterions which is explained by our theory previously reported.¹ In the case of PAA neutralized with large fractions of divalent ions, or with divalent ions only, the increment is very small because of reduction of the fluctuation by interaction between bound ions at the neighboring sites and reduction of the effective length of polyion probably due to chelation by divalent ions. There are some differences among the effects of Mg⁺⁺, Ca⁺⁺, and Ba⁺⁺ on dielectric increment which may result from affinity or chelating ability of these ions.     

Ions and Inhibitors in the Binding Site of HIV Protease: Comparison of Monte Carlo Simulations and the Linearized Poisson-Boltzmann Theory

Proteins can be influenced strongly by the electrolyte in which they are dissolved, and we wish to model, understand, and ultimately control such ionic effects. Relatively detailed Monte Carlo (MC) ion simulations are needed to capture biologically important properties of ion channels, but a simpler treatment of ions, the linearized Poisson-Boltzmann (LPB) theory, is often used to model processes such as binding and folding, even in settings where the LPB theory is expected to be inaccurate. This study uses MC simulations to assess the reliability of the LPB theory for such a system, the constrained, anionic active site of HIV protease. We study the distributions of ions in and around the active site, as well as the energetics of displacing ions when a protease inhibitor is inserted into the active site. The LPB theory substantially underestimates the density of counterions in the active site when divalent cations are present. It also underestimates the energy cost of displacing these counterions, but the error is not consequential because the energy cost is less than k_BT, according to the MC calculations. Thus, the LPB approach will often be suitable for studying energetics, but the more detailed MC approach is critical when ionic distributions and fluxes are at issue.     

Induced coalescence of cations through low-temperature Poisson-Boltzmann calculations

The computational determination of preferred binding regions of divalent counterions to nucleic acids is either inaccurate (standard Poisson-Boltzmann approaches) or extremely time-consuming (Monte Carlo or molecular dynamics simulations). A novel "selective low-temperature" Poisson-Boltzmann method is introduced that, although approximate in nature, qualitatively accounts for ion correlation and charge-transfer effects and allows for the rapid determination of such regions through an "induced coalescence" of divalent ions. The method is illustrated here for the binding of Mg(2+) to a double-helical sequence of B-form DNA (CGCGAATTCGCG) but the technique is readily applicable to locating divalent cations in other systems such as DNA-endonuclease complexes and ribozymes.