By introducing fractional Gaussian noise into the generalized Langevin equation, the subdiffusion of a particle can be described as a stationary Gaussian process with analytical tractability. This model is capable of explaining the equilibrium fluctuation of the distance between an electron transfer donor and acceptor pair within a protein that spans a broad range of time scales, and is in excellent agreement with a single-molecule experiment.