Modeling cell population dynamics

Markov field network model of multi-modal data predicts effects of immune system perturbations on intravenous BCG vaccination in macaques

Analysis of multi-modal datasets can identify multi-scale interactions underlying biological systems but can be beset by spurious connections due to indirect impacts propagating through an unmapped biological network. For example, studies in macaques have shown that Bacillus Calmette-Guerin (BCG) vaccination by an intravenous route protects against tuberculosis, correlating with changes across various immune data modes. To eliminate spurious correlations and identify critical immune interactions in a public multi-modal dataset (systems serology, cytokines, and cytometry) of vaccinated macaques, we applied Markov fields (MFs), a data-driven approach that explains vaccine efficacy and immune correlations via multivariate network paths, without requiring large numbers of samples (i.e., macaques) relative to multivariate features. We find that integrating multiple data modes with MFs helps remove spurious connections. Finally, we used the MF to predict outcomes of perturbations at various immune nodes, including an experimentally validated B cell depletion that induced network-wide shifts without reducing vaccine protection.

Markov Field network model of multi-modal data predicts effects of immune system perturbations on intravenous BCG vaccination in macaques

Analysis of multi-modal datasets can identify multi-scale interactions underlying biological systems, but can be beset by spurious connections due to indirect impacts propagating through an unmapped biological network. For example, studies in macaques have shown that BCG vaccination by an intravenous route protects against tuberculosis, correlating with changes across various immune data modes. To eliminate spurious correlations and identify critical immune interactions in a public multi-modal dataset (systems serology, cytokines, cytometry) of vaccinated macaques, we applied Markov Fields (MF), a data-driven approach that explains vaccine efficacy and immune correlations via multivariate network paths, without requiring large numbers of samples (i.e. macaques) relative to multivariate features. Furthermore, we find that integrating multiple data modes with MFs helps to remove spurious connections. Finally, we used the MF to predict outcomes of perturbations at various immune nodes, including a B-cell depletion that induced network-wide shifts without reducing vaccine protection, which we validated experimentally.

Magnetic field platform for experiments on well-mixed and spatially structured microbial populations

Magnetic fields have been shown to affect sensing, migration, and navigation in living organisms. However, the effects of magnetic fields on microorganisms largely remain to be elucidated. We develop an open-source, 3D-printed magnetic field exposure device to perform experiments on well-mixed and spatially structured microbial populations. This device is designed in AutoCAD, modeled in COMSOL, and validated using a Gaussmeter and experiments on the budding yeast Saccharomyces cerevisiae. We find that static magnetic field exposure slows the spatially structured expansion of yeast mats that expand in two dimensions, but not yeast mats that expand in three dimensions, across the surface of semi-solid yeast extract-peptone-dextrose agar media. We also find that magnetic fields do not affect the growth of planktonic yeast cells in well-mixed liquid yeast extract-peptone-dextrose media. This study provides an adaptable device for performing controlled magnetic field experiments on microbes and advances our understanding of the effects of magnetic fields on fungi.

Modelling the spatiotemporal dynamics of senescent cells in wound healing, chronic wounds, and fibrosis

Cellular senescence is known to drive age-related pathology through the senescence-associated secretory phenotype (SASP). However, it also plays important physiological roles such as cancer suppression, embryogenesis and wound healing. Wound healing is a tightly regulated process which when disrupted results in conditions such as fibrosis and chronic wounds. Senescent cells appear during the proliferation phase of the healing process where the SASP is involved in maintaining tissue homeostasis after damage. Interestingly, SASP composition and functionality was recently found to be temporally regulated, with distinct SASP profiles involved: a fibrogenic, followed by a fibrolytic SASP, which could have important implications for the role of senescent cells in wound healing. Given the number of factors at play a full understanding requires addressing the multiple levels of complexity, pertaining to the various cell behaviours, individually followed by investigating the interactions and influence each of these elements have on each other and the system as a whole. Here, a systems biology approach was adopted whereby a multi-scale model of wound healing that includes the dynamics of senescent cell behaviour and corresponding SASP composition within the wound microenvironment was developed. The model was built using the software CompuCell3D, which is based on a Cellular Potts modelling framework. We used an existing body of data on healthy wound healing to calibrate the model and validation was done on known disease conditions. The model provides understanding of the spatiotemporal dynamics of different senescent cell phenotypes and the roles they play within the wound healing process. The model also shows how an overall disruption of tissue-level coordination due to age-related changes results in different disease states including fibrosis and chronic wounds. Further specific data to increase model confidence could be used to explore senolytic treatments in wound disorders.

In silico model development and optimization of in vitro lung cell population growth

Tissue engineering predominantly relies on trial and error in vitro and ex vivo experiments to develop protocols and bioreactors to generate functional tissues. As an alternative, in silico methods have the potential to significantly reduce the timelines and costs of experimental programs for tissue engineering. In this paper, we propose a methodology to formulate, select, calibrate, and test mathematical models to predict cell population growth as a function of the biochemical environment and to design optimal experimental protocols for model inference of in silico model parameters. We systematically combine methods from the experimental design, mathematical statistics, and optimization literature to develop unique and explainable mathematical models for cell population dynamics. The proposed methodology is applied to the development of this first published model for a population of the airway-relevant bronchio-alveolar epithelial (BEAS-2B) cell line as a function of the concentration of metabolic-related biochemical substrates. The resulting model is a system of ordinary differential equations that predict the temporal dynamics of BEAS-2B cell populations as a function of the initial seeded cell population and the glucose, oxygen, and lactate concentrations in the growth media, using seven parameters rigorously inferred from optimally designed in vitro experiments.

Effects of CRISPR/Cas9 targeting of the myo-inositol biosynthesis pathway on hyper-osmotic tolerance of tilapia cells

Myo-inositol is an important compatible osmolyte in vertebrates. This osmolyte is produced by the myo-inositol biosynthesis (MIB) pathway composed of myo-inositol phosphate synthase and inositol monophosphatase. These enzymes are among the highest upregulated proteins in tissues and cell cultures from teleost fish exposed to hyperosmotic conditions indicating high importance of this pathway for tolerating this type of stress. CRISPR/Cas9 gene editing of tilapia cells produced knockout lines of MIB enzymes and control genes. Metabolic activity decreased significantly for MIB KO lines in hyperosmotic media. Trends of faster growth of the MIB knockout lines in isosmotic media and faster decline of MIB knockout lines in hyperosmotic media were also observed. These results indicate a decline in metabolic fitness but only moderate effects on cell survival when tilapia cells with disrupted MIB genes are exposed to hyperosmolality. Therefore MIB genes are required for full osmotolerance of tilapia cells.

Dynamical model of antibiotic responses linking expression of resistance genes to metabolism explains emergence of heterogeneity during drug exposures

Antibiotic responses in bacteria are highly dynamic and heterogeneous, with sudden exposure of bacterial colonies to high drug doses resulting in the coexistence of recovered and arrested cells. The dynamics of the response is determined by regulatory circuits controlling the expression of resistance genes, which are in turn modulated by the drug's action on cell growth and metabolism. Despite advances in understanding gene regulation at the molecular level, we still lack a framework to describe how feedback mechanisms resulting from the interdependence between expression of resistance and cell metabolism can amplify naturally occurring noise and create heterogeneity at the population level. To understand how this interplay affects cell survival upon exposure, we constructed a mathematical model of the dynamics of antibiotic responses that links metabolism and regulation of gene expression, based on the tetracycline resistancetetoperon inE. coli. We use this model to interpret measurements of growth and expression of resistance in microfluidic experiments, both in single cells and in biofilms. We also implemented a stochastic model of the drug response, to show that exposure to high drug levels results in large variations of recovery times and heterogeneity at the population level. We show that stochasticity is important to determine how nutrient quality affects cell survival during exposure to high drug concentrations. A quantitative description of how microbes respond to antibiotics in dynamical environments is crucial to understand population-level behaviors such as biofilms and pathogenesis.

Fitness effects of a demography-dispersal trade-off in expandingSaccharomyces cerevisiaemats

Fungi expand in space and time to form complex multicellular communities. The mechanisms by which they do so can vary dramatically and determine the life-history and dispersal traits of expanding populations. These traits influence deterministic and stochastic components of evolution, resulting in complex eco-evolutionary dynamics during colony expansion. We perform experiments on budding yeast strains genetically engineered to display rough-surface and smooth-surface phenotypes in colony-like structures called 'mats'. Previously, it was shown that the rough-surface strain has a competitive advantage over the smooth-surface strain when grown on semi-solid media. We experimentally observe the emergence and expansion of segments with a distinct smooth-surface phenotype during rough-surface mat development. We propose a trade-off between dispersal and local carrying capacity to explain the relative fitness of these two phenotypes. Using a modified stepping-stone model, we demonstrate that this trade-off gives the high-dispersing, rough-surface phenotype a competitive advantage from standing variation, but that it inhibits this phenotype's ability to invade a resident smooth-surface population via mutation. However, the trade-off improves the ability of the smooth-surface phenotype to invade in rough-surface mats, replicating the frequent emergence of smooth-surface segments in experiments. Together, these computational and experimental findings advance our understanding of the complex eco-evolutionary dynamics of fungal mat expansion.

Modularized Design and Construction of Tunable Microbial Consortia with Flexible Topologies

Complex and fluid bacterial community compositions are critical to diversity, stability, and function. However, quantitative and mechanistic descriptions of the dynamics of such compositions are still lacking. Here, we develop a modularized design framework that allows for bottom-up construction and the study of synthetic bacterial consortia with different topologies. We showcase the microbial consortia design and building process by constructing amensalism and competition consortia using only genetic circuit modules to engineer different strains to form the community. Functions of modules and hosting strains are validated and quantified to calibrate dynamic parameters, which are then directly fed into a full mechanistic model to accurately predict consortia composition dynamics for both amensalism and competition without further fitting. More importantly, such quantitative understanding successfully identifies the experimental conditions to achieve coexistence composition dynamics. These results illustrate the process of both computationally and experimentally building up bacteria consortia complexity and hence achieve robust control of such fluid systems.

Aging effects on osteoclast progenitor dynamics affect variability in bone turnover via feedback regulation

Bone turnover markers (BTMs) are commonly used in osteoporosis treatment as indicators of cell activities of bone-resorbing osteoclasts and bone-forming osteoblasts. The wide variability in their values due to multiple factors, such as aging and diseases, makes it difficult for physicians to utilize them for clinical decision-making. The progenitors of osteoclasts and osteoblasts are indispensable for a comprehensive interpretation of the variability in BTM values because these upstream progenitors strongly regulate the downstream cell activities of bone turnover. However, understanding the complex interactions among the multiple populations of bone cells is challenging. In this study, we aimed to gain a fundamental understanding of the mechanism by which the progenitor dynamics affect the variability in bone turnover through in silico experiments by exploring the cell dynamics with aging effects on osteoporosis. Negative feedback control driven by the consumptive loss of progenitors prevents rapid bone loss due to excessive bone turnover, and through feedback regulation, aging effects on osteoclast differentiation and osteoclast progenitor proliferation cause variability in the osteoclast and osteoblast activity balance and its temporal transition. By expressing the variability in the bone turnover status, our model describes the individualities of patients based on their clinical backgrounds. Therefore, our model could play a powerful role in assisting tailored treatment and has the potential to resolve the various health problems associated with osteoporosis worldwide.

Competitive inhibition and mutualistic growth in co-infections: deciphering Staphylococcus aureus-Acinetobacter baumannii interaction dynamics

Staphylococcus aureus (Sa) and Acinetobacter baumannii (Ab) are frequently co-isolated from polymicrobial infections that are severe and refractory to therapy. Here, we apply a combination of wet-lab experiments and in silico modeling to unveil the intricate nature of the Ab/Sa interaction using both, representative laboratory strains and strains co-isolated from clinical samples. This comprehensive methodology allowed uncovering Sa's capability to exert a partial interference on Ab by the expression of phenol-soluble modulins. In addition, we observed a cross-feeding mechanism by which Sa supports the growth of Ab by providing acetoin as an alternative carbon source. This study is the first to dissect the Ab/Sa interaction dynamics wherein competitive and cooperative strategies can intertwine. Through our findings, we illuminate the ecological mechanisms supporting their coexistence in the context of polymicrobial infections. Our research not only enriches our understanding but also opens doors to potential therapeutic avenues in managing these challenging infections.

Modelling the distinct roles of epithelial and stromal androgen receptor in the regulation of prostate epithelial dynamics

The prostate is an androgen-responsive organ, but the complex cellular and molecular interactions that mediate these responses remain incompletely defined. Here, I synthesise the existing literature to derive a simple conceptual framework describing the androgen-dependent regulation of prostate epithelial dynamics. In this framework, epithelial androgen receptor (AR) cell-autonomously controls luminal cell height, whereas stromal AR regulates the synthesis of growth factors that promote luminal cell survival and proliferation. With the additional aid of a reanalysis of single-cell RNA-seq data, I also propose that insulin-like growth factor 1 (IGF1) functions as a key androgen-dependent growth factor coordinating stromal-to-epithelial paracrine communication. A novel mathematical model based on this framework was able to quantitatively fit experimental data describing prostate regression and regeneration. Model analysis demonstrates how the luminal cell population can maintain a stable equilibrium size via competition for and degradation of stroma-derived IGF1 and how this population size can be controlled by androgen levels, without a requirement for distinct luminal cell subsets. Moreover, model simulations were able to qualitatively recapitulate experimental observations in inflammatory and cancerous states, thereby providing insights into potential disease mechanisms. This simple model could therefore serve as a foundation for more comprehensive modelling of both the healthy and diseased prostate.

Unified tumor growth mechanisms from multimodel inference and dataset integration

Mechanistic models of biological processes can explain observed phenomena and predict responses to a perturbation. A mathematical model is typically constructed using expert knowledge and informal reasoning to generate a mechanistic explanation for a given observation. Although this approach works well for simple systems with abundant data and well-established principles, quantitative biology is often faced with a dearth of both data and knowledge about a process, thus making it challenging to identify and validate all possible mechanistic hypothesis underlying a system behavior. To overcome these limitations, we introduce a Bayesian multimodel inference (Bayes-MMI) methodology, which quantifies how mechanistic hypotheses can explain a given experimental datasets, and concurrently, how each dataset informs a given model hypothesis, thus enabling hypothesis space exploration in the context of available data. We demonstrate this approach to probe standing questions about heterogeneity, lineage plasticity, and cell-cell interactions in tumor growth mechanisms of small cell lung cancer (SCLC). We integrate three datasets that each formulated different explanations for tumor growth mechanisms in SCLC, apply Bayes-MMI and find that the data supports model predictions for tumor evolution promoted by high lineage plasticity, rather than through expanding rare stem-like populations. In addition, the models predict that in the presence of cells associated with the SCLC-N or SCLC-A2 subtypes, the transition from the SCLC-A subtype to the SCLC-Y subtype through an intermediate is decelerated. Together, these predictions provide a testable hypothesis for observed juxtaposed results in SCLC growth and a mechanistic interpretation for tumor treatment resistance.

Identification of actionable targets for breast cancer intervention using a diversity outbred mouse model

HER2-targeted therapy has improved breast cancer survival, but treatment resistance and disease prevention remain major challenges. Genes that enable HER2/Neu oncogenesis are the next intervention targets. A bioinformatics discovery platform of HER2/Neu-expressing Diversity Outbred (DO) F1 Mice was established to identify cancer-enabling genes. Quantitative Trait Loci (QTL) associated with onset ages and growth rates of spontaneous mammary tumors were sought. Twenty-six genes in 3 QTL contain sequence variations unique to the genetic backgrounds that are linked to aggressive tumors and 21 genes are associated with human breast cancer survival. Concurrent identification of TSC22D3, a transcription factor, and its target gene LILRB4, a myeloid cell checkpoint receptor, suggests an immune axis for regulation, or intervention, of disease. We also investigated TIEG1 gene that impedes tumor immunity but suppresses tumor growth. Although not an actionable target, TIEG1 study revealed genetic regulation of tumor progression, forming the basis of the genetics-based discovery platform.

Basic Methods of Cell Cycle Analysis

Cellular growth and the preparation of cells for division between two successive cell divisions is called the cell cycle. The cell cycle is divided into several phases; the length of these particular cell cycle phases is an important characteristic of cell life. The progression of cells through these phases is a highly orchestrated process governed by endogenous and exogenous factors. For the elucidation of the role of these factors, including pathological aspects, various methods have been developed. Among these methods, those focused on the analysis of the duration of distinct cell cycle phases play important role. The main aim of this review is to guide the readers through the basic methods of the determination of cell cycle phases and estimation of their length, with a focus on the effectiveness and reproducibility of the described methods.

Time-resolved microfluidics unravels individual cellular fates during double-strand break repair

BACKGROUND

Double-strand break repair (DSBR) is a highly regulated process involving dozens of proteins acting in a defined order to repair a DNA lesion that is fatal for any living cell. Model organisms such as Saccharomyces cerevisiae have been used to study the mechanisms underlying DSBR, including factors influencing its efficiency such as the presence of distinct combinations of microsatellites and endonucleases, mainly by bulk analysis of millions of cells undergoing repair of a broken chromosome. Here, we use a microfluidic device to demonstrate in yeast that DSBR may be studied at a single-cell level in a time-resolved manner, on a large number of independent lineages undergoing repair.

RESULTS

We used engineered S. cerevisiae cells in which GFP is expressed following the successful repair of a DSB induced by Cas9 or Cpf1 endonucleases, and different genetic backgrounds were screened to detect key events leading to the DSBR efficiency. Per condition, the progenies of 80-150 individual cells were analyzed over 24 h. The observed DSBR dynamics, which revealed heterogeneity of individual cell fates and their contributions to global repair efficacy, was confronted with a coupled differential equation model to obtain repair process rates. Good agreement was found between the mathematical model and experimental results at different scales, and quantitative comparisons of the different experimental conditions with image analysis of cell shape enabled the identification of three types of DSB repair events previously not recognized: high-efficacy error-free, low-efficacy error-free, and low-efficacy error-prone repair.

CONCLUSIONS

Our analysis paves the way to a significant advance in understanding the complex molecular mechanism of DSB repair, with potential implications beyond yeast cell biology. This multiscale and multidisciplinary approach more generally allows unique insights into the relation between in vivo microscopic processes within each cell and their impact on the population dynamics, which were inaccessible by previous approaches using molecular genetics tools alone.

Non-genetic resistance facilitates survival while hindering the evolution of drug resistance due to intraspecific competition

Rising rates of resistance to antimicrobial drugs threaten the effective treatment of infections across the globe. Drug resistance has been established to emerge from non-genetic mechanisms as well as from genetic mechanisms. However, it is still unclear how non-genetic resistance affects the evolution of genetic drug resistance. We develop deterministic and stochastic population models that incorporate resource competition to quantitatively investigate the transition from non-genetic to genetic resistance during the exposure to static and cidal drugs. We find that non-genetic resistance facilitates the survival of cell populations during drug treatment while hindering the development of genetic resistance due to competition between the non-genetically and genetically resistant subpopulations. Non-genetic resistance in the presence of subpopulation competition increases the fixation times of drug resistance mutations, while increasing the probability of mutation before population extinction during cidal drug treatment. Intense intraspecific competition during drug treatment leads to extinction of susceptible and non-genetically resistant subpopulations. Alternating between drug and no drug conditions results in oscillatory population dynamics, increased resistance mutation fixation timescales, and reduced population survival. These findings advance our fundamental understanding of the evolution of resistance and may guide novel treatment strategies for patients with drug-resistant infections.

Benchmarking imputation methods for network inference using a novel method of synthetic scRNA-seq data generation

Background

Single cell RNA-sequencing (scRNA-seq) has very rapidly become the new workhorse of modern biology providing an unprecedented global view on cellular diversity and heterogeneity. In particular, the structure of gene-gene expression correlation contains information on the underlying gene regulatory networks. However, interpretation of scRNA-seq data is challenging due to specific experimental error and biases that are unique to this kind of data including drop-out (or technical zeros).

Methods

To deal with this problem several methods for imputation of zeros for scRNA-seq have been developed. However, it is not clear how these processing steps affect inference of genetic networks from single cell data. Here, we introduce Biomodelling.jl, a tool for generation of synthetic scRNA-seq data using multiscale modelling of stochastic gene regulatory networks in growing and dividing cells.

Results

Our tool produces realistic transcription data with a known ground truth network topology that can be used to benchmark different approaches for gene regulatory network inference. Using this tool we investigate the impact of different imputation methods on the performance of several network inference algorithms.

Conclusions

Biomodelling.jl provides a versatile and useful tool for future development and benchmarking of network inference approaches using scRNA-seq data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04778-9

Population balance modelling captures host cell protein dynamics in CHO cell cultures

Monoclonal antibodies (mAbs) have been extensively studied for their wide therapeutic and research applications. Increases in mAb titre has been achieved mainly by cell culture media/feed improvement and cell line engineering to increase cell density and specific mAb productivity. However, this improvement has shifted the bottleneck to downstream purification steps. The higher accumulation of the main cell-derived impurities, host cell proteins (HCPs), in the supernatant can negatively affect product integrity and immunogenicity in addition to increasing the cost of capture and polishing steps. Mathematical modelling of bioprocess dynamics is a valuable tool to improve industrial production at fast rate and low cost. Herein, a single stage volume-based population balance model (PBM) has been built to capture Chinese hamster ovary (CHO) cell behaviour in fed-batch bioreactors. Using cell volume as the internal variable, the model captures the dynamics of mAb and HCP accumulation extracellularly under physiological and mild hypothermic culture conditions. Model-based analysis and orthogonal measurements of lactate dehydrogenase activity and double-stranded DNA concentration in the supernatant show that a significant proportion of HCPs found in the extracellular matrix is secreted by viable cells. The PBM then served as a platform for generating operating strategies that optimise antibody titre and increase cost-efficiency while minimising impurity levels.

Network modeling-based identification of the switching targets between pyroptosis and secondary pyroptosis

The newly identified cell death type, pyroptosis plays crucial roles in various diseases. Most recently, mounting evidence accumulates that pyroptotic signaling is highly correlated with coronavirus disease 2019 (COVID-19). Thus, understanding the induction of the pyroptotic signaling and dissecting the detail molecular control mechanisms are urgently needed. Based on recent experimental studies, a core regulatory model of the pyroptotic signaling is constructed to investigate the intricate crosstalk dynamics between the two cell death types, i.e., pyroptosis and secondary pyroptosis. The model well reproduces the experimental observations under different conditions. Sensitivity analysis determines that only the expression level of caspase-1 or GSDMD has the potential to individually change death modes. The decrease of caspase-1 or GSDMD level switches cell death from pyroptosis to secondary pyroptosis. Besides, eight biochemical reactions are identified that can efficiently switch death modes. While from the viewpoint of bifurcation analysis, the expression level of caspase-3 is further identified and twelve biochemical reactions are obtained. The coexistence of pyroptosis and secondary pyroptosis is predicted to be observed not only within the bistable range, but also within proper monostable range, presenting two potential different control mechanisms. Combined with the landscape theory, we further explore the stochastic dynamic and global stability of the pyroptotic system, accurately quantifying how each component mediates the individual occurrence probability of pyroptosis and secondary pyroptosis. Overall, this study sheds new light on the intricate crosstalk of the pyroptotic signaling and uncovers the regulatory mechanisms of various stable state transitions, providing potential clues to guide the development for prevention and treatment of pyroptosis-related diseases.

Theoretical derivation and clinical dose-response quantification of a unified multi-activation (UMA) model of cell survival from a logistic equation

OBJECTIVE

To theoretically derive a unified multiactivation (UMA) model of cell survival after ionising radiation that can accurately assess doses and responses in radiotherapy and X-ray imaging.

METHODS

A unified formula with only two parameters in fitting of a cell survival curve (CSC) is first derived from an assumption that radiation-activated cell death pathways compose the first- and second-order reaction kinetics. A logit linear regression of CSC data is used for precise determination of the two model parameters. Intrinsic radiosensitivity, biologically effective dose (BED), equivalent dose to the traditional 2 Gy fractions (EQD2), tumour control probability, normal-tissue complication probability, BED₅₀ and steepness (Γ50) at 50% of tumour control probability (or normal-tissue complication probability) are analytical functions of the model and treatment (or imaging) parameters.

RESULTS

The UMA model has almost perfectly fit typical CSCs over the entire dose range with R²≥0.99. Estimated quantities for stereotactic body radiotherapy of early stage lung cancer and the skin reactions from X-ray imaging agree with clinical results.

CONCLUSION

The proposed UMA model has theoretically resolved the catastrophes of the zero slope at zero dose for multiple target model and the bending curve at high dose for the linear quadratic model. More importantly, it analytically predicts dose-responses to various dose-fraction schemes in radiotherapy and to low dose X-ray imaging based on these preclinical CSCs.

ADVANCES IN KNOWLEDGE

The discovery of a unified formula of CSC over the entire dose range may reveal a common mechanism of the first- and second-order reaction kinetics among multiple CD pathways activated by ionising radiation at various dose levels.

Quantifying the optimal strategy of population control of quorum sensing network in Escherichia coli

Biological functions of bacteria can be regulated by monitoring their own population density induced by the quorum sensing system. However, quantitative insight into the system’s dynamics and regulatory mechanism remain challenging. Here, we construct a comprehensive mathematical model of the synthetic quorum sensing circuit that controls population density in Escherichia coli. Simulations agree well with experimental results obtained under different ribosome-binding site (RBS) efficiencies. We present a quantitative description of the component dynamics and show how the components respond to isopropyl-β-D-1-thiogalactopyranoside (IPTG) induction. The optimal IPTG-induction range for efficiently controlling population density is quantified. The controllable area of population density by acyl-homoserine lactone (AHL) permeability is quantified as well, indicating that high AHL permeability should be treated with a high dose of IPTG, while low AHL permeability should be induced with low dose for efficiently controlling. Unexpectedly, an oscillatory behavior of the growth curve is observed with proper RBS-binding strengths and the oscillation is greatly restricted by the bacterial death induced by toxic metabolic by-products. Moreover, we identify that the mechanism underlying the emergence of oscillation is determined by the negative feedback loop structure within the signaling. Bifurcation analysis and landscape theory are further employed to study the stochastic dynamic and global stability of the system, revealing two faces of toxic metabolic by-products in controlling oscillatory behavior. Overall, our study presents a quantitative basis for understanding and new insights into the control mechanism of quorum sensing system, providing possible clues to guide the development of more rational control strategy.

Cell death modes are specified by the crosstalk dynamics within pyroptotic and apoptotic signaling

The crosstalk between pyroptosis and apoptosis pathways plays crucial roles in homeostasis, cancer, and other pathologies. However, its molecular regulatory mechanisms for cell death decision-making remain to be elucidated. Based on the recent experimental studies, we developed a core regulatory network model of the crosstalk between pyroptosis and apoptosis pathways. Sensitivity analysis and bifurcation analysis were performed to assess the death mode switching of the network. Both the approaches determined that only the level of caspase-1 or gasdermin D (GSDMD) has the potential to individually change death modes. The decrease of caspase-1 or GSDMD switches cell death from pyroptosis to apoptosis. Seven biochemical reactions among the 21 reactions in total that are essential for determining cell death modes are identified by using sensitivity analysis. While with bifurcation analysis of state transitions, nine reactions are suggested to be able to efficiently switch death modes. Monostability, bistability, and tristability are observed under different conditions. We found that only the reaction that caspase-1 activation induced by stimuli can trigger tristability. Six and two of the nine reactions are identified to be able to induce bistability and monostability, respectively. Moreover, the concurrence of pyroptosis and apoptosis is observed not only within proper bistable ranges, but also within tristable ranges, implying two potentially distinct regulatory mechanisms. Taken together, this work sheds new light on the crosstalk between pyroptosis and apoptosis and uncovers the regulatory mechanisms of various stable state transitions, which play important roles for the development of potential control strategies for disease prevention and treatment.

Does transcriptional heterogeneity facilitate the development of genetic drug resistance?

Non-genetic forms of antimicrobial (drug) resistance can result from cell-to-cell variability that is not encoded in the genetic material. Data from recent studies also suggest that non-genetic mechanisms can facilitate the development of genetic drug resistance. We speculate on how the interplay between non-genetic and genetic mechanisms may affect microbial adaptation and evolution during drug treatment. We argue that cellular heterogeneity arising from fluctuations in gene expression, epigenetic modifications, as well as genetic changes contribute to drug resistance at different timescales, and that the interplay between these mechanisms enhance pathogen resistance. Accordingly, developing a better understanding of the role of non-genetic mechanisms in drug resistance and how they interact with genetic mechanisms will enhance our ability to combat antimicrobial resistance. Also see the video abstract here: https://youtu.be/aefGpdh-bgU.

Seaweed cellulose scaffolds derived from green macroalgae for tissue engineering

Extracellular matrix (ECM) provides structural support for cell growth, attachments and proliferation, which greatly impact cell fate. Marine macroalgae species Ulva sp. and Cladophora sp. were selected for their structural variations, porous and fibrous respectively, and evaluated as alternative ECM candidates. Decellularization-recellularization approach was used to fabricate seaweed cellulose-based scaffolds for in-vitro mammalian cell growth. Both scaffolds were confirmed nontoxic to fibroblasts, indicated by high viability for up to 40 days in culture. Each seaweed cellulose structure demonstrated distinct impact on cell behavior and proliferation rates. The Cladophora sp. scaffold promoted elongated cells spreading along its fibers' axis, and a gradual linear cell growth, while the Ulva sp. porous surface, facilitated rapid cell growth in all directions, reaching saturation at week 3. As such, seaweed-cellulose is an environmentally, biocompatible novel biomaterial, with structural variations that hold a great potential for diverse biomedical applications, while promoting aquaculture and ecological agenda.

Coordination of gene expression noise with cell size: analytical results for agent-based models of growing cell populations

The chemical master equation and the Gillespie algorithm are widely used to model the reaction kinetics inside living cells. It is thereby assumed that cell growth and division can be modelled through effective dilution reactions and extrinsic noise sources. We here re-examine these paradigms through developing an analytical agent-based framework of growing and dividing cells accompanied by an exact simulation algorithm, which allows us to quantify the dynamics of virtually any intracellular reaction network affected by stochastic cell size control and division noise. We find that the solution of the chemical master equation-including static extrinsic noise-exactly agrees with the agent-based formulation when the network under study exhibits stochastic concentration homeostasis, a novel condition that generalizes concentration homeostasis in deterministic systems to higher order moments and distributions. We illustrate stochastic concentration homeostasis for a range of common gene expression networks. When this condition is not met, we demonstrate by extending the linear noise approximation to agent-based models that the dependence of gene expression noise on cell size can qualitatively deviate from the chemical master equation. Surprisingly, the total noise of the agent-based approach can still be well approximated by extrinsic noise models.

Do we really understand how drug eluted from stents modulates arterial healing?

The advent of drug-eluting stents (DES) has revolutionised the treatment of coronary artery disease. These devices, coated with anti-proliferative drugs, are deployed into stenosed or occluded vessels, compressing the plaque to restore natural blood flow, whilst simultaneously combating the evolution of restenotic tissue. Since the development of the first stent, extensive research has investigated how further advancements in stent technology can improve patient outcome. Mathematical and computational modelling has featured heavily, with models focussing on structural mechanics, computational fluid dynamics, drug elution kinetics and subsequent binding within the arterial wall; often considered separately. Smooth Muscle Cell (SMC) proliferation and neointimal growth are key features of the healing process following stent deployment. However, models which depict the action of drug on these processes are lacking. In this article, we start by reviewing current models of cell growth, which predominantly emanate from cancer research, and available published data on SMC proliferation, before presenting a series of mathematical models of varying complexity to detail the action of drug on SMC growth in vitro. Our results highlight that, at least for Sodium Salicylate and Paclitaxel, the current state-of-the-art nonlinear saturable binding model is incapable of capturing the proliferative response of SMCs across a range of drug doses and exposure times. Our findings potentially have important implications on the interpretation of current computational models and their future use to optimise and control drug release from DES and drug-coated balloons.

Study of Combinatorial Drug Synergy of Novel Acridone Derivatives With Temozolomide Using in-silico and in-vitro Methods in the Treatment of Drug-Resistant Glioma

Drug resistance is one of the critical challenges faced in the treatment of Glioma. There are only limited drugs available in the treatment of Glioma and among them Temozolomide (TMZ) has shown some effectiveness in treating Glioma patients, however, the rate of recovery remains poor due to the inability of this drug to act on the drug resistant tumor sub-populations. Hence, in this study three novel Acridone derivative drugs AC2, AC7, and AC26 have been proposed. These molecules when combined with TMZ show major tumor cytotoxicity that is effective in suppressing growth of cancer cells in both drug sensitive and resistant sub-populations of a tumor. In this study a novel mathematical model has been developed to explore the various drug combinations that may be useful for the treatment of resistant Glioma and show that the combinations of TMZ and Acridone derivatives have a synergistic effect. Also, acute toxicity studies of all three acridone derivatives were carried out for 14 days and were found safe for oral administration of 400 mg/kg body weight on albino Wistar rats. Molecular Docking studies of acridone derivatives with P-glycoprotein (P-gp), multiple resistant protein (MRP), and O6-methylguanine-DNA methyltransferase (MGMT) revealed different binding affinities to the transporters contributing to drug resistance. It is observed that while the Acridone derivatives bind with these drug resistance causing proteins, the TMZ can produce its cytotoxicity at a much lower concentration leading to the synergistic effect. The in silico analysis corroborate well with our experimental findings using TMZ resistant (T-98) and drug sensitive (U-87) Glioma cell lines and we propose three novel drug combinations (TMZ with AC2, AC7, and AC26) and dosages that show high synergy, high selectivity and low collateral toxicity for the use in the treatment of drug resistant Glioma, which could be future drugs in the treatment of Glioblastoma.

The influence of methylation status on a stochastic model of MGMT dynamics in glioblastoma: Phenotypic selection can occur with and without a downshift in promoter methylation status

Glioblastoma originates in the brain and is one of the most aggressive cancer types. Glioblastoma represents 15% of all brain tumours, with a median survival of 15 months. Although the current standard of care for such a tumour (the Stupp protocol) has shown positive results for the prognosis of patients, O-6-methylguanine-DNA methyltransferase (MGMT) driven drug resistance has been an issue of increasing concern and hence requires innovative approaches. In addition to the well established drug resistance factors such as tumour location and blood brain barriers, it is also important to understand how the genetic and epigenetic dynamics of the glioblastoma cells can play a role. One important aspect of this is the study of methylation status of MGMT following administration of temozolomide. In this paper, we extend our previously published model that simulated MGMT expression in glioblastoma cells to incorporate the promoter methylation status of MGMT. This methylation status has clinical significance and is used as a marker for patient outcomes. Using this model, we investigate the causative relationship between temozolomide treatment and the methylation status of the MGMT promoter in a population of cells. In addition by constraining the model to relevant biological data using Approximate Bayesian Computation, we were able to identify parameter regimes that yield different possible modes of resistances, namely, phenotypic selection of MGMT, a downshift in the methylation status of the MGMT promoter or both simultaneously. We analysed each of the parameter sets associated with the different modes of resistance, presenting representative solutions as well as discovering some similarities between them as well as unique requirements for each of them. Finally, we used them to devise optimal strategies for inhibiting MGMT expression with the aim of minimising live glioblastoma cell numbers.

MEDALT: single-cell copy number lineage tracing enabling gene discovery

We present a Minimal Event Distance Aneuploidy Lineage Tree (MEDALT) algorithm that infers the evolution history of a cell population based on single-cell copy number (SCCN) profiles, and a statistical routine named lineage speciation analysis (LSA), whichty facilitates discovery of fitness-associated alterations and genes from SCCN lineage trees. MEDALT appears more accurate than phylogenetics approaches in reconstructing copy number lineage. From data from 20 triple-negative breast cancer patients, our approaches effectively prioritize genes that are essential for breast cancer cell fitness and predict patient survival, including those implicating convergent evolution.The source code of our study is available at https://github.com/KChen-lab/MEDALT .

Inferring time-dependent population growth rates in cell cultures undergoing adaptation

Background

The population growth rate is an important characteristic of any cell culture. During sustained experiments, the growth rate may vary due to competition or adaptation. For instance, in presence of a toxin or a drug, an increasing growth rate indicates that the cells adapt and become resistant. Consequently, time-dependent growth rates are fundamental to follow on the adaptation of cells to a changing evolutionary landscape. However, as there are no tools to calculate the time-dependent growth rate directly by cell counting, it is common to use only end point measurements of growth rather than tracking the growth rate continuously.

Results

We present a computer program for inferring the growth rate over time in suspension cells using nothing but cell counts, which can be measured non-destructively. The program was tested on simulated and experimental data. Changes were observed in the initial and absolute growth rates, betraying resistance and adaptation.

Conclusions

For experiments where adaptation is expected to occur over a longer time, our method provides a means of tracking growth rates using data that is normally collected anyhow for monitoring purposes. The program and its documentation are freely available at https://github.com/Sandalmoth/ratrack under the permissive zlib license.

Influence of Culture Conditions and Medium Compositions on the Production of Bacteriocin-Like Inhibitory Substances by Lactococcus lactis Gh1

Antibacterial peptides or bacteriocins produced by many strains of lactic acid bacteria have been used as food preservatives for many years without any known adverse effects. Bacteriocin titres can be modified by altering the physiological and nutritional factors of the producing bacterium to improve the production in terms of yield and productivity. The effects of culture conditions (initial pH, inoculum age and inoculum size) and medium compositions (organic and inorganic nitrogen sources; carbon sources) were assessed for the production of bacteriocin-like inhibitory substances (BLIS) by Lactococcus lactis Gh1 in shake flask cultures. An inoculum of the mid-exponential phase culture at 1% (v/v) was the optimal age and size, while initial pH of culture media at alkaline and acidic state did not show a significant impact on BLIS secretion. Organic nitrogen sources were more favourable for BLIS production compared to inorganic sources. Production of BLIS by L. lactis Gh1 in soytone was 1.28-times higher as compared to that of organic nitrogen sources ((NH4)2SO4). The highest cell concentration (XmX = 0.69 ± 0.026 g·L-1) and specific growth rate (μmax = 0.14 h-1) were also observed in cultivation using soytone. By replacing carbon sources with fructose, BLIS production was increased up to 34.94% compared to BHI medium, which gave the biomass cell concentration and specific growth rate of 0.66 ± 0.002 g·L-1 and 0.11 h-1, respectively. It can be concluded that the fermentation factors have pronounced influences on the growth of L. lactis Gh1 and BLIS production. Results from this study could be used for subsequent application in process design and optimisation for improving BLIS production by L. lactis Gh1 at larger scale.

Advancing Antimicrobial Resistance Research Through Quantitative Modeling and Synthetic Biology

Antimicrobial resistance (AMR) is an emerging global health crisis that is undermining advances in modern medicine and, if unmitigated, threatens to kill 10 million people per year worldwide by 2050. Research over the last decade has demonstrated that the differences between genetically identical cells in the same environment can lead to drug resistance. Fluctuations in gene expression, modulated by gene regulatory networks, can lead to non-genetic heterogeneity that results in the fractional killing of microbial populations causing drug therapies to fail; this non-genetic drug resistance can enhance the probability of acquiring genetic drug resistance mutations. Mathematical models of gene networks can elucidate general principles underlying drug resistance, predict the evolution of resistance, and guide drug resistance experiments in the laboratory. Cells genetically engineered to carry synthetic gene networks regulating drug resistance genes allow for controlled, quantitative experiments on the role of non-genetic heterogeneity in the development of drug resistance. In this perspective article, we emphasize the contributions that mathematical, computational, and synthetic gene network models play in advancing our understanding of AMR to discover effective therapies against drug-resistant infections.

cuProCell: GPU-Accelerated Analysis of Cell Proliferation With Flow Cytometry Data

The investigation of cell proliferation can provide useful insights for the comprehension of cancer progression, resistance to chemotherapy and relapse. To this aim, computational methods and experimental measurements based on in vivo label-retaining assays can be coupled to explore the dynamic behavior of tumoral cells. ProCell is a software that exploits flow cytometry data to model and simulate the kinetics of fluorescence loss that is due to stochastic events of cell division. Since the rate of cell division is not known, ProCell embeds a calibration process that might require thousands of stochastic simulations to properly infer the parameterization of cell proliferation models. To mitigate the high computational costs, in this paper we introduce a parallel implementation of ProCell's simulation algorithm, named cuProCell, which leverages Graphics Processing Units (GPUs). Dynamic Parallelism was used to efficiently manage the cell duplication events, in a radically different way with respect to common computing architectures. We present the advantages of cuProCell for the analysis of different models of cell proliferation in Acute Myeloid Leukemia (AML), using data collected from the spleen of human xenografts in mice. We show that, by exploiting GPUs, our method is able to not only automatically infer the models' parameterization, but it is also 237× faster than the sequential implementation. This study highlights the presence of a relevant percentage of quiescent and potentially chemoresistant cells in AML in vivo, and suggests that maintaining a dynamic equilibrium among the different proliferating cell populations might play an important role in disease progression.

Phenotypic selection through cell death: stochastic modelling of O-6-methylguanine-DNA methyltransferase dynamics

Glioblastoma (GBM) is the most aggressive malignant primary brain tumour with a median overall survival of 15 months. To treat GBM, patients currently undergo a surgical resection followed by exposure to radiotherapy and concurrent and adjuvant temozolomide (TMZ) chemotherapy. However, this protocol often leads to treatment failure, with drug resistance being the main reason behind this. To date, many studies highlight the role of O-6-methylguanine-DNA methyltransferase (MGMT) in conferring drug resistance. The mechanism through which MGMT confers resistance is not well studied-particularly in terms of computational models. With only a few reasonable biological assumptions, we were able to show that even a minimal model of MGMT expression could robustly explain TMZ-mediated drug resistance. In particular, we showed that for a wide range of parameter values constrained by novel cell growth and viability assays, a model accounting for only stochastic gene expression of MGMT coupled with cell growth, division, partitioning and death was able to exhibit phenotypic selection of GBM cells expressing MGMT in response to TMZ. Furthermore, we found this selection allowed the cells to pass their acquired phenotypic resistance onto daughter cells in a stable manner (as long as TMZ is provided). This suggests that stochastic gene expression alone is enough to explain the development of chemotherapeutic resistance.

Modeling heterogeneous tumor growth dynamics and cell-cell interactions at single-cell and cell-population resolution

Cancer is a complex, dynamic disease that despite recent advances remains mostly incurable. Inter- and intratumoral heterogeneity are generally considered major drivers of therapy resistance, metastasis, and treatment failure. Recent advances in high-throughput experimentation have produced a wealth of data on tumor heterogeneity and researchers are increasingly turning to mathematical modeling to aid in the interpretation of these complex datasets. In this mini-review, we discuss three important classes of approaches for modeling cellular dynamics within heterogeneous tumors: agent-based models, population dynamics, and multiscale models. An important new focus, for which we provide an example, is the role of intratumoral cell-cell interactions.