We show that protein flexibility can be characterized using graph theory, from a single protein conformation. Covalent and hydrogen bonds are modeled by distance and angular constraints, and a map is constructed of the regions in this network that are flexible or rigid, based on whether their dihedral bonds remain rotatable or are locked by other interactions in the network. This analysis takes only a second on a typical PC, and interatomic potentials; the most time-consuming aspect of molecular dynamics calculations, are not required. Our preliminary work has shown that this approach identifies the experimentally observed, biologically important flexible regions in HIV protease and lysine-arginine-ornithine binding protein. Here we analyze three evolutionarily distant cytochromes c, and find strong similarity between their flexible regions, despite having only 39% sequence identity. Furthermore, we show how the structural flexibility increases as the weaker hydrogen bonds are removed, as would happen under thermal denaturation of the protein. This approach identifies the critical hydrogen bonds that cross-link the tertiary structure.