Data-driven modelling of biological multi-scale processes

Biological processes involve a variety of spatial and temporal scales. A holistic understanding of many biological processes therefore requires multi-scale models which capture the relevant properties on all these scales. In this manuscript we review mathematical modelling approaches used to describe the individual spatial scales and how they are integrated into holistic models. We discuss the relation between spatial and temporal scales and the implication of that on multi-scale modelling. Based upon this overview over state-of-the-art modelling approaches, we formulate key challenges in mathematical and computational modelling of biological multi-scale and multi-physics processes. In particular, we considered the availability of analysis tools for multi-scale models and model-based multi-scale data integration. We provide a compact review of methods for model-based data integration and model-based hypothesis testing. Furthermore, novel approaches and recent trends are discussed, including computation time reduction using reduced order and surrogate models, which contribute to the solution of inference problems. We conclude the manuscript by providing a few ideas for the development of tailored multi-scale inference methods.Comment: This manuscript will appear in the Journal of Coupled Systems and Multiscale Dynamics (American Scientific Publishers

Propagation of kinetic uncertainties through a canonical topology of the TLR4 signaling network in different regions of biochemical reaction space

<p>Abstract</p> <p>Background</p> <p>Signal transduction networks represent the information processing systems that dictate which dynamical regimes of biochemical activity can be accessible to a cell under certain circumstances. One of the major concerns in molecular systems biology is centered on the elucidation of the robustness properties and information processing capabilities of signal transduction networks. Achieving this goal requires the establishment of causal relations between the design principle of biochemical reaction systems and their emergent dynamical behaviors.</p> <p>Methods</p> <p>In this study, efforts were focused in the construction of a relatively well informed, deterministic, non-linear dynamic model, accounting for reaction mechanisms grounded on standard mass action and Hill saturation kinetics, of the canonical reaction topology underlying Toll-like receptor 4 (TLR4)-mediated signaling events. This signaling mechanism has been shown to be deployed in macrophages during a relatively short time window in response to lypopolysaccharyde (LPS) stimulation, which leads to a rapidly mounted innate immune response. An extensive computational exploration of the biochemical reaction space inhabited by this signal transduction network was performed via local and global perturbation strategies. Importantly, a broad spectrum of biologically plausible dynamical regimes accessible to the network in widely scattered regions of parameter space was reconstructed computationally. Additionally, experimentally reported transcriptional readouts of target pro-inflammatory genes, which are actively modulated by the network in response to LPS stimulation, were also simulated. This was done with the main goal of carrying out an unbiased statistical assessment of the intrinsic robustness properties of this canonical reaction topology.</p> <p>Results</p> <p>Our simulation results provide convincing numerical evidence supporting the idea that a canonical reaction mechanism of the TLR4 signaling network is capable of performing information processing in a robust manner, a functional property that is independent of the signaling task required to be executed. Nevertheless, it was found that the robust performance of the network is not solely determined by its design principle (topology), but this may be heavily dependent on the network's current position in biochemical reaction space. Ultimately, our results enabled us the identification of key rate limiting steps which most effectively control the performance of the system under diverse dynamical regimes.</p> <p>Conclusions</p> <p>Overall, our <it>in silico </it>study suggests that biologically relevant and non-intuitive aspects on the general behavior of a complex biomolecular network can be elucidated only when taking into account a wide spectrum of dynamical regimes attainable by the system. Most importantly, this strategy provides the means for a suitable assessment of the inherent variational constraints imposed by the structure of the system when systematically probing its parameter space.</p

A multiscale systems perspective on cancer, immunotherapy, and Interleukin-12

Monoclonal antibodies represent some of the most promising molecular targeted immunotherapies. However, understanding mechanisms by which tumors evade elimination by the immune system of the host presents a significant challenge for developing effective cancer immunotherapies. The interaction of cancer cells with the host is a complex process that is distributed across a variety of time and length scales. The time scales range from the dynamics of protein refolding (i.e., microseconds) to the dynamics of disease progression (i.e., years). The length scales span the farthest reaches of the human body (i.e., meters) down to the range of molecular interactions (i.e., nanometers). Limited ranges of time and length scales are used experimentally to observe and quantify changes in physiology due to cancer. Translating knowledge obtained from the limited scales observed experimentally to predict patient response is an essential prerequisite for the rational design of cancer immunotherapies that improve clinical outcomes. In studying multiscale systems, engineers use systems analysis and design to identify important components in a complex system and to test conceptual understanding of the integrated system behavior using simulation. The objective of this review is to summarize interactions between the tumor and cell-mediated immunity from a multiscale perspective. Interleukin-12 and its role in coordinating antibody-dependent cell-mediated cytotoxicity is used illustrate the different time and length scale that underpin cancer immunoediting. An underlying theme in this review is the potential role that simulation can play in translating knowledge across scales

Information Processing under Physiological Pulsatile Stimulation in a G-protein Coupled Signaling Pathway

The cellular microenvironment is often dynamic, and several physiological ligands are released in pulsatile bursts. The main hypothesis driving this study is that cells are able to discern these time-varying dynamic inputs and must have evolved to exploit the temporal information available in their microenvironment to their advantage. Taking Muscarinic M3 (a G-protein coupled receptor)-mediated signaling as an example, this thesis explores how information is processed under pulsatile stimulation. Several experimental and computational approaches techniques including microfluidics, real-time multi-color fluorescence imaging of single cells, reaction kinetics modeling and information and noise analysis are implemented to gain mechanistic insights into the signaling circuit architecture. A major finding of this thesis is that receptor-mediated signaling forms a low pass filter while downstream calcium-induced NFAT (Nuclear Factor of Activated T-Lymphocytes, a transcription factor) nuclear translocation forms a high pass filter. The combination acts as a band-pass filter optimized for intermediate frequencies of stimulation. Sensitivity analysis shows that receptor and downstream kinetics determine critical features of the band-pass and that the band-pass may be shifted for different receptors or NFAT dynamics. Another important finding in this thesis is that for weak physiological inputs, cells exhibit apparent stochastic responses that can be explained within a deterministic framework. Computational analysis suggests that cells may utilize apparent stochasticity to enhance selectivity in downstream responses. This thesis also demonstrates that pulsatile inputs enhance information transfer downstream in noisy biochemical pathways. Finally, a microfluidic experimental method is developed to measure two microfluidic observables in the same cell, similar to a ‘two-reporter’ system, to estimate biochemical noise. Analysis with this method suggests that effect of drug action increases with increasing biochemical noise. Although this thesis focuses on one particular receptor and ligand, the conclusions from this work may be applied to several signaling systems. Investigation of band-pass processing may lead to gaining mechanistic insights into hidden or unknown regulatory motifs in several signaling pathways that are poorly understood. Using pulsatility to modulate selectivity and sensitivity of signaling response amidst biochemical noise provides tools to synthetic biologists and pharmacologists for developing enhanced lab-on-chip devices and pharmaceutical interventions.PHDBiophysicsUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/135745/1/msumit_1.pd

An introduction to the maximum entropy approach and its application to inference problems in biology

A cornerstone of statistical inference, the maximum entropy framework is being increasingly applied to construct descriptive and predictive models of biological systems, especially complex biological networks, from large experimental data sets. Both its broad applicability and the success it obtained in different contexts hinge upon its conceptual simplicity and mathematical soundness. Here we try to concisely review the basic elements of the maximum entropy principle, starting from the notion of ‘entropy’, and describe its usefulness for the analysis of biological systems. As examples, we focus specifically on the problem of reconstructing gene interaction networks from expression data and on recent work attempting to expand our system-level understanding of bacterial metabolism. Finally, we highlight some extensions and potential limitations of the maximum entropy approach, and point to more recent developments that are likely to play a key role in the upcoming challenges of extracting structures and information from increasingly rich, high-throughput biological data

Computation of Steady-State Probability Distributions in Stochastic Models of Cellular Networks

Cellular processes are “noisy”. In each cell, concentrations of molecules are subject to random fluctuations due to the small numbers of these molecules and to environmental perturbations. While noise varies with time, it is often measured at steady state, for example by flow cytometry. When interrogating aspects of a cellular network by such steady-state measurements of network components, a key need is to develop efficient methods to simulate and compute these distributions. We describe innovations in stochastic modeling coupled with approaches to this computational challenge: first, an approach to modeling intrinsic noise via solution of the chemical master equation, and second, a convolution technique to account for contributions of extrinsic noise. We show how these techniques can be combined in a streamlined procedure for evaluation of different sources of variability in a biochemical network. Evaluation and illustrations are given in analysis of two well-characterized synthetic gene circuits, as well as a signaling network underlying the mammalian cell cycle entry

Dynamic optimization of metabolic networks coupled with gene expression

The regulation of metabolic activity by tuning enzyme expression levels is crucial to sustain cellular growth in changing environments. Metabolic networks are often studied at steady state using constraint-based models and optimization techniques. However, metabolic adaptations driven by changes in gene expression cannot be analyzed by steady state models, as these do not account for temporal changes in biomass composition. Here we present a dynamic optimization framework that integrates the metabolic network with the dynamics of biomass production and composition, explicitly taking into account enzyme production costs and enzymatic capacity. In contrast to the established dynamic flux balance analysis, our approach allows predicting dynamic changes in both the metabolic fluxes and the biomass composition during metabolic adaptations. We applied our algorithm in two case studies: a minimal nutrient uptake network, and an abstraction of core metabolic processes in bacteria. In the minimal model, we show that the optimized uptake rates reproduce the empirical Monod growth for bacterial cultures. For the network of core metabolic processes, the dynamic optimization algorithm predicted commonly observed metabolic adaptations, such as a diauxic switch with a preference ranking for different nutrients, re-utilization of waste products after depletion of the original substrate, and metabolic adaptation to an impending nutrient depletion. These examples illustrate how dynamic adaptations of enzyme expression can be predicted solely from an optimization principle

The 2019 Mathematical Oncology Roadmap.

Whether the nom de guerre is Mathematical Oncology, Computational or Systems Biology, Theoretical Biology, Evolutionary Oncology, Bioinformatics, or simply Basic Science, there is no denying that mathematics continues to play an increasingly prominent role in cancer research. Mathematical Oncology-defined here simply as the use of mathematics in cancer research-complements and overlaps with a number of other fields that rely on mathematics as a core methodology. As a result, Mathematical Oncology has a broad scope, ranging from theoretical studies to clinical trials designed with mathematical models. This Roadmap differentiates Mathematical Oncology from related fields and demonstrates specific areas of focus within this unique field of research. The dominant theme of this Roadmap is the personalization of medicine through mathematics, modelling, and simulation. This is achieved through the use of patient-specific clinical data to: develop individualized screening strategies to detect cancer earlier; make predictions of response to therapy; design adaptive, patient-specific treatment plans to overcome therapy resistance; and establish domain-specific standards to share model predictions and to make models and simulations reproducible. The cover art for this Roadmap was chosen as an apt metaphor for the beautiful, strange, and evolving relationship between Two Beasts: mathematics and cancer.NIH (R01CA16437, R01NS060752, U54CA210180, U54CA143970, U54193489, U01CA220378)James S. McDonnell FoundationBen & Catherine Ivy Foundatio

Aging and Health

Aging is a major risk factor for chronic diseases, which in turn can provide information about the aging of a biological system. This publication serves as an introduction to systems biology and its application to biological aging. Key pathways and processes that impinge on aging are reviewed, and how they contribute to health and disease during aging is discussed. The evolution of this situation is analyzed, and the consequences for the study of genetic effects on aging are presented. Epigenetic programming of aging, as a continuation of development, creates an interface between the genome and the environment. New research into the gut microbiome describes how this interface may operate in practice with marked consequences for a variety of disorders. This analysis is bolstered by a view of the aging organism as a whole, with conclusions about the mechanisms underlying resilience of the organism to change, and is expanded with a discussion of circadian rhythms in aging