The vibrational energy relaxation of dissociated carbon monoxide in the heme pocket of sperm whale myoglobin has been studied using equilibrium molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the delta- and epsilon-tautomers of the distal histidine, His64. Vibrational population relaxation times were estimated using the Landau-Teller model. For carbon monoxide (CO) in the myoglobin epsilon-tautomer, for a frequency of omega0 = 2131 cm-1 corresponding to the B1 state, T1epsilon(B1) = 640 +/- 185 ps, and for a frequency of omega0 = 2119 cm-1 corresponding to the B2 state, T1epsilon(B2) = 590 +/- 175 ps. Although the CO relaxation rates in both the epsilon- and delta-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the delta-tautomer. For CO in the myoglobin delta-tautomer, it was found that the relaxation times were identical within error for the two CO substate frequencies, T1delta(B1) = 335 +/- 115 ps and T1delta(B2) = 330 +/- 145 ps. These simulation results are in reasonable agreement with experimental results of Anfinrud and coworkers (unpublished results). Normal mode calculations were used to identify the dominant coupling between the protein and CO molecules. The calculations suggest that the residues of the myoglobin pocket, acting as a first solvation shell to the CO molecule, contribute the primary "doorway" modes in the vibrational relaxation of the oscillator.