Classical molecular dynamics simulations of the hydration thermodynamics, structure, and dynamics of water in hydration shells of charged buckminsterfullerenes are presented in this study. Charging of fullerenes leads to a structural transition in the hydration shell, accompanied by creation of a significant population of dangling O–H bonds pointing toward the solute. In contrast to the well accepted structure–function paradigm, this interfacial structural transition causes nearly no effect on either the dynamics of hydration water or on the solvation thermodynamics. Linear response to the solute charge is maintained despite significant structural changes in the hydration shell, and solvation thermodynamic potentials are nearly insensitive to the altering structure. Only solvation heat capacities, which are higher thermodynamic derivatives of the solvation free energy, indicate some sensitivity to the local hydration structure. We have separated the solvation thermodynamic potentials into direct solute–solvent interactions and restructuring of the hydration shell and analyzed the relative contributions of electrostatic and nonpolar interactions to the solvation thermodynamics.
more »
« less
Entropy connects water structure and dynamics in protein hydration layer
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
- Award ID(s):
- 1665157
- PAR ID:
- 10062698
- Date Published:
- Journal Name:
- Physical Chemistry Chemical Physics
- Volume:
- 20
- Issue:
- 21
- ISSN:
- 1463-9076
- Page Range / eLocation ID:
- 14765 to 14777
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Proteins and water couple dynamically over a wide range of time scales. Motivated by their central role in protein function, protein-water dynamics and thermodynamics have been extensively studied for structured proteins, where correspondence to structural features has been made. However, properties controlling intrinsically disordered protein (IDP)-water dynamics are not yet known. We report results of megahertz-to-terahertz dielectric spectroscopy and molecular dynamics simulations of a group of IDPs with varying charge content along with structured proteins of similar size. Hydration water around IDPs is found to exhibit more heterogeneous rotational and translational dynamics compared with water around structured proteins of similar size, yielding on average more restricted dynamics around individual residues of IDPs, charged or neutral, compared with structured proteins. The on-average slower water dynamics is found to arise from excess tightly bound water in the first hydration layer, which is related to greater exposure to charged groups. The more tightly bound water to IDPs correlates with the smaller hydration shell found experimentally, and affects entropy associated with protein-water interactions, the contribution of which we estimate based on the dielectric measurements and simulations. Water-IDP dynamic coupling at terahertz frequencies is characterized by the dielectric measurements and simulations.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
-
null (Ed.)Biological processes often take place at surfaces of proteins, where the dynamic and structural properties of aqueous solvents are modified. Information about solvent properties including hydration dynamics and structure, and protein collective motions can be obtained by measuring directly the dielectric response in the megahertz to terahertz frequencies of aqueous protein solutions. Due to the strong absorption of water in this frequency range, the experiment is challenging. Our home built dielectric spectrometer using a vector network analyzer together with frequency extenders allows us to perform the experiment in a wide range of frequency from megahertz to terahertz with a high dynamical range up to 120 dB. A detailed investigation of the dielectric response has revealed the hydration structure including the tightly, loosely bound layers and the number of water molecules in each hydration layer. These water molecules relax with different time constants at different temperatures. As a result, the dynamics of hydrated protein is also probed at different temperatures. Understanding the hydration structure and dynamics of lysozyme in biological conditions can explain the enzymatic activities of biomolecules.more » « less
An official website of the United States government

