Separable cooperative and localized translational motions of water confined by a chemically heterogeneous environment

Photo-Switching of Protein Dynamical Collectivity

We examine changes in the picosecond structural dynamics with irreversible photobleaching of red fluorescent proteins (RFP) mCherry, mOrange2 and TagRFP-T. Measurements of the protein dynamical transition using terahertz time-domain spectroscopy show in all cases an increase in the turn-on temperature in the bleached state. The result is surprising given that there is little change in the protein surface, and thus, the solvent dynamics held responsible for the transition should not change. A spectral analysis of the measurements guided by quasiharmonic calculations of the protein absorbance reveals that indeed the solvent dynamical turn-on temperature is independent of the thermal stability/photostate however the protein dynamical turn-on temperature shifts to higher temperatures. This is the first demonstration of switching the protein dynamical turn-on temperature with protein functional state. The observed shift in protein dynamical turn-on temperature relative to the solvent indicates an increase in the required mobile waters necessary for the protein picosecond motions, that is, these motions are more collective. Melting-point measurements reveal that the photobleached state is more thermally stable, and structural analysis of related RFP’s shows that there is an increase in internal water channels as well as a more uniform atomic root mean squared displacement. These observations are consistent with previous suggestions that water channels form with extended light excitation providing O2 access to the chromophore and subsequent fluorescence loss. We report that these same channels increase internal coupling enhancing thermal stability and collectivity of the picosecond protein motions. The terahertz spectroscopic characterization of the protein and solvent dynamical onsets can be applied generally to measure changes in collectivity of protein motions.

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.

Concentration dependence of hydration water in a model peptide

The molecular dynamics of aqueous solutions of a model amphiphilic peptide is studied as a function of concentration by broad-band light scattering experiments. Similarly to protein aqueous solutions, a considerable retardation, of about a factor 6-8, of hydration water dynamics with respect to bulk water is found, showing a slight dependence on solute concentration. Conversely, the average number of water molecules perturbed by the presence of peptide, i.e. the hydration number, appears to be strongly modified by adding solute. Its behaviour, decreasing upon increasing concentration, can be interpreted considering the random close-to-contact condition experienced by solute particles. Overall, the present findings support the view of a "long range" effect of peptides on the surrounding water, extending beyond the first two hydration shells.

Properties of water in the region between a tubulin dimer and a single motor head of kinesin

A kinesin is a molecular motor that can perform movement on a microtubule track in a stepping-like manner. This motion is connected with processes of association and dissociation of kinesin and tubulin. Water is an important participant in these kinds of molecular interactions. This is why we have decided to investigate the dynamical and structural properties of water in the region between the kinesin catalytic domain and the tubulin dimer. Using the molecular dynamics method, we found that these properties are different from the ones of bulk water. The changes in structure and dynamics are visible for water beyond the first solvation layers, even for the longest analyzed distance between proteins equal to 2.0 nm. However, these changes are not always enhanced compared to the situation when only one protein surface is present. One factor that distinguishes the investigated situation from the one with a single protein is the presence of an additional electric field originating from the second protein. The tendency of vectors of dipole moments of water molecules between the proteins to follow the vectors of electric field generated by the proteins causes a distortion of the water-water hydrogen bond network. It has been shown that this distortion affects the properties of water in this region: it induces structural changes in solvation water, and leads to increased water density and increased stiffness of the water structure.

Current status of the AMOEBA polarizable force field

Molecular force fields have been approaching a generational transition over the past several years, moving away from well-established and well-tuned, but intrinsically limited, fixed point charge models toward more intricate and expensive polarizable models that should allow more accurate description of molecular properties. The recently introduced AMOEBA force field is a leading publicly available example of this next generation of theoretical model, but to date, it has only received relatively limited validation, which we address here. We show that the AMOEBA force field is in fact a significant improvement over fixed charge models for small molecule structural and thermodynamic observables in particular, although further fine-tuning is necessary to describe solvation free energies of drug-like small molecules, dynamical properties away from ambient conditions, and possible improvements in aromatic interactions. State of the art electronic structure calculations reveal generally very good agreement with AMOEBA for demanding problems such as relative conformational energies of the alanine tetrapeptide and isomers of water sulfate complexes. AMOEBA is shown to be especially successful on protein-ligand binding and computational X-ray crystallography where polarization and accurate electrostatics are critical.

Effects of co-solvents on peptide hydration water structure and dynamics

We evaluate the molecular response of hydration water as a function of temperature and proximity to the surface of the peptide N-acetyl-leucine-methylamide (NALMA) when in the presence of the kosmotrope co-solvent glycerol or the chaotrope co-solvent dimethyl sulfoxide (DMSO), using molecular dynamics simulation with a polarizable force field. These detailed microscopic studies complement established thermodynamic analysis on the role of co-solvents in shifting the equilibrium for proteins away from or towards the native folded state. We find that the structure of the water at the peptide interfaces reflects an increase in hydration number in the glycerol solution and a decrease in hydration numbers in the DMSO solution. While the water dynamics around NALMA in the presence of both co-solvents is slower than that observed with the water solvent alone, in the DMSO mixture we no longer measure a separation in water motion time scales at low temperatures as is seen in the pure water solvent, but rather one single relaxation time. In the glycerol, however, we do observe a separation of time scales at low temperatures, supporting the hypothesis that hydration water near a hydrophobic solute evolves on a separate time scale than the extensive hydrogen-bonding network of more bulk-like water. Our simulation studies highlight the differences in the two co-solvent solutions due to the relative frequency of water contacts with the hydrophobic vs. hydrophilic peptide surface, and direct water interactions with the co-solvents.

Hydration water dynamics near biological interfaces

We performed classical molecular dynamics simulations using both fixed-charge and polarizable water and protein force fields to contrast the hydration dynamics near hydrophilic and amphiphilic peptides as a function of temperature. The high peptide concentrations we use serve as a model for the surface of folded proteins where hydration layers around each residue overlap significantly. Through simulation we determine that there are notable differences in the water dynamics analyzed from the outer and inner hydration layer regions of the amphiphilic peptide solution that explains the experimentally observed presence of two translational relaxations, while the hydrophilic peptide solution shows only a single non-Arrhenius translational process with no distinction between hydration layers. Given that water dynamics for the amphiphilic peptide system reproduces all known rotational and translational hydration dynamical anomalies exhibited by hydration water near protein surfaces, our analysis provides strong evidence that dynamical signatures near biological interfaces arises because of frustration in the hydration dynamics induced by chemical heterogeneity, as opposed to just topological roughness, of the protein surface.

Time scales of water dynamics at biological interfaces: peptides, proteins and cells

Water 2H and 17O spin relaxation is used to study water dynamics in the hydration layers of two small peptides, two globular proteins and in living cells of two microorganisms. The dynamical heterogeneity of hydration water is characterized by performing relaxation measurements over a wide temperature range, extending deeply into the supercooled regime, or by covering a wide frequency range. Protein hydration layers can be described by a power-law distribution of rotational correlation times with an exponent close to 2. This distribution comprises a small fraction of protein-specific hydration sites, where water rotation is strongly retarded, and a dominant fraction of generic hydration sites, where water rotation is as fast as in the hydration shells of small peptides. The generic dynamic perturbation factor is less than 2 at room temperature and exhibits a maximum near 260 K. The dynamic perturbation is induced by H-bond constraints that interfere with the cooperative mechanism that facilitates rotation in bulk water. Because these constraints are temperature-independent, hydration water does not follow the super-Arrhenius temperature dependence of bulk water. Water in living cells behaves as expected from studies of simpler model systems, the only difference being a larger fraction of secluded (strongly perturbed) hydration sites associated with the supramolecular organization in the cell. Intracellular water that is not in direct contact with biopolymers has essentially the same dynamics as bulk water. There is no significant difference in cell water dynamics between mesophilic and halophilic organisms, despite the high K+ and Na+ concentrations in the latter.

Single particle and collective hydration dynamics for hydrophobic and hydrophilic peptides

We have conducted extensive molecular dynamics simulations to study the single particle and collective dynamics of water in solutions of N-acetyl-glycine-methylamide, a model hydrophilic protein backbone, and N-acetyl-leucine-methylamide, a model (amphiphilic) hydrophobic peptide, as a function of peptide concentration. Various analytical models commonly used in the analysis of incoherent quasielastic neutron scattering (QENS), are tested against the translational and rotational intermediate scattering function, the mean square displacement of the water molecule center of mass, and fits to the second-order rotational correlation function of water evaluated directly from the simulation data. We find that while the agreement between the model-free analysis and analytical QENS models is quantitatively poor, the qualitative feature of dynamical heterogeneity due to caging is captured well by all approaches. The center of mass collective and single particle intermediate scattering functions of water calculated for these peptide solutions show that the crossover from collective to single particle-dominated motions occurs at a higher value of Q for high concentration solutions relative to low concentration because of the greater restriction in movement of water molecules due to confinement. Finally, we have shown that at the same level of confinement of the two peptides, the aqueous amphiphilic amino acid solution shows the strongest deviation between single particle and collective dynamics relative to the hydrophilic amino acid, indicating that chemical heterogeneity induces even greater spatial heterogeneity in the water dynamics.