Protein surface hydration mapped by site-specific mutations

Chasing weakly-bound biological water in aqueous environment near the peptide backbone by ultrafast 2D infrared spectroscopy

There has been a long-standing debate as to how many hydrogen bonds a peptide backbone amide can form in aqueous solution. Hydrogen-bonding structural dynamics of N-ethylpropionamide (a β-peptide model) in water was examined using infrared (IR) spectroscopy. Two amide-I sub bands arise mainly from amide C=O group that forms strong H-bonds with solvent water molecules (SHB state), and minorly from that involving one weak H-bond with water (WHB state). This picture is supported by molecular dynamics simulations and ab-initio calculations. Further, thermodynamics and kinetics of the SHB and WHB species were examined mainly by chemical-exchange two-dimensional IR spectroscopy, yielding an activation energy for the SHB-to-WHB exchange of 13.25 ± 0.52 kJ mol‒1, which occurs in half picosecond at room temperature. Our results provided experimental evidence of an unstable water molecule near peptide backbone, allowing us to gain more insights into the dynamics of the protein backbone hydration.

Directed ultrafast conformational changes accompany electron transfer in a photolyase as resolved by serial crystallography

Charge-transfer reactions in proteins are important for life, such as in photolyases which repair DNA, but the role of structural dynamics remains unclear. Here, using femtosecond X-ray crystallography, we report the structural changes that take place while electrons transfer along a chain of four conserved tryptophans in the Drosophila melanogaster (6-4) photolyase. At femto- and picosecond delays, photoreduction of the flavin by the first tryptophan causes directed structural responses at a key asparagine, at a conserved salt bridge, and by rearrangements of nearby water molecules. We detect charge-induced structural changes close to the second tryptophan from 1 ps to 20 ps, identifying a nearby methionine as an active participant in the redox chain, and from 20 ps around the fourth tryptophan. The photolyase undergoes highly directed and carefully timed adaptations of its structure. This questions the validity of the linear solvent response approximation in Marcus theory and indicates that evolution has optimized fast protein fluctuations for optimal charge transfer.

Explicit Projection of Stokes Shifts onto Solvent Motion in an Aqueous Liquid and Linear Response Theory

We investigate the molecular origin of the fluorescence Stokes shift in an aqueous liquid. By examining the speed of energy change, the solvation response function is explicitly projected onto the translational and rotational motions of water molecules for both nonequilibrium relaxation and equilibrium fluctuations. Molecular dynamics simulations of a tryptophan solution show that these two processes have highly consistent dynamics, not only for the total response function but also for the decomposed components in terms of specific molecular movements. We found that the rotational mode governs the relaxation of the Stokes shift, whereas the translational mode contributes non-negligibly with slower dynamics. This consistency implies the similarity of the underlying translational and rotational movements of water molecules as the system is far away from and at equilibrium, supporting the validity of the linear response theory at the molecular level. The decomposition methodology is also applicable to a rigid solvent.

The Structural and Dynamical Properties of the Hydration of SNase Based on a Molecular Dynamics Simulation

The dynamics of protein–water fluctuations are of biological significance. Molecular dynamics simulations were performed in order to explore the hydration dynamics of staphylococcal nuclease (SNase) at different temperatures and mutation levels. A dynamical transition in hydration water (at ~210 K) can trigger larger-amplitude fluctuations of protein. The protein–water hydrogen bonds lost about 40% in the total change from 150 K to 210 K, while the Mean Square Displacement increased by little. The protein was activated when the hydration water in local had a comparable trend in making hydrogen bonds with protein– and other waters. The mutations changed the local chemical properties and the hydration exhibited a biphasic distribution, with two time scales. Hydrogen bonding relaxation governed the local protein fluctuations on the picosecond time scale, with the fastest time (24.9 ps) at the hydrophobic site and slowest time (40.4 ps) in the charged environment. The protein dynamic was related to the water’s translational diffusion via the relaxation of the protein–water’s H-bonding. The structural and dynamical properties of protein–water at the molecular level are fundamental to the physiological and functional mechanisms of SNase.

Ultrafast Fluorescence Spectroscopy via Upconversion and Its Applications in Biophysics

In this review, the experimental set-up and functional characteristics of single-wavelength and broad-band femtosecond upconversion spectrophotofluorometers developed in our laboratory are described. We discuss applications of this technique to biophysical problems, such as ultrafast fluorescence quenching and solvation dynamics of tryptophan, peptides, proteins, reduced nicotinamide adenine dinucleotide (NADH), and nucleic acids. In the tryptophan dynamics field, especially for proteins, two types of solvation dynamics on different time scales have been well explored: ~1 ps for bulk water, and tens of picoseconds for “biological water”, a term that combines effects of water and macromolecule dynamics. In addition, some proteins also show quasi-static self-quenching (QSSQ) phenomena. Interestingly, in our more recent work, we also find that similar mixtures of quenching and solvation dynamics occur for the metabolic cofactor NADH. In this review, we add a brief overview of the emerging development of fluorescent RNA aptamers and their potential application to live cell imaging, while noting how ultrafast measurement may speed their optimization.

To unravel the connection between the non-equilibrium and equilibrium solvation dynamics of tryptophan: success and failure of the linear response theory of fluorescence Stokes shift

The connections between the non-equilibrium solvation dynamics upon optical transitions and the system's equilibrium fluctuations are explored in aqueous liquid. Linear response theory correlates time-dependent fluorescence with the equilibrium time correlation functions. In the previous work [T. Li, J. Chem. Theory Comput., 2017, 13, 1867], Stokes shift was explicitly decomposed into the contributions of various order time correlation functions on the excited state surface. Gaussian fluctuations of the solute-solvent interactions validate linear response theory. Correspondingly, the deviation of the Gaussian statistics causes the inefficiency of linear response evaluation. The above mechanism is thoroughly tested in this study. By employing molecular simulations, multiple non-equilibrium processes, not necessarily initiated from the ground state equilibrium minimum, were examined for tryptophan. Both the success and failure of linear response theory are found for this simple system and the mechanism is analyzed. These observations, assisted by the width dynamics, the initial state linear response approach, and the variation of the solvation structures, integrally verify the virtue of the excited state Gaussian statistics on the dynamics of Stokes shift.

Linear Response Theory for Decomposition Energies of Stokes Shift in Proteins

In aqueous solution, fluorescence Stokes shift experiments monitor the relaxation of the solute-solvent interactions upon photon excitation of the solute chromophore. Linear response (LR) theory expects the identical dynamics between the Stokes shift and the system's spontaneous fluctuations. Whether this identity guarantees similar dynamics between the nonequilibrated and equilibrium processes for the decomposition energy of the Stokes shift is the main focus of this study. In our previous work [Li, T. J. Chem. Theory Comput. 2017, 13, 1867-1873], Stokes shift is properly correlated with various order time-correlation functions. As a continuation, its decomposition energy from the subsystem is further represented as the full summation of all of the cross-time correlation functions between the decomposition energy and the total solute-solvent interactions. Gaussian statistics of the total solute-solvent interactions ensure the same decay rates among the odd orders not only for the time-correlation functions but also for the cross-time correlation functions, validating the LR of the Stokes shift and the decompositions, respectively. The above mechanism is verified by molecular dynamics simulations in the protein Staphylococcus nuclease and is robust even as the decomposed energy associated with an individual residue exhibits typical non-Gaussian properties. Further examinations reveal the consistent molecular motions for a specific residue over the nonequilibrium and equilibrium processes, which are responsible for the nonequilibrium dynamics of the associated decomposed energy. Our results show the appropriateness of LR on finer molecular scales.

Femtosecond Fluorescence Spectra of NADH in Solution: Ultrafast Solvation Dynamics

The ultrafast solvation dynamics of reduced nicotinamide adenine dinucleotide (NADH) free in solution has been investigated, using both a femtosecond upconversion spectrophotofluorometer and a picosecond time-correlated single-photon counting (TCSPC) apparatus. The familiar time constant of solvent relaxation originating in "bulk water" was found to be ∼1.4 ps, revealing ultrafast solvent reorientation upon excitation. We also found a slower spectral relaxation process with an apparent time of 27 ps, suggesting there could either be dissociable "biological water" hydration sites on the surface of NADH or internal dielectric rearrangements of the flexible solvated molecule on that timescale. In contrast, the femtosecond fluorescence anisotropy measurement revealed that rotational diffusion happened on two different timescales (3.6 ps (local) and 141 ps (tumbling)); thus, any dielectric rearrangement scenario for the 27 ps relaxation must occur without significant chromophore oscillator rotation. The coexistence of quasi-static self quenching (QSSQ) with the slower relaxation is also discussed.

Dipole-dipole interactions between tryptophan side chains and hydration water molecules dominate the observed dynamic stokes shift of lysozyme

The fluorescence intensity of tryptophan residues in hen egg-white lysozyme was measured up to 500 ps after the excitation by irradiation pulses at 290 nm. From the time-dependent variation of fluorescence intensity in a wavelength range of 320-370 nm, the energy relaxation in the dynamic Stokes shift was reconstructed as the temporal variation in wavenumber of the estimated fluorescence maximum. The relaxation was approximated by two exponential curves with decay constants of 1.2 and 26.7 ps. To interpret the relaxation, a molecular dynamics simulation of 75 ns was conducted for lysozyme immersed in a water box. From the simulation, the energy relaxation in the electrostatic interactions of each tryptophan residue was evaluated by using a scheme derived from the linear response theory. Dipole-dipole interactions between each of the Trp62 and Trp123 residues and hydration water molecules displayed an energy relaxation similar to that experimentally observed regarding time constants and magnitudes. The side chains of these residues were partly or fully exposed to the solvent. In addition, by inspecting the variation in dipole moments of the hydration water molecules around lysozyme, it was suggested that the observed relaxation could be attributed to the orientational relaxation of hydration water molecules participating in the hydrogen-bond network formed around each of the two tryptophan residues.

Surface Topography Effects of Globular Biomolecules on Hydration Water

Hydration water serves as a microscopic manifestation of structural stability and functions of biomolecules. To develop bio-nanomaterials in applications, it is important to study how the surface topography and heterogeneity of biomolecules result in their diversity of the hydration dynamics and energetics. We here performed molecular dynamics simulations combined with the steered molecular dynamics and umbrella sampling to investigate the dynamics and escape process associated with the free energy change of water molecules close to a globular biomolecule, i.e., hemoglobin (Hb) and G-quadruplex DNA (GDNA). The residence time, power of long-time tail, and dipole relaxation time were found to display drastic changes within the averaged hydration shell of 3.0-5.0 Å. Compared with bulk water, in the inner hydration shell, the water dipole moment displays a slower relaxation process and is more oriented toward GDNA than toward Hb, forming a hedgehog-like structure when it surrounds GDNA. In particular, a spine water structure is observed in the GDNA narrow groove. The water isotope effect not only prolongs the dynamic time scales of libration motion in the inner hydration shell and the dipole relaxation processes in the bulk but also strengthens the DNA spine water structure. The potential of the mean force profile reflects the integrity of the hydration shell structure and enables us to obtain detailed insights into the structures formed by water, such as the caged H-bond network and the edge bridge structures; it also reveals that local hydration shell free energy (LHSFE) depends on H-bond rupture processes and ranges from 0.2 to 4.2 kcal/mol. Our results demonstrate that the surface topography of a biomolecule influences the integrity of the hydration shell structure and LHSFE. Our studies are able to identify various further applications in the areas of microfluid devices and nano-dewetting on bioinspired surfaces.

Protein Solvent Shell Structure Provides Rapid Analysis of Hydration Dynamics

The solvation layer surrounding a protein is clearly an intrinsic part of protein structure-dynamics-function, and our understanding of how the hydration dynamics influences protein function is emerging. We have recently reported simulations indicating a correlation between regional hydration dynamics and the structure of the solvation layer around different regions of the enzyme Candida antarctica lipase B, wherein the radial distribution function (RDF) was used to calculate the pairwise entropy, providing a link between dynamics (diffusion) and thermodynamics (excess entropy) known as Rosenfeld scaling. Regions with higher RDF values/peaks in the hydration layer (the first peak, within 6 Å of the protein surface) have faster diffusion in the hydration layer. The finding thus hinted at a handle for rapid evaluation of hydration dynamics at different regions on the protein surface in molecular dynamics simulations. Such an approach may move the analysis of hydration dynamics from a specialized venture to routine analysis, enabling an informatics approach to evaluate the role of hydration dynamics in biomolecular function. This paper first confirms that the correlation between regional diffusive dynamics and hydration layer structure (via water center of mass around protein side-chain atom RDF) is observed as a general relationship across a set of proteins. Second, it seeks to devise an approach for rapid analysis of hydration dynamics, determining the minimum amount of information and computational effort required to get a reliable value of hydration dynamics from structural data in MD simulations based on the protein-water RDF. A linear regression model using the integral of the hydration layer in the water-protein RDF was found to provide statistically equivalent apparent diffusion coefficients at the 95% confidence level for a set of 92 regions within five different proteins. In summary, RDF analysis of 10 ns of data after simulation convergence is sufficient to accurately map regions of fast and slow hydration dynamics around a protein surface. Additionally, it is anticipated that a quick look at protein-water RDFs, comparing peak heights, will be useful to provide a qualitative ranking of regions of faster and slower hydration dynamics at the protein surface for rapid analysis when investigating the role of solvent dynamics in protein function.

Theoretical Insights on Nonlinear Response Theory of Fluorescence Spectroscopy in Liquids

Correlations between the nonequilibrium solvation dynamics upon the photon excitation of the chromophore and a system's equilibrium fluctuations are deeply studied. As the linear response of the solvent has been linked with Gaussian statistics of the energy fluctuations in the literature, we specifically explore the cases beyond the regime of the linear response theory due to deviation from Gaussian fluctuations. As a continuation of our previous work, an analytical formalism is presented to project the energy shift with various order moments, where the non-Gaussian statistics arise from the overlap of the energy basins on the perturbed potential energy surface. It is shown that the nonequilibrium dynamics still correlate with the spontaneous regressions at equilibrium and are controlled by the decay rates of those higher order components with the prevailing contributions to the energy shift. Molecular dynamics simulations were performed in the protein Staphylococcus nuclease, in which even the dynamics of the high order moments are available. The results further verify the above relationship. Our scheme is used to evaluate Stokes shift using the information on non-Gaussian statistics at equilibrium, thus presenting a broad picture on the correlation between the nonequilibrium process and equilibrium properties in liquids.

Distinguishing dynamical features of water inside protein hydration layer: Distribution reveals what is hidden behind the average

Since the pioneering works of Pethig, Grant, and Wüthrich on a protein hydration layer, many studies have been devoted to find out if there are any "general and universal" characteristic features that can distinguish water molecules inside the protein hydration layer from bulk. Given that the surface itself varies from protein to protein, and that each surface facing the water is heterogeneous, search for universal features has been elusive. Here, we perform an atomistic molecular dynamics simulation in order to propose and demonstrate that such defining characteristics can emerge if we look not at average properties but the distribution of relaxation times. We present results of calculations of distributions of residence times and rotational relaxation times for four different protein-water systems and compare them with the same quantities in the bulk. The distributions in the hydration layer are unusually broad and log-normal in nature due to the simultaneous presence of peptide backbones that form weak hydrogen bonds, hydrophobic amino acid side chains that form no hydrogen bond, and charged polar groups that form a strong hydrogen bond with the surrounding water molecules. The broad distribution is responsible for the non-exponential dielectric response and also agrees with large specific heat of the hydration water. Our calculations reveal that while the average time constant is just about 2-3 times larger than that of bulk water, it provides a poor representation of the real behaviour. In particular, the average leads to the erroneous conclusion that water in the hydration layer is bulk-like. However, the observed and calculated lower value of static dielectric constant of hydration layer remained difficult to reconcile with the broad distribution observed in dynamical properties. We offer a plausible explanation of these unique properties.

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.

Slow Dynamics around a Protein and Its Coupling to Solvent

Solvent is essential for protein dynamics and function, but its role in regulating the dynamics remains debated. Here, we employ saturation transfer electron spin resonance (ST-ESR) to explore the issue and characterize the dynamics on a longer (from μs to s) time scale than has been extensively studied. We first demonstrate the reliability of ST-ESR by showing that the dynamical changeovers revealed in the spectra agree to liquid-liquid transition (LLT) in the state diagram of the glycerol/water system. Then, we utilize ST-ESR with four different probes to systematically map out the variation in local (site-specific) dynamics around a protein surface at subfreezing temperatures (180-240 K) in 10 mol % glycerol/water mixtures. At highly exposed sites, protein and solvent dynamics are coupled, whereas they deviate from each other when temperature is greater than LLT temperature (∼190 K) of the solvent. At less exposed sites, protein however exhibits a dynamic, which is distinct from the bulk solvent, throughout the temperature range studied. Dominant dynamic components are thus revealed, showing that (from low to high temperatures) the overall structural fluctuation, rotamer dynamics, and internal side-chain dynamics, in turn, dominate the temperature dependence of spin-label motions. The structural fluctuation component is relatively slow, collective, and independent of protein structural segments, which is thus inferred to a fundamental dynamic component intrinsic to protein. This study corroborates that bulk solvent plasticizes protein and facilitates rather than slaves protein dynamics.

Deciphering complex dynamics of water counteraction around secondary structural elements of allosteric protein complex: Case study of SAP-SLAM system in signal transduction cascade

The first hydration shell of a protein exhibits heterogeneous behavior owing to several attributes, majorly local polarity and structural flexibility as revealed by solvation dynamics of secondary structural elements. We attempt to recognize the change in complex water counteraction generated due to substantial alteration in flexibility during protein complex formation. The investigation is carried out with the signaling lymphocytic activation molecule (SLAM) family of receptors, expressed by an array of immune cells, and interacting with SLAM-associated protein (SAP), composed of one SH2 domain. All atom molecular dynamics simulations are employed to the aqueous solutions of free SAP and SLAM-peptide bound SAP. We observed that water dynamics around different secondary structural elements became highly affected as well as nicely correlated with the SLAM-peptide induced change in structural rigidity obtained by thermodynamic quantification. A few instances of contradictory dynamic features of water to the change in structural flexibility are explained by means of occluded polar residues by the peptide. For βD, EFloop, and BGloop, both structural flexibility and solvent accessibility of the residues confirm the obvious contribution. Most importantly, we have quantified enhanced restriction in water dynamics around the second Fyn-binding site of the SAP due to SAP-SLAM complexation, even prior to the presence of Fyn. This observation leads to a novel argument that SLAM induced more restricted water molecules could offer more water entropic contribution during the subsequent Fyn binding and provide enhanced stability to the SAP-Fyn complex in the signaling cascade. Finally, SLAM induced water counteraction around the second binding site of the SAP sheds light on the allosteric property of the SAP, which becomes an integral part of the underlying signal transduction mechanism.

Spatially Heterogeneous Surface Water Diffusivity around Structured Protein Surfaces at Equilibrium

Hydration water on the surface of a protein is thought to mediate the thermodynamics of protein-ligand interactions. For hydration water to play a role beyond modulating global protein solubility or stability, the thermodynamic properties of hydration water must reflect on the properties of the heterogeneous protein surface and thus spatially vary over the protein surface. A potent read-out of local variations in thermodynamic properties of hydration water is its equilibrium dynamics spanning picosecond to nanosecond time scales. In this study, we employ Overhauser dynamic nuclear polarization (ODNP) to probe the equilibrium hydration water dynamics at select sites on the surface of Chemotaxis Y (CheY) in dilute solution. ODNP reports on site-specific hydration water dynamics within 5-10 Å of a label tethered to the biomolecular surface on two separate time scales of motion, corresponding to diffusive water (DW) and protein-water coupled motions, referred to as bound water (BW). We find DW dynamics to be highly heterogeneous across the surface of CheY. We identify a significant correlation between DW dynamics and the local hydropathy of the CheY protein surface, empirically determined by molecular dynamics (MD) simulations, and find the more hydrophobic sites to be hydrated with slower diffusing water. Furthermore, we compare the hydration water dynamics on different polypeptides and liposome surfaces and find the DW dynamics on globular proteins to be significantly more heterogeneous than on intrinsically disordered proteins (IDPs), peptides, and liposomes. The heterogeneity in the hydration water dynamics suggests that structured proteins have the capacity to encode information into the surrounding hydration shell.

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.

Efficient Criterion To Evaluate Linear Response Theory in Optical Transitions

The role of the Gaussian statistics on the solvation dynamics upon the photon excitation of the chromophore is deeply explored. The linear response theory for the fluorescence Stokes shift is investigated. An analytical formulism is presented to recast Stokes shift into the contributions of the equilibrium time correlation functions of the solute-solvent interactions on the excited-state surface, and the latter is further reformed and depicted by the time relaxation of the moment. As the first application of the formulism in the molecular dynamics simulations, it is verified that the efficiency of the linear response theory relies on the Gaussian characteristics of the dominant moments in terms of the Stokes shift, which is identified by the same relaxation dynamics between those moments and the linear order one. The comparisons between the above observations on the linearity of Stokes shift and the explanations in the literature are discussed. The key finding is the development of explicit criterion to measure the appropriateness of applying linear response theory.

Observation of the Global Dynamic Collectivity of a Hydration Shell around Apomyoglobin

Protein surface hydration is critical to the protein's structural properties and biological activities. However, it is still unknown whether the hydration shell is intrinsically connected and how its fluctuations dynamically interact with protein motion. Here, by selecting five site-specific locations with distinctly different environments around the surface of apomyoglobin, we used a tryptophan scan with femtosecond fluorescence spectroscopy and simultaneously detected hydration water dynamics and tryptophan side-chain relaxations with temperature dependence. We observed two types of relaxations for both interfacial hydration water and the tryptophan side chain. The former is always faster than the latter, and both motions show direct linear correlations with temperature changes, indicating one origin of their motions and hydration water driving of side-chain fluctuations. Significantly, we found the relaxation energy barriers are uniform across the entire protein surface, all less than 20 kJ/mol, strongly suggesting highly extended cooperative water networks and the nature of global dynamic collectivity of the entire hydration shell.

The opposite effects of sodium and potassium cations on water dynamics

Water rotational dynamics in NaSCN and KSCN solutions at a series of concentrations are investigated using femtosecond infrared spectroscopy and theory. Femtosecond infrared measurements, consistent with previous NMR observations, detect that sodium slows down while potassium accelerates the water O-H bond rotation. Results of reported neutron scattering measurements, on the other hand, suggested that these two cations have similar structure-breaking effects on water, and therefore should both accelerate water rotation through the presumably dominating large-amplitude angular jump component. To explain this discrepancy, theoretical studies with both classical and ab initio models were carried out, which indicate that both ions indeed accelerate the large-amplitude angular jump rotation of the water molecules, while the observed cation specific effect originates from the non-negligible opposite impact of the sodium and potassium cations on the diffusive rotation of water molecules.

Molecular Origin of Ultrafast Water-Protein Coupled Interactions

The fluctuations of hydration water and the protein are coupled together at the protein surface and often such water-protein dynamic interactions are controlled presumably by hydration water motions. However, direct evidence is scarce and it requires measuring the dynamics of hydration water and protein side chain simultaneously. Here, we use a unique protein with a single tryptophan to directly probe interfacial water and related side chain relaxations with temperature dependence. With systematic mutations to change local chemical identity and structural flexibility, we found that the side chain relaxations are always slower than hydration water motions and the two dynamic processes are linearly correlated with the same energy barriers, indicating the same origin of both relaxations. The charge mutations change the rates of hydration water relaxations but not the relaxation barriers. These results convincingly show that the water-protein relaxations are strongly coupled and the hydration water molecules govern such fluctuations on the picosecond time scales.

Dynamics and mechanism of ultrafast water-protein interactions

Protein hydration is essential to its structure, dynamics, and function, but water-protein interactions have not been directly observed in real time at physiological temperature to our awareness. By using a tryptophan scan with femtosecond spectroscopy, we simultaneously measured the hydration water dynamics and protein side-chain motions with temperature dependence. We observed the heterogeneous hydration dynamics around the global protein surface with two types of coupled motions, collective water/side-chain reorientation in a few picoseconds and cooperative water/side-chain restructuring in tens of picoseconds. The ultrafast dynamics in hundreds of femtoseconds is from the outer-layer, bulk-type mobile water molecules in the hydration shell. We also found that the hydration water dynamics are always faster than protein side-chain relaxations but with the same energy barriers, indicating hydration shell fluctuations driving protein side-chain motions on the picosecond time scales and thus elucidating their ultimate relationship.

Proteins Take up Water Before Unfolding

Proteins perform specific biological functions that strongly depend on their three-dimensional structure. This three-dimensional structure, i.e. the way the protein folds, is strongly determined by the interaction between the protein and the water solvent. We study the dynamics of water in aqueous solutions of several globular proteins at different degrees of urea-induced unfolding, using polarization-resolved femtosecond infrared spectroscopy. We observe that a fraction of the water molecules is strongly slowed down by their interaction with the protein surface. By monitoring the slow water fraction we can directly probe the amount of water-exposed protein surface. We find that at mild denaturing conditions, the water-exposed surface increases by almost 50%, while the secondary structure is still intact. This finding indicates that protein unfolding starts with the protein structure becoming less tight, thereby allowing water to enter.

Determination of Protein Surface Hydration by Systematic Charge Mutations

Protein surface hydration is critical to its structural stability, flexibility, dynamics, and function. Recent observations of surface solvation on picosecond time scales have evoked debate on the origin of such relatively slow motions, from hydration water or protein charged side chains, especially with molecular dynamics simulations. Here we used a unique nuclease with a single tryptophan as a local probe and systematically mutated three neighboring charged residues to differentiate the contributions from hydration water and charged side chains. By various mutations of one, two, and all three charged residues, we observed slight increases in the total tryptophan Stokes shifts with fewer neighboring charged residue(s) and found insensitivity of charged side chains to the relaxation patterns. The dynamics is correlated with hydration water relaxation with the slowest time in a dense charged environment and the fastest time at a hydrophobic site. On such picosecond time scales, the protein surface motion is restricted. The total Stokes shifts are dominantly from hydration water relaxation and the slow dynamics is from water-driven relaxation, coupled to local protein fluctuations.

Protein dynamics: from rattling in a cage to structural relaxation

We present an overview of protein dynamics based mostly on results of neutron scattering, dielectric relaxation spectroscopy and molecular dynamics simulations. We identify several major classes of protein motions on the time scale from faster than picoseconds to several microseconds, and discuss the coupling of these processes to solvent dynamics. Our analysis suggests that the microsecond backbone relaxation process might be the main structural relaxation of the protein that defines its glass transition temperature, while faster processes present some localized secondary relaxations. Based on the overview, we formulate a general picture of protein dynamics and discuss the challenges in this field.

Dynamics and mechanisms of DNA repair by photolyase

Photolyases, a class of flavoproteins, use blue light to repair two types of ultraviolet-induced DNA damage, a cyclobutane pyrimidine dimer (CPD) and a pyrimidine-pyrimidone (6-4) photoproduct (6-4PP). In this perspective, we review the recent progress in the repair dynamics and mechanisms of both types of DNA restoration by photolyases. We first report the spectroscopic characterization of flavin in various redox states and the active-site solvation dynamics in photolyases. We then systematically summarize the detailed repair dynamics of damaged DNA by photolyases and a biomimetic system through resolving all elementary steps on ultrafast timescales, including multiple intermolecular electron- and proton-transfer reactions and bond-breaking and -making processes. We determined the unique electron tunneling pathways, identified the key functional residues and revealed the molecular origin of high repair efficiency, and thus elucidate the molecular mechanisms and repair photocycles at the most fundamental level. We finally conclude that the active sites of photolyases, unlike the aqueous solution for the biomimetic system, provide a unique electrostatic environment and local flexibility and thus a dedicated synergy for all elementary dynamics to maximize the repair efficiency. This repair photomachine is the first enzyme that the entire functional evolution is completely mapped out in real time.

Picosecond fluorescence dynamics of tryptophan and 5-fluorotryptophan in monellin: slow water-protein relaxation unmasked

Time dependent fluorescence Stokes (emission wavelength) shifts (TDFSS) from tryptophan (Trp) following sub-picosecond excitation are increasingly used to investigate protein dynamics, most recently enabling active research interest into water dynamics near the surface of proteins. Unlike many fluorescence probes, both the efficiency and the wavelength of Trp fluorescence in proteins are highly sensitive to microenvironment, and Stokes shifts can be dominated by the well-known heterogeneous nature of protein structure, leading to what we call pseudo-TDFSS: shifts that arise from differential decay rates of subpopulations. Here we emphasize a novel, general method that obviates pseudo-TDFSS by replacing Trp by 5-fluorotryptophan (5Ftrp), a fluorescent analogue with higher ionization potential and greatly suppressed electron-transfer quenching. 5FTrp slows and suppresses pseudo-TDFSS, thereby providing a clearer view of genuine relaxation caused by solvent and protein response. This procedure is applied to the sweet-tasting protein monellin which has uniquely been the subject of ultrafast studies in two different laboratories (Peon, J.; et al. Proc. Natl. Acad. Sci. U.S.A. 2002, 99, 10964; Xu, J.; et al. J. Am. Chem. Soc. 2006, 128, 1214) that led to disparate interpretations of a 20 ps transient. They differed because of the pseudo-TDFSS present. The current study exploiting special properties of 5FTrp strongly supports the conclusion that both lifetime heterogeneity-based TDFSS and environment relaxation-based TDFSS are present in monellin and 5FTrp-monellin. The original experiments on monellin were most likely dominated by pseudo-TDFSS, whereas, in the present investigation of 5FTrp-monellin, the TDFSS is dominated by relaxation and any residual pseudo-TDFSS is overwhelmed and/or slowed to irrelevance.

Ultrafast water dynamics at the interface of the polymerase-DNA binding complex

DNA polymerases slide on DNA during replication, and the interface must be mobile for various conformational changes. The role of lubricant interfacial water is not understood. In this report, we systematically characterize the water dynamics at the interface and in the active site of a tight binding polymerase (pol β) in its binary complex and ternary state using tryptophan as a local optical probe. Using femtosecond spectroscopy, we observed that upon DNA recognition the surface hydration water is significantly confined and becomes bound water at the interface, but the dynamics are still ultrafast and occur on the picosecond time scale. These interfacial water molecules are not trapped but are mobile in the heterogeneous binding nanospace. Combining our findings with our previous observation of ultrafast water motions at the interface of a loose binding polymerase (Dpo4), we conclude that the binding interface is dynamic and the water molecules in various binding clefts, channels, and caves are mobile and even fluid with different levels of mobility for loose or tight binding polymerases. Such a dynamic interface should be general to all DNA polymerase complexes to ensure the biological function of DNA synthesis.

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.

Picosecond-resolved fluorescent probes at functionally distinct tryptophans within a thermophilic alcohol dehydrogenase: relationship of temperature-dependent changes in fluorescence to catalysis

Two single-tryptophan variants were generated in a thermophilic alcohol dehydrogenase with the goal of correlating temperature-dependent changes in local fluorescence with the previously demonstrated catalytic break at ca. 30 °C (Kohen et al., Nature 1999, 399, 496). One tryptophan variant, W87in, resides at the active site within van der Waals contact of bound alcohol substrate; the other variant, W167in, is a remote-site surface reporter located >25 Å from the active site. Picosecond-resolved fluorescence measurements were used to analyze fluorescence lifetimes, time-dependent Stokes shifts, and the extent of collisional quenching at Trp87 and Trp167 as a function of temperature. A subnanosecond fluorescence decay rate constant has been detected for W87in that is ascribed to the proximity of the active site Zn(2+) and shows a break in behavior at 30 °C. For the remainder of the reported lifetime measurements, there is no detectable break between 10 and 50 °C, in contrast with previously reported hydrogen/deuterium exchange experiments that revealed a temperature-dependent break analogous to catalysis (Liang et al., Proc. Natl. Acad. Sci. U.S.A. 2004, 101, 9556). We conclude that the motions that lead to the rigidification of ht-ADH below 30 °C are likely to be dominated by global processes slower than the picosecond to nanosecond motions measured herein. In the case of collisional quenching of fluorescence by acrylamide, W87in and W167in behave in a similar manner that resembles free tryptophan in water. Stokes shift measurements, by contrast, show distinctive behaviors in which the active-site tryptophan relaxation is highly temperature-dependent, whereas the solvent-exposed tryptophan's dynamics are temperature-independent. These data are concluded to reflect a significantly constrained environment surrounding the active site Trp87 that both increases the magnitude of the Stokes shift and its temperature-dependence. The results are discussed in the context of spatially distinct differences in enthalpic barriers for protein conformational sampling that may be related to catalysis.

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.

Charge invariant protein-water relaxation in GB1 via ultrafast tryptophan fluorescence

The protein–water interface is a critical determinant of protein structure and function, yet the precise nature of dynamics in this complex system remains elusive. Tryptophan fluorescence has become the probe of choice for such dynamics on the picosecond time scale (especially via fluorescence “upconversion”). In the absence of ultrafast (“quasi-static”) quenching from nearby groups, the TDFSS (time-dependent fluorescence Stokes shift) for exposed Trp directly reports on dipolar relaxation near the interface (both water and polypeptide). The small protein GB1 contains a single Trp (W43) of this type, and its structure is refractory to pH above 3. Thus, it can be used to examine the dependence of dipolar relaxation upon charge reconfiguration with titration. Somewhat surprisingly, the dipolar dynamics in the 100 fs to 100 ps range were unchanged with pH, although nanosecond yield, rates, and access all changed. These results were rationalized with the help of molecular dynamics (including QM-MM) simulations that reveal a balancing, sometimes even countervailing influence of protein and water dipoles. Interestingly, these simulations also showed the dominant influence of water molecules which are associated with the protein interface for up to 30 ps yet free to rotate at approximately “bulk” water rates.

Dynamic insight into protein structure utilizing red edge excitation shift

Proteins are considered the workhorses in the cellular machinery. They are often organized in a highly ordered conformation in the crowded cellular environment. These conformations display characteristic dynamics over a range of time scales. An emerging consensus is that protein function is critically dependent on its dynamics. The subtle interplay between structure and dynamics is a hallmark of protein organization and is essential for its function. Depending on the environmental context, proteins can adopt a range of conformations such as native, molten globule, unfolded (denatured), and misfolded states. Although protein crystallography is a well established technique, it is not always possible to characterize various protein conformations by X-ray crystallography due to transient nature of these states. Even in cases where structural characterization is possible, the information obtained lacks dynamic component, which is needed to understand protein function. In this overall scenario, approaches that reveal information on protein dynamics are much appreciated. Dynamics of confined water has interesting implications in protein folding. Interfacial hydration combines the motion of water molecules with the slow moving protein molecules. The red edge excitation shift (REES) approach becomes relevant in this context. REES is defined as the shift in the wavelength of maximum fluorescence emission toward higher wavelengths, caused by a shift in the excitation wavelength toward the red edge of absorption spectrum. REES arises due to slow rates (relative to fluorescence lifetime) of solvent relaxation (reorientation) around an excited state fluorophore in organized assemblies such as proteins. Consequently, REES depends on the environment-induced motional restriction imposed on the solvent molecules in the immediate vicinity of the fluorophore. In the case of a protein, the confined water in the protein creates a dipolar field that acts as the solvent for a fluorophore in the protein. In this Account, we focus on REES to monitor organization and dynamics of soluble and membrane proteins utilizing intrinsic protein fluorescence. We discuss here the application of REES in various conformations of proteins. While application of REES to proteins in native conformation has been in use for a long time, our work highlights the potential of this approach in case of molten globule and denatured conformations. For example, we have demonstrated the presence of residual structure, that could not be detected using other methods, by REES of denatured spectrin. Given the functional relevance of such residual structures, these results are very far reaching. We discuss here the application of REES to molten globule conformation and to the green fluorescent protein (GFP). The case of GFP is particularly interesting since the dipolar field in this case is provided by the protein matrix itself and not confined water. We envision that future applications of REES in proteins will involve generating a dynamic hydration map of the protein, which would allow us to explore protein function in terms of local dynamics and hydration.