Search for: All records

We report a Monte Carlo simulation study of length-scale dependent density fluctuations in cavities in the coarse-grained mW representation of water at ambient conditions. Specifically, we use a combination of test particle insertion and umbrella sampling techniques to examine the full range of water occupation states in spherical cavities up to 6.3 Å in radius in water. As has previously been observed, water density fluctuations are found to be effectively Gaussian in nature for atomic-scale cavities, but as the cavities get larger they exhibit a non-Gaussian “fat-tail” distribution for lower occupancy states. We introduce a new statistical thermodynamic approach to analyze non-Gaussian fluctuations based on the radial distribution of waters about cavities with varying numbers of waters within its boundaries. It is shown that the onset on these non-Gaussian fluctuations is a result of the formation of a bubble within the cavity as it is emptied that is accompanied by the adsorption of waters onto its interior surface. We revisit a theoretical framework we previously introduced to describe Gaussian fluctuations within cavities to now incorporate bubble formation by including surface tension contributions. This modified theory accurately describes density fluctuations within both atomic and meso-scale cavities. Moreover, the theory predicts the transition from Gaussian to non-Gaussian fluctuations at a specific cavity occupancy in excellent agreement with simulation observations. 

Hydrophobic interactions drive the binding of nonpolar ligands to the oily pockets of proteins and supramolecular species in aqueous solution. As such, the wetting of host pockets is expected to play a critical role in determining the thermodynamics of guest binding. Here we use molecular simulations to examine the impact of pressure on the wetting and dewetting of the nonpolar pockets of a series of deep-cavity cavitands in water. The portals to the cavitand pockets are functionalized with both nonpolar (methyl) and polar (hydroxyl) groups oriented pointing either upward or inward toward the pocket. We find wetting of the pocket is favored by the hydroxyl groups and dewetting is favored by the methyl groups. The distribution of waters in the pocket is found to exhibit a two-state-like equilibrium between wet and dry states with a free energy barrier between the two states. Moreover, we demonstrate that the pocket hydration of the cavitands can be collapsed onto a unified adsorption isotherm by assuming the effective pressures within each cavitand pocket differ by a shift pressure that depends on the chemical identity and number of functional groups placed about the portal. These observations support the development of a twostate capillary evaporation model that accurately describes the equilibrium between states and naturally gives rise to the effective shift pressures observed from simulation. This work demonstrates that the hydration of host pockets can be tuned following simple design rules that in turn are expected to impact the thermodynamics of guest complexation. 

There are many open questions regarding the hydration of solvent-exposed, non-polar tracts and pockets in proteins. Although water is predicted to de-wet purely repulsive surfaces and evacuate crevices, the extent of de-wetting is unclear when ubiquitous van der Waals interactions are in play. The structural simplicity of synthetic supramolecular hosts imbues them with considerable potential to address this issue. To this end, here we detail a combination of densimetry and molecular dynamics simulations of three cavitands, coupled with calorimetric studies of their complexes with short-chain carboxylates. Our results reveal the range of wettability possible within the ostensibly identical cavitand pockets — which differ only in the presence/position of the methyl groups that encircle the portal to their non-polar pockets. The results demonstrate the ability of macrocycles to template water cavitation within their binding sites and show how the orientation of methyl groups can trigger the drying of non-polar pockets in liquid water, suggesting new avenues to control guest complexation. 

Science still does not have the ability to accurately predict the affinity that ligands have for proteins. In an attempt to address this, the Statistical Assessment of Modeling of Proteins and Ligands (SAMPL) series of blind predictive challenges is a community-wide exercise aimed at advancing computational techniques as standard predictive tools in rational drug design. In each cycle, a range of biologically relevant systems of different levels of complexity are selected to test the latest modeling methods. As part of this on-going exercise, and as a step towards understanding the important factors in context dependent guest binding, we challenged the computational community to determine the affinity of a series of negatively and positively charged guests to two constitutionally isomeric cavitand hosts: octa-acid 1 , and exo -octa acid 2 . Our affinity determinations, combined with molecular dynamics simulations, reveal asymmetries in affinities between host–guest pairs that cannot alone be explained by simple coulombic interactions, but also point to the importance of host–water interactions. Our work reveals the key facets of molecular recognition in water, emphasizes where improvements need to be made in modelling, and shed light on the complex problem of ligand-protein binding in the aqueous realm.