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.     

DNA and its counterions: a molecular dynamics study

The behaviour of mobile counterions, Na⁺ and K⁺, was analysed around a B-DNA double helix with the sequence CCATGCGCTGAC in aqueous solution during two 50 ns long molecular dynamics trajectories. The movement of both monovalent ions remains diffusive in the presence of DNA. Ions sample the complete space available during the simulation time, although individual ions sample only about one-third of the simulation box. Ions preferentially sample electronegative sites around DNA, but direct binding to DNA bases remains a rather rare event, with highest site occupancy values of <13%. The location of direct binding sites depends greatly on the nature of the counterion. While Na⁺ binding in both grooves is strongly sequence-dependent with the preferred binding site in the minor groove, K⁺ mainly visits the major groove and binds close to the centre of the oligomer. The electrostatic potential of an average DNA structure therefore cannot account for the ability of a site to bind a given cation; other factors must also play a role. An extensive analysis of the influence of counterions on DNA conformation showed no evidence of minor groove narrowing upon ion binding. A significant difference between the conformations of the double helix in the different simulations can be attributed to extensive α/γ transitions in the phosphate backbone during the simulation with Na⁺. These transitions, with lifetimes over tens of nanoseconds, however, appear to be correlated with ion binding to phosphates. The ion-specific conformational properties of DNA, hitherto largely overlooked, may play an important role in DNA recognition and binding.     

Adsorption of divalent cations on DNA

The distribution of divalent ions in semidilute solutions of high-molecular-mass DNA containing both sodium chloride and strontium chloride in near-physiological conditions is studied by small-angle x-ray scattering and by small-angle neutron scattering. Both small-angle neutron scattering and small-angle x-ray scattering reveal a continuous increase in the scattering intensity at low q with increasing divalent ion concentration, while at high q the scattering curves converge. The best fit to the data is found for a configuration in which DNA strands of cross-sectional radius 10 angstroms are surrounded by a counterion sheath of outer radius approximately 13.8 angstroms, independent of the strontium chloride concentration. When the strontium chloride is replaced by calcium chloride, similar results are obtained, but the thickness of the sheath increases when the divalent salt concentration decreases. These results correspond in both cases to partial localization of the counterions within a layer that is thinner than the effective Debye screening length.     

The release of monovalent counterions by addition of divalent counterions in coulombic interaction system

The additivity rule of counterion activity or osmotic pressure in rodlike polyelectrolyte solutions has been discussed on the basis of the Fokker-Planck and Poisson equations in relation to the fluctuation of counterion distribution. This new theory has concluded that the additivity rule of counterion activity is less applicable than that of osmotic pressure due to the electric expansion force acting on the free-volume surface resulting from the fluctuation of counterion distribution. The theory has introduced an approximate relation between the counterion activities in the mixture solution of divalent and monovalent counterions, such that Deltaa+ = DeltaC++ - Deltaa++, in which Deltaa+ represents the increase of activity of monovalent counter-ions resulting from the addition of divalent counterionsDeltaC++, (in molar) to the solution, and Deltaa++ means the increase of the divalent counterion activity (in molar) in this process. This relation has been experimentally examined for Na-PSS solutions in the process of Cu2+ ion addition by the use of Na+ and Cu2+ sensitive electrodes, and it has been turned out that the relation is established in the low charge state of polyion.     

Differential stability of DNA crossovers in solution mediated by divalent cations

The assembly of DNA duplexes into higher-order structures plays a major role in many vital cellular functions such as recombination, chromatin packaging and gene regulation. However, little is currently known about the molecular structure and stability of direct DNA–DNA interactions that are required for such functions. In nature, DNA helices minimize electrostatic repulsion between double helices in several ways. Within crystals, B-DNA forms either right-handed crossovers by groove–backbone interaction or left-handed crossovers by groove–groove juxtaposition. We evaluated the stability of such crossovers at various ionic concentrations using large-scale atomistic molecular dynamics simulations. Our results show that right-handed DNA crossovers are thermodynamically stable in solution in the presence of divalent cations. Attractive forces at short-range stabilize such crossover structures with inter-axial separation of helices less than 20 Å. Right-handed crossovers, however, dissociate swiftly in the presence of monovalent ions only. Surprisingly, left-handed crossovers, assembled by sequence-independent juxtaposition of the helices, appear unstable even at the highest concentration of Mg²⁺studied here. Our study provides new molecular insights into chiral association of DNA duplexes and highlights the unique role divalent cations play in differential stabilization of crossover structures. These results may serve as a rational basis to understand the role DNA crossovers play in biological processes.     

Phosphatidic acid domains in membranes: effect of divalent counterions

Phosphatidic acid (PA) is emerging as a key phospholipid in a wide range of biological processes such as signal transduction, secretion, or membrane fusion. In most cases, the biological functionality of PA is associated with the presence of micromolar to millimolar calcium concentrations. It has been argued that PA can create defects in the packing of lipids in membranes due to lateral phase separation by divalent ions, which in turn aggregate proteins with high affinity for PA. In this article, we present a detailed investigation of the properties of PA domains in the presence of divalent ions by a combination of molecular dynamics simulations and theoretical methods. Our results show that PA is extremely effective in binding divalent ions through its oxygen atoms, with a broad distribution of binding constants and exhibiting the phenomenon of charge inversion (a total number of bound counterion charges that exceeds the negative PA charge). We predict that a PA-rich domain undergoes a drastic reorganization when divalent cations reach micromolar concentrations (i.e., typical physiological conditions), as PA lipids become doubly charged by releasing their protons. We also present a detailed investigation of the properties of interfacial water, which determine the binding of proteins or other molecules. We conclude with a discussion of the implications of our results in the context of recent experimental studies in model systems and in real cells.     

Hydration of counterions interacting with DNA double helix: a molecular dynamics study

In the present work, molecular dynamics simulations have been carried out to study the dependence of counterion distribution around the DNA double helix on the character of ion hydration. The simulated systems consisted of DNA fragment d(CGCGAATTCGCG) in water solution with the counterions Na⁺, K⁺, Cs⁺ or Mg²⁺. The characteristic binding sites of the counterions with DNA and the changes in their hydration shell have been determined. The results show that due to the interaction with DNA at least two hydration shells of the counterions undergo changes. The first hydration shell of Na⁺, K⁺, Cs⁺, and Mg²⁺ counterions in the bulk consists of six, seven, ten, and six water molecules, respectively, while the second one has several times higher values. The Mg²⁺ and Na⁺ counterions, constraining water molecules of the first hydration shell, mostly form with DNA water-mediated contacts. In this case the coordination numbers of the first hydration shell do not change, while the coordination numbers of the second one decrease about twofold. The Cs⁺ and K⁺ counterions that do not constrain surrounding water molecules may be easily dehydrated, and when interacting with DNA their first hydration shell may be decreased by three and five water molecules, respectively. Due to the dehydration effect, these counterions can squeeze through the hydration shell of DNA to the bottom of the double helix grooves. The character of ion hydration establishes the correlation between the coordination numbers of the first and the second hydration shells.

Open image in new window

Graphical Abstract Hydration of counterions interacting with DNA double helix

     

The release of monovalent counterions by addition of divalent counterions in coulombic interaction system of polyions (II)

In this work, the increase in activity. Deltaa+ of a monovalent counterion (Na+) resulting from the addition or a divalent ion (Cu2+) was experimentally investigated in sulfate or sulfonate polyelectrolyte solutions. It was found that in the case of polyvinylsulfonate, a+ was expressed by Deltaa+ = Delta(C++ -a++). where C(++) and a(++) represent the molar concentration and activity of added Cu2+ respectively, while in the case of dextran sulfate, the expression or Deltaa+ became Deltaa+ =2Delta(C++-a++). Considering that interactions between these cations and sulfate or sulfonale groups are Durely electrostatic, these results were analyzed in relation to the conformation of the polyelectrolytes.     

Transmission of viruses to mosquito larvae mediated by divalent cations

The two major groups of pathogenic viruses in mosquitoes are the occluded viruses, represented by baculoviruses and cypoviruses, and the non-occluded viruses, represented by the densoviruses and the iridoviruses. Baculoviruses, densoviruses, and iridoviruses are DNA viruses, while cypoviruses are the major group of RNA viruses reported from mosquitoes. Research on mosquito pathogenic viruses has been limited, in part, due to the inability to effectively transmit them to the larval mosquito host. Recently, there have been tremendous advancements in the ability to transmit mosquito baculoviruses and cypoviruses with the finding that transmission is mediated by divalent cations. Oral transmission of both baculoviruses and cypoviruses to mosquito larvae is enhanced by magnesium and inhibited by calcium ions. The current status of transmission for each of the major groups is reviewed with emphasis on the common role of divalent cations in transmission of the distantly related baculoviruses and cypoviruses.     

Mechanistic clues to the mutagenicity of alkylated DNA bases: a theoretical study

Experiment indicates that the N7-guanine site in DNA is not "promutagenic" (mutation-inducing) on alkylation, while the O6-guanine and O4-thymine sites are so. These differences in nucleic acid template activity are attributed to corresponding differences in acidity of the Watson-Crick hydrogen bonding protons. Mechanistic indicators for ease of Watson-Crick proton loss are calculated using molecular orbital theory for DNA bases alkylated at the N7-guanine, O6-guanine and O4-thymine sites. Their values point to a definite favouring of the proton loss for the O-alkylated bases compared to the N7-alkylguanines. This may suggest the possibility that, at biological pH, the O-alkylated bases deprotonate readily while the N7-alkylguanines do not, thus accounting for observed differences in promutagenicity and nucleic acid template activity.     

The release of monovalent counterions by addition of divalent countemons to aqueous solutions of maleic acid copolymer

Divalent cation binding and the release of monovalent cations accompanying the cation binding were experimentally studied by ion-selective electrode methods in aqueous solutions of copolymer of maleic acid and ethyl vinyl ether. It was found that in the process of Ca2+ addition, all the Ca2+ added was bound to polyions and the initially condensed Na+ was released in proportion to the concentration of the added Ca2+ up to the critical concentration of added Ca2+ at which the condensation of Ca2+ ceases. Values of the structural charge density parameter xi(s), were determined from the end-points of condensation of Ca2+. The process of Na+ release by adding Ca2+ was analyzed on the basis of the counterion condensation theory by using these xi(s) values. In addition, the relationship between the activity coefficient gamma-- of Ca2+ and degree of neutralization alpha in salt-free solutions was obtained from the Manning theory. Agreement between the calculated and experimental values was excellent in both cases.     

Nonlinear mechanical property of tracheal cartilage: a theoretical and experimental study

BACKGROUND: Despite being the stiffest airway of the bronchial tree, the trachea undergoes significant deformation due to intrathoracic pressure during breathing. The mechanical properties of the trachea affect the flow in the airway and may contribute to the biological function of the lung. METHOD: A Fung-type strain energy density function was used to investigate the nonlinear mechanical behavior of tracheal cartilage. A bending test on pig tracheal cartilage was performed and a mathematical model for analyzing the deformation of tracheal cartilage was developed. The constants included in the strain energy density function were determined by fitting the experimental data. RESULT: The experimental data show that tracheal cartilage is a nonlinear material displaying higher strength in compression than in tension. When the compression forces varied from -0.02 to -0.03N and from -0.03 to -0.04N, the deformation ratios were 11.03+/-2.18% and 7.27+/-1.59%, respectively. Both were much smaller than the deformation ratios (20.01+/-4.49%) under tension forces of 0.02 to 0.01N. The Fung-type strain energy density function can capture this nonlinear behavior very well, whilst the linear stress-strain relation cannot. It underestimates the stability of trachea by exaggerating the displacement in compression. This study may improve our understanding of the nonlinear behavior of tracheal cartilage and it may be useful for the future study on tracheal collapse behavior under physiological and pathological conditions.     

Denaturation maps of DNA: experimental and theoretical maps of phiX174 DNA.

A formaldehyde denaturation map of the replicative form of phiX174 DNA is obtained. The RFI DNA was converted into a linear state by restriction endonuclease pst I which introduces into this DNA a single double-stranded break. The map has four clear-cut peaks. Their positions excellently correlate with the peak positions on the map of equilibrium denaturation theoretically obtained earlier from the known nucleotide sequence of phiX174 DNA. The sequence is also used for a calculation of the maps of smoothed AT-content. The maxima on these maps correlate well with the peaks on the denaturation maps. To reveal the causes of a good correlation between the experimental formaldehyde and theoretical equilibrium denaturation maps, the theoretical formaldehyde denaturation maps are calculated for different conditions (temperature, formaldehyde concentration) using the detailed theory of DNA interaction with formaldehyde developed earlier.     

Left-handed Z form in superhelical DNA: a theoretical study

This is a comprehensive statistical mechanical treatment of the Z form formation in purinepyrimidine stretches of different length inserted into superhelical DNA. The B-Z transition for short inserts is shown to follow the "all-or-none" principle. Over some critical value of the insert length n, the B-Z transition in the insert proceeds in two stages. The flipping of m base pairs into the Z form is followed by a gradual growth of the Z-form stretch until it occupies the whole insert. By fitting the theoretical transition curves to experimental ones the fundamental thermodynamic parameters of the B-Z transition have been determined: the B-Z junction energy Fj = 4-5kcal.mol-1 and the free energy change delta FBZ = 0.5-7.0 kcal.mol-1 under standard salt conditions. Calculations show that the B-Z transition in short purinepyrimidine inserts may be seriously affected by cruciform formation in the carrier DNA.     

The stereochemistry of a four-way DNA junction: a theoretical study.

The stereochemical conformation of the four-way helical junction in DNA (the Holliday junction; the postulated central intermediate of genetic recombination) has been analysed, using molecular mechanical computer modelling. A version of the AMBER program package was employed, that had been modified to include the influence of counterions and a global optimisation procedure. Starting from an extended planar structure, the conformation was varied in order to minimise the energy, and we discuss three structures obtained by this procedure. One structure is closely related to a square-planar cross, in which there is no stacking interaction between the four double helical stems. This structure is probably closely similar to that observed experimentally in the absence of cations. The remaining two structures are based on related, yet distinct, conformations, in which there is pairwise coaxial stacking of neighbouring stems. In these structures, the four DNA stems adopt the form of two quasi-continuous helices, in which base stacking is very similar to that found in standard B-DNA geometry. The two stacked helices so formed are not aligned parallel to each other, but subtend an angle of approximately 60 degrees. The strands that exchange between one stacked helix and the other are disposed about the smaller angle of the cross (i.e. 60 degrees rather than 120 degrees), generating an approximately antiparallel alignment of DNA sequences. This structure is precisely the stacked X-structure proposed on the basis of experimental data. The calculations indicate distortions from standard B-DNA conformation that are required to adopt the stacked X-structure; a widening of the minor groove at the junction, and reorientation of the central phosphate groups of the exchanging strands. An important feature of the stacked X-structure is that it presents two structurally distinct sides. These may be recognised differently by enzymes, providing a rationalisation for the points of cleavage by Holliday resolvases.     

Adsorption of divalent cations to bilayer membranes containing phosphatidylserine

The Stern equation, a combination of the Langmuir adsorption isotherm, the Boltzmann relation, and the Grahame equation from the theory of the diffuse double layer, provides a simple theoretical framework for describing the adsorption of charged molecules to surfaces. The ability of this equation to describe the adsorption of divalent cations to membranes containing brain phosphatidylserine (PS) was tested in the following manner. Charge reversal measurements were first made to determine the intrinsic 1:1 association constants of the divalent cations with the anionic PS molecules: when the net charge of a PS vesicle is zero one-half of the available sites are occupied by divalent cations. The intrinsic association constant, therefore, is equal to the reciprocal of the divalent cation concentration at which the mobility of a PS vesicle reverses sign. The Stern equation with this association constant is capable of accurately describing both the zeta potential data obtained with PS vesicles at other concentrations of the divalent cations and the data obtained with with vesicles formed from mixtures of PS and zwitterionic phospholipids. Independent measurements of the number of ions adsorbed to sonicated PS vesicles were made with a calcium-sensitive electrode. The results agreed with the zeta potential results obtained with multilamellar vesicles. When membranes are formed at 20 degrees C in 0.1 M NaCl, the intrinsic 1:1 association constants of Ni, Co, Mn, Ba, Sr, Ca, and Mg with PS are 40, 28, 25, 20, 14, 12, and 8 M-1, respectively.     

Contributions of long-range electrostatic interactions to 4-chlorobenzoyl-CoA dehalogenase catalysis: a combined theoretical and experimental study

It is well established that electrostatic interactions play a vital role in enzyme catalysis. In this work, we report theory-guided mutation experiments that identified strong electrostatic contributions of a remote residue, namely, Glu232 located on the adjacent subunit, to 4-chlorobenzoyl-CoA dehalogenase catalysis. The Glu232Asp mutant was found to bind the substrate analogue 4-methylbenzoyl-CoA more tightly than does the wild-type dehalogenase. In contrast, the kcat for 4-chlorobenzoyl-CoA conversion to product was reduced 10000-fold in the mutant. UV difference spectra measured for the respective enzyme-ligand complexes revealed an approximately 3-fold shift in the equilibrium of the two active site conformers away from that inducing strong pi-electron polarization in the ligand benzoyl ring. Increased substrate binding, decreased ring polarization, and decreased catalytic efficiency indicated that the repositioning of the point charge in the Glu232Asp mutant might affect the orientation of the Asp145 carboxylate with respect to the substrate aromatic ring. The time course for formation and reaction of the arylated enzyme intermediate during a single turnover was measured for wild-type and Glu232Asp mutant dehalogenases. The accumulation of arylated enzyme in the wild-type dehalogenase was not observed in the mutant. This indicates that the reduced turnover rate in the mutant is the result of a slow arylation of Asp145, owing to decreased efficiency in substrate nucleophilic attack by Asp145. To rationalize the experimental observations, a theoretical model is proposed, which computes the potential of mean force for the nucleophilic aromatic substitution step using a hybrid quantum mechanical/molecular mechanical method. To this end, the removal or reorientation of the side chain charge of residue 232, modeled respectively by the Glu232Gln and Glu232Asp mutants, is shown to increase the rate-limiting energy barrier. The calculated 23.1 kcal/mol free energy barrier for formation of the Meisenheimer intermediate in the Glu232Asp mutant represents an increase of 6 kcal/mol relative to that of the wild-type enzyme, consistent with the 5.6 kcal/mol increase calculated from the difference in experimentally determined rate constants. On the basis of the combination of the experimental and theoretical evidence, we hypothesize that the Glu232(B) residue contributes to catalysis by providing an electrostatic force that acts on the Asp145 nucleophile.     

Adsorption ofRhizobium meliloti to alfalfa roots: Dependence on divalent cations and pH

Adsorption ofRhizobium meliloti L5-30 in low numbers to alfalfa (Medicago sativa L.) roots was dependent on the presence of divalent cations, and required neutral pH. Adsorption was proportional to Ca and/or Mg concentrations up to 1.5 mM. Ca was not substituted by Sr, Ba or Mn. Adsorption was abolished and viability decreased at pH6. When lowering pH, higher Ca concentrations were required to attain similar adsorption levels, indicating a marked interactive effect between Ca and H ions. Pretreatment of the roots with Ca and low pH did not affect subsequent adsorption of the bacteria. However, Ca pretreatment ofR. meliloti sustained further adsorption at low Ca levels and low pH substantially affected their ability to adsorb. Low pH appears to affect the stability of binding causing desorption of the previously bound bacteria. The presence of saturating concentrations of heterologousR. leguminosarum bv.trifolii A118, did not prevent the expression of divalent cations and pH requirements, as well as their interaction. Our results suggest that rhizobial binding to the root surface already shows the Ca and pH dependence of alfalfa nodulation, which was generally associated to some event prior to rhizobial penetration of root hairs.