Water Dynamics at Protein Interfaces: Ultrafast Optical Kerr Effect Study

Dynamics of Acrylamide Hydrogels, Polymers, and Monomers in Water Measured with Optical Heterodyne-Detected Optical Kerr Effect Spectroscopy

The ultrafast dynamics of acrylamide monomers (AAm), polyacrylamide (PAAm), and polyacrylamide hydrogels (PAAm-HG) in water were studied using optical heterodyne-detected optical Kerr effect (OHD-OKE) spectroscopy. Previous ultrafast infrared (IR) measurements of the water dynamics showed that at the same concentration of the acrylamide moiety, AAm, PAAm, and PAAm-HG exhibited identical water dynamics and that these dynamics slowed with increasing concentration. In contrast to the IR measurements, OHD-OKE experiments measure the dynamics of both the water and the acrylamide species, which occur on different time scales. In this study, the dynamics of all the acrylamide systems slowed with increasing concentration. We found that AAm exhibits tetraexponential decays, the longest component of which followed Debye-Stokes-Einstein behavior except for the highest concentration, 40% (w/v). Low concentrations of PAAm followed a single power law decay, while high concentrations of PAAm and all concentrations of PAAm-HG decayed with two power laws. The highest concentrations, 25% and 40%, of PAAm and PAAm-HG showed nearly identical dynamics. We interpreted this result as reflecting a similar extent of chain-chain interactions. At low concentrations, PAAm displays non-Markovian, single-chain dynamics (single power law), but PAAm displays entangled chain-chain interactions at high concentrations (two power laws). PAAm-HG has chain-chain interactions at all concentrations that arise from the cross-linking. At high concentrations, the dynamics of the entangled of PAAm become identical within error as those of the cross-linked PAAm-HG.

Protein-Ligand Binding Molecular Details Revealed by Terahertz Optical Kerr Spectroscopy: A Simulation Study

Picosecond fast motions and their involvement in the biochemical processes such as protein-ligand binding has engaged significant attention. Terahertz optical Kerr spectroscopy (OKE) has the superior potential to probe these fast motions directly. Application of OKE in protein-ligand binding study is, however, limited by the difficulty of quantitative atomistic interpretation, and the calculation of Kerr spectrum for entire solvated protein complex was considered not yet feasible, due to the lack of one consistent polarizable model for both configuration sampling and polarizability calculation. Here, we analyzed the biochemical relevance of OKE to the lysozyme-triacetylchitotriose binding based on the first OKE simulation using one consistent Drude polarizable model. An analytical multipole and induced dipole scheme was employed to calculate the off-diagonal Drude polarizability more efficiently and accurately. Further theoretical analysis revealed how the subtle twisting and stiffening of aromatic protein residues' spatial arrangement as well as the confinement of small water clusters between ligand and protein cavity due to the ligand binding can be examined using Kerr spectroscopy. Comparison between the signals of bound complex and that of uncorrelated protein/ligand demonstrated that binding action alone has reflection in the OKE spectrum. Our study indicated OKE as a powerful terahertz probe for protein-ligand binding chemistry and dynamics.

Local Mutations Can Serve as a Game Changer for Global Protein Solvent Interaction

Although it is well-known that limited local mutations of enzymes, such as matrix metalloproteinases (MMPs), may change enzyme activity by orders of magnitude as well as its stability, the completely rational design of proteins is still challenging. These local changes alter the electrostatic potential and thus local electrostatic fields, which impacts the dynamics of water molecules close the protein surface. Here we show by a combined computational design, experimental, and molecular dynamics (MD) study that local mutations have not only a local but also a global effect on the solvent: In the specific case of the matrix metalloprotease MMP14, we found that the nature of local mutations, coupled with surface morphology, have the ability to influence large patches of the water hydrogen-bonding network at the protein surface, which is correlated with stability. The solvent contribution can be experimentally probed via terahertz (THz) spectroscopy, thus opening the door to the exciting perspective of rational protein design in which a systematic tuning of hydration water properties allows manipulation of protein stability and enzymatic activity.

Toward theoretical terahertz spectroscopy of glassy aqueous solutions: partially frozen solute-solvent couplings of glycine in water

The molecular-level understanding of THz spectra of aqueous solutions under ambient conditions has been greatly advanced in recent years. Here, we go beyond previous analyses by performing ab initio molecular dynamics simulations of glycine in water with artificially frozen solute or solvent molecules, respectively, while computing the total THz response as well as its decomposition into mode-specific resonances based on the "supermolecular solvation complex" technique. Clamping the water molecules and keeping glycine moving breaks the coupling of glycine to the structural dynamics of the solvent, however, the polarization and dielectric solvation effects in the static solvation cage are still at work since the full electronic structure of the quenched solvent is taken into account. The complementary approach of fixing glycine reveals both the dynamical and electronic response of the solvation cage at the level of its THz response. Moreover, to quantitatively account for the electronic contribution solely due to solvent embedding, the solute species is "vertically desolvated", thus preserving the fully coupled solute-solvent motion in terms of the solute's structural dynamics in solution, while its electronic structure is no longer subject to solute-solvent polarization and charge transfer effects. When referenced to the free simulation of Gly(aq), this three-fold approach allows us to decompose the THz spectral contributions due to the correlated solute-solvent dynamics into entirely structural and purely electronic effects. Beyond providing hitherto unknown insights, the observed systematic changes of THz spectra in terms of peak shifts and lineshape modulations due to conformational freezing and frozen solvation cages might be useful to investigate the solvation of molecules in highly viscous H-bonding solvents such as ionic liquids and even in cryogenic ices as relevant to polar stratospheric and dark interstellar clouds.

Picosecond orientational dynamics of water in living cells

Cells are extremely crowded, and a central question in biology is how this affects the intracellular water. Here, we use ultrafast vibrational spectroscopy and dielectric-relaxation spectroscopy to observe the random orientational motion of water molecules inside living cells of three prototypical organisms: Escherichia coli, Saccharomyces cerevisiae (yeast), and spores of Bacillus subtilis. In all three organisms, most of the intracellular water exhibits the same random orientational motion as neat water (characteristic time constants ~9 and ~2 ps for the first-order and second-order orientational correlation functions), whereas a smaller fraction exhibits slower orientational dynamics. The fraction of slow intracellular water varies between organisms, ranging from ~20% in E. coli to ~45% in B. subtilis spores. Comparison with the water dynamics observed in solutions mimicking the chemical composition of (parts of) the cytosol shows that the slow water is bound mostly to proteins, and to a lesser extent to other biomolecules and ions.

The cytoplasm’s crowdedness leads one to expect that cell water is different from bulk water. By measuring the rotational motion of water molecules in living cells, Tros et al. find that apart from a small fraction of water solvating biomolecules, cell water has the same dynamics as bulk water.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

Perspective: Structure and ultrafast dynamics of biomolecular hydration shells

The structure and function of biomolecules can be strongly influenced by their hydration shells. A key challenge is thus to determine the extent to which these shells differ from bulk water, since the structural fluctuations and molecular excitations of hydrating water molecules within these shells can cover a broad range in both space and time. Recent progress in theory, molecular dynamics simulations, and ultrafast vibrational spectroscopy has led to new and detailed insight into the fluctuations of water structure, elementary water motions, and electric fields at hydrated biointerfaces. Here, we discuss some central aspects of these advances, focusing on elementary molecular mechanisms and processes of hydration on a femto- to picosecond time scale, with some special attention given to several issues subject to debate.

Hydrogen Bond Network of Water around Protein Investigated with Terahertz and Infrared Spectroscopy

The dynamical and structural properties of water at protein interfaces were characterized on the basis of the broadband complex dielectric constant (0.25 to 400 THz) of albumin aqueous solutions. Our analysis of the dielectric responses between 0.25 and 12 THz first revealed hydration water with retarded reorientational dynamics extending ∼8.5 Å (corresponding to three to four layers) out from the albumin surface. Second, the number of nonhydrogen-bonded water was decreased in the presence of the albumin solute, indicating protein inhibits the fragmentation of the water hydrogen-bond network. Finally, water molecules at the albumin interface were found to form a distorted hydrogen-bond structure due to topological and energetic disorder of the protein surface. In addition, the intramolecular O-H stretching vibration of water (∼100 THz), which is sensitive to hydrogen-bond environment, pointed to a trend that hydration water has a larger population of strongly hydrogen-bonded water molecules compared with that of bulk water. From these experimental results, we concluded that the "strengthened" water hydrogen bonds at the protein interface dynamically slow down the reorientational motion of water and form the less-defective hydrogen-bond network by inhibiting the fragmentation of water-water hydrogen bonds. Nevertheless, such a strengthened water hydrogen-bond network is composed of heterogeneous hydrogen-bond distances and angles, and thus characterized as structurally "distorted."

New Method to Study the Vibrational Modes of Biomolecules in the Terahertz Range Based on a Single-Stage Raman Spectrometer

The low-frequency vibrational (LFV) modes of biomolecules reflect specific intramolecular and intermolecular thermally induced fluctuations that are driven by external perturbations, such as ligand binding, protein interaction, electron transfer, and enzymatic activity. Large efforts have been invested over the years to develop methods to access the LFV modes due to their importance in the studies of the mechanisms and biological functions of biomolecules. Here, we present a method to measure the LFV modes of biomolecules based on Raman spectroscopy that combines volume holographic filters with a single-stage spectrometer, to obtain high signal-to-noise-ratio spectra in short acquisition times. We show that this method enables LFV mode characterization of biomolecules even in a hydrated environment. The measured spectra exhibit distinct features originating from intra- and/or intermolecular collective motion and lattice modes. The observed modes are highly sensitive to the overall structure, size, long-range order, and configuration of the molecules, as well as to their environment. Thus, the LFV Raman spectrum acts as a fingerprint of the molecular structure and conformational state of a biomolecule. The comprehensive method we present here is widely applicable, thus enabling high-throughput study of LFV modes of biomolecules.

Molecular properties of aqueous solutions: a focus on the collective dynamics of hydration water

When a solute is dissolved in water, their mutual interactions determine the molecular properties of the solute on one hand, and the structure and dynamics of the surrounding water particles (the so-called hydration water) on the other. The very existence of soft matter and its peculiar properties are largely due to the wide variety of possible water-solute interactions. In this context, water is not an inert medium but rather an active component, and hydration water plays a crucial role in determining the structure, stability, dynamics, and function of matter. This review focuses on the collective dynamics of hydration water in terms of retardation with respect to the bulk, and of the number of molecules whose dynamics is perturbed. Since water environments are in a dynamic equilibrium, with molecules continuously exchanging from around the solute towards the bulk and vice versa, we examine the ability of different techniques to measure the water dynamics on the basis of the explored time scales and exchange rates. Special emphasis is given to the collective dynamics probed by extended depolarized light scattering and we discuss whether and to what extent the results obtained in aqueous solutions of small molecules can be extrapolated to the case of large biomacromolecules. In fact, recent experiments performed on solutions of increasing complexity clearly indicate that a reductionist approach is not adequate to describe their collective dynamics. We conclude this review by presenting current ideas that are being developed to describe the dynamics of water interacting with macromolecules.

Water Determines the Structure and Dynamics of Proteins

Water is an essential participant in the stability, structure, dynamics, and function of proteins and other biomolecules. Thermodynamically, changes in the aqueous environment affect the stability of biomolecules. Structurally, water participates chemically in the catalytic function of proteins and nucleic acids and physically in the collapse of the protein chain during folding through hydrophobic collapse and mediates binding through the hydrogen bond in complex formation. Water is a partner that slaves the dynamics of proteins, and water interaction with proteins affect their dynamics. Here we provide a review of the experimental and computational advances over the past decade in understanding the role of water in the dynamics, structure, and function of proteins. We focus on the combination of X-ray and neutron crystallography, NMR, terahertz spectroscopy, mass spectroscopy, thermodynamics, and computer simulations to reveal how water assist proteins in their function. The recent advances in computer simulations and the enhanced sensitivity of experimental tools promise major advances in the understanding of protein dynamics, and water surely will be a protagonist.

Salt effects on the picosecond dynamics of lysozyme hydration water investigated by terahertz time-domain spectroscopy and an insight into the Hofmeister series for protein stability and solubility

The addition of salts into protein aqueous solutions causes changes in protein solubility and stability, the ability of which is known to be ordered in the Hofmeister series.

Femtosecond optical Kerr effect setup with signal "live view" for measurements in the solid, liquid, and gas phases

We present the experimental setup constructed in our laboratory for measurement of the femtosecond optical Kerr effect. The setup allows measurements with high temporal resolution and acquisition speed. The high signal to noise ratio is obtained with use of a homemade balanced detector. Due to the high acquisition speed and good signal to noise ratio, it is possible to have a "live view" of the signal and to easily tune the sample position and orientation before the measurement. We show the example results obtained in the solid, liquid, and the gas phases and we use them in order to check on the precision of our setup. As the samples we have used a YAG crystal, liquid acetone, and atmospheric air. In the latter two cases, a good agreement with the literature data has been found. The measurements in the gas phase confirm that our setup, although utilizing low energy pulses from the sapphire oscillator, is able to acquire high quality rotational signal in a low density sample.

Comparative study of hydration shell dynamics around a hyperactive antifreeze protein and around ubiquitin

The hydration layer surrounding a protein plays an essential role in its biochemical function and consists of a heterogeneous ensemble of water molecules with different local environments and different dynamics. What determines the degree of dynamical heterogeneity within the hydration shell and how this changes with temperature remains unclear. Here, we combine molecular dynamics simulations and analytic modeling to study the hydration shell structure and dynamics of a typical globular protein, ubiquitin, and of the spruce budworm hyperactive antifreeze protein over the 230-300 K temperature range. Our results show that the average perturbation induced by both proteins on the reorientation dynamics of water remains moderate and changes weakly with temperature. The dynamical heterogeneity arises mostly from the distribution of protein surface topographies and is little affected by temperature. The ice-binding face of the antifreeze protein induces a short-ranged enhancement of water structure and a greater slowdown of water reorientation dynamics than the non-ice-binding faces whose effect is similar to that of ubiquitin. However, the hydration shell of the ice-binding face remains less tetrahedral than the bulk and is not "ice-like". We finally show that the hydrogen bonds between water and the ice-binding threonine residues are particularly strong due to a steric confinement effect, thereby contributing to the strong binding of the antifreeze protein on ice crystals.

Hydration shells of proteins probed by depolarized light scattering and dielectric spectroscopy: orientational structure is significant, positional structure is not

Water interfacing hydrated proteins carry properties distinct from those of the bulk and is often described as a separate entity, a "biological water." We address here the question of which dynamical and structural properties of hydration water deserve this distinction. The study focuses on different aspects of the density and orientational fluctuations of hydration water and the ability to separate them experimentally by combining depolarized light scattering with dielectric spectroscopy. We show that the dynamics of the density fluctuations of the hydration shells reflect the coupled dynamics of the solute and solvent and do not require a special distinction as "biological water." The orientations of shell water molecules carry dramatically different physics and do require a separation into a sub-ensemble. Depending on the property considered, the perturbation of water's orientational structure induced by the protein propagates 3-5 hydration shells into the bulk at normal temperature.

One- and two-photon absorption in solution: the effects of a passive auxiliary beam

The efficiencies of one- and two-photon absorption by chromophores in solution may be significantly modified by a sufficiently intense beam of off-resonant light. A molecular analysis based on quantum electrodynamics (QED) fully accounts for this phenomenon of laser-modified absorption. A time-dependent perturbation-theory treatment describes the process in terms of stimulated forward Rayleigh-scattering of the auxiliary beam occurring simultaneously with the absorption interaction(s). Our formulation accommodates media modifications to the basic character of light-matter interactions, taking into account the refractive and dispersive properties of a solution-phase environment. This introduces the bulk refractive index of the solvent directly into the QED framework. The measurable electronic response of molecules freely rotating in solution is defined by an average of all orientations. We explicitly derive fixed-orientation and rotationally averaged calculations for the Fermi-rule rate of laser-modified one- and two-photon absorption. For a given beam polarization geometry, the solution-phase molecular response is expressible as a set of natural invariant scalars. These results reveal details of the dependence on the beam polarisations and on the rotationally averaged molecular response: we illustrate the breadth of variation available via geometric manipulation of beam polarization, and raise new possibilities for quantum weak measurements of laser states.

Water dynamics in protein hydration shells: the molecular origins of the dynamical perturbation

Protein hydration shell dynamics play an important role in biochemical processes including protein folding, enzyme function, and molecular recognition. We present here a comparison of the reorientation dynamics of individual water molecules within the hydration shell of a series of globular proteins: acetylcholinesterase, subtilisin Carlsberg, lysozyme, and ubiquitin. Molecular dynamics simulations and analytical models are used to access site-resolved information on hydration shell dynamics and to elucidate the molecular origins of the dynamical perturbation of hydration shell water relative to bulk water. We show that all four proteins have very similar hydration shell dynamics, despite their wide range of sizes and functions, and differing secondary structures. We demonstrate that this arises from the similar local surface topology and surface chemical composition of the four proteins, and that such local factors alone are sufficient to rationalize the hydration shell dynamics. We propose that these conclusions can be generalized to a wide range of globular proteins. We also show that protein conformational fluctuations induce a dynamical heterogeneity within the hydration layer. We finally address the effect of confinement on hydration shell dynamics via a site-resolved analysis and connect our results to experiments via the calculation of two-dimensional infrared spectra.

Terahertz underdamped vibrational motion governs protein-ligand binding in solution

Low-frequency collective vibrational modes in proteins have been proposed as being responsible for efficiently directing biochemical reactions and biological energy transport. However, evidence of the existence of delocalized vibrational modes is scarce and proof of their involvement in biological function absent. Here we apply extremely sensitive femtosecond optical Kerr-effect spectroscopy to study the depolarized Raman spectra of lysozyme and its complex with the inhibitor triacetylchitotriose in solution. Underdamped delocalized vibrational modes in the terahertz frequency domain are identified and shown to blue-shift and strengthen upon inhibitor binding. This demonstrates that the ligand-binding coordinate in proteins is underdamped and not simply solvent-controlled as previously assumed. The presence of such underdamped delocalized modes in proteins may have significant implications for the understanding of the efficiency of ligand binding and protein-molecule interactions, and has wider implications for biochemical reactivity and biological function.

Implementing electrostatic polarization cannot fill the gap between experimental and theoretical measurements for the ultrafast fluorescence decay of myoglobin

Over the past few years, time-dependent ultrafast fluorescence spectroscopy method has been applied to the study of protein dynamics. However, observations from these experiments are in a controversy with other experimental studies. Participating of theoretical methods in this debate has not reconciled the contradiction, because the predicted initial relaxation from computer simulations is one-order faster than the ultrafast fluorescence spectroscopy experiment. In those simulations, pairwise force fields are employed, which have been shown to underestimate the roughness of the free energy landscape. Therefore, the relaxation rate of protein and water molecules under pairwise force fields is falsely exaggerated. In this work, we compared the relaxations of tryptophan/environment interaction under linear response approximation employing pairwise, polarized, and polarizable force fields. Results show that although the relaxation can be slowed down to a certain extent, the large gap between experiment and theory still cannot be filled.

Biopharmaceutical liquid formulation: a review of the science of protein stability and solubility in aqueous environments

Formulation scientists employed in the biopharmaceutical industry face the challenge of creating liquid aqueous formulations for proteins that never had evolutionary pressure to be exceptionally stable or soluble. Yet commercial products usually need a shelf life of 2 years to be economically viable. The research done in this field is dominated by physical chemists who have developed theories like preferential interaction, preferential hydration and excluded volume to explain the mechanisms for the interaction between salt, small organic molecules and proteins. This review aims to translate the research findings on protein stability and solubility produced by the physical chemists and make it accessible to formulation scientists working within the biopharmaceutical industry.

Crowding induced collective hydration of biological macromolecules over extended distances

Ultrafast two-dimensional infrared (2D-IR) spectroscopy reveals picosecond protein and hydration dynamics of crowded hen egg white lysozyme (HEWL) labeled with a metal-carbonyl vibrational probe covalently attached to a solvent accessible His residue. HEWL is systematically crowded alternatively with polyethylene glycol (PEG) or excess lysozyme in order to distinguish the chemically inert polymer from the complex electrostatic profile of the protein crowder. The results are threefold: (1) A sharp dynamical jamming-like transition is observed in the picosecond protein and hydration dynamics that is attributed to an independent-to-collective hydration transition induced by macromolecular crowding that slows the hydration dynamics up to an order of magnitude relative to bulk water. (2) The interprotein distance at which the transition occurs suggests collective hydration of proteins over distances of 30-40 Å. (3) Comparing the crowding effects of PEG400 to our previously reported experiments using glycerol exposes fundamental differences between small and macromolecular crowding agents.

Water lubricates hydrogen-bonded molecular machines

The mechanical behaviour of molecular machines differs greatly from that of their macroscopic counterparts. This applies particularly when considering concepts such as friction and lubrication, which are key to optimizing the operation of macroscopic machinery. Here, using time-resolved vibrational spectroscopy and NMR-lineshape analysis, we show that for molecular machinery consisting of hydrogen-bonded components the relative motion of the components is accelerated strongly by adding small amounts of water. The translation of a macrocycle along a thread and the rotation of a molecular wheel around an axle both accelerate significantly on the addition of water, whereas other protic liquids have much weaker or opposite effects. We tentatively assign the superior accelerating effect of water to its ability to form a three-dimensional hydrogen-bond network between the moving parts of the molecular machine. These results may indicate a more general phenomenon that helps explain the function of water as the 'lubricant of life'.

The effect of protein composition on hydration dynamics

Water dynamics at the surface of two homologous proteins with different thermal resistances is found to be unaffected by the different underlying amino-acid compositions, and when proteins are folded it responds similarly to temperature variations. Upon unfolding the water dynamics slowdown with respect to bulk decreases by a factor of two. Our findings are explained by the dominant topological perturbation induced by the protein on the water hydrogen bond dynamics.