The aim of this study is to construct invariant regions in which we can establish global existence of classical solutions for reaction-diffusion systems with a general full matrix of diffusion coefficients. Our techniques are based on invariant regions and Lyapunov functional methods. The nonlinear reaction term has been supposed to be of exponential growth.