Hydration dependence of the mass fractal dimension and anomalous diffusion of vibrational energy in proteins

Locating dynamic contributions to allostery via determining rates of vibrational energy transfer

Determining rates of energy transfer across non-covalent contacts for different states of a protein can provide information about dynamic and associated entropy changes during transitions between states. We investigate the relationship between rates of energy transfer across polar and nonpolar contacts and contact dynamics for the β₂-adrenergic receptor, a rhodopsin-like G-protein coupled receptor, in an antagonist-bound inactive state and agonist-bound active state. From structures sampled during molecular dynamics (MD) simulations, we find the active state to have, on average, a lower packing density, corresponding to generally more flexibility and greater entropy than the inactive state. Energy exchange networks (EENs) are computed for the inactive and active states from the results of the MD simulations. From the EENs, changes in the rates of energy transfer across polar and nonpolar contacts are found for contacts that remain largely intact during activation. Change in dynamics of the contact, and entropy associated with the dynamics, can be estimated from the change in rates of energy transfer across the contacts. Measurement of change in the rates of energy transfer before and after the transition between states thereby provides information about dynamic contributions to activation and allostery.

Dependence of Vibrational Energy Transfer on Distance in a Four-Helix Bundle Protein: Equidistant Increments with the Periodicity of α Helices

Vibrational energy exchanges between various degrees of freedom are critical to barrier-crossing processes in proteins. Heme proteins are highly suitable for studies of the vibrational energy exchanges in proteins. The migration of excess energy released by heme in a protein moiety can be observed using time-resolved anti-Stokes ultraviolet resonance Raman spectroscopy. The anti-Stokes resonance Raman intensity of a tryptophan residue is an excellent probe for the excess energy and the spatial resolution of a single amino acid residue can be achieved. Here, we studied dependence of vibrational energy transfer on the distance in cytochrome b₅₆₂, which is a heme-containing, four-helix bundle protein. The vibrational energy transfer from the heme group to a single tryptophan residue introduced by site-directed mutagenesis was examined for different heme-tryptophan distances by a quasi-constant length with the periodicity of α helices. Taken together with structural data obtained by molecular dynamics simulations, the energy transfer could be well described by the model of classical thermal diffusion, which suggests that continuum media provide a good approximation of the protein interior, of which the atomic packing density is very high.

Energy Transfer across Nonpolar and Polar Contacts in Proteins: Role of Contact Fluctuations

Molecular dynamics simulations of the villin headpiece subdomain HP36 have been carried out to examine relations between rates of vibrational energy transfer across non-covalently bonded contacts and equilibrium structural fluctuations, with focus on van der Waals contacts. Rates of energy transfer across van der Waals contacts vary inversely with the variance of the contact length, with the same constant of proportionality for all nonpolar contacts of HP36. A similar relation is observed for hydrogen bonds, but the proportionality depends on contact pairs, with hydrogen bonds stabilizing the α-helices all exhibiting the same constant of proportionality, one that is distinct from those computed for other polar contacts. Rates of energy transfer across van der Waals contacts are found to be up to 2 orders of magnitude smaller than rates of energy transfer across polar contacts.

Dynamics of water clusters confined in proteins: a molecular dynamics simulation study of interfacial waters in a dimeric hemoglobin

Water confined in proteins exhibits dynamics distinct from the dynamics of water in the bulk or near the surface of a biomolecule. We examine the water dynamics at the interface of the two globules of the homodimeric hemoglobin from Scapharca inaequivalvis (HbI) by molecular dynamics (MD) simulations, with focus on water-protein hydrogen bond lifetimes and rotational anisotropy of the interfacial waters. We find that relaxation of the waters at the interface of both deoxy- and oxy-HbI, which contain a cluster of 17 and 11 interfacial waters, respectively, is well described by stretched exponentials with exponents from 0.1 to 0.6 and relaxation times of tens to thousands of picoseconds. The interfacial water molecules of oxy-HbI exhibit slower rotational relaxation and hydrogen bond rearrangement than those of deoxy-HbI, consistent with an allosteric transition from unliganded to liganded conformers involving the expulsion of several water molecules from the interface. Though the interfacial waters are translationally and rotationally static on the picosecond time scale, they contribute to fast communication between the globules via vibrations. We find that the interfacial waters enhance vibrational energy transport across the interface by ≈10%.

Water and backbone dynamics in a hydrated protein

Rotational immobilization of proteins permits characterization of the internal peptide and water molecule dynamics by magnetic relaxation dispersion spectroscopy. Using different experimental approaches, we have extended measurements of the magnetic field dependence of the proton-spin-lattice-relaxation rate by one decade from 0.01 to 300 MHz for (1)H and showed that the underlying dynamics driving the protein (1)H spin-lattice relaxation is preserved over 4.5 decades in frequency. This extension is critical to understanding the role of (1)H(2)O in the total proton-spin-relaxation process. The fact that the protein-proton-relaxation-dispersion profile is a power law in frequency with constant coefficient and exponent over nearly 5 decades indicates that the characteristics of the native protein structural fluctuations that cause proton nuclear spin-lattice relaxation are remarkably constant over this wide frequency and length-scale interval. Comparison of protein-proton-spin-lattice-relaxation rate constants in protein gels equilibrated with (2)H(2)O rather than (1)H(2)O shows that water protons make an important contribution to the total spin-lattice relaxation in the middle of this frequency range for hydrated proteins because of water molecule dynamics in the time range of tens of ns. This water contribution is with the motion of relatively rare, long-lived, and perhaps buried water molecules constrained by the confinement. The presence of water molecule reorientational dynamics in the tens of ns range that are sufficient to affect the spin-lattice relaxation driven by (1)H dipole-dipole fluctuations should make the local dielectric properties in the protein frequency dependent in a regime relevant to catalytically important kinetic barriers to conformational rearrangements.

Effects of surface water on protein dynamics studied by a novel coarse-grained normal mode approach

Normal mode analysis (NMA) has received much attention as a direct approach to extract the collective motions of macromolecules. However, the stringent requirement of computational resources by classical all-atom NMA limits the size of the macromolecules to which the method is normally applied. We implemented a novel coarse-grained normal mode approach based on partitioning the all-atom Hessian matrix into relevant and nonrelevant parts. It is interesting to note that, using classical all-atom NMA results as a reference, we found that this method generates more accurate results than do other coarse-grained approaches, including elastic network model and block normal mode approaches. Moreover, this new method is effective in incorporating the energetic contributions from the nonrelevant atoms, including surface water molecules, into the coarse-grained protein motions. The importance of such improvements is demonstrated by the effect of surface water to shift vibrational modes to higher frequencies and by an increase in overlap of the coarse-grained eigenvector space (the motion directions) with that obtained from molecular dynamics simulations of solvated protein in a water box. These results not only confirm the quality of our method but also point out the importance of incorporating surface structural water in studying protein dynamics.