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.
Quantifying the Molecular Polarization Response of Liquid Water Interfaces at Heterogeneously Charged Surfaces
The hydration shells of proteins mediate interactions, such as small molecule binding, that are vital to their biological function or in some cases their dysfunction. However, even when the structure of a protein is known, the properties of its hydration environment cannot be easily predicted due to the complex interplay between protein surface heterogeneity and the collective structure of water’s hydrogen bonding network. This manuscript presents a theoretical study of the influence of surface charge heterogeneity on the polarization response of the liquid water interface. We focus our attention on classical point charge models of water, where the polarization response is limited to molecular reorientation. We introduce a new computational method for analyzing simulation data that is capable of quantifying water’s collective polarization response and determining the effective surface charge distribution of hydrated surfaces over atomistic length scales. To illustrate the utility of this method, we present the results of molecular dynamics simulations of liquid water in contact with a heterogeneous model surface and the CheY protein.
more »
« less
- Award ID(s):
- 1654415
- NSF-PAR ID:
- 10402744
- Date Published:
- Journal Name:
- Journal of Chemical Theory and Computation
- ISSN:
- 1549-9618
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
Dr. Sudipta Maiti (Ed.)IA3 is a 68 amino acid natural peptide/protein inhibitor of yeast aspartic proteinase A (YPRA) that is intrinsically dis-ordered in solution with induced N-terminal helicity when in the protein complex with YPRA. Based upon the intrinsical-ly disordered proteins (IDPs) parameters of fractional net charge (FNC), of net charge density per residue (NCPR) and of charge patterning (), the two domains of IA3 are defined to occupy different domains within conformationally based subclasses of IDPs; thus, making IA3 a bimodal-domain IDP. Site-directed spin-labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy and low-field Overhauser dynamic nuclear polarization (ODNP) spectroscopy results show that these two domains possess different degrees of compaction and hydration diffusivity behavior. This work suggests that SDSL EPR line shapes – analyzed in terms of their local tumbling volume (VL) – provide insight into the compaction of the unstructured IDP ensemble in solution and that protein sequence and net charge distribution pat-terns within a conformational subclass can impact bound water hydration dynamics; thus, possibly offering an alter-native thermodynamic property that can encode conforma-tional binding and behavior of IDPs and liquid-liquid phase separations.more » « less
-
na (Ed.)The surface hydration diffusivity of Bacillus subtilis Lipase A (BSLA) has been characterized by low-field Overhauser dynamic nuclear polarization (ODNP) relaxometry using a series of spin-labeled constructs. Sites for spin-label incorporation were previously designed via an atomistic computational approach that screened for surface exposure, reflective of the surface hydration comparable to other proteins studied by this method, as well as minimal impact on protein function, dynamics, and structure of BSLA by excluding any surface site that participated in greater than 30% occupancy of a hydrogen bonding network within BSLA. Experimental ODNP relaxometry coupling factor results verify the overall surface hydration behavior for these BSLA spin-labeled sites similar to other globular proteins. Here, by plotting the ODNP parameters of relative diffusive water versus the relative bound water, we introduce an effective "phase-space" analysis, which provides a facile visual comparison of the ODNP parameters of various biomolecular systems studied to date. We find notable differences when comparing BSLA to other systems, as well as when comparing different clusters on the surface of BSLA. Specifically, we find a grouping of sites that correspond to the spin-label surface location within the two main hydrophobic core clusters of the branched aliphatic amino acids isoleucine, leucine, and valine cores observed in the BSLA crystal structure. The results imply that hydrophobic clustering may dictate local surface hydration properties, perhaps through modulation of protein conformations and samplings of the unfolded states, providing insights into how the dynamics of the hydration shell is coupled to protein motion and fluctuations.more » « less
-
DNA functions only in aqueous environments and adopts different conformations depending on the hydration level. The dynamics of hydration water and hydrated DNA leads to rotating and oscillating dipoles that, in turn, give rise to a strong megahertz to terahertz absorption. Investigating the impact of hydration on DNA dynamics and the spectral features of water molecules influenced by DNA, however, is extremely challenging because of the strong absorption of water in the megahertz to terahertz frequency range. In response, we have employed a high-precision megahertz to terahertz dielectric spectrometer, assisted by molecular dynamics simulations, to investigate the dynamics of water molecules within the hydration shells of DNA as well as the collective vibrational motions of hydrated DNA, which are vital to DNA conformation and functionality. Our results reveal that the dynamics of water molecules in a DNA solution is heterogeneous, exhibiting a hierarchy of four distinct relaxation times ranging from ~8 ps to 1 ns, and the hydration structure of a DNA chain can extend to as far as ~18 A from its surface. The low-frequency collective vibrational modes of hydrated DNA have been identified and found to be sensitive to environmental conditions including temperature and hydration level. The results reveal critical information on hydrated DNA dynamics and DNA-water interfaces, which impact the biochemical functions and reactivity of DNA.more » « less
-
Zwitterionic materials are an important class of antifouling biomaterials for various applications. Despite such desirable antifouling properties, molecular-level understanding of the structure–property relationship associated with surface chemistry/topology/hydration and antifouling performance still remains to be elucidated. In this work, we computationally studied the packing structure, surface hydration, and antifouling property of three zwitterionic polymer brushes of poly(carboxybetaine methacrylate) (pCBMA), poly(sulfobetaine methacrylate) (pSBMA), and poly((2-(methacryloyloxy)ethyl)phosporylcoline) (pMPC) brushes and a hydrophilic PEG brush using a combination of molecular mechanics (MM), Monte Carlo (MC), molecular dynamics (MD), and steered MD (SMD) simulations. We for the first time determined the optimal packing structures of all polymer brushes from a wide variety of unit cells and chain orientations in a complex energy landscape. Under the optimal packing structures, MD simulations were further conducted to study the structure, dynamics, and orientation of water molecules and protein adsorption on the four polymer brushes, while SMD simulations to study the surface resistance of the polymer brushes to a protein. The collective results consistently revealed that the three zwitterionic brushes exhibited stronger interactions with water molecules and higher surface resistance to a protein than the PEG brush. It was concluded that both the carbon space length between zwitterionic groups and the nature of the anionic groups have a distinct effect on the antifouling performance, leading to the following antifouling ranking of pCBMA > pMPC > pSBMA. This work hopefully provides some structural insights into the design of new antifouling materials beyond traditional PEG-based antifouling materials.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1021/acs.jctc.2c01256
