AgentCell: a digital single-cell assay for bacterial chemotaxis

MotGen: a closed-loop bacterial motility control framework using generative adversarial networks

MOTIVATION

Many organisms' survival and behavior hinge on their responses to environmental signals. While research on bacteria-directed therapeutic agents has increased, systematic exploration of real-time modulation of bacterial motility remains limited. Current studies often focus on permanent motility changes through genetic alterations, restricting the ability to modulate bacterial motility dynamically on a large scale. To address this gap, we propose a novel real-time control framework for systematically modulating bacterial motility dynamics.

RESULTS

We introduce MotGen, a deep learning approach leveraging Generative Adversarial Networks to analyze swimming performance statistics of motile bacteria based on live cell imaging data. By tracking objects and optimizing cell trajectory mapping under environmentally altered conditions, we trained MotGen on a comprehensive statistical dataset derived from real image data. Our experimental results demonstrate MotGen's ability to capture motility dynamics from real bacterial populations with low mean absolute error in both simulated and real datasets. MotGen allows us to approach optimal swimming conditions for desired motility statistics in real-time. MotGen's potential extends to practical biomedical applications, including immune response prediction, by providing imputation of bacterial motility patterns based on external environmental conditions. Our short-term, in-situ interventions for controlling motility behavior offer a promising foundation for the development of bacteria-based biomedical applications.

AVAILABILITY AND IMPLEMENTATION

MotGen is presented as a combination of Matlab image analysis code and a machine learning workflow in Python. Codes are available at https://github.com/bgmseo/MotGen, for cell tracking and implementation of trained models to generate bacterial motility statistics.

Learning black- and gray-box chemotactic PDEs/closures from agent based Monte Carlo simulation data

We propose a machine learning framework for the data-driven discovery of macroscopic chemotactic Partial Differential Equations (PDEs)-and the closures that lead to them- from high-fidelity, individual-based stochastic simulations of Escherichia coli bacterial motility. The fine scale, chemomechanical, hybrid (continuum-Monte Carlo) simulation model embodies the underlying biophysics, and its parameters are informed from experimental observations of individual cells. Using a parsimonious set of collective observables, we learn effective, coarse-grained "Keller-Segel class" chemotactic PDEs using machine learning regressors: (a) (shallow) feedforward neural networks and (b) Gaussian Processes. The learned laws can be black-box (when no prior knowledge about the PDE law structure is assumed) or gray-box when parts of the equation (e.g. the pure diffusion part) is known and "hardwired" in the regression process. More importantly, we discuss data-driven corrections (both additive and functional), to analytically known, approximate closures.

Agent-Based Modeling of Microbial Communities

Microbial communities are complex living systems that populate the planet with diverse functions and are increasingly harnessed for practical human needs. To deepen the fundamental understanding of their organization and functioning as well as to facilitate their engineering for applications, mathematical modeling has played an increasingly important role. Agent-based models represent a class of powerful quantitative frameworks for investigating microbial communities because of their individualistic nature in describing cells, mechanistic characterization of molecular and cellular processes, and intrinsic ability to produce emergent system properties. This review presents a comprehensive overview of recent advances in agent-based modeling of microbial communities. It surveys the state-of-the-art algorithms employed to simulate intracellular biomolecular events, single-cell behaviors, intercellular interactions, and interactions between cells and their environments that collectively serve as the driving forces of community behaviors. It also highlights three lines of applications of agent-based modeling, namely, the elucidation of microbial range expansion and colony ecology, the design of synthetic gene circuits and microbial populations for desired behaviors, and the characterization of biofilm formation and dispersal. The review concludes with a discussion of existing challenges, including the computational cost of the modeling, and potential mitigation strategies.

Individual-Based Modeling of Spatial Dynamics of Chemotactic Microbial Populations

One important direction of synthetic biology is to establish desired spatial structures from microbial populations. Underlying this structural development process are different driving factors, among which bacterial motility and chemotaxis serve as a major force. Here, we present an individual-based, biophysical computational framework for mechanistic and multiscale simulation of the spatiotemporal dynamics of motile and chemotactic microbial populations. The framework integrates cellular movement with spatial population growth, mechanical and chemical cellular interactions, and intracellular molecular kinetics. It is validated by a statistical comparison of single-cell chemotaxis simulations with reported experiments. The framework successfully captures colony range expansion of growing isogenic populations and also reveals chemotaxis-modulated, spatial patterns of a two-species amensal community. Partial differential equation-based models subsequently validate these simulation findings. This study provides a versatile computational tool to uncover the fundamentals of microbial spatial ecology as well as to facilitate the design of synthetic consortia for desired spatial patterns.

Emerging computational paradigms to address the complex role of gut microbial metabolism in cardiovascular diseases

The human gut microbiota and its associated perturbations are implicated in a variety of cardiovascular diseases (CVDs). There is evidence that the structure and metabolic composition of the gut microbiome and some of its metabolites have mechanistic associations with several CVDs. Nevertheless, there is a need to unravel metabolic behavior and underlying mechanisms of microbiome-host interactions. This need is even more highlighted when considering that microbiome-secreted metabolites contributing to CVDs are the subject of intensive research to develop new prevention and therapeutic techniques. In addition to the application of high-throughput data used in microbiome-related studies, advanced computational tools enable us to integrate omics into different mathematical models, including constraint-based models, dynamic models, agent-based models, and machine learning tools, to build a holistic picture of metabolic pathological mechanisms. In this article, we aim to review and introduce state-of-the-art mathematical models and computational approaches addressing the link between the microbiome and CVDs.

Steady-state running rate sets the speed and accuracy of accumulation of swimming bacteria

We study the chemotaxis of a population of genetically identical swimming bacteria undergoing run and tumble dynamics driven by stochastic switching between clockwise and counterclockwise rotation of the flagellar rotary system, where the steady-state rate of the switching changes in different environments. Understanding chemotaxis quantitatively requires that one links the measured steady-state switching rates of the rotary system, as well as the directional changes of individual swimming bacteria in a gradient of chemoattractant/repellant, to the efficiency of a population of bacteria in moving up/down the gradient. Here we achieve this by using a probabilistic model, parametrized with our experimental data, and show that the response of a population to the gradient is complex. We find the changes to the steady-state switching rate in the absence of gradients affect the average speed of the swimming bacterial population response as well as the width of the distribution. Both must be taken into account when optimizing the overall response of the population in complex environments.

BioDynaMo: a modular platform for high-performance agent-based simulation

MOTIVATION

Agent-based modeling is an indispensable tool for studying complex biological systems. However, existing simulation platforms do not always take full advantage of modern hardware and often have a field-specific software design.

RESULTS

We present a novel simulation platform called BioDynaMo that alleviates both of these problems. BioDynaMo features a modular and high-performance simulation engine. We demonstrate that BioDynaMo can be used to simulate use cases in: neuroscience, oncology and epidemiology. For each use case, we validate our findings with experimental data or an analytical solution. Our performance results show that BioDynaMo performs up to three orders of magnitude faster than the state-of-the-art baselines. This improvement makes it feasible to simulate each use case with one billion agents on a single server, showcasing the potential BioDynaMo has for computational biology research.

AVAILABILITY AND IMPLEMENTATION

BioDynaMo is an open-source project under the Apache 2.0 license and is available at www.biodynamo.org. Instructions to reproduce the results are available in the supplementary information.

SUPPLEMENTARY INFORMATION

Available at https://doi.org/10.5281/zenodo.5121618.

Rules of Engagement: A Guide to Developing Agent-Based Models

Agent-based models (ABM), also called individual-based models, first appeared several decades ago with the promise of nearly real-time simulations of active, autonomous individuals such as animals or objects. The goal of ABMs is to represent a population of individuals (agents) interacting with one another and their environment. Because of their flexible framework, ABMs have been widely applied to study systems in engineering, economics, ecology, and biology. This chapter is intended to guide the users in the development of an agent-based model by discussing conceptual issues, implementation, and pitfalls of ABMs from first principles. As a case study, we consider an ABM of the multi-scale dynamics of cellular interactions in a microbial community. We develop a lattice-free agent-based model of individual cells whose actions of growth, movement, and division are influenced by both their individual processes (cell cycle) and their contact with other cells (adhesion and contact inhibition).

Chemotactic migration of bacteria in porous media

Chemotactic migration of bacteria-their ability to direct multicellular motion along chemical gradients-is central to processes in agriculture, the environment, and medicine. However, current understanding of migration is based on studies performed in bulk liquid, despite the fact that many bacteria inhabit tight porous media such as soils, sediments, and biological gels. Here, we directly visualize the chemotactic migration of Escherichia coli populations in well-defined 3D porous media in the absence of any other imposed external forcing (e.g., flow). We find that pore-scale confinement is a strong regulator of migration. Strikingly, cells use a different primary mechanism to direct their motion in confinement than in bulk liquid. Furthermore, confinement markedly alters the dynamics and morphology of the migrating population-features that can be described by a continuum model, but only when standard motility parameters are substantially altered from their bulk liquid values to reflect the influence of pore-scale confinement. Our work thus provides a framework to predict and control the migration of bacteria, and active matter in general, in complex environments.

Tracing the footsteps of autophagy in computational biology

Autophagy plays a crucial role in maintaining cellular homeostasis through the degradation of unwanted materials like damaged mitochondria and misfolded proteins. However, the contribution of autophagy toward a healthy cell environment is not only limited to the cleaning process. It also assists in protein synthesis when the system lacks the amino acids’ inflow from the extracellular environment due to diet consumptions. Reduction in the autophagy process is associated with diseases like cancer, diabetes, non-alcoholic steatohepatitis, etc., while uncontrolled autophagy may facilitate cell death. We need a better understanding of the autophagy processes and their regulatory mechanisms at various levels (molecules, cells, tissues). This demands a thorough understanding of the system with the help of mathematical and computational tools. The present review illuminates how systems biology approaches are being used for the study of the autophagy process. A comprehensive insight is provided on the application of computational methods involving mathematical modeling and network analysis in the autophagy process. Various mathematical models based on the system of differential equations for studying autophagy are covered here. We have also highlighted the significance of network analysis and machine learning in capturing the core regulatory machinery governing the autophagy process. We explored the available autophagic databases and related resources along with their attributes that are useful in investigating autophagy through computational methods. We conclude the article addressing the potential future perspective in this area, which might provide a more in-depth insight into the dynamics of autophagy.

A Multi-Scale Approach to Modeling E. coli Chemotaxis

The degree to which we can understand the multi-scale organization of cellular life is tied to how well our models can represent this organization and the processes that drive its evolution. This paper uses Vivarium-an engine for composing heterogeneous computational biology models into integrated, multi-scale simulations. Vivarium's approach is demonstrated by combining several sub-models of biophysical processes into a model of chemotactic E. coli that exchange molecules with their environment, express the genes required for chemotaxis, swim, grow, and divide. This model is developed incrementally, highlighting cross-compartment mechanisms that link E. coli to its environment, with models for: (1) metabolism and transport, with transport moving nutrients across the membrane boundary and metabolism converting them to useful metabolites, (2) transcription, translation, complexation, and degradation, with stochastic mechanisms that read real gene sequence data and consume base pairs and ATP to make proteins and complexes, and (3) the activity of flagella and chemoreceptors, which together support navigation in the environment.

An on-lattice agent-based Monte Carlo model simulating the growth kinetics of multicellular tumor spheroids

PURPOSE

To develop an on-lattice agent-based model describing the growth of multicellular tumor spheroids using simple Monte Carlo tools.

METHODS

Cells are situated on the vertices of a cubic grid. Different cell states (proliferative, hypoxic or dead) and cell evolution rules, driven by 10 parameters, and the effects of the culture medium are included. About twenty spheroids of MCF-7 human breast cancer were cultivated and the experimental data were used for tuning the model parameters.

RESULTS

Simulated spheroids showed adequate sizes of the necrotic nuclei and of the hypoxic and proliferative cell phases as a function of the growth time, mimicking the overall characteristics of the experimental spheroids. The relation between the radii of the necrotic nucleus and the whole spheroid obtained in the simulations was similar to the experimental one and the number of cells, as a function of the spheroid volume, was well reproduced. The statistical variability of the Monte Carlo model described the whole volume range observed for the experimental spheroids. Assuming that the model parameters vary within Gaussian distributions it was obtained a sample of spheroids that reproduced much better the experimental findings.

CONCLUSIONS

The model developed allows describing the growth of in vitro multicellular spheroids and the experimental variability can be well reproduced. Its flexibility permits to vary both the agents involved and the rules that govern the spheroid growth. More general situations, such as, e. g., tumor vascularization, radiotherapy effects on solid tumors, or the validity of the tumor growth mathematical models can be studied.

From indication to decision: A hierarchical approach to model the chemotactic behavior of Escherichia coli

Reducing the complex behavior of living entities to its underlying physical and chemical processes is a formidable task in biology. Complex behaviors can be characterized as decision making: the ability to process the incoming information via an intracellular network and act upon this information to choose appropriate strategies. Motility is one such behavior that has been the focus many modeling efforts in the past. Our aim is to reduce the chemotactic behavior in Escherichia coli to its molecular constituents in order to paint a comprehensive and end-to-end picture of this intricate behavior. We utilize a hierarchical approach, consisting of three layers, to achieve this goal: at the first level, chemical reactions involved in chemotaxis are simulated. In the second level, the chemical reactions give rise to the mechanical movement of six independent flagella. At the last layer, the two lower layers are combined to allow a digital bacterium to receive information from its environment and swim through it with verve. Our results are in concert with the experimental studies concerning the motility of E.coli cells. In addition, we show that our detailed model of chemotaxis is reducible to a non-homogeneous Markov process.

Modelling bacterial chemotaxis for indirectly binding attractants

In bacterial chemotaxis, chemoattractant molecules may bind either directly or indirectly with receptors within the cell periplasmic space. The indirect binding mechanism, which involves an intermediate periplasmic binding protein, has been reported to increase sensitivity to dilute attractant concentrations as well as range of response. Current mathematical models for bacterial chemotaxis at the population scale do not appear to take the periplasmic binding protein (BP) concentration or the indirect binding mechanics into account. We formulate an indirect binding extension to the existing Rivero equation for chemotactic velocity based on fundamental reversible enzyme kinetics. The formulated indirect binding expression accounts for the periplasmic BP concentration and the dissociation constants for binding between attractant and periplasmic BP, as well as between BP and chemoreceptor. We validate the indirect-binding model using capillary assay simulations of the chemotactic responses of E. coli to the indirectly-binding attractants maltose and AI-2. The predicted response agrees well with experimental data from a number of maltose capillary assay studies conducted in previous literature. The model is also able to achieve good agreement with AI-2 capillary assay data of one study out of two tested. The chemotactic response of E. coli towards AI-2 appears to be of higher complexity due to reports of variable periplasmic BP concentration as well as the low concentration of periplasmic BP relative to the total receptor concentration. Our current model is thus suitable for indirect binding chemotactic response systems with constant periplasmic BP concentration that is significantly larger than the total receptor concentration, such as the response of E. coli towards maltose. Further considerations may be taken into account to model the chemotactic response towards AI-2 with greater accuracy.

On the efficiency of chemotactic pursuit - Comparing blind search with temporal and spatial gradient sensing

In chemotaxis, cells are modulating their migration patterns in response to concentration gradients of a guiding substance. Immune cells are believed to use such chemotactic sensing for remotely detecting and homing in on pathogens. Considering that immune cells may encounter a multitude of targets with vastly different migration properties, ranging from immobile to highly mobile, it is not clear which strategies of chemotactic pursuit are simultaneously efficient and versatile. We tackle this problem theoretically and define a tunable response function that maps temporal or spatial concentration gradients to migration behavior. The seven free parameters of this response function are optimized numerically with the objective of maximizing search efficiency against a wide spectrum of target cell properties. Finally, we reverse-engineer the best-performing parameter sets to uncover strategies of chemotactic pursuit that are efficient under different biologically realistic boundary conditions. Although strategies based on the temporal or spatial sensing of chemotactic gradients are significantly more efficient than unguided migration, such ‘blind search’ turns out to work surprisingly well, in particular if the immune cells are fast and directionally persistent. The resulting simulated data can be used for the design of chemotaxis experiments and for the development of algorithms that automatically detect and quantify goal oriented behavior in measured immune cell trajectories.

GutLogo: Agent-based modeling framework to investigate spatial and temporal dynamics in the gut microbiome

Knowledge of the spatial and temporal dynamics of the gut microbiome is essential to understanding the state of human health, as over a hundred diseases have been correlated with changes in microbial populations. Unfortunately, due to the complexity of the microbiome and the limitations of in vivo and in vitro experiments, studying spatial and temporal dynamics of gut bacteria in a biological setting is extremely challenging. Thus, in silico experiments present an excellent alternative for studying such systems. In consideration of these issues, we have developed a user-friendly agent-based model, GutLogo, that captures the spatial and temporal development of four representative bacterial genera populations in the ileum. We demonstrate the utility of this model by simulating population responses to perturbations in flow rate, nutrition, and probiotics. While our model predicts distinct changes in population levels due to these perturbations, most of the simulations suggest that the gut populations will return to their original steady states once the disturbance is removed. We hope that, in the future, the GutLogo model is utilized and customized by interested parties, as GutLogo can serve as a basic modeling framework for simulating a variety of physiological scenarios and can be extended to capture additional complexities of interest.

High-Performance Agent-Based Modeling Applied to Vocal Fold Inflammation and Repair

Fast and accurate computational biology models offer the prospect of accelerating the development of personalized medicine. A tool capable of estimating treatment success can help prevent unnecessary and costly treatments and potential harmful side effects. A novel high-performance Agent-Based Model (ABM) was adopted to simulate and visualize multi-scale complex biological processes arising in vocal fold inflammation and repair. The computational scheme was designed to organize the 3D ABM sub-tasks to fully utilize the resources available on current heterogeneous platforms consisting of multi-core CPUs and many-core GPUs. Subtasks are further parallelized and convolution-based diffusion is used to enhance the performance of the ABM simulation. The scheme was implemented using a client-server protocol allowing the results of each iteration to be analyzed and visualized on the server (i.e., in-situ) while the simulation is running on the same server. The resulting simulation and visualization software enables users to interact with and steer the course of the simulation in real-time as needed. This high-resolution 3D ABM framework was used for a case study of surgical vocal fold injury and repair. The new framework is capable of completing the simulation, visualization and remote result delivery in under 7 s per iteration, where each iteration of the simulation represents 30 min in the real world. The case study model was simulated at the physiological scale of a human vocal fold. This simulation tracks 17 million biological cells as well as a total of 1.7 billion signaling chemical and structural protein data points. The visualization component processes and renders all simulated biological cells and 154 million signaling chemical data points. The proposed high-performance 3D ABM was verified through comparisons with empirical vocal fold data. Representative trends of biomarker predictions in surgically injured vocal folds were observed.

Behavioral Variability and Phenotypic Diversity in Bacterial Chemotaxis

Living cells detect and process external signals using signaling pathways that are affected by random fluctuations. These variations cause the behavior of individual cells to fluctuate over time (behavioral variability) and generate phenotypic differences between genetically identical individuals (phenotypic diversity). These two noise sources reduce our ability to predict biological behavior because they diversify cellular responses to identical signals. Here, we review recent experimental and theoretical advances in understanding the mechanistic origin and functional consequences of such variation in Escherichia coli chemotaxis-a well-understood model of signal transduction and behavior. After briefly summarizing the architecture and logic of the chemotaxis system, we discuss determinants of behavior and chemotactic performance of individual cells. Then, we review how cell-to-cell differences in protein abundance map onto differences in individual chemotactic abilities and how phenotypic variability affects the performance of the population. We conclude with open questions to be addressed by future research.

Computational Systems Biology in Cancer: Modeling Methods and Applications

In recent years it has become clear that carcinogenesis is a complex process, both at the molecular and cellular levels. Understanding the origins, growth and spread of cancer, therefore requires an integrated or system-wide approach. Computational systems biology is an emerging sub-discipline in systems biology that utilizes the wealth of data from genomic, proteomic and metabolomic studies to build computer simulations of intra and intercellular processes. Several useful descriptive and predictive models of the origin, growth and spread of cancers have been developed in an effort to better understand the disease and potential therapeutic approaches. In this review we describe and assess the practical and theoretical underpinnings of commonly-used modeling approaches, including ordinary and partial differential equations, petri nets, cellular automata, agent based models and hybrid systems. A number of computer-based formalisms have been implemented to improve the accessibility of the various approaches to researchers whose primary interest lies outside of model development. We discuss several of these and describe how they have led to novel insights into tumor genesis, growth, apoptosis, vascularization and therapy.

Multi-fractal characterization of bacterial swimming dynamics: a case study on real and simulated Serratia marcescens

To add to the current state of knowledge about bacterial swimming dynamics, in this paper, we study the fractal swimming dynamics of populations of Serratia marcescens bacteria both in vitro and in silico, while accounting for realistic conditions like volume exclusion, chemical interactions, obstacles and distribution of chemoattractant in the environment. While previous research has shown that bacterial motion is non-ergodic, we demonstrate that, besides the non-ergodicity, the bacterial swimming dynamics is multi-fractal in nature. Finally, we demonstrate that the multi-fractal characteristic of bacterial dynamics is strongly affected by bacterial density and chemoattractant concentration.

Agent-based modelling in synthetic biology

Biological systems exhibit complex behaviours that emerge at many different levels of organization. These span the regulation of gene expression within single cells to the use of quorum sensing to co-ordinate the action of entire bacterial colonies. Synthetic biology aims to make the engineering of biology easier, offering an opportunity to control natural systems and develop new synthetic systems with useful prescribed behaviours. However, in many cases, it is not understood how individual cells should be programmed to ensure the emergence of a required collective behaviour. Agent-based modelling aims to tackle this problem, offering a framework in which to simulate such systems and explore cellular design rules. In this article, I review the use of agent-based models in synthetic biology, outline the available computational tools, and provide details on recently engineered biological systems that are amenable to this approach. I further highlight the challenges facing this methodology and some of the potential future directions.

Enhancement of Swimming Speed Leads to a More-Efficient Chemotactic Response to Repellent

Negative chemotaxis refers to the motion of microorganisms away from regions with high concentrations of chemorepellents. In this study, we set controlled gradients of NiCl2, a chemorepellent, in microchannels to quantify the motion of Escherichia coli over a broad range of concentrations. The experimental technique measured the motion of the bacteria in space and time and further related the motion to the local concentration profile of the repellent. Results show that the swimming speed of bacteria increases with an increasing concentration of repellent, which in turn enhances the drift velocity. The contribution of the increased swimming speed to the total drift velocity was in the range of 20 to 40%, with the remaining contribution coming from the modulation of the tumble frequency. A simple model that incorporates receptor dynamics, including adaptation, intracellular signaling, and swimming speed variation, was able to qualitatively capture the observed trend in drift velocity.

Understanding the link between single cell and population scale responses of Escherichia coli in differing ligand gradients

We formulate an agent-based population model of Escherichia coli cells which incorporates a description of the chemotaxis signalling cascade at the single cell scale. The model is used to gain insight into the link between the signalling cascade dynamics and the overall population response to differing chemoattractant gradients. Firstly, we consider how the observed variation in total (phosphorylated and unphosphorylated) signalling protein concentration affects the ability of cells to accumulate in differing chemoattractant gradients. Results reveal that a variation in total cell protein concentration between cells may be a mechanism for the survival of cell colonies across a wide range of differing environments. We then study the response of cells in the presence of two different chemoattractants. In doing so we demonstrate that the population scale response depends not on the absolute concentration of each chemoattractant but on the sensitivity of the chemoreceptors to their respective concentrations. Our results show the clear link between single cell features and the overall environment in which cells reside.

Established Microbial Colonies Can Survive Type VI Secretion Assault

Type VI secretion (T6S) is a cell-to-cell injection system that can be used as a microbial weapon. T6S kills vulnerable cells, and is present in close to 25% of sequenced Gram-negative bacteria. To examine the ecological role of T6S among bacteria, we competed self-immune T6S+ cells and T6S-sensitive cells in simulated range expansions. As killing takes place only at the interface between sensitive and T6S+ strains, while growth takes place everywhere, sufficiently large domains of sensitive cells can achieve net growth in the face of attack. Indeed T6S-sensitive cells can often outgrow their T6S+ competitors. We validated these findings through in vivo competition experiments between T6S+ Vibrio cholerae and T6S-sensitive Escherichia coli. We found that E. coli can survive and even dominate so long as they have an adequate opportunity to form microcolonies at the outset of the competition. Finally, in simulated competitions between two equivalent and mutually sensitive T6S+ strains, the more numerous strain has an advantage that increases with the T6S attack rate. We conclude that sufficiently large domains of T6S-sensitive individuals can survive attack and potentially outcompete self-immune T6S+ bacteria.

Modeling evolution of spatially distributed bacterial communities: a simulation with the haploid evolutionary constructor

Background

Multiscale approaches for integrating submodels of various levels of biological organization into a single model became the major tool of systems biology. In this paper, we have constructed and simulated a set of multiscale models of spatially distributed microbial communities and study an influence of unevenly distributed environmental factors on the genetic diversity and evolution of the community members.

Results

Haploid Evolutionary Constructor software http://evol-constructor.bionet.nsc.ru/ was expanded by adding the tool for the spatial modeling of a microbial community (1D, 2D and 3D versions). A set of the models of spatially distributed communities was built to demonstrate that the spatial distribution of cells affects both intensity of selection and evolution rate.

Conclusion

In spatially heterogeneous communities, the change in the direction of the environmental flow might be reflected in local irregular population dynamics, while the genetic structure of populations (frequencies of the alleles) remains stable. Furthermore, in spatially heterogeneous communities, the chemotaxis might dramatically affect the evolution of community members.

Fast, high-throughput measurement of collective behaviour in a bacterial population

Swimming bacteria explore their environment by performing a random walk, which is biased in response to, for example, chemical stimuli, resulting in a collective drift of bacterial populations towards 'a better life'. This phenomenon, called chemotaxis, is one of the best known forms of collective behaviour in bacteria, crucial for bacterial survival and virulence. Both single-cell and macroscopic assays have investigated bacterial behaviours. However, theories that relate the two scales have previously been difficult to test directly. We present an image analysis method, inspired by light scattering, which measures the average collective motion of thousands of bacteria simultaneously. Using this method, a time-varying collective drift as small as 50 nm s(-1) can be measured. The method, validated using simulations, was applied to chemotactic Escherichia coli bacteria in linear gradients of the attractant α-methylaspartate. This enabled us to test a coarse-grained minimal model of chemotaxis. Our results clearly map the onset of receptor methylation, and the transition from linear to logarithmic sensing in the bacterial response to an external chemoeffector. Our method is broadly applicable to problems involving the measurement of collective drift with high time resolution, such as cell migration and fluid flows measurements, and enables fast screening of tactic behaviours.

Methodologies for the modeling and simulation of biochemical networks, illustrated for signal transduction pathways: a primer

Biochemical networks depict the chemical interactions that take place among elements of living cells. They aim to elucidate how cellular behavior and functional properties of the cell emerge from the relationships between its components, i.e. molecules. Biochemical networks are largely characterized by dynamic behavior, and exhibit high degrees of complexity. Hence, the interest in such networks is growing and they have been the target of several recent modeling efforts. Signal transduction pathways (STPs) constitute a class of biochemical networks that receive, process, and respond to stimuli from the environment, as well as stimuli that are internal to the organism. An STP consists of a chain of intracellular signaling processes that ultimately result in generating different cellular responses. This primer presents the methodologies used for the modeling and simulation of biochemical networks, illustrated for STPs. These methodologies range from qualitative to quantitative, and include structural as well as dynamic analysis techniques. We describe the different methodologies, outline their underlying assumptions, and provide an assessment of their advantages and disadvantages. Moreover, publicly and/or commercially available implementations of these methodologies are listed as appropriate. In particular, this primer aims to provide a clear introduction and comprehensive coverage of biochemical modeling and simulation methodologies for the non-expert, with specific focus on relevant literature of STPs.

A microfluidic digital single-cell assay for the evaluation of anticancer drugs

Digital single-cell assays hold high potentials for the analysis of cell apoptosis and the evaluation of chemotherapeutic reagents for cancer therapy. In this paper, a microfluidic hydrodynamic trapping system was developed for digital single-cell assays with the capability of monitoring cellular dynamics over time. The microfluidic chip was designed with arrays of bypass structures for trapping individual cells without the need for surface modification, external electric force, or robotic equipment. After optimization of the bypass structure by both numerical simulations and experiments, a single-cell trapping efficiency of ∼90 % was achieved. We demonstrated the method as a digital single-cell assay for the evaluation of five clinically established chemotherapeutic reagents. As a result, the half maximal inhibitory concentration (IC50) values of these compounds could be conveniently determined. We further modeled the gradual decrease of active drugs over time which was often observed in vivo after an injection to investigate cell apoptosis against chemotherapeutic reagents. The developed method provided a valuable means for cell apoptotic analysis and evaluation of anticancer drugs.

Cooperative target tracking by a mobile bionanosensor network

This paper describes a mobile bionanosensor network designed for target tracking. The mobile bionanosensor network is composed of bacterium-based autonomous biosensors that coordinate their movement through the use of two types of signaling molecules, repellents and attractants. In search of a target, the bacterium-based autonomous biosensors release repellents to quickly spread over the environment, while, upon detecting a target, they release attractants to recruit other biosensors in the environment toward the location around the target. A mobility model of bacterium-based autonomous biosensors is first developed based on the rotational diffusion model of bacterial chemotaxis, and from this their collective movement to track a moving target is demonstrated. In simulation experiments, the mobile bionanosensor network is evaluated based on the mean tracking time. Simulation results show a set of parameter values that can optimize the mean tracking time, providing an insight into how bacterium-based autonomous biosensors may be designed and engineered for target tracking.

Toward an integrated software platform for systems pharmacology

Understanding complex biological systems requires the extensive support of computational tools. This is particularly true for systems pharmacology, which aims to understand the action of drugs and their interactions in a systems context. Computational models play an important role as they can be viewed as an explicit representation of biological hypotheses to be tested. A series of software and data resources are used for model development, verification and exploration of the possible behaviors of biological systems using the model that may not be possible or not cost effective by experiments. Software platforms play a dominant role in creativity and productivity support and have transformed many industries, techniques that can be applied to biology as well. Establishing an integrated software platform will be the next important step in the field. © 2013 The Authors. Biopharmaceutics & Drug Disposition published by John Wiley & Sons, Ltd.

Escherichia coli swimming is robust against variations in flagellar number

Bacterial chemotaxis is a paradigm for how environmental signals modulate cellular behavior. Although the network underlying this process has been studied extensively, we do not yet have an end-to-end understanding of chemotaxis. Specifically, how the rotational states of a cell’s flagella cooperatively determine whether the cell ‘runs’ or ‘tumbles’ remains poorly characterized. Here, we measure the swimming behavior of individual E. coli cells while simultaneously detecting the rotational states of each flagellum. We find that a simple mathematical expression relates the cell’s run/tumble bias to the number and average rotational state of its flagella. However, due to inter-flagellar correlations, an ‘effective number’ of flagella—smaller than the actual number—enters into this relation. Data from a chemotaxis mutant and stochastic modeling suggest that fluctuations of the regulator CheY-P are the source of flagellar correlations. A consequence of inter-flagellar correlations is that run/tumble behavior is only weakly dependent on number of flagella.

DOI:http://dx.doi.org/10.7554/eLife.01916.001

Escherichia coli is a rod-shaped bacterium commonly found in the lower intestines of humans and other warm-blooded animals. While most strains of E. coli are harmless, including most of those found in the human gut, some can cause diseases such as food poisoning. Due to its close association with humans and the fact that it is easy to grow and work with in the laboratory, E. coli has been intensively studied for over 60 years.

Many bacteria are capable of ‘swimming’ by using one or more flagella. These rotating whip-like structures are each driven by a reversible motor, and they act a bit like a propeller on a boat. While some bacteria have only a single flagellum, others, such as E. coli, have multiple flagella distributed over the cell surface. Rotating all their flagella in a counterclockwise direction allows the bacterium to swim—and it has been proposed that the clockwise movement of at least one flagellum will cause the bacterium cell to stop swimming and start tumbling.

E. coli is able to control the time it spends swimming or tumbling to move towards a nutrient, such as glucose, or away from certain harmful chemicals. However, the details of how the number of flagella and the direction of rotation (clockwise or counterclockwise) influence the motion of the bacterium are not fully understood.

Now, Mears et al. have used ‘optical tweezers’ to immobilize individual E. coli cells under a microscope, and then track both their swimming behavior and the movements of their flagella. This revealed that the individual flagella on the same cell tend to move in a coordinated way. Therefore, whilst tumbling could be caused by a single flagellum stopping swimming behavior, it often involved a concerted effort by many of the cell’s flagella.

After observing that E. coli cells with more flagella spent less time tumbling than would be predicted if a single flagella always ‘vetoed’ swimming, Mears et al. propose a new mathematical relationship between the number of flagella on the cell, the direction of rotation, and the resulting probability that the cell will tumble. This work shows that swimming behavior in bacteria is less affected by variations in the number of flagella than previously thought—and this phenomenon may provide evolutionary advantages to E. coli. The next step is to explore the mechanism by which bacteria coordinate their flagella.

DOI:http://dx.doi.org/10.7554/eLife.01916.002

Evaluation of bacterial motility from non-Gaussianity of finite-sample trajectories using the large deviation principle

Motility of bacteria is usually recognized in the trajectory data and compared with Brownian motion, but the diffusion coefficient is insufficient to evaluate it. In this paper, we propose a method based on the large deviation principle. We show that it can be used to evaluate the non-Gaussian characteristics of model Escherichia coli motions and to distinguish combinations of the mean running duration and running speed that lead to the same diffusion coefficient. Our proposed method does not require chemical stimuli to induce the chemotaxis in a specific direction, and it is applicable to various types of self-propelling motions for which no a priori information of, for example, threshold parameters for run and tumble or head/tail direction is available. We also address the issue of the finite-sample effect on the large deviation quantities, but we propose to make use of it to characterize the nature of motility.

Investigating biocomplexity through the agent-based paradigm

Capturing the dynamism that pervades biological systems requires a computational approach that can accommodate both the continuous features of the system environment as well as the flexible and heterogeneous nature of component interactions. This presents a serious challenge for the more traditional mathematical approaches that assume component homogeneity to relate system observables using mathematical equations. While the homogeneity condition does not lead to loss of accuracy while simulating various continua, it fails to offer detailed solutions when applied to systems with dynamically interacting heterogeneous components. As the functionality and architecture of most biological systems is a product of multi-faceted individual interactions at the sub-system level, continuum models rarely offer much beyond qualitative similarity. Agent-based modelling is a class of algorithmic computational approaches that rely on interactions between Turing-complete finite-state machines--or agents--to simulate, from the bottom-up, macroscopic properties of a system. In recognizing the heterogeneity condition, they offer suitable ontologies to the system components being modelled, thereby succeeding where their continuum counterparts tend to struggle. Furthermore, being inherently hierarchical, they are quite amenable to coupling with other computational paradigms. The integration of any agent-based framework with continuum models is arguably the most elegant and precise way of representing biological systems. Although in its nascence, agent-based modelling has been utilized to model biological complexity across a broad range of biological scales (from cells to societies). In this article, we explore the reasons that make agent-based modelling the most precise approach to model biological systems that tend to be non-linear and complex.

Multi-agent model of hepatitis C virus infection

OBJECTIVES

The objective of this study is to design a method for modeling hepatitis C virus (HCV) infection using multi-agent simulation and to verify it in practice.

METHODS AND MATERIALS

In this paper, first, the modeling of HCV infection using a multi-agent system is compared with the most commonly used model type, which is based on differential equations. Then, the implementation and results of the model using a multi-agent simulation is presented. To find the values of the parameters used in the model, a method using inverted simulation flow and genetic algorithm is proposed. All of the data regarding HCV infection are taken from the paper describing the model based on the differential equation to which the proposed method is compared.

RESULTS

Important advantages of the proposed method are noted and demonstrated: these include flexibility, clarity, re-usability and the possibility to model more complex dependencies. Then, the simulation framework that uses the proposed approach is successfully implemented in C++ and is verified by comparing it to the approach based on differential equations. The verification proves that an objective function that performs the best is the function that minimizes the maximal differences in the data. Finally, an analysis of one of the already known models is performed, and it is proved that it incorrectly models a decay in the hepatocytes number by 40%.

CONCLUSIONS

The proposed method has many advantages in comparison to the currently used model types and can be used successfully for analyzing HCV infection. With almost no modifications, it can also be used for other types of viral infections.

Data-driven quantification of the robustness and sensitivity of cell signaling networks

Robustness and sensitivity of responses generated by cell signaling networks has been associated with survival and evolvability of organisms. However, existing methods analyzing robustness and sensitivity of signaling networks ignore the experimentally observed cell-to-cell variations of protein abundances and cell functions or contain ad hoc assumptions. We propose and apply a data-driven maximum entropy based method to quantify robustness and sensitivity of Escherichia coli (E. coli) chemotaxis signaling network. Our analysis correctly rank orders different models of E. coli chemotaxis based on their robustness and suggests that parameters regulating cell signaling are evolutionary selected to vary in individual cells according to their abilities to perturb cell functions. Furthermore, predictions from our approach regarding distribution of protein abundances and properties of chemotactic responses in individual cells based on cell population averaged data are in excellent agreement with their experimental counterparts. Our approach is general and can be used to evaluate robustness as well as generate predictions of single cell properties based on population averaged experimental data in a wide range of cell signaling systems.

Biomimicry of quorum sensing using bacterial lifecycle model

BACKGROUND

Recent microbiologic studies have shown that quorum sensing mechanisms, which serve as one of the fundamental requirements for bacterial survival, exist widely in bacterial intra- and inter-species cell-cell communication. Many simulation models, inspired by the social behavior of natural organisms, are presented to provide new approaches for solving realistic optimization problems. Most of these simulation models follow population-based modelling approaches, where all the individuals are updated according to the same rules. Therefore, it is difficult to maintain the diversity of the population.

RESULTS

In this paper, we present a computational model termed LCM-QS, which simulates the bacterial quorum-sensing (QS) mechanism using an individual-based modelling approach under the framework of Agent-Environment-Rule (AER) scheme, i.e. bacterial lifecycle model (LCM). LCM-QS model can be classified into three main sub-models: chemotaxis with QS sub-model, reproduction and elimination sub-model and migration sub-model. The proposed model is used to not only imitate the bacterial evolution process at the single-cell level, but also concentrate on the study of bacterial macroscopic behaviour. Comparative experiments under four different scenarios have been conducted in an artificial 3-D environment with nutrients and noxious distribution. Detailed study on bacterial chemotatic processes with quorum sensing and without quorum sensing are compared. By using quorum sensing mechanisms, artificial bacteria working together can find the nutrient concentration (or global optimum) quickly in the artificial environment.

CONCLUSIONS

Biomimicry of quorum sensing mechanisms using the lifecycle model allows the artificial bacteria endowed with the communication abilities, which are essential to obtain more valuable information to guide their search cooperatively towards the preferred nutrient concentrations. It can also provide an inspiration for designing new swarm intelligence optimization algorithms, which can be used for solving the real-world problems.

A multi-paradigm modeling framework to simulate dynamic reciprocity in a bioreactor

Despite numerous technology advances, bioreactors are still mostly utilized as functional black-boxes where trial and error eventually leads to the desirable cellular outcome. Investigators have applied various computational approaches to understand the impact the internal dynamics of such devices has on overall cell growth, but such models cannot provide a comprehensive perspective regarding the system dynamics, due to limitations inherent to the underlying approaches. In this study, a novel multi-paradigm modeling platform capable of simulating the dynamic bidirectional relationship between cells and their microenvironment is presented. Designing the modeling platform entailed combining and coupling fully an agent-based modeling platform with a transport phenomena computational modeling framework. To demonstrate capability, the platform was used to study the impact of bioreactor parameters on the overall cell population behavior and vice versa. In order to achieve this, virtual bioreactors were constructed and seeded. The virtual cells, guided by a set of rules involving the simulated mass transport inside the bioreactor, as well as cell-related probabilistic parameters, were capable of displaying an array of behaviors such as proliferation, migration, chemotaxis and apoptosis. In this way the platform was shown to capture not only the impact of bioreactor transport processes on cellular behavior but also the influence that cellular activity wields on that very same local mass transport, thereby influencing overall cell growth. The platform was validated by simulating cellular chemotaxis in a virtual direct visualization chamber and comparing the simulation with its experimental analogue. The results presented in this paper are in agreement with published models of similar flavor. The modeling platform can be used as a concept selection tool to optimize bioreactor design specifications.

Agent-based models of cellular systems

Software agents are particularly suitable for engineering models and simulations of cellular systems. In a very natural and intuitive manner, individual software components are therein delegated to reproduce "in silico" the behavior of individual components of alive systems at a given level of resolution. Individuals' actions and interactions among individuals allow complex collective behavior to emerge. In this chapter we first introduce the readers to software agents and multi-agent systems, reviewing the evolution of agent-based modeling of biomolecular systems in the last decade. We then describe the main tools, platforms, and methodologies available for programming societies of agents, possibly profiting also of toolkits that do not require advanced programming skills.

BSim: an agent-based tool for modeling bacterial populations in systems and synthetic biology

Large-scale collective behaviors such as synchronization and coordination spontaneously arise in many bacterial populations. With systems biology attempting to understand these phenomena, and synthetic biology opening up the possibility of engineering them for our own benefit, there is growing interest in how bacterial populations are best modeled. Here we introduce BSim, a highly flexible agent-based computational tool for analyzing the relationships between single-cell dynamics and population level features. BSim includes reference implementations of many bacterial traits to enable the quick development of new models partially built from existing ones. Unlike existing modeling tools, BSim fully considers spatial aspects of a model allowing for the description of intricate micro-scale structures, enabling the modeling of bacterial behavior in more realistic three-dimensional, complex environments. The new opportunities that BSim opens are illustrated through several diverse examples covering: spatial multicellular computing, modeling complex environments, population dynamics of the lac operon, and the synchronization of genetic oscillators. BSim is open source software that is freely available from http://bsim-bccs.sf.net and distributed under the Open Source Initiative (OSI) recognized MIT license. Developer documentation and a wide range of example simulations are also available from the website. BSim requires Java version 1.6 or higher.

Stochastic coordination of multiple actuators reduces latency and improves chemotactic response in bacteria

Individual neuronal, signal transduction, and regulatory pathways often control multiple stochastic downstream actuators, which raises the question of how coordinated response to a single input can be achieved when individual actuators fluctuate independently. In Escherichia coli, the bacterial chemotaxis pathway controls the activity of multiple flagellar motors to generate the run-and-tumble motion of the cell. High-resolution microscopy experiments have identified the key conformational changes adopted by individual flagella during this process. By incorporating these observations into a stochastic model of the flagellar bundle, we demonstrate that the presence of multiple motors imposes a trade-off on chemotactic performance. Multiple motors reduce the latency of the response below the time scale of the stochastic switching of a single motor, which improves performance on steep gradients of attractants. However, the uncoordinated switching of multiple motors interrupts and shortens cell runs, which thereby reduces signal detection and performance on shallow gradients. Remarkably, when slow fluctuations generated by the adaptation mechanism of the chemotaxis system are incorporated in the model at levels measured in experiments, the chemotactic sensitivity and performance in shallow gradients is partially restored with marginal effects for steep gradients. The noise is beneficial because it simultaneously generates long events in the statistics of individual motors and coordinates the motors to generate a long tail in the run length distribution of the cell. Occasional long runs are known to enhance exploration of random walkers. Here we show that they have the additional benefit of enhancing the sensitivity of the bacterium to very shallow gradients.

Modeling formalisms in Systems Biology

Systems Biology has taken advantage of computational tools and high-throughput experimental data to model several biological processes. These include signaling, gene regulatory, and metabolic networks. However, most of these models are specific to each kind of network. Their interconnection demands a whole-cell modeling framework for a complete understanding of cellular systems. We describe the features required by an integrated framework for modeling, analyzing and simulating biological processes, and review several modeling formalisms that have been used in Systems Biology including Boolean networks, Bayesian networks, Petri nets, process algebras, constraint-based models, differential equations, rule-based models, interacting state machines, cellular automata, and agent-based models. We compare the features provided by different formalisms, and discuss recent approaches in the integration of these formalisms, as well as possible directions for the future.

Software for systems biology: from tools to integrated platforms

Understanding complex biological systems requires extensive support from software tools. Such tools are needed at each step of a systems biology computational workflow, which typically consists of data handling, network inference, deep curation, dynamical simulation and model analysis. In addition, there are now efforts to develop integrated software platforms, so that tools that are used at different stages of the workflow and by different researchers can easily be used together. This Review describes the types of software tools that are required at different stages of systems biology research and the current options that are available for systems biology researchers. We also discuss the challenges and prospects for modelling the effects of genetic changes on physiology and the concept of an integrated platform.

Toward a comprehensive language for biological systems

Rule-based modeling has become a powerful approach for modeling intracellular networks, which are characterized by rich molecular diversity. Truly comprehensive models of cell behavior, however, must address spatial complexity at both the intracellular level and at the level of interacting populations of cells, and will require richer modeling languages and tools. A recent paper in BMC Systems Biology represents a signifcant step toward the development of a unified modeling language and software platform for the development of multi-level, multiscale biological models.

See research article: http://www.biomedcentral.com/1752-0509/5/166

Multi-scale modeling in biology: how to bridge the gaps between scales?

Human physiological functions are regulated across many orders of magnitude in space and time. Integrating the information and dynamics from one scale to another is critical for the understanding of human physiology and the treatment of diseases. Multi-scale modeling, as a computational approach, has been widely adopted by researchers in computational and systems biology. A key unsolved issue is how to represent appropriately the dynamical behaviors of a high-dimensional model of a lower scale by a low-dimensional model of a higher scale, so that it can be used to investigate complex dynamical behaviors at even higher scales of integration. In the article, we first review the widely-used different modeling methodologies and their applications at different scales. We then discuss the gaps between different modeling methodologies and between scales, and discuss potential methods for bridging the gaps between scales.

Discrete diffusion models to study the effects of Mg2+ concentration on the PhoPQ signal transduction system

Background

The challenge today is to develop a modeling and simulation paradigm that integrates structural, molecular and genetic data for a quantitative understanding of physiology and behavior of biological processes at multiple scales. This modeling method requires techniques that maintain a reasonable accuracy of the biological process and also reduces the computational overhead. This objective motivates the use of new methods that can transform the problem from energy and affinity based modeling to information theory based modeling. To achieve this, we transform all dynamics within the cell into a random event time, which is specified through an information domain measure like probability distribution. This allows us to use the “in silico” stochastic event based modeling approach to find the molecular dynamics of the system.

Results

In this paper, we present the discrete event simulation concept using the example of the signal transduction cascade triggered by extra-cellular Mg²⁺ concentration in the two component PhoPQ regulatory system of Salmonella Typhimurium. We also present a model to compute the information domain measure of the molecular transport process by estimating the statistical parameters of inter-arrival time between molecules/ions coming to a cell receptor as external signal. This model transforms the diffusion process into the information theory measure of stochastic event completion time to get the distribution of the Mg²⁺ departure events. Using these molecular transport models, we next study the in-silico effects of this external trigger on the PhoPQ system.

Conclusions

Our results illustrate the accuracy of the proposed diffusion models in explaining the molecular/ionic transport processes inside the cell. Also, the proposed simulation framework can incorporate the stochasticity in cellular environments to a certain degree of accuracy. We expect that this scalable simulation platform will be able to model more complex biological systems with reasonable accuracy to understand their temporal dynamics.

Chemotaxis of Escherichia coli to L-serine

A novel experimental technique was used to quantify the motion of E. coli to varying serine concentrations and gradients so as to capture the spatial and temporal variation of the chemotactic response. The average run speed and the cell diffusivity are found to be dependent on the serine concentration. The measured diffusivities were in the range of 1.2-2.5 x 10 (-10) m(2) s(-1). The study revealed that the rotational diffusivity of the cells, induced by the extracellular environment, also varies with the serine concentration. The drift velocity increased with serine gradients reaching a maximum value of approximately 5.5 microm s(-1) at 1.6 microM microm(-1) after which it decreased. Experimental analysis demonstrated the interdependence of run speed, rotational diffusivity and drift velocity that characterizes the motion. Further, the motion was found to critically depend on the oxygen concentration and energy level of the cells.