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

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

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

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.

Prolonged Association between Water Molecules under Hydrophobic Nanoconfinement

We present an investigation of the dynamics of water confined among rigid carbon rods and between parallel graphene sheets with molecular dynamics simulations. Diffusion coefficients, activation energy of diffusion, and residence-time correlation functions as a function of confinement geometry reveal a retardation of water dynamics under hydrophobic confinement compared to bulk water. In fact, water under various confinements possesses longer associations with its neighbors and exhibits diffusion dynamics characteristic of a lower temperature. Analysis of the residence-time correlation functions reveals long and short residence times, which we relate to the diffusion coefficient and activation energy of diffusion, respectively. Additional investigations reveal how the level of confining surface hydrophobicity affects water dynamics, further broadening our understanding of water diffusion inside diverse media. Overall, this study sheds light on the physical origin of retarded water dynamics under hydrophobic confinement and the close relationship between residence times and diffusion behavior.

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.

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.

Unraveling the Role of Monoolein in Fluidity and Dynamical Response of a Mixed Cationic Lipid Bilayer

The maintenance of cell membrane fluidity is of critical importance for various cellular functions. At lower temperatures when membrane fluidity decreases, plants and cyanobacteria react by introducing unsaturation in the lipids, so that the membranes return to a more fluidic state. To probe how introduction of unsaturation leads to reduced membrane fluidity, a model cationic lipid dioctadecyldimethylammonium bromide (DODAB) has been chosen, and the effects of an unsaturated lipid monoolein (MO) on the structural dynamics and phase behavior of DODAB have been monitored by quasielastic neutron scattering and time-resolved fluorescence measurements. In the coagel phase, fluidity of the lipid bilayer increases significantly in the presence of MO relative to pure DODAB vesicles and becomes manifest in significantly enhanced dynamics of the constituent lipids along with faster hydration and orientational relaxation dynamics of a fluorophore. On the contrary, MO restricts both lateral and internal motions of the lipid molecules in the fluid phase (>330 K), which is consistent with relatively slow hydration and orientational relaxation dynamics of the fluorophore embedded in the mixed lipid bilayer. The present study illustrates how incorporation of an unsaturated lipid at lower temperatures (below the phase transition) assists the model lipid (DODAB) in regulating fluidity via enhancement of dynamics of the constituent lipids.

Probing relaxation dynamics of a cationic lipid based non-viral carrier: a time-resolved fluorescence study

Lipid vesicles composed of cationic dioctadecyldimethylammonium bromide (DODAB) and neutral 1-monooleoyl-rac-glycerol (MO) are promising non-viral carriers of nucleic acids for delivery into cells. Among them, higher cell transfection efficiency was displayed by DODAB-rich vesicles than those enriched with MO. Structural relaxation of these mixed lipid vesicles plays a key role in their cell transfection efficiency because structural organization of the DODAB-rich vesicles are different from that of the MO-rich vesicles. Polarization-gated anisotropy in conjunction with picosecond resolved emission transients of a novel fluorophore 6-acetyl-(2-((4-hydroxycyclohexyl)(methyl)amino)naphthalene) (ACYMAN) has been employed to probe relaxation dynamics in pure DODAB vesicles, and in mixed vesicles of DODAB with varying content of MO. Both orientational relaxation of ACYMAN and relaxation dynamics of its local environment are retarded significantly in mixed lipid vesicles with increasing MO content, from a mole fraction (χ_MO) of 0.2 to that of 0.8 which is consistent with increased rigidity of the MO-rich (χ_MO > 0.5) vesicles relative to the DODAB-rich (χ_MO < 0.5) vesicles. Therefore, higher structural rigidity of the MO-rich vesicles (χ_MO > 0.5) gives rise to their lower cell transfection efficiency than the more flexible DODAB-rich (χ_MO < 0.5) vesicles as observed in previous in vivo studies (Biochim. Biophys. Acta, Biomembr., 2014, 1838, 2555–2567).

Lipid vesicles composed of cationic dioctadecyldimethylammonium bromide (DODAB) and neutral 1-monooleoyl-rac-glycerol (MO) are promising non-viral carriers of nucleic acids for delivery into cells.

Fundamental limitations of the time-dependent Stokes shift for investigating protein hydration dynamics

In the TDSS measured in protein systems, large protein contributions fully obscure hydration dynamics.

Changes in protein hydration dynamics by encapsulation or crowding of ubiquitin: strong correlation between time-dependent Stokes shift and intermolecular nuclear Overhauser effect

The local changes in protein hydration dynamics upon encapsulation of the protein or macromolecular crowding are essential to understand protein function in cellular environments. We were able to obtain a spatially-resolved picture of the influence of confinement and crowding on the hydration dynamics of the protein ubiquitin by analyzing the time-dependent Stokes shift (TDSS), as well as the intermolecular Nuclear Overhauser Effect (NOE) at different sites of the protein by large-scale computer simulation of single and multiple proteins in water and confined in reverse micelles. Besides high advanced space resolved information on hydration dynamics we found a strong correlation of the change in NOE upon crowding or encapsulation and the change in the integral TDSS relaxation times in all investigated systems relative to the signals in a diluted protein solution.

Changes in local protein hydration dynamics caused by encapsulation or crowding are reflected in the TDSS and the intermolecular NOE alike.

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.

On the coupling between the dynamics of protein and water

The solvation entropy of flexible protein regions is higher than that of rigid regions and contributes differently to the overall thermodynamic stability.

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.

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.

On the Coupling between the Collective Dynamics of Proteins and Their Hydration Water

Picosecond time scale dynamics of hydrated proteins has been connected with the onset of biological activity as it coincides with solvent-solute hydrogen bond rearrangements and amino acid rotational relaxation time scales. The presence and fluctuations of protein hydration water (PHW) largely influence protein motions that are believed to be slaved to those of the solvent, yet to date, how protein and hydration water dynamics are coupled remains unclear. Here, we provide a significant advance in characterizing this coupling; we present the first full study of both the longitudinal and transverse coherent collective motions in a protein-solvent system. The data show unexpectedly the presence in the water dynamics of collective modes belonging to the protein. The properties of these modes, in particular, their propagation velocities and amplitudes, indicate a strengthening of the interactions and a higher rigidity of the network of solvent molecules close to the protein surface. Accordingly, the present study presents the most compelling and clear evidence of a very strong dynamical coupling between a protein and its hydration water, previously suggested by studies using various experimental techniques.

Communication maps of vibrational energy transport through Photoactive Yellow Protein

We calculate communication maps for Photoactive Yellow Protein (PYP) from the purple phototropic eubacterium Halorhodospira halophile and use them to elucidate energy transfer pathways from the chromophore through the rest of the protein in the ground and excited state. The calculations reveal that in PYP excess energy from the chromophore flows mainly to regions of the surrounding residues that hydrogen bond to the chromophore. In addition, quantum mechanics/molecular mechanics and molecular dynamics (MD) simulations of the dielectric response of the protein and solvent environment due to charge rearrangement on the chromophore following photoexcitation are also presented, with both approaches yielding similar time constants for the response. Results of MD simulations indicate that the residues hydrogen bonding to the chromophore make the largest contribution to the response.

Anomalous water diffusion in salt solutions

The dynamics of water exhibits anomalous behavior in the presence of different electrolytes. Recent experiments [Kim JS, Wu Z, Morrow AR, Yethiraj A, Yethiraj A (2012) J Phys Chem B 116(39):12007-12013] have found that the self-diffusion of water (Dw) can either be enhanced or suppressed around CsI and NaCl, respectively, relative to that of neat water. Here we show that unlike classical empirical potentials, ab initio molecular dynamics simulations successfully reproduce the qualitative trends observed experimentally. These types of phenomena have often been rationalized in terms of the "structure-making" or "structure-breaking" effects of different ions on the solvent, although the microscopic origins of these features have remained elusive. Rather than disrupting the network in a significant manner, the electrolytes studied here cause rather subtle changes in both structural and dynamical properties of water. In particular, we show that water in the ab initio molecular dynamics simulations is characterized by dynamic heterogeneity, which turns out to be critical in reproducing the experimental trends.

Protein electron transfer: Dynamics and statistics

Electron transfer between redox proteins participating in energy chains of biology is required to proceed with high energetic efficiency, minimizing losses of redox energy to heat. Within the standard models of electron transfer, this requirement, combined with the need for unidirectional (preferably activationless) transitions, is translated into the need to minimize the reorganization energy of electron transfer. This design program is, however, unrealistic for proteins whose active sites are typically positioned close to the polar and flexible protein-water interface to allow inter-protein electron tunneling. The high flexibility of the interfacial region makes both the hydration water and the surface protein layer act as highly polar solvents. The reorganization energy, as measured by fluctuations, is not minimized, but rather maximized in this region. Natural systems in fact utilize the broad breadth of interfacial electrostatic fluctuations, but in the ways not anticipated by the standard models based on equilibrium thermodynamics. The combination of the broad spectrum of static fluctuations with their dispersive dynamics offers the mechanism of dynamical freezing (ergodicity breaking) of subsets of nuclear modes on the time of reaction/residence of the electron at a redox cofactor. The separation of time-scales of nuclear modes coupled to electron transfer allows dynamical freezing. In particular, the separation between the relaxation time of electro-elastic fluctuations of the interface and the time of conformational transitions of the protein caused by changing redox state results in dynamical freezing of the latter for sufficiently fast electron transfer. The observable consequence of this dynamical freezing is significantly different reorganization energies describing the curvature at the bottom of electron-transfer free energy surfaces (large) and the distance between their minima (Stokes shift, small). The ratio of the two reorganization energies establishes the parameter by which the energetic efficiency of protein electron transfer is increased relative to the standard expectations, thus minimizing losses of energy to heat. Energetically efficient electron transfer occurs in a chain of conformationally quenched cofactors and is characterized by flattened free energy surfaces, reminiscent of the flat and rugged landscape at the stability basin of a folded protein.

Direct probing of solvent accessibility and mobility at the binding interface of polymerase (Dpo4)-DNA complex

Water plays essential structural and dynamical roles in protein-DNA recognition through contributing to enthalpic or entropic stabilization of binding complex and by mediating intermolecular interactions and fluctuations for biological function. These interfacial water molecules are confined by the binding partners in nanospace, but in many cases they are highly mobile and exchange with outside bulk solution. Here, we report our studies of the interfacial water dynamics in the binary and ternary complexes of a polymerase (Dpo4) with DNA and an incoming nucleotide using a site-specific tryptophan probe with femtosecond resolution. By systematic comparison of the interfacial water motions and local side chain fluctuations in the apo, binary, and ternary states of Dpo4, we observed that the DNA binding interface and active site are dynamically solvent accessible and the interfacial water dynamics are similar to the surface hydration water fluctuations on picosecond time scales. Our molecular dynamics simulations also show the binding interface full of water molecules and nonspecific weak interactions. Such a fluid binding interface facilitates the polymerase sliding on DNA for fast translocation whereas the spacious and mobile hydrated active site contributes to the low fidelity of the lesion-bypass Y-family DNA polymerase.

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.

Surface Polarity and Nanoscale Solvation

Linear solvation theories are well established to describe electrostatic hydration of small solutes when the hydration free energy is dominated by the electrostatic free energy of the solute multipole. In contrast, hydration of nanometer solutes is driven by surface hydration. We address the question of whether the linear-response thermodynamics established for small multipolar solutes applies to surface hydration. To this end, molecular dynamics simulations are carried out on a model C180 solute that carries no global multipole, but the surface of which is decorated with radially pointing dipoles. Linear response is dramatically violated in this case. Further, two crossovers in the solvation thermodynamics are discovered as the surface polarity is increased. Both transformations produce strongly nonlinear solvation response. The second, more collective, crossover leads to a dramatic slowing down of the interfacial dynamics, reaching the time-scales of nanoseconds. Our picture offers the possibility of flipping water domains at interfaces of nanoparticles and biomolecules.

Non-Gaussian statistics and nanosecond dynamics of electrostatic fluctuations affecting optical transitions in proteins

We show that electrostatic fluctuations of the protein-water interface are globally non-Gaussian. The electrostatic component of the optical transition energy (energy gap) in a hydrated green fluorescent protein is studied here by classical molecular dynamics simulations. The distribution of the energy gap displays a high excess in the breadth of electrostatic fluctuations over the prediction of the Gaussian statistics. The energy gap dynamics include a nanosecond component. When simulations are repeated with frozen protein motions, the statistics shifts to the expectations of linear response and the slow dynamics disappear. We therefore suggest that both the non-Gaussian statistics and the nanosecond dynamics originate largely from global, low-frequency motions of the protein coupled to the interfacial water. The non-Gaussian statistics can be experimentally verified from the temperature dependence of the first two spectral moments measured at constant-volume conditions. Simulations at different temperatures are consistent with other indicators of the non-Gaussian statistics. In particular, the high-temperature part of the energy gap variance (second spectral moment) scales linearly with temperature and extrapolates to zero at a temperature characteristic of the protein glass transition. This result, violating the classical limit of the fluctuation-dissipation theorem, leads to a non-Boltzmann statistics of the energy gap and corresponding non-Arrhenius kinetics of radiationless electronic transitions, empirically described by the Vogel-Fulcher-Tammann law.

Local order and mobility of water molecules around ambivalent helices

Water on a protein surface plays a key role in determining the structure and dynamics of proteins. Compared to the properties of bulk water, many aspects of the structure and dynamics of the water surrounding the proteins are less understood. It is interesting therefore to explore how the properties of the water within the solvation shell around the peptide molecule depend on its specific secondary structure. In this work we investigate the orientational order and residence times of the water molecules to characterize the structure, energetics, and dynamics of the hydration shell water around ambivalent peptides. Ambivalent sequences are identical sequences which display multiple secondary structures in different proteins. Molecular dynamics simulations of representative proteins containing variable helix, variable nonhelix, and conserved helix are also used to explore the local structure and mobility of water molecules in their vicinity. The results, for the first time, depict a different water distribution pattern around the conserved and variable helices. The water molecules surrounding the helical segments in variable helices are found to possess a less locally ordered structure compared to those around their corresponding nonhelical counterparts and conserved helices. The long conserved helices exhibit extremely high local residence times compared to the helical conformations of the variable helices, whereas the residence times of the nonhelical conformations of the variable helices are comparable to those of the short conserved helices. This differential pattern of the structure and dynamics of water molecules in the vicinity of conserved/variable helices may lend valuable insights for understanding the role of solvent effects in determining sequence ambivalency and help in improving the accuracy of water models used in the simulations of proteins.

An AIMD study of CPD repair mechanism in water: role of solvent in ring splitting

In this paper, we continue to explore the repair mechanisms of the cyclobutane pyrimidine dimer. We find that a full description of both C5-C5' and C6-C6' bond splitting requires a multidimensional treatment involving a solvent coordinate in addition to changes in internal dimer coordinates. Nonequilibrium effects are likely to be important as well, although the initial conditions following forward electron transfer to the dimer, beyond the scope of this study, will ultimately determine the importance of these effects. Throughout the splitting of C5-C5' and C6-C6' bonds, a significant amount of excess charge is delocalized onto the solvent. We have verified that this is not an artifact of the electronic density functional theory (DFT) method used for this anionic system with Schrödinger equation-based quantum chemical cluster calculations. The amount and variability of charge delocalization changes with the course of the reaction. The splitting of the C6-C6' bond is accompanied by both an increase in electron density on the C6 and C6' carbon atoms and an increase in the water density near those atoms. These features are observed both in our equilibrium umbrella sampling simulations and nonequilibrium trajectories.

Probing deuterium isotope effect on structure and solvation dynamics of human serum albumin

The deuterium isotopic effect on the structure and solvation dynamics of the protein, human serum albumin (HSA), has been studied by using circular dichroism (CD), femtosecond up-conversion, FRET, and single-molecule spectroscopy. The CD spectra suggest that D(2)O affects the structure of HSA, leading to a 20% decrease in the helical structure. The FRET study indicates that the distance of C153 from the lone tryptophan residue of HSA is quite similar (≈21 Å) in H(2)O and D(2)O, and hence, the location of the probe in the protein remains the same in the two solvents. The single-molecule study suggests that coumarin 153 (C153) binds almost exclusively (>96%) to one site of HSA. Solvation dynamics of C153 in HSA is found to be markedly retarded in D(2)O compared with H(2)O. In H(2)O, the solvation of C153 bound to HSA is found to be biexponential with one component of 7 ps (30%) and a long component of 350 ps (70%). In D(2)O, we detected a short component of 4 ps (41%) and a long component of 950 ps (59%). Thus, the ultraslow component of the solvation dynamics of C153 bound to HSA in D(2)O (950 ps) is 2.5-fold slower than that in H(2)O (350 ps). The marked deuterium isotope effect has been ascribed to water molecules confined in the protein environment and to a lesser extent to the structural modification of protein by D(2)O.

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.

Relaxation dynamics of nucleosomal DNA

Recent experimental and theoretical evidence demonstrates that proteins and water in the hydration layer can follow complex stretched exponential or power law relaxation dynamics. Here, we report on a 50 ns all atom molecular dynamics (MD) simulation of the yeast nucleosome, where the interactions between DNA, histones, surrounding water and ions are explicitly included. DNA interacts with the histone core in 14 locations, approximately every 10.4 base pairs. We demonstrate that all sites of interaction exhibit anomalously slow power law relaxation, extending up to 10 ns, while fast exponential relaxation dynamics of hundreds of picoseconds applies to DNA regions outside these locations. The appearance of 1/f(alpha) noise or pink noise in DNA dynamics is ubiquitous. For histone-bound nucleotide dynamics alpha --> 1 and is a signature of complexity of the protein-DNA interactions. For control purposes two additional DNA simulations free of protein are conducted. Both utilize the same sequence of DNA, as found the in the nucleosome. In one simulation the initial conformation of the double helix is a straight B-form. In the other, the initial conformation is super helical. Neither of these simulations exhibits the variation of alpha as a function of position, the measure of power law for dynamical behavior, which we observe in the nucleosome simulation. The unique correspondence (high alpha to DNA-histone interaction sites, low alpha to free DNA sites), suggests that alpha may be an important and new quantification of protein-DNA interactions for future experiments.