Distinct Protein Hydration Water Species Defined by Spatially Resolved Spectra of Intermolecular Vibrations

FTIR spectroscopy and molecular level insight of diluted aqueous solutions of acetic acid

Aqueous solutions of acetic acid (AA) have been intensively explored for decades with a particular attention addressed to the hydrogen bond network generated by COOH group at different concentrations. In majority of studies conducted so far the envelope originated from νCO is decomposed into two bands assigned to differently hydrated monomers: the one presumably to AA···H2O, and another one to AA···(H2O)2. In order to examine if species other than the mentioned monomers produce this spectral signature, we performed computational and FTIR spectroscopic study of AA in aqueous solutions. Dilute solutions of deuterated acetic acid (CD3COOD) in D2O and in C2Cl4 as a reference were prepared (c0 = 0.001, 0.01 and 0.1 mol dm-3) as well as of deuterated sodium acetate (CD3COONa) in D2O. CD3COOD in 0.1 mol dm-3 solution in D2O displays a feature that separated in two signals with maxima at 1706 cm-1 and 1687 cm-1. A combined DFT and molecular dynamics study performed in this work showed the assignation of those spectral bands to be a more complex problem than previously thought, with syn-anti isomerism and hydration contributing to the experimentally observed broad νCO envelope.

Molecular Dynamics Simulation of the Influence of Temperature and Salt on the Dynamic Hydration Layer in a Model Polyzwitterionic Polymer PAEDAPS

We investigate the hydration of poly(3-[2-(acrylamido) ethyldimethylammonio] propanesulfonate) over a range of temperatures in pure water and with the inclusion of 0.1 mol/L NaCl using atomistic molecular dynamics simulation. Drawing on concepts drawn from the field of glass-forming liquids, we use the Debye-Waller parameter () for describing the water mobility gradient around the polybetaine backbone extending to an overall distance ≈18 Å. The water mobility in this layer is defined through the mean-square water molecule displacement at a time on the order of water's β-relaxation time. The brushlike topology of polybetaines leads to two regions in the dynamic hydration layer. The inner region of ≈10.5 Å is explored by pendant group conformational motions, and the outer region of ≈7.5 Å represents an extended layer of reduced water mobility relative to bulk water. The dynamic hydration layer extends far beyond the static hydration layer, adjacent to the polymer.

Model-Dependent Solvation of the K-18 Domain of the Intrinsically Disordered Protein Tau

A known imbalance between intra-protein and protein-water interactions in many empirical force fields results in collapsed conformational ensembles of intrinsically disordered proteins in explicit solvent simulations that disagree with experiments. Multiple strategies have been introduced in the literature to modify protein-water interactions, which improve agreement between experiments and simulations. In this work, we combine simulations with standard and modified force fields with a spatially resolved analysis of solvation free energy contributions and compare the consequences of each strategy. We find that enhanced Lennard-Jones (LJ) interactions between protein atoms and water oxygens primarily improve the solvation of nonpolar functional groups of the protein. In contrast, modified electrostatics in the water model or strengthened LJ interactions between the protein and water hydrogens mainly affect the hydration of polar functional groups. Modified electrostatics further impact the average orientation of water molecules in the hydration shell. As a result, protein-water interactions with the first hydration layers are strengthened, while interactions with water molecules in higher hydration shells are weakened. Hence, distinct strategies to balance intra-protein and protein-water interactions in simulations have qualitatively different effects on protein solvation. These differences are not necessarily captured by comparisons to experiments that report on global parameters describing protein conformational ensembles, e.g., the radius of gyration, but will influence the tendency of a protein to form aggregates or phase-separated droplets.

Solvation free energy arithmetic for small organic molecules

Solvent-mediated interactions contribute to ligand binding affinities in computational drug design and provide a challenge for theoretical predictions. In this study, we analyze the solvation free energy of benzene derivatives in water to guide the development of predictive models for solvation free energies and solvent-mediated interactions. We use a spatially resolved analysis of local solvation free energy contributions and define solvation free energy arithmetic, which enable us to construct additive models to describe the solvation of complex compounds. The substituents analyzed in this study are carboxyl and nitro-groups due to their similar sterical requirements but distinct interactions with water. We find that nonadditive solvation free energy contributions are primarily attributed to electrostatics, which are qualitatively reproduced with computationally efficient continuum models. This suggests a promising route for the development of efficient and accurate models for the solvation of complex molecules with varying substitution patterns using solvation arithmetic.

Dynamic Electric Field Alignment Determines the Water Rotational Motion around Protein

Water rotational dynamics in biomolecular solution is crucial to evaluating and controlling biomolecule stability. In this molecular dynamics simulation (MD) study on lysozyme solutions, we present how the exerted internal electric field determines water rotational dynamics. We find that the relaxation time of water rotation is equivalent to that of the reorientation of the exerted overall electric field for every single water molecule, regardless of its translation mode. Namely, water molecular rotation synchronizes with the exerted field reorientation. We also map the reorientation process of the electric field at fixed points relative to protein in the solution, which displays the local hydration dynamics commensurate with the reported time-dependent fluorescence Stokes shift (TDFSS) measurements. Comparing the spatial distribution of local field reorientation relaxation time with that of rotational relaxation time, we further suggest that water rotation dynamics are subject to the reorientation of the local overall field within the hydration layer. While outside the hydration layer, the relaxation time of the local electric field reorientation is short enough (subpicosecond) to assume the δ function, showing the electric force with randomly changing orientation is applied to each water molecule.

Understanding dynamics in coarse-grained models. I. Universal excess entropy scaling relationship

Coarse-grained (CG) models facilitate an efficient exploration of complex systems by reducing the unnecessary degrees of freedom of the fine-grained (FG) system while recapitulating major structural correlations. Unlike structural properties, assessing dynamic properties in CG modeling is often unfeasible due to the accelerated dynamics of the CG models, which allows for more efficient structural sampling. Therefore, the ultimate goal of the present series of articles is to establish a better correspondence between the FG and CG dynamics. To assess and compare dynamical properties in the FG and the corresponding CG models, we utilize the excess entropy scaling relationship. For Paper I of this series, we provide evidence that the FG and the corresponding CG counterpart follow the same universal scaling relationship. By carefully reviewing and examining the literature, we develop a new theory to calculate excess entropies for the FG and CG systems while accounting for entropy representability. We demonstrate that the excess entropy scaling idea can be readily applied to liquid water and methanol systems at both the FG and CG resolutions. For both liquids, we reveal that the scaling exponents remain unchanged from the coarse-graining process, indicating that the scaling behavior is universal for the same underlying molecular systems. Combining this finding with the concept of mapping entropy in CG models, we show that the missing entropy plays an important role in accelerating the CG dynamics.

Total free energy analysis of fully hydrated proteins

The total free energy of a hydrated biomolecule and its corresponding decomposition of energy and entropy provides detailed information about regions of thermodynamic stability or instability. The free energies of four hydrated globular proteins with different net charges are calculated from a molecular dynamics simulation, with the energy coming from the system Hamiltonian and entropy using multiscale cell correlation. Water is found to be most stable around anionic residues, intermediate around cationic and polar residues, and least stable near hydrophobic residues, especially when more buried, with stability displaying moderate entropy-enthalpy compensation. Conversely, anionic residues in the proteins are energetically destabilized relative to singly solvated amino acids, while trends for other residues are less clear-cut. Almost all residues lose intraresidue entropy when in the protein, enthalpy changes are negative on average but may be positive or negative, and the resulting overall stability is moderate for some proteins and negligible for others. The free energy of water around single amino acids is found to closely match existing hydrophobicity scales. Regarding the effect of secondary structure, water is slightly more stable around loops, of intermediate stability around β strands and turns, and least stable around helices. An interesting asymmetry observed is that cationic residues stabilize a residue when bonded to its N-terminal side but destabilize it when on the C-terminal side, with a weaker reversed trend for anionic residues.

Stochastic Analysis of Molecular Dynamics Reveals the Rotation Dynamics Distribution of Water around Lysozyme

Water dynamics is essential to biochemical processes by mediating all such reactions, including biomolecular degeneration in solutions. To disentangle the molecular-scale distribution of water dynamics around a solute biomolecule, we investigated here the rotational dynamics of water around lysozyme by combining molecular dynamics (MD) simulations and broadband dielectric spectroscopy (BDS). A statistical analysis using the relaxation times and trajectories of every single water molecule was proposed, and the two-dimensional probability distribution of water at a distance from the lysozyme surface with a rotational relaxation time was given. For the observed lysozyme solutions of 34-284 mg/mL, we discovered that the dielectric relaxation time obtained from this distribution agrees well with the measured γ relaxation time, which suggests that rotational self-correlation of water molecules underlies the gigahertz domain of the dielectric spectra. Regardless of protein concentration, water rotational relaxation time versus the distance from the lysozyme surface revealed that the water rotation is severely retarded within 3 Å from the lysozyme surface and is nearly comparable to pure water when farther than 10 Å. The dimension of the first hydration layer was subsequently identified in terms of the relationship between the acceleration of water rotation and the distance from the protein surface.

Spatially Resolved Hydration Thermodynamics in Biomolecular Systems

Water is essential for the structure, dynamics, energetics, and thus the function of biomolecules. It is a formidable challenge to elicit, in microscopic detail, the role of the solvation-related driving forces of biomolecular processes, such as the enthalpy and entropy contributions to the underlying free-energy landscape. In this Perspective, we discuss recent developments and applications of computational methods that provide a spatially resolved map of hydration thermodynamics in biomolecular systems and thus yield atomic-level insights to guide the interpretation of experimental observations. An emphasis is on the challenge of quantifying the hydration entropy, which requires characterization of both the motions of the biomolecules and of the water molecules in their surrounding.

Hydrolysis Mechanism of the Linkers by Matrix Metalloproteinase-9 Using QM/MM Calculations

Activatable cell-penetrating peptides (ACPPs) are known to be able to decrease the cytotoxicity of cell-penetrating peptide (CPP)-based drug delivery systems. Furthermore, they can improve the targeting of CPPs when specifically recognized and hydrolyzed by characteristic proteases. A comprehensive and profound understanding of the recognition and hydrolysis process will provide a better design of the ACPP-based drug delivery system. Previous studies have clearly described how ACPPs are recognized and bound by MMPs. However, the hydrolysis mechanism of ACPPs is still unsolved. This work focuses on a proteinase-sensitive cleavable linker of ACPPs (PLGLAG), the key structure for recognition and hydrolysis, trying to determine the mechanism by which MMP-9 hydrolyzes its substrate PLGLAG. The quantum mechanics/molecular mechanics (QM/MM) calculations herein show that MMP-9 proteolysis is a water-mediated four-step reaction. More specifically, it consists of (i) nucleophilic attack, (ii) hydrogen-bond rearrangement, (iii) proton transfer, and finally (iv) amide bond rupture. Considering the reversibility of multistep reaction, the second step (i.e., hydrogen-bond rearrangement) has the highest barrier and is the rate-limiting step in the hydrolysis of PLGLAG. The possible design and improvement of the key P1 and P1' sites are also explored through mutations. The present results indicate that, while the mutations affect the reaction energy barriers and the rate-limiting steps, all mutants considered could be hydrolyzed by MMP-9. To provide further insights, the hydrolysis mechanism of MMP-2, which has a similar hydrolysis process to that of MMP-9 but with different reaction barriers, is also studied and compared. As a result, this work provides detailed insights into the hydrolysis mechanism of ACPPs by MMP-9 and, thus, also possible insights for the development of new strategies for ACPP-based delivery systems.

Spectrally Resolved Estimation of Water Entropy in the Active Site of Human Carbonic Anhydrase II

A major challenge in understanding ligand binding to biomacromolecules lies in dissecting the underlying thermodynamic driving forces at the atomic level. Quantifying the contributions of water molecules is often especially demanding, although they can play important roles in biomolecular recognition and binding processes. One example is human carbonic anhydrase II, whose active site harbors a conserved network of structural water molecules that are essential for enzymatic catalysis. Inhibitor binding disrupts this water network and changes the hydrogen-bonding patterns in the active site. Here, we use atomistic molecular dynamics simulations to compute the absolute entropy of the individual water molecules confined in the active site of hCAII using a spectrally resolved estimation (SRE) approach. The entropy decrease of water molecules that remain in the active site upon binding of a dorzolamide inhibitor is caused by changes in hydrogen bonding and stiffening of the hydrogen-bonding network. Overall, this entropy decrease is overcompensated by the gain due to the release of three water molecules from the active site upon inhibitor binding. The spectral density calculations enable the assignment of the changes to certain vibrational modes. In addition, the range of applicability of the SRE approximation is systematically explored by exploiting the gradually changing degree of immobilization of water molecules as a function of the distance to a phospholipid bilayer surface, which defines an "entropy ruler". These results demonstrate the applicability of SRE to biomolecular solvation, and we expect it to become a useful method for entropy calculations in biomolecular systems.

Protein flexibility reduces solvent-mediated friction barriers of ligand binding to a hydrophobic surface patch

Solvent fluctuations have been explored in detail for idealized and rigid hydrophobic model systems, but so far it has remained unclear how internal protein motions and their coupling to the surrounding solvent affect the dynamics of ligand binding to biomolecular surfaces. Here, molecular dynamics simulations were used to elucidate the solvent-mediated binding of a model ligand to the hydrophobic surface patch of ubiquitin. The ligand's friction profiles reveal pronounced long-time correlations and enhanced friction in the vicinity of the protein, similar to idealized hydrophobic surfaces. Interestingly, these effects are shaped by internal protein motions. Protein flexibility modulates water density fluctuations near the hydrophobic surface patch and smooths out the friction profile of ligand binding.

Destabilisation of the structure of transthyretin is driven by Ca2

Tetrameric transthyretin (TTR) transports thyroid hormones and retinol in plasma and cerebrospinal fluid and performs protective functions under stress conditions. Ageing and mutations result in TTR destabilisation and the formation of the amyloid deposits that dysregulate Ca²⁺ homeostasis. Our aim was to determine whether Ca²⁺ affects the structural stability of TTR. We show, using multiple techniques, that Ca²⁺ does not induce prevalent TTR dissociation and/or oligomerisation. However, in the presence of Ca²⁺, TTR exhibits altered conformational flexibility and different interactions with the solvent molecules. These structural changes lead to the formation of the sub-populations of non-native TTR conformers and to the destabilisation of the structure of TTR. Moreover, the sub-population of TTR molecules undergoes fragmentation that is augmented by Ca²⁺. We postulate that Ca²⁺ constitutes the structural and functional switch between the native and non-native forms of TTR, and therefore tip the balance towards age-dependent pathological calcification.

Comparison of Free-Energy Methods to Calculate the Barriers for the Nucleophilic Substitution of Alkyl Halides by Hydroxide

Calculating the free-energy barriers of liquid-phase chemical reactions with explicit solvent is a considerable challenge. Most studies use the energy and entropy of minimized single-point geometries of the reactants and transition state in implicit solvent using normal mode analysis (NMA). Explicit-solvent methods instead make use of the potential of mean force (PMF). Here, we propose a new energy-entropy (EE) method to calculate the Gibbs free energy of reactants and transition states in explicit solvent by combining quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations with multiscale cell correlation (MCC). We apply it to six nucleophilic substitution reactions of the hydroxide transfer to methyl and ethyl halides in water, where the halides are F, Cl, and Br. We compare EE-MCC Gibbs free energy barriers using two Hamiltonians, self-consistent charge density functional based tight-binding (SCC-DFTB) and B3LYP/6-31+G* density functional theory (DFT) with respective PMF values, EE-NMA values using B3LYP/6-31+G* and M06/6-31+G* DFT in implicit solvent and experimental values derived via transition state theory. The barriers using SCC-DFTB are found to agree well with the PMF and experiment and previous computational studies, being slightly higher but improving on the lower values obtained for the implicit solvent. Achieving convergence over many degrees of freedom remains a challenge for EE-MCC in explicit-solvent QM/MM systems, particularly for the more expensive B3LYP/6-31+G* and M06/6-31+G* DFT methods, but the insightful decomposition of entropy over all degrees of freedom should make EE-MCC a valuable tool for deepening the understanding of chemical reactions.

Dielectrophoresis of proteins: experimental data and evolving theory

The ability to selectively move and trap proteins is core to their effective use as building blocks and for their characterization. Analytical and preparative strategies for proteins have been pursued and modeled for nearly a hundred years, with great advances and success. Core to all of these studies is the separation, isolation, purification, and concentration of pure homogeneous fractions of a specific protein in solution. Processes to accomplish this useful solution include biphasic equilibrium (chromatographies, extractions), mechanical, bulk property, chemical equilibria, and molecular recognition. Ultimately, the goal of all of these is to physically remove all non-like protein molecules-to the finest detail: all atoms in the full three-dimensional structure being identical down the chemical bond and bulk structure chirality. One strategy which has not been effectively pursued is exploiting the higher order subtle electrical properties of the protein-solvent system. The advent of microfluidic systems has enabled the use of very high electric fields and well-defined gradients such that extremely high resolution separations of protein mixtures are possible. These advances and recognition of these capabilities have caused a re-evaluation of the underlying theoretical models and they were found to be inadequate. New theoretical descriptions are being considered which align more closely to the total forces present and the subtlety of differences between similar proteins. These are focused on the interfacial area between the protein and hydrating solvent molecules, as opposed to the macroscale assumptions of homogeneous solutions and particles. This critical review examines all data which has been published that place proteins in electric field gradients which induce collection of those proteins, demonstrating a force greater than dispersive effects or countering forces. Evolving theoretical constructs are presented and discussed, and a general estimate of future capabilities using the higher order effects and the high fields and precise gradients of microfluidic systems is discussed. Graphical abstract.

Mechanical Unfolding of Spectrin Repeats Induces Water-Molecule Ordering

Mechanical processes are involved at many stages of the development of living cells, and often external forces applied to a biomolecule result in its unfolding. Although our knowledge of the unfolding mechanisms and the magnitude of the forces involved has evolved, the role that water molecules play in the mechanical unfolding of biomolecules has not yet been fully elucidated. To this end, we investigated with steered molecular dynamics simulations the mechanical unfolding of dystrophin's spectrin repeat 1 and related the changes in the protein's structure to the ordering of the surrounding water molecules. Our results indicate that upon mechanically induced unfolding of the protein, the solvent molecules become more ordered and increase their average number of hydrogen bonds. In addition, the unfolded structures originating from mechanical pulling expose an increasing amount of the hydrophobic residues to the solvent molecules, and the uncoiled regions adapt a convex surface with a small radius of curvature. As a result, the solvent molecules reorganize around the protein's small protrusions in structurally ordered waters that are characteristic of the so-called "small-molecule regime," which allows water to maintain a high hydrogen bond count at the expense of an increased structural order. We also determined that the response of water to structural changes in the protein is localized to the specific regions of the protein that undergo unfolding. These results indicate that water plays an important role in the mechanically induced unfolding of biomolecules. Our findings may prove relevant to the ever-growing interest in understanding macromolecular crowding in living cells and their effects on protein folding, and suggest that the hydration layer may be exploited as a means for short-range allosteric communication.

Role of hydration water in the onset of protein structural dynamics

Proteins are the molecular workhorses in a living organism. Their 3D structures are animated by a multitude of equilibrium fluctuations and specific out-of-equilibrium motions that are required for proteins to be biologically active. When studied as a function of temperature, functionally relevant dynamics are observed at and above the so-called protein dynamical transition (~240 K) in hydrated, but not in dry proteins. In this review we present and discuss the main experimental and computational results that provided evidence for the dynamical transition, with a focus on the role of hydration water dynamics in sustaining functional protein dynamics. The coupling and mutual influence of hydration water dynamics and protein dynamics are discussed and the hypotheses illustrated that have been put forward to explain the physical origin of their onsets.

Heterogeneity of water structure and dynamics at the protein-water interface

In this molecular dynamics simulation study, we analyze the local structural and dynamic properties of water hydrating the protein ubiquitin on a spatial grid with 1 Å resolution. This allows for insights into the spatial distribution of water number densities, molecular orientations, translations, and rotations as a function of distance from the protein surface. Water molecule orientations follow a heterogeneous distribution with preferred local orientations of water dipoles and O-H bond vectors up to 10-15 Å distances from the protein, while local variations of the water number density converge to homogeneous bulk-like values within less than 8 Å. Interestingly, we find that the long-ranged orientational structure of water does not impact either the translational or rotational dynamics of water. Instead, heterogeneous distributions of local dynamical parameters and averaged dynamical retardation factors are only found close to the protein surface and follow a distance dependence comparable to heterogeneities in the local water number density. This study shows that the formation of nanodomains of preferred water orientations far from the protein does not significantly impact dynamical processes probed as a non-local average in most experiments.

Structural Relaxation Processes and Collective Dynamics of Water in Biomolecular Environments

In this simulation study, we investigate the influence of biomolecular confinement on dynamical processes in water. We compare water confined in a membrane protein nanopore at room temperature to pure liquid water at low temperatures with respect to structural relaxations, intermolecular vibrations, and the propagation of collective modes. We observe distinct potential energy landscapes experienced by water molecules in the two environments, which nevertheless result in comparable hydrogen bond lifetimes and sound propagation velocities. Hence, we show that a viscoelastic argument that links slow rearrangements of the water-hydrogen bond network to ice-like collective properties applies to both, the pure liquid and biologically confined water, irrespective of differences in the microscopic structure.

Overhauser Dynamic Nuclear Polarization for the Study of Hydration Dynamics, Explained

We outline the physical properties of hydration water that are captured by Overhauser Dynamic Nuclear Polarization (ODNP) relaxometry and explore the insights that ODNP yields about the water and the surface that this water is coupled to. As ODNP relies on the pairwise cross-relaxation between the electron spin of a spin probe and a proton nuclear spin of water, it captures the dynamics of single-particle diffusion of an ensemble of water molecules moving near the spin probe. ODNP principally utilizes the same physics as other nuclear magnetic resonance (NMR) relaxometry (i.e., relaxation measurement) techniques. However, in ODNP, electron paramagnetic resonance (EPR) excites the electron spins probes and their high net polarization acts as a signal amplifier. Furthermore, it renders ODNP parameters highly sensitive to water moving at rates commensurate with the EPR frequency of the spin probe (typically 10GHz). Also, ODNP selectively enhances the NMR signal contributions of water moving within close proximity to the spin label. As a result, ODNP can capture ps-ns movements of hydration waters with high sensitivity and locality, even in samples with protein concentrations as dilute as 10 µM. To date, the utility of the ODNP technique has been demonstrated for two major applications: the characterization of the spatial variation in the properties of the hydration layer of proteins or other surfaces displaying topological diversity, and the identification of structural properties emerging from highly disordered proteins and protein domains. The former has been shown to correlate well with the properties of hydration water predicted by MD simulations and has been shown capable of evaluating the hydrophilicity or hydrophobicity of a surface. The latter has been demonstrated for studies of an interhelical loop of proteorhodopsin, the partial structure of α-synuclein embedded at the lipid membrane surface, incipient structures adopted by tau proteins en route to fibrils, and the structure and hydration profile of a transmembrane peptide. This chapter focuses on offering a mechanistic understanding of the ODNP measurement and the molecular dynamics encoded in the ODNP parameters. In particular, it clarifies how the electron-nuclear dipolar coupling encodes information about the molecular dynamics in the nuclear spin self-relaxation and, more importantly, the electron-nuclear spin cross-relaxation rates. The clarification of the molecular dynamics underlying ODNP should assist in establishing a connection to theory and computer simulation that will offer far richer interpretations of ODNP results in future studies.

Macromolecular crowding effects in flexible polymer solutions

Biological environments are often “crowded” due to high concentrations (300–400[Formula: see text]g/L) of macromolecules. Computational modeling approaches like Molecular Dynamics (MD), rigid-body Brownian Dynamics and Monte Carlo simulations have recently emerged, which allow to study the effects macromolecular crowding at a microscopic level and to provide complementary information to experiments. Here, we use a recently introduced multiple-conformation Monte Carlo (mcMC) approach in order to study the influence of intermolecular interactions on the structural equilibrium of flexible polyethylene glycol (PEG) polymers under self-crowding conditions. The large conformational space accessible to PEG polymers allows us to evaluate the general applicability of the mcMC approach, which describes the intramolecular degrees of freedom by a finite-size ensemble of discrete conformations. Despite the simplicity of the approach, we show that influences of intermolecular interactions on the intramolecular free energy surface can be described qualitatively using mcMC. By varying the magnitude of distinct terms in the intermolecular potential, we can further study the compensating effects of repulsive and nonspecific attractive intermolecular interactions, which favor compact and extended polymer states, respectively. We use our simulation results to derive an analytical model that describes the effects of intermolecular interactions on the stability of PEG polymer conformations as a function of the radius of gyration and the corresponding solvent accessible surface. We use this model to confirm the role of molecular surfaces for attractive interactions that can counteract excluded volume effects. Extrapolation of the model further allows for the analysis of scenarios that are not easily accessible to direct simulations as described here.

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.