Identification of age-structured models: cell cycle phase transitions

A comprehensive review of computational cell cycle models in guiding cancer treatment strategies

This article reviews the current knowledge and recent advancements in computational modeling of the cell cycle. It offers a comparative analysis of various modeling paradigms, highlighting their unique strengths, limitations, and applications. Specifically, the article compares deterministic and stochastic models, single-cell versus population models, and mechanistic versus abstract models. This detailed analysis helps determine the most suitable modeling framework for various research needs. Additionally, the discussion extends to the utilization of these computational models to illuminate cell cycle dynamics, with a particular focus on cell cycle viability, crosstalk with signaling pathways, tumor microenvironment, DNA replication, and repair mechanisms, underscoring their critical roles in tumor progression and the optimization of cancer therapies. By applying these models to crucial aspects of cancer therapy planning for better outcomes, including drug efficacy quantification, drug discovery, drug resistance analysis, and dose optimization, the review highlights the significant potential of computational insights in enhancing the precision and effectiveness of cancer treatments. This emphasis on the intricate relationship between computational modeling and therapeutic strategy development underscores the pivotal role of advanced modeling techniques in navigating the complexities of cell cycle dynamics and their implications for cancer therapy.

DNA-structured mathematical model of cell-cycle progression in cyclic hypoxia

New experimental data have shown how the periodic exposure of cells to low oxygen levels (i.e., cyclic hypoxia) impacts their progress through the cell-cycle. Cyclic hypoxia has been detected in tumours and linked to poor prognosis and treatment failure. While fluctuating oxygen environments can be reproduced in vitro, the range of oxygen cycles that can be tested is limited. By contrast, mathematical models can be used to predict the response to a wide range of cyclic dynamics. Accordingly, in this paper we develop a mechanistic model of the cell-cycle that can be combined with in vitro experiments, to better understand the link between cyclic hypoxia and cell-cycle dysregulation. A distinguishing feature of our model is the inclusion of impaired DNA synthesis and cell-cycle arrest due to periodic exposure to severely low oxygen levels. Our model decomposes the cell population into five compartments and a time-dependent delay accounts for the variability in the duration of the S phase which increases in severe hypoxia due to reduced rates of DNA synthesis. We calibrate our model against experimental data and show that it recapitulates the observed cell-cycle dynamics. We use the calibrated model to investigate the response of cells to oxygen cycles not yet tested experimentally. When the re-oxygenation phase is sufficiently long, our model predicts that cyclic hypoxia simply slows cell proliferation since cells spend more time in the S phase. On the contrary, cycles with short periods of re-oxygenation are predicted to lead to inhibition of proliferation, with cells arresting from the cell-cycle in the G2 phase. While model predictions on short time scales (about a day) are fairly accurate (i.e, confidence intervals are small), the predictions become more uncertain over longer periods. Hence, we use our model to inform experimental design that can lead to improved model parameter estimates and validate model predictions.

A checkpoint-oriented cell cycle simulation model

Modeling and in silico simulations are of major conceptual and applicative interest in studying the cell cycle and proliferation in eukaryotic cells. In this paper, we present a cell cycle checkpoint-oriented simulator that uses agent-based simulation modeling to reproduce the dynamics of a cancer cell population in exponential growth. Our in silico simulations were successfully validated by experimental in vitro supporting data obtained with HCT116 colon cancer cells. We demonstrated that this model can simulate cell confluence and the associated elongation of the G1 phase. Using nocodazole to synchronize cancer cells at mitosis, we confirmed the model predictivity and provided evidence of an additional and unexpected effect of nocodazole on the overall cell cycle progression. We anticipate that this cell cycle simulator will be a potential source of new insights and research perspectives.

Quantifying the length and variance of the eukaryotic cell cycle phases by a stochastic model and dual nucleoside pulse labelling

A fundamental property of cell populations is their growth rate as well as the time needed for cell division and its variance. The eukaryotic cell cycle progresses in an ordered sequence through the phases and and is regulated by environmental cues and by intracellular checkpoints. Reflecting this regulatory complexity, the length of each phase varies considerably in different kinds of cells but also among genetically and morphologically indistinguishable cells. This article addresses the question of how to describe and quantify the mean and variance of the cell cycle phase lengths. A phase-resolved cell cycle model is introduced assuming that phase completion times are distributed as delayed exponential functions, capturing the observations that each realization of a cycle phase is variable in length and requires a minimal time. In this model, the total cell cycle length is distributed as a delayed hypoexponential function that closely reproduces empirical distributions. Analytic solutions are derived for the proportions of cells in each cycle phase in a population growing under balanced growth and under specific non-stationary conditions. These solutions are then adapted to describe conventional cell cycle kinetic assays based on pulse labelling with nucleoside analogs. The model fits well to data obtained with two distinct proliferating cell lines labelled with a single bromodeoxiuridine pulse. However, whereas mean lengths are precisely estimated for all phases, the respective variances remain uncertain. To overcome this limitation, a redesigned experimental protocol is derived and validated in silico. The novelty is the timing of two consecutive pulses with distinct nucleosides that enables accurate and precise estimation of both the mean and the variance of the length of all phases. The proposed methodology to quantify the phase length distributions gives results potentially equivalent to those obtained with modern phase-specific biosensor-based fluorescent imaging.

Among the important characteristics of dividing cell populations is the time necessary for cells to complete each of the cell cycle phases, that is, to increase the cell's mass, to duplicate and repair its genome, to properly segregate its chromosomes, and to make decisions whether to continue dividing or enter a quiescent state. The cycle phase times also determine the maximal rate at which a dividing cell population can grow in size. Cell cycle phase completion times largely differ between cell types, cellular environments as well as metabolic stages, and can thus be considered as part of the phenotype of a given cell. Our article advances the methods to quantitatively characterize this phenotype. We introduce a novel phase-resolved cell cycle progression model and use it to estimate the mean and variance of the cycle phase completion times based on nucleoside analog pulse labelling experiments. This classic workhorse of cell cycle kinetic studies is revamped by our approach to potentially rival in accuracy and precision with modern phase-specific biosensor-based fluorescent imaging, while superseding the latter in its application scope.

Population Balance Modeling: Current Status and Future Prospects

Population balance modeling is undergoing phenomenal growth in its applications, and this growth is accompanied by multifarious reviews. This review aims to fortify the model's fundamental base, as well as point to a variety of new applications, including modeling of crystal morphology, cell growth and differentiation, gene regulatory processes, and transfer of drug resistance. This is accomplished by presenting the many faces of population balance equations that arise in the foregoing applications.

From single cells to microbial population dynamics: modelling in biotechnology based on measurements of individual cells

The development of dynamic modelling of microbial populations in bioprocesses is reviewed. In the 1960s Arnold Fredrickson established the theoretical basis of such models, and other researchers have subsequently advanced them. This review explores the relationships that describe cell proliferation and evaluates the importance of the application of flow cytometry to the fundamental parameterisation of the models for their use in bioprocess engineering. The section "Individual-Based Modelling" discusses recent theoretical developments. Delay-differential equations are demonstrated to describe accurately temporal variation of the cell proliferation cycle and specialised approaches and related iconography are applied to stochastic and deterministic modelling of stages of cellular development. Synchronised cultures of the bacterium Cupriavidus necator were prepared and monitored using a flow cytometer. The data obtained demonstrate that cell proliferation could be simulated quantitatively using the developed model.

Deconstructing stem cell population heterogeneity: single-cell analysis and modeling approaches

Isogenic stem cell populations display cell-to-cell variations in a multitude of attributes including gene or protein expression, epigenetic state, morphology, proliferation and proclivity for differentiation. The origins of the observed heterogeneity and its roles in the maintenance of pluripotency and the lineage specification of stem cells remain unclear. Addressing pertinent questions will require the employment of single-cell analysis methods as traditional cell biochemical and biomolecular assays yield mostly population-average data. In addition to time-lapse microscopy and flow cytometry, recent advances in single-cell genomic, transcriptomic and proteomic profiling are reviewed. The application of multiple displacement amplification, next generation sequencing, mass cytometry and spectrometry to stem cell systems is expected to provide a wealth of information affording unprecedented levels of multiparametric characterization of cell ensembles under defined conditions promoting pluripotency or commitment. Establishing connections between single-cell analysis information and the observed phenotypes will also require suitable mathematical models. Stem cell self-renewal and differentiation are orchestrated by the coordinated regulation of subcellular, intercellular and niche-wide processes spanning multiple time scales. Here, we discuss different modeling approaches and challenges arising from their application to stem cell populations. Integrating single-cell analysis with computational methods will fill gaps in our knowledge about the functions of heterogeneity in stem cell physiology. This combination will also aid the rational design of efficient differentiation and reprogramming strategies as well as bioprocesses for the production of clinically valuable stem cell derivatives.

Optimizing cancer pharmacotherapeutics using mathematical modeling and a systems biology approach

Research in mathematics and in mathematical biology on cancer and its treatments has been soaring in the past 10 years at an unprecedented speed. Such thriving is likely due as much to new findings in fundamental biology as to an emerging general interest from mathematicians and engineers towards applications in biology and medicine and to their subsequently designed representations and predictions of tumor processes that are now allowed by modern means of computation and simulation. This article, which does not claim the status of an extended review paper on mathematical models of cancer and its treatment, is focused on modeling in a systems biology perspective. I will list here the most necessary mathematical methods, in my opinion, which, while enforcing already existing methods, should be further developed towards designing and applying optimized individualized treatments of cancer in the clinic.

Experimental analysis and modelling of in vitro HUVECs proliferation in the presence of various types of drugs

OBJECTIVES

This study focuses on experimental analysis and corresponding mathematical simulation of in vitro HUVECs (human umbilical vein endothelial cells) proliferation in the presence of various types of drugs.

MATERIALS AND METHODS

HUVECs, once seeded in Petri dishes, were expanded to confluence. Temporal profiles of total count obtained by classic haemocytometry and cell size distribution measured using an electronic Coulter counter, are quantitatively simulated by a suitable model based on the population balance approach. Influence of drugs on cell proliferation is also properly simulated by accounting for suitable kinetic equations.

RESULTS AND DISCUSSION

The models' parameters have been determined by comparison with experimental data related to cell population expansion and cell size distribution in the absence of drugs. Inhibition constant for each type of drug has been estimated by comparing the experimental data with model results concerning temporal profiles of total cell count. The reliability of the model and its predictive capability have been tested by simulating cell size distribution for experiments performed in the presence of drugs. The proposed model will be useful in interpreting effects of selected drugs on expansion of readily available human cells.

Scalable stirred-suspension bioreactor culture of human pluripotent stem cells

Advances in stem cell biology have afforded promising results for the generation of various cell types for therapies against devastating diseases. However, a prerequisite for realizing the therapeutic potential of stem cells is the development of bioprocesses for the production of stem cell progeny in quantities that satisfy clinical demands. Recent reports on the expansion and directed differentiation of human embryonic stem cells (hESCs) in scalable stirred-suspension bioreactors (SSBs) demonstrated that large-scale production of therapeutically useful hESC progeny is feasible with current state-of-the-art culture technologies. Stem cells have been cultured in SSBs as aggregates, in microcarrier suspension and after encapsulation. The various modes in which SSBs can be employed for the cultivation of hESCs and human induced pluripotent stem cells (hiPSCs) are described. To that end, this is the first account of hiPSC cultivation in a microcarrier stirred-suspension system. Given that cultured stem cells and their differentiated progeny are the actual products used in tissue engineering and cell therapies, the impact of bioreactor's operating conditions on stem cell self-renewal and commitment should be considered. The effects of variables specific to SSB operation on stem cell physiology are discussed. Finally, major challenges are presented which remain to be addressed before the mainstream use of SSBs for the large-scale culture of hESCs and hiPSCs.

In vitro ovine articular chondrocyte proliferation: experiments and modelling

This study focuses on analysis of in vitro cultures of chondrocytes from ovine articular cartilage. Isolated cells were seeded in Petri dishes, then expanded to confluence and phenotypically characterized by flow cytometry. The sigmoidal temporal profile of total counts was obtained by classic haemocytometry and corresponding cell size distributions were measured electronically using a Coulter Counter. A mathematical model recently proposed (1) was adopted for quantitative interpretation of these experimental data. The model is based on a 1-D (that is, mass-structured), single-staged population balance approach capable of taking into account contact inhibition at confluence. The model's parameters were determined by fitting measured total cell counts and size distributions. Model reliability was verified by predicting cell proliferation counts and corresponding size distributions at culture times longer than those used when tuning the model's parameters. It was found that adoption of cell mass as the intrinsic characteristic of a growing chondrocyte population enables sigmoidal temporal profiles of total counts in the Petri dish, as well as cell size distributions at 'balanced growth', to be adequately predicted.

Experimental analysis and modelling of in vitro proliferation of mesenchymal stem cells

OBJECTIVES

Stem cell therapies based on differentiation of adult or embryonic stem cells into specialized ones appear to be effective for treating several human diseases. This work addresses the mathematical simulation of proliferation kinetics of stem cells.

MATERIALS AND METHODS

Sheep bone marrow mesenchymal stem cells (phenotype characterized by flow cytometry analysis) seeded at different initial concentrations in Petri dishes were expanded to confluence. Sigmoid temporal profiles of total counts obtained through classic haemocytometry were quantitatively interpreted by both a phenomenological logistic equation and a novel model based on a one-dimensional, single-staged population balance approach capable of taking into account contact inhibition at confluence. The models' parameters were determined by comparison with experimental data on population expansion starting from single seeding concentration. Reliability of the models was tested by predicting cell proliferation carried out starting from different seeding concentrations.

RESULTS AND DISCUSSION

It was found that the proposed population balance modelling approach was successful in predicting the experimental data over the whole range of initial cell numbers investigated, while prediction capability of phenomenological logistic equation was more limited.