The distributions of the quantum vibrational energy levels of the protein molecular chain are found by the discretely nonlinear Schrodinger equation appropriate to protein obtained from the Davydov theory. The results calculated by this method are basically consistent with the experimental values. Furthermore, the energy spectra at high excited states have also been obtained for the molecular chain which is helpful in researching the properties of infrared absorption and Raman scattering of the protein molecules.