Cooperative hydrophobic/hydrophilic interactions in the hydration of dimethyl ether

Hydration Free Energies of Polypeptides from Popular Implicit Solvent Models versus All-Atom Simulation Results Based on Molecular Quasichemical Theory

Calculating the hydration free energy of a macromolecule in all-atom simulations has long remained a challenge, necessitating the use of models wherein the effect of the solvent is captured without explicit account of solvent degrees of freedom. This situation has changed with developments in the molecular quasi-chemical theory (QCT)─an approach that enables calculation of the hydration free energy of macromolecules within all-atom simulations at the same resolution as is possible for small molecular solutes. The theory also provides a rigorous and physically transparent framework to conceptualize and model interactions in molecular solutions and thus provides a convenient framework to investigate the assumptions in implicit solvent models. In this study, we compare the results using molecular QCT versus predictions from EEF1, ABSINTH, and GB/SA implicit solvent models for polyglycine and polyalanine solutes covering a range of chain lengths and conformations. The hydration free energies or the differences in hydration free energies between conformers obtained from the implicit solvent models do not agree with explicit solvent results, with the deviations being largest for the group additive EEF1 and ABSINTH models. GB/SA does better in capturing the qualitative trends seen in explicit solvent results. Analysis founded on QCT reveals the critical importance of the cooperativity of hydration that is inherent in the hydrophilic and hydrophobic contributions to hydration─physics that is not well captured in additive models but somewhat better accounted for by means of a dielectric in the GB/SA approach.

Thermodynamics of Hydration from the Perspective of the Molecular Quasichemical Theory of Solutions

The quasichemical organization of the potential distribution theorem, molecular quasichemical theory (QCT), enables practical calculations and also provides a conceptual framework for molecular hydration phenomena. QCT can be viewed from multiple perspectives: (a) as a way to regularize an ill-conditioned statistical thermodynamic problem; (b) as an introduction of and emphasis on the neighborship characteristics of a solute of interest; or (c) as a way to include accurate electronic structure descriptions of near-neighbor interactions in defensible statistical thermodynamics by clearly defining neighborship clusters. The theory has been applied to solutes of a wide range of chemical complexity, ranging from ions that interact with water with both long-ranged and chemically intricate short-ranged interactions, to solutes that interact with water solely through traditional van der Waals interations, and including water itself. The solutes range in variety from monatomic ions to chemically heterogeneous macromolecules. A notable feature of QCT is that, in applying the theory to this range of solutes, the theory itself provides guidance on the necessary approximations and simplifications that can facilitate the calculations. In this Perspective, we develop these ideas and document them with examples that reveal the insights that can be extracted using the QCT formulation.

Importance of Hydrophilic Hydration and Intramolecular Interactions in the Thermodynamics of Helix-Coil Transition and Helix-Helix Assembly in a Deca-Alanine Peptide

For a model deca-alanine peptide the cavity (ideal hydrophobic) contribution to hydration favors the helix state over extended states and the paired helix bundle in the assembly of two helices. The energetic contributions of attractive protein-solvent interactions are separated into quasi-chemical components consisting of a short-range part arising from interactions with solvent in the first hydration shell and the remaining long-range part that is well described by a Gaussian. In the helix-coil transition, short-range attractive protein-solvent interactions outweigh hydrophobic hydration and favor the extended coil states. Analysis of enthalpic effects shows that it is the favorable hydration of the peptide backbone that favors the unfolded state. Protein intramolecular interactions favor the helix state and are decisive in favoring folding. In the pairing of two helices, the cavity contribution outweighs the short-range attractive protein-water interactions. However, long-range, protein-solvent attractive interactions can either enhance or reverse this trend depending on the mutual orientation of the helices. In helix-helix assembly, change in enthalpy arising from change in attractive protein-solvent interactions favors disassembly. In helix pairing as well, favorable protein intramolecular interactions are found to be as important as hydration effects. Overall, hydrophilic protein-solvent interactions and protein intramolecular interactions are found to play a significant role in the thermodynamics of folding and assembly in the system studied.

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.