Resolving anisotropic distributions of correlated vibrational motion in protein hydration water

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.

Frequency-Selective Anharmonic Mode Analysis of Thermally Excited Vibrations in Proteins

Low-frequency molecular vibrations at far-infrared frequencies are thermally excited at room temperature. As a consequence, thermal fluctuations are not limited to the immediate vicinity of local minima on the potential energy surface, and anharmonic properties cannot be ignored. The latter is particularly relevant in molecules with multiple conformations, such as proteins and other biomolecules. However, existing theoretical and computational frameworks for the analysis of molecular vibrations have so far been limited by harmonic or quasi-harmonic approximations, which are ill-suited to describe anharmonic low-frequency vibrations. Here, we introduce a fully anharmonic analysis of molecular vibrations based on a time correlation formalism that eliminates the need for harmonic or quasi-harmonic approximations. We use molecular dynamics simulations of a small protein to demonstrate that this new approach, in contrast to harmonic and quasi-harmonic normal modes, correctly identifies the collective degrees of freedom associated with molecular vibrations at any given frequency. This allows us to unambiguously characterize the anharmonic character of low-frequency vibrations in the far-infrared spectrum.

Molecular dynamics study of collective water vibrations in a DNA hydration shell

The structure of DNA double helix is stabilized by water molecules and metal counterions that form the ion-hydration shell around the macromolecule. Understanding the role of the ion-hydration shell in the physical mechanisms of the biological functioning of DNA requires detailed studies of its structure and dynamics at the atomistic level. In the present work, the study of collective vibrations of water molecules around the DNA double helix was performed within the framework of classical all-atom molecular dynamics methods. Calculating the vibrational density of states, the vibrations of water molecules in the low-frequency spectra ranged from [Formula: see text]30 to [Formula: see text]300 [Formula: see text] were analyzed for the case of different regions of the DNA double helix (minor groove, major groove, and phosphate groups). The analysis revealed significant differences in the collective vibrations behavior of water molecules in the DNA hydration shell, compared to the vibrations of bulk water. All low-frequency modes of the DNA ion-hydration shell are shifted by about 15-20 [Formula: see text] towards higher frequencies, which is more significant for water molecules in the minor groove region of the double helix. The interactions of water molecules with the atoms of the macromolecule induce intensity decrease of the mode of hydrogen-bond symmetrical stretching near 150 [Formula: see text], leading to the disappearance of this mode in the DNA spectra. The obtained results can provide an interpretation of the experimental data for DNA low-frequency spectra and may be important for the understanding of the processes of indirect protein-nucleic recognition.

Knowledge-based prediction of DNA hydration using hydrated dinucleotides as building blocks

Water plays an important role in stabilizing the structure of DNA and mediating its interactions. Here, the hydration of DNA was analyzed in terms of dinucleotide fragments from an ensemble of 2727 nonredundant DNA chains containing 41 853 dinucleotides and 316 265 associated first-shell water molecules. The dinucleotides were classified into categories based on their 16 sequences and the previously determined structural classes known as nucleotide conformers (NtCs). The construction of hydrated dinucleotide building blocks allowed dinucleotide hydration to be calculated as the probability of water density distributions. Peaks in the water densities, known as hydration sites (HSs), uncovered the interplay between base and sugar-phosphate hydration in the context of sequence and structure. To demonstrate the predictive power of hydrated DNA building blocks, they were then used to predict hydration in an independent set of crystal and NMR structures. In ten tested crystal structures, the positions of predicted HSs and experimental waters were in good agreement (more than 40% were within 0.5 Å) and correctly reproduced the known features of DNA hydration, for example the `spine of hydration' in B-DNA. Therefore, it is proposed that hydrated building blocks can be used to predict DNA hydration in structures solved by NMR and cryo-EM, thus providing a guide to the interpretation of experimental data and computer models. The data for the hydrated building blocks and the predictions are available for browsing and visualization at the website https://watlas.datmos.org/watna/.

Protein flexibility reduces solvent-mediated friction barriers of ligand binding to a hydrophobic surface patch

Solvent fluctuations have been explored in detail for idealized and rigid hydrophobic model systems, but so far it has remained unclear how internal protein motions and their coupling to the surrounding solvent affect the dynamics of ligand binding to biomolecular surfaces. Here, molecular dynamics simulations were used to elucidate the solvent-mediated binding of a model ligand to the hydrophobic surface patch of ubiquitin. The ligand's friction profiles reveal pronounced long-time correlations and enhanced friction in the vicinity of the protein, similar to idealized hydrophobic surfaces. Interestingly, these effects are shaped by internal protein motions. Protein flexibility modulates water density fluctuations near the hydrophobic surface patch and smooths out the friction profile of ligand binding.

Dynamics of hydration water in gelatin and hyaluronic acid hydrogels

We employed broadband dielectric spectroscopy (BDS), for the investigation of the water dynamics in partially hydrated hyaluronic acid (HA), and gelatin (Gel), enzymatically crosslinked hydrogels, in the water fraction ranges [Formula: see text]. Our results indicate that at low hydrations ([Formula: see text]), where the dielectric response of the hydrogels is identical during cooling and heating, water plasticizes strongly the polymeric matrix and is organized in clusters giving rise to [Formula: see text]-process, secondary water relaxation and to an additional slower relaxation process. This later process has been found to be related with the dc charge conductivity and can be described in terms of the conduction current relaxation mechanism. At slightly higher hydrations, however, always below the hydration level where ice is formed during cooling, we have recorded in HA hydrogel a strong water dielectric relaxation process, [Formula: see text], which has Arrhenius-like temperature dependence and large time scale resembling relaxation processes recorded in bulk low density amorphous solid water structures. This relaxation process shows a strong-to-fragile transition at [Formula: see text]C and our data suggest that the VTF-like process recorded at [Formula: see text]C is controlled by the same molecular process like long range charge transport. In addition, our data imply that the crossover temperature is related with the onset of structural rearrangements (increase in configurational entropy) of the macromolecules. In partially crystallized hydrogels ([Formula: see text]) HA exhibits at low temperatures the ice dielectric process consistent with the bulk hexagonal ice, whereas Gel hydrogel exhibits as main low temperature process a slow relaxation process that refers to open tetrahedral structures of water similar to low density amorphous ice structures and to bulk cubic ice. Regarding the water secondary relaxation processes, we have shown that the [Formula: see text]-process and the [Formula: see text] process are activated in water hydrogen bond networks with different structures.

Atomistic characterization of collective protein–water–membrane dynamics

Water mediates correlated vibrations of atoms of protein and membrane bilayer surfaces.

Hydration-mediated stiffening of collective membrane dynamics by cholesterol

Hydration water governs the cholesterol-induced changes in collective headgroup dynamics in lipid bilayers.

Aqueous solvation from the water perspective

The response of water re-solvating a charge-transfer dye (deprotonated Coumarin 343) after photoexcitation has been measured by means of transient THz spectroscopy. Two steps of increasing THz absorption are observed, a first ∼10 ps step on the time scale of Debye relaxation of bulk water and a much slower step on a 3.9 ns time scale, the latter of which reflecting heating of the bulk solution upon electronic relaxation of the dye molecules from the S₁ back into the S₀ state. As an additional reference experiment, the hydroxyl vibration of water has been excited directly by a short IR pulse, establishing that the THz signal measures an elevated temperature within ∼1 ps. This result shows that the first step upon dye excitation (10 ps) is not limited by the response time of the THz signal; it rather reflects the reorientation of water molecules in the solvation layer. The apparent discrepancy between the relatively slow reorientation time and the general notion that water is among the fastest solvents with a solvation time in the sub-picosecond regime is discussed. Furthermore, non-equilibrium molecular dynamics simulations have been performed, revealing a close-to-quantitative agreement with experiment, which allows one to disentangle the contribution of heating to the overall THz response from that of water orientation.

Pressure effects on collective density fluctuations in water and protein solutions

Neutron Brillouin scattering and molecular dynamics simulations have been used to investigate protein hydration water density fluctuations as a function of pressure. Our results show significant differences between the pressure and density dependence of collective dynamics in bulk water and in concentrated protein solutions. Pressure-induced changes in the tetrahedral order of the water HB network have direct consequences for the high-frequency sound velocity and damping coefficients, which we find to be a sensitive probe for changes in the HB network structure as well as the wetting of biomolecular surfaces.

On the ability of molecular dynamics simulation and continuum electrostatics to treat interfacial water molecules in protein-protein complexes

Interfacial waters are increasingly appreciated as playing a key role in protein-protein interactions. We report on a study of the prediction of interfacial water positions by both Molecular Dynamics and explicit solvent-continuum electrostatics based on the Dipolar Poisson-Boltzmann Langevin (DPBL) model, for three test cases: (i) the barnase/barstar complex (ii) the complex between the DNase domain of colicin E2 and its cognate Im2 immunity protein and (iii) the highly unusual anti-freeze protein Maxi which contains a large number of waters in its interior. We characterize the waters at the interface and in the core of the Maxi protein by the statistics of correctly predicted positions with respect to crystallographic water positions in the PDB files as well as the dynamic measures of diffusion constants and position lifetimes. Our approach provides a methodology for the evaluation of predicted interfacial water positions through an investigation of water-mediated inter-chain contacts. While our results show satisfactory behaviour for molecular dynamics simulation, they also highlight the need for improvement of continuum methods.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Protein Structural Memory Influences Ligand Binding Mode(s) and Unbinding Rates

The binding of small molecules (e.g., natural ligands, metabolites, and drugs) to proteins governs most biochemical pathways and physiological processes. Here, we use molecular dynamics to investigate the unbinding of dimethyl sulfoxide (DMSO) from two distinct states of a small rotamase enzyme, the FK506-binding protein (FKBP). These states correspond to the FKBP protein relaxed with and without DMSO in the active site. Since the time scale of ligand unbinding (2-20 ns) is faster than protein relaxation (100 ns), a novel methodology is introduced to relax the protein without having to introduce an artificial constraint. The simulation results show that the unbinding time is an order of magnitude longer for dissociation from the DMSO-bound state (holo-relaxed). That is, the actual rate of unbinding depends on the state of the protein, with the protein having a long-lived memory. The rate thus depends on the concentration of the ligand as the apo and holo states reflect low and high concentrations of DMSO, respectively. Moreover, there are multiple binding modes in the apo-relaxed state, while a single binding mode dominates the holo-relaxed state in which DMSO acts as hydrogen bond acceptor from the backbone NH of Ile56, as in the crystal structure of the DMSO/FKBP complex. The solvent relaxes very fast (∼1 ns) close to the NH of Ile56 and with the same time scale of the protein far away from the active site. These results have implications for high-throughput docking, which makes use of a rigid structure of the protein target.