Mathematical models for cell-matrix interactions during dermal wound healing

This paper contains a review of our recent work on the mathematical modeling of cell interaction with extracellular matrix components during the process of dermal wound healing. The models are of partial differential equation type and allow us to investigate in detail how various mechanochemical effects may be responsible for certain wound healing disorders such as fibrocontractive and fibroproliferative diseases. We also present a model for wound healing angiogenesis. The latter has several features in common with angiogenesis during cancer tumour growth and spread so a deeper understanding of the phenomenon in the context of wound healing may also help in the treatment of certain cancers

Computational and Mathematical Modelling of the EGF Receptor System

This chapter gives an overview of computational and mathematical modelling of the EGF receptor system. It begins with a survey of motivations for producing such models, then describes the main approaches that are taken to carrying out such modelling, viz. differential equations and individual-based modelling. Finally, a number of projects that applying modelling and simulation techniques to various aspects of the EGF receptor system are described

Systems biology in animal sciences

Systems biology is a rapidly expanding field of research and is applied in a number of biological disciplines. In animal sciences, omics approaches are increasingly used, yielding vast amounts of data, but systems biology approaches to extract understanding from these data of biological processes and animal traits are not yet frequently used. This paper aims to explain what systems biology is and which areas of animal sciences could benefit from systems biology approaches. Systems biology aims to understand whole biological systems working as a unit, rather than investigating their individual components. Therefore, systems biology can be considered a holistic approach, as opposed to reductionism. The recently developed ‘omics’ technologies enable biological sciences to characterize the molecular components of life with ever increasing speed, yielding vast amounts of data. However, biological functions do not follow from the simple addition of the properties of system components, but rather arise from the dynamic interactions of these components. Systems biology combines statistics, bioinformatics and mathematical modeling to integrate and analyze large amounts of data in order to extract a better understanding of the biology from these huge data sets and to predict the behavior of biological systems. A ‘system’ approach and mathematical modeling in biological sciences are not new in itself, as they were used in biochemistry, physiology and genetics long before the name systems biology was coined. However, the present combination of mass biological data and of computational and modeling tools is unprecedented and truly represents a major paradigm shift in biology. Significant advances have been made using systems biology approaches, especially in the field of bacterial and eukaryotic cells and in human medicine. Similarly, progress is being made with ‘system approaches’ in animal sciences, providing exciting opportunities to predict and modulate animal traits

Lattice and Continuum Modelling of a Bioactive Porous Tissue Scaffold

A contemporary procedure to grow artificial tissue is to seed cells onto a porous biomaterial scaffold and culture it within a perfusion bioreactor to facilitate the transport of nutrients to growing cells. Typical models of cell growth for tissue engineering applications make use of spatially homogeneous or spatially continuous equations to model cell growth, flow of culture medium, nutrient transport, and their interactions. The network structure of the physical porous scaffold is often incorporated through parameters in these models, either phenomenologically or through techniques like mathematical homogenization. We derive a model on a square grid lattice to demonstrate the importance of explicitly modelling the network structure of the porous scaffold, and compare results from this model with those from a modified continuum model from the literature. We capture two-way coupling between cell growth and fluid flow by allowing cells to block pores, and by allowing the shear stress of the fluid to affect cell growth and death. We explore a range of parameters for both models, and demonstrate quantitative and qualitative differences between predictions from each of these approaches, including spatial pattern formation and local oscillations in cell density present only in the lattice model. These differences suggest that for some parameter regimes, corresponding to specific cell types and scaffold geometries, the lattice model gives qualitatively different model predictions than typical continuum models. Our results inform model selection for bioactive porous tissue scaffolds, aiding in the development of successful tissue engineering experiments and eventually clinically successful technologies.Comment: 38 pages, 16 figures. This version includes a much-expanded introduction, and a new section on nonlinear diffusion in addition to polish throughou

Bridging the gap between the micro- and the macro-world of tumors

At present it is still quite difficult to match the vast knowledge on the behavior of individual tumor cells with macroscopic measurements on clinical tumors. On the modeling side, we already know how to deal with many molecular pathways and cellular events, using systems of differential equations and other modeling tools, and ideally, we should be able to extend such a mathematical description up to the level of large tumor masses. An extended model should thus help us forecast the behavior of large tumors from our basic knowledge of microscopic processes. Unfortunately, the complexity of these processes makes it very difficult -- probably impossible -- to develop comprehensive analytical models. We try to bridge the gap with a simulation program which is based on basic biochemical and biophysical processes -- thereby building an effective computational model -- and in this paper we describe its structure, endeavoring to make the description sufficiently detailed and yet understandable.Comment: 24 pages, 10 figures. Accepted for publication in AIP Advances, in the special issue on the physics of cance

Computational analysis of dynamic light exposure of unicellular algal cells in a flat-panel photobioreactor to support light-induced CO2 bioprocess development

Cyanobacterial cell factories trace a vibrant pathway to climate change neutrality and sustainable development owing to their ability to turn carbon dioxide-rich waste into a broad portfolio of renewable compounds, which are deemed valuable in green chemistry cross-sectorial applications. Cell factory design requires to define the optimal operational and cultivation conditions. The paramount parameter in biomass cultivation in photobioreactors is the light intensity since it impacts cellular physiology and productivity. Our modeling framework provides a basis for the predictive control of light-limited, light-saturated, and light-inhibited growth of the Synechocystis sp. PCC 6803 model organism in a flat-panel photobioreactor. The model here presented couples computational fluid dynamics, light transmission, kinetic modeling, and the reconstruction of single cell trajectories in differently irradiated areas of the photobioreactor to relate key physiological parameters to the multi-faceted processes occurring in the cultivation environment. Furthermore, our analysis highlights the need for properly constraining the model with decisive qualitative and quantitative data related to light calibration and light measurements both at the inlet and outlet of the photobioreactor in order to boost the accuracy and extrapolation capabilities of the model

A mass action model of a fibroblast growth factor signaling pathway and its simplification

We consider a kinetic law of mass action model for Fibroblast Growth Factor (FGF) signaling, focusing on the induction of the RAS-MAP kinase pathway via GRB2 binding. Our biologically simple model suffers a combinatorial explosion in the number of differential equations required to simulate the system. In addition to numerically solving the full model, we show that it can be accurately simplified. This requires combining matched asymptotics, the quasi-steady state hypothesis, and the fact subsets of the equations decouple asymptotically. Both the full and simplified models reproduce the qualitative dynamics observed experimentally and in previous stochastic models. The simplified model also elucidates both the qualitative features of GRB2 binding and the complex relationship between SHP2 levels, the rate SHP2 induces dephosphorylation and levels of bound GRB2. In addition to providing insight into the important and redundant features of FGF signaling, such work further highlights the usefulness of numerous simplification techniques in the study of mass action models of signal transduction, as also illustrated recently by Borisov and co-workers (Borisov et al. in Biophys. J. 89, 951–66, 2005, Biosystems 83, 152–66, 2006; Kiyatkin et al. in J. Biol. Chem. 281, 19925–9938, 2006). These developments will facilitate the construction of tractable models of FGF signaling, incorporating further biological realism, such as spatial effects or realistic binding stoichiometries, despite a more severe combinatorial explosion associated with the latter