Recent computational methodologies, such as individual-based modelling, pave the way to the search for explanatory insight into the collective behaviour of molecules. Many reviews offer an up-to-date perspective about such methodologies, but little is discussed about the practical information requirements involved. The biological information used as input should be easily and routinely determined in the laboratory, publicly available and, preferably, organized in programmatically accessible databases.

This review is the first to provide a systematic and comprehensive overview of available resources for the modelling of metabolic events at the molecular scale. The glycolysis pathway of Escherichia coli , which is one of the most studied pathways in Microbiology, serves as case study. This curation addressed structural information about E. coli (i.e. defining the simulation environment), the reactions forming the glycolysis pathway including the enzymes and the metabolites (i.e. the molecules to be represented), the kinetics of each reaction (i.e. behavioural logic of the molecules) and diffusion parameters for all enzymes and metabolites (i.e. molecule movement in the environment).

Furthermore, the interpretation of relevant biological features, such as molecular diffusion and enzyme kinetics, and the connection of experimental determination and simulation validation are detailed. Notably, the information from classical theories, such as enzymatic rates and diffusion coefficients, is translated to simulation parameters, such as collision efficiency and particle velocity.

This work presents a molecular-scale agent-based model for the simulation of enzymatic reactions at experimentally measured concentrations. The model incorporates stochasticity and spatial dependence, using diffusing and reacting particles with physical dimensions. We developed strategies to adjust and validate the enzymatic rates and diffusion coefficients to the information required by the computational agents, i.e., collision efficiency, interaction logic between agents, the time scale associated with interactions (e.g., kinetics), and agent velocity. Also, we tested the impact of molecular location (a source of biological noise) in the speed at which the reactions take place.

Simulations were conducted for experimental data on the 2-hydroxymuconate tautomerase (EC 5.3.2.6, UniProt ID Q01468) and the Steroid Delta-isomerase (EC 5.3.3.1, UniProt ID P07445). Obtained results demonstrate that our approach is in accordance to existing experimental data and long-term biophysical and biochemical assumptions.

Agent-based simulations are increasingly popular in exploring and understanding cellular systems, but the natural complexity of these systems and the desire to grasp different modelling levels demand cost-effective simulation strategies and tools.

In this context, the present paper introduces novel sequential and distributed approaches for the three-dimensional agent-based simulation of individual molecules in cellular events. These approaches are able to describe the dimensions and position of the molecules with high accuracy and thus, study the critical effect of spatial distribution on cellular events. Moreover, two of the approaches allow multi-thread high performance simulations, distributing the three-dimensional model in a platform independent and computationally efficient way.

Evaluation addressed the reproduction of molecular scenarios and different scalability aspects of agent creation and agent interaction. The three approaches simulate common biophysical and biochemical laws faithfully. The distributed approaches show improved performance when dealing with large agent populations while the sequential approach is better suited for small to medium size agent populations.

Evaluation addressed the reproduction of molecular scenarios and different scalability aspects of agent creation and agent interaction. The three approaches simulate common biophysical and biochemical laws faithfully. The distributed approaches show improved performance when dealing with large agent populations while the sequential approach is better suited for small to medium size agent populations.

Experimental incapacity to track microbe-microbe interactions in structures like biofilms, and the complexity inherent to the mathematical modelling of those interactions, raises the need for feasible, alternative modelling approaches.

This work proposes an agent-based representation of the diffusion of N-acyl homoserine lactones (AHL) in a multicellular environment formed by Pseudomonas aeruginosa and Candida albicans. Depending on the spatial location, C. albicans cells were variably exposed to AHLs, an observation that might help explain why phenotypic switching of individual cells in biofilms occurred at different time points. The simulation and algebraic results were similar for simpler scenarios, although some statistical differences could be observed (p < 0.05).

The model was also successfully applied to a more complex scenario representing a small multicellular environment containing C. albicans and P. aeruginosa cells encased in a 3-D matrix. Further development of this model may help create a predictive tool to depict biofilm heterogeneity at the single-cell level.