Effect of molecular shape and electrostatic interactions on the water layer around polar and apolar groups in solution

Charge fluctuations in charge-regulated systems: dependence on statistical ensemble

We investigate charge regulation of nanoparticles in concentrated suspensions, focusing on the effect of different statistical ensembles. We find that the choice of ensemble does not affect the mean charge of nanoparticles, but significantly alters the magnitude of its fluctuation. Specifically, we compared the behaviors of colloidal charge fluctuations in the semi-grand canonical and canonical ensembles and identified significant differences between the two. The choice of ensemble-whether the system is isolated or is in contact with a reservoir of acid and salt-will, therefore, affect the Kirkwood-Shumaker fluctuation-induced force inside concentrated suspensions. Our results emphasize the importance of selecting an appropriate ensemble that accurately reflects the experimental conditions when studying fluctuation-induced forces between polyelectrolytes, proteins, and colloidal particles in concentrated suspensions.

Hydrophobic Cluster Formation in Aqueous Ethanol Solutions Probed by Soft X-ray Absorption Spectroscopy

Hydrophobic cluster structures in aqueous ethanol solutions at different concentrations have been investigated by soft X-ray absorption spectroscopy (XAS). In the O K-edge XAS, we have found that hydrogen bond structures among water molecules are enhanced in the middle-concentration region by the hydrophobic interaction of the ethyl groups in ethanol. In the C K-edge XAS, the lower energy features arise from a transition from the terminal methyl C 1s electron to an unoccupied orbital of 3s Rydberg character, which is sensitive to the nearest-neighbor intermolecular interactions. From the comparison of C K-edge XAS with the inner-shell calculations, we have found that ethanol clusters are easily formed in the middle-concentration region due to the hydrophobic interaction of the ethyl group in ethanol, resulting in the enhancement of the hydrogen bond structures among water molecules. This behavior is different from aqueous methanol solutions, where the methanol-water mixed clusters are more predominant in the middle-concentration region due to the relatively weak hydrophobic interactions of the methyl group in methanol.

Molecular photoswitches in aqueous environments

Molecular photoswitches enable dynamic control of processes with high spatiotemporal precision, using light as external stimulus, and hence are ideal tools for different research areas spanning from chemical biology to smart materials. Photoswitches are typically organic molecules that feature extended aromatic systems to make them responsive to (visible) light. However, this renders them inherently lipophilic, while water-solubility is of crucial importance to apply photoswitchable organic molecules in biological systems, like in the rapidly emerging field of photopharmacology. Several strategies for solubilizing organic molecules in water are known, but there are not yet clear rules for applying them to photoswitchable molecules. Importantly, rendering photoswitches water-soluble has a serious impact on both their photophysical and biological properties, which must be taken into consideration when designing new systems. Altogether, these aspects pose considerable challenges for successfully applying molecular photoswitches in aqueous systems, and in particular in biologically relevant media. In this review, we focus on fully water-soluble photoswitches, such as those used in biological environments, in both in vitro and in vivo studies. We discuss the design principles and prospects for water-soluble photoswitches to inspire and enable their future applications.

Highly Selective Supported Graphene Oxide Membranes for Water-Ethanol Separation

A polyethersulfone (PES)-supported graphene oxide (GO) membrane has been developed by a simple casting approach. This stable membrane is applied for ethanol/water separation at different temperatures. The 5.0 µm thick GO film coated on PES support membrane showed a long-term stability over a testing period of one month and excellent water/ethanol selectivity at elevated temperatures. The water/ethanol selectivity is dependent on ethanol weight percentage in water/ethanol feed mixtures and on operating temperature. The water/ethanol selectivity was enhanced with an increase of ethanol weight percentage in water/ethanol mixtures, from below 100 at RT to close to 874 at a 90 °C for 90% ethanol/10% water mixture. Molecular dynamics simulation of water-ethanol mixtures in graphene bilayers, that are considered to play a key role in transport, revealed that molecular transport is negligible for layer spacing below 1 nm. The differences in the diffusion of ethanol and water in the bilayer are not consistent with the large selectivity value experimentally observed. The entry of water and ethanol into the interlayer space may be the crucial step controlling the selectivity.

Induction of Au-methotrexate conjugates by sugar molecules: production, assembly mechanism, and bioassay studies

Au-methotrexate (Au-MTX) conjugates induced by sugar molecules were produced by a simple, one-pot, hydrothermal growth method. Herein, the Au(III)-MTX complexes were used as the precursors to form Au-MTX conjugates. Addition of different types of sugar molecules with abundant hydroxyl groups resulted in the formation of Au-MTX conjugates featuring distinct characteristics that could be explained by the diverse capping mechanisms of sugar molecules. That is, the instant-capping mechanism of glucose favored the generation of peanut-like Au-MTX conjugates with high colloidal stability while the post-capping mechanism of dextran and sucrose resulted in the production of Au-MTX conjugates featuring excellent near-infrared (NIR) optical properties with a long-wavelength plasmon resonance near 630-760 nm. Moreover, in vitro bioassays showed that cancer cell viabilities upon incubation with free MTX, Au-MTX conjugates doped with glucose, dextran and sucrose for 48 h were 74.6%, 55.0%, 62.0%, and 63.1%, respectively. Glucose-doped Au-MTX conjugates exhibited a higher anticancer activity than those doped with dextran and sucrose, therefore potentially presenting a promising treatment platform for anticancer therapy. Based on the present study, this work may provide the first example of using biocompatible sugars as regulating agents to effectively guide the shape and assembly behavior of Au-MTX conjugates. Potentially, the synergistic strategy of drug molecules and sugar molecules may offer the possibility to create more gold-based nanocarriers with new shapes and beneficial features for advanced anticancer therapy.

Compressed-tube pressure cell for optical studies at ocean pressures: Application to glucose mutarotation kinetics

A self-contained compressed-tube pressure cell is tested to 25 MPa. The cell is very simple to construct and offers stable pressure control with optical access to fluid samples. The physical path length of light through the cell is large enough to measure optical activity. The entire system is relatively small and portable, and it is vibration-free, since a compressor is not used. Operation of the cell is demonstrated by measuring the mutarotation rate of aqueous glucose solutions at 25 °C. A logarithmic plot of the rate constant vs. pressure yields an activation volume for mutarotation of -22 cm³/mol, approximately twice the value measured previously at higher pressures.

Physical origin underlying the entropy loss upon hydrophobic hydration

The hydrophobic effect (HE) is commonly associated with the demixing of oil and water at ambient conditions and plays the leading role in determining the structure and stability of biomolecular assembly in aqueous solutions. On the molecular scale HE has an entropic origin. It is believed that hydrophobic particles induce order in the surrounding water by reducing the volume of configuration space available for hydrogen bonding. Here we show with computer simulation results that this traditional picture, based on average structural features of hydration water, configurational properties of single water molecules, and up to pairwise correlations, is not correct. Analyzing collective fluctuations in water clusters we are able to provide a fundamentally new picture of HE based on pronounced many-body correlations affecting the switching of hydrogen bonds (HBs) between molecules. These correlations emerge as a nonlocal compensation of reduced fluctuations of local electrostatic fields in the presence of an apolar solute. We propose an alternative view which may also be formulated as a maximization principle: The electrostatic noise acting on water molecules is maximized under the constraint that each water molecule on average maintains as many HBs as possible. In the presence of the solute the maximized electrostatic noise is a result of nonlocal fluctuations in the labile HB network giving rise to strong correlations among at least up to four water molecules.

Aqueous solvation of methane from first principles

Structural, dynamical, bonding, and electronic properties of water molecules around a soluted methane molecule are studied from first principles. The results are compatible with experiments and qualitatively support the conclusions of recent classical molecular dynamics simulations concerning the controversial issue on the presence of "immobilized" water molecules around hydrophobic groups: the hydrophobic solute slightly reduces (by a less than 2 factor) the mobility of many surrounding water molecules rather than immobilizing just the few ones which are closest to methane, similarly to what was obtained by previous first-principles simulations of soluted methanol. Moreover, the rotational slowing down is compatible with the one predicted on the basis of the excluded volume fraction, which leads to a slower hydrogen bond exchange rate. The analysis of simulations performed at different temperatures suggests that the target temperature of the soluted system must be carefully chosen, in order to avoid artificial slowing-down effects. By generating maximally localized Wannier functions, a detailed description of the polarization effects in both solute and solvent molecules is obtained, which better characterizes the solvation process.

Simulation and Neutron Diffraction Studies of Small Biomolecules in Water

Modern biophysics has benefited greatly from the use of X-ray and neutron diffraction from ordered single crystals of proteins and other macromolecules to give highly detailed pictures of these molecules in the solid state. However, the most biologically relevant environments for these molecules are liquid solutions, and their liquid state properties are sensitive to details of the liquid structuring. The best experimental method for studying such structuring is also neutron diffraction, but of course, the inherent disorder of the liquid state means that these experiments cannot hope to achieve the level of informational detail available from single crystal diffraction. Nonetheless, recent advances in neutron beam intensity, beam stability, and detector sensitivity mean that it should be possible, at least in principle, to use such measurements to extract information about structuring in much more complex systems than have previously been studied. We describe a series of neutron diffraction studies of isotopically labeled molecules in aqueous solution which, when combined with results from computer simulations, can be used to extract conformational information of the hydration of the molecules themselves, essentially opening up new avenues of investigation in structural biology.

Observation of pyridine aggregation in aqueous solution using neutron scattering experiments and MD simulations

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.

Solvation of Biomolecules by the Soft Sticky Dipole-Quadrupole-Octupole Water Model

The soft sticky dipole-quadrupole-octupole (SSDQO) potential energy function represents a water molecule by a single site with a van der Waals sphere and point multipoles. Previously, SSDQO was shown to give good properties for liquid water and solvation of simple ions and is faster than three point models. Here, SSDQO is assessed for solvating biologically relevant molecules having a multi-site, partial charge description. Monte Carlo simulations of ethanol, benzene, and N-methylacetamide in SSDQO with SPC/E moments showed the water structure was as good as in SPC/E. Thus, SSDQO is potentially useful for simulations of biological macromolecules in aqueous solution.

Are there immobilized water molecules around hydrophobic groups? Aqueous solvation of methanol from first principles

Structural, dynamical, bonding, and electronic properties of water molecules around a soluted methanol molecule are studied from first principles. The results are compatible with experiments and qualitatively support the conclusions of recent classical molecular dynamics simulations concerning the controversial issue on the presence of "immobilized" water molecules around hydrophobic groups: the hydrophobic solute slightly reduces the mobility of many surrounding water molecules rather than immobilizing just the few ones which are closest to the methyl group. By generating maximally localized Wannier functions, a detailed description of the polarization effects in both solute and solvent molecules is obtained, which better elucidates the solvation process.

Molecular dynamics studies of the conformation of sorbitol

Molecular dynamics simulations of a 3 molal aqueous solution of D-sorbitol (also called D-glucitol) have been performed at 300 K, as well as at two elevated temperatures to promote conformational transitions. In principle, sorbitol is more flexible than glucose since it does not contain a constraining ring. However, a conformational analysis revealed that the sorbitol chain remains extended in solution, in contrast to the bent conformation found experimentally in the crystalline form. While there are 243 staggered conformations of the backbone possible for this open-chain polyol, only a very limited number were found to be stable in the simulations. Although many conformers were briefly sampled, only eight were significantly populated in the simulation. The carbon backbones of all but two of these eight conformers were completely extended, unlike the bent crystal conformation. These extended conformers were stabilized by a quite persistent intramolecular hydrogen bond between the hydroxyl groups of carbon C-2 and C-4. The conformational populations were found to be in good agreement with the limited available NMR data except for the C-2-C-3 torsion (spanned by the O-2-O-4 hydrogen bond), where the NMR data support a more bent structure.

The hydration of glucose: the local configurations in sugar-water hydrogen bonds

The hydration of a simple sugar is an essential model for understanding interactions between hydrophilic groups and interfacial water molecules. Here I perform first-principles molecular dynamics simulations on a glucose-water system and investigate how individual hydroxyl groups are locally hydrated. I demonstrate that the hydroxyl groups are less hydrated and more incompatible with a locally tetrahedral network of hydrogen bonds than previously thought. The results suggest that the hydroxyl groups form roughly two hydrogen bonds. Further, I find that the local hydration of the hydroxyl groups is sensitively affected by seemingly small variations in the local electronic structure and bond polarity of the groups. My findings offer insight into an atomic-level understanding of sugar-water interactions.

Conformation, dynamics, solvation and relative stabilities of selected β-hexopyranoses in water: a molecular dynamics study with the gromos 45A4 force field

The present article reports long timescale (200 ns) simulations of four beta-D-hexopyranoses (beta-D-glucose, beta-D-mannose, beta-D-galactose and beta-D-talose) using explicit-solvent (water) molecular dynamics and vacuum stochastic dynamics simulations together with the GROMOS 45A4 force field. Free-energy and solvation free-energy differences between the four compounds are also calculated using thermodynamic integration. Along with previous experimental findings, the present results suggest that the formation of intramolecular hydrogen-bonds in water is an 'opportunistic' consequence of the close proximity of hydrogen-bonding groups, rather than a major conformational driving force promoting this proximity. In particular, the conformational preferences of the hydroxymethyl group in aqueous environment appear to be dominated by 1,3-syn-diaxial repulsion, with gauche and solvation effects being secondary, and intramolecular hydrogen-bonding essentially negligible. The rotational dynamics of the exocyclic hydroxyl groups, which cannot be probed experimentally, is found to be rapid (10-100 ps timescale) and correlated (flip-flop hydrogen-bonds interconverting preferentially through an asynchronous disrotatory pathway). Structured solvent environments are observed between the ring and lactol oxygen atoms, as well as between the 4-OH and hydroxymethyl groups. The calculated stability differences between the four compounds are dominated by intramolecular effects, while the corresponding differences in solvation free energies are small. An inversion of the stereochemistry at either C(2) or C(4) from equatorial to axial is associated with a raise in free energy. Finally, the particularly low hydrophilicity of beta-D-talose appears to be caused by the formation of a high-occurrence hydrogen-bonded bridge between the 1,3-syn-diaxial 2-OH and 4-OH groups. Overall, good agreement is found with available experimental and theoretical data on the structural, dynamical, solvation and energetic properties of these compounds. However, this detailed comparison also reveals some discrepancies, suggesting the need (and providing a solid basis) for further refinement.

Nanometer-scale ion aggregates in aqueous electrolyte solutions: guanidinium carbonate

Neutron diffraction with isotopic substitution (NDIS) experiments and molecular dynamics (MD) simulations have been used to characterize the structure of aqueous guanidinium carbonate (Gdm2CO3) solutions. The MD simulations found very strong hetero-ion pairing in Gdm2CO3 solution and were used to determine the best structural experiment to demonstrate this ion pairing. The NDIS experiments confirm the most significant feature of the MD simulation, which is the existence of strong hetero-ion pairing between the Gdm+ and CO3(2-) ions. The neutron structural data also support the most interesting feature of the MD simulation, that the hetero-ion pairing is sufficiently strong as to lead to nanometer-scale aggregation of the ions. The presence of such clustering on the nanometer length scale was then confirmed using small-angle neutron scattering experiments. Taken together, the experiment and simulation suggest a molecular-level explanation for the contrasting denaturant properties of guanidinium salts in solution.

Structure of aqueous glucose solutions as determined by neutron diffraction with isotopic substitution experiments and molecular dynamics calculations

Neutron diffraction with isotopic substitution (NDIS) experiments and molecular dynamics (MD) simulations have been used to examine the structuring of solvent around d-glucose in aqueous solution. As expected, no significant tendency for glucose molecules to aggregate was found in either the experiments or the simulation. To the extent that solute pairing does occur as the result of the high concentration, it was found to take place through hydroxyl-hydroxyl hydrogen bonds, in competition with water molecules for the same hydrogen-bonding sites. A detailed analysis of the hydrogen-bonding patterns occurring in the simulations found that the sugar hydroxyl groups are more efficient hydrogen bond donors than acceptors. From the comparison of the MD and NDIS data, it was found that while the modeling generally does a satisfactory job in reproducing the experimental data the force fields may produce sugar rings that are too rigid and thus may require future revisions.

Spatial hydration structures and dynamics of phenol in sub- and supercritical water

The hydration structures and dynamics of phenol in aqueous solution at infinite dilution are investigated using molecular-dynamics simulation technique. The simulations are performed at several temperatures along the coexistence curve of water up to the critical point, and above the critical point with density fixed at 0.3 g/cm3. The hydration structures of phenol are characterized using the radial, cylindrical, and spatial distribution functions. In particular, full spatial maps of local atomic (solvent) density around a solute molecule are presented. It is demonstrated that in addition to normal H bonds with hydroxyl group of phenol, water forms pi-type complexes with the center of the benzene ring, in which H2O molecules act as H-bond donor. At ambient conditions phenol is solvated by 38 water molecules, which make up a large hydrophobic cavity, and forms on average 2.39 H bonds (1.55 of which are due to the hydroxyl group-water interactions and 0.84 are due to the pi complex) with its hydration shell. As temperature increases, the hydration structure of phenol undergoes significant changes. The disappearance of the pi-type H bonding is observed near the critical point. Self-diffusion coefficients of water and phenol are also calculated. Dramatic increase in the diffusivity of phenol in aqueous solution is observed near the critical point of simple point-charge-extended water and is related to the changes in water structure at these conditions.

Nanometer-Scale Ion Aggregates in Aqueous Electrolyte Solutions: Guanidinium Sulfate and Guanidinium Thiocyanate

Neutron diffraction experiments and molecular dynamics simulations are used to study the structure of aqueous solutions of two electrolytes: guanidinium sulfate (a mild protein conformation stabilizer) and guanidinium thiocyanate (a powerful denaturant). The MD simulations find the unexpected result that in the Gdm2SO4 solution the ions aggregated into mesoscopic (nanometer-scale) clusters, while no such aggregation is found in the GdmSCN solution. The neutron diffraction studies, the most direct experimental probe of solution structure, provide corroborating evidence that the predicted very strong ion pairing does occur in solutions of 1.5 m Gdm2SO4 but not in 3 m solutions of GdmSCN. A mechanism is proposed as to how this mesoscopic solution structure affects solution denaturant properties and suggests an explanation for the Hofmeister ordering of these solutions in terms of this ion pairing and the ability of sulfate to reverse the denaturant power of guanidinium.

Neutron Diffraction and Computer Simulation Studies of d-Xylose

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.

On the Intactness of Hydrogen Bonds around Nonpolar Solutes Dissolved in Water

Angell developed a simple two-state model of hydrogen bonds with the aim to describe some properties of pure water. Muller extended the two-state model description to treat the unusual thermodynamics of hydrophobic hydration. We show here that, to correctly reproduce a qualitative feature of the temperature dependence of the hydration heat capacity change of nonpolar solutes by means of the two-state Muller's model, the hydrogen bonds in the hydration shell have to be more broken than those in bulk water. This contrasts with the suggestion in the literature that more hydrogen bonds form around a nonpolar solute in water.

The Structure of Aqueous Guanidinium Chloride Solutions

The combination of neutron diffraction with isotopic substitution (NDIS) experiments and molecular dynamics (MD) simulations to characterize the structuring in an aqueous solution of the denaturant guanidinium chloride is described. The simulations and experiments were carried out at a concentration of 3 m at room temperature, allowing for an examination of any propensity for ion association in a realistic solution environment. The simulations satisfactorily reproduced the principal features of the neutron scattering and indicate a bimodal hydration of the guanidinium ions, with the N-H groups making well-ordered hydrogen bonds in the molecular plane, but with the planar faces relatively deficient in interactions with water. The most striking feature of these solutions is the rich ion-ion ordering observed around the guanidinium ion in the simulations. The marked tendency of the guanidinium ions to stack parallel to their water-deficient surfaces indicates that the efficiency of this ion as a denaturant is due to its ability to simultaneously interact favorably with both water and hydrophobic side chains of proteins.