Pavel Jungwirth

Hofmeister series and specific interactions of charged headgroups with aqueous ions

In this paper, we propose a Hofmeister-like ordering of charged headgroups. To this purpose we review various literature data and complete them with some new experimental and computational results on interactions of ions with alkyl sulfates and carboxylates. We further combine the proposed headgroup ordering with the law of matching water affinities in order to obtain a general description and predictions of ion-headgroup interactions. Examples from colloidal chemistry and from biological systems are provided to illustrate the power of this approach.

Water surface is acidic

Water autoionization reaction 2H2O → H3O− + OH− is a textbook process of basic importance, resulting in pH = 7 for pure water. However, pH of pure water surface is shown to be significantly lower, the reduction being caused by proton stabilization at the surface. The evidence presented here includes ab initio and classical molecular dynamics simulations of water slabs with solvated H3O+ and OH− ions, density functional studies of (H2O)48H+ clusters, and spectroscopic isotopic-exchange data for D2O substitutional impurities at the surface and in the interior of ice nanocrystals. Because H3O+ does, but OH− does not, display preference for surface sites, the H2O surface is predicted to be acidic with pH < 4.8. For similar reasons, the strength of some weak acids, such as carbonic acid, is expected to increase at the surface. Enhanced surface acidity can have a significant impact on aqueous surface chemistry, e.g., in the atmosphere.

Biomolecular simulations of membranes: Physical properties from different force fields

Phospholipid force fields are of ample importance for the simulation of artificial bilayers, membranes, and also for the simulation of integral membrane proteins. Here, we compare the two most applied atomic force fields for phospholipids, the all-atom CHARMM27 and the united atom Berger force field, with a newly developed all-atom generalized AMBER force field (GAFF) for dioleoylphosphatidylcholine molecules. Only the latter displays the experimentally observed difference in the order of the C2 atom between the two acyl chains. The interfacial water dynamics is smoothly increased between the lipid carbonyl region and the bulk water phase for all force fields; however, the water order and with it the electrostatic potential across the bilayer showed distinct differences between the force fields. Both Berger and GAFF underestimate the lipid self-diffusion. GAFF offers a consistent force field for the atomic scale simulation of biomembranes.

Molecular Mechanisms of Ion-Specific Effects on Proteins

The specific binding sites of Hofmeister ions with an uncharged 600-residue elastin-like polypeptide, (VPGVG)(120), were elucidated using a combination of NMR and thermodynamic measurements along with molecular dynamics simulations. It was found that the large soft anions such as SCN(-) and I(-) interact with the polypeptide backbone via a hybrid binding site that consists of the amide nitrogen and the adjacent α-carbon. The hydrocarbon groups at these sites bear a slight positive charge, which enhances anion binding without disrupting specific hydrogen bonds to water molecules. The hydrophobic side chains do not contribute significantly to anion binding or the corresponding salting-in behavior of the biopolymer. Cl(-) binds far more weakly to the amide nitrogen/α-carbon binding site, while SO(4)(2-) is repelled from both the backbone and hydrophobic side chains of the polypeptide. The Na(+) counterions are also repelled from the polypeptide. The identification of these molecular-level binding sites provides new insights into the mechanism of peptide-anion interactions.

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.

The orientation and charge of water at the hydrophobic oil droplet-water interface.

We established the charge and structure of the oil/water interface by combining ζ-potential measurements, sum frequency scattering (SFS) and molecular dynamics simulations. The SFS experiments show that the orientation of water molecules can be followed on the oil droplet/water interface. The average water orientation on a neat oil droplet/water interface is the same as the water orientation on a negatively charged interface. pH dependent experiments show, however, that there is no sign of selective adsorption of hydroxide ions. Molecular dynamics simulations, both with and without intermolecular charge transfer, show that the balance of accepting and donating hydrogen bonds is broken in the interfacial layer, leading to surface charging. This can account for the negative surface charge that is found in experiments.

Effects of Alkali Cations and Halide Anions on the DOPC Lipid Membrane

By means of molecular dynamics simulations with an all-atom force field, we investigated the affinities of alkali cations and halide anions for the dioleoylphosphatidylcholine lipid membrane in aqueous salt solutions. In addition, changes in phospholipid lateral diffusion and in headgroup mobility upon adding NaCl were observed using fluorescence spectroscopy. The simulations revealed that sodium is attracted to the headgroup region with its concentration being maximal in the vicinity of the phosphate groups. Potassium and cesium, however, do not preferentially adsorb to the membrane. Similarly, halide anions do not exhibit a strong affinity for the lipid headgroups but merely compensate for the positive charge of the sodium countercations. Nevertheless, larger halides such as bromide and iodide penetrate deeper into the headgroup region toward the boundary with the hydrophobic alkyl chain, this effect being likely underestimated within the present nonpolarizable force field. Addition of alkali halide salts modifies physical properties of the bilayer including the electronic density profiles, the electrostatic potential, and the area per lipid headgroup.

Cation-Specific Interactions with Carboxylate in Amino Acid and Acetate Aqueous Solutions: X-ray Absorption and ab initio Calculations

Relative interaction strengths between cations (X = Li (+), Na (+), K (+), NH 4 (+)) and anionic carboxylate groups of acetate and glycine in aqueous solution are determined. These model systems mimic ion pairing of biologically relevant cations with negatively charged groups at protein surfaces. With oxygen 1s X-ray absorption spectroscopy, we can distinguish between spectral contributions from H 2O and carboxylate, which allows us to probe the electronic structure changes of the atomic site of the carboxylate group being closest to the countercation. From the intensity variations of the COO (-) aq O 1s X-ray absorption peak, which quantitatively correlate with the change in the local partial density of states from the carboxylic site, interactions are found to decrease in the sequence Na (+) > Li (+) > K (+) > NH 4 (+). This ordering, as well as the observed bidental nature of the -COO (-) aq and X (+) aq interaction, is supported by combined ab initio and molecular dynamics calculations.

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.

The Molecular Origin of Like-Charge Arginine−Arginine Pairing in Water

Molecular dynamics simulations show significant like-charge pairing of guanidinium side chains in aqueous poly-arginine, while this effect is absent in aqueous poly-lysine containing ammonium-terminated side chains. This behavior of the guanidinium group is revealed also by protein database searches, having important biochemical implications. Combination of molecular dynamics simulations with explicit solvent and ab initio calculations employing a polarizable continuum model of water allows one to rationalize the formation of contact ion pairs between guanidinium cations in terms of individual interactions at the molecular level.