We develop methodology for modelling dynamical biochemical reaction networks, e.g. for subcellular signalling pathways of nerve cells. The methodology includes parameter estimation, uncertainty quantification and sensitivity analysis (UQSA, see
https://doi.org/10.48550/arXiv.2308.05527). The model that we are currently considering use stochastic simulation and require likelihood-free methods for uncertainty quantification.