We have investigated the relaxational dynamics for a protein model at various temperatures. Theoretical analysis of this model in conjunction with numerical simulations suggests several relaxation regimes, including a single exponential, a power law and a logarithmic time dependence. Even though a stretched exponential form gives a good fit to the simulation results in the crossover regime between a single exponential and a power law decay, we have not been able to directly deduce this form from the theoretical analysis.