Towards Flexibility and Interpretability of Gaussian Process State-Space Model