Ionic force field optimization based on single-ion and ion-pair solvation properties.

Molecular dynamics simulations of ionic solutions depend sensitively on the force fields employed for the ions. To resolve the fine differences between ions of the same valence and roughly similar size and in particular to correctly describe ion-specific effects, it is clear that accurate force fields are necessary. In the past, optimization strategies for ionic force fields either considered single-ion properties (such as the solvation free energy at infinite dilution or the ion-water structure) or ion-pair properties (in the form of ion-ion distribution functions). In this paper we investigate strategies to optimize ionic force fields based on single-ion and ion-pair properties simultaneously. To that end, we simulate five different salt solutions, namely, CsCl, KCl, NaI, KF, and CsI, at finite ion concentration. The force fields of these ions are systematically varied under the constraint that the single-ion solvation free energy matches the experimental value, which reduces the two-dimensional {sigma,epsilon} parameter space of the Lennard-Jones interaction to a one dimensional line for each ion. From the finite-concentration simulations, the pair potential is extracted and the osmotic coefficient is calculated, which is compared to experimental data. We find a strong dependence of the osmotic coefficient on the force field, which is remarkable as the single-ion solvation free energy and the ion-water structure remain invariant under the parameter variation. Optimization of the force field is achieved for the cations Cs(+) and K(+), while for the anions I(-) and F(-) the experimental osmotic coefficient cannot be reached. This suggests that in the long run, additional parameters might have to be introduced into the modeling, for example by modified mixing rules.

Ion specific effects of sodium and potassium on the catalytic activity of HIV-1 protease.

The behavior of HIV-1 protease in aqueous NaCl and KCl solutions is investigated by kinetic measurements and molecular dynamics simulations. Experiments show cation-specific effects on the enzymatic activity. The initial velocity of peptide substrate hydrolysis increases with salt concentration more dramatically in potassium than in sodium chloride solutions. Furthermore, significantly higher catalytic efficiencies (k(cat)/K(M)) are observed in the presence of K+ compared to Na+ at comparable salt concentrations. Molecular dynamics simulations provide insight into this ion-specific behavior. Sodium is attracted more strongly than potassium to the protein surface primarily due to stronger interactions with carboxylate side chain groups of aspartates and glutamates. These effects are of particular importance for acidic amino acid residues at or near the active site of the enzyme, including a pair of aspartates at the entrance to the reaction cavity. We infer that the presence of more Na+ than K+ at the active site leads to a lower increase in enzymatic activity with increasing salt concentration in the presence of Na+, likely due to the ability of the alkali cations at the active site to lower the efficiency of substrate binding.

Specificity of ion-protein interactions: complementary and competitive effects of tetrapropylammonium, guanidinium, sulfate, and chloride ions.

The interactions of ions with a model peptide (a single melittin alpha-helix) in solutions of tetrapropylammonium sulfate or guanidinium chloride were examined by molecular dynamics simulations. The tetrapropylammonium cation shares the geometrical property of essentially flat faces with the previously examined guanidinium cation, and it was found that that this geometry leads to a strong preference for tetrapropylammonium to interact in a similar stacking-type fashion with flat nonpolar groups such as the indole side chain of tryptophan. In contrast to guanidinium, however, tetrapropylammonium does not exhibit strong ion pairing or clustering with sulfate counterions in the solution. Sulfate was found to interact almost exclusively and strongly with the cationic groups of the peptide, such that, already in a 0.1 m solution of tetrapropylammonium sulfate, the 6+ charge of the peptide is effectively locally neutralized. In combination with previous simulations, neutron scattering studies, and experiments on the conformational stability of model peptides, the present results suggest that the Hofmeister series can be explained in higher detail by splitting ions according to the effect they have on hydrogen bonding, salt bridges, and hydrophobic interactions in the protein and how these effects are altered by the counterion.

Ion-specific thermodynamics of multicomponent electrolytes: a hybrid HNC/MD approach.

Using effective infinite dilution ion-ion interaction potentials derived from explicit-water molecular dynamics (MD) computer simulations in the hypernetted-chain (HNC) integral equation theory we calculate the liquid structure and thermodynamic properties, namely, the activity and osmotic coefficients of various multicomponent aqueous electrolyte mixtures. The electrolyte structure expressed by the ion-ion radial distribution functions is for most ions in excellent agreement with MD and implicit solvent Monte Carlo (MC) simulation results. Calculated thermodynamic properties are also represented consistently among these three methods. Our versatile HNC/MD hybrid method allows for a quick prediction of the thermodynamics of multicomponent electrolyte solutions for a wide range of concentrations and an efficient assessment of the validity of the employed MD force-fields with possible implications in the development of thermodynamically consistent parameter sets.

Specific ion binding to nonpolar surface patches of proteins.

Employing detailed atomistic modeling we study the mechanisms behind ion binding to proteins and other biomolecules and conclude that (1) small, hard ions bind via direct ion pairing to charged surface groups and (2) large, soft ions bind to nonpolar groups via a solvent assisted attraction. Our predictions are in qualitative agreement with bulk solution data and may provide an important clue for the basic understanding of ion-specific effects in biological systems.

Quantification and rationalization of the higher affinity of sodium over potassium to protein surfaces.

For a series of different proteins, including a structural protein, enzyme, inhibitor, protein marker, and a charge-transfer system, we have quantified the higher affinity of Na+ over K+ to the protein surface by means of molecular dynamics simulations and conductivity measurements. Both approaches show that sodium binds at least twice as strongly to the protein surface than potassium does with this effect being present in all proteins under study. Different parts of the protein exterior are responsible to a varying degree for the higher surface affinity of sodium, with the charged carboxylic groups of aspartate and glutamate playing the most important role. Therefore, local ion pairing is the key to the surface preference of sodium over potassium, which is further demonstrated and quantified by simulations of glutamate and aspartate in the form of isolated amino acids as well as short oligopeptides. As a matter of fact, the effect is already present at the level of preferential pairing of the smallest carboxylate anions, formate or acetate, with Na+ versus K+, as shown by molecular dynamics and ab initio quantum chemical calculations. By quantifying and rationalizing the higher preference of sodium over potassium to protein surfaces, the present study opens a way to molecular understanding of many ion-specific (Hofmeister) phenomena involving protein interactions in salt solutions.

Homogeneous freezing of water starts in the subsurface.

Molecular dynamics simulations of homogeneous ice nucleation in extended aqueous slabs show that freezing preferentially starts in the subsurface. The top surface layer remains disordered during the freezing process. The subsurface accommodates better than the bulk the increase of volume connected with freezing. It also experiences strong electric fields caused by oriented surface water molecules, which can enhance ice nucleation. Our computational results shed new light on the experimental controversy concerning the bulk vs surface origin of homogeneous ice nucleation in water droplets. This has important atmospheric implications for the microphysics of formation of high altitude clouds.

Selected biologically relevant ions at the air/water interface: a comparative molecular dynamics study.

Interfacial behavior of selected biologically and technologically relevant ions is studied using molecular dynamics simulations employing polarizable potentials. Propensities of choline, tetraalkylammonium (TAA), and sodium cations, and sulfate and chloride anions for the air/water interface are analyzed by means of density profiles. Affinity of TAA ions for the interface increases with their increasing hydrophobicity. Tetramethylammonium favors bulk solvation, whereas cations with propyl and butyl chains behave as surfactants. The choice of counter-anions has only a weak effect on the behavior of these cations. For choline, sodium, chloride and sulfate, the behavior at the air/water interface was compared to the results of our recent study on the segregation of these ions at protein surfaces. No analogy between these two interfaces in terms of ion segregation is found.

Specific ion effects at protein surfaces: a molecular dynamics study of bovine pancreatic trypsin inhibitor and horseradish peroxidase in selected salt solutions.

The distribution of sodium, choline, sulfate, and chloride ions around two proteins, horseradish peroxidase (HRP) and bovine pancreatic trypsin inhibitor (BPTI), is investigated by means of molecular dynamics simulations with the aim to elucidate ion adsorption at the protein surface. Although the two proteins under investigation are very different from each other, the ion distributions around them are remarkably similar. Sulfate is always strongly attached to the proteins, choline shows a significant, but unspecific, propensity for the protein surfaces, and sodium ions have a weak surface affinity, while chloride has virtually no preference for the protein surface. In mixtures of all four ion species in protein solutions, the resulting distributions are almost a superposition of the distributions of sodium sulfate and choline chloride, except that sodium partially replaces choline close to the proteins. The present simulations support a picture of ions interacting with individual ionic and polar amino acid groups rather than with an averaged protein surface. The results thus show how subtle the so-called Hofmeister and electroselectivity effects are in salt solution of proteins, making all simplified interaction models questionable.

Aqueous ionic and complementary zwitterionic soluble surfactants: molecular dynamics simulations and sum frequency generation spectroscopy of the surfaces.

An aqueous ionic surfactant, 1-dodecyl-4-(dimethylamino)pyridinium (DMP) bromide, and the corresponding zwitterion 2-[4-(dimethylamino)pyridinio]dodecanoate (DPN) were explored by means of molecular dynamics (MD) simulations and, for the ionic system, by infrared-visible sum frequency generation (IR-vis SFG). The molecular structure of the interfacial layer was investigated for the ionic and zwitterionic systems as a function of surfactant concentration, both in water and in salt (KF or KBr) solutions, by MD simulations in a slab geometry. The buildup of the surface monolayer and a sublayer was monitored, and density and orientational profiles of the surfactants were evaluated. The difference between the ionic and zwitterionic systems and the effect of the added salt were analyzed at the molecular level. The results of MD simulations were compared to those of nonlinear optical spectroscopy measurements. IR-vis SFG was employed to study the DMP ionic surfactant in water and upon addition of simple salts. The influence of added salts on the different molecular moieties at the interface was quantified in detail experimentally.

Brine rejection from freezing salt solutions: a molecular dynamics study.

The atmospherically and technologically very important process of brine rejection from freezing salt solutions is investigated with atomic resolution using molecular dynamics simulations. The present calculations allow us to follow the motion of each water molecule and salt ion and to propose a microscopic mechanism of brine rejection, in which a fluctuation (reduction) of the ion density in the vicinity of the ice front is followed by the growth of a new ice layer. The presence of salt slows down the freezing process, which leads to the formation of an almost neat ice next to a disordered brine layer.

Interior and interfacial aqueous solvation of benzene dicarboxylate dianions and their methylated analogues: A combined molecular dynamics and photoelectron spectroscopy study.

Aqueous solvation of benzene dicarboxylate dianions (BCD(2-)) was studied by means of photoelectron spectroscopy and molecular dynamics simulations. Photoelectron spectra of hydrated o- and p-BCD(2-) with up to 25 water molecules were obtained. An even-odd effect was observed for the p-BCD(2-) system as a result of the alternate solvation of the two negative charges. However, the high polarizability of the benzene ring makes the two carboxylate groups interact with each other in p-BCD(2-), suppressing the strength of this even-odd effect compared with the linear dicarboxylate dianions linked by an aliphatic chain. No even-odd effect was observed for the o-BCD(2-) system, because each solvent molecule can interact with the two carboxylate groups at the same time due to their proximity. For large solvated clusters, the spectral features of the solute decreased while the solvent features became dominant, suggesting that both o- and p-BCD(2-) are situated in the center of the solvated clusters. Molecular dynamics simulations with both nonpolarizable and polarizable force fields confirmed that all three isomers (o-, m-, and p-BCD(2-)) solvate in the aqueous bulk. However, upon methylation the hydrophobic forces overwhelm electrostatic interactions and, as a result, the calculations predict that the tetramethyl-o-BCD(2-) is located at the water surface with the carboxylate groups anchored in the liquid and the methylated benzene ring tilted away from the aqueous phase.