Synergy in protein-osmolyte mixtures

Solvent accessible surface area-assessed molecular basis of osmolyte-induced protein stability

In solvent-modulated protein folding, under certain physiological conditions, an equilibrium exists between the unfolded and folded states of the protein without any need to break or make a covalent bond. In this process, interactions between various protein groups (peptides) and solvent molecules are known to play a major role in determining the directionality of the chemical reaction. However, an understanding of the mechanism of action of the co(solvent) by a generic theoretical underpinning is lacking. In this study, a generic solvation model is developed based on statistical mechanics and the thermodynamic transfer free energy model by considering the microenvironment polarity of the interacting co(solvent)-protein system. According to this model, polarity and the fractional solvent-accessible surface areas contribute to the interaction energies. The present model includes various orientations of participating interactant solvent surfaces of suitable areas. As model systems, besides the backbone we consider naturally occurring amino acid residues solvated in ten different osmolytes, small organic compounds known to modulate protein stability. The present model is able to predict the correct trend of the osmolyte-peptide interactions ranging from stabilizing to destabilizing not only for the backbone but also for side chains. Our model predicts Asn, Gln, Asp, Glu, Arg and Pro to be highly stable in most of the protecting osmolytes while Ala, Val, Ile, Leu, Thr, Met, Lys, Phe, Trp and Tyr are predicted to be moderately stable, and Ser, Cys and Histidine are predicted to be least stable. However, in denaturing solvents, both backbone and side chain models show similar stabilities in urea and guanidine. One of the important aspects of this model is that it is essentially parameter-free and consistent with the electrostatics of the interaction partners that make this model suitable for estimating any solute-solvent interaction energies.

Synergistic Solvation as the Enhancement of Local Mixing

Mixing two solvents can sometimes make a much better solvent than expected from their weighted mean. This phenomenon, called synergistic solvation, has commonly been explained via the Hildebrand and Hansen solubility parameters, yet their inability in other solubilization phenomena, most notably hydrotropy, necessitates an alternative route to elucidating solubilization. While, recently, the universal theory of solubilization was founded on the statistical thermodynamic fluctuation theory (as a generalization of the Kirkwood-Buff theory), its demand for experimental data processing has been a hindrance for its wider application. This can be overcome by the solubility isotherm theory, which is founded on the fluctuation theory yet reduces experimental data processing significantly to the level of isotherm analysis in sorption. The isotherm analysis identifies the driving force of synergistic solvation as the enhancement of solvent mixing around the solute, opposite in behavior to hydrotropy (characterized by the enhancement of demixing or self-association around the solute). Thus, the fluctuation theory, including its solubility isotherms, provides a universal language for solubilization across the historic subcategorization of solubilizers, for which different (and often contradictory) mechanistic models have been proposed.

Preferential Solvation Phenomena as Solute-Induced Structure-Making/Breaking Processes: Linking Thermodynamic Preferential Interaction Parameters to Fundamental Structure Making/Breaking Functions

In this work, we identify the explicit macroscopic-to-microscopic rigorous links between existing thermodynamic preferential interaction parameters Γ Q α Q β ( χ i ) and microstructural descriptors based on total correlation function integrals, leading to their unambiguous characterization in terms of fundamental structure making/breaking functions S α β . First, we provide the statistics mechanical framework to identify a universal molecular-based signature for the preferential solvation P S phenomenon involving solutes at infinite dilution in mixed-solvent environments and discuss its fundamental properties. Then, we characterize the S α β functions relevant to the P S process, identify the microscopic markers for the existing preferential interaction parameters Γ Q α Q β ( χ i ) in terms of the S α β functions, and test their compliance with a pair of essential microstructural constraints linked to the properties of the universal P S signature. Moreover, we illustrate the analysis by probing the behavior of a representative ternary system comprising the solubility of methane in aqueous 1,4-dioxane mixed-solvent environments under ambient conditions. Finally, we discuss some relevant issues surrounding the statistical mechanical (microstructural) interpretation of the thermodynamic (macroscopic) preferential interaction parameters, review some pitfalls in their evaluation from molecular simulation, and provide an outlook.

Counteracting Effects of Trimethylamine N-Oxide against Urea in Aqueous Solutions: Insights from Theoretical Two-Dimensional Infrared Spectroscopy

The study of small osmolytes in their aqueous solutions has gained significant attention because of their relevance to structure and thermodynamics of proteins in aqueous media. Special attention has been given to the binary and ternary aqueous solutions of urea and trimethylamine N-oxide (TMAO). Urea is a well-known protein denaturant, while TMAO protects proteins in their native states. Interestingly, TMAO counteracts urea's ability to denature proteins when present in solutions with approximately half of the concentration of urea. Vibrational spectroscopy can improve our understanding of the molecular origin of this counteracting effect because of its sensitivity to local structure and dynamics. We present results of theoretical linear vibrational and two-dimensional infrared (2DIR) spectroscopy of water in the binary and ternary aqueous solutions of TMAO and urea. The 2DIR spectra are calculated using the electronic structure/molecular dynamics approach. The non-Condon effects in spectral transitions are incorporated in the theoretical calculations of 2DIR spectra. It is found that TMAO disrupts the local structure of water, while urea leaves it essentially unaffected. The 2DIR results show that both TMAO and urea slow down the dynamics of spectral diffusion of water. The extent of slowing down is found to be particularly significant for both hydration and bulk water in the presence of TMAO which can be attributed to strong hydrogen bonds between the water and TMAO molecules. The water molecules present in the hydration layer of the solutes in the ternary solutions are found to relax at even slower rates compared to that in their binary solutions in water. The hydrogen bonds between TMAO and urea are found to be not stable. Thus, the counteracting effect of TMAO against urea is seen to take place mainly through water-mediated interactions. Such TMAO-induced effects giving rise to more structured and slower hydrogen-bonded network are successfully captured through 2DIR spectroscopic calculations.

Deep Eutectic Solvents for Efficient Drug Solvation: Optimizing Composition and Ratio for Solubility of β-Cyclodextrin

Deep eutectic solvents (DESs) show promise in pharmaceutical applications, most prominently as excellent solubilizers. Yet, because DES are complex multi-component mixtures, it is challenging to dissect the contribution of each component to solvation. Moreover, deviations from the eutectic concentration lead to phase separation of the DES, making it impractical to vary the ratios of components to potentially improve solvation. Water addition alleviates this limitation as it significantly decreases the melting temperature and stabilizes the DES single-phase region. Here, we follow the solubility of β-cyclodextrin (β-CD) in DES formed by the eutectic 2:1 mole ratio of urea and choline chloride (CC). Upon water addition to DES, we find that at almost all hydration levels, the highest β-CD solubility is achieved at DES compositions that are shifted from the 2:1 ratio. At higher urea to CC ratios, due to the limited solubility of urea, the optimum composition allowing the highest β-CD solubility is reached at the DES solubility limit. For mixtures with higher CC concentration, the composition allowing optimal solvation varies with hydration. For example, β-CD solubility at 40 wt% water is enhanced by a factor of 1.5 for a 1:2 urea to CC mole ratio compared with the 2:1 eutectic ratio. We further develop a methodology allowing us to link the preferential accumulation of urea and CC in the vicinity of β-CD to its increased solubility. The methodology we present here allows a dissection of solute interactions with DES components that is crucial for rationally developing improved drug and excipient formulations.

A simple theory for interfacial properties of dilute solutions

Recent studies suggest that cosolute mixtures may exert significant non-additive effects upon protein stability. The corresponding liquid–vapor interfaces may provide useful insight into these non-additive effects. Accordingly, in this work, we relate the interfacial properties of dilute multicomponent solutions to the interactions between solutes. We first derive a simple model for the surface excess of solutes in terms of thermodynamic observables. We then develop a lattice-based statistical mechanical perturbation theory to derive these observables from microscopic interactions. Rather than adopting a random mixing approximation, this dilute solution theory (DST) exactly treats solute–solute interactions to lowest order in perturbation theory. Although it cannot treat concentrated solutions, Monte Carlo (MC) simulations demonstrate that DST describes the interactions in dilute solutions with much greater accuracy than regular solution theory. Importantly, DST emphasizes a fundamental distinction between the “intrinsic” and “effective” preferences of solutes for interfaces. DST predicts that three classes of solutes can be distinguished by their intrinsic preference for interfaces. While the surface preference of strong depletants is relatively insensitive to interactions, the surface preference of strong surfactants can be modulated by interactions at the interface. Moreover, DST predicts that the surface preference of weak depletants and weak surfactants can be qualitatively inverted by interactions in the bulk. We also demonstrate that DST can be extended to treat surface polarization effects and to model experimental data. MC simulations validate the accuracy of DST predictions for lattice systems that correspond to molar concentrations.

A green and facile method to prepare graphitic carbon nitride nanosheets with outstanding photocatalytic H2O2 production ability via NaClO hydrothermal treatment

Transportation and transformation of photogenerated carriers during the photocatalytic process of graphitic carbon nitride (g-C₃N₄) are restricted by the low surface area and insufficient charge separation efficiency.

Combinations of Osmolytes, Including Monosaccharides, Disaccharides, and Sugar Alcohols Act in Concert During Cryopreservation to Improve Mesenchymal Stromal Cell Survival

There is demand for non-dimethyl sulfoxide (DMSO) cryoprotective agents that maintain cell viability without causing poor postthaw function or systemic toxicity. The focus of this investigation involves expanding our understanding of multicomponent osmolyte solutions and their ability to preserve cell viability during freezing. Controlled cooling rate freezing, Raman microscopy, and differential scanning calorimetry (DSC) were utilized to evaluate the differences in recovery and ice crystal formation behavior for solutions containing multiple cryoprotectants, including sugars, sugar alcohols, and small molecule additives. Postthaw recovery of mesenchymal stem cells (MSCs) in solutions containing multiple osmolytes have been shown to be comparable or better than that of MSCs frozen in 10% DMSO at 1°C/min when the solution composition is optimized. Maximum postthaw recovery was observed in these multiple osmolyte solutions with incubation times of up to 2 h before freezing. Raman images demonstrate large ice crystal formation in cryopreserved cells incubated for shorter periods of time (∼30 min), suggesting that longer permeation times are needed for these solutions. Recovery was dependent upon the concentration of each component in solution, and was not strongly correlated with osmolarity. It is noteworthy that the postthaw recovery varied significantly with the composition of solutions containing the same three components and this variation exhibited an inverted U-shape behavior, indicating that there may be a "sweet spot" for different combinations of osmolytes. Raman images of freezing behavior in different solution compositions were consistent with the observed postthaw recovery. Phase change behavior (solidification patterns and glass-forming tendency) did not differ for solutions with similar osmolarity, but differences in postthaw recovery suggest that biological, not physical, methods of protection are at play. Lastly, molecular substitution of glucose (a monosaccharide) for sucrose (a disaccharide) resulted in a significant drop in recovery. Taken together, the information from these studies increases our understanding of non-DMSO multicomponent cryoprotective solutions and the manner by which they enhance postthaw recovery.