In this study, it is demonstrated that the motion of hinges in single protein molecules can be modeled as general fractional Langevin dynamics of harmonically bound particles driven by non-local Gaussian noise with a power-law correlation. This conclusion was justified by comparing the theoretical predictions of the proposed model to the existing molecular dynamics simulations performed on the M6I mutant of phage T4 lysozyme and E. coli ribonuclease H1.