Patrick Varilly

Fluctuations of Water near Extended Hydrophobic and Hydrophilic Surfaces

We use molecular dynamics simulations of the SPC-E model of liquid water to derive probability distributions for water density fluctuations in probe volumes of different shapes and sizes, both in the bulk as well as near hydrophobic and hydrophilic surfaces. Our results are obtained with a biased sampling of coarse-grained densities that is easily combined with molecular dynamics integration algorithms. Our principal result is that the probability for density fluctuations of water near a hydrophobic surface, with or without surface water attractions, is akin to density fluctuations at the water-vapor interface. Specifically, the probability of density depletion near the surface is significantly larger than that in the bulk, and this enhanced probability is responsible for hydrophobic forces of assembly. In contrast, we find that the statistics of water density fluctuations near a model hydrophilic surface are similar to that in the bulk.

Sitting at the Edge: How Biomolecules use Hydrophobicity to Tune Their Interactions and Function

Water near extended hydrophobic surfaces is like that at a liquid-vapor interface, where fluctuations in water density are substantially enhanced compared to those in bulk water. Here we use molecular simulations with specialized sampling techniques to show that water density fluctuations are similarly enhanced, even near hydrophobic surfaces of complex biomolecules, situating them at the edge of a dewetting transition. Consequently, water near these surfaces is sensitive to subtle changes in surface conformation, topology, and chemistry, any of which can tip the balance toward or away from the wet state and thus significantly alter biomolecular interactions and function. Our work also resolves the long-standing puzzle of why some biological surfaces dewet and other seemingly similar surfaces do not.

Extended surfaces modulate hydrophobic interactions of neighboring solutes.

Amish J. Patel, Patrick Varilly, Sumanth N. Jamadagni, Hari Acharya, Shekhar Garde, ∗ and David Chandler ∗ Howard P. Isermann Department of Chemical & Biological Engineering, and Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute, Troy, NY 12180, USA Department of Chemistry, University of California, Berkeley, CA 94720, USA Abstract Interfaces are a most common motif in complex systems. To understand how the presence of interfaces affect hydrophobic phenomena, we use molecular simulations and theory to study hydration of solutes at interfaces. The solutes range in size from sub-nanometer to a few nanometers. The interfaces are self-assembled monolayers with a range of chemistries, from hydrophilic to hydrophobic. We show that the driving force for assembly in the vicinity of a hydrophobic surface is weaker than that in bulk water, and decreases with increasing temperature, in contrast to that in the bulk. We explain these distinct features in terms of an interplay between interfacial fluctuations and excluded volume effects—the physics encoded in Lum-Chandler-Weeks theory [J. Phys. Chem. B 103 4570–4577 (1999)]. Our results suggest a catalytic role for hydrophobic interfaces in the unfolding of proteins, for example, in the interior of chaperonins and in amyloid formation.Interfaces are a most common motif in complex systems. To understand how the presence of interfaces affects hydrophobic phenomena, we use molecular simulations and theory to study hydration of solutes at interfaces. The solutes range in size from subnanometer to a few nanometers. The interfaces are self-assembled monolayers with a range of chemistries, from hydrophilic to hydrophobic. We show that the driving force for assembly in the vicinity of a hydrophobic surface is weaker than that in bulk water and decreases with increasing temperature, in contrast to that in the bulk. We explain these distinct features in terms of an interplay between interfacial fluctuations and excluded volume effects—the physics encoded in Lum–Chandler–Weeks theory [Lum K, Chandler D, Weeks JD (1999) J Phys Chem B 103:4570–4577]. Our results suggest a catalytic role for hydrophobic interfaces in the unfolding of proteins, for example, in the interior of chaperonins and in amyloid formation.

An improved coarse-grained model of solvation and the hydrophobic effect

We present a coarse-grained lattice model of solvation thermodynamics and the hydrophobic effect that implements the ideas of Lum-Chandler-Weeks theory [J. Phys. Chem. B 134, 4570 (1999)] and improves upon previous lattice models based on it. Through comparison with molecular simulation, we show that our model captures the length-scale and curvature dependence of solvation free energies with near-quantitative accuracy and 2-3 orders of magnitude less computational effort, and further, correctly describes the large but rare solvent fluctuations that are involved in dewetting, vapor tube formation, and hydrophobic assembly. Our model is intermediate in detail and complexity between implicit-solvent models and explicit-water simulations.

A general theory of DNA-mediated and other valence-limited colloidal interactions

We present a general theory for predicting the interaction potentials between DNA-coated colloids, and more broadly, any particles that interact via valence-limited ligand-receptor binding. Our theory correctly incorporates the configurational and combinatorial entropic factors that play a key role in valence-limited interactions. By rigorously enforcing self-consistency, it achieves near-quantitative accuracy with respect to detailed Monte Carlo calculations. With suitable approximations and in particular geometries, our theory reduces to previous successful treatments, which are now united in a common and extensible framework. We expect our tools to be useful to other researchers investigating ligand-mediated interactions. A complete and well-documented Python implementation is freely available at http://github.com/patvarilly/DNACC.

Water evaporation: a transition path sampling study.

We use transition path sampling to study evaporation in the SPC/E model of liquid water. On the basis of thousands of evaporation trajectories, we characterize the members of the transition state ensemble (TSE), which exhibit a liquid-vapor interface with predominantly negative mean curvature at the site of evaporation. We also find that after evaporation is complete, the distributions of translational and angular momenta of the evaporated water are Maxwellian with a temperature equal to that of the liquid. To characterize the evaporation trajectories in their entirety, we find that it suffices to project them onto just two coordinates: the distance of the evaporating molecule to the instantaneous liquid-vapor interface and the velocity of the water along the average interface normal. In this projected space, we find that the TSE is well-captured by a simple model of ballistic escape from a deep potential well, with no additional barrier to evaporation beyond the cohesive strength of the liquid. Equivalently, they are consistent with a near-unity probability for a water molecule impinging upon a liquid droplet to condense. These results agree with previous simulations and with some, but not all, recent experiments.

Transcendental obstructions to weak approximation on general K3 surfaces

Abstract We construct an explicit K3 surface over the field of rational numbers that has geometric Picard rank one, and for which there is a transcendental Brauer–Manin obstruction to weak approximation. To do so, we exploit the relationship between polarized K3 surfaces endowed with particular kinds of Brauer classes and cubic fourfolds.

Effect of Inert Tails on the Thermodynamics of DNA Hybridization

The selective hybridization of DNA is of key importance for many practical applications such as gene detection and DNA-mediated self-assembly. These applications require a quantitative prediction of the hybridization free energy. Existing methods ignore the effects of non-complementary ssDNA tails beyond the first unpaired base. We use experiments and simulations to show that the binding strength of complementary ssDNA oligomers is altered by these sequences of non-complementary nucleotides. Even a small number of non-binding bases are enough to raise the hybridization free energy by approximately 1 kcal/mol at physiological salt concentrations. We propose a simple analytical expression that accounts quantitatively for this variation as a function of tail length and salt concentration.

Predicting DNA-mediated colloidal pair interactions

Recently, Rogers and Crocker (1) proposed a method to predict the interaction between colloids coated with two kinds of ssDNA, A and B. A key step in ref. 1 was to estimate the average number of DNA bonds, 〈N〉, assuming local chemical equilibrium (LCE) between hybridized and unhybridized sticky end concentrations:

Consistent Treatment of Hydrophobicity in Protein Lattice Models Accounts for Cold Denaturation.

The hydrophobic effect stabilizes the native structure of proteins by minimizing the unfavorable interactions between hydrophobic residues and water through the formation of a hydrophobic core. Here, we include the entropic and enthalpic contributions of the hydrophobic effect explicitly in an implicit solvent model. This allows us to capture two important effects: a length-scale dependence and a temperature dependence for the solvation of a hydrophobic particle. This consistent treatment of the hydrophobic effect explains cold denaturation and heat capacity measurements of solvated proteins.