Octa-Coordination and the Aqueous Ba2+ Ion - American Chemical

Jun 18, 2015 - Octa-coordination resembles the crystal structure of Ba2+ and. K+ bound in potassium ion channels, but differs from the local hydration...
0 downloads 0 Views 4MB Size
Subscriber access provided by NEW YORK UNIV

Article 2+

Octa-coordination and the Aqueous Ba Ion Mangesh I. Chaudhari, Marielle Soniat, and Susan B. Rempe J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.5b03050 • Publication Date (Web): 18 Jun 2015 Downloaded from http://pubs.acs.org on June 21, 2015

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry B is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Octa-coordination and the Aqueous Ba2+ Ion Mangesh I. Chaudhari,∗,† Marielle Soniat,∗,‡ and Susan B. Rempe∗,† Center for Biological and Material Sciences, Sandia National Laboratories, Albuquerque, NM, 87185, and Department of Chemistry, University of New Orleans, New Orleans, LA, 70148 E-mail: [email protected]; [email protected]; [email protected]

Abstract The hydration structure of Ba2+ ion is important for understanding blocking mechanisms in potassium ion channels. Here, we combine statistical mechanical theory, ab initio molecular dynamics simulations, and electronic structure methods to calculate the hydration free energy and local hydration structure of Ba2+ (aq). The predicted hydration free energy (- 304±1 kcal/mol) agrees with the experimental value (-303 kcal/mol) when a maximally occupied, unimodal inner solvation shell is treated. In the local environment defined by the first shell of hydrating waters, Ba2+ is directly and stably coordinated by eight (8) waters. Octacoordination resembles the crystal structure of Ba2+ and K+ bound in potassium ion channels, but differs from the local hydration structure of K+ (aq) determined earlier.

Introduction

ter and ionic solutions, Bernal and Fowler proposed an unusual octa-coordinated Ba2+ hydration structure as an exception to the anticipated hexa-coordination of smaller divalent ions. 24 Structural data to test that proposition directly is sparse. One reason is that Ba2+ extensively absorbs X-rays, leading to unfavorable conditions for structural studies. 25 Consequently, only a few experimental groups have targeted Ba2+ ion, producing discrepant results. X-ray diffraction studies reported a hydration number greater than eight (9.5), 26 while extended X-ray absorption fine structure (EXAFS) spectroscopy experiments reported hydration numbers of eight (8.1 and 7.8). 27,28 Similarly, theoretical analysis of Ba2+ hydration using ab initio methods has been limited, partly due to the large number of electrons involved. 29 No purely ab initio studies of hydration structure in aqueous solution have been reported. One study combined an ab initio quantum mechanical approach with classical molecular mechanics simulations (QM/MM) and re-

Barium (Ba2+ ) ion is about the same size as potassium (K+ ) based on crystal ionic radii that fall within 0.03 ˚ A. 1 Both ions partition between water and octa-coordinated binding sites in potassium ion channels. 2–6 Thus, Ba2+ can act as an analogue of K+ . Octa-ligation by oxygens along the protein backbone of K channels is widely believed to facilitate K+ permeation by mimicking that ion’s local hydration structure. 5,7,8 While K+ permeates rapidly, Ba2+ instead sticks and blocks permeation of other ions. That inhibitory behavior has been used since the 1970’s to probe the mechanism of K channel function. 3,8–20 Yet, recent works still debate the blocking mechanism and structure of the blocking sites in various K channels. 21–23 To help clarify the debate, we analyze here the physical chemical properties of Ba2+ in aqueous solution, the reference environment for ion channel block. In a 1933 landmark theoretical paper on wa-

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ported a hydration number of nine (9.3). 29 That structural result was not substantiated by a prediction of hydration free energy. Classical simulations have produced values ranging from eight 30,31 to nine. 28 Those simulations are typically not designed to provide sole determination of such properties, but they may shed light on issues determining the hydration number of ions. A statistical mechanical theory developed earlier permits computation of solvation free energy based on local structural results determined for systems treated with ab initio models. 32–36 The coupling of structure with thermodynamic predictions provides an advantage for validating the results compared to more standard studies of structure alone. Also, structural data obtained by molecular simulation contains information about the spatial distributions of each neighboring solvent molecule. Although a neighborship analysis is typically unresolvable for experimental data and seldom applied to simulation data, the results can help clarify how many solvent molecules define the local hydration structure around an ion. The free energy theory has been coupled with structural studies to obtain new insights about K+ 37–39 and its monovalent analogue, rubidium (Rb+ ). 40 Here, we apply this approach for the first time to Ba2+ (aq) to resolve fundamental questions about that ion’s interactions in bulk liquid water. We base the structural analysis on the first purely ab initio studies of Ba2+ (aq). At the same time, we also investigate how choices in the theoretical analysis affect the free energy predictions. Our studies suggest that the best hydration free energy predictions result from analysis of a maximally occupied, unimodal inner solvation shell, as determined by a neighborship analysis. Those restrictions on inner-shell occupation and boundary arise because of the dielectric continuum approximation used here to treat water outside the inner shell. Contributions from water within and outside the inner shell are both important for predicting the total hydration free energy. Barium ion directly coordinates with n=8 waters in aqueous solution, in contrast to lower hydration numbers

Page 2 of 24

reported earlier in experimental and theoretical studies of K+ 37–39,41–44 and its analogue, Rb+ . 40,45–47 An interesting observation is that octa-coordinated ion binding sites in K channels appear to mimic the local hydration structure of Ba2+ , a blocking ion, not the permeant ions, K+ and Rb+ .

Theory Quasi-chemical theory (QCT) 32–36,39 divides the excess chemical potential of Ba2+ hydra(ex) tion, µBa2+ , into three contributions, (ex)

µBa2+ = −kT lnKn(0) ρnH2 O + kT lnpBa2+ (n) (ex)

+ (µBa(H

(ex)

2+ 2 O)n

− nµH2 O ). (1)

The first term represents an equilibrium ratio (0) Kn for Ba2+ -water association reactions (Eq. 2) treated as in an ideal gas phase, hence the superscript (0): Ba2+ + nH2 O ⇋ Ba(H2 O)2+ n .

(2)

The association reactions occur within a chosen observation volume. While any shape may be selected, 35,36 we choose a sphere of radius r = λ centered around the ion for convenience. Other works have also chosen a spherical boundary. 37–39,44,48–55 Solvent molecules within λ, measured here by the location of water oxygen atoms, are local to the ion and constitute inner-shell solvent. Accordingly, the reactions form inner-shell solvation complexes, Ba(H2 O)2+ n . The density of water in solution, ρH2 O , accounts for the availability of (0) water ligands. Both Kn and ρH2 O are made dimensionless by scaling with the density of water equivalent to an ideal gas at ambient pressure, p=1 atm. In the second term, pBa2+ (n) is the probability of observing n waters around Ba2+ within the inner solvation shell of radius λ. This term accounts for fluctuations in the inner-shell solvent population in aqueous solution. The population fluctuation term will contribute zero if the inner shell is strictly occupied by a single

ACS Paragon Plus Environment

2

Page 3 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Computational Methods

coordination number, n, so that p(n)=1. The probability is readily evaluated from ab initio molecular dynamics (AIMD) simulations. The third term represents solvation of a Ba2+ water inner-shell complex by the outer solvation environment, and removal of the water ligands from the same environment. Note that the environment accounts for the compositions of bulk solutions studied experimentally. That combi(ex) (ex) nation, µX(H O)2+ − nµH2 O , balances the free en2 n ergy for the association reaction when the solute, X, is only weakly solvated. Gas molecules provide example cases. 36,56–59 A balance is not anticipated here because Ba2+ hydration free (ex) energy, µBa2+ , is large. The Boltzmann factor, k, and absolute temperature, T , set the energy scale. (ex) The overall hydration free energy, µBa2+ , is independent of the choices made for λ and n in QCT (Eq. 1). In fact, Eq. 1 applies to any single n that occurs within the chosen innershell boundary, λ. Nevertheless, specific n or λ may be more effective than others in practical applications that evaluate the free energy contributions approximately. One goal pursued here is to determine how different choices for the radius (λ) and occupancy (n) of the inner solvation shell affect QCT predictions when the terms (Eq. 1) are evaluated by electronic structure, AIMD, and dielectric continuum models. To achieve that goal, we analyze AIMD simulation data to inform the selection of the inner shell boundary and occupancy. In place of the single λ value considered previously, our study compares hydration free energies computed with seven λ values (2.9-3.5 in increments of 0.1 ˚ A). Those values define inner-shell regions that include specific structural features found in the simulated Ba2+ -oxygen radial distribution function, g(r). We also consider n that span the full range of possible occupancy distributed exclusively within the first hydration shell of Ba2+ (aq). Those occupancies are identified by a neighborship analysis of g(r).

The main challenge in predicting local hydration structure and hydration free energy for ions is to represent the broad range of molecular interactions involved in ion complexation reactions. Previous studies suggest that treatment of multi-body interactions is important for those predictions. 60–62 Accordingly, we chose to model Ba2+ and near water neighbors using density functional theory (DFT) since that approach naturally accounts for interactions between pairs and larger groups of atoms. Ab initio molecular dynamics (AIMD) simulations on a single Ba2+ solvated by 64 waters were performed using the VASP AIMD simulation package 63 to investigate local hydration structure. Simulations were carried out for a total time of 41 ps. A system with 64 waters extends well beyond the first hydration shell of Ba2+ , which is the region of interest in this study. The ion/water ratio was chosen to match a 0.8 M ion concentration studied experimentally. 27 The simulated AIMD system was defined by a cubic box of 12.417 ˚ A with periodic boundary conditions to mimic bulk liquid conditions. In that box volume, the water density matches the experimental density of liquid water at standard conditions of room temperature and pressure used in experimental structural studies. 27,28 The boundaries contain a background charge to neutralize the overall charge of the system. We utilized the PW91 generalized gradient approximation (GGA) to the electron density 64 to facilitate comparisons with earlier work on K+ 43,44 and Rb+ 40 hydration. As a pure GGA functional, PW91 enables a significant computational savings over more specialized functionals such as the hybrid GGA, B3LYP, and metaGGA, TPSS. We described the core-valence interactions using the projector augmented-wave (PAW) method, 65 expanded the valence electronic orbitals in plane waves with a high kinetic energy cut-off of 36.75 Ry (500 eV), and used 10−6 eV as the convergence criteria for the electronic structure self-consistent iterations. We also used a time step of 0.5 fs to resolve water vibrational motions. Prior to the pro-

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 4 of 24

Page 5 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

fied to match the chosen λ value. Additionally, default parameters from the PCM model were used for hydrogen and oxygen radii in the clusters and for individual waters. The radii define cavities for the PCM solvation model based on a set of overlapping spheres. The dielectric constant of the outer-shell environment was set to mimic water at room temperature (78). The ion-water complexes were re-optimized in the presence of the environment. Hydration free energies were calculated at standard conditions of T =298 K and p=1 atm by adding the individual contributions of Eq. 1. The density in the first term was subsequently adjusted to account for the actual concentration of water ligands in liquid water, ρH2 O =1 g cm−3 , to match experimental conditions. If this density is tracked as an adjustment to the ideal gas pressure, then it corresponds to a pressure factor of 1354 atm. 32,89 Predicted values are compared with the absolute hydration free energy compiled by Marcus. 85 Due to a sign error, the compiled value was adjusted by 2 × −1.9 kcal/mol to convert to the standard state pertaining to hydration (1M Ba2+ in gas phase and aqueous solution). 32

determined gas phase free energies for sequential addition of waters to Ba2+ under standard conditions of room temperature and ambient pressure provided data 74 for selecting the exchange-correlation functional and basis sets (FIG. 1). A comparison of density functionals led to selection of the TPSS exchange-correlation density functional with the aug-cc-pvtz (oxygen) and cc-pvtz (hydrogen) basis sets 83 and the MWB46 84 (Ba) effective core potential and corresponding basis set. The slightly smaller basis on hydrogen atoms results in negligible changes to structure and electronic energy, but represents as much as a five-fold savings in computational time compared with the aug-cc-pvtz basis set. As anticipated, the exchange-correlation functionals produced nearly identical lowestenergy structures. Note also that differences in binding free energies among functionals are small relative to total ion hydration free energy (≈ -300 kcal/mol 85 ), particularly between TPSS, PW91, and PBE (FIG. 1). That agreement means that the final results from combining terms in Eq. 1 will not depend sensitively on the level of electronic structure theory. Cluster conformations were exhaustively sampled from the ab initio simulation trajectories and then submitted to the G09 software to obtain optimized geometries and electronic energies. Tight convergence criteria on the optimization (10−5 a.u.) and energy (10−8 a.u.), along with an ultra-fine integration grid, facilitated the optimization procedure. Vibrational frequency analysis based on the normal modes 86,87 were performed to obtain thermal corrections to the electronic energy under the same standard conditions studied experimentally. All vibrational frequencies were positive, confirming that optimized cluster configurations represent minimum-energy geometries. To evaluate the outer-shell solvation term of Eq. 1, the outer environment was treated implicitly as a reaction field with a polarizable continuum model (PCM). 88 The PCM model includes contributions to the outer-shell solvation free energy from electrostatic, packing, and dispersion terms. 88 The radius of the barium atom within each Ba(H2 O)2+ n cluster was modi-

Results and Discussion The hydration structure determined by the Ba2+ -oxygen radial distribution function, g(r), shows a distinct division between water neighbors occupying the first and second hydration shells (FIG. 2). Analysis of the near neighbor distributions reveals that n=1-6 waters fill in the principal maximum at rmax ≈ 2.8 ˚ A. That position agrees with experimental results of rmax =2.81 ˚ A 27 and 2.78 ˚ A. 28 Those waters, and the n=7th and n=8th near neighbors, directly coordinate the ion and contribute to the first peak in a unimodal way. An inflection point on the running coordination number hn(r)i at the first minimum in g(r), rmin ≈ 3.5 ˚ A, defines the first hydration shell and confirms h¯ ni = 8 as the average and most probable hydration number, in agreement with the average hydration structure from EXAFS data of hni=8.1 27 and 7.8. 28 The radial distri-

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 6 of 24

Page 7 of 24

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 8 of 24

Page 9 of 24

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

permeant K+ , but not the smaller impermeant sodium ion (Na+ ). 39,44,96 An interesting observation is that ion binding sites in K channels appear to mimic the local hydration structure of a blocking ion (Ba2+ ), not the permeant ion (K+ ). This result establishes a foundation for future molecular studies of the blocking mechanism and blocking sites of K channels.

Page 10 of 24

and barium binding in NaK2K. J. Gen. Phys. 2014, 144, 181–192. (6) Koepfer, D.; Song, C.; Gruene, T.; Sheldrick, G.; Zachariae, U.; de Groot, B. L. Ion permeation in K+ channels occurs by direct Coulomb knock-on. Science 2014, 346, 352 – 355. (7) Zhou, Y.; Morais-Cabral, J.; Kaufman, A.; MacKinnon, R. Chemistry of ion coordination and hydration revealed by a K+ channel-Fab complex at 2.0 ˚ A resolution. Nature 2001, 414, 43–48.

Acknowledgements We thank Lawrence R. Pratt, Sameer Varma, and Dubravko Sabo for helpful discussions. Sandia National Laboratories is a multiprogram laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy’s National Nuclear Security Administration under Contract DE-AC0494AL8500. This work was supported by Sandia’s LDRD program (M. I. C. and S. B. R.) and the State of Louisiana Board of Regents (M. S.).

(8) Piasta, K. N.; Theobald, D. L.; Miller, C. Potassium-selective block of barium permeation through single KcsA channels. J. Gen. Phys. 2011, 138, 421–436. (9) Hagiwara, S.; Miyazaki, S.; Moody, W.; Patlak, J. Blocking effects of barium and hydrogen ions on the potassium current during anomalous rectification in the starfish egg. J. Physiol. 1978, 279, 167– 185.

References

(10) Eaton, D. C.; Brodwick, M. S. Effects of barium on the potassium conductance of squid axon. J. Gen. Phys. 1980, 75, 727– 750.

(1) Shannon, R. D. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides. Acta Cryst. 1976, A32, 751–767.

(11) Armstrong, C. M.; Swenson, R. P.; Taylor, S. R. Block of squid axon K-channels by internally and externally applied barium ions. J. Gen. Phys. 1982, 80, 663– 682.

(2) Jiang, Y.; MacKinnon, R. The barium site in a potassium channel by x-ray crystallography. J. Gen. Phys. 2000, 115, 269. (3) Lockless, S.; Zhou, M.; MacKinnon, R. Structural and thermodynamic properties of selective ion binding in a K+ channel. PLoS Bio. 2007, 5, e121.

(12) Miller, C.; Latorre, R.; Reisin, I. Coupling of voltage-dependent gating and Ba++ block in the high- conductance, Ca++ -activated K+ channel. J. Gen. Phys. 1987, 90, 427–449.

(4) Guo, R.; Zeng, W.; Cui, H.; Chen, L.; Ye, S. Ionic interactions of Ba2+ blockades in the MthK K+ channel. J. Gen. Phys. 2014, 144, 193–200.

(13) Neyton, J.; Miller, C. Discrete Ba2+ block as a probe of ion occupancy and pore structure in the high-conductance Ca2+ -activated K+ channel. J. Gen. Phys. 1988, 92, 569–586.

(5) Lam, Y. L.; Zeng, W.; Sauer, D.; Jiang, Y. The conserved potassium channel filter can have distinct ion binding profiles: Structural analysis of rubidium, cesium,

ACS Paragon Plus Environment

10

Page 11 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(14) Harris, R. E.; Larsson, H. P.; Isacoff, E. Y. A permeant ion binding site located between two gates of the Shaker K+ channel. Biophys. J. 1998, 74, 1808–1820.

(22) Rossi, M.; Tkatchenko, A.; Rempe, S. B.; Varma, S. Role of methyl-induced polarization in ion binding. Proc. Nat. Ac. Sci. 2013, 110, 12978–12983.

(15) Vergara, C.; Alvarez, O.; Latorre, R. Localization of the K+ lock-in and the Ba2+ binding sites in a voltage-gated calciummodulated channel. J. Gen. Phys. 1999, 114, 365–376.

(23) Rowley, C. N.; Roux, B. A computational study of barium blockades in the KcsA potassium channel based on multi-ion potential of mean force calculations and free energy perturbation. J. Gen. Phys. 2013, 142, 451–463.

(16) Alagem, N.; Dvir, M.; Reuveny, E. Mechanism of Ba2+ block of a mouse inwardly rectifying K+ channel: differential contribution by two discrete residues. J. Phys. 2001, 534, 381–393.

(24) Bernal, J. D.; Fowler, R. H. A theory of water and ionic solution, with particular reference to hydrogen and hydroxyl ions. J. Chem. Phys. 1933, 1, 515–548. (25) Ohtaki, H.; Radnai, T. Structure and dynamics of hydrated ions. Chem Rev 1993, 93, 1157–1204.

(17) Proks, P.; Antcliff, J. F.; Ashcroft, F. M. The ligand-sensitive gate of a potassium channel lies close to the selectivity filter. EMBO Reports 2003, 4, 70–75.

(26) Albright, J. N. X-Ray Diffraction Studies of Aqueous Alkaline-Earth Chloride Solutions. J. Chem. Phys. 1972, 56, 3783.

(18) Krishnan, M. N.; Bingham, J.-P.; Lee, S. H.; Trombley, P.; Moczydlowski, E. Functional role and affinity of inorganic cations in stabilizing the tetrameric structure of the KcsA K+ channel. J. Gen. Phys. 2005, 126, 271–283.

(27) Persson, I.; Sandstrom, M.; Yokoyama, H.; Chaudhry, M. Structure of the solvated strontium and barium ions in aqueous, dimethyl-sulfoxide and pyridine solution, and crystal-structure of strontium and barium hydroxide octahydrate. Zeitschr. Naturforsch. Sec. A - J. Phys. Sci. 1995, 50, 21–37.

(19) Krishnan, M. N.; Trombley, P.; Moczydlowski, E. G. Thermal stability of the K+ channel tetramer: cation interactions and the conserved threonine residue at the innermost site (S4) of the KcsA selectivity filter. Biochem. 2008, 47, 5354–5367.

(28) D’Angelo, P.; Pavel, N.; Roccatano, D.; Nolting, H. F. Multielectron excitations at the L edges of barium in aqueous solution. Phys. Rev. B 1996, 54, 12129–12138.

(20) Chatelain, F. C.; Gazzarrini, S.; Fujiwara, Y.; Arrigoni, C.; Domigan, C.; Ferrara, G.; Pantoja, C.; Thiel, G.; Moroni, A.; Minor, D. L. Selection of inhibitor-resistant viral potassium channels identifies a selectivity filter site that affects barium and amantadine block. PloS One 2009, 4, e7496.

(29) Hofer, T. S.; Rode, B. M.; Randolf, B. R. Structure and dynamics of solvated Ba(II) in dilute aqueous solution – an ab initio QM/MM MD approach. Chem. Phys. 2005, 312, 81–88. (30) Stack, A. G.; Rustad, J. R. Structure and dynamics of water on aqueous barium ion and the 001 barite surface. J. Phys. Chem. C 2007, 111, 16387–16391.

(21) Kim, I.; Allen, T. W. On the selective ion binding hypothesis for potassium channels. Proc. Nat. Ac. Sci. 2011, 108, 17963–17968.

(31) Yu, H.; Whitfield, T.; Harder, E.; Lamoureux, G.; Vorobyov, I.; Anisimov, V.;

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 24

changes in competition for ligand binding. J. Am. Chem. Soc. 2008, 130, 15405– 15419.

Jr., A. M.; Roux, B. Simulating Monovalent and Divalent Ions in Aqueous Solution Using a Drude Polarizable Force Field. J. Chem. Theo. and Comp. 2010, 6, 774–786.

(40) Sabo, D.; Jiao, D.; Varma, S.; Pratt, L. R.; Rempe, S. B. Case study of Rb+ (aq), quasi-chemical theory of ion hydration, and the no split occupancies rule. Ann. Rep. Prog. Chem, Sect. C (Phys. Chem.) 2013, 109, 266–278.

(32) Pratt, L. R.; Rempe, S. B. In Simulation and Theory of Electrostatic Interactions in Solution; Hummer, G., Pratt, L. R., Eds.; AIP: New York, 1999; pp 177–201.

(41) Neilson, G. W.; Skipper, N. T. K+ coordination in aqueous solution. Chem. Phys. Lett. 1985, 114, 35.

(33) Beck, T. L.; Paulaitis, M. E.; Pratt, L. R. The Potential Distribution Theorem and Models of Molecular Solutions; Cambridge University Press, 2006.

(42) Glezakou, V.; Chen, Y.; Fulton, J.; Schenter, G.; Dang, L. X. Electronic structure, statistical mechanical simulations, and EXAFS spectroscopy of aqueous potassium. Theor. Chem. Acc. 2006, 115, 86– 99.

(34) Asthagiri, D.; Dixit, P. D.; Merchant, S.; Paulaitis, M. E.; Pratt, L. R.; Rempe, S. B.; Varma, S. Ion selectivity from local configurations of ligands in solutions and ion channels. Chem. Phys. Lett. 2010, 485, 1–7.

(43) Varma, S.; Rempe, S. B. Coordination numbers of alkali metal ions in aqueous solutions. Biophys. Chem. 2006, 124, 192– 199.

(35) Rogers, D. M.; Rempe, S. B. Probing the thermodynamics of competitive ion binding using minimum energy structures. J. Phys. Chem. B 2011, 115, 9116–9129.

(44) Varma, S.; Rempe, S. B. Tuning ion coordination preferences to enable selective permeation. Biophys. J. 2007, 93, 1093– 99.

(36) Rogers, D. M.; Jiao, D.; Pratt, L.; Rempe, S. B. Structural models and molecular thermodynamics of hydration of ions and small molecules. Ann. Rep. Comp. Chem. 2013, 8, 71–128.

(45) Ramos, S.; Barnes, A. C.; Neilson, G. W.; Capitan, M. Anomalous X-ray diffraction studies of hydration effects in concentrated aqueous electrolyte solutions. Chem. Phys. 2000, 258, 171.

(37) Rempe, S. B.; Asthagiri, D.; Pratt, L. R. Inner shell definition and absolute hydration free energy of K+ (aq) on the basis of quasi-chemical theory and ab initio molecular dynamics. Phys. Chem. Chem. Phys. 2004, 6, 1966.

(46) Fulton, J. L.; Pfund, D. M.; Wallen, S. L.; Newville, M.; Stern, E. A.; Ma, Y. Rubidium ion hydration in ambient and supercritical water. J. Chem. Phys. 1996, 105, 2161.

(38) Asthagiri, D.; Pratt, L. R.; Paulaitis, M. E.; Rempe, S. B. Hydration structure and free energy of biomolecularly specific aqueous dications, including Zn2+ and first transition row metals. J. Am. Chem. Soc. 2004, 126, 1285–1289.

(47) Filipponi, A.; Panfilis, S. D.; Oliva, C.; Ricci, M. A.; D’Angelo, P.; Bowron, D. T. Ion hydration under pressure. Phys. Rev. Lett. 2003, 91, 165505.

(39) Varma, S.; Rempe, S. B. Structural transitions in ion coordination driven by

(48) Rempe, S. B.; Pratt, L. R.; Hummer, G.; Kress, J. D.; Martin, R. L.; Redondo, A.

ACS Paragon Plus Environment

12

Page 13 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The hydration number of Li+ in liquid water. J. Am. Chem. Soc. 2000, 122, 966– 967.

(58) Sabo, D.; Varma, S.; Martin, M. G.; Rempe, S. B. Studies of the thermodynamic properties of hydrogen gas in bulk water. J. Phys. Chem. B 2008, 112, 867– 876.

(49) Rempe, S. B.; Pratt, L. R. The hydration number of Na+ in liquid water. Fl. Ph. Equ. 2001, 121–132.

(59) Chaudhari, M. I.; Sabo, D.; Pratt, L. R.; Rempe, S. B. Hydation of Kr(aq) in dilute and concentrated solutions. J. Phys. Chem. B 2014, doi: 10.1021/jp508866h.

(50) Ashbaugh, H. S.; Asthagiri, D.; Pratt, L. R.; Rempe, S. B. Hydration of krypton and consideration of clathrate models of hydrophobic effects from the perspective of quasi-chemical theory. Biophys. Chem. 2003, 105, 323–338.

(60) Whitfield, T. W.; Varma, S.; Harder, E.; Lamoureux, G.; Rempe, S. B.; Roux, B. Theoretical study of aqueous solvation of K+ comparing ab initio, polarizable, and fixed-charge models. J. Chem. Theo. Comp. 2007, 3, 2068–2082.

(51) Asthagiri, D.; Pratt, L. R. Quasi-chemical study of Be2+ (aq) speciation. Chem. Phys. Lett. 2003, 371, 613–619.

(61) Varma, S.; Rempe, S. B. Multibody effects in ion binding and selectivity. Biophys. J. 2010, 99, 3394–3401.

(52) Jiao, D.; Leung, K.; Rempe, S. B.; Nenoff, T. M. First principles calculations of atomic bickel redox potentials and dimerization free energies: A study of metal nanoparticle growth. J. Chem. Theo. Comp. 2011, 7, 485–495.

(62) Bostick, D. L.; Brooks, C. L. Selective complexation of K+ and Na+ in simple polarizable ion-ligating systems. J. Am. Chem. Soc. 2010, 132, 13185–13187. (63) Kresse, G.; Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 1993, 47, 558–561.

(53) Clawson, J. S.; Cygan, R. T.; Alam, T. M.; Leung, K.; Rempe, S. B. Ab initio study of hydrogen storage in water clathrates. J. Comput. and Theor. Nanosci. 2010, 7, 2602–2606.

(64) Perdew, J. P.; Chevary, J. A.; Vosko, S. H.; Jackson, K. A.; Pederson, M. R.; Singh, D. J.; Fiolhais, C. Atoms, molecules, solids, and surfaces - applications of the generalized gradient approximation for exchange and correlation. Phys. Rev. B 1992, 46, 6671–6687.

(54) Varma, S.; Sabo, D.; Rempe, S. B. K+ /Na+ selectivity in K channels and valinomycin: Over-coordination versus cavity-size constraints. J. Mol. Bio. 2008, 376, 13–22. (55) Jiao, D.; Rempe, S. B. Combined DFT and continuum calculations of pKa in carbonic anhydrase. Biochem 2012, 51, 5979–5989.

(65) Bl¨ochl, P. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953– 17979. (66) Alam, T. M.; Hart, D.; Rempe, S. L. B. Computing the 7 Li NMR chemical shielding of hydrated Li+ using cluster calculations and time-averaged configurations from ab initio molecular dynamics simulations. Phys. Chem. Chem. Phys. 2011, 13, 13629–13637.

(56) Jiao, D.; Rempe, S. B. CO2 solvation free energy using quasi-chemical theory. J. Chem. Phys. 2011, 134, 224506. (57) Sabo, D.; Rempe, S. B.; Greathouse, J. A.; Martin, M. G. Molecular studies of the structural properties of hydrogen gas in bulk water. Mol. Sim. 2006, 32, 269–278.

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(67) Ding, Y.; Hassanali, A. A.; Parrinello, M. Anomalous water diffusion in salt solutions. Proc. Nat. Acad. Sci. U.S.A. 2014, 111, 3310–3315.

Page 14 of 24

(76) Vosko, S. H.; Wilk, L.; Nusair, M. Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis. Can. J. Phys. 1980, 58, 1200–1211.

(68) Bankura, A.; Carnevale, V.; Klein, M. L. Hydration structure of Na+ and K+ from ab initio molecular dynamics based on modern density functional theory. Mol. Phys. 2014, 112, 1448–1456.

(77) Zhao, Y.; Truhlar, D. G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215 – 241.

(69) Gaiduk, A. P.; Zhang, C.; Gygi, F.; Galli, G. Structural and electronic properties of aqueous NaCl solutions from ab initio molecular dynamics simulations with hybrid density functionals. Chem. Phys. Lett. 2014, 604, 89–96.

(78) Ernzerhof, M.; Perdew, J. P. Generalized gradient approximation to the angle- and system-averaged exchange hole. J. Chem. Phys. 1998, 109, 3313 – 3320.

(70) VandeVondele, J.; Mohamed, F.; Krack, M.; Hutter, J.; Sprik, M.; Parrinello, M. The influence of temperature and density functional models in ab initio molecular dynamics simulation of liquid water. J. Chem. Phys. 2005, 122, 014515.

(79) Yanai, T.; Tew, D. P.; Handy, N. C. A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP). Chem. Phys. Lett. 2004, 393, 51–57.

(71) Rempe, S. B.; Mattsson, T. R.; Leung, K. On ”the complete basis set limit” and plane-wave methods in first-principles simulations of water. Phys. Chem. Chem. Phys. 2008, 10, 4685–4687.

(80) Tao, J. M.; Perdew, J. P.; Staroverov, V. N.; Scuseria, G. E. Climbing the density functional ladder: Non-empirical meta-generalized gradient approximation designed for molecules and solids. Phys. Rev. Lett. 2003, 91, 146401.

(72) Berendsen, H. J. C.; Grigera, J. R.; Straatsma, T. P. The missing term in effective pair potentials. J. Phys. Chem. 1987, 91, 6269–6271.

(81) Frisch, M. J. et al. Gaussian 09 Revision A.1. Gaussian Inc. Wallingford CT 2009.

(73) Mamatkulov, S.; Fyta, M.; Netz, R. Force fields for divalent cations based on singleion and ion-pair properties. J. Chem. Phys. 2013, 138, 024505.

(82) Soniat, M.; Rogers, D. M.; Rempe, S. B. Dispersion- and exchange-corrected density functional theory for sodium ion hydration. J. Chem. Theory Comput. 2015,

(74) Peschke, M.; Blades, A. T.; Kebarle, P. Hydration energies and entropies for Mg2+ , Ca2+ , Sr2+ , and Ba2+ from gasphase ion-water molecule equilibria determinations. J. Phys. Chem. A 1998, 102, 9978–9985.

(83) Kendall, R.; Jr., T. D.; Harrison, R. Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions. J. Chem. Phys. 1992, 96, 6796 – 806. (84) Bergner, A.; Dolg, M.; Kuechle, W.; Stoll, H.; Preuss, H. Ab initio energyadjusted pseudopotetials for elements of

(75) Becke, A. D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652.

ACS Paragon Plus Environment

14

Page 15 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(94) Kathmann, S. Mundy, C. J.; derstanding the water. J. Phys. 4369.

groups 13-17 . Mol. Phys. 1993, 80, 1431 – 1441. (85) Marcus, Y. A simple empirical model describing the thermodynamics of hydration of ions of widely varying charges, sizes, and shapes. Biophys. Chem. 1994, 51, 111–127.

M.; Kuo, I. W.; Schenter, G. K. Unsurface potential of Chem. B 2011, 115,

(95) Beck, T. The influence of water interfacial potentials on ion hydration in bulk water and near interfaces. Chem. Phys. Lett. 2013, 561-562, 1–13.

(86) Rempe, S. B.; J´onsson, H. A computational exercise illustrating molecular vibrations and normal modes. Chem. Ed. 1998, 3, 1–17.

(96) Varma, S.; Rogers, D. M.; Pratt, L. R.; Rempe, S. B. Design principles for K+ selectivity in membrane transport. J. Gen. Phys. 2011, 137, 479–488.

(87) Rempe, S. B.; Watts, R. The exact quantum mechanical kinetic energy operator in internal coordinates for vibration of a hexatomic molecule. J. Chem. Phys. 1998, 108, 10084. (88) Tomasi, J.; Mennucci, B.; Cammi, R. Quantum mechanical continuum solvation models. Chem. Rev. 2005, 105, 2999– 3093. (89) Grabowski, P.; Riccardi, D.; Gomez, M. A.; Asthagiri, D.; Pratt, L. R. Quasi-Chemical Theory and the Standard Free Energy of H+ (aq). J. Phys. Chem. A 2002, 106, 9145– 48. (90) Shah, J. K.; Asthagiri, D.; Pratt, L. R.; Paulaitis, M. E. Balancing local order and long-ranged interactions in the molecular theory of liquid water. J. Chem. Phys. 2007, 127, 144508. (91) Asthagiri, D.; Pratt, L. R.; Ashbaugh, H. Absolute hydration free energies of ions, ion–water clusters, and quasichemical theory. J. Chem. Phys. 2003, 119, 2702. (92) Leung, K.; Rempe, S. B.; von Lilienfeld, O. A. Ab initio molecular dynamics calculations of ion hydration free energies. J. Chem. Phys. 2009, 130, 204507– 204517. (93) Leung, K. Surface potential at the airwater interface computed using density functional theory. J. Chem. Phys. Lett. 2010, 1, 496–499.

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 16 of 24

Page 17The of 24 Journal of Physical Chemistry

−kT lnK(n−1)→n (kcal/mol)

−5

(0)

1 2 3 4 −10 B3LYP 5 M06 6 7 PW91 8 PBE0 9 −15 10 CAM − B3LYP 11 TPSS 12 Experiment ACS Paragon Plus Environment 13 −20 14 5 6 7 8 15

n

gBa2+ O (r)

1 2 3 4 5 6 7 8 9 10 11 12 13

10

The Journal of Physical Chemistry Page 18 of 24

8

8

6

6 1-6

4

4 1 23 45

2

2

67

9

8

0

ACS Paragon Plus Environment

2.5

3.0 r (˚ A)

3.5

0 4.0

hn(r)i

10

Page 19 of 24

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 20 of 24

kT lnpBa2+ (nλ = n)

Page6 21 of 24The Journal 2.9 ˚ A of Physical Chemistry

3.0 3.1 3.2

4

1 2 2 3 4 5 0 6 7 −2 8 9 −4 10 11 12−6 0 13

˚ A ˚ A ˚ A

ACS Paragon Plus Environment

2

4

n

6

8

10

−200

The Journal of Physical ChemistryPage 22 of 24

λ = 3.5 ˚ A −225

(ex)

µBa2+ (kcal/mol)

1 2 3 4 −250 5 6 7 −275 8 9 10 −300 11 12 1 13

λ = 2.9 ˚ A ACS Paragon Plus Environment

2

3

4

n

5

6

7

8

0

Page 23 of 24The Journal of Physicalassociation Chemistry free energy

kcal/mol

1 2 3 −100 4 5 6 7 8 9 −200 10 11 12 13 14 15−300 16 1 17

with ligand density outer − shell hydration free energy

ACS Paragon Plus Environment

2

3

4

n

5

6

7

8

−200

The Journal of Physical ChemistryPage 24 of 24

λ = 3.5 ˚ A −225

(ex)

µBa2+ (kcal/mol)

1 2 3 4 −250 5 6 7 −275 8 9 10 −300 11 12 1 13

λ = 2.9 ˚ A ACS Paragon Plus Environment

2

3

4

n

5

6

7

8