Training signaling pathway maps to biochemical data with constrained fuzzy logic: quantitative analysis of liver cell responses to inflammatory stimuli

Predictive understanding of cell signaling network operation based on general prior knowledge but consistent with empirical data in a specific environmental context is a current challenge in computational biology. Recent work has demonstrated that Boolean logic can be used to create context-specific network models by training proteomic pathway maps to dedicated biochemical data; however, the Boolean formalism is restricted to characterizing protein species as either fully active or inactive. To advance beyond this limitation, we propose a novel form of fuzzy logic sufficiently flexible to model quantitative data but also sufficiently simple to efficiently construct models by training pathway maps on dedicated experimental measurements. Our new approach, termed constrained fuzzy logic (cFL), converts a prior knowledge network (obtained from literature or interactome databases) into a computable model that describes graded values of protein activation across multiple pathways. We train a cFL-converted network to experimental data describing hepatocytic protein activation by inflammatory cytokines and demonstrate the application of the resultant trained models for three important purposes: (a) generating experimentally testable biological hypotheses concerning pathway crosstalk, (b) establishing capability for quantitative prediction of protein activity, and (c) prediction and understanding of the cytokine release phenotypic response. Our methodology systematically and quantitatively trains a protein pathway map summarizing curated literature to context-specific biochemical data. This process generates a computable model yielding successful prediction of new test data and offering biological insight into complex datasets that are difficult to fully analyze by intuition alone.     

From proteomics toward systems biology: integration of different types of proteomics data into network models

Living organisms are comprised of various systems at different levels, i.e., organs, tissues, and cells. Each system carries out its diverse functions in response to environmental and genetic perturbations, by utilizing biological networks, in which nodal components, such as, DNA, mRNAs, proteins, and metabolites, closely interact with each other. Systems biology investigates such systems by producing comprehensive global data that represent different levels of biological information, i.e., at the DNA, mRNA, protein, or metabolite levels, and by integrating this data into network models that generate coherent hypotheses for given biological situations. This review presents a systems biology framework, called the 'Integrative Proteomics Data Analysis Pipeline' (IPDAP), which generates mechanistic hypotheses from network models reconstructed by integrating diverse types of proteomic data generated by mass spectrometry-based proteomic analyses. The devised framework includes a serial set of computational and network analysis tools. Here, we demonstrate its functionalities by applying these tools to several conceptual examples.     

Body size in ecological networks

Body size determines a host of species traits that can affect the structure and dynamics of food webs, and other ecological networks, across multiple scales of organization. Measuring body size provides a relatively simple means of encapsulating and condensing a large amount of the biological information embedded within an ecological network. Recently, important advances have been made by incorporating body size into theoretical models that explore food web stability, the patterning of energy fluxes, and responses to perturbations. Because metabolic constraints underpin body-size scaling relationships, metabolic theory offers a potentially useful new framework within which to develop novel models to describe the structure and functioning of ecological networks and to assess the probable consequences of biodiversity change.     

Studies of folding and misfolding using simplified models

Computer simulations are as vital to our studies of biological systems as experiments. They bridge and rationalize experimental observations, extend the experimental "field of view", which is often limited to a specific time or length scale, and, most importantly, provide novel insights into biological systems, offering hypotheses about yet-to-be uncovered phenomena. These hypotheses spur further experimental discoveries. Simplified molecular models have a special place in the field of computational biology. Branded as less accurate than all-atom protein models, they have offered what all-atom molecular dynamics simulations could not--the resolution of the length and time scales of biological phenomena. Not only have simplified models proven to be accurate in explaining or reproducing several biological phenomena, they have also offered a novel multiscale computational strategy for accessing a broad range of time and length scales upon integration with traditional all-atom simulations. Recent computer simulations of simplified models have shaken or advanced the established understanding of biological phenomena. It was demonstrated that simplified models can be as accurate as traditional molecular dynamics approaches in identifying native conformations of proteins. Their application to protein structure prediction yielded phenomenal accuracy in recapitulating native protein conformations. New studies that utilize the synergy of simplified protein models with all-atom models and experiments yielded novel insights into complex biological processes, such as protein folding, aggregation and the formation of large protein complexes.     

Mechanistic systems models of cell signaling networks: a case study of myocyte adrenergic regulation

Signal transduction networks coordinate a wide variety of cellular functions, including gene expression, metabolism, and cell fate processes. Understanding biological networks quantitatively is a major challenge to post-genomic biology, and mechanistic systems models will be crucial for this task. Here, we review approaches towards developing mechanistic systems models of established cell signaling networks. The ability of mechanistic system models to generate testable biological hypotheses and experimental strategies is discussed. As a case study of model development and analysis, we examined the functional roles of phospholamban, the L-type calcium channel, the ryanodine receptor, and troponin I phosphorylation upon β-adrenergic stimulation in the rat ventricular myocyte. Model analysis revealed that while protein kinase A-mediated phosphorylation of the ryanodine receptor greatly increases its calcium sensitivity, calcium autoregulation may adapt quickly by negating potential increases in contractility. Systematic combinations of in silico perturbations supported the conclusion that phospholamban phosphoregulation is the primary mechanism for increased sarcoplasmic reticulum load and calcium relaxation rate during β-adrenergic stimulation, while both phospholamban and the L-type calcium channel contribute to increased systolic calcium. Combined with detailed experimental studies, mechanistic systems models will be valuable for developing a quantitative understanding of cell signaling networks.     

Linear fuzzy gene network models obtained from microarray data by exhaustive search

Background  

Recent technological advances in high-throughput data collection allow for experimental study of increasingly complex systems on the scale of the whole cellular genome and proteome. Gene network models are needed to interpret the resulting large and complex data sets. Rationally designed perturbations (e.g., gene knock-outs) can be used to iteratively refine hypothetical models, suggesting an approach for high-throughput biological system analysis. We introduce an approach to gene network modeling based on a scalable linear variant of fuzzy logic: a framework with greater resolution than Boolean logic models, but which, while still semi-quantitative, does not require the precise parameter measurement needed for chemical kinetics-based modeling.     

A novel modeling approach for the “end-to-end” analysis of marine ecosystems

There is a growing demand for “end-to-end” models, which are modeling tools used to analyze and understand the fundamental complexities of marine ecosystems and processes emerging from the interaction of individuals from different trophic groups with respect to the physical environment and, even, human activity. These models are valuable quantitative tools for ecosystem-based management. To explore potential answers to complex questions regarding ecosystems using these models, it is necessary to incorporate classical ontogenic changes through the life cycle of target individuals, in addition to inherited behavioral strategies, as an additional differentiating aspect, particularly when the behavior has a direct impact on the ecosystem phenomena under study. However, it is difficult to combine different fine scale time and spatial granularities to infer animal behavior and ontogenic development. This complexity has kept these two levels of analysis separated, because most current tools do not have the required computational resources and advanced software architecture. To address this issue, we propose an individual-based modeling framework that is capable of handling and unifying the two experimental categories with a comprehensive biological and behavioral model that strictly adheres to the physiological functions of ingestion, growth, and metabolism of organisms. In addition, this model incorporates the exchange and transfer of mass and energy through local interactions at all trophic levels (lower to higher), the physical environment, and anthropogenic activity. For the framework to model short time events, such as classical predator–prey interactions, while also generating long-term ecosystem emergent properties, a special interleaving scheduling engine and physical space computer model was devised, which optimizes memory and processing resources. The framework was tested through several experiments with a three-population ecosystem containing up to 40 thousand organisms evolving inside a 200,000 m² simulation environment during 12,000 model-hours; yet, requiring only a few hours of program execution on a regular personal computer. The model included various environmental physical elements, such as several hundred shelters, the number of which can be easily modified in each experiment to simulate substrate degradation and its impact on populations. With the aid of the quantitative and qualitative tools provided by the model, it was possible to observe a coupling between prey and predator population dynamics. In conclusion, we confirmed that the end-to-end model developed here could successfully generate detailed specific hypotheses about fish behavior and quantify impacts on population dynamics.     

The Millennial model: in search of measurable pools and transformations for modeling soil carbon in the new century

Soil organic carbon (SOC) can be defined by measurable chemical and physical pools, such as mineral-associated carbon, carbon physically entrapped in aggregates, dissolved carbon, and fragments of plant detritus. Yet, most soil models use conceptual rather than measurable SOC pools. What would the traditional pool-based soil model look like if it were built today, reflecting the latest understanding of biological, chemical, and physical transformations in soils? We propose a conceptual model—the Millennial model—that defines pools as measurable entities. First, we discuss relevant pool definitions conceptually and in terms of the measurements that can be used to quantify pool size, formation, and destabilization. Then, we develop a numerical model following the Millennial model conceptual framework to evaluate against the Century model, a widely-used standard for estimating SOC stocks across space and through time. The Millennial model predicts qualitatively similar changes in total SOC in response to single factor perturbations when compared to Century, but different responses to multiple factor perturbations. We review important conceptual and behavioral differences between the Millennial and Century modeling approaches, and the field and lab measurements needed to constrain parameter values. We propose the Millennial model as a simple but comprehensive framework to model SOC pools and guide measurements for further model development.     

A computational approach for functional mapping of quantitative trait loci that regulate thermal performance curves

Whether and how thermal reaction norm is under genetic control is fundamental to understand the mechanistic basis of adaptation to novel thermal environments. However, the genetic study of thermal reaction norm is difficult because it is often expressed as a continuous function or curve. Here we derive a statistical model for dissecting thermal performance curves into individual quantitative trait loci (QTL) with the aid of a genetic linkage map. The model is constructed within the maximum likelihood context and implemented with the EM algorithm. It integrates the biological principle of responses to temperature into a framework for genetic mapping through rigorous mathematical functions established to describe the pattern and shape of thermal reaction norms. The biological advantages of the model lie in the decomposition of the genetic causes for thermal reaction norm into its biologically interpretable modes, such as hotter-colder, faster-slower and generalist-specialist, as well as the formulation of a series of hypotheses at the interface between genetic actions/interactions and temperature-dependent sensitivity. The model is also meritorious in statistics because the precision of parameter estimation and power of QTLdetection can be increased by modeling the mean-covariance structure with a small set of parameters. The results from simulation studies suggest that the model displays favorable statistical properties and can be robust in practical genetic applications. The model provides a conceptual platform for testing many ecologically relevant hypotheses regarding organismic adaptation within the Eco-Devo paradigm.     

Heterogeneity,Mixing, and the Spatial Scales of Mosquito-Borne Pathogen Transmission

The Ross-Macdonald model has dominated theory for mosquito-borne pathogen transmission dynamics and control for over a century. The model, like many other basic population models, makes the mathematically convenient assumption that populations are well mixed; i.e., that each mosquito is equally likely to bite any vertebrate host. This assumption raises questions about the validity and utility of current theory because it is in conflict with preponderant empirical evidence that transmission is heterogeneous. Here, we propose a new dynamic framework that is realistic enough to describe biological causes of heterogeneous transmission of mosquito-borne pathogens of humans, yet tractable enough to provide a basis for developing and improving general theory. The framework is based on the ecological context of mosquito blood meals and the fine-scale movements of individual mosquitoes and human hosts that give rise to heterogeneous transmission. Using this framework, we describe pathogen dispersion in terms of individual-level analogues of two classical quantities: vectorial capacity and the basic reproductive number, . Importantly, this framework explicitly accounts for three key components of overall heterogeneity in transmission: heterogeneous exposure, poor mixing, and finite host numbers. Using these tools, we propose two ways of characterizing the spatial scales of transmission—pathogen dispersion kernels and the evenness of mixing across scales of aggregation—and demonstrate the consequences of a model''s choice of spatial scale for epidemic dynamics and for estimation of , both by a priori model formulas and by inference of the force of infection from time-series data.     

A general structured model for a sequential hermaphrodite population

This paper introduces and analyzes a model of sequential hermaphroditism in the framework of continuously structured population models with sexual reproduction. The model is general in the sense that the birth, transition (from one sex to the other) and death processes of the population are given by arbitrary functions according to a biological meaningful hypotheses. The system is reduced to a single equation introducing the intrinsic sex-ratio subspace. The steady states are analyzed and illustrated for several cases. In particular, neglecting the competition for resources we have explicitly found a unique non-trivial equilibrium which is unstable.     

The emerging impact of tRNA modifications in the brain and nervous system

A remarkable number of neurodevelopmental disorders have been linked to defects in tRNA modifications. These discoveries place tRNA modifications in the spotlight as critical modulators of gene expression pathways that are required for proper organismal growth and development. Here, we discuss the emerging molecular and cellular functions of the diverse tRNA modifications linked to cognitive and neurological disorders. In particular, we describe how the structure and location of a tRNA modification influences tRNA folding, stability, and function. We then highlight how modifications in tRNA can impact multiple aspects of protein translation that are instrumental for maintaining proper cellular proteostasis. Importantly, we describe how perturbations in tRNA modification lead to a spectrum of deleterious biological outcomes that can disturb neurodevelopment and neurological function. Finally, we summarize the biological themes shared by the different tRNA modifications linked to cognitive disorders and offer insight into the future questions that remain to decipher the role of tRNA modifications. This article is part of a Special Issue entitled: mRNA modifications in gene expression control edited by Dr. Soller Matthias and Dr. Fray Rupert.     

Anticipatory Functions, Digital-Analog Forms and Biosemiotics: Integrating the Tools to Model Information and Normativity in Autonomous Biological Agents

We argue that living systems process information such that functionality emerges in them on a continuous basis. We then provide a framework that can explain and model the normativity of biological functionality. In addition we offer an explanation of the anticipatory nature of functionality within our overall approach. We adopt a Peircean approach to Biosemiotics, and a dynamical approach to Digital-Analog relations and to the interplay between different levels of functionality in autonomous systems, taking an integrative approach. We then apply the underlying biosemiotic logic to a particular biological system, giving a model of the B-Cell Receptor signaling system, in order to demonstrate how biosemiotic concepts can be used to build an account of biological information and functionality. Next we show how this framework can be used to explain and model more complex aspects of biological normativity, for example, how cross-talk between different signaling pathways can be avoided. Overall, we describe an integrated theoretical framework for the emergence of normative functions and, consequently, for the way information is transduced across several interconnected organizational levels in an autonomous system, and we demonstrate how this can be applied in real biological phenomena. Our aim is to open the way towards realistic tools for the modeling of information and normativity in autonomous biological agents.     

Modelling nutritional mutualisms: challenges and opportunities for data integration

Nutritional mutualisms are ancient, widespread, and profoundly influential in biological communities and ecosystems. Although much is known about these interactions, comprehensive answers to fundamental questions, such as how resource availability and structured interactions influence mutualism persistence, are still lacking. Mathematical modelling of nutritional mutualisms has great potential to facilitate the search for comprehensive answers to these and other fundamental questions by connecting the physiological and genomic underpinnings of mutualisms with ecological and evolutionary processes. In particular, when integrated with empirical data, models enable understanding of underlying mechanisms and generalisation of principles beyond the particulars of a given system. Here, we demonstrate how mathematical models can be integrated with data to address questions of mutualism persistence at four biological scales: cell, individual, population, and community. We highlight select studies where data has been or could be integrated with models to either inform model structure or test model predictions. We also point out opportunities to increase model rigour through tighter integration with data, and describe areas in which data is urgently needed. We focus on plant‐microbe systems, for which a wealth of empirical data is available, but the principles and approaches can be generally applied to any nutritional mutualism.     

Information-Theoretic Analysis of the Dynamics of an Executable Biological Model

To facilitate analysis and understanding of biological systems, large-scale data are often integrated into models using a variety of mathematical and computational approaches. Such models describe the dynamics of the biological system and can be used to study the changes in the state of the system over time. For many model classes, such as discrete or continuous dynamical systems, there exist appropriate frameworks and tools for analyzing system dynamics. However, the heterogeneous information that encodes and bridges molecular and cellular dynamics, inherent to fine-grained molecular simulation models, presents significant challenges to the study of system dynamics. In this paper, we present an algorithmic information theory based approach for the analysis and interpretation of the dynamics of such executable models of biological systems. We apply a normalized compression distance (NCD) analysis to the state representations of a model that simulates the immune decision making and immune cell behavior. We show that this analysis successfully captures the essential information in the dynamics of the system, which results from a variety of events including proliferation, differentiation, or perturbations such as gene knock-outs. We demonstrate that this approach can be used for the analysis of executable models, regardless of the modeling framework, and for making experimentally quantifiable predictions.