In this work, the Gaussian Network Model (GNM) and Anisotropic Network Model (ANM) approaches are applied to describe the dynamics of protein structure graphs built from calculated promolecular electron density (ED) distribution functions. A first set of analyses is carried out on results obtained from ED maxima calculated at various smoothing levels. A second set is achieved for ED networks whose edges are weighted by ED overlap integral values. Results are compared with those obtained through the classical GNM and ANM approaches applied to networks of C, atoms. It is shown how the network model and the consideration of crystal packing as well as of the side chains may lead to various improvements dependent upon the structure under study. The selected protein structures are Crambin and Pancreatic Trypsin Inhibitor because of their small size and numerous dynamical data obtained by other authors.