On the nature of the molecular ordering of water in aqueous DMSO mixtures

Conformational Landscape of NADH and Ion Binding in Water/DMSO Mixtures via 31P NMR Spectroscopy

This study reports on the conformational states of nicotinamide adenine dinucleotide (NADH) in water/DMSO mixtures and examines the influence of ion binding. We observe evidence of conformational changes through a series of NMR techniques, including 31P NMR relaxation (R1 and R2), chemical exchange saturation transfer (CEST), and diffusion-ordered spectroscopy (DOSY) experiments. The observed variation of the conformational states is indicative of the solvent's influence on NADH's structural flexibility. The experimental findings, in combination with viscosity data, are shown to be in line with findings from earlier molecular dynamics studies. The reported observations emphasize the critical role of the solvent environment in determining the conformational states of NADH and similar molecules with relevance for the biophysiological context. The results found herein can help in studying the biophysical behavior of NADH and similar biomolecules and their associated metabolic pathways under various solvent conditions.

Concentration Fluctuation/Microheterogeneity Duality Illustrated with Aqueous 1,4-Dioxane Mixtures

The structural properties of aqueous 1-4 dioxane mixtures are studied by computer simulations of different water and dioxane force field models, from the perspective of illustrating the link between structural properties at the molecular level and measurable properties such as radiation scattering intensities and Kirkwood-Buff integrals (KBIs). A strategy to consistently correct the KBI obtained from simulations is proposed, which allows us to obtain the genuine KBI corresponding to a given pair of molecular species, in the entire concentration range, and without necessitating excessively large system sizes. The application of this method to the aqueous dioxane mixtures, with an all-atom CHARMM dioxane model and 2 water models, namely, SPC/E and TIP3P, allows one to understand the differences in the structure of the corresponding mixtures at the molecular level, particularly concerning the role of the water aggregates and its model dependence. This study allows us to characterize the dual role played by the concentration fluctuations and the domain segregation, particularly in what concerns the calculated X-ray spectra.

Controlling Interfacial Hydrogen Bonding at a Gold Surface: The Effect of Organic Cosolvents

Water often serves as both a reactant and solvent in electrocatalytic reactions. Interfacial water networks can affect the transport and kinetics of these reactions, e.g., hydrogen evolution reaction and CO2 reduction reaction. Adding cosolvents that influence the hydrogen-bonding (H-bonding) environment, such as dimethyl sulfoxide (DMSO), has the potential to tune the reactivity of these important electrocatalytic reactions by regulating the interfacial local environment and water network. We investigate interfacial H-bonding networks in water-DMSO cosolvent mixtures on gold surfaces by using surface-enhanced infrared absorption spectroscopy and molecular dynamics simulations. Experiments and simulations show that the gold surface is enriched with dehydrated DMSO molecules and the mixture phase-separates to form water clusters. Simulations show a "buckled" water conformation at the surface, further constraining interfacial H-bonding. The small size of these water clusters and the energetically unfavorable H-bond conformations might inhibit H-bonding with bulk water, suppressing the proton diffusion required for efficient hydrogen evolution reaction processes.

Supramolecular Protein Stabilization with Zwitterionic Polypeptide-Cucurbit[7]uril Conjugates

Protein aggregation is an obstacle for the development of new biopharmaceuticals, presenting challenges in shipping and storage of vital therapies. Though a variety of materials and methods have been explored, the need remains for a simple material that is biodegradable, nontoxic, and highly efficient at stabilizing protein therapeutics. In this work, we investigated zwitterionic polypeptides prepared using a rapid and scalable polymerization technique and conjugated to a supramolecular macrocycle host, cucurbit[7]uril, for the ability to inhibit aggregation of model protein therapeutics insulin and calcitonin. The polypeptides are based on the natural amino acid methionine, and zwitterion sulfonium modifications were compared to analogous cationic and neutral structures. Each polymer was end-modified with a single cucurbit[7]uril macrocycle to afford supramolecular recognition and binding to terminal aromatic amino acids on proteins. Only conjugates prepared from zwitterionic structures of sufficient chain lengths were efficient inhibitors of insulin aggregation and could also inhibit aggregation of calcitonin. This polypeptide exhibited no cytotoxicity in human cells even at concentrations that were five-fold of the intended therapeutic regime. We explored treatment of the zwitterionic polypeptides with a panel of natural proteases and found steady biodegradation as expected, supporting eventual clearance when used as a protein formulation additive.

Nanostructured Temperature Indicator for Cold Chain Logistics

Food, chemicals, agricultural products, drugs, and vaccines should be transported and stored within an appropriate low-temperature range, following cold chain logistics. Violations of the required temperature regime are generally reported by time-temperature indicators; however, current sensors do not cover a sufficiently broad low-temperature range and may lack thermal and photostability. Here, we report a nanostructured solvatochromic temperature indicator formed from cellulose nanocrystals decorated with carbon dots (C-dots). The indicator utilizes a strong nonlinear dependence of photoluminescence of C-dots on the composition of water/dimethyl sulfoxide (DMSO) solvent and a composition-dependent variation of the melting temperature of the water/DMSO mixture. Exceeding the temperature of the frozen mixed solvent above a designated threshold value results in solvent melting, flow, and impregnation of the nanostructured film, thus causing an irreversible change in the intensity and wavelength of photoluminescence emission of the film, which is reported both qualitatively and quantitatively. The indicator covers a temperature range from -68 to +19 °C and is cost-efficient, portable and photo- and thermostable.

Lifetime distribution of clusters in binary mixtures involving hydrogen bonding liquids

Hydrogen bonded liquids are associated liquids and tend to exhibit local inhomogeneity in the form of clusters and segregated sub-nano domains. It is an open question as to whether Hbonded clusters in pure water have common features with the water segregated pockets observed in various aqueous binary mixtures, such as water–alcohol mixtures, for example. In the present study, we demonstrate through classical molecular dynamics studies of the lifetime distributions of the hydrogen bonds in different types of binary mixtures, that these lifetimes exhibit the same universal features in the case of the pure liquids, independently of the species concentrations. The same types of three distinct lifetimes are observed, all of them in the sub picosecond regime. The primary lifetime concerns that of Hbonded dimers, and strongly depends on Hbonding criteria such as the bonding distance. The two others are independent of bonding criteria and appear as universal accross many liquids and mixtures. The secondary lifetime (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau _1 \approx 20$$\end{document}τ1≈20 fs) concerns Hbonded cluster lifetimes, while the tertiary lifetime (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau _2 \approx 50$$\end{document}τ2≈50 fs) concerns the topology of these clusters, such as chains or globules, for example. This surprizing separation in three distinct lifetimes suggests the existence of associated three distinct kinetic mechanisms in the very short sub-picosecond time scales, with, in addition, an appealing connection to the concepts of local energy and entropy.

Molecular heterogeneity in aqueous cosolvent systems

Aqueous cosolvent systems (ACoSs) are mixtures of small polar molecules such as amides, alcohols, dimethyl sulfoxide, or ions in water. These liquids have been the focus of fundamental studies due to their complex intermolecular interactions as well as their broad applications in chemistry, medicine, and materials science. ACoSs are fully miscible at the macroscopic level but exhibit nanometer-scale spatial heterogeneity. ACoSs have recently received renewed attention within the chemical physics community as model systems to explore the relationship between intermolecular interactions and microscopic liquid-liquid phase separation. In this perspective, we provide an overview of ACoS spatial segregation, dynamic heterogeneity, and multiscale relaxation dynamics. We describe emerging approaches to characterize liquid microstructure, H-bond networks, and dynamics using modern experimental tools combined with molecular dynamics simulations and network-based analysis techniques.

Microscopic origin of breakdown of Stokes-Einstein relation in binary mixtures: Inherent structure analysis

Aqueous binary mixtures often exhibit dramatic departure from the predicted hydrodynamic behavior when transport properties are plotted against composition. We show by inherent structure (IS) analysis that this sharp composition dependent breakdown of the Stokes-Einstein relation can be attributed to the non-monotonic variation in the average inherent structure energy of these mixtures. Further IS analysis reveals the existence of a unique ground state, stabilized by both the formation of an optimum number of H-bonds and a favorable hydrophobic interaction at this composition. The surprisingly sharp turnaround behavior observed in the effective hydrodynamic radius also owes its origin to the same combination of these two factors. Interestingly, the temperature dependence of isothermal compressibility shows a minimum at the particular composition. Extensive studies on water-dimethyl sulfoxide and water-ethanol mixtures using two different force-fields of water reveal many features that are nearly universal. A justification of this quasi-universal behavior is provided in terms of a mode-coupling theory (MCT) of viscosity, which can serve as the starting point of a remarkable correlation observed with the nearest neighbor structure, as captured by the first peaks of the radial distribution function, and the slowdown in the intermediate scattering function at intermediate wavenumbers. Therefore, the formation of the local structure captured through IS analysis can be correlated with the MCT.

Liquid-Liquid Phase Separation Produces Fast H-Bond Dynamics in DMSO-Water Mixtures

Liquid-liquid phase separation is common in complex mixtures, but the behavior of nanoconfined liquids is poorly understood from a physical perspective. Dimethyl sulfoxide (DMSO) is an amphiphilic molecule with unique concentration-dependent bulk properties in mixtures with water. Here, we use ultrafast two-dimensional infrared (2D IR) spectroscopy to measure the H-bond dynamics of two probe molecules with different polarities: formamide (FA) and dimethylformamide (DMF). Picosecond H-bond dynamics are fastest in the intermediate concentration regime (20-50 mol % DMSO), because such confined water exhibits bulk-like dynamics. Each vibrational probe experiences a unique microscopic environment as a result of nanoscale phase separation. Molecular dynamics simulations show that the dynamics span multiple time scales, from femtoseconds to nanoseconds. Our studies suggest a previously unknown liquid environment, which we label "local bulk", in which despite the local heterogeneity, the ultrafast H-bond dynamics are similar to bulk water.

Molecular Mechanism of Water Reorientation Dynamics in Dimethyl Sulfoxide Aqueous Mixtures

Nonmonotonic composition dependence is often observed for numerous properties in the aqueous mixtures of small amphiphilic molecules. The molecular picture underlying this structure-activity relationship, however, remains largely elusive. We herein studied water reorientation dynamics in the aqueous mixture of dimethyl sulfoxide (DMSO), which has a significant nonmonotonic composition dependence, using molecular dynamic simulation and an extended molecular jump model. The analysis indicates that this nonideal behavior is driven by the collective frame diffusion component of water reorientation, which decelerates in the water-rich regime because of the strengthened hydrogen bonds and accelerates in the water-poor regime as the hydrogen bonding network is broken into smaller aggregates. The current work therefore connects the microheterogeneity in the solvation structure of DMSO-water with its nonmonotonic hydration dynamics and sheds new light on how microsegregation leads to the multiscale hydration nonideality in general.

Empirical S=O stretch vibrational frequency map

Dimethyl sulfoxide (DMSO) water mixtures have been widely studied due to their unique concentration-dependent bulk properties. Here, we present an empirical spectroscopic map for the sulfinyl (S=O) stretching mode. The model can be used to interpret infrared (IR) absorption and ultrafast two-dimensional infrared (2D IR) spectra and quantify hydrogen bond populations and lifetimes by directly connecting spectroscopic measurements with structures and dynamics from molecular dynamics simulations. The electrostatic map is directly parameterized against experimental absorption spectra in the S=O stretching region (980-1100 cm^-1) of dilute DMSO in water. A comparison of center peak frequencies shows that the map performs well across the entire DMSO concentration range, accurately reproducing the ∼10 cm^-1 red-shift per hydrogen bond observed in the experiments. We further benchmark the map by comparing experimental and simulated 2D IR spectra generated by direct numerical integration of the Schrödinger equation. We expect that this empirical frequency map will provide a quantitative platform for investigating intermolecular interactions, microscopic heterogeneity, and ultrafast dynamics in complex liquid mixtures containing DMSO.

Inversion of thermodiffusive properties of ionic colloidal dispersions in water-DMSO mixtures probed by forced Rayleigh scattering

Thermodiffusion properties at room temperature of colloidal dispersions of hydroxyl-coated nanoparticles (NPs) are probed in water, in dimethyl sulfoxide (DMSO) and in mixtures of water and DMSO at various proportions of water, [Formula: see text]. In these polar solvents, the positive NPs superficial charge imparts the systems with a strong electrostatic interparticle repulsion, slightly decreasing from water to DMSO, which is here probed by Small Angle Neutron Scattering and Dynamic Light Scattering. However if submitted to a gradient of temperature, the NPs dispersed in water with ClO₄^- counterions present a thermophilic behavior, the same NPs dispersed in DMSO with the same counterions present a thermophobic behavior. Mass diffusion coefficient [Formula: see text] and Ludwig-Soret coefficient [Formula: see text] are measured as a function of NP volume fraction [Formula: see text] at various [Formula: see text]. The [Formula: see text]-dependence of [Formula: see text] is analyzed in terms of thermoelectric and thermophoretic contributions as a function of [Formula: see text]. Using two different models for evaluating the Eastman entropy of transfer of the co- and counterions in the mixtures, the single-particle thermophoretic contribution (the NP's Eastman entropy of transfer) is deduced. It is found to evolve from negative in water to positive in DMSO. It is close to zero on a large range of [Formula: see text] values, meaning that in this [Formula: see text]-range [Formula: see text] largely depends on the thermoelectric effect of free co- and counterions.

Microscopic origin of the scattering pre-peak in aqueous propylamine mixtures: X-ray and neutron experiments versus simulations

The structure of aqueous propylamine mixtures is investigated through X-ray and neutron scattering experiments, and the scattered intensities compared with computer simulation data.

How do glycerol and dimethyl sulphoxide affect local tetrahedral structure of water around a nonpolar solute at low temperature? Importance of preferential interaction

Glycerol and dimethyl sulphoxide (DMSO) have vital roles in cryoprotection of living cells, tissues, etc. The above action has been directly linked with disruption of hydrogen (H-) bond structure and dynamics of water by these cosolvents at bulk region and around various complex units, such as peptide, amino acid, protein, and lipid membrane. However, the disruption of the local structure of the water solvent around a purely hydrophobic solute is still not studied extensively. The latter is also important in the context of stabilization of protein from cold denaturation. Through all-atom molecular dynamics simulation, we have investigated the comparative effect of glycerol and DMSO on the orientational order of water around a nonpolar solute at -5 °C. A steady reduction of the tetrahedral order of water is observed at bulk (>10 Å distance from the solute) and solute interface (<5.5 Å distance from the solute) with increasing the cosolvent concentration. Contrasting roles of glycerol and DMSO have been evidenced. While DMSO affects the H-bond structure of the interfacial water more than that of the bulk water, glycerol affects the water structure almost uniformly at all regions around the solute. Furthermore, while glycerol helps to retain water molecules at the interface, DMSO significantly reduces the water content in that region. We have put forward a plausible mechanism for these contrasting roles of these cosolvents. The solute-cosolvent hydrophobic-interaction-induced orientational alignment of an interfacial cosolvent molecule determines whether the involvement of the cosolvent molecules in H-bonding with solvent water in the interface is akin to the bulk region or not.

Preferential solvation and ion association properties in aqueous dimethyl sulfoxide solutions

We study the solvation and the association properties of ion pairs in aqueous dimethyl sulfoxide (DMSO) solution by atomistic molecular dynamics (MD) simulations. The ion pair is composed of two lithium and a single sulfonated diphenyl sulfone ion whose properties are studied under the influence of different DMSO concentrations. For increasing mole fractions of DMSO, we observe a non-ideal behavior of the solution as indicated by the derivatives of the chemical activity. Our findings are complemented by dielectric spectra, which also verify a complex DMSO-water mixing behavior. In agreement with these results, further simulation outcomes reveal an aqueous homoselective solvation of the ion species which fosters the occurrence of pronounced ion association constants at higher DMSO mole fractions. The consequences of this finding are demonstrated by lower ionic conductivities for increasing concentrations of DMSO.

Role of Hydrogen Bonding of Cyclodextrin-Drug Complexes Probed by Thermodiffusion

Temperature gradient-induced migration of biomolecules, known as thermophoresis or thermodiffusion, changes upon ligand binding. In recent years, this effect has been used to determine protein-ligand binding constants. The mechanism through which thermodiffusive properties change when complexes are formed, however, is not understood. An important contribution to thermodiffusive properties originates from the thermal response of hydrogen bonds. Because there is a considerable difference between the degree of solvation of the protein-ligand complex and its isolated components, ligand-binding is accompanied by a significant change in hydration. The aim of the present work is therefore to investigate the role played by hydrogen bonding on the change in thermodiffusive behavior upon ligand-binding. As a model system, we use cyclodextrins (CDs) and acetylsalicylic acid (ASA), where quite a significant change in hydration is expected and where no conformational changes occur when a CD/ASA complex is formed in aqueous solution. Thermophoresis was investigated in the temperature range of 10-50 °C by infrared thermal diffusion forced Rayleigh scattering. Nuclear magnetic resonance measurements were performed at 25 °C to obtain information about the structure of the complexes. All CD/ASA complexes show a stronger affinity toward regions of lower temperature compared to the free CDs. We found that the temperature sensitivity of thermophoresis correlates with the 1-octanol/water partition coefficient. This observation not only establishes the relation between thermodiffusion and degree of hydrogen bonding but also opens the possibility to relate thermodiffusive properties of complexes to their partition coefficient, which cannot be determined otherwise. This concept is especially interesting for protein-ligand complexes where the protein undergoes a conformational change, different from the CD/ASA complexes, giving rise to additional changes in their hydrophilicity.

Charge ordering and scattering pre-peaks in ionic liquids and alcohols

The structural properties of ionic liquids and alcohols are viewed under the charge ordering process as a common basis to explain the peculiarity of their radiation scattering properties, namely the presence, or absence, of a scattering pre-peak.

A new approach for estimating the density of liquids

We propose a novel approach with which to estimate the density of liquids. The approach is based on the assumption that the systems would be structurally similar when viewed at around the length scale (inverse wavenumber) of the first peak of the structure factor, unless their thermodynamic states differ significantly. The assumption was implemented via a similarity transformation to the radial distribution function to extract the density from the structure factor of a reference state with a known density. The method was first tested using two model liquids, and could predict the densities within an error of several percent unless the state in question differed significantly from the reference state. The method was then applied to related real liquids, and satisfactory results were obtained for predicted densities. The possibility of applying the method to amorphous materials is discussed.

Accumulation of formamide in hydrothermal pores to form prebiotic nucleobases

Formamide is one of the important compounds from which prebiotic molecules can be synthesized, provided that its concentration is sufficiently high. For nucleotides and short DNA strands, it has been shown that a high degree of accumulation in hydrothermal pores occurs, so that temperature gradients might play a role in the origin of life [Baaske P, et al. (2007)Proc Natl Acad Sci USA104(22):9346-9351]. We show that the same combination of thermophoresis and convection in hydrothermal pores leads to accumulation of formamide up to concentrations where nucleobases are formed. The thermophoretic properties of aqueous formamide solutions are studied by means of Infrared Thermal Diffusion Forced Rayleigh Scattering. These data are used in numerical finite element calculations in hydrothermal pores for various initial concentrations, ambient temperatures, and pore sizes. The high degree of formamide accumulation is due to an unusual temperature and concentration dependence of the thermophoretic behavior of formamide. The accumulation fold in part of the pores increases strongly with increasing aspect ratio of the pores, and saturates to highly concentrated aqueous formamide solutions of ∼85 wt% at large aspect ratios. Time-dependent studies show that these high concentrations are reached after 45-90 d, starting with an initial formamide weight fraction of[Formula: see text]wt % that is typical for concentrations in shallow lakes on early Earth.

From solutions to molecular emulsions

Concentration fluctuations play an important role in the statistical description of the stability of liquids, particularly in the neighborhood of phase transitions. Classical thermodynamics is blind to fluctuations, and statistical thermodynamics is required to fully understand quantities such as the isothermal compressibility or heat capacity, by linking them to fluctuations of appropriate statistical microscopic quantities and showing that they are response functions. This is illustrated by the seminal Kirkwood–Buff theory of solutions. However, the existence of micro-heterogeneous structures, particularly in aqueous mixtures, which leads to large Kirkwood–Buff integrals, suggest that micro-heterogeneity is a form of concentration fluctuation. This interpretation becomes difficult to accept when extrapolated to larger micro-heterogeneous structures such as micellar aggregates in micro-emulsions. By analyzing how different methods, experimental, computer experiments and theoretical approaches deal with the underlying duality behind these two physical manifestations, we put in evidence the need to reconsider the description of liquids by incorporating the description of emergent “objects”, such as the micro-heterogeneous structures from a molecular point of view. On this path, the concept of “molecular emulsion” allows to describe in a unified way all type of disordered liquids, from solutions to the organized liquids of soft matter.

Micro-heterogeneity versus clustering in binary mixtures of ethanol with water or alkanes

Snapshots of the difference in complex disorder, with analogy with direct (ethanol–water) and inverse (ethanol–alkanes) emulsions.

Comparison of tetrel bonds and halogen bonds in complexes of DMSO with ZF3X (Z = C and Si; X = halogen)

A theoretical study of the complexes formed by dimethylsulfoxide (DMSO) with ZF₃X (Z = C and Si; X = halogen) has been performed at the MP2/aug-cc-pVTZ level.