Molecular fluctuations of the native conformation of c-AMP dependent protein kinase (cAPK) have been investigated with three different approaches. The first approach is the full atomic normal mode analysis (NMA) with empirical force fields. The second and third approaches are based on a coarse-grained model with a single single-parameter- harmonic potential between close residues in the crystal structure of the molecule without any residue specificity. The second method calculates only the magnitude of fluctuations whereas the third method is developed to find the directionality of the fluctuations which are essential to understand the functional importance of biological molecules. The aim, in this study, is to determine whether using such coarse-grained models are appropriate for elucidating the global dynamic characteristics of large proteins which reduces the size of the system at least by a factor of ten. The mean-square fluctuations of Cα atoms and the residue cross-correlations are obtained by three approaches. These results are then compared to test the results of coarse grained models on the overall collective motions. AH three of the approaches show that highly flexible regions correspond to the activation and solvent exposed loops, whereas the conserved residues (especially in substrate binding regions) exhibit almost no flexibility, adding stability to the structure. The anti-correlated motions of the two lobes of the catalytic core provide flexibility to the molecule. High similarities among the results of these methods indicate that the slowest modes governing the most global motions are preserved in the coarse grained models for proteins. This finding may suggest that the general shapes of the structures are representative of their dynamic characteristics and the dominant motions of protein structures are robust at coarse-grained levels.