A Molecular Dynamics Study of Lys-Trp-Lys: Structure and Dynamics in Solution Following Photoexcitation

Examining the origins of observed terahertz modes from an optically pumped atomistic model protein in aqueous solution

The microscopic origins of terahertz (THz) vibrational modes in biological systems are an active and open area of current research. Recent experiments [Phys Rev X. 8, 031061 (2018)] have revealed the presence of a pronounced mode at ∼0.3 THz in fluorophore-decorated bovine serum albumin (BSA) protein in aqueous solution under nonequilibrium conditions induced by optical pumping. This result was heuristically interpreted as a collective elastic fluctuation originating from the activation of a low-frequency phonon mode. In this work, we show that the sub-THz spectroscopic response emerges in a statistically significant manner (> 2 σ ) from such collective behavior, illustrating how photoexcitation can alter specific THz vibrational modes. We revisit the theoretical analysis with proof-of-concept molecular dynamics that introduce optical excitations into the simulations. Using information theory techniques, we show that these excitations can give rise to a multiscale response involving two optically excited chromophores (tryptophans), other amino acids in the protein, ions, and water. Our results motivate new experiments and fully nonequilibrium simulations to probe these phenomena, as well as the refinement of atomistic models of Fröhlich condensates that are fundamentally determined by nonlinear interactions in biology.

Investigation of the Effect of Substituents on Electronic and Charge Transport Properties of Benzothiazole Derivatives

A series of new benzothiazole-derived donor-acceptor-based compounds (Comp1-4) were synthesized and characterized with the objective of tuning their multifunctional properties, i.e., charge transport, electronic, and optical. All the proposed structural formulations (Comp1-4) were commensurate using FTIR, 1H NMR, 13C NMR, ESI-mass, UV-vis, and elemental analysis techniques. The effects of the electron-donating group (-CH3) and electron-withdrawing group (-NO2) on the optoelectronic and charge transfer properties were studied. The substituent effect on absorption was calculated at the TD-B3LYP/6-31+G** level in the gas and solvent phases. The effect of solvent polarity on the absorption spectra using various polar and nonpolar solvents, i.e., ethanol, acetone, DMF, and DMSO was investigated. Light was shed on the charge transport in benzothiazole compounds by calculating electron affinity, ionization potential, and reorganization energies. Furthermore, the synthesized compounds were used to prepare thin films on the FTO substrate to evaluate the charge carrier mobility and other related device parameters with the help of I-V characteristic measurements.

Tuning the Electronic and Charge Transport Properties of Schiff Base Compounds by Electron Donor and/or Acceptor Groups

Organic semiconductors have gained substantial interest as active materials in electronic devices due to their advantages over conventional semiconductors. We first designed four Schiff base compounds, then the effect of electron donor/acceptor groups (methyl/nitro) was studied on the compounds' electronic and transport nature. The absorption spectra (λ_abs) were computed by time-dependent DFT at TD-B3LYP/6-31+G** level. The effect of different solvents (ethanol, DMF, DMSO, and acetone) was investigated on the λ_abs. The substitution of the -NO₂ group to the furan moiety at the 5th position in Compound 3 leads to a red-shift in the absorption spectrum. A smaller hole reorganization energy value in Compound 3 would be beneficial to get the hole's intrinsic mobility. In contrast, a reduced-electron reorganization energy value of Compound 4 than hole may result in enhanced electron charge transfer capabilities. The reorganization energies of compounds 1 and 2 exposed balanced hole/electron transport probability. The optical, electronic, and charge transport properties at the molecular level indicate that Compound 3 is suitable for organic electronic device applications.

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.

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.

Solvation dynamics and rotational relaxation of coumarin 153 in mixed micelles of Triton X-100 and cationic gemini surfactants: effect of composition and spacer chain length of gemini surfactants

To demonstrate simultaneously how the solvation dynamics and rotational relaxation in nonionic micelles change with the composition of a gemini surfactant and how this change depends on spacer chain length of gemini surfactants.

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.

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.

Insensitivity of tryptophan fluorescence to local charge mutations

The steady state fluorescence spectral maximum (λmax) for tryptophan 140 of Staphylococcal nuclease remains virtually unchanged when nearby charged groups are removed by mutation, even though large electrostatic effects on λmax might be expected. To help understand the underlying mechanism of this curious result, we have modeled λmax with three sets of 50-ns molecular dynamics simulations in explicit water, equilibrated with excited state and with ground state charges. Semiempirical quantum mechanics and independent electrostatic analysis for the wild-type protein and four charge-altering mutants were performed on the chromophore using the coordinates from the simulations. Electrostatic contributions from the nearby charged lysines by themselves contribute 30-90 nm red shifts relative to the gas phase, but in each case, contributions from water create compensating blue shifts that bring the predicted λmax within 2 nm of the experimental value, 332 ± 0.5 nm for all five proteins. Although long-range collective interactions from ordered water make large blue shifts, crucial for determining the steady state λmax for absorption and fluorescence, such blue shifts do not contribute to the amplitude of the time dependent Stokes shift following excitation, which comes from nearby charges and only ∼6 waters tightly networked with those charges. We therefore conclude that for STNase, water and protein effects on the Stokes shift are not separable.

Ultrafast dynamics of nonequilibrium resonance energy transfer and probing globular protein flexibility of myoglobin

Protein structural plasticity is critical to many biological activities and accurate determination of its temporal and spatial fluctuations is challenging and difficult. Here, we report our extensive characterization of global flexibility of a globular heme protein of myoglobin using resonance energy transfer as a molecular ruler. With site-directed mutagenesis, we use a tryptophan scan to examine local structural fluctuations from B to H helices utilizing 10 tryptophan-heme energy transfer pairs with femtosecond resolution. We observed ultrafast resonance energy transfer dynamics by following a nearly single exponential behavior in 10-100 ps, strongly indicating that the globular structure of myoglobin is relatively rigid, with no observable static or slow dynamic conformational heterogeneity. The observation is against our molecular dynamics simulations, which show large local fluctuations and give multiple exponential energy transfer behaviors, suggesting too flexible of the global structure and thus raising a serious issue of the force fields used in simulations. Finally, these ultrafast energy transfer dynamics all occur on the similar time scales of local environmental relaxations (solvation), leading to nonexponential processes caused by energy relaxations, not structural fluctuations. Our analyses of such processes reveal an intrinsic compressed- and/or stretched-exponential behaviors and elucidate the nature of inherent nonequilibrium of ultrafast resonance energy transfer in proteins. This new concept of compressed nonequilibrium transfer dynamics should be applied to all protein studies by time-resolved Förster resonance energy transfer (FRET).

Relaxation dynamics of tryptophan in water: A UV fluorescence up-conversion and molecular dynamics study

We report on an ultrafast experimental and simulations study of the early relaxation events of photoexcited tryptophan in water. Experimentally, we used fluorescence up-conversion in both polychromatic and single wavelength detection modes in the 300-480 nm range with polarization dependence. We report on the time evolution of the Stokes shift, bandwidth, and anisotropy from tens of femtoseconds to picoseconds. These observables contain signatures of the simultaneous occurrence of intramolecular and solvent-molecule interactions, which we disentangle with the help of nonequilibrium molecular dynamics simulations. We also observe a breakdown of the linear response approximation to describe our results.

Picosecond protein dynamics: the origin of the time-dependent spectral shift in the fluorescence of the single Trp in the protein GB1

How a biological system responds to a charge shift is a challenging question directly relevant to biological function. Time-resolved fluorescence of a tryptophan residue reflects protein and solvent response to the difference in pi-electron density between the excited and the ground state. In this study we use molecular dynamics to calculate the time-dependent spectral shift (TDSS) in the fluorescence of Trp-43 in GB1 protein. A new computational method for separating solvent, protein, and fluorophore contributions to TDSS is applied to 100 nonequilibrium trajectories for GB1 in TIP3P water. The results support several nontrivial conclusions. Both longitudinal and transverse relaxation modes of bulk solvent contribute to the TDSS in proteins. All relaxation components slower than the transverse relaxation of bulk solvent have significant contributions from both protein and solvent, with a negative correlation between them. Five exponential terms in the TDSS of GB1 are well separated by their relaxation times. A 0.036 ps term is due to both solvent (60%) and protein (40%). Two exponential terms represent longitudinal (tau(L) approximately = 0.4 ps) and transverse (tau(D) approximately = 5.6 ps) relaxation modes of TIP3P water. A 131 ps term is attributable to a small change in the tertiary structure, with the alpha-helix moving 0.2 A away from the beta-strand containing Trp-43. A 2580 ps term is due to the change in the conformation of the Glu-42 side chain that brings its carboxyl group close to the positively charged end of the excited fluorophore. Interestingly, water cancels 60% of the TDSS resulting from this conformational change.

Protein hydration dynamics and molecular mechanism of coupled water-protein fluctuations

Protein surface hydration is fundamental to its structural stability and flexibility, and water-protein fluctuations are essential to biological function. Here, we report a systematic global mapping of water motions in the hydration layer around a model protein of apomyoglobin in both native and molten globule states. With site-directed mutagenesis, we use intrinsic tryptophan as a local optical probe to scan the protein surface one at a time with single-site specificity. With femtosecond resolution, we examined 16 mutants in two states and observed two types of water-network relaxation with distinct energy and time distributions. The first water motion results from the local collective hydrogen-bond network relaxation and occurs in a few picoseconds. The initial hindered motions, observed in bulk water in femtoseconds, are highly suppressed and drastically slow down due to structured water-network collectivity in the layer. The second water-network relaxation unambiguously results from the lateral cooperative rearrangements in the inner hydration shell and occurs in tens to hundreds of picoseconds. Significantly, this longtime dynamics is the coupled interfacial water-protein motions and is the direct measurement of such cooperative fluctuations. These local protein motions, although highly constrained, are necessary to assist the longtime water-network relaxation. A series of correlations of hydrating water dynamics and coupled fluctuations with local protein's chemical and structural properties were observed. These results are significant and reveal various water behaviors in the hydration layer with wide heterogeneity. We defined a solvation speed and an angular speed to quantify the water-network rigidity and local protein flexibility, respectively. We also observed that the dynamic hydration layer extends to more than 10 A. Finally, from native to molten globule states, the hydration water networks loosen up, and the protein locally becomes more flexible with larger global plasticity and partial unfolding.

Site-dependent excited-state dynamics of a fluorescent probe bound to avidin and streptavidin

The excited-state dynamics of biotin-spacer-Lucifer-Yellow (LY) constructs bound to avidin (Avi) and streptavidin (Sav) was investigated using femtosecond spectroscopy. Two different locations in the proteins, identified by molecular dynamics simulations of Sav, namely the entrance of the binding pocket and the protein surface, were probed by varying the length of the spacer. A reduction of the excited-state lifetime, stronger in Sav than in Avi, was observed with the long spacer construct. Transient absorption measurements show that this effect originates from an electron transfer quenching of LY, most probably by a nearby tryptophan residue. The local environment of the LY chromophore could be probed by measuring the time-dependent polarisation anisotropy and Stokes shift of the fluorescence. Substantial differences in both dynamics were observed. The fluorescence anisotropy decays analysed by using the wobbling-in-a-cone model reveal a much more constrained environment of the chromophore with the short spacer. Moreover, the dynamic Stokes shift is multiphasic in all cases, with a approximately 1 ps component that can be ascribed to diffusive motion of bulk-like water molecules, and with slower components with time constants varying not only with the spacer, but with the protein as well. These slow components, which depend strongly on the local environment of the probe, are ascribed to the motion of the hydration layer coupled to the conformational dynamics of the protein.

Origin of slow relaxation following photoexcitation of W7 in myoglobin and the dynamics of its hydration layer

Molecular dynamics simulations are used to calculate the time-dependent Stokes shift following photoexcitation of Trp-7 (W7) in myoglobin. In agreement with experiment, a long time (approximately 60 ps) component is observed. Since the long time Stokes shift component is absent when we repeat the calculation with protein frozen at the instant of photoexcitation, we firmly establish that protein flexibility is required to observe slow Stokes shift dynamics in this case. A transition between sub-states near the middle of a 30 ns ground-state trajectory gave us an opportunity to compare solvation dynamics in two different environments. While some of the superficial features are different, we find that the underlying dynamics are shared by the two isomers. It is necessary to look beyond a decomposition of the Stokes shift into protein and water contributions and probe the underlying dynamics of protein side groups, backbone, and water dynamics to obtain a full picture of the relaxation process. We analyze water residence times, diffusion, and reorientation dynamics in the hydration layer. We find slow components in each of these quantities and critically examine their origin and how they affect the observed Stokes shift.

Dielectric relaxation in proteins: the computational perspective

In photoexcitation and electron transfer, a new dipole or charge is introduced, and the structure is adjusted. This adjustment represents dielectric relaxation, which is the focus of this review. We concentrate on a few selected topics. We discuss linear response theory, as a unifying framework and a tool to describe non-equilibrium states. We review recent, molecular dynamics simulation studies that illustrate the calculation of dynamic and thermodynamic properties, such as Stokes shifts or reorganization free energies. We then turn to the macroscopic, continuum electrostatic view. We recall the physical definition of a dielectric constant and revisit the decomposition of the free energy into a reorganization and a static term. We review some illustrative continuum studies and discuss some difficulties that can arise with the continuum approach. In conclusion, we consider recent developments that will increase the accuracy and broaden the scope of all these methods.

Model for the Water−Amorphous Silica Interface: The Undissociated Surface

The physical and chemical properties of the amorphous silica-water interface are of crucial importance for a fundamental understanding of electrochemical and electrokinetic phenomena, and for various applications including chromatography, sensors, metal ion extraction, and the construction of micro- and nanoscale devices. A model for the undissociated amorphous silica-water interface reported here is a step toward a practical microscopic model of this important system. We have extended the popular BKS and SPC/E models for bulk silica and water to describe the hydrated, hydroxylated amorphous silica surface. The parameters of our model were determined using ab initio quantum chemical studies on small fragments. Our model will be useful in empirical potential studies, and as a starting point for ab initio molecular dynamics calculations. At this stage, we present a model for the undissociated surface. Our calculated value for the heat of immersion, 0.3 J x m(-2), falls within the range of reported experimental values of 0.2-0.8 J x m(-2). We also study the perturbation of water properties near the silica-water interface. The disordered surface is characterized by regions that are hydrophilic and hydrophobic, depending on the statistical variations in silanol group density.

Hydration dynamics and time scales of coupled water-protein fluctuations

We report experimental and theoretical studies on water and protein dynamics following photoexcitation of apomyoglobin. Using site-directed mutation and with femtosecond resolution, we experimentally observed relaxation dynamics with a biphasic distribution of time scales, 5 and 87 ps, around the site Trp7. Theoretical studies using both linear response and direct nonequilibrium molecular dynamics (MD) calculations reproduced the biphasic behavior. Further constrained MD simulations with either frozen protein or frozen water revealed the molecular mechanism of slow hydration processes and elucidated the role of protein fluctuations. Observation of slow water dynamics in MD simulations requires protein flexibility, regardless of whether the slow Stokes shift component results from the water or protein contribution. The initial dynamics in a few picoseconds represents fast local motions such as reorientations and translations of hydrating water molecules, followed by slow relaxation involving strongly coupled water-protein motions. We observed a transition from one isomeric protein configuration to another after 10 ns during our 30 ns ground-state simulation. For one isomer, the surface hydration energy dominates the slow component of the total relaxation energy. For the other isomer, the slow component is dominated by protein interactions with the chromophore. In both cases, coupled water-protein motion is shown to be necessary for observation of the slow dynamics. Such biologically important water-protein motions occur on tens of picoseconds. One significant discrepancy exists between theory and experiment, the large inertial relaxation predicted by simulations but clearly absent in experiment. Further improvements required in the theoretical model are discussed.

A Bioassay Based on the Ultrafast Response of a Reporter Molecule

The capability of using ultrafast detection technologies for a fast analysis of biomolecular reactions has been explored. As an example, the ultrafast response of tetramethylrhodamine (TMR)—labeled bovine serum albumin (BSA) as a function of different extents in proteolytic cleavage was investigated. The authors compared 4 samples of masses differing over several orders of magnitude: untreated, TMR-labeled BSA (66 kDa), TMR-labeled BSA treated with elastase (6-33 kDa) and with subtilisin (< 3 kDa), and the pure label TMR (0.4 kDa). A direct comparison with gel electrophoresis revealed that various ultrafast parameters give robust information about the progress of the proteolytic cleavage. The authors found the ratio of the transient absorption signal observed at 0 psec and 50 psec after excitation (λPump = 540 nm, λProbe = 570 nm) to be the most precise parameter for determining the cleavage. This parameter allowed determining the mass accurately within 1 sec (Z' factor of 0.83) or 600 msec (Z' factor of 0.64), measuring time per sample. This indicates that many of the known ultrafast detection technologies might be used for monitoring biochemical reactions, probably even without any labeling procedure. The authors also discuss briefly which ultrafast processes contribute to the signals and how they are affected by changes in the biomolecular environment. ( Journal of Biomolecular Screening 2007:341-350)

Exploring DNA groove water dynamics through hydrogen bond lifetime and orientational relaxation

Dynamics of water molecules in the grooves of DNA are of great interest both for practical (functionality of DNA) and fundamental (as examples of confined systems) interest. Here the authors employ atomistic molecular dynamics simulations to understand varying water dynamics at the minor and the major grooves of a 38 base-pair long DNA duplex in water. In order to understand and quantify the diversity in the nature of hydrogen bond due to many hydrogen bond donors and acceptors present in the four bases, they have undertaken study of hydrogen bond lifetime (HBLT) correlation functions of all the specific hydrogen bonds between the base atoms and water molecules. They find that the HBLT correlation functions are in general multiexponential, with the average lifetime depending significantly on the specificity and may thus be biologically relevant. The average hydrogen bond lifetime is longer in the minor groove than that in the major groove by almost a factor of 2. Analysis further shows that water hydrogen bonds with phosphate oxygen have substantially shorter lifetimes than those with the groove atoms. They also compute two different orientational time correlation functions (OTCFs) of the water molecules present at the major and the minor grooves and attempt to correlate OTCF with HBLT correlation function. The OTCFs in the minor groove exhibit three time scales, with the time constant of the slowest component one to two orders of magnitude longer than what is observed for bulk water. A slow component is also present for the major groove water but with shorter time constant. Interestingly, correlation between reformations allowed HBLT correlation function [C(HB)(t)] and the OTCF markedly deviates from each other in the grooves, indicating enhanced rigidity of water molecules in the grooves.

Protein surface hydration mapped by site-specific mutations

Water motion at protein surfaces is fundamental to protein structure, stability, dynamics, and function. By using intrinsic tryptophans as local optical probes, and with femtosecond resolution, it is possible to probe surface-water motions in the hydration layer. Here, we report our studies of local hydration dynamics at the surface of the enzyme Staphylococcus nuclease using site-specific mutations. From these studies of the WT and four related mutants, which change local charge distribution and structure, we are able to ascertain the contribution to solvation by protein side chains as relatively insignificant. We determined the time scales of hydration to be 3-5 ps and 100-150 ps. The former is the result of local librational/rotational motions of water near the surface; the latter is a direct measure of surface hydration assisted by fluctuations of the protein. Experimentally, these hydration dynamics of the WT and the four mutants are also consistent with results of the total dynamic Stokes shifts and fluorescence emission maxima and are correlated with their local charge distribution and structure. We discuss the role of protein fluctuation on the time scale of labile hydration and suggest reexamination of recent molecular dynamics simulations.