  The concentration gradient of the Bicoid morphogen, which is established
during the early stages of a Drosophila melanogaster embryonic development,
determines the differential spatial patterns of gene expression and subsequent
cell fate determination. This is mainly achieved by diffusion elicited by the
different concentrations of the Bicoid protein in the embryo. Such chemical
dynamic progress can be simulated by stochastic models, particularly the
Gillespie alogrithm. However, as with various modelling approaches in biology,
each technique involves drawing assumptions and reducing the model complexity
sometimes limiting the model capability. This is mainly due to the complexity
of the software modelling approaches to construct these models. Agent-based
modelling is a technique which is becoming increasingly popular for modelling
the behaviour of individual molecules or cells in computational biology.
  This paper attempts to compare these two popular modelling techniques of
stochastic and agent-based modelling to show how the model can be studied in
detail using the different approaches. This paper presents how to use these
techniques with the advantages and disadvantages of using either of these.
Through various comparisons, such as computation complexity and results
obtained, we show that although the same model is implemented, both approaches
can give varying results. The results of the paper show that the stochastic
model is able to give smoother results compared to the agent-based model which
may need further analysis at a later stage. We discuss the reasons for these
results and how these could be rectified in systems biology research.
