Interactions between proteins lie at the heart of numerous biological processes and are essential for the proper functioning of the cell. Although the importance of hydrophobic residues in driving protein interactions is universally accepted, a characterization of protein hydrophobicity, which informs its interactions, has remained elusive. The challenge lies in capturing the collective response of the protein hydration waters to the nanoscale chemical and topographical protein patterns, which determine protein hydrophobicity. To address this challenge, here, we employ specialized molecular simulations wherein water molecules are systematically displaced from the protein hydration shell; by identifying protein regions that relinquish their waters more readily than others, we are then able to uncover the most hydrophobic protein patches. Surprisingly, such patches contain a large fraction of polar/charged atoms and have chemical compositions that are similar to the more hydrophilic protein patches. Importantly, we also find a striking correspondence between the most hydrophobic protein patches and regions that mediate protein interactions. Our work thus establishes a computational framework for characterizing the emergent hydrophobicity of amphiphilic solutes, such as proteins, which display nanoscale heterogeneity, and for uncovering their interaction interfaces.
more »
« less
Water dynamics around T vs R4 of hemoglobin from local hydrophobicity analysis
The local hydration around tetrameric hemoglobin (Hb) in its T0 and R4 conformational substates is analyzed based on molecular dynamics simulations. Analysis of the local hydrophobicity (LH) for all residues at the α1β2 and α2β1 interfaces, responsible for the quaternary T → R transition, which is encoded in the Monod–Wyman–Changeux model, as well as comparison with earlier computations of the solvent accessible surface area, makes clear that the two quantities measure different aspects of hydration. Local hydrophobicity quantifies the presence and structure of water molecules at the interface, whereas “buried surface” reports on the available space for solvent. For simulations with Hb frozen in its T0 and R4 states, the correlation coefficient between LH and buried surface is 0.36 and 0.44, respectively, but it increases considerably if the 95% confidence interval is used. The LH with Hb frozen and flexible changes little for most residues at the interfaces but is significantly altered for a few select ones: Thr41α, Tyr42α, Tyr140α, Trp37β, Glu101β (for T0) and Thr38α, Tyr42α, Tyr140α (for R4). The number of water molecules at the interface is found to increase by ∼25% for T0 → R4, which is consistent with earlier measurements. Since hydration is found to be essential to protein function, it is clear that hydration also plays an essential role in allostery.
more »
« less
- Award ID(s):
- 1654415
- PAR ID:
- 10487345
- Publisher / Repository:
- Journal of Chemical Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 2
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The hydrophobicity of an interface determines the magnitude of hydrophobic interactions that drive numerous biological and industrial processes. Chemically heterogeneous interfaces are abundant in these contexts; examples include the surfaces of proteins, functionalized nanomaterials, and polymeric materials. While the hydrophobicity of nonpolar solutes can be predicted and related to the structure of interfacial water molecules, predicting the hydrophobicity of chemically heterogeneous interfaces remains a challenge because of the complex, non-additive contributions to hydrophobicity that depend on the chemical identity and nanoscale spatial arrangements of polar and nonpolar groups. In this work, we utilize atomistic molecular dynamics simulations in conjunction with enhanced sampling and data-centric analysis techniques to quantitatively relate changes in interfacial water structure to the hydration free energy (a thermodynamically well-defined descriptor of hydrophobicity) of chemically heterogeneous interfaces. We analyze a large data set of 58 self-assembled monolayers (SAMs) composed of ligands with nonpolar and polar end groups of different chemical identity (amine, amide, and hydroxyl) in five mole fractions, two spatial patterns, and with scaled partial charges. We find that only five features of interfacial water structure are required to accurately predict hydration free energies. Examination of these features reveals mechanistic insights into the interfacial hydrogen bonding behaviors that distinguish different surface compositions and patterns. This analysis also identifies the probability of highly coordinated water structures as a unique signature of hydrophobicity. These insights provide a physical basis to understand the hydrophobicity of chemically heterogeneous interfaces and connect hydrophobicity to experimentally accessible perturbations of interfacial water structure.more » « less
-
Resonantly enhanced second harmonic generation (SHG) spectra of Coumarin 152 (C152) adsorbed at the water-silica interface show that C152 experiences a local dielectric environment slightly more polar than that of bulk water. This result stands in contrast to recently reported time-resolved fluorescence experiments and simulations that suggest a alkane-like permittivity for interfacial water at strongly associating, hydrophilic solid surfaces. Taken together, these results imply that while the static electric field across the aqueous-silica interface may be large, restricted water dynamics lead to apparent nonpolar solvation behavior similar to that experienced by solutes in confinement. Resonance-enhanced SHG spectra and time-resolved fluorescence of C152 adsorbed to aqueous-hydrophobic silica surfaces show that when water’s ability to hydrogen bond with the silica surface is eliminated, a solute’s interfacial solvation and corresponding ability to photoisomerize converge to an intermediate limit similar to that experienced in bulk acetone or methanol. While water structure and dynamics at solid-liquid interfaces have received considerable attention, results presented below show how strong solvent-substrate interactions can create conflicting pictures of solute reactivity across buried interfaces.more » « less
-
The enzyme Candida Antarctica lipase B (CALB) serves here as a model for understanding connections among hydration layer dynamics, solvation shell structure, and protein surface structure. The structure and dynamics of water molecules in the hydration layer were characterized for regions of the CALB surface, divided around each α-helix, β-sheet, and loop structure. Heterogeneous hydration dynamics were observed around the surface of the enzyme, in line with spectroscopic observations of other proteins. Regional differences in the structure of the biomolecular hydration layer were found to be concomitant with variations in dynamics. In particular, it was seen that regions of higher density exhibit faster water dynamics. This is analogous to the behavior of bulk water, where dynamics (diffusion coefficients) are connected to water structure (density and tetrahedrality) by excess (or pair) entropy, detailed in the Rosenfeld scaling relationship. Additionally, effects of protein surface topology and hydrophobicity on water structure and dynamics were evaluated using multiregression analysis, showing that topology has a somewhat larger effect on hydration layer structure–dynamics. Concave and hydrophobic protein surfaces favor a less dense and more tetrahedral solvation layer, akin to a more ice-like structure, with slower dynamics. Results show that pairwise entropies of local hydration layers, calculated from regional radial distribution functions, scale logarithmically with local hydration dynamics. Thus, the Rosenfeld relationship describes the heterogeneous structure–dynamics of the hydration layer around the enzyme CALB. These findings raise the question of whether this may be a general principle for understanding the structure–dynamics of biomolecular solvation.more » « less
-
Determination of the surface hydrophobicity or wettability of nanomaterials and nanoparticles (NPs) is often challenged by the heterogeneous properties of NPs that vary with particle size, shape, surface charge, aggregation states, and surface sorption or coating. This study first summarized inherent limitations of the water contact angle, octanol–water partition coefficient ( K ow ) and surface adsorption of probe molecules in probing nanomaterial hydrophobicity. Then, we demonstrated the principle of a scanning probe method based on atomic force microscopy (AFM) for the local surface hydrophobicity measurement. Specifically, we measured the adhesion forces between functionalized AFM tips and self-assembled monolayers (SAMs) to establish a linear relationship between the adhesion forces and water contact angles based on the continuum thermodynamic approach (CTA). This relationship was used to determine the local surface hydrophobicity of seven different NPs ( i.e. , TiO 2 , ZnO, SiO 2 , CuO, CeO 2 , α-Fe 2 O 3 , and Ag), which agreed well with bulk contact angles of these NPs. Some discrepancies were observed for Fe 2 O 3 , CeO 2 and SiO 2 NPs, probably because of surface hydration and roughness effects. Moreover, the solution pH and ionic strength had negligible effects on the adhesion forces between the AFM tip and MWCNTs or C 60 , indicating that the hydrophobicity of carbonaceous nanomaterials is not influenced by pH or ionic strength (IS). By contrast, natural organic matter (NOM) appreciably decreased the hydrophobicity of MWCNTs and C 60 due to surface coating of hydrophilic NOM. This scanning probe method has been proved to be reliable and robust toward the accurate measurement of the nanoscale hydrophobicity of individual NPs or nanomaterials in liquid environments.more » « less
An official website of the United States government

