Before revision:   To validate our approach we employ the Gillespiemethod to generate benchmark data. The numerical simulation shows that the data produced by both methods agree quite well .
After revision:  To validate our approach we employ Gillespie's stochastic simulation algorithm (SSA) to generate (forward) benchmark data, according to predefined gene network models. Numerical simulations show that the BSDE method is able to correctly infer the protein level distributions that preceded a known final condition, obtained originally from the forward SSA. This makes the BSDE method a powerful systems biology tool for time reversed simulations, allowing, for example, the assessment of the biological conditions (e.g. protein concentrations) that preceded an experimentally measured event of interest (e.g. mitosis, apoptosis, etc.) .