Cellular adaptation to external signals is essential for biological functions, and it is an important field of interest in systems biology. This study examines the impact of cooperativity on the adaptation response of the Incoherent Feedforward Loop (IFFL) network motif to various signal profiles. Through comprehensive simulations, we studied how the IFFL motif responds to constant and pulse-type signals under varying levels of cooperativity. The results of our study demonstrate that positive cooperativity generally enhances the system's ability to adapt to different signal profiles. Nevertheless, given specific signal profiles, higher levels of cooperativity may decrease the system's adaptability. On the other hand, the adaptive response breaks down for negative cooperativity. For constant signals, increased positive cooperativity leads to a response with higher amplitude, and it accelerates the response time but delays the return time required to settle back down to the pre-stimulus state. Upon signal cessation, high positive cooperativity not only slows the system's response and return times but, in some cases, can lead to a complete temporary halt in response. For the pulse-like signal, cooperativity increases the maximum amplitude of the oscillatory response. These insights highlight the delicate balance between cooperativity and signal profile in cellular adaptation mechanisms involving the IFFL network motif.
Single cell transcriptomics has recently seen a surge in popularity, leading to the need for data analysis pipelines that are reproducible, modular, and interoperable across different systems and institutions.To meet this demand, we introduce scAN1.0, a processing pipeline for analyzing 10X single cell RNA sequencing data. scAN1.0 is built using the Nextflow DSL2 and can be run on most computational systems. The modular design of Nextflow pipelines enables easy integration and evaluation of different blocks for specific analysis steps.We demonstrate the usefulness of scAN1.0 by showing its ability to examine the impact of the mapping step during the analysis of two datasets: (i) a 10X scRNAseq of a human pituitary gonadotroph tumor dataset and (ii) a murine 10X scRNAseq acquired on CD8 T cells during an immune response.
Darwin's and Wallace's 1859 explanation that novel speciation resulted from natural variants that had been subjected to selection was refined over the next 150 years as genetic inheritance and the importance of mutation-induced change were discovered, the quantitative theory of evolutionary population genetics was produced, the speed of genetic change in small populations became apparent and the ramifications of the DNA revolution became clear. This paper first discusses the modern view of speciation in its historical context. It then uses systems-biology approaches to consider the many complex processes that underpin the production of a new species; these extend in scale from genes to populations with the processes of variation, selection and speciation being affected by factors that range from mutation to climate change. Here, events at a particular scale level (e.g. protein network activity) are activated by the output of the level immediately below (i.e. gene expression) and generate a new output that activates the layer above (e.g. embryological development), with this change often being modulated by feedback from higher and lower levels. The analysis shows that activity at each level in the evolution of a new species is marked by stochastic activity, with mutation of course being the key step for variation. The paper examines events at each of these scale levels and particularly considers how the pathway by which mutation leads to phenotypic variants and the wide range of factors that drive selection can be investigated computationally. It concludes that, such is the complexity of speciation, most steps in the process are currently difficult to model and that predictions about future speciation will, apart from a few special cases, be hard to make. The corollary is that opportunities for novel variants to form are maximised.
CABs (Concave Actin Bundles) are oriented against the scaffold transversally in a manner different from traditional longitudinal F-actin bundles. CABs are present in a specific area, and do not exist in random areas. Biologically, CABs are developed to attach cells to fibers firmly so that CABs are found near cells. Based on this knowledge, we closely examined 3D confocal microcopy images containing fiber scaffolds, actin, and cells. Then, we assumed that the areas containing high values of compactness of fiber, compactness of actin, and density of cells would have many numbers of CABs.In this research, we wanted to prove this assumption. We first incorporated a two-point correlation function to define a measure of compactness. Then, we used the Bayes' theorem to prove the above assumption. As the assumption, our results verified that CABs exist in an area of high compactness of a fiber network, high compactness of actin distribution, and high density of cells. Thus, we concluded that CABs are developed to attach cells to a fibrillar scaffold firmly. This finding may be further verified mathematically in future studies.
Differentiation can be modeled at the single cell level as a stochastic process resulting from the dynamical functioning of an underlying Gene Regulatory Network (GRN), driving stem or progenitor cells to one or many differentiated cell types. Metastability seems inherent to differentiation process as a consequence of the limited number of cell types. Moreover, mRNA is known to be generally produced by bursts, which can give rise to highly variable non-Gaussian behavior, making the estimation of a GRN from transcriptional profiles challenging. In this article, we present CARDAMOM (Cell type Analysis from scRna-seq Data achieved from a Mixture MOdel), a new algorithm for inferring a GRN from timestamped scRNA-seq data, which crucially exploits these notions of metastability and transcriptional bursting. We show that such inference can be seen as the successive resolution of as many regression problem as timepoints, after a preliminary clustering of the whole set of cells with regards to their associated bursts frequency. We demonstrate the ability of CARDAMOM to infer a reliable GRN from in silico expression datasets, with good computational speed. To the best of our knowledge, this is the first description of a method which uses the concept of metastability for performing GRN inference.
Many biological ecosystems exhibit chaotic behavior, demonstrated either analytically using parameter choices in an associated dynamical systems model or empirically through analysis of experimental data. In this paper, we use existing software tools (COPASI, R) to explore dynamical systems and uncover regions with positive Lyapunov exponents where thus chaos exists. We evaluate the ability of the software’s optimization algorithms to find these positive values with several dynamical systems used to model biological populations. The algorithms have been able to identify parameter sets which lead to positive Lyapunov exponents, even when those exponents lie in regions with small support. For one of the examined systems, we observed that positive Lyapunov exponents were not uncovered when executing a search over the parameter space with small spacings between values of the independent variables.
Large-scale visualization and analysis of HPIs involved in microbial CVDs can provide crucial insights into the mechanisms of pathogenicity. The comparison of CVD associated HPIs with the entire set of HPIs can identify the pathways specific to CVDs. Therefore, topological properties of HPI networks in CVDs and all pathogens was studied using Cytoscape3.5.1. Ontology and pathway analysis were done using KOBAS 3.0. HPIs of Papilloma, Herpes, Influenza A virus as well as Yersinia pestis and Bacillus anthracis among bacteria were predominant in the whole (wHPI) and the CVD specific (cHPI) network. The central viral and secretory bacterial proteins were predicted virulent. The central viral proteins had higher number of interactions with host proteins in comparison with bacteria. Major fraction of central and essential host proteins interacts with central viral proteins. Alpha-synuclein, Ubiquitin ribosomal proteins, TATA-box-binding protein, and Polyubiquitin-C &B proteins were the top interacting proteins specific to CVDs. Signaling by NGF, Fc epsilon receptor, EGFR and ubiquitin mediated proteolysis were among the top enriched CVD specific pathways. DEXDc and HELICc were enriched host mimicry domains that may help in hijacking of cellular machinery by pathogens. This study provides a system level understanding of cardiac damage in microbe induced CVDs.
Vascular endothelial growth factor (VEGF) has been known as a key mediator of angiogenesis in cancer. Bevacizumab is anti-VEGF monoclonal antibody that has been approved by the FDA as a first-line treatment in many types of cancer. In this paper, we extend a previously validated multiscale tumor model to comprehensively include the multiple roles of VEGF during the course of angiogenesis and its binding mechanism with bevacizumab. We use the model to simulate tumor system response under various bevacizumab concentrations, both in stand-alone treatment and in combination with chemotherapy. Our simulation indicates that periodic administration of bevacizumab with lower concentration can achieve greater efficacy than a single treatment with higher concentration. The simulation of the combined therapy also shows that the continuous administration of bevacizumab during the maintenance phase can lead to antitumor activity which further suppresses its growth. Agreement with experimental results indicates the potential of the model in predicting the efficacy of anti-VEGF therapies and could therefore contribute to developing prospective clinical trials.
Yeasts exist in communities that expand over space and time to form complex structures and patterns. We developed a lattice-based framework to perform spatial-temporal Monte Carlo simulations of budding yeast colonies exposed to different nutrient and magnetic field conditions. The budding patterns of haploid and diploid yeast cells were incorporated into the framework, as well as the filamentous growth that occurs in yeast colonies under nutrient limiting conditions. Simulation of the framework predicted that magnetic fields decrease colony growth rate, solidity, and roundness. Magnetic field simulations further predicted that colony elongation and boundary fluctuations increase in a nutrient- and ploidy-dependent manner. These in-silico predictions are an important step towards understanding the effects of the physico-chemical environment on microbial colonies and for informing bioelectromagnetic experiments on yeast colony biofilms and fungal pathogens.
The role of the immune system in tumor development increasingly includes the idea of cancer immunoediting. It comprises three phases: elimination, equilibrium, and escape. In the first phase, elimination, transformed cells are recognized and destroyed by immune system. The rare tumor cells that are not destroyed in this phase may then enter the equilibrium phase, where their growth is prevented by immunity mechanisms. The escape phase represents the final phase of this process, where cancer cells begin to grow unconstrained by the immune system. In this study, we describe and analyze an evolutionary game theoretical model of proliferating, quiescent, and immune cells interactions for the first time. The proposed model is evaluated with constant and dynamic approaches. Population dynamics and interactions between the immune system and cancer cells are investigated. Stability of equilibria or critical points are analyzed by applying algebraic analysis. This model allows us to understand the process of cancer development and might help us design better treatment strategies to account for immunoediting.
To develop vaccines it is mandatory yet challenging to account for inter-individual variability during immune responses. Even in laboratory mice, T cell responses of single individuals exhibit a high heterogeneity that may come from genetic backgrounds, intra-specific processes (e.g. antigen-processing and presentation) and immunization protocols. To account for inter-individual variability in CD8 T cell responses in mice, we propose a dynamical model coupled to a statistical, nonlinear mixed effects model. Average and individual dynamics during a CD8 T cell response are characterized in different immunization contexts (vaccinia virus and tumor). On one hand, we identify biological processes that generate inter-individual variability (activation rate of naive cells, the mortality rate of effector cells, and dynamics of the immunogen). On the other hand, introducing categorical covariates to analyze two different immunization regimens, we highlight the steps of the response impacted by immunogens (priming, differentiation of naive cells, expansion of effector cells and generation of memory cells). The robustness of the model is assessed by confrontation to new experimental data. Our approach allows to investigate immune responses in various immunization contexts, when measurements are scarce or missing, and contributes to a better understanding of inter-individual variability in CD8 T cell immune responses.
Micro-Tissue Engineered Neural Networks (Micro-TENNs) are living three-dimensional constructs designed to replicate the neuroanatomy of white matter pathways in the brain, and are being developed as implantable microtissue for axon tract reconstruction or as anatomically-relevant in vitro experimental platforms. Micro-TENNs are composed of discrete neuronal aggregates connected by bundles of long-projecting axonal tracts within miniature tubular hydrogels. In order to help design and optimize micro-TENN performance, we have created a new computational model including geometric and functional properties. The model is built upon the three-dimensional diffusion equation and incorporates large-scale uni- and bi-directional growth that simulates realistic neuron morphologies. The model captures unique features of 3D axonal tract development that are not apparent in planar outgrowth, and may be insightful for how white matter pathways form during brain development. The processes of axonal outgrowth, branching, turning and aggregation/bundling from each neuron are described through functions built on concentration equations and growth time distributed across the growth segments. Once developed we conducted multiple parametric studies to explore the applicability of the method and conducted preliminary validation via comparisons to experimentally grown micro-TENNs for a range of growth conditions. Using this framework, this model can be applied to study micro-TENN growth processes and functional characteristics using spiking network or compartmental network modeling. This model may be applied to improve our understanding of axonal tract development and functionality, as well as to optimize the fabrication of implantable tissue engineered brain pathways for nervous system reconstruction and/or modulation.
A dynamical model of the pathophysiological behaviors of IL18 and IL10 cytokines with their receptors is tested against data for the case of early sepsis. The proposed approach considers the surroundings (organs and bone marrow) and the different subsystems (cells and cyctokines). The interactions between blood cells, cytokines and the surroundings are described via mass balances. Cytokines are adsorbed onto associated receptors at the cell surface. The adsorption is described by the Langmuir model and gives rise to the production of more cytokines and associated receptors inside the cell. The quantities of pro and anti-inflammatory cytokines present in the body are combined to give global information via an inflammation level function which describes the patient's state. Data for parameter estimation comes from the Sepsis 48 H database. Comparisons between patient data and simulations are presented and are in good agreement. For the IL18/IL10 cytokine pair, 5 key parameters have been found. They are linked to pro-inflammatory IL18 cytokine and show that the early sepsis is driven by components of inflammatory character.
Moonlighting refers to a protein with at least two unrelated, mechanistically different functions. As a concept, moonlighting describes a large and diverse group of proteins which have been discovered in a multitude of organisms. As of today, a systematized view on these proteins is missing. Here, we propose a classification of moonlighting proteins by two classifiers. We use the function of the protein as a first classifier: activating - activating (Type I), activating - inhibiting (Type II), inhibiting - activating (Type III) and inhibiting - inhibiting (Type IV). To further specify the type of moonlighting protein, we used a second classifier based on the character of the factor that switches the function of the protein: external factor affecting the protein (Type A), change in the first pathway (Type B), change in the second pathway (Type C), equal competition between both pathways (Type D). Using a small two-pathway model we simulated these types of moonlighting proteins to elucidate possible behaviors of the types of moonlighting proteins. We find that, using the results of our simulations, we can classify the behavior of the moonlighting types into Blinker, Splitter andSwitch.
The in vivo erythropoiesis, which is the generation of mature red blood cells in the bone marrow of whole organisms, has been described by a variety of mathematical models in the past decades. However, the in vitro erythropoiesis, which produces red blood cells in cultures, has received much less attention from the modelling community. In this paper, we propose the first mathematical model of in vitro erythropoiesis. We start by formulating different models and select the best one at fitting experimental data of in vitro erythropoietic differentiation obtained from chicken erythroid progenitor cells. It is based on a set of linear ODE, describing 3 hypothetical populations of cells at different stages of differentiation. We then compute confidence intervals for all of its parameters estimates, and conclude that our model is fully identifiable. Finally, we use this model to compute the effect of a chemical drug called Rapamycin, which affects all states of differentiation in the culture, and relate these effects to specific parameter variations. We provide the first model for the kinetics of in vitro cellular differentiation which is proven to be identifiable. It will serve as a basis for a model which will better account for the variability which is inherent to the experimental protocol used for the model calibration.
Angiogenesis, a formation of blood vessels from an existing vasculature, plays a key role in tumor growth and its progression into cancer. The lining of blood vessels consists of endothelial cells (ECs) which proliferate and migrate, allowing the capillaries to sprout towards the tumor to deliver the needed oxygen. Various treatments aiming to suppress or even inhibit angiogenesis have been explored. Mesenchymal stem cells (MSCs) have recently been undergoing development in cell-based therapy for cancer due to their ability to migrate towards the capillaries and induce the apoptosis of the ECs, causing capillary degeneration. However, further investigations in this direction are needed as it is usually difficult to preclinically assess the efficacy of such therapy. We develop a hybrid multiscale model that integrates molecular, cellular, tissue and extracellular components of tumor system to investigate angiogenesis and tumor growth under MSC-mediated therapy. Our simulations produce angiogenesis and vascular tumor growth profiles as observed in the experiments. Furthermore, the simulations show that the effectiveness of MSCs in inducing EC apoptosis is density dependent and its full effect is reached within several days after MSCs application. Quantitative agreements with experimental data indicate the predictive potential of our model for evaluating the efficacy of cell-based therapies targeting angiogenesis.
