A cluster size distribution theory to study the thermodynamics and phase behavior of multi-bonding single site solutes in patchy colloidal mixtures.

We study binary mixtures of multi-bonding single site solute particles in a solvent comprising patchy colloid particles. The particles in the mixture interact by very short-ranged attraction and hard-sphere repulsion. The attractive patch on the solute can bond with multiple solvent particles, whereas the patch on the solvent is restricted to bond only once. From a quasi-chemical analysis of association, in the hard-sphere reference we develop an accurate multi-body correlation information for the distribution of solvent particles over the patch region of the solute. We use this information within Wertheim's multi-density formalism to develop a cluster size distribution theory that is capable of capturing the physics of multi-body association for any geometry of association sites on the solute. We use this general framework to study a mixture containing Janus solutes and one- or two-patch solvent particles over a range of concentration of the solute and association strengths. We find that a mixture of two-patch solvent (with both patches of the same kind) and multi-bonding solutes with different patch geometries can have a vapor-liquid equilibrium, although the pure components themselves cannot phase separate. The liquid state occurs at very low densities, forming a so-called empty liquid. For the relative association strengths studied in this work, we observe that the vapor-liquid coexistence curve broadens as the concentration of the patchy solvent particles in the liquid phase is increased. The pressure-composition phase equilibrium curves show negative azeotropes for these mixtures. We also observe that, for these mixtures, as the size of the patch on the solute particles is decreased, the critical temperature and the critical packing fraction decreases.

Solvophobic and solvophilic contributions in the water-to-aqueous guanidinium chloride transfer free energy of model peptides.

We study the solvation free energy of two different conformations (helix and extended) of two different peptides (deca-alanine and deca-glycine) in two different solvents (water and aqueous guanidinium chloride, GdmCl). The free energies are obtained using the quasichemical organization of the potential distribution theorem, an approach that naturally provides the repulsive (solvophobic or cavity) and attractive (solvophilic) contributions to solvation. The solvophilic contribution is further parsed into a chemistry contribution arising from solute interaction with the solvent in the first solvation shell and a long-range contribution arising from non-specific interactions between the solute and the solvent beyond the first solvation shell. The cavity contribution is obtained for two different envelopes, ΣSE, which theory helps identify as the solvent excluded volume, and ΣG, a larger envelope beyond which solute-solvent interactions are Gaussian. The ΣSE envelope is independent of the solvent, as expected on the basis of the insensitivity to the solvent type of the distance of closest approach between protein heavy atoms and solvent heavy atoms, but contrary to the intuition based on treating solvent constituents as spheres of some effective radii. For both envelopes, the cavity contribution in water is proportional to the surface area of the envelope. The same does not hold for GdmCl(aq), revealing the limitation of using molecular area to assess solvation energetics. The ΣG-cavity contribution predicts that GdmCl(aq) should favor the more compact state, contrary to the role of GdmCl in unfolding proteins. The chemistry contribution attenuates this effect, but still the net local (chemistry plus ΣG-packing) contribution is inadequate in capturing the role of GdmCl. With the inclusion of the long-range contribution, which is dominated by van der Waals interaction, aqueous GdmCl favors the extended conformation over the compact conformation. Our finding emphasizes the importance of weak, but attractive, long-range dispersion interactions in protein solution thermodynamics.

NMR spin-rotation relaxation and diffusion of methane.

The translational diffusion-coefficient and the spin-rotation contribution to the 1H NMR relaxation rate for methane (CH4) are investigated using MD (molecular dynamics) simulations, over a wide range of densities and temperatures, spanning the liquid, supercritical, and gas phases. The simulated diffusion-coefficients agree well with measurements, without any adjustable parameters in the interpretation of the simulations. A minimization technique is developed to compute the angular velocity for non-rigid spherical molecules, which is used to simulate the autocorrelation function for spin-rotation interactions. With increasing diffusivity, the autocorrelation function shows increasing deviations from the single-exponential decay predicted by the Langevin theory for rigid spheres, and the deviations are quantified using inverse Laplace transforms. The 1H spin-rotation relaxation rate derived from the autocorrelation function using the "kinetic model" agrees well with measurements in the supercritical/gas phase, while the relaxation rate derived using the "diffusion model" agrees well with measurements in the liquid phase. 1H spin-rotation relaxation is shown to dominate over the MD-simulated 1H-1H dipole-dipole relaxation at high diffusivity, while the opposite is found at low diffusivity. At high diffusivity, the simulated spin-rotation correlation time agrees with the kinetic collision time for gases, which is used to derive a new expression for 1H spin-rotation relaxation, without any adjustable parameters.

Role of internal motions and molecular geometry on the NMR relaxation of hydrocarbons.

The role of internal motions and molecular geometry on 1H NMR relaxation rates in liquid-state hydrocarbons is investigated using MD (molecular dynamics) simulations of the autocorrelation functions for intramolecular and intermolecular 1H-1H dipole-dipole interactions. The effects of molecular geometry and internal motions on the functional form of the autocorrelation functions are studied by comparing symmetric molecules such as neopentane and benzene to corresponding straight-chain alkanes n-pentane and n-hexane, respectively. Comparison of rigid versus flexible molecules shows that internal motions cause the intramolecular and intermolecular correlation-times to get significantly shorter, and the corresponding relaxation rates to get significantly smaller, especially for longer-chain n-alkanes. Site-by-site simulations of 1H's across the chains indicate significant variations in correlation times and relaxation rates across the molecule, and comparison with measurements reveals insights into cross-relaxation effects. Furthermore, the simulations reveal new insights into the relative strength of intramolecular versus intermolecular relaxation as a function of internal motions, as a function of molecular geometry, and on a site-by-site basis across the chain.

Quasichemical theory and the description of associating fluids relative to a reference: Multiple bonding of a single site solute.

We derive an expression for the chemical potential of an associating solute in a solvent relative to the value in a reference fluid using the quasichemical organization of the potential distribution theorem. The fraction of times the solute is not associated with the solvent, the monomer fraction, is expressed in terms of (a) the statistics of occupancy of the solvent around the solute in the reference fluid and (b) the Widom factors that arise because of turning on solute-solvent association. Assuming pair-additivity, we expand the Widom factor into a product of Mayer f-functions and the resulting expression is rearranged to reveal a form of the monomer fraction that is analogous to that used within the statistical associating fluid theory (SAFT). The present formulation avoids all graph-theoretic arguments and provides a fresh, more intuitive, perspective on Wertheim's theory and SAFT. Importantly, multi-body effects are transparently incorporated into the very foundations of the theory. We illustrate the generality of the present approach by considering examples of multiple solvent association to a colloid solute with bonding domains that range from a small patch on the sphere to a Janus particle to a solute whose entire surface is available for association.

Electrostatic and induction effects in the solubility of water in alkanes.

Experiments show that at 298 K and 1 atm pressure, the transfer free energy, µex, of water from its vapor to liquid normal alkanes CnH2n+2 (n=5…12) is negative. Earlier it was found that with the united-atom TraPPE model for alkanes and the SPC/E model for water, one had to artificially enhance the attractive alkane-water cross interaction to capture this behavior. Here we revisit the calculation of µex using the polarizable AMOEBA and the non-polarizable Charmm General (CGenFF) forcefields. We test both the AMOEBA03 and AMOEBA14 water models; the former has been validated with the AMOEBA alkane model while the latter is a revision of AMOEBA03 to better describe liquid water. We calculate µex using the test particle method. With CGenFF, µex is positive and the error relative to experiments is about 1.5 kBT. With AMOEBA, µex is negative and deviations relative to experiments are between 0.25 kBT (AMOEBA14) and 0.5 kBT (AMOEBA03). Quantum chemical calculations in a continuum solvent suggest that zero point effects may account for some of the deviation. Forcefield limitations notwithstanding, electrostatic and induction effects, commonly ignored in consideration of water-alkane interactions, appear to be decisive in the solubility of water in alkanes.

Thermodynamics of mixtures of patchy and spherical colloids of different sizes: A multi-body association theory with complete reference fluid information.

We present a theory to predict the structure and thermodynamics of mixtures of colloids of different diameters, building on our earlier work [A. Bansal et al., J. Chem. Phys. 145, 074904 (2016)] that considered mixtures with all particles constrained to have the same size. The patchy, solvent particles have short-range directional interactions, while the solute particles have short-range isotropic interactions. The hard-sphere mixture without any association site forms the reference fluid. An important ingredient within the multi-body association theory is the description of clustering of the reference solvent around the reference solute. Here we account for the physical, multi-body clusters of the reference solvent around the reference solute in terms of occupancy statistics in a defined observation volume. These occupancy probabilities are obtained from enhanced sampling simulations, but we also present statistical mechanical models to estimate these probabilities with limited simulation data. Relative to an approach that describes only up to three-body correlations in the reference, incorporating the complete reference information better predicts the bonding state and thermodynamics of the physical solute for a wide range of system conditions. Importantly, analysis of the residual chemical potential of the infinitely dilute solute from molecular simulation and theory shows that whereas the chemical potential is somewhat insensitive to the description of the structure of the reference fluid, the energetic and entropic contributions are not, with the results from the complete reference approach being in better agreement with particle simulations.

Structure and thermodynamics of a mixture of patchy and spherical colloids: A multi-body association theory with complete reference fluid information.

A mixture of solvent particles with short-range, directional interactions and solute particles with short-range, isotropic interactions that can bond multiple times is of fundamental interest in understanding liquids and colloidal mixtures. Because of multi-body correlations, predicting the structure and thermodynamics of such systems remains a challenge. Earlier Marshall and Chapman [J. Chem. Phys. 139, 104904 (2013)] developed a theory wherein association effects due to interactions multiply the partition function for clustering of particles in a reference hard-sphere system. The multi-body effects are incorporated in the clustering process, which in their work was obtained in the absence of the bulk medium. The bulk solvent effects were then modeled approximately within a second order perturbation approach. However, their approach is inadequate at high densities and for large association strengths. Based on the idea that the clustering of solvent in a defined coordination volume around the solute is related to occupancy statistics in that defined coordination volume, we develop an approach to incorporate the complete information about hard-sphere clustering in a bulk solvent at the density of interest. The occupancy probabilities are obtained from enhanced sampling simulations but we also develop a concise parametric form to model these probabilities using the quasichemical theory of solutions. We show that incorporating the complete reference information results in an approach that can predict the bonding state and thermodynamics of the colloidal solute for a wide range of system conditions.

Solvation free energy of the peptide group: its model dependence and implications for the additive-transfer free-energy model of protein stability.

The group-additive decomposition of the unfolding free energy of a protein in an osmolyte solution relative to that in water poses a fundamental paradox: whereas the decomposition describes the experimental results rather well, theory suggests that a group-additive decomposition of free energies is, in general, not valid. In a step toward resolving this paradox, here we study the peptide-group transfer free energy. We calculate the vacuum-to-solvent (solvation) free energies of (Gly)n and cyclic diglycine (cGG) and analyze the data according to experimental protocol. The solvation free energies of (Gly)n are linear in n, suggesting group additivity. However, the slope interpreted as the free energy of a peptide unit differs from that for cGG scaled by a factor of half, emphasizing the context dependence of solvation. However, the water-to-osmolyte transfer free energies of the peptide unit are relatively independent of the peptide model, as observed experimentally. To understand these observations, a way to assess the contribution to the solvation free energy of solvent-mediated correlation between distinct groups is developed. We show that linearity of solvation free energy with n is a consequence of uniformity of the correlation contributions, with apparent group-additive behavior in the water-to-osmolyte transfer arising due to their cancellation. Implications for inferring molecular mechanisms of solvent effects on protein stability on the basis of the group-additive transfer model are suggested.

Separating the role of protein restraints and local metal-site interaction chemistry in the thermodynamics of a zinc finger protein.

We express the effective Hamiltonian of an ion-binding site in a protein as a combination of the Hamiltonian of the ion-bound site in vacuum and the restraints of the protein on the site. The protein restraints are described by the quadratic elastic network model. The Hamiltonian of the ion-bound site in vacuum is approximated as a generalized Hessian around the minimum energy configuration. The resultant of the two quadratic Hamiltonians is cast into a pure quadratic form. In the canonical ensemble, the quadratic nature of the resultant Hamiltonian allows us to express analytically the excess free energy, enthalpy, and entropy of ion binding to the protein. The analytical expressions allow us to separate the roles of the dynamic restraints imposed by the protein on the binding site and the temperature-independent chemical effects in metal-ligand coordination. For the consensus zinc-finger peptide, relative to the aqueous phase, the calculated free energy of exchanging Zn(2+) with Fe(2+), Co(2+), Ni(2+), and Cd(2+) are in agreement with experiments. The predicted excess enthalpy of ion exchange between Zn(2+) and Co(2+) also agrees with the available experimental estimate. The free energy of applying the protein restraints reveals that relative to Zn(2+), the Co(2+), and Cd(2+)-site clusters are more destabilized by the protein restraints. This leads to an experimentally testable hypothesis that a tetrahedral metal binding site with minimal protein restraints will be less selective for Zn(2+) over Co(2+) and Cd(2+) compared to a zinc finger peptide. No appreciable change is expected for Fe(2+) and Ni(2+). The framework presented here may prove useful in protein engineering to tune metal selectivity.

The role of bulk protein in local models of ion-binding to proteins: comparative study of KcsA, its semisynthetic analog with a locked-in binding site, and valinomycin.

In studying ion-selectivity in biomaterials, it is common to study ion-protein interactions within a local neighborhood around the ion. This local system analysis for the S(2) site of KcsA, its semisynthetic analog, and valinomycin yields the free energy change in exchanging K(+) with Na(+) in quantitative agreement with the value obtained by considering ion-interactions with the entire system. But the energetics of ion binding in the local system and in the entire system differ significantly and lead to different conclusions regarding the physical basis of ion selectivity. For configurations sampled from an all-atom simulation, we show that the selectivity free energy can be decomposed into a contribution arising from interactions of the ion with its local neighborhood, ΔW(local), and a term arising from the field imposed on the ion and the binding site by the rest of the medium, ΔW(ϕ). The local contribution ΔW(local) is numerically close to the actual free energy difference because the field contribution is small. The field contribution is small because of cancellation of inversely related ion-medium and site-medium interactions. Our analysis presents a rigorous foundation for the numerical success of the local system analysis and shows that its implications do not always hold for the entire protein.

Communication: regularizing binding energy distributions and thermodynamics of hydration: theory and application to water modeled with classical and ab initio simulations.

The high-energy tail of the distribution of solute-solvent interaction energies is poorly characterized for condensed systems, but this tail region is of principal interest in determining the excess free energy of the solute. We introduce external fields centered on the solute to modulate the short-range repulsive interaction between the solute and solvent. This regularizes the binding energy distribution and makes it easy to calculate the free energy of the solute with the field. Together with the work done to apply the field in the presence and absence of the solute, we calculate the excess chemical potential of the solute. We present the formal development of this idea and apply it to study liquid water.

Water coordination structures and the excess free energy of the liquid.

We assess the contribution of each coordination state to the hydration free energy of a distinguished water molecule, the solute water. We define a coordination sphere, the inner-shell, and separate the hydration free energy into packing, outer-shell, and local, solute-specific (chemical) contributions. The coordination state is defined by the number of solvent water molecules within the coordination sphere. The packing term accounts for the free energy of creating a solute-free coordination sphere in the liquid. The outer-shell contribution accounts for the interaction of the solute with the fluid outside the coordination sphere and it is accurately described by a Gaussian model of hydration for coordination radii greater than the minimum of the oxygen-oxygen pair-correlation function: theory helps identify the length scale to parse chemical contributions from bulk, nonspecific contributions. The chemical contribution is recast as a sum over coordination states. The nth term in this sum is given by the probability p(n) of observing n water molecules inside the coordination sphere in the absence of the solute water times a factor accounting for the free energy, W(n), of forming an n-water cluster around the solute. The p(n) factors thus reflect the intrinsic properties of the solvent while W(n) accounts for the interaction between the solute and inner-shell solvent ligands. We monitor the chemical contribution to the hydration free energy by progressively adding solvent ligands to the inner-shell and find that four-water molecules are needed to fully account for the chemical term. For a chemically meaningful coordination radius, we find that W(4) ≈ W(1) and thus the interaction contribution is principally accounted for by the free energy for forming a one-water cluster, and intrinsic occupancy factors alone account for over half of the chemical contribution. Our study emphasizes the need to acknowledge the intrinsic solvent properties in interpreting the hydration structure of any solute, with particular care in cases where the solute-solvent interaction strength is similar to that between the solvent molecules.

Ion-water clusters, bulk medium effects, and ion hydration.

Thermochemistry of gas-phase ion-water clusters together with estimates of the hydration free energy of the clusters and the water ligands are used to calculate the hydration free energy of the ion. Often the hydration calculations use a continuum model of the solvent. The primitive quasichemical approximation to the quasichemical theory provides a transparent framework to anchor such efforts. Here we evaluate the approximations inherent in the primitive quasichemical approach and elucidate the different roles of the bulk medium. We find that the bulk medium can stabilize configurations of the cluster that are usually not observed in the gas phase, while also simultaneously lowering the excess chemical potential of the ion. This effect is more pronounced for soft ions. Since the coordination number that minimizes the excess chemical potential of the ion is identified as the optimal or most probable coordination number, for such soft ions the optimum cluster size and the hydration thermodynamics obtained with and without account of the bulk medium on the ion-water clustering reaction can be different. The ideas presented in this work are expected to be relevant to experimental studies that translate thermochemistry of ion-water clusters to the thermodynamics of the hydrated ion and to evolving theoretical approaches that combine high-level calculations on clusters with coarse-grained models of the medium.

Ion selectivity from local configurations of ligands in solutions and ion channels.

Probabilities of numbers of ligands proximal to an ion lead to simple, general formulae for the free energy of ion selectivity between different media. That free energy does not depend on the definition of an inner shell for ligand-counting, but other quantities of mechanistic interest do. If analysis is restricted to a specific coordination number, then two distinct probabilities are required to obtain the free energy in addition. The normalizations of those distributions produce partition function formulae for the free energy. Quasi-chemical theory introduces concepts of chemical equilibrium, then seeks the probability that is simplest to estimate, that of the most probable coordination number. Quasi-chemical theory establishes the utility of distributions of ligand-number, and sharpens our understanding of quasi-chemical calculations based on electronic structure methods. This development identifies contributions with clear physical interpretations, and shows that evaluation of those contributions can establish a mechanistic understanding of the selectivity in ion channels.

Communication: Thermodynamics of water modeled using ab initio simulations.

We regularize the potential distribution framework to calculate the excess free energy of liquid water simulated with the BLYP-D density functional. Assuming classical statistical mechanical simulations at 350 K model the liquid at 298 K, the calculated free energy is found in fair agreement with experiments, but the excess internal energy and hence also the excess entropy are not. The utility of thermodynamic characterization in understanding the role of high temperatures to mimic nuclear quantum effects and in evaluating ab initio simulations is noted.

Molecular packing and chemical association in liquid water simulated using ab initio hybrid Monte Carlo and different exchange-correlation functionals.

In the free energy of hydration of a solute, the chemical contribution is given by the free energy required to expel water molecules from the coordination sphere and the packing contribution is given by the free energy required to create the solute-free coordination sphere (the observation volume) in bulk water. With the simple point charge/extended (SPC/E) water model as a reference, we examine the chemical and packing contributions in the free energy of water simulated using different electron density functionals. The density is fixed at a value corresponding to that for SPC/E water at a pressure of 1 bar. The chemical contribution shows that water simulated at 300 K with BLYP is somewhat more tightly bound than water simulated at 300 K with revised PBE (revPBE) functional or at 350 K with the BLYP and BLYP-D functionals. The packing contribution for various radii of the observation volume is studied. In the size range where the distribution of water molecules in the observation volume is expected to be Gaussian, the packing contribution is expected to scale with the volume of the observation sphere. Water simulated at 300 K with the revPBE and at 350 K with BLYP-D or BLYP conforms to this expectation, but the results suggest an earlier onset of system size effects in the BLYP 350 K and revPBE 300 K systems than that observed for either BLYP-D 350 K or SPC/E. The implication of this observation for constant pressure simulations is indicated. For water simulated at 300 K with BLYP, in the size range where Gaussian distribution of occupation is expected, we instead find non-Gaussian behavior, and the packing contribution scales with surface area of the observation volume, suggesting the presence of heterogeneities in the system.

NMR 1H-1H Dipole Relaxation in Fluids: Relaxation of Individual 1H-1H Pairs versus Relaxation of Molecular Modes.

The intramolecular 1H NMR dipole-dipole relaxation of molecular fluids has traditionally been interpreted within the Bloembergen-Purcell-Pound (BPP) theory of NMR intramolecular relaxation. The BPP theory draws upon Debye's theory for describing the rotational diffusion of the 1H-1H pair and predicts a monoexponential decay of the 1H-1H dipole-dipole autocorrelation function between distinct spin pairs. Using molecular dynamics (MD) simulations, we show that for both n-heptane and water this is not the case. In particular, the autocorrelation function of individual 1H-1H intramolecular pairs itself evinces a rich stretched-exponential behavior, implying a distribution in rotational correlation times. However, for the high-symmetry molecule neopentane, the individual 1H-1H intramolecular pairs do conform to the BPP description, suggesting an important role of molecular symmetry in aiding agreement with the BPP model. The intermolecular autocorrelation functions for n-heptane, water, and neopentane also do not admit a monoexponential behavior of individual 1H-1H intermolecular pairs at distinct initial separations. We suggest expanding the autocorrelation function in terms of modes, provisionally termed molecular modes, that do have an exponential relaxation behavior. With care, the resulting Fredholm integral equation of the first kind can be inverted to recover the probability distribution of the molecular modes. The advantages and limitations of this approach are noted.