In Silico Biology最新文献

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.

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.

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 micro-tissue 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, the 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.

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.

 Quantitative modeling is quickly becoming an integral part of biology, due to the ability of mathematical models and computer simulations to generate insights and predict the behavior of living systems. Single-cell models can be incapable or misleading for inferring population dynamics, as they do not consider the interactions between cells via metabolites or physical contact, nor do they consider competition for limited resources such as nutrients or space. Here we examine methods that are commonly used to model and simulate cell populations. First, we cover simple models where analytic solutions are available, and then move on to more complex scenarios where computational methods are required. Overall, we present a summary of mathematical models used to describe cell population dynamics, which may aid future model development and highlights the importance of population modeling in biology.

Global level network analysis of molecular links is necessary for systems level view of complex diseases like cancer. Using genome-wide expression datasets, we constructed and compared gene co-expression based specific networks of pre-cancerous tumors (adenoma) and cancerous tumors (carcinoma) with paired normal networks to assess for any possible changes in network connectivity. Previously, loss of connectivity was reported as a characteristics of cancer samples. Here, we observed that pre-cancerous conditions also had significantly less connections than paired normal samples. We observed a loss of connectivity trend for colorectal adenoma, aldosterone producing adenoma and uterine leiomyoma. We also showed that the loss of connectivity trend is not specific to positive or negative correlation based networks. Differential hub genes, which were the most highly differentially less connected genes in tumor, were mostly different between different datasets. No common gene list could be defined which underlies the lower connectivity of tumor specific networks. Connectivity of colorectal cancer methylation targets was different from other genes. Extracellular space related terms were enriched in negative correlation based differential hubs and common methylation targets of colorectal carcinoma. Our results indicate a systems level change of lower connectivity as cells transform to not only cancer but also pre-cancerous conditions. This systems level behavior could not be attributed to a group of genes.

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.

A variety of mathematical models is used to describe and simulate the multitude of natural processes examined in life sciences. In this paper we present a scalable and adjustable foundation for the simulation of natural systems. Based on neighborhood relations in graphs and the complex interactions in cellular automata, the model uses recurrence relations to simulate changes on a mesoscopic scale. This implicit definition allows for the manipulation of every aspect of the model even during simulation. The definition of value rules ω facilitates the accumulation of change during time steps. Those changes may result from different physical, chemical or biological phenomena. Value rules can be combined into modules, which in turn can be used to create baseline models. Exemplarily, a value rule for the diffusion of chemical substances was designed and its applicability is demonstrated. Finally, the stability and accuracy of the solutions is analyzed.