Divalent cations and the electrostatic potential around DNA: Monte Carlo and Poisson-Boltzmann calculations

Modeling Microtubule Counterion Distributions and Conductivity Using the Poisson-Boltzmann Equation

Microtubules are highly negatively charged proteins which have been shown to behave as bio-nanowires capable of conducting ionic currents. The electrical characteristics of microtubules are highly complicated and have been the subject of previous work; however, the impact of the ionic concentration of the buffer solution on microtubule electrical properties has often been overlooked. In this work we use the non-linear Poisson Boltzmann equation, modified to account for a variable permittivity and a Stern Layer, to calculate counterion concentration profiles as a function of the ionic concentration of the buffer. We find that for low-concentration buffers ([KCl] from 10 μM to 10 mM) the counterion concentration is largely independent of the buffer's ionic concentration, but for physiological-concentration buffers ([KCl] from 100 to 500 mM) the counterion concentration varies dramatically with changes in the buffer's ionic concentration. We then calculate the conductivity of microtubule-counterion complexes, which are found to be more conductive than the buffer when the buffer's ionic concentrations is less than ≈100 mM and less conductive otherwise. These results demonstrate the importance of accounting for the ionic concentration of the buffer when analyzing microtubule electrical properties both under laboratory and physiological conditions. We conclude by calculating the basic electrical parameters of microtubules over a range of ionic buffer concentrations applicable to nanodevice and medical applications.

Ionic effects in semi-dilute biopolymer solutions: A small angle scattering study

Systematic investigations using neutron and X-ray small angle scattering in near-physiological salt solutions were made to reveal the effect of polymer concentration, pH, and calcium ion concentration on the structure of semi-dilute solutions of four model biopolymers [polyaspartic acid, DNA, chondroitin sulfate, and hyaluronic acid (HA)] representing typical backbone structures. In the low q range (<0.01 Å-1), the scattering response I(q) is dominated by scattering from large clusters. In the intermediate q range, I(q) varies approximately as q -1, exposing the linear nature of the scatterers. In these polyelectrolyte solutions, the correlation length L displays a power law dependence on the polymer concentration c that resembles that of neutral polymer solutions. L increases with increasing calcium chloride concentration and with decreasing pH. The effect of the different divalent cations, Ba, Mg, Ca, Sr, and Mn, on the structure of DNA solutions is practically identical. However, in mixed salt conditions at the same ionic strength, the combined effect of mono- and divalent counter-ions on the structure of the polymer solutions deviates significantly from additivity. Anomalous small angle X-ray scattering observations on both DNA and HA solutions reveal that the divalent strontium counter-ions form a tight sheath around the polymer chain. The shape of the divalent ion cloud is similar in these two systems.

On the energy components governing molecular recognition in the framework of continuum approaches

Molecular recognition is a process that brings together several biological macromolecules to form a complex and one of the most important characteristics of the process is the binding free energy. Various approaches exist to model the binding free energy, provided the knowledge of the 3D structures of bound and unbound molecules. Among them, continuum approaches are quite appealing due to their computational efficiency while at the same time providing predictions with reasonable accuracy. Here we review recent developments in the field emphasizing on the importance of adopting adequate description of physical processes taking place upon the binding. In particular, we focus on the efforts aiming at capturing some of the atomistic details of the binding phenomena into the continuum framework. When possible, the energy components are reviewed independently of each other. However, it is pointed out that rigorous approaches should consider all energy contributions on the same footage. The two major schemes for utilizing the individual energy components to predict binding affinity are outlined as well.

Effect of sodium ions on RNA duplex stability

The standard sodium concentration for RNA optical melting experiments is 1.021 M. Algorithms that predict Tm, ΔG°37, and secondary structure from sequence generally rely on parameters derived from optical melting experiments performed in 1.021 M sodium. Physiological monovalent cation concentrations are much lower than 1.021 M. In fact, many molecular biology techniques require buffers containing monovalent cation concentrations other than 1.021 M. Predictions based on the 1.021 M Na(+) parameters may not be accurate when the monovalent cation concentration is not 1.021 M. Here, we report thermodynamic data from optical melting experiments for a set of 18 RNA duplexes, each melted over a wide range of sodium ion concentrations (71, 121, 221, and 621 mM). Using these data and previously published data for the same sequences melted in 1.021 M Na(+), we report Tm and ΔG°37 correction factors to scale the standard 1.021 M Na(+) RNA parameters to other sodium ion concentrations. The recommended Tm correction factor predicts the melting temperature within 0.7 °C, and the recommended ΔG°37 correction factor predicts the free energy within 0.14 kcal/mol. These correction factors can be incorporated into prediction algorithms that predict RNA secondary structure from sequence and provide Tm and ΔG°37 values for RNA duplexes.

Bullied no more: when and how DNA shoves proteins around

The predominant protein-centric perspective in protein-DNA-binding studies assumes that the protein drives the interaction. Research focuses on protein structural motifs, electrostatic surfaces and contact potentials, while DNA is often ignored as a passive polymer to be manipulated. Recent studies of DNA topology, the supercoiling, knotting, and linking of the helices, have shown that DNA has the capability to be an active participant in its transactions. DNA topology-induced structural and geometric changes can drive, or at least strongly influence, the interactions between protein and DNA. Deformations of the B-form structure arise from both the considerable elastic energy arising from supercoiling and from the electrostatic energy. Here, we discuss how these energies are harnessed for topology-driven, sequence-specific deformations that can allow DNA to direct its own metabolism.

Electrostatics of nucleic acid folding under conformational constraint

RNA folding is enabled by interactions between the nucleic acid and its ion atmosphere, the mobile sheath of aqueous ions that surrounds and stabilizes it. Understanding the ion atmosphere requires the interplay of experiment and theory. However, even an apparently simple experiment to probe the ion atmosphere, measuring the dependence of DNA duplex stability upon ion concentration and identity, suffers from substantial complexity, because the unfolded ensemble contains many conformational states that are difficult to treat accurately with theory. To minimize this limitation, we measured the unfolding equilibrium of a DNA hairpin using a single-molecule optical trapping assay, in which the unfolded state is constrained to a limited set of elongated conformations. The unfolding free energy increased linearly with the logarithm of monovalent cation concentration for several cations, such that smaller cations tended to favor the folded state. Mg(2+) stabilized the hairpin much more effectively at low concentrations than did any of the monovalent cations. Poisson-Boltzmann theory captured trends in hairpin stability measured for the monovalent cation titrations with reasonable accuracy, but failed to do so for the Mg(2+) titrations. This finding is consistent with previous work, suggesting that Poisson-Boltzmann and other mean-field theories fail for higher valency cations where ion-ion correlation effects may become significant. The high-resolution data herein, because of the straightforward nature of both the folded and the unfolded states, should serve as benchmarks for the development of more accurate electrostatic theories that will be needed for a more quantitative and predictive understanding of nucleic acid folding.

Counterion condensation and shape within Poisson-Boltzmann theory

An analytical approximation to the nonlinear Poisson-Boltzmann (PB) equation is applied to charged macromolecules that possess one-dimensional symmetry and can be modeled by a plane, infinite cylinder, or sphere. A functional substitution allows the nonlinear PB equation subject to linear boundary conditions to be transformed into an approximate linear (Debye-Hückel-type) equation subject to nonlinear boundary conditions. A simple analytical result for the surface potential of such polyelectrolytes follows, leading to expressions for the amount of condensed (or renormalized) charge and the electrostatic Helmholtz energy for polyelectrolytes. Analytical high-charge/low-salt and low-charge/high-salt limits are shown to be similar to results obtained by others based on PB or counterion condensation theory. Several important general observations concerning polyelectrolytes treated within the context of PB theory can be made including: (1) all charged surfaces display some counterion condensation for finite electrolyte concentration, (2) the effect of surface geometry is described primarily by the sum of the Debye constant and the mean curvature of the surface, (3) two surfaces with the same surface charge density and mean curvature condense approximately identical fractions of counterions, (4) the amount of condensation is not determined by a predefined "condensation distance" although such a distance can be determined uniquely from it, and (5) substantial condensation occurs if the Debye constant of the electrolyte is much less than the mean curvature of a highly charged polyelectrolyte.

Critical assessment of nucleic acid electrostatics via experimental and computational investigation of an unfolded state ensemble

Electrostatic forces, acting between helices and modulated by the presence of the ion atmosphere, are key determinants in the energetic balance that governs RNA folding. Previous studies have employed Poisson-Boltzmann (PB) theory to compute the energetic contribution of these forces in RNA folding. However, the complex interaction of these electrostatic forces with RNA features such as tertiary contact formation, specific ion-binding, and complex interhelical junctions present in prior studies precluded a rigorous evaluation of PB theory, especially in physiologically important Mg(2+) solutions. To critically assess PB theory, we developed a model system that isolates these electrostatic forces. The model system, composed of two DNA duplexes tethered by a polyethylene glycol junction, is an analog for the unfolded state of canonical helix-junction-helix motifs found in virtually all structured RNAs. This model system lacks the complicating features that have precluded a critical assessment of PB in prior studies, ensuring that interhelical electrostatic forces dominate the behavior of the system. The system's simplicity allows PB predictions to be directly compared with small-angle X-ray scattering experiments over a range of monovalent and divalent ion concentrations. These comparisons indicate that PB is a reasonable description of the underlying electrostatic energies for monovalent ions, but large deviations are observed for divalent ions. The validation of PB for monovalent solutions allows analysis of the change in the conformational ensemble of this simple motif as salt concentration is changed. Addition of ions allows the motif to sample more compact microstates, increasing its conformational entropy. The increase of conformational entropy presents an additional barrier to folding by stabilizing the unfolded state. Neglecting this effect will adversely impact the accuracy of folding analyses and models.

Mono- and trivalent ions around DNA: a small-angle scattering study of competition and interactions

The presence of small numbers of multivalent ions in DNA-containing solutions results in strong attractive forces between DNA strands. Despite the biological importance of this interaction, e.g., DNA condensation, its physical origin remains elusive. We carried out a series of experiments to probe interactions between short DNA strands as small numbers of trivalent ions are included in a solution containing DNA and monovalent ions. Using resonant (anomalous) and nonresonant small angle x-ray scattering, we coordinated measurements of the number and distribution of each ion species around the DNA with the onset of attractive forces between DNA strands. DNA-DNA interactions occur as the number of trivalent ions increases. Surprisingly good agreement is found between data and size-corrected numerical Poisson-Boltzmann predictions of ion competition for non- and weakly interacting DNAs. We also obtained an estimate for the minimum number of trivalent ions needed to initiate DNA-DNA attraction.

Predicting stability of DNA duplexes in solutions containing magnesium and monovalent cations

Accurate predictions of DNA stability in physiological and enzyme buffers are important for the design of many biological and biochemical assays. We therefore investigated the effects of magnesium, potassium, sodium, Tris ions, and deoxynucleoside triphosphates on melting profiles of duplex DNA oligomers and collected large melting data sets. An empirical correction function was developed that predicts melting temperatures, transition enthalpies, entropies, and free energies in buffers containing magnesium and monovalent cations. The new correction function significantly improves the accuracy of predictions and accounts for ion concentration, G-C base pair content, and length of the oligonucleotides. The competitive effects of potassium and magnesium ions were characterized. If the concentration ratio of [Mg (2+)] (0.5)/[Mon (+)] is less than 0.22 M (-1/2), monovalent ions (K (+), Na (+)) are dominant. Effects of magnesium ions dominate and determine duplex stability at higher ratios. Typical reaction conditions for PCR and DNA sequencing (1.5-5 mM magnesium and 20-100 mM monovalent cations) fall within this range. Conditions were identified where monovalent and divalent cations compete and their stability effects are more complex. When duplexes denature, some of the Mg (2+) ions associated with the DNA are released. The number of released magnesium ions per phosphate charge is sequence dependent and decreases surprisingly with increasing oligonucleotide length.

Continuum molecular electrostatics, salt effects, and counterion binding--a review of the Poisson-Boltzmann theory and its modifications

This work is a review of the Poisson-Boltzmann (PB) continuum electrostatics theory and its modifications, with a focus on salt effects and counterion binding. The PB model is one of the mesoscopic theories that describes the electrostatic potential and equilibrium distribution of mobile ions around molecules in solution. It serves as a tool to characterize electrostatic properties of molecules, counterion association, electrostatic contributions to solvation, and molecular binding free energies. We focus on general formulations which can be applied to large molecules of arbitrary shape in all-atomic representation, including highly charged biomolecules such as nucleic acids. These molecules present a challenge for theoretical description, because the conventional PB model may become insufficient in those cases. We discuss the conventional PB equation, the corresponding functionals of the electrostatic free energy, including a connection to DFT, simple empirical extensions to this model accounting for finite size of ions, the modified PB theory including ionic correlations and fluctuations, the cell model, and supplementary methods allowing to incorporate site-bound ions in the PB calculations.

Quantitative and comprehensive decomposition of the ion atmosphere around nucleic acids

The ion atmosphere around nucleic acids critically affects biological and physical processes such as chromosome packing, RNA folding, and molecular recognition. However, the dynamic nature of the ion atmosphere renders it difficult to characterize. The basic thermodynamic description of this atmosphere, a full accounting of the type and number of associated ions, has remained elusive. Here we provide the first complete accounting of the ion atmosphere, using buffer equilibration and atomic emission spectroscopy (BE-AES) to accurately quantitate the cation association and anion depletion. We have examined the influence of ion size and charge on ion occupancy around simple, well-defined DNA molecules. The relative affinity of monovalent and divalent cations correlates inversely with their size. Divalent cations associate preferentially over monovalent cations; e.g., with Na+ in 4-fold excess of Mg2+ (20 vs 5 mM), the ion atmosphere nevertheless has 3-fold more Mg2+ than Na+. Further, the dicationic polyamine putrescine2+ does not compete effectively for association relative to divalent metal ions, presumably because of its lower charge density. These and other BE-AES results can be used to evaluate and guide the improvement of electrostatic treatments. As a first step, we compare the BE-AES results to predictions from the widely used nonlinear Poisson Boltzmann (NLPB) theory and assess the applicability and precision of this theory. In the future, BE-AES in conjunction with improved theoretical models, can be applied to complex binding and folding equilibria of nucleic acids and their complexes, to parse the electrostatic contribution from the overall thermodynamics of important biological processes.

Molecular dynamics studies of ion distributions for DNA duplexes and DNA clusters: salt effects and connection to DNA melting

We present extensive molecular dynamics simulations of the ion distributions for DNA duplexes and DNA clusters using the Amber force field with implicit water. The distribution of ions and the electrostatic energy of ions around an isolated DNA duplex and clusters of DNA duplexes in different salt (NaCl) concentrations over the range 0.2-1.0 mol/L are determined on the basis of the simulation results. Using the electrostatic energy profile, we determine a local net charge fraction phi, which is found to increase with increasing of salt concentration. For DNA clusters containing two DNA duplexes (DNA pair) or four DNA duplexes, phi increases as the distance between the duplexes decreases. Combining this result with experimental results for the dependence of the DNA melting temperature on bulk salt concentration, we conclude that for a pair of DNA duplexes the melting temperature increases by 5-10 K for interaxis separations of 25-40 A. For a cluster of four DNA duplexes, an even larger melting temperature increase should occur. We argue that this melting temperature increase in dense DNA clusters is responsible for the cooperative melting mechanism in DNA-linked nanoparticle aggregates and DNA-linked polymer aggregates.

Preferential interaction between DNA and small ions in mixed-size counterion systems: Monte Carlo simulation and density functional study

Competitive binding between counterions around DNA molecule is characterized using the preferential interaction coefficient of individual ion in single and mixed electrolyte solutions. The canonical Monte Carlo (MC) simulation, nonlinear Poisson-Boltzmann (PB) equation, and density functional theory (DFT) proposed in our previous work [Wang, Yu, Gao, and Luo, J. Chem. Phys. 123, 234904 (2005)] are utilized to calculate the preferential interaction coefficients. The MC simulations and theoretical results show that for single electrolyte around DNA, the preferential interaction coefficient of electrolyte decreases as the cation size is increased, indicating that the larger cation has less accumulation ability in the vicinity of DNA. For the mixed electrolyte solution, it is found that cation diameter has a significant effect on the competitive ability while anion diameter has a negligible effect. It proves that the preferential interaction coefficients of all ions decrease as the total ionic concentration is increased. The DFT generally has better performance than the PB equation does when compared to the MC simulation data. The DFT behaves quite well for the real ionic solutions such as the KCl-NaCl-H2O, NaCl-CaCl2-H2O, and CaCl2-MgCl2-H2O systems.

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

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

Density-functional theory and Monte Carlo simulation study on the electric double layer around DNA in mixed-size counterion systems

A density-functional approach and canonical Monte Carlo simulations are presented for describing the ionic microscopic structure around the DNA molecule immersed in mixed-size counterion solutions. In the density-functional approach, the hard-sphere contribution to the Helmholtz energy functional is obtained from the modified fundamental measure theory [Y.-X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)], and the electrostatic contribution is evaluated through a quadratic functional Taylor expansion. The new theory is suitable to the systems containing ions of arbitrary sizes and valences. In the established canonical Monte Carlo simulation, an iterative self-consistent method is used to evaluate the long-range energy, and another iterative algorithm is adopted to obtain desired bulk ionic concentrations. The ion distributions from the density-functional theory (DFT) are in good agreement with those from the corresponding Monte Carlo (MC) simulations. It is found that the ratio of the bulk concentrations of two species of counterions (cations) makes significant contribution to the ion distributions in the vicinity of DNA. Comparisons with the electrostatic potential profiles from the MC simulations show that the accuracy of the DFT becomes low when a small divalent cation exists. Both the DFT and MC simulation results illustrate that the electrostatic potential at the surface of DNA increases as the anion diameter or the total cation concentration is increased and decreases as the diameter of one cation species is increased. The calculation of electrostatic potential using real ion diameters shows that the accuracy of DFT predictions for divalent ions is also acceptable.

Calculation of the binding free energy for magnesium-RNA interactions

The nature of the interaction between nucleic acids and divalent ions in solution is complex. It includes long-range electrostatic and short-range nonelectrostatic forces. Water molecules can be in an inner coordination shell that intervenes between the ion and its binding site. This work describes a method for calculating the binding free energy and applies it to a specific Mg-RNA system in the presence of monovalent salt. The approach combines high-level ab initio theory with Poisson-Boltzmann calculations and provides an accurate description of all terms of the binding free energy for magnesium ions located at the RNA surface (including nonelectrostatic interactions). Some alternative macroscopic approaches are also discussed.

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.

"Like-charge attraction" between anionic polyelectrolytes: molecular dynamics simulations

"Like-charge attraction" is a phenomenon found in many biological systems containing DNA or proteins, as well as in polyelectrolyte systems of industrial importance. "Like-charge attraction" between polyanions is observed in the presence of mobile multivalent cations. At a certain limiting concentration of cations, the negatively charged macroions cease to repel each other and even an attractive force between the anions is found. With classical molecular dynamics simulations it is possible to elucidate the processes that govern the attractive behavior with atomistic resolution. As an industrially relevant example we study the interaction of negatively charged carboxylate groups of sodium polyacrylate molecules with divalent cationic Ca2+ counterions. Here we show that Ca2+ ions initially associate with single chains of polyacrylates and strongly influence sodium ion distribution; shielded polyanions approach each other and eventually "stick" together (precipitate), contrary to the assumption that precipitation is initially induced by intermolecular Ca2+ bridging.

Counterion distribution around DNA probed by solution X-ray scattering

Counterion atmospheres condensed onto charged biopolymers strongly affect their physical properties and biological functions, but have been difficult to quantify experimentally. Here, monovalent and divalent counterion atmospheres around DNA double helices in solution are probed using small-angle x-ray scattering techniques. Modulation of the ion scattering factors by anomalous (resonant) x-ray scattering and by interchanging ion identities yields direct measurements of the scattering signal due to the spatial correlation of surrounding ions to the DNA. The quality of the data permit, for the first time, quantitative tests of extended counterion distributions calculated from atomic-scale models of biologically relevant molecules.

The wall and multivalent counterion effects on the electrostatic force between like-charged spherical particles confined in a charged pore

The effect of wall confinement (wall charge and wall-sphere separation distance) on the electrostatic force between two charged spheres confined in a long charged pore in symmetric and asymmetric electrolytes have been quantified by solving the nonlinear Poisson-Boltzmann equation (PBE), using adaptive finite elements combined with error minimization techniques. The computed force indicated the strong effect of the wall potential on the reduction of the repulsive force for all type of electrolytes. The influence of the wall effect was reduced when the valence of the electrolyte was increased. A significant reduction in the repulsive force between the two spheres was also observed when the distance between the pore wall and the sphere surface was reduced. A smaller long-range repulsive interaction was observed between spheres when the solutions contained multivalent counterions as compared with a monovalent solution. However, at short ranges of separation distances multivalent counterions increase the electrostatic repulsive force between the spheres. The effect of the dimensionless radius of the spheres on the electrostatic force between them has been determined and a significant reduction observed as the dimensionless radius was reduced.

The linkage between magnesium binding and RNA folding

Understanding the linkage between Mg(2+) binding and RNA folding requires a proper theoretical model describing the energetics of Mg(2+) binding to the folded and unfolded states of RNA. Our current understanding of Mg(2+) binding to these different RNA states derives from empirical thermodynamic models that depend on a number of unjustified assumptions. We present a rigorous theoretical model describing the linkage between RNA folding and magnesium ion binding. In this model, based on the non-linear Poisson-Boltzmann (NLPB) equation, the stabilization of RNA by Mg(2+) arises from two distinct binding modes, diffuse binding and site binding. Diffusely bound Mg(2+) are described as an ensemble of hydrated ions that are attracted to the negative charge of the RNA. Site-bound Mg(2+) are partially desolvated ions that are attracted to electronegative pockets on the RNA surface. We explore two systems, yeast tRNA(Phe) and a 58-nucleotide rRNA fragment, with different Mg(2+) binding properties. The NLPB equation accurately describes both the stoichiometric and energetic linkage between Mg(2+) binding and RNA folding for both of these systems without requiring any fitted parameters in the calculation. Moreover, the NLPB model presents a well-defined physical description of how Mg(2+) binding helps fold an RNA. For both of the molecules studied here, the relevant unfolded state is a disordered intermediate state (I) that contains stable helical secondary structure without any tertiary contacts. Diffusely bound Mg(2+) interact with these secondary structure elements to stabilize the I state. The secondary structural elements of the I state fold into a compact, native tertiary structure (the N state). Diffuse binding plays a dominant role in stabilizing the N state for both RNAs studied. However, for the rRNA fragment, site-binding to a location with extraordinarily high electrostatic potential is also coupled to folding. Our results suggest that much experimental data measuring the linkage between Mg(2+) binding and RNA folding must be reinterpreted.

Effective interaction between helical biomolecules

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

The interpretation of Mg(2+) binding isotherms for nucleic acids using Poisson-Boltzmann theory

Magnesium ions play a crucial role in the structural integrity and biological activity of nucleic acids. Experimental thermodynamic descriptions of Mg(2+) interactions with nucleic acids in solution have generally relied on the analyses of binding polynomials to estimate the energetic contributions of diffuse and site-bound ions. However, since ion binding is dominated by long-range electrostatic forces, such models provide only a phenomenological description of the experimental Mg(2+) binding data and provide little insight into the actual mechanism of the binding equilibria. Here, we present a rigorous theoretical framework based on the non-linear Poisson-Boltzmann (NLPB) equation for understanding diffuse ion interactions that can be used to interpret experimental Mg(2+) binding isotherms. As intuitively expected, in the NLPB model binding is simply the total accumulation of the ion around the nucleic acid. Comparing the experimental data to the calculated curves shows that the NLPB equation provides a remarkably accurate description of Mg(2+) binding to linear polynucleotides like DNA and poly(A x U) without any fitted parameters. In particular, the NLPB model explains two general features of magnesium binding; the strong dependence on univalent salt concentration, and its substantial anticooperativity. Each of these effects can be explained by changes in the Mg(2+) distribution around the polyion under different solution conditions. In order to more fully understand these different aspects of magnesium binding, the free energy of Mg(2+) binding, DeltaGMg, is calculated and partitioned into several salt-dependent contributions: the change in the electrostatic interaction free energy of the charges, DeltaDeltaGE.D (including Mg(2+)-phosphate, Mg(2+)-Mg(2+), Mg(2+)-Na(+), Na(+)-Na(+), Na(+)-phosphate interactions, and similar contributions for Cl(-)) and the cratic free energies of (re)organizing the MgCl2 and NaCl atmospheres, DeltaG(Mg)org and DeltaDeltaG(Na)org, respectively. For the systems studied here, DeltaGMg is strongly influenced by entropic free energy changes in the distributions of both NaCl and MgCl2, DeltaG(Mg)org and DeltaDeltaG(Na)org. From this analysis, we also raise the possibility that coions added with the magnesium salt might play an important role in the overall stability of nucleic acids under some conditions.