Origin of Asymmetric Solvation Effects for Ions in Water and Organic Solvents Investigated Using Molecular Dynamics Simulations: The Swain Acity-Basity Scale Revisited

A charge-dependent long-ranged force drives tailored assembly of matter in solution

The interaction between charged objects in solution is generally expected to recapitulate two central principles of electromagnetics: (1) like-charged objects repel, and (2) they do so regardless of the sign of their electrical charge. Here we demonstrate experimentally that the solvent plays a hitherto unforeseen but crucial role in interparticle interactions, and importantly, that interactions in the fluid phase can break charge-reversal symmetry. We show that in aqueous solution, negatively charged particles can attract at long range while positively charged particles repel. In solvents that exhibit an inversion of the net molecular dipole at an interface, such as alcohols, we find that the converse can be true: positively charged particles may attract whereas negatives repel. The observations hold across a wide variety of surface chemistries: from inorganic silica and polymeric particles to polyelectrolyte- and polypeptide-coated surfaces in aqueous solution. A theory of interparticle interactions that invokes solvent structuring at an interface captures the observations. Our study establishes a nanoscopic interfacial mechanism by which solvent molecules may give rise to a strong and long-ranged force in solution, with immediate ramifications for a range of particulate and molecular processes across length scales such as self-assembly, gelation and crystallization, biomolecular condensation, coacervation, and phase segregation.

Efficient Computation of the Electrostatic Component of Solvation Free Energy via a Two-Point Padé Approximation

We develop an efficient method to compute the electrostatic component of the solvation free energy via the two-point Padé approximation. The Padé approximant uses four parameters to describe the electrostatic free energy change of the solvation process, which could be readily determined from four thermodynamic properties obtained in two simulations, namely, the first- and second-order free energy gradients of any two states. Therefore, instead of sampling at multiple intermediate states, only two states, e.g., electrostatically fully solvated and desolvated, are needed to determine the Padé approximant and compute the corresponding free energy contribution. Applications to several model systems, including both neutral and charged species, show that the method can accurately produce electrostatic solvation free energy. The method would be very useful to save computational cost in applications in which accurate but expensive energy functions like quantum mechanics are used.

Solvation Thermodynamics of Solutes in Water and Ionic Liquids Using the Multiscale Solvation-Layer Interface Condition Continuum Model

Molecular assembly processes are generally driven by thermodynamic properties in solutions. Atomistic modeling can be very helpful in designing and understanding complex systems, except that bulk solvent is very inefficient to treat explicitly as discrete molecules. In this work, we develop and assess two multiscale solvation models for computing solvation thermodynamic properties. The new SLIC/CDC model combines continuum solvent electrostatics based on the solvent layer interface condition (SLIC) with new statistical thermodynamic models for hydrogen bonding and nonpolar modes: cavity formation, dispersion interactions, combinatorial mixing (CDC). Given the structures of 500 solutes, the SLIC/CDC model predicts Gibbs energies of solvation in water with an average accuracy better than 1 kcal/mol, when compared to experimental measurements, and better than 0.8 kcal/mol, when compared to explicit-solvent molecular dynamics simulations. The individual SLIC/CDC energy mode values agree quantitatively with those computed from explicit-solvent molecular dynamics. The previously published SLIC/SASA multiscale model combines the SLIC continuum electrostatic model with the solvent-accessible surface area (SASA) nonpolar energy mode. With our new, improved parametrization method, the SLIC/SASA model now predicts Gibbs energies of solvation with better than 1.4 kcal/mol average accuracy in aqueous systems, compared to experimental and explicit-solvent molecular dynamics, and better than 1.6 kcal/mol average accuracy in ionic liquids, compared to explicit-solvent molecular dynamics. Both models predict solvation entropies, and are the first implicit-solvation models capable of predicting solvation heat capacities.

Measurements and Utilization of Consistent Gibbs Energies of Transfer of Single Ions: Towards a Unified Redox Potential Scale for All Solvents

Utilizing the "ideal" ionic liquid salt bridge to measure Gibbs energies of transfer of silver ions between the solvents water, acetonitrile, propylene carbonate and dimethylformamide results in a consistent data set with a precision of 0.6 kJ mol-1 over 87 measurements in 10 half-cells. This forms the basis for a coherent experimental thermodynamic framework of ion solvation chemistry. In addition, we define the solvent independent pe abs H 2 O - and the E abs H 2 O values that account for the electronating potential of any redox system similar to the pH abs H 2 O value of a medium that accounts for its protonating potential. This E abs H 2 O scale is thermodynamically well-defined enabling a straightforward comparison of the redox potentials (reducities) of all media with respect to the aqueous redox potential scale, hence unifying all conventional solvents' redox potential scales. Thus, using the Gibbs energy of transfer of the silver ion published herein, one can convert and unify all hitherto published redox potentials measured, for example, against ferrocene, to the E abs H 2 O scale.

The Role of Surface Chemistry in the Orientational Behavior of Water at an Interface

Molecular dynamics studies have demonstrated that molecular water at an interface, with either a gas or a solid, displays anisotropic orientational behavior in contrast to its bulk counterpart. This effect has been recently implicated in the like-charge attraction problem for colloidal particles in solution. Here, negatively charged particles in solution display a long-ranged attraction where continuum electrostatic theory predicts monotonically repulsive interactions, particularly in solutions with monovalent salt ions at low ionic strength. Anisotropic orientational behavior of solvent molecules at an interface gives rise to an excess interfacial electrical potential which we suggest generates an additional solvation contribution to the total free energy that is traditionally overlooked in continuum descriptions of interparticle interactions in solution. In the present investigation we perform molecular dynamics simulation based calculations of the interfacial potential using realistic surface models representing various chemistries as well as different solvents. Similar to previous work that focused on simple model surfaces constructed by using oxygen atoms, we find that solvents at more realistic model surfaces exhibit substantial anisotropic orientational behavior. We explore the dependence of the interfacial solvation potential on surface properties such as surface group chemistry and group density at silica and carboxylated polystyrene interfaces. For water, we note surprisingly good agreement between results obtained for a simple O-atom wall and more complex surface models, suggesting a general qualitative consistency of the interfacial solvation effect for surfaces in contact with water. In contrast, for an aprotic solvent such as DMSO, surface chemistry appears to exert a stronger influence on the sign and magnitude of the interfacial solvation potential. The study carries broad implications for molecular-scale interactions and may find relevance in explaining a range of phenomena in soft-matter physics and cell biology.

Nonmonotonic Pair Potentials in the Interaction of Like-Charged Objects in Solution

We consider the long-standing like-charge attraction problem, wherein under certain conditions, similarly charged spheres suspended in aqueous electrolyte have been observed to display a minimum in their interaction potential, contrary to the intuitively expected monotonically varying repulsion. Recently, we described an interfacial mechanism invoking the molecular nature of the solvent that explains this anomalous experimental observation. In our model for the interaction of negatively charged particles in water, the minimum in the pair potential results from the superposition of competing contributions to the total free energy. One of these contributions is the canonical repulsive electrostatic term, whereas the other is a solvation-induced attractive contribution. We find that whereas both contributions grow approximately exponentially with decreasing interparticle separation, the occurrence of a stable, long-ranged minimum in the pair potential arises from differences in the precise interparticle separation dependence of the two terms. Specifically, the interfacial solvation term exhibits a more gradual decay with distance than the electrostatic repulsion, permitting the attractive contribution to dominate the interaction at large distances. Importantly, these disparities become evident in quantities calculated from exact numerical solutions to the governing nonlinear Poisson-Boltzmann (PB) equation for the spatial electrical potential distribution in the system. In marked contrast, we find that the linearized PB equation, applicable in the regime of low surface electrical potentials, does not support nonmonotonic trends in the total interaction free energy within the present model. Our results point to the importance of exact descriptions of electrostatic interactions in real systems that most often do not subscribe to particular mathematical limits where analytical approximations may provide a sufficiently accurate description of the problem.

The effect of solvent on the binding of anions and ion-pairs with a neutral [2]rotaxane

In this work we report the binding properties of rotaxane 1 towards a series of tetraalkylammonium salts of Cl^-, OCN^- and NO₃^- anions in acetone and a CHCl₃/MeOH solvent mixture. We used ¹H NMR titrations and Isothermal Titration Calorimetry (ITC) experiments to monitor and analyze the binding processes. We compared the obtained results with those previously described by us in chloroform solution. In acetone solution, the determined binding constants for the 1 : 1 complexes were 1 to 3 orders of magnitude larger than those measured in chloroform, a less competitive solvent for hydrogen-bonding. The thermodynamic signatures of the binding processes in acetone, determined by ITC experiments, revealed favorable enthalpic and entropic contributions having similar magnitudes. These results suggested that solvation/desolvation processes in acetone play a significant role in the binding processes. Conversely, the addition of just 5% of methanol to chloroform solutions of 1 significantly reduces the magnitude of the binding constants of all studied ion-pairs. In this solvent mixture, the entropy term is also favorable but it does not compensate the experienced loss of binding enthalpy. Moreover, in acetone solution, the addition of the Cl^- and OCN^- tetraalkylammonium salts in excess (more than 1 equiv.) led to the immediate appearance of 2 : 1 complexes. Related high-stoichiometry complexes are not observed in the solvent mixture (CHCl₃/MeOH 95/5). In chloroform, a large excess of the salt (> 6 equiv.) is required for its formation. From the analysis of the obtained binding data we infer that, in acetone and in CHCl₃/MeOH mixture, the formed complexes are mainly anionic.

Efficient Alchemical Intermediate States in Free Energy Calculations Using λ-Enveloping Distribution Sampling

Alchemical free energy calculations generally require intermediate states along a coupling parameter λ to establish sufficient phase space overlap for obtaining converged results. Such intermediate states can also be engineered to lower the energy barriers and, consequently, reduce the required sampling time. The recently introduced λ-enveloping distribution sampling (λ-EDS) scheme combines the properties of the minimum variance pathway and the EDS methods to improve sampling and allow for larger steps along the alchemical pathway compared to conventional approaches. This scheme also eliminates the need for soft-core potentials and retains the behavior of conventional λ-intermediate states as a limiting case. In this study, an automated procedure is developed to select the parameters of λ-EDS for optimal performance. The underlying theory is illustrated based on simulations of simple test systems (bond length changes in harmonic oscillators, mutations of dihedral angles, and charge creation in water), as well as on the calculation of the absolute hydration free energies of 12 small organic molecules.

Enhancing the promiscuity of a member of the Caspase protease family by rational design

The N-terminal cleavage of fusion tags to restore the native N-terminus of recombinant proteins is a challenging task and up to today, protocols need to be optimized for different proteins individually. Within this work, we present a novel protease that was designed in-silico to yield enhanced promiscuity toward different N-terminal amino acids. Two mutations in the active-site amino acids of human Caspase-2 were determined to increase the recognition of branched amino-acids, which show only poor binding capabilities in the unmutated protease. These mutations were determined by sequential and structural comparisons of Caspase-2 and Caspase-3 and their effect was additionally predicted using free-energy calculations. The two mutants proposed in the in-silico studies were expressed and in-vitro experiments confirmed the simulation results. Both mutants showed not only enhanced activities toward branched amino acids, but also smaller, unbranched amino acids. We believe that the created mutants constitute an important step toward generalized procedures to restore original N-termini of recombinant fusion proteins.

Hybrid Solvation Model with First Solvation Shell for Calculation of Solvation Free Energy

We present a hybrid solvation model with first solvation shell to calculate solvation free energies. This hybrid model combines the quantum mechanics and molecular mechanics methods with the analytical expression based on the Born solvation model to calculate solvation free energies. Based on calculated free energies of solvation and reaction profiles in gas phase, we set up a unified scheme to predict reaction profiles in solution. The predicted solvation free energies and reaction barriers are compared with experimental results for twenty bimolecular nucleophilic substitution reactions. These comparisons show that our hybrid solvation model can predict reliable solvation free energies and reaction barriers for chemical reactions of small molecules in aqueous solution.

Interfacial solvation can explain attraction between like-charged objects in aqueous solution

Over the past few decades, the experimental literature has consistently reported observations of attraction between like-charged colloidal particles and macromolecules in aqueous solution. Examples include nucleic acids and colloidal particles in the bulk solution and under confinement, and biological liquid-liquid phase separation. This observation is at odds with the intuitive expectation of an interparticle repulsion that decays monotonically with distance. Although attraction between like-charged particles can be rationalized theoretically in the strong-coupling regime, e.g., in the presence of multivalent counterions, recurring accounts of long-range attraction in aqueous solution containing monovalent ions at low ionic strength have posed an open conundrum. Here, we show that the behavior of molecular water at an interface-traditionally disregarded in the continuum electrostatics picture-provides a mechanism to explain the attraction between like-charged objects in a broad spectrum of experiments. This basic principle will have important ramifications in the ongoing quest to better understand intermolecular interactions in solution.

The Born model can accurately describe electrostatic ion solvation

The solvation free energies of ions in water are consistent with the Born linear response model if the centre on which the ion–water repulsion force acts is moved from the oxygen atom towards the hydrogens.

Single-Ion Thermodynamics from First Principles: Calculation of the Absolute Hydration Free Energy and Single-Electrode Potential of Aqueous Li+ Using ab Initio Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations

A recently proposed thermodynamic integration (TI) approach formulated in the framework of quantum mechanical/molecular mechanical molecular dynamics (QM/MM MD) simulations is applied to study the structure, dynamics, and absolute intrinsic hydration free energy Δ_s G_M⁺,wat^◦ of the Li⁺ ion at a correlated ab initio level of theory. Based on the results, standard values (298.15 K, ideal gas at 1 bar, ideal solute at 1 molal) for the absolute intrinsic hydration free energy [Formula: see text] of the proton, the surface electric potential jump χ_wat^◦ upon entering bulk water, and the absolute single-electrode potential [Formula: see text] of the reference hydrogen electrode are calculated to be -1099.9 ± 4.2 kJ·mol^-1, 0.13 ± 0.08 V, and 4.28 ± 0.04 V, respectively, in excellent agreement with the standard values recommended by Hünenberger and Reif on the basis of an extensive evaluation of the available experimental data (-1100 ± 5 kJ·mol^-1, 0.13 ± 0.10 V, and 4.28 ± 0.13 V). The simulation results for Li⁺ are also compared to those for Na⁺ and K⁺ from a previous study in terms of relative hydration free energies ΔΔ_s G_M⁺,wat^◦ and relative electrode potentials [Formula: see text]. The calculated values are found to agree extremely well with the experimental differences in standard conventional hydration free energies ΔΔ_s G_M⁺,wat^• and redox potentials [Formula: see text]. The level of agreement between simulation and experiment, which is quantitative within error bars, underlines the substantial accuracy improvement achieved by applying a highly demanding QM/MM approach at the resolution-of-identity second-order Møller-Plesset perturbation (RIMP2) level over calculations relying on purely molecular mechanical or density functional theory (DFT) descriptions. A detailed analysis of the structural and dynamical properties of the Li⁺ hydrate indicates that a correct description of the solvation structure and dynamics is achieved as well at this level of theory. Consideration of the QM/MM potential-energy components also shows that the partitioning into QM and MM zones does not induce any significant energetic artifact for the system considered.

Interfacial structural crossover and hydration thermodynamics of charged C60 in water

Classical molecular dynamics simulations of the hydration thermodynamics, structure, and dynamics of water in hydration shells of charged buckminsterfullerenes are presented in this study. Charging of fullerenes leads to a structural transition in the hydration shell, accompanied by creation of a significant population of dangling O-H bonds pointing toward the solute. In contrast to the well accepted structure-function paradigm, this interfacial structural transition causes nearly no effect on either the dynamics of hydration water or on the solvation thermodynamics. Linear response to the solute charge is maintained despite significant structural changes in the hydration shell, and solvation thermodynamic potentials are nearly insensitive to the altering structure. Only solvation heat capacities, which are higher thermodynamic derivatives of the solvation free energy, indicate some sensitivity to the local hydration structure. We have separated the solvation thermodynamic potentials into direct solute-solvent interactions and restructuring of the hydration shell and analyzed the relative contributions of electrostatic and nonpolar interactions to the solvation thermodynamics.

Understanding the scale of the single ion free energy: A critical test of the tetra-phenyl arsonium and tetra-phenyl borate assumption

The tetra-phenyl arsonium and tetra-phenyl borate (TATB) assumption is a commonly used extra-thermodynamic assumption that allows single ion free energies to be split into cationic and anionic contributions. The assumption is that the values for the TATB salt can be divided equally. This is justified by arguing that these large hydrophobic ions will cause a symmetric response in water. Experimental and classical simulation work has raised potential flaws with this assumption, indicating that hydrogen bonding with the phenyl ring may favor the solvation of the TB^- anion. Here, we perform ab initio molecular dynamics simulations of these ions in bulk water demonstrating that there are significant structural differences. We quantify our findings by reproducing the experimentally observed vibrational shift for the TB^- anion and confirm that this is associated with hydrogen bonding with the phenyl rings. Finally, we demonstrate that this results in a substantial energetic preference of the water to solvate the anion. Our results suggest that the validity of the TATB assumption, which is still widely used today, should be reconsidered experimentally in order to properly reference single ion solvation free energy, enthalpy, and entropy.

Basic Remarks on Acidity

This Review provides a unified view on Brønsted acidity. For this purpose, a brief overview of the concepts acidity, acid strengths, and pH value is given, including problems, proposed solutions, and the use of the pH_abs /pHabsH2O scale as a unifying concept. Thereafter, some examples of the accessibility and application of unified pH_abs values are given. The Review is rounded off with the analogy of acid-base chemistry to redox chemistry with the introduction of the unified redox scale pe_abs . The combination of pH_abs and pe_abs values in the protoelectric potential map (PPM), as elaborated in ongoing studies on the thermochemistry of single ions, provides a means to classify and to compare all possible acid-base/redox reactions in a medium-independent and, thus, unified fashion.

Effect of Sodium and Chloride Binding on a Lecithin Bilayer. A Molecular Dynamics Study

The effect of ion binding on the structural, mechanical, dynamic and electrostatic properties of a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer in a 0.5 M aqueous NaCl solution is investigated using classical atomistic molecular dynamics simulation with different force-field descriptions for ion-ion and ion-lipid interactions. Most importantly, the repulsive Lennard-Jones parameters for the latter were modified, such that approximately similar binding of cations and anions to the lipid membrane is achieved. This was done to qualitatively improve the apparent ion-lipid binding constants obtained from simulations with the original force field (Berger lipids and GROMOS87 ions in combination with the SPC water model) in comparison to experimental data. Furthermore, various parameters characterizing membrane structure, elasticity, order and dynamics are analyzed. It is found that ion binding as observed in simulations involving the modified in comparison to the original force-field description leads to: (i) a smaller salt-induced change in the area per lipid, which is in closer agreement with the experiment; (ii) a decrease in the area compressibility and bilayer thickness to values comparable to a bilayer in pure water; (iii) lipid deuterium order parameters and lipid diffusion coefficients on nanosecond timescales that are very similar to the values for a membrane in pure water. In general, salt effects on the structural properties of a POPC bilayer in an aqueous sodium-chloride solution appear to be reproduced reasonably well by the new force-field description. An analysis of membrane-membrane disjoining pressure suggests that the smaller salt-induced change in area per lipid induced by the new force-field description is not due to the alteration of membrane-associated net charge, but must rather be understood as a consequence of ion-specific effects on the arrangement of lipid molecules.

The influence of hydrogen bonding on partition coefficients

This Perspective explores how consideration of hydrogen bonding can be used to both predict and better understand partition coefficients. It is shown how polarity of both compounds and substructures can be estimated from measured alkane/water partition coefficients. When polarity is defined in this manner, hydrogen bond donors are typically less polar than hydrogen bond acceptors. Analysis of alkane/water partition coefficients in conjunction with molecular electrostatic potential calculations suggests that aromatic chloro substituents may be less lipophilic than is generally believed and that some of the effect of chloro-substitution stems from making the aromatic π-cloud less available to hydrogen bond donors. Relationships between polarity and calculated hydrogen bond basicity are derived for aromatic nitrogen and carbonyl oxygen. Aligned hydrogen bond acceptors appear to present special challenges for prediction of alkane/water partition coefficients and this may reflect 'frustration' of solvation resulting from overlapping hydration spheres. It is also shown how calculated hydrogen bond basicity can be used to model the effect of aromatic aza-substitution on octanol/water partition coefficients.

Accuracy limit of rigid 3-point water models

Classical 3-point rigid water models are most widely used due to their computational efficiency. Recently, we introduced a new approach to constructing classical rigid water models [S. Izadi et al., J. Phys. Chem. Lett. 5, 3863 (2014)], which permits a virtually exhaustive search for globally optimal model parameters in the sub-space that is most relevant to the electrostatic properties of the water molecule in liquid phase. Here we apply the approach to develop a 3-point Optimal Point Charge (OPC3) water model. OPC3 is significantly more accurate than the commonly used water models of same class (TIP3P and SPCE) in reproducing a comprehensive set of liquid bulk properties, over a wide range of temperatures. Beyond bulk properties, we show that OPC3 predicts the intrinsic charge hydration asymmetry (CHA) of water - a characteristic dependence of hydration free energy on the sign of the solute charge - in very close agreement with experiment. Two other recent 3-point rigid water models, TIP3PFB and H2ODC, each developed by its own, completely different optimization method, approach the global accuracy optimum represented by OPC3 in both the parameter space and accuracy of bulk properties. Thus, we argue that an accuracy limit of practical 3-point rigid non-polarizable models has effectively been reached; remaining accuracy issues are discussed.