The evolution of vibrational wave packets built from the normal modes of cytochrome c, myoglobin and green fluorescent protein is investigated. Vibrational energy flow in these proteins is found to exhibit anomalous subdiffusion, a consequence of trapping of energy by spatially localized normal modes contained in the wave packet. Anomalous subdiffusion is characterized by an exponent, nu, that is related to the spectral dimension, (d) over bar, and fractal dimension, D, of the protein. The dispersion relation describing variation of the protein's normal mode frequencies with wave number is also characterized by an exponent, a, that is related to (d) over bar and D. Values of the exponent, a, computed for the three proteins are consistent with the computed values for nu. The values of D obtained from nu, a, and (d) over bar for each protein are the same within computational error, and close to the mass fractal dimension computed for each protein, all values falling in the range D=2.3+/-0.2. We find also that relaxation of the center of energy of a wave packet computed in terms of normal modes follows stretched exponential kinetics with an exponent equal to 2nu. (C) 2003 American Institute of Physics.