Rahul Godawat

Characterizing hydrophobicity of interfaces by using cavity formation, solute binding, and water correlations

Hydrophobicity is often characterized macroscopically by the droplet contact angle. Molecular signatures of hydrophobicity have, however, remained elusive. Successful theories predict a drying transition leading to a vapor-like region near large hard-sphere solutes and interfaces. Adding attractions wets the interface with local density increasing with attractions. Here we present extensive molecular simulation studies of hydration of realistic surfaces with a wide range of chemistries from hydrophobic (−CF3, −CH3) to hydrophilic (−OH, −CONH2). We show that the water density near weakly attractive hydrophobic surfaces (e.g., −CF3) can be bulk-like or larger, and provides a poor quantification of surface hydrophobicity. In contrast, the probability of cavity formation or the free energy of binding of hydrophobic solutes to interfaces correlates quantitatively with the macroscopic wetting properties and serves as an excellent signature of hydrophobicity. Specifically, the probability of cavity formation is enhanced in the vicinity of hydrophobic surfaces, and water–water correlations correspondingly display characteristics similar to those near a vapor–liquid interface. Hydrophilic surfaces suppress cavity formation and reduce the water–water correlation length. Our results suggest a potentially robust approach for characterizing hydrophobicity of more complex and heterogeneous surfaces of proteins and biomolecules, and other nanoscopic objects.

Hydrophobicity of proteins and interfaces: insights from density fluctuations.

Macroscopic characterizations of hydrophobicity (e.g., contact angle measurements) do not extend to the surfaces of proteins and nanoparticles. Molecular measures of hydrophobicity of such surfaces need to account for the behavior of hydration water. Theory and state-of-the-art simulations suggest that water density fluctuations provide such a measure; fluctuations are enhanced near hydrophobic surfaces and quenched with increasing surface hydrophilicity. Fluctuations affect conformational equilibria and dynamics of molecules at interfaces. Enhanced fluctuations are reflected in enhanced cavity formation, more favorable binding of hydrophobic solutes, increased compressibility of hydration water, and enhanced water-water correlations at hydrophobic surfaces. These density fluctuation-based measures can be used to develop practical methods to map the hydrophobicity/philicity of heterogeneous surfaces including those of proteins. They highlight that the hydrophobicity of a group is context dependent and is significantly affected by its environment (e.g., chemistry and topography) and especially by confinement. The ability to include information about hydration water in mapping hydrophobicity is expected to significantly impact our understanding of protein-protein interactions as well as improve drug design and discovery methods and bioseparation processes.

Unfolding of Hydrophobic Polymers in Guanidinium Chloride Solutions

Guanidinium chloride (GdmCl) is a widely used chemical denaturant that unfolds proteins. Its effects on hydrophobic interactions are, however, not fully understood. We quantify the effects of GdmCl on various manifestations of hydrophobicity--from solvation and interactions of small solutes to folding-unfolding of hydrophobic polymers--in water and in concentrated GdmCl solutions. For comparison, we also perform similar calculations in solutions of NaCl and CsCl in water. Like NaCl and CsCl, GdmCl increases the surface tension of water, decreases the solubility of small hydrophobic solutes, and enhances the strength of hydrophobic interactions at the pair level. However, unlike NaCl and CsCl, GdmCl destabilizes folded states of hydrophobic polymers. We show that Gdm(+) ions preferentially coat the hydrophobic polymer, and it is the direct van der Waals interaction between Gdm(+) ions and the polymer that contributes to the destabilization of folded states. Interestingly, the temperature dependence of the free energy of unfolding of the hydrophobic polymer in water is protein-like, with signatures of both heat and cold denaturation. Addition of GdmCl shifts the cold denaturation temperature higher, into the experimentally accessible region. Finally, translational as well as conformational dynamics of the polymer are slower in GdmCl and correlate with dynamics of water molecules in solution.

How Surface Wettability Affects the Binding, Folding, and Dynamics of Hydrophobic Polymers at Interfaces

We present an extensive molecular simulation study of the behavior of a flexible hydrophobic 25-mer polymer at interfaces presenting a range of chemistries from hydrophobic (-CH(3)) to hydrophilic (-CONH(2)). We quantify the free energy of adsorption, conformational equilibria, and translational and conformational dynamics of the polymer at these diverse interfaces. Water-mediated interactions drive the polymer to adsorb strongly at a hydrophobic interface and repel it from hydrophilic ones. At hydrophilic surfaces, van der Waals interactions between the polymer and the surface mitigate this water-mediated repulsion, leading to weak adsorption of the polymer. Although the polymer is strongly adsorbed to hydrophobic surfaces, it is also most dynamic there. Translational diffusion and conformational dynamics are faster at hydrophobic surfaces compared to those at hydrophilic ones. In bulk water, the polymer collapses into compact globular shapes, whereas the thermodynamic stability of folded polymers is significantly lowered at hydrophobic surfaces. The polymer spreads into pancake-like 2D conformations at hydrophobic surfaces and gradually beads up into globular shapes as the surface is made more hydrophilic. Interestingly, the binding thermodynamics and dynamics correlate with macroscopic droplet contact angles that characterize the wetting properties of the different interfaces.

Mechanistic studies of displacer–protein binding in chemically selective displacement systems using NMR and MD simulations

A parallel batch screening technique was employed to identify chemically selective displacers which exhibited exclusive separation behavior for the protein pair α‐chymotrypsin/ribonuclease A on a strong cation exchange resin. Two selective displacers, 1‐(4‐chlorobenzyl)piperidin‐3‐aminesulfate and N′1′‐(4‐methyl‐quinolin‐2‐yl)‐ethane‐1,2‐diamine dinitrate, and one non‐selective displacer, spermidine, were selected as model systems to investigate the mechanism of chemically selective displacement chromatography. Saturation transfer difference (STD) NMR was used to directly evaluate displacer–protein binding. The results indicated that while binding occurred between the two chemically selective displacers and the more hydrophobic protein, α‐chymotrypsin, no binding was observed with ribonuclease A. Further, the non‐selective displacer, spermidine, was not observed to bind to either protein. Importantly, the binding event was observed to occur primarily on the aromatic portion of the selective displacers. Extensive molecular dynamic simulations of protein–displacer–water solution were also carried out. The MD results corroborated the NMR findings demonstrating that the binding of selective displacers occurred primarily on hydrophobic surface patches of α‐chymotrypsin, while no significant long term binding to ribonuclease A was observed. The non‐selective displacer did not show significant binding to either of the proteins. MD simulations also indicated that the charged amine group of the selective displacers in the bound state was primarily oriented towards the solvent, potentially facilitating their interaction with a resin surface. These results directly confirm that selective binding between a protein and displacer is the mechanism by which chemically selective displacement occurs. This opens up many possibilities for future molecular design of selective displacers for a range of applications. Biotechnol. Bioeng. 2009;102: 1428–1437.