The model was encoded and simulated using COPASI 4.24 and the figure was plotted using R 3.5.1.