Hydration structure and dynamics of B- and Z-DNA in the presence of counterions via molecular dynamics simulations

Following our previous attempts at understanding the structural and dynamical properties of water and counterions hydrating nucleic acids, we have performed molecular dynamics simulations for B- and Z-DNA. In these simulations, the nucleic acids were held rigid. In the case of B-DNA, one turn of B-DNA double helix was considered in the presence of 1500 water molecules and 20 counterions (K⁺). The simulations were performed for 4.0 ps after equilibrating the system. For Z-DNA, we considered one turn of the double helix in the presence of 1851 water molecules and 24 counterions (K⁺). The simulations were carried out for 3.5 ps after equilibration. The average temperature of these simulations was ～ 360 K for Z-DNA and ～ 345 K for B-DNA. In these simulations the hydrogen atoms were explicitly taken into account. For both simulations, a fifth-order predictor-corrector was used for solving the translational equations of motion. The rotational motion of the water molecules was represented in terms of quaternion algebra and the rotational equations of motion were solved with a second-order quaternion method using a sixth-order predictor-corrector method. A time step of 0.5 · 10^?15 s was used in these simulations. The structural and the dynamical properties of water solvating the counterions, and the phosphate groups of the DNA, were computed to understand the hydration structure. Diffusion coefficients and velocity correlation functions were calculated for both ions and the water molecules. The velocity correlation functions for the ions exhibit a caged behavior. The dipole correlation functions for the water molecules indicate that the water molecules close to the helix retain the memory of their initial orientations for longer periods of time than those away from the helix. During the time period of our simulation (3–4 ps) the ion probability distributions show a well-defined pattern and suggest limited mobility for the ions, being close to the helix.     

Structural rigidity of paranemic crossover and juxtapose DNA nanostructures

Crossover motifs are integral components for designing DNA-based nanostructures and nanomechanical devices due to their enhanced rigidity compared to the normal B-DNA. Although the structural rigidity of the double helix B-DNA has been investigated extensively using both experimental and theoretical tools, to date there is no quantitative information about structural rigidity and the mechanical strength of parallel crossover DNA motifs. We have used fully atomistic molecular dynamics simulations in explicit solvent to get the force-extension curve of parallel DNA nanostructures to characterize their mechanical rigidity. In the presence of monovalent Na⁺ ions, we find that the stretch modulus (γ₁) of the paranemic crossover and its topoisomer JX DNA structure is significantly higher (∼30%) compared to normal B-DNA of the same sequence and length. However, this is in contrast to the original expectation that these motifs are almost twice as rigid compared to the double-stranded B-DNA. When the DNA motif is surrounded by a solvent with Mg²⁺ counterions, we find an enhanced rigidity compared to Na⁺ environment due to the electrostatic screening effects arising from the divalent nature of Mg²⁺ ions. To our knowledge, this is the first direct determination of the mechanical strength of these crossover motifs, which can be useful for the design of suitable DNA for DNA-based nanostructures and nanomechanical devices with improved structural rigidity.     

Structure of aggregates of water and Li+, Na+, or K+ counterions with nucleic acid in solution

The position and orientation of water molecules hydrating fragments of DNA in the B and Z conformations are analyzed with the help of computer simulations. Monte Carlo studies are carried out at room temperature, high relative humidity (500 water molecules per pitch) and in the presence of counterions such as Li⁺, Na⁺, and K⁺. Differences in hydration patterns and in the counterionic structures were found by compairing B-DNA with Z-DNA double helices and B-DNA helices with different base-pair distributions. The present extension of our similations to Z-DNA and to Li⁺ and K⁺ counterions permits some general conclusions concerning nucleic acids in solution.     

Quantum calculations on water in the KcsA channel cavity with permeant and non-permeant ions

Different ions in the pore of the KcsA channel behave differently, and we relate this to their solvation. We show that the selectivity is dependent, in part, on the solvation in the cavity (sometimes referred to as the vestibule, it is the region containing water molecules between the intracellular gate and the selectivity filter at the extracellular end of the pore). We have shown earlier that potassium is more dependent at the upper end of the cavity region on solvation by the threonines there, while sodium ion has more water molecules as ligands. In addition, sodium ion is placed asymmetrically, while potassium is nearly exactly symmetric with respect to the four-fold symmetry of the channel. We have now extended these calculations to rubidium and cesium ions, and find that rubidium solvation resembles that of potassium (and both are permeant ions), while cesium resembles sodium (and both are non-permeant), in terms of the geometry of up to eight hydrating, and four non-hydrating, water molecules. In each case, a maximum of 12 water molecules are relevant to the calculation. The placement of the water molecules in the two cases is essentially the same as found from the electron density in the X-ray structure of Zhou and MacKinnon. For Na⁺ and K⁺, we show that energy decreases from bulk to the cavity to the lowest position in the selectivity filter (accurate energy could not be calculated for the heavier ions). A separate calculation shows that fixing the Na⁺ ion at the position of the K⁺ minimum, followed by re-optimization produced a significantly modified system, not something that could be produced by thermal fluctuations. Moving the K⁺ into the Na⁺ position in the upper cavity led to a small increase in energy, ≈ 3 k_BT, but was accompanied by large shifts in the positions of hydrating waters, which would create a major kinetic barrier. Therefore, thermal fluctuations could not invalidate the conclusions of the main calculations.     

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

     

Simulations of the solvent structure for macromolecules. I. Solvation of B-DNA double helix at T = 300 K

Monte Carlo simulations are reported for a system of 447 water molecules enclosing a B-DNA double-helix fragment with 12 base pairs and the corresponding sugar and phosphate units. From a detailed analysis on the interaction energies and probability distributions (at a simulated temperature of 300 K), the water molecules can be partitioned into clusters strongly interacting with (1) the phosphates, (2) the sugars, (3) the sugars and the bases, and (4) the base pairs. In addition, transgroove and interphosphate filament of hydrogen-bonded water molecules have been detected. From simulations performed with variable numbers of water molecules, a theoretical isotherm has been obtained, with the characteristic sigmoidal shape, known from absorption–desorption experiments on related systems. The expected main features for the structure of water molecules solvating B-DNA with Na⁺ counterions are briefly discussed at the end of the paper.     

Na+ uptake through the brush border membranes of intestine from fresh water and sea-water adapted trout (Salmo gairdneri,R.)

Summary A comparative study of the mechanisms of Na⁺ absorption through brush border membranes of enterocytes from freshwater (FW) and seawater (SW) adapted trout were carried out using purified vesicle preparations. In contrast to FW trout, SW trout were found to possess a Na⁺–K⁺–Cl^– cotransport process. This finding is regarded as a major adaptation to SW since this cotransport allows an increase of ions and water absorption. Both FW and SW trout were equipped with a Na⁺–H⁺ exchange. In FW, the intestine of the trout had both a Na⁺–Na⁺ exchange and a Na⁺ conductance which may be responsible for enterocyte Na⁺ uptake along the potential gradient.     

Molecular Dynamics Simulation of Limiting Conductance for Na2+, Cl2−, Na°, and Cl° in Supercritical Water

Abstract

We report results of molecular dynamics simulations of the limiting conductance of Na²⁺, Cl^2?, Na°, and Cl° in supercritical water using the SPC/E model for water in conjuction with our previous study (Lee et al., Chem. Phys. Lett. 293, 289 (1998)). The behavior of the limiting conductances of Na²⁺ and Cl^2? in the whole range of water density shows almost the same trend as those of Na⁺ and Cl^?, but the deviation from the assumed linear dependence of limiting conductances of Na²⁺ and Cl^2? on the water density is smaller than that of Na⁺ and Cl^?. The ratio of the limiting conductance of the divalentions to that of the corresponding monovalentions over the whole range of water density is almost constant. In the cases of Na²⁺ and Cl^2?, the dominating factor of the number of hydration water molecules around ions in the higher-density region and the dominating factor of the interaction strength between the ions and the hydration water molecules in the lower-density region are also found as was the cases for Na⁺ and Cl^?. These factors, however, are not so strong as for the corresponding monovalent ions because the change in the energetics, structure, and dynamics are very small mainly due to the strong Coulomb interaction of the divalent ions with the hydration water molecules. The diffusion coefficient of Na° and Cl° monotonically increases with decreasing water density over the whole range of water density. The increase of the diffusion coefficient with decreasing water density is attributed only to the dramatic decrease of the hydration number of water in the first solvation shell around the uncharged species. Among the two important competing factors in the limiting conductance of Na⁺ and Cl^?, the effect of the number of hydration water molecules around the uncharged species is the only existing factor over the whole range of water density since the interaction strength between the uncharged species and the hydration water molecules very small through the LJ interaction. This result has confirmed the dominating factor of the number of hydration water molecules around ions in the higher-density region in the explanation of the limiting conductance of Na⁺ and Cl^? in supercritical water at 673 K.     

Identification of the Third Na Site and the Sequence of Extracellular Binding Events in the Glutamate Transporter

The transport cycle in the glutamate transporter (GlT) is catalyzed by the cotransport of three Na⁺ ions. However, the positions of only two of these ions (Na1 and Na2 sites) along with the substrate have been captured in the crystal structures reported for both the outward-facing and the inward-facing states of Glt_ph. Characterizing the third ion binding site (Na3) is necessary for structure-function studies attempting to investigate the mechanism of transport in GlTs at an atomic level, particularly for the determination of the sequence of the binding events during the transport cycle. In this study, we report a series of molecular dynamics simulations performed on various bound states of Glt_ph (the apo state, as well as in the presence of Na⁺, the substrate, or both), which have been used to identify a putative Na3 site. The calculated trajectories have been used to determine the water accessibility of potential ion-binding residues in the protein, as a prerequisite for their ion binding. Combined with conformational analysis of the key regions in the protein in different bound states and several additional independent simulations in which a Na⁺ ion was randomly introduced to the interior of the transporter, we have been able to characterize a putative Na3 site and propose a plausible binding sequence for the substrate and the three Na⁺ ions to the transporter during the extracellular half of the transport cycle. The proposed Na3 site is formed by a set of highly conserved residues, namely, Asp³¹², Thr⁹², and Asn³¹⁰, along with a water molecule. Simulation of a fully bound state, including the substrate and the three Na⁺ ions, reveals a stable structure—showing closer agreement to the crystal structure when compared to previous models lacking an ion in the putative Na3 site. The proposed sequence of binding events is in agreement with recent experimental models suggesting that two Na⁺ ions bind before the substrate, and one after that. Our results, however, provide additional information about the sites involved in these binding events.     

Hydration valve controlled non-selective conduction of Na and K in the NaK channel

The Na⁺ and K⁺ channels are essential to neural signaling, but our current knowledge at the atomic level is mainly limited to the conducting mechanism of K⁺. Unlike a K⁺ channel having four equivalent K⁺-binding sites in its selectivity filter, a NaK channel has a vestibule in the middle part of its selectivity filter, and can conduct both Na⁺ and K⁺ ions. However, the underlying mechanism for non-selective ion conduction in NaK remains elusive. Here we find four small grottos connecting with the vestibule of the NaK selectivity filter, which form a vestibule-grotto complex perpendicular to the filter pore with a few water molecules within it. It is shown that two or more of the water molecules coming to the vestibule to coordinate the cation are necessary for conducting both Na⁺ and K⁺ ions, while only one water molecule in the vestibule will obstruct ion permeation. Thus, the complex with the aid of interior water movement forms a dynamic hydration valve which is flexible in conveying different cations through the vestibule. Similar exquisite hydration valve mechanisms are expected to be utilized by other non-selective cation channels, and the results should shed new light on the importance of water in neural signaling.     

Circular dichroism of flow-oriented nucleic acids. II. Comparison between observed and calculated spectra.

The ultraviolet circular dichroism (CD) of oriented DNA and RNA molecules is calculated by an extension of Johnson and Tinoco's theory [(1969) Biopolymers 7 , 727–749] for unoriented molecules. The calculations are carried out for molecular models of A-DNA, B-DNA, planar B-DNA, C-DNA, and RNA-11-α. The calculated curves are compared to measured spectra [(1975) J. Mol. Biol. 92 , 449–466] Chung and Holzwarth, for oriented solutions of DNA in buffer, DNA in 6 M LiCl or in ethylene glycol, and double-stranded viral RNA. The calculation, which considers only base–base interactions, predicts that the CD of B-DNA, measured with light propagating parallel to the helix axis, should be large and semiconservative, whereas the CD for light propagating perpendicular to the helix axis should be nonconservative. These predictions agree qualitatively with the experimental observations for DNA in buffer; agreement improves if one assumes the bases to be exactly perpendicular to the helix axis. For the other geometries, agreement is less satisfactory, but qualitative agreement with experiment is obtained and the signs of the specific CD spectra are in accord with observations.     

A Nanosecond Molecular Dynamics Study of Antiparallel d(G)7 Quadruplex Structures: Effect of the Coordinated Cations

Abstract

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G₇) quadruplex structures with different coordinated ions, namely Na⁺ and K⁺ ion, water and Na⁺ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of ～ 1.5 Å from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K⁺ ions are quite immobile. MD studies at 400K indicate that K⁺ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K⁺ ion coordinated structure is more favourable than that for the Na⁺ ion coordinated antiparallel quadruplex structure.     

Electrostatic free energy landscapes for DNA helix bending

Nucleic acids are highly charged polyanionic molecules; thus, the ionic conditions are crucial for nucleic acid structural changes such as bending. We use the tightly bound ion theory, which explicitly accounts for the correlation and ensemble effects for counterions, to calculate the electrostatic free energy landscapes for DNA helix bending. The electrostatic free energy landscapes show that DNA bending energy is strongly dependent on ion concentration, valency, and size. In a Na⁺ solution, DNA bending is electrostatically unfavorable because of the strong charge repulsion on backbone. With the increase of the Na⁺ concentration, the electrostatic bending repulsion is reduced and thus the bending becomes less unfavorable. In contrast, in an Mg²⁺ solution, ion correlation induces a possible attractive force between the different parts of the helical strands, resulting in bending. The electrostatically most favorable and unfavorable bending directions are toward the major and minor grooves, respectively. Decreasing the size of the divalent ions enhances the electrostatic bending attraction, causing an increased bending angle, and shifts the most favorable bending to the direction toward the minor groove. The microscopic analysis on ion-binding distribution reveals that the divalent ion-induced helix bending attraction may come from the correlated distribution of the ions across the grooves in the bending direction.     

First-principles study of water desorption from montmorillonite surface

Knowledge about water desorption is important to give a full picture of water diffusion in montmorillonites (MMT), which is a driving factor in MMT swelling. The desorption paths and energetics of water molecules from the surface of MMT with trapped Li⁺, Na⁺ or K⁺ counterions were studied using periodic density functional theory calculations. Two paths—surface and vacuum desorption—were designed for water desorption starting from a stationary structure in which water bonds with both the counterion and the MMT surface. Surface desorption is energetically more favorable than vacuum desorption due to water–surface hydrogen bonds that help stabilize the intermediate structure of water released from the counterion. The energy barriers of water desorption are in the order of Li⁺?>?Na⁺?>?K⁺, which can be attributed to the short ionic radius of Li⁺, which favors strong binding with the water molecule. The temperature dependence of water adsorption and desorption rates were compared based on the computed activation energies. Our calculations reveal that the water desorption on the MMT surface has a different mechanism from water adsorption, which results from surface effects favoring stabilization of water conformers during the desorption process.     

Hyaluronic acid: molecular conformations and interactions in two sodium salts.

A detailed structure for the tetragonal form (a = b = 0.989 nm, c, fibre axis, = 3.394 nm) of sodium hyaluronate has been obtained by analysing X-ray fibre diffraction data using new molecular modelling techniques. Two polysaccharide chains pass through each unit cell, one at the corner and one at the centre. The chains are anti-parallel to one another. Each chain is a left-handed, 4-fold helix of disaccharide units. There are intramolecular hydrogen bonds stabilising each glycosidic linkage. Octahedrally co-ordinated sodium ions link, by O … Na⁺ … O bridges, neighbouring polysaccharide chains that are further linked by hydrogen bonds. No double-helix model (as originally proposed for this structure) has been found to be free of unacceptable non-bonded contacts or to fit the diffraction intensities as closely.The tetragonal form, which is stable at zero relative humidity, contains no detectable water molecules. At higher relative humidities a related orthorhombic form is observed in which only the a dimension of the lattice is different (a = 1.153 nm, b = 0.989 nm, c = 3.386 nm). In this form the hyaluronate helix is 2-fold with tetrasaccharide units conformationally similar to the 4-fold helix of the tetragonal form. The Na⁺ … O binding and hydrogen bonds lost on expansion of the tetragonal lattice are all replaced in the orthorhombic structure by bridges through water molecules, four of which associated with each tetrasaccharide.     

Molecular Dynamics Simulations Elucidate the Mechanism of Proton Transport in the Glutamate Transporter EAAT3

The uptake of glutamate in nerve synapses is carried out by the excitatory amino acid transporters (EAATs), involving the cotransport of a proton and three Na⁺ ions and the countertransport of a K⁺ ion. In this study, we use an EAAT3 homology model to calculate the pK_a of several titratable residues around the glutamate binding site to locate the proton carrier site involved in the translocation of the substrate. After identifying E374 as the main candidate for carrying the proton, we calculate the protonation state of this residue in different conformations of EAAT3 and with different ligands bound. We find that E374 is protonated in the fully bound state, but removing the Na2 ion and the substrate reduces the pK_a of this residue and favors the release of the proton to solution. Removing the remaining Na⁺ ions again favors the protonation of E374 in both the outward- and inward-facing states, hence the proton is not released in the empty transporter. By calculating the pK_a of E374 with a K⁺ ion bound in three possible sites, we show that binding of the K⁺ ion is necessary for the release of the proton in the inward-facing state. This suggests a mechanism in which a K⁺ ion replaces one of the ligands bound to the transporter, which may explain the faster transport rates of the EAATs compared to its archaeal homologs.     

Molecular Dynamics Simulations Elucidate the Mechanism of Proton Transport in the Glutamate Transporter EAAT3

The uptake of glutamate in nerve synapses is carried out by the excitatory amino acid transporters (EAATs), involving the cotransport of a proton and three Na⁺ ions and the countertransport of a K⁺ ion. In this study, we use an EAAT3 homology model to calculate the pK_a of several titratable residues around the glutamate binding site to locate the proton carrier site involved in the translocation of the substrate. After identifying E374 as the main candidate for carrying the proton, we calculate the protonation state of this residue in different conformations of EAAT3 and with different ligands bound. We find that E374 is protonated in the fully bound state, but removing the Na2 ion and the substrate reduces the pK_a of this residue and favors the release of the proton to solution. Removing the remaining Na⁺ ions again favors the protonation of E374 in both the outward- and inward-facing states, hence the proton is not released in the empty transporter. By calculating the pK_a of E374 with a K⁺ ion bound in three possible sites, we show that binding of the K⁺ ion is necessary for the release of the proton in the inward-facing state. This suggests a mechanism in which a K⁺ ion replaces one of the ligands bound to the transporter, which may explain the faster transport rates of the EAATs compared to its archaeal homologs.     

Solvation free energies of glutamate and its metal complexes: A computer simulation study

Most biological ion channels demonstrate a high degree of selectivity for one type of ion more than others, and in many cases, how they control attaining this is still not clear. So we have studied on some metal ion compounds of glutamate. The Glutamate and its meal ion compounds (Ca²⁺, Na⁺, K⁺ and Li⁺) were first modeled by ab initio calculations and then Monte Carlo simulation was used to calculate solvation free energies and also the complexes free energies for the related structures. The results indicated that Glutamate-Ca²⁺ have more stability in water than other metal ion. Also, it was found out that the more movement in ions; less stability of the structure would result. This trend can be seen both in gas and liquid phase.     

Improved measurements of Na+ fluxes in plants using calixarene-based microelectrodes

Ion-selective microelectrodes are a powerful tool in studying adaptive responses of plant cells and tissues to various abiotic stresses. However, application of this technique in Na⁺ flux measurements was limited due to poor selectivity for Na⁺ ions of commercially available Na⁺ cocktails. Often, these cocktails cannot discriminate between Na⁺ and other interfering ions such as K⁺ and Ca²⁺, leading to inaccurate measurements of Na⁺ concentration and, consequently, inaccurate Na⁺ flux calculations. To overcome this problem, three Na⁺-selective cocktail mixtures were prepared using tetramethoxyethyl ester derivative of p-t-butyl calix[4]arene. These cocktail mixtures were compared with commercially available ETH 227-based Na⁺ cocktail for selectivity for Na⁺ ions over other ions (particularly K⁺ and Ca²⁺). Among the three calixarene-based Na⁺ cocktails tested, cocktail 2 [in % w/w: Na⁺ ionophore (4-tert-butylcalix[4]arene-tetra acetic acid tetraethyl ester) 3.5, the plasticizer (2-nitrophenyl octyl ether) 95.9 and lipophilic anion (potassium tetrakis (4-chlorophenyl) borate) 0.6] showed the best selectivity for Na⁺ ions over K⁺ and Ca²⁺ ions and was highly stable over time (up to 10 h). Na⁺ flux measurements under a wide range of NaCl concentrations (25-150 mM) using Na⁺ cocktail 2 established a clear dose-response relationship between severity of salt stress and magnitude of Na⁺ influx at the distal elongation and mature zones of Arabidopsis thaliana roots. Furthermore, Na⁺ cocktail 2 was compared with commercially available ETH 227-based Na⁺ cocktail by measuring Na⁺ fluxes at the two Arabidopsis root zones in response to 100 mM NaCl treatment. With calixarene-based Na⁺ cocktail 2, a large decreasing Na⁺ influx (0-15 min) followed by small Na⁺ influx (15-45 min) was measured, whereas with ETH-based Na⁺ cocktail Na⁺ influx was short-lived (1-3 min) and was followed by Na⁺ efflux (3-45 min) that might have been due to K⁺ and Ca²⁺ efflux measured together with Na⁺ influx. In conclusion, Na⁺-selective calixarene-based microelectrodes have excellent potential to be used in real-time Na⁺ flux measurements in plants.     

Interaction of phospholipids (Lysophosphatidylethanolamines) with water and sodium cation

We have performed detailed ab initio SCF calculations on the intermolecular interaction energies for one Na⁺ ion and one water molecule with two molecular fragments, one exemplifying a phospholipid (PL) head (PLHD) and the other, a phospholipid tail (PLTL). A 6-12-1 atom-atom pair potential for the interaction of a Na⁺ ion and water with a lysophosphatidyl-ethanolamine (LPEA) was derived from these results by a fitting procedure. This fitted potential was used to obtain isoenergy maps that provide energy profiles of the Na⁺ ion and the water around the phospholipids. The interaction of the Na⁺ ion with PL, as well as the interaction of water with the PL, can be visualized from these maps, which, as expected, show regions of hydrophilicity and hydrophobicity for the water and indicate a very strong binding site for the Na⁺ ion on the phosphate. It appears to be a stationary site that would limit the Na⁺ ion mobility. This binding site is located near the double-bonded oxygen atom of the phosphate group; its binding energy for Na⁺ is 67 kcal/mol. On the other hand the NH⁺ group of PLHD ahows strong electrostatic repulsion of Na⁺ while interacting with water with a binding energy of 13 kcal/mol. This potential energy well region for water is separated from another of similar depth near the phosphate by a barrier and both regions are expected to act as binding sites for water.