Evidence of protein collective motions on the picosecond timescale

We investigate the presence of structural collective motions on a picosecond timescale for the heme protein, cytochrome c, as a function of oxidation and hydration, using terahertz (THz) time domain spectroscopy and molecular dynamics simulations. The THz response dramatically increases with oxidation, with the largest increase for lowest hydrations, and highest frequencies. For both oxidation states the THz response rapidly increases with hydration saturating above ∼25% (g H₂O/g protein). Quasiharmonic vibrational modes and dipole-dipole correlation functions were calculated from molecular dynamics trajectories. The collective mode density of states alone reproduces the measured hydration dependence, providing strong evidence of the existence of these motions. The large oxidation dependence is reproduced only by the dipole-dipole correlation function, indicating the contrast arises from diffusive motions consistent with structural changes occurring in the vicinity of buried internal water molecules. This source for the observed oxidation dependence is consistent with the lack of an oxidation dependence in nuclear resonant vibrational spectroscopy measurements.