Molecular-scale hydrophobic interactions between hard-sphere reference solutes are attractive and endothermic

How do water-mediated interactions and osmotic second virial coefficients vary with particle size?

We examine quantitatively the solute-size dependences of the effective interactions between nonpolar solutes in water and in a simple liquid. The potential w(r) of mean force and the osmotic second virial coefficients B are calculated with high accuracy from molecular dynamics simulations. As the solute diameter increases from methane's to C60's with the solute-solute and solute-solvent attractive interaction parameters fixed to those for the methane-methane and methane-water interactions, the first minimum of w(r) lowers from -1.1 to -4.7 in units of the thermal energy kT. Correspondingly, the magnitude of B (<0) increases proportional to σα with some power close to 6 or 7, which reinforces the solute-size dependence of B found earlier for a smaller range of σ [H. Naito, R. Okamoto, T. Sumi and K. Koga, J. Chem. Phys., 2022, 156, 221104]. We also demonstrate that the strength of the attractive interactions between solute and solvent molecules can qualitatively change the characteristics of the effective pair interaction between solute particles, both in water and in a simple liquid. If the solute-solvent attractive force is set to be weaker (stronger) than a threshold, the effective interaction becomes increasingly attractive (repulsive) with increasing solute size.

In Silico Study: Combination of α-Mangostin and Chitosan Conjugated with Trastuzumab against Human Epidermal Growth Factor Receptor 2

Breast cancer is a type of cancer with the highest prevalence worldwide. Almost 10–30% of breast cancer cases are diagnosed as positive for HER2 (human epidermal growth factor receptor 2). The currently available treatment methods still exhibit many shortcomings such as a high incidence of side effects and treatment failure due to resistance. This in silico study aims to simulate α-mangostin and chitosan combination conjugated to trastuzumab formulation against HER2 as an effort to improve breast cancer patient therapy. This molecular docking simulation was done through using PatchDock Server. The materials used including the two-dimensional structure of α-mangostin, chitosan, and sodium tripolyphosphate from the PubChem database; trastuzumab FASTA sequence from the DrugBank database; and HER2 structure obtained from a crystal complex with PDB ID: 1N8Z. The results indicated that the particle of α-mangostin and chitosan combinations interacted mostly with the crystallizable fragment (Fc region) of trastuzumab in the conjugation process. The conjugation of trastuzumab to the particle of a combination of α-mangostin and chitosan resulted in the greatest increase in the binding score of the smallest-sized particles (50 Å) with an increase in the score of 3828 and also gave the most similar mode of interaction with trastuzumab. However, the conjugation of trastuzumab eliminated the similarity of the mode of interaction and increased the value of atomic contact energy. Thus, a cominbation of α-mangostin and chitosan conjugated to a trastuzumab formulation was predicted can increase the effectiveness of breast cancer therapy at a relatively small particle size but with the consequence of decreasing atomic contact energy.

The behavior of methane-water mixtures under elevated pressures from simulations using many-body potentials

Non-polarizable empirical potentials have been proven to be incapable of capturing the mixing of methane-water mixtures at elevated pressures. Although density functional theory-based ab initio simulations may circumvent this discrepancy, they are limited in terms of the relevant time and length scales associated with mixing phenomena. Here, we show that the many-body MB-nrg potential, designed to reproduce methane-water interactions with coupled cluster accuracy, successfully captures this phenomenon up to 3 GPa and 500 K with varying methane concentrations. Two-phase simulations and long time scales that are required to fully capture the mixing, affordable due to the speed and accuracy of the MBX software, are assessed. Constructing the methane-water equation of state across the phase diagram shows that the stable mixtures are denser than the sum of their parts at a given pressure and temperature. We find that many-body polarization plays a central role, enhancing the induced dipole moments of methane by 0.20 D during mixing under pressure. Overall, the mixed system adopts a denser state, which involves a significant enthalpic driving force as elucidated by a systematic many-body energy decomposition analysis.

Hydrophobic hydration and pairwise hydrophobic interaction of Lennard-Jones and Mie particles in different water models

The study provides a deep computational analysis of the thermodynamic and structural features associated with the hydration of xenon, Xe, and its pairwise hydrophobic interaction (i.e., the potential of mean force, PMF), over a large temperature range. Xe is described both as a Lennard-Jones particle, LJ-Xe, and as a Mie particle, Mie-Xe (pseudo hard sphere). Three different water models are used: TIP3P-Ew, SPCE and TIP4P-2005. Mie-Xe is more hydrophobic than LJ-Xe due to the lack of the attractive energetic interactions with water molecules; its hydration, around room temperature, is opposed by a large and negative entropy change and a positive enthalpy change. The PMF of Mie-Xe is characterized by a deep minimum at contact distance whose depth increases with temperature, and whose magnitude is significantly larger than that obtained for LJ-Xe. The contact minimum configuration of Mie-Xe is favoured by a large positive entropy change and contrasted by a positive enthalpy change. These results are qualitatively the same regardless of the water model used. There is no clear connection between the values determined for the thermodynamic functions and the structural features of the hydration shells surrounding the single Mie-Xe and the couple of Mie-Xe particles in the contact minimum configuration. This confirms that the structural reorganization of water associated with such processes is characterized by an almost complete enthalpy-entropy compensation.

Short solvent model for ion correlations and hydrophobic association

Coulomb interactions play a major role in determining the thermodynamics, structure, and dynamics of condensed-phase systems, but often present significant challenges. Computer simulations usually use periodic boundary conditions to minimize corrections from finite cell boundaries but the long range of the Coulomb interactions generates significant contributions from distant periodic images of the simulation cell, usually calculated by Ewald sum techniques. This can add significant overhead to computer simulations and hampers the development of intuitive local pictures and simple analytic theory. In this paper, we present a general framework based on local molecular field theory to accurately determine the contributions from long-ranged Coulomb interactions to the potential of mean force between ionic or apolar hydrophobic solutes in dilute aqueous solutions described by standard classical point charge water models. The simplest approximation leads to a short solvent (SS) model, with truncated solvent-solvent and solute-solvent Coulomb interactions and long-ranged but screened Coulomb interactions only between charged solutes. The SS model accurately describes the interplay between strong short-ranged solute core interactions, local hydrogen-bond configurations, and long-ranged dielectric screening of distant charges, competing effects that are difficult to capture in standard implicit solvent models.

Characterizing Solvent Density Fluctuations in Dynamical Observation Volumes

Hydrophobic effects drive diverse aqueous assemblies, such as micelle formation or protein folding, wherein the solvent plays an important role. Consequently, characterizing the free energetics of solvent density fluctuations can lead to important insights into these processes. Although techniques such as the indirect umbrella sampling (INDUS) method can be used to characterize solvent fluctuations in static observation volumes of various sizes and shapes, characterizing how the solvent mediates inherently dynamic processes, such as self-assembly or conformational change, remains a challenge. In this work, we generalize the INDUS method to facilitate the enhanced sampling of solvent fluctuations in dynamical observation volumes, whose positions and shapes can evolve. We illustrate the usefulness of this generalization by characterizing water density fluctuations in dynamical volumes pertaining to the hydration of flexible solutes, the assembly of small hydrophobes, and conformational transitions in a model peptide. We also use the method to probe the dynamics of hard spheres.

Connections between the Anomalous Volumetric Properties of Alcohols in Aqueous Solution and the Volume of Hydrophobic Association

The partial molar volumes of alcohols in water exhibit a non-monotonic dependence on concentration at room temperature, initially decreasing with increasing concentration before passing through a minimum and rising to the pure liquid plateau. This anomalous behavior is associated with hydrophobic interactions. We report molecular simulations of short chain alcohols and alkanes in water to examine the volumetric properties of these mixtures at infinite dilution over a range of temperatures. Our simulations find this anomaly disappears at a crossover temperature, above which the solute volume only varies monotonically with concentration. A Voronoi volume analysis of solution configurations finds that solutes in clusters take up less space than individual solutes at low temperature and more space at elevated temperatures. These changes in cluster volumes are subsequently shown to correlate with the derivative of the solute partial molar volume with respect to solute concentration. The changes in solute volume upon nonpolar solute association impact the response of molecular-scale hydrophobic interactions for assembly with increasing pressure.

Drug Partitioning in Micellar Media and Its Implications in Rational Drug Design: Insights with Streptomycin

Oral bioavailability of a drug molecule requires its effective delivery to the target site. In general, majority of synthetically developed molecular entities have high hydrophobic nature as well as low bioavailability, therefore the need for suitable delivery vehicles arises. Self-assembled structures such as micelles, niosomes, and liposomes have been used as effective delivery vehicles and studied extensively. However, the information available in literature is mostly qualitative in nature. We have quantitatively investigated the partitioning of antibiotic drug streptomycin into cationic, nonionic, and a mixture of cationic and nonionic surfactant micelles and its interaction with the transport protein serum albumin upon subsequent delivery. A combination of calorimetry and spectroscopy has been used to obtain the thermodynamic signatures associated with partitioning and interaction with the protein and the resulting conformational changes in the latter. The results have been correlated with other class of drugs of different nature to understand the role of molecular features in the partitioning process. These studies are oriented toward understanding the physical chemistry of partitioning of a variety of drug molecules into suitable delivery vehicles and hence establishing structure-property-energetics relationships. Such studies provide general guidelines toward a broader goal of rational drug design.

Molecular Shape and the Hydrophobic Effect

This review focuses on papers published since 2000 on the topic of the properties of solutes in water. More specifically, it evaluates the state of the art of our understanding of the complex relationship between the shape of a hydrophobe and the hydrophobic effect. To highlight this, we present a selection of references covering both empirical and molecular dynamics studies of small (molecular-scale) solutes. These include empirical studies of small molecules, synthetic hosts, crystalline monolayers, and proteins, as well as in silico investigations of entities such as idealized hard and soft spheres, small solutes, hydrophobic plates, artificial concavity, molecular hosts, carbon nanotubes and spheres, and proteins.

Hydrophobicity of proteins and nanostructured solutes is governed by topographical and chemical context

Hydrophobic interactions drive many important biomolecular self-assembly phenomena. However, characterizing hydrophobicity at the nanoscale has remained a challenge due to its nontrivial dependence on the chemistry and topography of biomolecular surfaces. Here we use molecular simulations coupled with enhanced sampling methods to systematically displace water molecules from the hydration shells of nanostructured solutes and calculate the free energetics of interfacial water density fluctuations, which quantify the extent of solute-water adhesion, and therefore solute hydrophobicity. In particular, we characterize the hydrophobicity of curved graphene sheets, self-assembled monolayers (SAMs) with chemical patterns, and mutants of the protein hydrophobin-II. We find that water density fluctuations are enhanced near concave nonpolar surfaces compared with those near flat or convex ones, suggesting that concave surfaces are more hydrophobic. We also find that patterned SAMs and protein mutants, having the same number of nonpolar and polar sites but different geometrical arrangements, can display significantly different strengths of adhesion with water. Specifically, hydroxyl groups reduce the hydrophobicity of methyl-terminated SAMs most effectively not when they are clustered together but when they are separated by one methyl group. Hydrophobin-II mutants show that a charged amino acid reduces the hydrophobicity of a large nonpolar patch when placed at its center, rather than at its edge. Our results highlight the power of water density fluctuations-based measures to characterize the hydrophobicity of nanoscale surfaces and caution against the use of additive approximations, such as the commonly used surface area models or hydropathy scales for characterizing biomolecular hydrophobicity and the associated driving forces of assembly.

Synergistic Inhibition of Protein Fibrillation by Proline and Sorbitol: Biophysical Investigations

We report here interesting synergistic effects of proline and sorbitol, two well-known chemical chaperones, in the inhibition of fibrillation of two proteins, insulin and lysozyme. A combination of many biophysical techniques has been used to understand the structural morphology and modes of interaction of the chaperones with the proteins during fibrillation. Both the chaperones establish stronger polar interactions in the elongation and saturation stages of fibrillation compared to that in the native stage. However, when presented as a mixture, we also see contribution of hydrophobic interactions. Thus, a co-operative adjustment of polar and hydrophobic interactions between the chaperones and the protein surface seems to drive the synergistic effects in the fibrillation process. In insulin, this synergy is quantitatively similar in all the stages of the fibrillation process. These observations would have significant implications for understanding protein folding concepts, in general, and for designing combination therapies against protein fibrillation, in particular.

Statistical Analyses of Hydrophobic Interactions: A Mini-Review

This review focuses on the striking recent progress in solving for hydrophobic interactions between small inert molecules. We discuss several new understandings. First, the inverse temperature phenomenology of hydrophobic interactions, i.e., strengthening of hydrophobic bonds with increasing temperature, is decisively exhibited by hydrophobic interactions between atomic-scale hard sphere solutes in water. Second, inclusion of attractive interactions associated with atomic-size hydrophobic reference cases leads to substantial, nontrivial corrections to reference results for purely repulsive solutes. Hydrophobic bonds are weakened by adding solute dispersion forces to treatment of reference cases. The classic statistical mechanical theory for those corrections is not accurate in this application, but molecular quasi-chemical theory shows promise. Finally, because of the masking roles of excluded volume and attractive interactions, comparisons that do not discriminate the different possibilities face an interpretive danger.

Necessity of capillary modes in a minimal model of nanoscale hydrophobic solvation

Modern theories of the hydrophobic effect highlight its dependence on length scale, emphasizing the importance of interfaces in the vicinity of sizable hydrophobes. We recently showed that a faithful treatment of such nanoscale interfaces requires careful attention to the statistics of capillary waves, with significant quantitative implications for the calculation of solvation thermodynamics. Here, we show that a coarse-grained lattice model like that of Chandler [Chandler D (2005)Nature437(7059):640-647], when informed by this understanding, can capture a broad range of hydrophobic behaviors with striking accuracy. Specifically, we calculate probability distributions for microscopic density fluctuations that agree very well with results of atomistic simulations, even many SDs from the mean and even for probe volumes in highly heterogeneous environments. This accuracy is achieved without adjustment of free parameters, because the model is fully specified by well-known properties of liquid water. As examples of its utility, we compute the free-energy profile for a solute crossing the air-water interface, as well as the thermodynamic cost of evacuating the space between extended nanoscale surfaces. These calculations suggest that a highly reduced model for aqueous solvation can enable efficient multiscale modeling of spatial organization driven by hydrophobic and interfacial forces.

Hydrophobic Ambivalence: Teetering on the Edge of Randomness

Processes ranging from oil-water phase separation to the formation of solid clathrate hydrates send mixed messages regarding whether oil molecules hate or love to be surrounded by water. Recent experimental and theoretical results help decipher these mixed messages by illuminating the conditions under which the stability of a hydrophobic contact is expected to exceed thermal energy fluctuations - thus facilitating hydrophobic self-assembly and the emergence of structure from randomness. Important open questions remain regarding the dependence of hydrophobic interactions on molecular size and temperature, as well as the balance of direct and water-mediated interactions.