Electron-vibration coupling in zinc substituted myoglobin has been calculated using a quantum mechanical/molecular mechanical method. The methodology has been tested by a direct comparison of the calculated optical observables, the steady-state optical spectra and three-pulse-photon-echo-peak-shift (3PEPS) function, to those experimentally measured showing a qualitative agreement. A range of experiments and calculations were performed to explain the discrepancies, which lead to the conclusion that the discrepancy originates from adiabatic coupling of the two nearly degenerate electronic transitions.