Probability diffusion in non-Markhovian, non-Gaussian molecular ensembles: A theoretical analysis and computer simulation