Surrogate Modeling with Gaussian Processes for Molecular Models of Viscoelasticity