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.