Search

We employed density functional theory-based ab initio molecular dynamics simulations to examine the hydration structure of several common alkali and alkali earth metal cations. We found that the commonly used atom pairwise dispersion correction scheme D3, which assigns dispersion coefficients based on the neutral form of the atom rather than its actual oxidation state, leads to inaccuracies in the hydration structures of these cations. We evaluated this effect for lithium, sodium, potassium, and calcium and found that the inaccuracies are particularly pronounced for sodium and potassium compared to the experiment. To remedy this issue, we propose disabling the D3 correction specifically for all cation-including pairs, which leads to a much better agreement with experimental data., Main text: 6 pages, 5 figures. Supporting information: 3 pages, 3 figures

We report valence band photoelectron spectroscopy measurements of gas-phase and liquid-phase benzene as well as those of benzene dissolved in liquid ammonia, complemented by electronic structure calculations. The origins of the sizable gas-to-liquid-phase shifts in electron binding energies deduced from the benzene valence band spectral features are quantitatively characterized in terms of the Born-Haber solvation model. This model also allows to rationalize the observation of almost identical shifts in liquid ammonia and benzene despite the fact that the former solvent is polar while the latter is not. For neutral solutes like benzene, it is the electronic polarization response determined by the high frequency dielectric constant of the solvent, which is practically the same in the two liquids, that primarily determines the observed gas-to-liquid shifts.

The change in number densities of aqueous solutions of alkali chlorides should be qualitatively predictable. Typically, as cations get bigger the number density of the solution decreases. However, aqueous solutions of lithium and sodium chloride exhibit at ambient conditions practically identical number densities at equal molalities despite different ionic sizes. Here, we provide an atomistic interpretation of this experimentally observed anomalous behavior using molecular dynamics simulations. The obtained results show that rigidity of the Li+ first and second solvation shells and associated compromised hydrogen bonding result in practically equal average water densities in the local hydration regions for Li+ and Na+ despite different sizes of the cations. In addition, in more distant regions from the cations, the water densities of these two solutions also coincide. These findings thus provide an atomistic interpretation for matching number densities of LiCl and NaCl solutions. In contrast, the number density differences between NaCl and KCl solutions, as well as between LiCl and KCl solutions behave in a regular fashion with lower number densities of solutions observed for larger cations.

Neutron diffraction with isotopic substitution (NDIS) experiments were done on both natural abundance potassium and isotopically labeled 41KCl heavy water solutions to characterize the solvent structuring around the potassium ion in water. Preliminary measurements suggested that the literature value for the coherent neutron scattering length (2.69 fm) for 41K was significantly in error. This value was remeasured using a neutron powder diffractometer and found to be 2.40 fm. This revision increases significantly the contrast between the natural abundance K and 41K by about 30% (from 1.0 to 1.3 fm). The experimentally determined structure factor of the potassium ion was then compared to that calculated from molecular dynamics (MD) simulations. Previous neutron scattering measurements of potassium gave a solvation number of 5.5 (see below). In this study, the NDIS and MD results are in good agreement and allowed us to derive a coordination number of 6.1 for water molecules and 0.8 for chloride ions around each K+ ion in 4 molal aqueous KCl solution.

Nature 595(7869), 673 - 676 (2021). doi:10.1038/s41586-021-03646-5, Insulating materials can in principle be made metallic by applying pressure. In the case of pure water, this is estimated1 to require a pressure of 48 megabar, which is beyond current experimental capabilities and may only exist in the interior of large planets or stars2,3,4. Indeed, recent estimates and experiments indicate that water at pressures accessible in the laboratory will at best be superionic with high protonic conductivity5, but not metallic with conductive electrons1. Here we show that a metallic water solution can be prepared by massive doping with electrons upon reacting water with alkali metals. Although analogous metallic solutions of liquid ammonia with high concentrations of solvated electrons have long been known and characterized6,7,8,9, the explosive interaction between alkali metals and water10,11 has so far only permitted the preparation of aqueous solutions with low, submetallic electron concentrations12,13,14. We found that the explosive behaviour of the water–alkali metal reaction can be suppressed by adsorbing water vapour at a low pressure of about 10$^{−4}$ millibar onto liquid sodium–potassium alloy drops ejected into a vacuum chamber. This set-up leads to the formation of a transient gold-coloured layer of a metallic water solution covering the metal alloy drops. The metallic character of this layer, doped with around 5 × 10$^{21}$ electrons per cubic centimetre, is confirmed using optical reflection and synchrotron X-ray photoelectron spectroscopies., Published by Nature Publ. Group, London [u.a.]

Hydrogen to deuterium isotopic substitution has only a minor effect on physical and chemical properties of water and, as such, is not supposed to influence its neutral taste. Here we conclusively demonstrate that humans are, nevertheless, able to distinguish D2O from H2O by taste. Indeed, highly purified heavy water has a distinctly sweeter taste than same-purity normal water and can add to perceived sweetness of sweeteners. In contrast, mice do not prefer D2O over H2O, indicating that they are not likely to perceive heavy water as sweet. HEK 293T cells transfected with the TAS1R2/TAS1R3 heterodimer and chimeric G-proteins are activated by D2O but not by H2O. Lactisole, which is a known sweetness inhibitor acting via the TAS1R3 monomer of the TAS1R2/TAS1R3, suppresses the sweetness of D2O in human sensory tests, as well as the calcium release elicited by D2O in sweet taste receptor-expressing cells. The present multifaceted experimental study, complemented by homology modelling and molecular dynamics simulations, resolves a long-standing controversy about the taste of heavy water, shows that its sweet taste is mediated by the human TAS1R2/TAS1R3 taste receptor, and opens way to future studies of the detailed mechanism of action., Ben Abu, Mason and colleagues use molecular dynamics, cell-based experiments, mouse models, and human subjects to determine that, unlike ordinary water, heavy water tastes sweet to humans, but not mice. Mechanistically, this effect is mediated by the human TAS1R/TAS1R3 sweet taste receptor.

A versatile, temperature controlled apparatus is presented, which generates deeply cooled liquid microjets of condensed gases, expelling them via a small aperture into vacuum for use in photoelectron spectroscopy (PES). The functionality of the design is demonstrated by temperature- and concentration-dependent PES measurements of liquid ammonia and solutions of KI and NH4I in liquid ammonia. The experimental setup is not limited to the usage of liquid ammonia solutions solely.

Significance The passive translocation mechanism of arginine-rich cell-penetrating peptides has puzzled the scientific community for more than 20 y. In this study we propose a hitherto unrecognized mechanism of passive cell entry involving fusion of multilamellar structures generated by the cell-penetrating peptides. The geometry of entry for this mechanism is completely different from previously suggested direct translocation mechanisms, leading to another paradigm for designing molecular carriers for drug delivery to the cell., Arginine-rich cell-penetrating peptides do not enter cells by directly passing through a lipid membrane; they instead passively enter vesicles and live cells by inducing membrane multilamellarity and fusion. The molecular picture of this penetration mode, which differs qualitatively from the previously proposed direct mechanism, is provided by molecular dynamics simulations. The kinetics of vesicle agglomeration and fusion by an iconic cell-penetrating peptide—nonaarginine—are documented via real-time fluorescence techniques, while the induction of multilamellar phases in vesicles and live cells is demonstrated by a combination of electron and fluorescence microscopies. This concert of experiments and simulations reveals that the identified passive cell penetration mechanism bears analogy to vesicle fusion induced by calcium ions, indicating that the two processes may share a common mechanistic origin.

Insight into the molecular interactions of homotactic and heterotactic association of caffeine and pyridine in aqueous solution is given on the basis of both experimental and simulation studies. Caffeine is about 5 times more soluble in a 3 m aqueous pyridine solution than it is in pure water (an increase from ∼0.1 m to 0.5 m). At this elevated concentration the system becomes suitable for neutron scattering study. Caffeine-pyridine interactions were studied by neutron scattering and molecular dynamics simulations, allowing a detailed characterization of the spatial and orientational structure of the solution. It was found that while pyridine-caffeine interactions are not as strong as caffeine-caffeine interactions, the pyridine-caffeine interactions still significantly disrupted caffeine-caffeine stacking. The alteration of the caffeine-caffeine stacking, occasioned by the presence of pyridine molecules in solution and the consequent formation of heterotactic interactions, leads to the experimentally detected increase in caffeine solubility.

The molecular structure and strength of a model salt bridge between a guanidinium cation as the charged side chain group of arginine and the carboxylic group of acetate in an aqueous solutions is characterized by a combination of neutron diffraction with isotopic substitution and molecular dynamics simulations. Being able to recover the second order difference signal, the present neutron scattering experiments provide direct information about ion pairing in the investigated solution. At the same time, these measurements serve as benchmarks for assessing the quality of the force field employed in the simulation. We show that a standard non-polarizable force field, which tends to overestimate the strength of salt bridges, does not reproduce the structural features from neutron scattering pertinent to ion pairing. In contrast, a quantitative agreement with experiment is obtained when electronic polarization effects are accounted for in a mean-field way via charge scaling. Such simulations are then used to quantify the weak character of a fully hydrated salt bridge. Finally, on top of the canonical hydrogen-bonding binding mode between guanidinium and acetate, these simulations also point to another interaction motif involving an out-of-plane hydrophobic contact of the methyl group of acetate with the guanidinium cation.

Photoelectron spectroscopy of microjets expanded into vacuum allows access to orbital energies for solute or solvent molecules in the liquid phase. Microjets of water, acetonitrile and alcohols have previously been studied; however, it has been unclear whether jets of low temperature molecular solvents could be realized. Here we demonstrate a stable 20 μm jet of liquid ammonia (-60 °C) in a vacuum, which we use to record both valence and core-level band photoelectron spectra using soft X-ray synchrotron radiation. Significant shifts from isolated ammonia in the gas-phase are observed, as is the liquid-phase photoelectron angular anisotropy. Comparisons with spectra of ammonia in clusters and the solid phase, as well as spectra for water in various phases potentially reveal how hydrogen bonding is reflected in the condensed phase electronic structure.

We present a combination of force field and ab initio molecular dynamics simulations together with neutron scattering experiments with isotopic substitution that aim at characterizing ion hydration and pairing in aqueous calcium chloride and formate/acetate solutions. Benchmarking against neutron scattering data on concentrated solutions together with ion pairing free energy profiles from ab initio molecular dynamics allows us to develop an accurate calcium force field which accounts in a mean-field way for electronic polarization effects via charge rescaling. This refined calcium parameterization is directly usable for standard molecular dynamics simulations of processes involving this key biological signaling ion.

It is a textbook knowledge that charges of the same polarity repel each other. For two monovalent ions in the gas phase at a close contact this repulsive interaction amounts to hundreds of kilojoules per mole. In aqueous solutions, however, this Coulomb repulsion is strongly attenuated by a factor equal to the dielectric constant of the medium. The residual repulsion, which now amounts only to units of kilojoules per mole, may be in principle offset by attractive interactions. Probably the smallest cationic pair, where a combination of dispersion and cavitation forces overwhelms the Coulomb repulsion, consists of two guanidinium ions in water. Indeed, by a combination of molecular dynamics with electronic structure calculations and electrophoretic, as well as spectroscopic, experiments, we have demonstrated that aqueous guanidinium cations form (weakly) thermodynamically stable like- charge ion pairs. The importance of pairing of guanidinium cations in aqueous solutions goes beyond a mere physical curiosity, since it has significant biochemical implications. Guanidinium chloride is known to be an efficient and flexible protein denaturant. This is due to the ability of the orientationally amphiphilic guanidinium cations to disrupt various secondary structural motifs of proteins by pairing promiscuously with both hydrophobic and hydrophilic groups, including guanidinium- containing side chains of arginines. The fact that the cationic guanidinium moiety forms the dominant part of the arginine side chain implies that the like-charge ion pairing may also play a role for interactions between peptides and proteins. Indeed, arginine– arginine pairing has been frequently found in structural protein databases. In particular, when strengthened by a presence of negatively charged glutamate, aspartate, or C-terminal carboxylic groups, this binding motif helps to stabilize peptide or protein dimers and is also found in or near active sites of several enzymes. The like-charge pairing of the guanidinium side-chain groups may also hold the key to the understanding of the arginine “magic”, that is, the extraordinary ability of arginine-rich polypeptides to passively penetrate across cellular membranes. Unlike polylysines, which are also highly cationic but lack the ease in crossing membranes, polyarginines do not exhibit mutual repulsion. Instead, they accumulate at the membrane, weaken it, and might eventually cross in a concerted, “train-like” manner. This behavior of arginine-rich cell penetrating peptides can be exploited when devising smart strategies how to deliver in a targeted way molecular cargos into the cell.

Modeled ions, described by nonpolarizable force fields, can suffer from unphysical ion pairing and clustering in aqueous solutions well below their solubility limit. The electronic continuum correction takes electronic polarization effects of the solvent into account in an effective way by scaling the charges on the ions, resulting in a much better description of the ionic behavior. Here, we present parameters for the sodium ion consistent with this effective polarizability approach and in agreement with experimental data from neutron scattering, which could be used for simulations of complex aqueous systems where polarization effects are important.

Neutron diffraction with isotopic substitution was performed on aqueous solutions of isopropyl alcohol and isopropylamine. The difference between these two measurements primarily contains information about the different hydration of the alcohol and amino group. This data is used as a test of the accuracy of molecular dynamic simulations of the same systems. Having established the level of accuracy of the modeling, it is employed as an interpretive tool for the experimental data. Even though the alcohol and the amine possess comparable hydrogen bonding capabilities, consisting respectively of either two hydrogen bond acceptors and one donor, or two hydrogen bond donors and one acceptor, we find significant differences in the hydration of the hydroxyl and amino groups.

Alkali metals can react explosively with water and it is textbook knowledge that this vigorous behaviour results from heat release, steam formation and ignition of the hydrogen gas that is produced. Here we suggest that the initial process enabling the alkali metal explosion in water is, however, of a completely different nature. High-speed camera imaging of liquid drops of a sodium/potassium alloy in water reveals submillisecond formation of metal spikes that protrude from the surface of the drop. Molecular dynamics simulations demonstrate that on immersion in water there is an almost immediate release of electrons from the metal surface. The system thus quickly reaches the Rayleigh instability limit, which leads to a ‘coulomb explosion’ of the alkali metal drop. Consequently, a new metal surface in contact with water is formed, which explains why the reaction does not become self-quenched by its products, but can rather lead to explosive behaviour. The explosion of alkali metals in water is a typical high-school chemistry experiment, but its mechanism is not fully understood. Using high-speed cameras and molecular simulations it is now shown that a key early step in this reaction is the migration of electrons from the alkali metal into water, leading to a charging of the metal's surface and subsequent Coulomb explosion.

New results derived from the experimental method of neutron diffraction and isotopic substitution (NDIS) are presented for the hydration structure of the lithium cation (Li(+)) in aqueous solutions of lithium chloride in heavy water (D2O) at concentrations of 6, 3, and 1 m and at 1.5 m lithium sulfate. By introducing new and more-accurate data reduction procedures than in our earlier studies (I. Howell and G. W. Neilson, J. Phys: Condens. Matter, 1996, 8, 4455-4463), we find, in the first hydration shell of Li(+), ∼4.3(2) water molecules at 6 m, 4.9(3) at 3 m, 4.8(3) at 1 m in the LiCl solutions, and 5.0(3) water molecules in the case of Li2SO4 solution. The general form of the first hydration shell is similar in all four solutions, with the correlations for Li-O and Li-D sited at 1.96 (0.02) Å and 2.58 (0.02) Å, respectively. The results resemble those presented in 1996, in terms of ion-water distances and local coordination, but the hydration number is significantly lower for the case at 1 m than the 6.5 (1.0) given at that time. Thus, experimental and theoretical results now agree that lithium is hydrated by a small number of water molecules (4-5) in the nearest coordination shell.

Recent advances in computational and experimental techniques have allowed for accurate description of ion pairing in aqueous solutions. Free energy methods based on ab initio molecular dynamics, as well as on force fields accounting effectively for electronic polarization, can provide quantitative information about the structures and occurrences of individual types of ion pairs. When properly benchmarked against electronic structure calculations for model systems and against structural experiments, in particular neutron scattering, such force field simulations represent a powerful tool for elucidating interactions of salt ions in complex biological aqueous environments.

Magnesium and zinc dications possess the same charge and have an almost identical size, yet they behave very differently in aqueous solutions and play distinct biological roles. It is thus crucial to identify the origins of such different behaviors and to assess to what extent they can be captured by force-field molecular dynamics simulations. In this work, we combine neutron scattering experiments in a specific mixture of H

Magnesium and zinc dications possess the same charge and have an almost identical size, yet they behave very differently in aqueous solutions and play distinct biological roles. It is thus crucial to identify the origins of such different behaviors and to assess to what extent they can be captured by force-field molecular dynamics simulations. In this work, we combine neutron scattering experiments in a specific mixture of H2O and D2O (the so-called null water) with ab initio molecular dynamics simulations to probe the difference in the hydration structure and ion-pairing properties of chloride solutions of the two cations. The obtained data are used as a benchmark to develop a scaled-charge force field for Mg2+ that includes electronic polarization in a mean field way. We show that using this electronic continuum correction we can describe aqueous magnesium chloride solutions well. However, in aqueous zinc chloride specific interaction terms between the ions need to be introduced to capture ion pairing q...

International audience; The imidazole motif is widely encountered in biomolecules, and its biological role, for instance as a proton relay, is often linked to its ability to form hydrogen bonds with water molecules. Detailed characterization of the hydration pattern of imidazole and its changes upon protonation is thus of high interest. Here, we combine neutron scattering experiments with force field simulations to provide an unprecedented characterization of the neutral and protonated imidazole solvation at the atomistic level. We show that neutron diffraction data can be used to assess the quality of the imidazole force field in molecular simulations. Simulations using the CHARMM general force field for imidazole are in excellent agreement with the experimental neutron scattering data and we use them to provide an atomic scale interpretation of the neutron scattering patterns. Upon protonation, we clearly identify the signature of the reorganization in the hydration pattern caused by the change from one H-bond donor and one H-bond acceptor group for imidazole to two H-bond donor groups for imidazolium. We also point the limits of the experiment, which is rather insensitive to details of the H-bond geometry at the deprotonated nitrogen of imidazole and further complement the description of the hydration structure with ab initio molecular dynamics simulations.

Neutron scattering experiments were performed on 6 m LiCl solutions in order to obtain the solvation structure around the chloride ion. Molecular dynamics simulations on systems mirroring the concentrated electrolyte conditions of the experiment were carried out with a variety of chloride force-fields. In each case the simulations were run with both full ionic charges and employing the electronic continuum correction (implemented through charge scaling) to account effectively for electronic polarisation. The experimental data were then used to assess the successes and shortcomings of the investigated force-fields. We found that due to the very good signal-to-noise ratio in the experimental data, they provide a very narrow window for the position of the first hydration shell of the chloride ion. This allowed us to establish the importance of effectively accounting for electronic polarisation, as well as adjusting the ionic size, for obtaining a force-field which compares quantitatively to the experimental ...

The structure of aqueous solutions of methyl β-D-ribofuranoside was investigated by coupling molecular dynamics (MD) simulations and neutron scattering measurements with isotopic substitution. Using a sample of the sugar isotopically-labeled at a single unique position, neutron scattering structure factors and radial distribution functions can be compared with MD simulations constrained to different conformations to determine which conformer best fits the experimental results. Three different simulations were performed with the methyl ether group of the sugar unconstrained and constrained in each of its staggered orientations. The results of the unconstrained simulation showed that the methyl ester group occupied predominantly the 300° position, which is in agreement with the diffraction experimental results. This result suggests that the molecular mechanics force field used in the simulation adequately describes the conformation of the 1-methyl ether group in the methyl β-D-ribofuranoside.

Molecular dynamics simulations of concentrated aqueous solutions of LiCl and Li2SO4 were conducted in order to provide molecular insight into recent neutron scattering data. The structures predicted from the molecular dynamics simulations using standard nonpolarizable force fields provided a very poor fit to the experiment; therefore, refinement was needed. The electronic polarizability of the medium was effectively accounted for by implementing the electronic continuum correction, which practically means rescaling the ionic charges. Consistent with previous studies, we found that this approach in each case provided a significantly improved fit to the experimental data, which was further enhanced by slightly adjusting the radius of the lithium ion. The polarization effect was particularly pronounced in the Li2SO4 solution where the ions in the nonpolarizable simulations tended to cluster unphysically. With the above alterations, the employed force field displayed an excellent fit to the neutron scattering data and provided a useful interpretative framework for the experimental measurements. At the same time, the present study underlines the importance of solvent polarization effects in hydration of ions with high charge density.

Alkali metals in water are always at the brink of explosion. Herein, we show that this vigorous reaction can be kept in a non-exploding regime, revealing a fascinating richness of hitherto unexplored chemical processes. A combination of high-speed camera imaging and visible/near-infrared/infrared spectroscopy allowed us to catch and characterize the system at each stage of the reaction. After gently placing a drop of a sodium/potassium alloy on water under an inert atmosphere, the production of solvated electrons became so strong that their characteristic blue color could be observed with the naked eye. The exoergic reaction leading to the formation of hydrogen and hydroxide eventually heated the alkali metal drop such that it became glowing red, and part of the metal evaporated. As a result of the reaction, a perfectly transparent drop consisting of molten hydroxide was temporarily stabilized on water through the Leidenfrost effect, bursting spectacularly after it had cooled sufficiently.

Gating charges in voltage-sensing domains (VSD) of voltage-sensitive ion channels and enzymes are carried on arginine side chains rather than lysine. This arginine preference may result from the unique hydration properties of the side chain guanidinium group which facilitates its movement through a hydrophobic plug that seals the center of the VSD, as suggested by molecular dynamics simulations. To test for side chain interactions implicit in this model we inspected interactions of the side chains of arginine and lysine with each of the 19 non-glycine amino acids in proteins in the protein data bank. The arginine guanidinium interacts with non-polar aromatic and aliphatic side chains above and below the guanidinium plane while hydrogen bonding with polar side chains is restricted to in-plane positions. In contrast, non-polar side chains interact largely with the aliphatic part of the lysine side chain. The hydration properties of arginine and lysine are strongly reflected in their respective interactions with non-polar and polar side chains as observed in protein structures and in molecular dynamics simulations and likely underlie the preference for arginine as a mobile charge carrier in VSD.

The behavior of guanidinium chloride at the surface of aqueous solutions is investigated using classical molecular dynamics simulation. It is found that the population of guanidinium ions oriented parallel to the interface is greater in the surface region than in bulk. The opposite is true for ions in other orientations. Overall, guanidinium chloride is depleted in the surface region, in agreement with the fact that addition of guanidinium chloride increases the surface tension of water. The orientational dependence of the surface affinity of the guanidinium cation appears to be related to its anisotropic hydration. To bring the ion to the surface in the parallel orientation does not require hydrogen bonds to be broken, in contrast to other orientations. The surface enrichment of parallelly oriented guanidinium indicates that its solvation is more favorable near the surface than in bulk solution for this orientation. The dependence of the bulk and surface properties of guanidinium on the force field parameters, including polarizability, is also investigated. Despite significant quantitative differences between the force fields, the surface behavior is qualitatively robust. The implications of the above findings for the behaviorof guanidinium groups near biological interfaces are demonstrated on the example of arginine side chain orientations in lysozyme.

Molecular dynamics simulations were carried out on a system of eight independent caffeine molecules in a periodic box of water at 300 K, representing a solution near the solubility limit for caffeine at room temperature, using a newly-developed CHARMM-type force field for caffeine in water. Simulations were also conducted for single caffeine molecules in water using two different water models (TIP3P and TIP4P). Water was found to structure in a complex fashion around the planar caffeine molecules, which was not sensitive to the water model used. As expected, extensive aggregation of the caffeine molecules was observed, with the molecules stacking their flat faces against one another like coins, with their methylene groups staggered to avoid steric clashes. A dynamic equilibrum was observed between large n-mers, including stacks with all eight solute molecules, and smaller clusters, with the calculated osmotic coefficient being in acceptable agreement with the experimental value. The insensitivity of the results to water model and the congruence with experimental thermodynamic data suggest that the observed stacking interactions are a realistic representation of the actual association mechanism in aqueous caffeine solutions.

Surface affinity of aqueous guanidinium chloride (GdmCl) is compared to that of aqueous tetrapropylammonium chloride (TPACl) upon addition of sodium chloride (NaCl) or disodium sulfate (Na2SO4). The experimental results have been acquired using the surface sensitive technique X-ray photoelectron spectroscopy on a liquid jet. Molecular dynamics (MD) simulations have been used to produce radial distribution functions and surface density plots. The surface affinities of both TPA+ and Gdm+ increase upon adding NaCl to the solution. With the addition of Na2SO4 the surface affinity of TPA+ increases while that of Gdm+ decreases. From the results of MD simulations it is seen that Gdm+ and SO4^2- ion pairs. This finding can be used to explain the decreased surface affinity of Gdm+ when co-dissolved with SO4^2- ions. Since SO4^2- ions avoid the surface due to the double charge and strong water interaction, the Gdm+-SO4^2- ion pair resides deeper in the solutions' bulk than the Gdm+ ions. Since TPA+ does not form ion pairs with SO4^2- the TPA+ ions are instead enriched at the surface.

Neutron diffraction with isotopic substitution (NDIS) experiments have been used to examine the structuring of aqueous solutions of pyridine. A new method is described for extracting the structure factors relating to intermolecular correlations from neutron scattering experiments on liquid solutions of complex molecular species. This approach performs experiments at different concentrations and exploits the intramolecular coordination number concentration invariance (ICNCI) to separate the internal and intermolecular contributions to the total intensities. The ability of this method to deconvolute molecular and intermolecular correlations is tested and demonstrated using simulated neutron scattering results predicted from molecular dynamics simulations of aqueous solutions of the polyatomic solute pyridine in which the inter- and intramolecular terms are known. The method is then implemented using neutron scattering measurements on solutions of pyridine. The results confirm that pyridine shows a significant propensity to aggregate in solution and demonstrate the prospects for the future application of the ICNCI approach to the study of large polyatomic solutes using next-generation neutron sources and detectors.

Small angle neutron scattering (SANS) and molecular dynamics (MD) simulations were used to characterize the long range structuring (aggregation) of aqueous solutions of isopropanol (IPA) and pyridine, and the effect on structuring of guanidinium chloride (GdmCl). These solutes serve as highly soluble analogs of the non-polar aliphatic (IPA), and aromatic (pyridine) side chains of proteins. SANS data showed that isopropanol and pyridine both form clusters in water resulting from interaction between non-polar groups of the solutes, with pyridine aggregation producing longer-range structuring than isopropanol in 3 molal solutions. Addition of GdmCl at 3 molal concentration considerably reduced pyridine aggregation, but had no effect on isopropanol aggregation. MD simulations of these solutions support the conclusion that long range structuring involves hydrophobic solute interactions, and that Gdm+ interacts with the planar pyridine group to suppress pyridine-pyridine interactions in solution. Hydrophobic interactions involving the aliphatic groups of isopropanol were unaffected by GdmCl, indicating that the planar and weakly hydrated Gdm+ cation cannot make productive interactions with the highly curved or “lumpy” aliphatic groups of this solute. These observations support the conclusion that the effects of Gdm+ ions on protein-stabilizing interactions involving aromatic amino acid side chains make significant contributions to the denaturant activity of GdmCl, whereas interactions with the “lumpy” aliphatic side chains are likely to be less important.

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.

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.

Guanidinium (Gdm+) chloride is a powerful protein denaturant, whereas the sulfate dianion (SO42-) is a strong stabilizer of folded protein states; Gdm2SO4 is effectively neutral in its effects on protein stability. While the "neutralizing" effects of protein-stabilizing solutes on the activity of denaturants can be broadly interpreted in terms of additive effects of the solutes, recent experimental and simulation studies support a role for hetero-ion interactions in the effect of sulfate on Gdm+ denaturation [Mason, P. E.; et al. J. Phys. Chem. B 2005, 109, 24185-24196]. Here we describe an experimental strategy for testing this mechanism that involves spectroscopic analysis of the separate effects of alkali metal sulfates (Na2SO4, Rb2SO4), GdmCl, and Gdm2SO4 on the folded populations of several peptides chosen to dissect specific noncovalent contributions to the conformational stability of proteins [alanine-based helical peptides stabilized by hydrogen bonds, tryptophan zipper (trpzip) peptides stabilized largely by cross-strand indole-indole interactions]. While the trpzip peptides are highly sensitive to GdmCl denaturation, they are unaffected by NaCl, Na2SO4, or Gdm2SO4, indicating that the reversal of the denaturant activity of Gdm+ by sulfate in this case is not due to competing stabilizing (sulfate) and destabilizing (Gdm+) interactions. Gdm2SO4 was found to retain considerable denaturant activity against alanine-based alpha-helical peptides. The differences in the effects of Gdm2SO4 on the two peptide types can be understood in terms of the different mechanisms of Gdm+ denaturation of trpzip peptides and helical peptides, respectively, and the specific nature of Gdm+ and SO42- ionic "clustering" that differentially affects the ability of Gdm+ to make the molecular interactions with the peptides that underlie its denaturant activity.

Molecular dynamics simulations of pure liquid water under ambient conditions using four common empirical water models have been analyzed to determine how well the oxygen-oxygen radial distribution function, g(OO)(r), used as the sole criterion of congruence with experiment, captures variations in the actual anisotropic collective structuring for these models. The largest systematic deviations from tetrahedrality were found to be due to deformations of the angle between the two closet hydrogen bond donor neighbors, but for intrinsic geometric reasons, these were found to contribute less to g(OO)(r) than deformations of the angles between one hydrogen bond donor neighbor and one hydrogen bond acceptor neighbor. Relying exclusively on a qualitative characterization of the second peak in g(OO)(r) seems to overemphasize the differences between the structuring in some of these models.

The presence of ions in water provides a rich and varied environment in which many natural processes occur with important consequences in biology, geology and chemistry. This article will focus on the structural properties of ions in water and it will be shown how the 'difference' methods of neutron diffraction with isotopic substitution (NDIS) and anomalous X-ray diffraction (AXD) can be used to obtain direct information regarding the radial pair distribution functions of many cations and anions in solution. This information can subsequently be used to calculate coordination numbers and to determine ion-water conformation in great detail. As well as enabling comparisons to be made amongst ions in particular groups in the periodic table, such information can also be contrasted with results provided by molecular dynamics (MD) simulation techniques. To illustrate the power of these 'difference' methods, reference will be made to the alkali group of ions, all of which have been successfully investigated by the above methods, with the exception of the radioactive element francium. Additional comments will be made on how NDIS measurements are currently being combined with MD simulations to determine the structure around complex ions and molecules, many of which are common in biological systems.

Neutron diffraction with isotopic substitution (NDIS) experiments and molecular dynamics (MD) simulations have been used to study the structuring in aqueous solution of two cesium salts, cesium carbonate, and cesium nitrate. As was previously found for guanidinium salts of carbonate, mesoscopic-scale clusters were seen to form in the Cs2CO3 solution both in the MD simulations and in the diffraction experiments. No such large scale ion clusters were found in the CsNO3 solutions in either the modeling or experiments. The results are dominated by the strength and geometry of the direct first-neighbor interactions, which explain the differences in the clustering behavior between the two solutions without need to refer to longer-range water-water structuring.

A method of neutron diffraction is described which enables the first hydration zone of small cations to be investigated at atomic resolution. It is shown that the cation structures of aqueous electrolyte solutions dissolved in a 'null' mixture of water (H(2)O) and heavy water (D(2)O), can be calculated directly from the neutron scattering patterns. The hitherto unresolved structure around Na(+) is used to illustrate the power of this method, the accuracy of which is discussed formally with reference to standard nickel chloride solutions. Possible applications to a variety of other systems and at different thermodynamic states are proposed.

Neutron diffraction with isotopic substitution (NDIS) experiments and molecular dynamics (MD) simulations have been used to examine the pentose D-xylose in aqueous solution. By specifically labeling D-xylose molecules with a deuterium atom at the nonexchangeable hydrogen position on C4, it was possible to extract information about the atomic structuring around just that specific position. The MD simulations were found to give satisfactory agreement with the experimental NDIS results and could be used to help interpret the scattering data in terms of the solvent structuring as well as the intramolecular hydroxyl conformations. Although the experiment is challenging and on the limit of modern instrumentation, it is possible by careful analysis, in conjunction with MD studies, to show that the conformation trans to H4 at 180 degrees is strongly disfavored, in excellent agreement with the MD results. This is the first attempt to use NDIS experiments to determine the rotameric conformation of a hydroxyl group.

Atomic ordering in network glasses on length scales longer than nearest-neighbour length scales has long been a source of controversy(1-6). Detailed experimental information is therefore necessary to understand both the network properties and the fundamentals of glass formation. Here we address the problem by investigating topological and chemical ordering in structurally disordered AX2 systems by applying the method of isotopic substitution in neutron diffraction to glassy ZnCl2. This system may be regarded as a prototypical ionic network forming glass, provided that ion polarization effects are taken into account(7), and has thus been the focus of much attention(8-14). By experiment, we show that both the topological and chemical ordering are described by two length scales at distances greater than nearest-neighbour length scales. One of these is associated with the intermediate range, as manifested by the appearance in the measured diffraction patterns of a first sharp diffraction peak at 1.09( 3) angstrom(-1); the other is associated with an extended range, which shows ordering in the glass out to 62( 4) angstrom. We also find that these general features are characteristic of glassy GeSe2, a prototypical covalently bonded network material(15,16). The results therefore offer structural insight into those length scales that determine many important aspects of supercooled liquid and glass phenomenology(11).

Understanding the molecular basis for protein denaturation by urea and guanidinium chloride (GdmCl) should accommodate the observation that, on a molar basis, GdmCl is generally 2-2.5-fold more effective as a protein denaturant than urea. Previous studies [Smith, J. S., and Scholtz, J. M. (1996) Biochemistry 35, 7292-7297] have suggested that the effects of GdmCl on the stability of alanine-based helical peptides can be separated into denaturant and salt effects, since adding equimolar NaCl to urea enhanced urea-induced unfolding to an extent that was close to that of Gdm. We reinvestigated this observation using an alanine-based helical peptide (alahel) that lacks side chain electrostatic contributions to stability, and compared the relative denaturant sensitivities of this peptide with that of tryptophan zipper peptides (trpzip) whose native conformations are stabilized largely by cross-strand indole ring interactions. In contrast to the observations of Smith and Scholtz, GdmCl was only slightly more powerful as a denaturant of alahel than urea in salt-free buffer (the denaturant m value m(GdmCl)/m(urea) ratio = 1.4), and the denaturation of alahel by urea exhibited only a small dependence on NaCl or KCl. The trpzip peptides were much more sensitive to GdmCl than to urea (m(GdmCl)/m(urea) = 3.5-4). These observations indicate that the m(GdmCl)/m(urea) ratio of 2-2.5 for proteins results from a combination of effects on the multiple contributions to protein stability, for which GdmCl may be only slightly more effective than urea (e.g., hydrogen bonds) or considerably more effective than urea (e.g., indole-indole interactions).

Neutron diffraction experiments were carried out at ambient temperature on aqueous solutions of glucose at concentrations of 1, 3, and 5 molal. The difference methods of neutron diffraction and isotopic substitution were used to determine various structure factors associated with all the exchangeable hydrogen atoms of the solutions. The results imply that the hydration of the glucose molecule is not concentration sensitive. It is also found that the average water structure itself is relatively unperturbed by the presence of glucose.

Neutron diffraction experiments were carried out on aqueous solutions containing either guanidinium or thiocyanate ions. The first-order difference method of neutron diffraction and isotopic substitution was applied, and the hydration structures of two of nature's strongest denaturant ions were determined. Each ion is shown to interact weakly with water: Guanidinium has no recognizable hydration shell and is one of the most weakly hydrated cations yet characterized. Hydration of thiocyanate is characterized by a low coordination number involving around one hydrogen-bonded water molecule and approximately two water molecules weakly interacting through “hydration bonds.” The weak hydration of these denaturant ions strongly supports suggestions that a major contribution to the denaturant effect is the preferential interaction of the denaturant with the protein surface. By contrast, solute species such as many sugars and related polyols that stabilize proteins are strongly hydrated and are thus preferentially retained in the bulk solvent and excluded from the protein surface.

Calcium is one of the biologically most important ions; however, its accurate description by classical molecular dynamics simulations is complicated by strong electrostatic and polarization interactions with surroundings due to its divalent nature. Here, we explore the recently suggested approach for effectively accounting for polarization effects via ionic charge rescaling and develop a new and accurate parametrization of the calcium dication. Comparison to neutron scattering and viscosity measurements demonstrates that our model allows for an accurate description of concentrated aqueous calcium chloride solutions. The present model should find broad use in efficient and accurate modeling of calcium in aqueous environments, such as those encountered in biological and technological applications.

The presence of ions and/or apolar species in water provides a rich and varied environment in which many natural processes occur. This review provides results of recent structural studies of aqueou...

[formula: see text] Inspired by folded, nonpseudorotaxane complexes of bis(m-phenylene)-32-crown-10 systems, we synthesized a new bicyclic crown ether containing two 1,3,5-phenylene units linked by three tetra(ethyleneoxy) units. The new cryptand forms a "pseudorotaxane-like" inclusion complex with N,N'-dimethyl-4,4'-bipyridinium bis(hexafluorophosphate) with association constant Ka = 6.1 x 10(4) M-1, 100-fold greater than that of an analogous simple crown ether.

First, the complexation of rigid bis(paraquat)-based cyclophane host 1 with linear MDI-based polyurethane guests 2−5 prepared from tetra(ethylene glycol) was studied. Through detailed analysis of 1H NMR data, we have shown a relationship between the length of the linear polyurethanes and the fraction of all threaded cyclophanes rapidly exchanging with solution upon the 1H NMR time scale. In other words, the molecular weight of the polymer can be estimated from the proportion of cyclic threaded but exchanging rapidly, i.e., those near the ends, versus those that are threaded but exchanging slowly because they do not have ready access to the ends of the macromolecule; this is essentially an end group analysis. And then in a similar study of branched polyurethanes 7−11 derived from tetra(ethylene glycol) and glycerol, it was shown that the cyclophane can be used to determine the (accessible) fraction of the polymer between a branch point and a terminus and the average length between a branch point and a term...