Optimization in computational systems biology

Nutrient-epigenome interactions: Implications for personalized nutrition against aging-associated diseases

Aging is a multifaceted process involving genetic and environmental interactions often resulting in epigenetic changes, potentially leading to aging-related diseases. Various strategies, like dietary interventions and calorie restrictions, have been employed to modify these epigenetic landscapes. A burgeoning field of interest focuses on the role of microbiota in human health, emphasizing system biology and computational approaches. These methods help decipher the intricate interplay between diet and gut microbiota, facilitating the creation of personalized nutrition strategies. In this review, we analysed the mechanisms related to nutritional interventions while highlighting the influence of dietary strategies, like calorie restriction and intermittent fasting, on microbial composition and function. We explore how gut microbiota affects the efficacy of interventions using tools like multi-omics data integration, network analysis, and machine learning. These tools enable us to pinpoint critical regulatory elements and generate individualized models for dietary responses. Lastly, we emphasize the need for a deeper comprehension of nutrient-epigenome interactions and the potential of personalized nutrition informed by individual genetic and epigenetic profiles. As knowledge and technology advance, dietary epigenetics stands on the cusp of reshaping our strategy against aging and related diseases.

CHOmpact: A reduced metabolic model of Chinese hamster ovary cells with enhanced interpretability

Metabolic modeling has emerged as a key tool for the characterization of biopharmaceutical cell culture processes. Metabolic models have also been instrumental in identifying genetic engineering targets and developing feeding strategies that optimize the growth and productivity of Chinese hamster ovary (CHO) cells. Despite their success, metabolic models of CHO cells still present considerable challenges. Genome-scale metabolic models (GeMs) of CHO cells are very large (>6000 reactions) and are difficult to constrain to yield physiologically consistent flux distributions. The large scale of GeMs also makes the interpretation of their outputs difficult. To address these challenges, we have developed CHOmpact, a reduced metabolic network that encompasses 101 metabolites linked through 144 reactions. Our compact reaction network allows us to deploy robust, nonlinear optimization and ensure that the computed flux distributions are physiologically consistent. Furthermore, our CHOmpact model delivers enhanced interpretability of simulation results and has allowed us to identify the mechanisms governing shifts in the anaplerotic consumption of asparagine and glutamate as well as an important mechanism of ammonia detoxification within mitochondria. CHOmpact, thus, addresses key challenges of large-scale metabolic models and will serve as a platform to develop dynamic metabolic models for the control and optimization of biopharmaceutical cell culture processes.

Bayesian Optimization for Design of Multiscale Biological Circuits

Recent advances in synthetic biology have enabled the construction of molecular circuits that operate across multiple scales of cellular organization, such as gene regulation, signaling pathways, and cellular metabolism. Computational optimization can effectively aid the design process, but current methods are generally unsuited for systems with multiple temporal or concentration scales, as these are slow to simulate due to their numerical stiffness. Here, we present a machine learning method for the efficient optimization of biological circuits across scales. The method relies on Bayesian optimization, a technique commonly used to fine-tune deep neural networks, to learn the shape of a performance landscape and iteratively navigate the design space toward an optimal circuit. This strategy allows the joint optimization of both circuit architecture and parameters, and provides a feasible approach to solve a highly nonconvex optimization problem in a mixed-integer input space. We illustrate the applicability of the method on several gene circuits for controlling biosynthetic pathways with strong nonlinearities, multiple interacting scales, and using various performance objectives. The method efficiently handles large multiscale problems and enables parametric sweeps to assess circuit robustness to perturbations, serving as an efficient in silico screening method prior to experimental implementation.

Automatic Generation of SBML Kinetic Models from Natural Language Texts Using GPT

Kinetic modeling is an essential tool in systems biology research, enabling the quantitative analysis of biological systems and predicting their behavior. However, the development of kinetic models is a complex and time-consuming process. In this article, we propose a novel approach called KinModGPT, which generates kinetic models directly from natural language text. KinModGPT employs GPT as a natural language interpreter and Tellurium as an SBML generator. We demonstrate the effectiveness of KinModGPT in creating SBML kinetic models from complex natural language descriptions of biochemical reactions. KinModGPT successfully generates valid SBML models from a range of natural language model descriptions of metabolic pathways, protein-protein interaction networks, and heat shock response. This article demonstrates the potential of KinModGPT in kinetic modeling automation.

Approximating scaffold printability utilizing computational methods

Bioprinting facilitates the generation of complex, three-dimensional (3D), cell-based constructs for various applications. Although multiple bioprinting technologies have been developed, extrusion-based systems have become the dominant technology due to the diversity of materials (bioinks) that can be utilized, either individually or in combination. However, each bioink has unique material properties and extrusion characteristics that affect bioprinting utility, accuracy, and precision. Here, we have extended our previous work to achieve high precision (i.e. repeatability) and printability across samples by optimizing bioink-specific printing parameters. Specifically, we hypothesized that a fuzzy inference system (FIS) could be used as a computational method to address the imprecision in 3D bioprinting test data and uncover the optimal printing parameters for a specific bioink that result in high accuracy and precision. To test this hypothesis, we have implemented a FIS model consisting of four inputs (bioink concentration, printing flow rate, speed, and temperature) and two outputs to quantify the precision (scaffold bioprinted linewidth variance) and printability. We validate our use of the bioprinting precision index with both standard and normalized printability factors. Finally, we utilize optimized printing parameters to bioprint scaffolds containing up to 30 × 106cells ml-1with high printability and precision. In total, our results indicate that computational methods are a cost-efficient measure to improve the precision and robustness of extrusion 3D bioprinting.

PhysiCOOL: A generalized framework for model Calibration and Optimization Of modeLing projects

In silico models of biological systems are usually very complex and rely on a large number of parameters describing physical and biological properties that require validation. As such, parameter space exploration is an essential component of computational model development to fully characterize and validate simulation results. Experimental data may also be used to constrain parameter space (or enable model calibration) to enhance the biological relevance of model parameters. One widely used computational platform in the mathematical biology community is PhysiCell, which provides a standardized approach to agent-based models of biological phenomena at different time and spatial scales. Nonetheless, one limitation of PhysiCell is the lack of a generalized approach for parameter space exploration and calibration that can be run without high-performance computing access. Here, we present PhysiCOOL, an open-source Python library tailored to create standardized calibration and optimization routines for PhysiCell models.

ASAS-NANP symposium: Mathematical Modeling in Animal Nutrition: The power of identifiability analysis for dynamic modeling in animal science:a practitioner approach

Constructing dynamic mathematical models of biological systems requires estimating unknown parameters from available experimental data, usually using a statistical fitting procedure. This procedure is usually called parameter identification, parameter estimation, model fitting, or model calibration. In animal science, parameter identification is often performed without analytic considerations on the possibility of determining unique values of the model parameters. These analytical studies are related to the mathematical property of structural identifiability, which refers to the theoretical ability to recover unique values of the model parameters from the measures defined in an experimental setup and use the model structure as the sole basis. The structural identifiability analysis is a powerful tool for model construction because it informs whether the parameter identification problem is well-posed (i.e., the problem has a unique solution). Structural identifiability analysis is helpful to determine which actions (e.g., model reparameterization, choice of new data measurements, and change of the model structure) are needed to render the model parameters identifiable (when possible). The mathematical technicalities associated with structural identifiability analysis are very sophisticated. However, the development of dedicated, freely available software tools enables the application of identifiability analysis without needing to be an expert in mathematics and computer programming. We refer to such a non-expert user as a practitioner for hands-on purposes. However, a practitioner should be familiar with the model construction and software implementation process. In this paper, we propose to adopt a practitioner approach that takes advantage of available software tools to integrate identifiability analysis in the modeling practice in the animal science field. The application of structural identifiability implies switching our regard of the parameter identification problem as a downstream process (after data collection) to an upstream process (before data collection) where experiment design is applied to guarantee identifiability. This upstream approach will substantially improve the workflow of model construction toward robust and valuable models in animal science. Illustrative examples with different levels of complexity support our work. The source codes of the examples were provided for learning purposes and to promote open science practices.

MLAGO: machine learning-aided global optimization for Michaelis constant estimation of kinetic modeling

BACKGROUND

Kinetic modeling is a powerful tool for understanding the dynamic behavior of biochemical systems. For kinetic modeling, determination of a number of kinetic parameters, such as the Michaelis constant (K_m), is necessary, and global optimization algorithms have long been used for parameter estimation. However, the conventional global optimization approach has three problems: (i) It is computationally demanding. (ii) It often yields unrealistic parameter values because it simply seeks a better model fitting to experimentally observed behaviors. (iii) It has difficulty in identifying a unique solution because multiple parameter sets can allow a kinetic model to fit experimental data equally well (the non-identifiability problem).

RESULTS

To solve these problems, we propose the Machine Learning-Aided Global Optimization (MLAGO) method for K_m estimation of kinetic modeling. First, we use a machine learning-based K_m predictor based only on three factors: EC number, KEGG Compound ID, and Organism ID, then conduct a constrained global optimization-based parameter estimation by using the machine learning-predicted K_m values as the reference values. The machine learning model achieved relatively good prediction scores: RMSE = 0.795 and R² = 0.536, making the subsequent global optimization easy and practical. The MLAGO approach reduced the error between simulation and experimental data while keeping K_m values close to the machine learning-predicted values. As a result, the MLAGO approach successfully estimated K_m values with less computational cost than the conventional method. Moreover, the MLAGO approach uniquely estimated K_m values, which were close to the measured values.

CONCLUSIONS

MLAGO overcomes the major problems in parameter estimation, accelerates kinetic modeling, and thus ultimately leads to better understanding of complex cellular systems. The web application for our machine learning-based K_m predictor is accessible at https://sites.google.com/view/kazuhiro-maeda/software-tools-web-apps , which helps modelers perform MLAGO on their own parameter estimation tasks.

Optimization of synthetic molecular reporters for a mesenchymal glioblastoma transcriptional program by integer programing

MOTIVATION

A recent approach to perform genetic tracing of complex biological problems involves the generation of synthetic deoxyribonucleic acid (DNA) probes that specifically mark cells with a phenotype of interest. These synthetic locus control regions (sLCRs), in turn, drive the expression of a reporter gene, such as fluorescent protein. To build functional and specific sLCRs, it is critical to accurately select multiple bona fide cis-regulatory elements from the target cell phenotype cistrome. This selection occurs by maximizing the number and diversity of transcription factors (TFs) within the sLCR, yet the size of the final sLCR should remain limited.

RESULTS

In this work, we discuss how optimization, in particular integer programing, can be used to systematically address the construction of a specific sLCR and optimize pre-defined properties of the sLCR. Our presented instance of a linear optimization problem maximizes the activation potential of the sLCR such that its size is limited to a pre-defined length and a minimum number of all TFs deemed sufficiently characteristic for the phenotype of interest is covered. We generated an sLCR to trace the mesenchymal glioblastoma program in patients by solving our corresponding linear program with the software optimizer Gurobi. Considering the binding strength of transcription factor binding sites (TFBSs) with their TFs as a proxy for activation potential, the optimized sLCR scores similarly to an sLCR experimentally validated in vivo, and is smaller in size while having the same coverage of TFBSs.

AVAILABILITY AND IMPLEMENTATION

We provide a Python implementation of the presented framework in the Supplementary Material with which an optimal selection of cis-regulatory elements can be calculated once the target set of TFs and their binding strength with their TFBSs is known.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Understanding mitochondria and the utility of optimization as a canonical framework for identifying and modeling mitochondrial pathways

The goal of this paper is to provide an overview of our current understanding of mitochondrial function as a framework to motivate the hypothesis that mitochondrial behavior is governed by optimization principles that are constrained by the laws of the physical and biological sciences. Then, mathematical optimization tools can generally be useful to model some of these processes under reasonable assumptions and limitations. We are specifically interested in optimizations via variational methods, which are briefly summarized. Within such an optimization framework, we suggest that the numerous mechanical instigators of cell and intracellular functioning can be modeled utilizing some of the principles of mechanics that govern engineered systems, as well as by the frequently observed feedback and feedforward mechanisms that coordinate the multitude of processes within cells. These mechanical aspects would need to be coupled to governing biochemical rules. Of course, biological systems are significantly more complex than engineered systems, and require considerably more experimentation to ascertain and characterize parameters and subsequent behavior. That complexity requires well-defined limitations and assumptions for any derived models. Optimality is being motivated as a framework to help us understand how cellular decisions are made, especially those that transition between physiological behaviors and dysfunctions along pathophysiological pathways. We elaborate on our interpretation of optimality and cellular decision making within the body of this paper, as we revisit these ideas in the numerous different contexts of mitochondrial functions.

Mathematical Modelling in Plant Synthetic Biology

Mathematical modelling techniques are integral to current research in plant synthetic biology. Modelling approaches can provide mechanistic understanding of a system, allowing predictions of behaviour and thus providing a tool to help design and analyse biological circuits. In this chapter, we provide an overview of mathematical modelling methods and their significance for plant synthetic biology. Starting with the basics of dynamics, we describe the process of constructing a model over both temporal and spatial scales and highlight crucial approaches, such as stochastic modelling and model-based design. Next, we focus on the model parameters and the techniques required in parameter analysis. We then describe the process of selecting a model based on tests and criteria and proceed to methods that allow closer analysis of the system's behaviour. Finally, we highlight the importance of uncertainty in modelling approaches and how to deal with a lack of knowledge, noisy data, and biological variability; all aspects that play a crucial role in the cooperation between the experimental and modelling components. Overall, this chapter aims to illustrate the importance of mathematical modelling in plant synthetic biology, providing an introduction for those researchers who are working with or working on modelling techniques.

Optimization and feedback control system of dilution rate for 1,3-propanediol in two-stage fermentation: A theoretical study

In utilizing glycerol to produce 1,3-propanediol by microbial fermentation, the problems of low utilization rate and poor production performance need to be addressed. Based on the analysis of a mathematical model for 1,3-propanediol production from glycerol by Klebsiella pneumoniae, this study theoretically investigated the effects of the dilution rate and the initial glycerol concentration in a two-stage fermentation process and the feasibility of applying the feedback control methods. First, the optimal operation conditions of initial glycerol concentration and dilution rate were obtained. Through the use of feedback control theory, a control strategy for dilution rate was designed and optimized to shorten the settling time (time required for fermentation to reach stability) from 60.92 to 36.68 h for the first reactor, and from 53.66 to 22.68 h for the second reactor. In addition, the yield of 1,3-propanediol in both two reactors reached up to 0.5 g·g^-1 . The simulation results indicated that the feedback control strategy for dilution rate increased the product concentration, reduced the residual glycerol in the fermentation broth, and greatly improved the performance of the fermentation. A feeding strategy of automatic control for dilution rate has been established and will be applied as an effective guiding scheme in automatic continuous fermentations for production of 1,3-propanediol.

Optimal proteome allocation and the temperature dependence of microbial growth laws

Although the effect of temperature on microbial growth has been widely studied, the role of proteome allocation in bringing about temperature-induced changes remains elusive. To tackle this problem, we propose a coarse-grained model of microbial growth, including the processes of temperature-sensitive protein unfolding and chaperone-assisted (re)folding. We determine the proteome sector allocation that maximizes balanced growth rate as a function of nutrient limitation and temperature. Calibrated with quantitative proteomic data for Escherichia coli, the model allows us to clarify general principles of temperature-dependent proteome allocation and formulate generalized growth laws. The same activation energy for metabolic enzymes and ribosomes leads to an Arrhenius increase in growth rate at constant proteome composition over a large range of temperatures, whereas at extreme temperatures resources are diverted away from growth to chaperone-mediated stress responses. Our approach points at risks and possible remedies for the use of ribosome content to characterize complex ecosystems with temperature variation.

Using optimal control to understand complex metabolic pathways

BACKGROUND

Optimality principles have been used to explain the structure and behavior of living matter at different levels of organization, from basic phenomena at the molecular level, up to complex dynamics in whole populations. Most of these studies have assumed a single-criteria approach. Such optimality principles have been justified from an evolutionary perspective. In the context of the cell, previous studies have shown how dynamics of gene expression in small metabolic models can be explained assuming that cells have developed optimal adaptation strategies. Most of these works have considered rather simplified representations, such as small linear pathways, or reduced networks with a single branching point, and a single objective for the optimality criteria.

RESULTS

Here we consider the extension of this approach to more realistic scenarios, i.e. biochemical pathways of arbitrary size and structure. We first show that exploiting optimality principles for these networks poses great challenges due to the complexity of the associated optimal control problems. Second, in order to surmount such challenges, we present a computational framework which has been designed with scalability and efficiency in mind, including mechanisms to avoid the most common pitfalls. Third, we illustrate its performance with several case studies considering the central carbon metabolism of S. cerevisiae and B. subtilis. In particular, we consider metabolic dynamics during nutrient shift experiments.

CONCLUSIONS

We show how multi-objective optimal control can be used to predict temporal profiles of enzyme activation and metabolite concentrations in complex metabolic pathways. Further, we also show how to consider general cost/benefit trade-offs. In this study we have considered metabolic pathways, but this computational framework can also be applied to analyze the dynamics of other complex pathways, such as signal transduction or gene regulatory networks.

Parameterization of mechanistic models from qualitative data using an efficient optimal scaling approach

Quantitative dynamical models facilitate the understanding of biological processes and the prediction of their dynamics. These models usually comprise unknown parameters, which have to be inferred from experimental data. For quantitative experimental data, there are several methods and software tools available. However, for qualitative data the available approaches are limited and computationally demanding. Here, we consider the optimal scaling method which has been developed in statistics for categorical data and has been applied to dynamical systems. This approach turns qualitative variables into quantitative ones, accounting for constraints on their relation. We derive a reduced formulation for the optimization problem defining the optimal scaling. The reduced formulation possesses the same optimal points as the established formulation but requires less degrees of freedom. Parameter estimation for dynamical models of cellular pathways revealed that the reduced formulation improves the robustness and convergence of optimizers. This resulted in substantially reduced computation times. We implemented the proposed approach in the open-source Python Parameter EStimation TOolbox (pyPESTO) to facilitate reuse and extension. The proposed approach enables efficient parameterization of quantitative dynamical models using qualitative data.

Calibration of ionic and cellular cardiac electrophysiology models

Cardiac electrophysiology models are among the most mature and well-studied mathematical models of biological systems. This maturity is bringing new challenges as models are being used increasingly to make quantitative rather than qualitative predictions. As such, calibrating the parameters within ion current and action potential (AP) models to experimental data sets is a crucial step in constructing a predictive model. This review highlights some of the fundamental concepts in cardiac model calibration and is intended to be readily understood by computational and mathematical modelers working in other fields of biology. We discuss the classic and latest approaches to calibration in the electrophysiology field, at both the ion channel and cellular AP scales. We end with a discussion of the many challenges that work to date has raised and the need for reproducible descriptions of the calibration process to enable models to be recalibrated to new data sets and built upon for new studies. This article is categorized under: Analytical and Computational Methods > Computational Methods Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Cellular Models.

Efficient exact inference for dynamical systems with noisy measurements using sequential approximate Bayesian computation

MOTIVATION

Approximate Bayesian computation (ABC) is an increasingly popular method for likelihood-free parameter inference in systems biology and other fields of research, as it allows analyzing complex stochastic models. However, the introduced approximation error is often not clear. It has been shown that ABC actually gives exact inference under the implicit assumption of a measurement noise model. Noise being common in biological systems, it is intriguing to exploit this insight. But this is difficult in practice, as ABC is in general highly computationally demanding. Thus, the question we want to answer here is how to efficiently account for measurement noise in ABC.

RESULTS

We illustrate exemplarily how ABC yields erroneous parameter estimates when neglecting measurement noise. Then, we discuss practical ways of correctly including the measurement noise in the analysis. We present an efficient adaptive sequential importance sampling-based algorithm applicable to various model types and noise models. We test and compare it on several models, including ordinary and stochastic differential equations, Markov jump processes and stochastically interacting agents, and noise models including normal, Laplace and Poisson noise. We conclude that the proposed algorithm could improve the accuracy of parameter estimates for a broad spectrum of applications.

AVAILABILITY AND IMPLEMENTATION

The developed algorithms are made publicly available as part of the open-source python toolbox pyABC (https://github.com/icb-dcm/pyabc).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Systems Genetics for Mechanistic Discovery in Heart Diseases

Cardiovascular diseases are the leading cause of death worldwide. Complex diseases with highly heterogenous disease progression among patient populations, cardiovascular diseases feature multifactorial contributions from both genetic and environmental stressors. Despite significant effort utilizing multiple approaches from molecular biology to genome-wide association studies, the genetic landscape of cardiovascular diseases, particularly for the nonfamilial forms of heart failure, is still poorly understood. In the past decade, systems-level approaches based on omics technologies have become an important approach for the study of complex traits in large populations. These advances create opportunities to integrate genetic variation with other biological layers to identify and prioritize candidate genes, understand pathogenic pathways, and elucidate gene-gene and gene-environment interactions. In this review, we will highlight some of the recent progress made using systems genetics approaches to uncover novel mechanisms and molecular bases of cardiovascular pathophysiological manifestations. The key technology and data analysis platforms necessary to implement systems genetics will be described, and the current major challenges and future directions will also be discussed. For complex cardiovascular diseases, such as heart failure, systems genetics represents a powerful strategy to obtain mechanistic insights and to develop individualized diagnostic and therapeutic regiments, paving the way for precision cardiovascular medicine.

Medusa: Software to build and analyze ensembles of genome-scale metabolic network reconstructions

Uncertainty in the structure and parameters of networks is ubiquitous across computational biology. In constraint-based reconstruction and analysis of metabolic networks, this uncertainty is present both during the reconstruction of networks and in simulations performed with them. Here, we present Medusa, a Python package for the generation and analysis of ensembles of genome-scale metabolic network reconstructions. Medusa builds on the COBRApy package for constraint-based reconstruction and analysis by compressing a set of models into a compact ensemble object, providing functions for the generation of ensembles using experimental data, and extending constraint-based analyses to ensemble scale. We demonstrate how Medusa can be used to generate ensembles and perform ensemble simulations, and how machine learning can be used in conjunction with Medusa to guide the curation of genome-scale metabolic network reconstructions. Medusa is available under the permissive MIT license from the Python Packaging Index (https://pypi.org) and from github (https://github.com/opencobra/Medusa), and comprehensive documentation is available at https://medusa.readthedocs.io/en/latest.

The FLAME-accelerated signalling tool (FaST) for facile parallelisation of flexible agent-based models of cell signalling

Agent-based modelling is particularly adept at modelling complex features of cell signalling pathways, where heterogeneity, stochastic and spatial effects are important, thus increasing our understanding of decision processes in biology in such scenarios. However, agent-based modelling often is computationally prohibitive to implement. Parallel computing, either on central processing units (CPUs) or graphical processing units (GPUs), can provide a means to improve computational feasibility of agent-based applications but generally requires specialist coding knowledge and extensive optimisation. In this paper, we address these challenges through the development and implementation of the FLAME-accelerated signalling tool (FaST), a software that permits easy creation and parallelisation of agent-based models of cell signalling, on CPUs or GPUs. FaST incorporates validated new agent-based methods, for accurate modelling of reaction kinetics and, as proof of concept, successfully converted an ordinary differential equation (ODE) model of apoptosis execution into an agent-based model. We finally parallelised this model through FaST on CPUs and GPUs resulting in an increase in performance of 5.8× (16 CPUs) and 53.9×, respectively. The FaST takes advantage of the communicating X-machine approach used by FLAME and FLAME GPU to allow easy alteration or addition of functionality to parallel applications, but still includes inherent parallelisation optimisation. The FaST, therefore, represents a new and innovative tool to easily create and parallelise bespoke, robust, agent-based models of cell signalling.

Global optimization of the Michaelis-Menten parameters using physiologically-based pharmacokinetic (PBPK) modeling and chloroform vapor uptake data in F344 rats

Objective: To quantify metabolism, a physiologically based pharmacokinetic (PBPK) model for a volatile compound can be calibrated with the closed chamber (i.e. vapor uptake) inhalation data. Here, we introduce global optimization as a novel component of the predictive process and use it to illustrate a procedure for metabolic parameter estimation.Materials and methods: Male F344 rats were exposed in vapor uptake chambers to initial concentrations of 100, 500, 1000, and 3000 ppm chloroform. Chamber time-course data from these experiments, in combination with optimization using a chemical-specific PBPK model, were used to estimate Michaelis-Menten metabolic constants. Matlab^® simulation software was used to integrate the mass balance equations and to perform the global optimizations using MEIGO (MEtaheuristics for systems biology and bIoinformatics Global Optimization - Version 64 bit, R2016A), a toolbox written for Matlab®. The cost function used the chamber time-course data and least squares to minimize the difference between data and simulation values.Results and discussion: The final values estimated for V_max (maximum metabolic rate) and K_m (affinity constant) were 1.2 mg/h and a range between 0.0005 and 0.6 mg/L, respectively. Also, cost function plots were used to analyze the dose-dependent capacity to estimate V_max and K_m within the experimental range used. Sensitivity analysis was used to assess identifiability for both parameters and show these kinetic data may not be sufficient to identify K_m.Conclusion: In summary, this work should help toxicologists interested in optimization techniques understand the overall process employed when calibrating metabolic parameters in a PBPK model with inhalation data.

Trends in mathematical modeling of host-pathogen interactions

Pathogenic microorganisms entail enormous problems for humans, livestock, and crop plants. A better understanding of the different infection strategies of the pathogens enables us to derive optimal treatments to mitigate infectious diseases or develop vaccinations preventing the occurrence of infections altogether. In this review, we highlight the current trends in mathematical modeling approaches and related methods used for understanding host-pathogen interactions. Since these interactions can be described on vastly different temporal and spatial scales as well as abstraction levels, a variety of computational and mathematical approaches are presented. Particular emphasis is placed on dynamic optimization, game theory, and spatial modeling, as they are attracting more and more interest in systems biology. Furthermore, these approaches are often combined to illuminate the complexities of the interactions between pathogens and their host. We also discuss the phenomena of molecular mimicry and crypsis as well as the interplay between defense and counter defense. As a conclusion, we provide an overview of method characteristics to assist non-experts in their decision for modeling approaches and interdisciplinary understanding.

Mathematical Modeling of Plant Metabolism in a Changing Temperature Regime

Changes in environmental temperature regimes significantly affect plant growth, development and reproduction. Within a multigenic process termed acclimation, many plant species of the temperate region are able to adjust their metabolism to low and high temperature. Temperature-induced metabolic reprogramming is a nonlinear process affecting numerous enzyme kinetic reactions and pathways. The analysis of metabolic reprogramming during temperature acclimation is essentially supported by mathematical modeling which enables the study of nonlinear enzyme kinetics in context of metabolic networks and pathway regulation. This chapter introduces mathematical modeling of plant metabolism during a dynamic environmental temperature regime. A focus is laid on kinetic modeling and thermodynamic constraints.

Kinetic model optimization and its application to mitigating the Warburg effect through multiple enzyme alterations

Pathway engineering is a powerful tool in biotechnological and clinical applications. However, many phenomena cannot be rewired with a single enzyme change, and in a complex network like energy metabolism, the selection of combinations of targets to engineer is a daunting task. To facilitate this process, we have developed an optimization framework and applied it to a mechanistic kinetic model of energy metabolism. We then identified combinations of enzyme alternations that led to the elimination of the Warburg effect seen in the metabolism of cancer cells and cell lines, a phenomenon coupling rapid proliferation to lactate production. Typically, optimization approaches use integer variables to achieve the desired flux redistribution with a minimum number of altered genes. This framework uses convex penalty terms to replace these integer variables and improve computational tractability. Optimal solutions are identified which substantially reduce or eliminate lactate production while maintaining the requirements for cellular proliferation using three or more enzymes.

Ensemble optimization of microbial conversion of glycerol into 1, 3-propanediol by Klebsiella pneumoniae

Using mathematical model and computer simulation to predict biological processes and optimize the target production is an important strategy for optimizing fermentation process. However, the inherent uncertainty of the kinetic model severely limits the predictive capability. In this study, optimize target production, such as productivity and yield of 1, 3-propanediol produced by Klebsiella pneumoniae using glycerol as substrate, the ensemble modeling approach was used to reduce the model's uncertainty for fermentation process as much as possible, and effectively improve its prediction performance. Firstly, through sensitivity analysis, the parameters having significant influence on the model were determined as the adjustable parameters for the ensemble modeling. After comparison, the appropriate threshold coefficient of the model error was determined, and the sampling method was used to generate as many equivalent parameter sets as possible. Each set of parameters was separately applied for the simulation, and all the predicted values were integrated for the weighted average. Therefore, the expected value of the prediction was obtained. Compared with the traditional simulation using single parameter set, the ensemble modeling method achieved the lower relative error between the prediction and the experimental value and the greatly improved model prediction performance. Moreover, the optimal productivity and yield of 1, 3-propanediol and the corresponding operating conditions were obtained, respectively. The ensemble modeling approach effectively compensates for the uncertainties of the model, making its prediction performance more practical, which is important for computer simulations to predict and guide the actual production process.

Identification of regulatory variables for state transition of biological networks

Attractors represent steady states of biological networks. Recent studies have shown that regulatory variables can be used to steer a network state transition from an undesired attractor, such as a cancerous state, to a desired healthy one. Therefore, it is important to identify the regulatory variables and determine their time-dependent profile for state transition of a given network. However, this is a challenging task since regulatory variables have to be identified among numerous candidates in a large-scale biological network. In this study, we developed a new method for identifying regulatory variables in large-scale biological networks for the purpose of state transition. As a result, a set of optimal regulatory variables can be determined based on formulating and solving a mixed-integer nonlinear dynamic optimization problem. A relaxation scheme is used to overcome the difficulties in solving this complex problem containing a large number of binary variables. The solution to this problem simultaneously identifies the optimal regulatory variables, provides strength of regulatory interactions, and obtains the minimal control time to realize the required state transition. In addition, by adjusting the objective function, various combinations of the strength of regulatory interactions and the transition time can be achieved according to the requirement for disease therapy. Results of three case studies (a myeloid differentiation regulatory network, a cancer gene regulatory network, and a T-LGL signaling network) demonstrate the efficacy of the proposed approach. Therefore, this study establishes an appropriate framework for identifying the regulatory variables for state transition of complex biological networks.

Scalable Inference of Ordinary Differential Equation Models of Biochemical Processes

Ordinary differential equation models have become a standard tool for the mechanistic description of biochemical processes. If parameters are inferred from experimental data, such mechanistic models can provide accurate predictions about the behavior of latent variables or the process under new experimental conditions. Complementarily, inference of model structure can be used to identify the most plausible model structure from a set of candidates, and, thus, gain novel biological insight. Several toolboxes can infer model parameters and structure for small- to medium-scale mechanistic models out of the box. However, models for highly multiplexed datasets can require hundreds to thousands of state variables and parameters. For the analysis of such large-scale models, most algorithms require intractably high computation times. This chapter provides an overview of the state-of-the-art methods for parameter and model inference, with an emphasis on scalability.

Regulation of Gene Expression and Signaling Pathway Activity in Mammalian Cells by Automated Microfluidics Feedback Control

Gene networks and signaling pathways display complex topologies and, as a result, complex nonlinear behaviors. Accumulating evidence shows that both static (concentration) and dynamical (rate-of-change) features of transcription factors, ligands and environmental stimuli control downstream processes and ultimately cellular functions. Currently, however, methods to generate stimuli with the desired features to probe cell response are still lacking. Here, combining tools from Control Engineering and Synthetic Biology (cybergenetics), we propose a simple and cost-effective microfluidics-based platform to precisely regulate gene expression and signaling pathway activity in mammalian cells by means of real-time feedback control. We show that this platform allows (i) to automatically regulate gene expression from inducible promoters in different cell types, including mouse embryonic stem cells; (ii) to precisely regulate the activity of the mTOR signaling pathway in single cells; (iii) to build a biohybrid oscillator in single embryonic stem cells by interfacing biological parts with virtual in silico counterparts. Ultimately, this platform can be used to probe gene networks and signaling pathways to understand how they process static and dynamic features of specific stimuli, as well as for the rapid prototyping of synthetic circuits for biotechnology and biomedical purposes.

Misinterpretation risks of global stochastic optimisation of kinetic models revealed by multiple optimisation runs

One of use cases for metabolic network optimisation of biotechnologically applied microorganisms is the in silico design of new strains with an improved distribution of metabolic fluxes. Global stochastic optimisation methods (genetic algorithms, evolutionary programing, particle swarm and others) can optimise complicated nonlinear kinetic models and are friendly for unexperienced user: they can return optimisation results with default method settings (population size, number of generations and others) and without adaptation of the model. Drawbacks of these methods (stochastic behaviour, undefined duration of optimisation, possible stagnation and no guaranty of reaching optima) cause optimisation result misinterpretation risks considering the very diverse educational background of the systems biology and synthetic biology research community. Different methods implemented in the COPASI software package are tested in this study to determine their ability to find feasible solutions and assess the convergence speed to the best value of the objective function. Special attention is paid to the potential misinterpretation of results. Optimisation methods are tested with additional constraints that can be introduced to ensure the biological feasibility of the resulting optimised design: (1) total enzyme activity constraint (called also amino acid pool constraint) to limit the sum of enzyme concentrations and (2) homeostatic constraint limiting steady state metabolite concentration corridor around the steady state concentrations of metabolites in the original model. Impact of additional constraints on the performance of optimisation methods and misinterpretation risks is analysed.

Optimal control of bacterial growth for the maximization of metabolite production

Microorganisms have evolved complex strategies for controlling the distribution of available resources over cellular functions. Biotechnology aims at interfering with these strategies, so as to optimize the production of metabolites and other compounds of interest, by (re)engineering the underlying regulatory networks of the cell. The resulting reallocation of resources can be described by simple so-called self-replicator models and the maximization of the synthesis of a product of interest formulated as a dynamic optimal control problem. Motivated by recent experimental work, we are specifically interested in the maximization of metabolite production in cases where growth can be switched off through an external control signal. We study various optimal control problems for the corresponding self-replicator models by means of a combination of analytical and computational techniques. We show that the optimal solutions for biomass maximization and product maximization are very similar in the case of unlimited nutrient supply, but diverge when nutrients are limited. Moreover, external growth control overrides natural feedback growth control and leads to an optimal scheme consisting of a first phase of growth maximization followed by a second phase of product maximization. This two-phase scheme agrees with strategies that have been proposed in metabolic engineering. More generally, our work shows the potential of optimal control theory for better understanding and improving biotechnological production processes.

Dynamic Optimization with Particle Swarms (DOPS): a meta-heuristic for parameter estimation in biochemical models

BACKGROUND

Mathematical modeling is a powerful tool to analyze, and ultimately design biochemical networks. However, the estimation of the parameters that appear in biochemical models is a significant challenge. Parameter estimation typically involves expensive function evaluations and noisy data, making it difficult to quickly obtain optimal solutions. Further, biochemical models often have many local extrema which further complicates parameter estimation. Toward these challenges, we developed Dynamic Optimization with Particle Swarms (DOPS), a novel hybrid meta-heuristic that combined multi-swarm particle swarm optimization with dynamically dimensioned search (DDS). DOPS uses a multi-swarm particle swarm optimization technique to generate candidate solution vectors, the best of which is then greedily updated using dynamically dimensioned search.

RESULTS

We tested DOPS using classic optimization test functions, biochemical benchmark problems and real-world biochemical models. We performed [Formula: see text] = 25 trials with [Formula: see text] = 4000 function evaluations per trial, and compared the performance of DOPS with other commonly used meta-heuristics such as differential evolution (DE), simulated annealing (SA) and dynamically dimensioned search (DDS). On average, DOPS outperformed other common meta-heuristics on the optimization test functions, benchmark problems and a real-world model of the human coagulation cascade.

CONCLUSIONS

DOPS is a promising meta-heuristic approach for the estimation of biochemical model parameters in relatively few function evaluations. DOPS source code is available for download under a MIT license at http://www.varnerlab.org .

A Review of Dynamic Modeling Approaches and Their Application in Computational Strain Optimization for Metabolic Engineering

Mathematical modeling is a key process to describe the behavior of biological networks. One of the most difficult challenges is to build models that allow quantitative predictions of the cells' states along time. Recently, this issue started to be tackled through novel in silico approaches, such as the reconstruction of dynamic models, the use of phenotype prediction methods, and pathway design via efficient strain optimization algorithms. The use of dynamic models, which include detailed kinetic information of the biological systems, potentially increases the scope of the applications and the accuracy of the phenotype predictions. New efforts in metabolic engineering aim at bridging the gap between this approach and other different paradigms of mathematical modeling, as constraint-based approaches. These strategies take advantage of the best features of each method, and deal with the most remarkable limitation—the lack of available experimental information—which affects the accuracy and feasibility of solutions. Parameter estimation helps to solve this problem, but adding more computational cost to the overall process. Moreover, the existing approaches include limitations such as their scalability, flexibility, convergence time of the simulations, among others. The aim is to establish a trade-off between the size of the model and the level of accuracy of the solutions. In this work, we review the state of the art of dynamic modeling and related methods used for metabolic engineering applications, including approaches based on hybrid modeling. We describe approaches developed to undertake issues regarding the mathematical formulation and the underlying optimization algorithms, and that address the phenotype prediction by including available kinetic rate laws of metabolic processes. Then, we discuss how these have been used and combined as the basis to build computational strain optimization methods for metabolic engineering purposes, how they lead to bi-level schemes that can be used in the industry, including a consideration of their limitations.

Heparan Sulfate Induces Necroptosis in Murine Cardiomyocytes: A Medical-In silico Approach Combining In vitro Experiments and Machine Learning

Life-threatening cardiomyopathy is a severe, but common, complication associated with severe trauma or sepsis. Several signaling pathways involved in apoptosis and necroptosis are linked to trauma- or sepsis-associated cardiomyopathy. However, the underling causative factors are still debatable. Heparan sulfate (HS) fragments belong to the class of danger/damage-associated molecular patterns liberated from endothelial-bound proteoglycans by heparanase during tissue injury associated with trauma or sepsis. We hypothesized that HS induces apoptosis or necroptosis in murine cardiomyocytes. By using a novel Medical-In silico approach that combines conventional cell culture experiments with machine learning algorithms, we aimed to reduce a significant part of the expensive and time-consuming cell culture experiments and data generation by using computational intelligence (refinement and replacement). Cardiomyocytes exposed to HS showed an activation of the intrinsic apoptosis signal pathway via cytochrome C and the activation of caspase 3 (both p < 0.001). Notably, the exposure of HS resulted in the induction of necroptosis by tumor necrosis factor α and receptor interaction protein 3 (p < 0.05; p < 0.01) and, hence, an increased level of necrotic cardiomyocytes. In conclusion, using this novel Medical-In silico approach, our data suggest (i) that HS induces necroptosis in cardiomyocytes by phosphorylation (activation) of receptor-interacting protein 3, (ii) that HS is a therapeutic target in trauma- or sepsis-associated cardiomyopathy, and (iii) indicate that this proof-of-concept is a first step toward simulating the extent of activated components in the pro-apoptotic pathway induced by HS with only a small data set gained from the in vitro experiments by using machine learning algorithms.

Modelling the host-pathogen interactions of macrophages and Candida albicans using Game Theory and dynamic optimization

The release of fungal cells following macrophage phagocytosis, called non-lytic expulsion, is reported for several fungal pathogens. On one hand, non-lytic expulsion may benefit the fungus in escaping the microbicidal environment of the phagosome. On the other hand, the macrophage could profit in terms of avoiding its own lysis and being able to undergo proliferation. To analyse the causes of non-lytic expulsion and the relevance of macrophage proliferation in the macrophage–Candida albicans interaction, we employ Evolutionary Game Theory and dynamic optimization in a sequential manner. We establish a game-theoretical model describing the different strategies of the two players after phagocytosis. Depending on the parameter values, we find four different Nash equilibria and determine the influence of the systems state of the host upon the game. As our Nash equilibria are a direct consequence of the model parameterization, we can depict several biological scenarios. A parameter region, where the host response is robust against the fungal infection, is determined. We further apply dynamic optimization to analyse whether macrophage mitosis is relevant in the host–pathogen interaction of macrophages and C. albicans. For this, we study the population dynamics of the macrophage–C. albicans interactions and the corresponding optimal controls for the macrophages, indicating the best macrophage strategy of switching from proliferation to attacking fungal cells.

Mathematical Modeling Approaches in Plant Metabolomics

The experimental analysis of a plant metabolome typically results in a comprehensive and multidimensional data set. To interpret metabolomics data in the context of biochemical regulation and environmental fluctuation, various approaches of mathematical modeling have been developed and have proven useful. In this chapter, a general introduction to mathematical modeling is presented and discussed in context of plant metabolism. A particular focus is laid on the suitability of mathematical approaches to functionally integrate plant metabolomics data in a metabolic network and combine it with other biochemical or physiological parameters.

An improved hybrid of particle swarm optimization and the gravitational search algorithm to produce a kinetic parameter estimation of aspartate biochemical pathways

Mathematical modelling is fundamental to understand the dynamic behavior and regulation of the biochemical metabolisms and pathways that are found in biological systems. Pathways are used to describe complex processes that involve many parameters. It is important to have an accurate and complete set of parameters that describe the characteristics of a given model. However, measuring these parameters is typically difficult and even impossible in some cases. Furthermore, the experimental data are often incomplete and also suffer from experimental noise. These shortcomings make it challenging to identify the best-fit parameters that can represent the actual biological processes involved in biological systems. Computational approaches are required to estimate these parameters. The estimation is converted into multimodal optimization problems that require a global optimization algorithm that can avoid local solutions. These local solutions can lead to a bad fit when calibrating with a model. Although the model itself can potentially match a set of experimental data, a high-performance estimation algorithm is required to improve the quality of the solutions. This paper describes an improved hybrid of particle swarm optimization and the gravitational search algorithm (IPSOGSA) to improve the efficiency of a global optimum (the best set of kinetic parameter values) search. The findings suggest that the proposed algorithm is capable of narrowing down the search space by exploiting the feasible solution areas. Hence, the proposed algorithm is able to achieve a near-optimal set of parameters at a fast convergence speed. The proposed algorithm was tested and evaluated based on two aspartate pathways that were obtained from the BioModels Database. The results show that the proposed algorithm outperformed other standard optimization algorithms in terms of accuracy and near-optimal kinetic parameter estimation. Nevertheless, the proposed algorithm is only expected to work well in small scale systems. In addition, the results of this study can be used to estimate kinetic parameter values in the stage of model selection for different experimental conditions.

A systematic model identification method for chemical transformation pathways - the case of heroin biomarkers in wastewater

This study presents a novel statistical approach for identifying sequenced chemical transformation pathways in combination with reaction kinetics models. The proposed method relies on sound uncertainty propagation by considering parameter ranges and associated probability distribution obtained at any given transformation pathway levels as priors for parameter estimation at any subsequent transformation levels. The method was applied to calibrate a model predicting the transformation in untreated wastewater of six biomarkers, excreted following human metabolism of heroin and codeine. The method developed was compared to parameter estimation methods commonly encountered in literature (i.e., estimation of all parameters at the same time and parameter estimation with fix values for upstream parameters) by assessing the model prediction accuracy, parameter identifiability and uncertainty analysis. Results obtained suggest that the method developed has the potential to outperform conventional approaches in terms of prediction accuracy, transformation pathway identification and parameter identifiability. This method can be used in conjunction with optimal experimental designs to effectively identify model structures and parameters. This method can also offer a platform to promote a closer interaction between analytical chemists and modellers to identify models for biochemical transformation pathways, being a prominent example for the emerging field of wastewater-based epidemiology.

A parallel metaheuristic for large mixed-integer dynamic optimization problems, with applications in computational biology

Background

We consider a general class of global optimization problems dealing with nonlinear dynamic models. Although this class is relevant to many areas of science and engineering, here we are interested in applying this framework to the reverse engineering problem in computational systems biology, which yields very large mixed-integer dynamic optimization (MIDO) problems. In particular, we consider the framework of logic-based ordinary differential equations (ODEs).

Methods

We present saCeSS2, a parallel method for the solution of this class of problems. This method is based on an parallel cooperative scatter search metaheuristic, with new mechanisms of self-adaptation and specific extensions to handle large mixed-integer problems. We have paid special attention to the avoidance of convergence stagnation using adaptive cooperation strategies tailored to this class of problems.

Results

We illustrate its performance with a set of three very challenging case studies from the domain of dynamic modelling of cell signaling. The simpler case study considers a synthetic signaling pathway and has 84 continuous and 34 binary decision variables. A second case study considers the dynamic modeling of signaling in liver cancer using high-throughput data, and has 135 continuous and 109 binaries decision variables. The third case study is an extremely difficult problem related with breast cancer, involving 690 continuous and 138 binary decision variables. We report computational results obtained in different infrastructures, including a local cluster, a large supercomputer and a public cloud platform. Interestingly, the results show how the cooperation of individual parallel searches modifies the systemic properties of the sequential algorithm, achieving superlinear speedups compared to an individual search (e.g. speedups of 15 with 10 cores), and significantly improving (above a 60%) the performance with respect to a non-cooperative parallel scheme. The scalability of the method is also good (tests were performed using up to 300 cores).

Conclusions

These results demonstrate that saCeSS2 can be used to successfully reverse engineer large dynamic models of complex biological pathways. Further, these results open up new possibilities for other MIDO-based large-scale applications in the life sciences such as metabolic engineering, synthetic biology, drug scheduling.

Identification of optimal strategies for state transition of complex biological networks

Complex biological networks typically contain numerous parameters, and determining feasible strategies for state transition by parameter perturbation is not a trivial task. In the present study, based on dynamical and structural analyses of the biological network, we optimized strategies for controlling variables in a two-node gene regulatory network and a T-cell large granular lymphocyte signaling network associated with blood cancer by using an efficient dynamic optimization method. Optimization revealed the critical value for each decision variable to steer the system from an undesired state into a desired attractor. In addition, the minimum time for the state transition was determined by defining and solving a time-optimal control problem. Moreover, time-dependent variable profiles for state transitions were achieved rather than constant values commonly adopted in previous studies. Furthermore, the optimization method allows multiple controls to be simultaneously adjusted to drive the system out of an undesired attractor. Optimization improved the results of the parameter perturbation method, thus providing a valuable guidance for experimental design.

Search for a Minimal Set of Parameters by Assessing the Total Optimization Potential for a Dynamic Model of a Biochemical Network

Selecting an efficient small set of adjustable parameters to improve metabolic features of an organism is important for a reduction of implementation costs and risks of unpredicted side effects. In practice, to avoid the analysis of a huge combinatorial space for the possible sets of adjustable parameters, experience-, and intuition-based subsets of parameters are often chosen, possibly leaving some interesting counter-intuitive combinations of parameters unrevealed. The combinatorial scan of possible adjustable parameter combinations at the model optimization level is possible; however, the number of analyzed combinations is still limited. The total optimization potential (TOP) approach is proposed to assess the full potential for increasing the value of the objective function by optimizing all possible adjustable parameters. This seemingly unpractical combination of adjustable parameters allows assessing the maximum attainable value of the objective function and stopping the combinatorial space scanning when the desired fraction of TOP is reached and any further increase in the number of adjustable parameters cannot bring any reasonable improvement. The relation between the number of adjustable parameters and the reachable fraction of TOP is a valuable guideline in choosing a rational solution for industrial implementation. The TOP approach is demonstrated on the basis of two case studies.

Optimality principles reveal a complex interplay of intermediate toxicity and kinetic efficiency in the regulation of prokaryotic metabolism

A precise and rapid adjustment of fluxes through metabolic pathways is crucial for organisms to prevail in changing environmental conditions. Based on this reasoning, many guiding principles that govern the evolution of metabolic networks and their regulation have been uncovered. To this end, methods from dynamic optimization are ideally suited since they allow to uncover optimality principles behind the regulation of metabolic networks. We used dynamic optimization to investigate the influence of toxic intermediates in connection with the efficiency of enzymes on the regulation of a linear metabolic pathway. Our results predict that transcriptional regulation favors the control of highly efficient enzymes with less toxic upstream intermediates to reduce accumulation of toxic downstream intermediates. We show that the derived optimality principles hold by the analysis of the interplay between intermediate toxicity and pathway regulation in the metabolic pathways of over 5000 sequenced prokaryotes. Moreover, using the lipopolysaccharide biosynthesis in Escherichia coli as an example, we show how knowledge about the relation of regulation, kinetic efficiency and intermediate toxicity can be used to identify drug targets, which control endogenous toxic metabolites and prevent microbial growth. Beyond prokaryotes, we discuss the potential of our findings for the development of antifungal drugs.

Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks

Mechanistic mathematical modeling of biochemical reaction networks using ordinary differential equation (ODE) models has improved our understanding of small- and medium-scale biological processes. While the same should in principle hold for large- and genome-scale processes, the computational methods for the analysis of ODE models which describe hundreds or thousands of biochemical species and reactions are missing so far. While individual simulations are feasible, the inference of the model parameters from experimental data is computationally too intensive. In this manuscript, we evaluate adjoint sensitivity analysis for parameter estimation in large scale biochemical reaction networks. We present the approach for time-discrete measurement and compare it to state-of-the-art methods used in systems and computational biology. Our comparison reveals a significantly improved computational efficiency and a superior scalability of adjoint sensitivity analysis. The computational complexity is effectively independent of the number of parameters, enabling the analysis of large- and genome-scale models. Our study of a comprehensive kinetic model of ErbB signaling shows that parameter estimation using adjoint sensitivity analysis requires a fraction of the computation time of established methods. The proposed method will facilitate mechanistic modeling of genome-scale cellular processes, as required in the age of omics.

In this manuscript, we introduce a scalable method for parameter estimation for genome-scale biochemical reaction networks. Mechanistic models for genome-scale biochemical reaction networks describe the behavior of thousands of chemical species using thousands of parameters. Standard methods for parameter estimation are usually computationally intractable at these scales. Adjoint sensitivity based approaches have been suggested to have superior scalability but any rigorous evaluation is lacking. We implement a toolbox for adjoint sensitivity analysis for biochemical reaction network which also supports the import of SBML models. We show by means of a set of benchmark models that adjoint sensitivity based approaches unequivocally outperform standard approaches for large-scale models and that the achieved speedup increases with respect to both the number of parameters and the number of chemical species in the model. This demonstrates the applicability of adjoint sensitivity based approaches to parameter estimation for genome-scale mechanistic model. The MATLAB toolbox implementing the developed methods is available from http://ICB-DCM.github.io/AMICI/.

Unearthing the transition rates between photoreceptor conformers

Background

Obtaining accurate estimates of biological or enzymatic reaction rates is critical in understanding the design principles of a network and how biological processes can be experimentally manipulated on demand. In many cases experimental limitations mean that some enzymatic rates cannot be measured directly, requiring mathematical algorithms to estimate them. Here, we describe a methodology that calculates rates at which light-regulated proteins switch between conformational states. We focus our analysis on the phytochrome family of photoreceptors found in cyanobacteria, plants and many optogenetic tools. Phytochrome proteins change between active (P_A) and inactive (P_I) states at rates that are proportional to photoconversion cross-sections and influenced by light quality, light intensity, thermal reactions and dimerisation. This work presents a method that can accurately calculate these photoconversion cross-sections in the presence of multiple non-light regulated reactions.

Results

Our approach to calculating the photoconversion cross-sections comprises three steps: i) calculate the thermal reversion reaction rate(s); ii) develop search spaces from which all possible sets of photoconversion cross-sections exist, and; iii) estimate extinction coefficients that describe our absorption spectra. We confirm that the presented approach yields accurate results through the use of simulated test cases. Our test cases were further expanded to more realistic scenarios where noise, multiple thermal reactions and dimerisation are considered. Finally, we present the photoconversion cross-sections of an Arabidopsis phyB N-terminal fragment commonly used in optogenetic tools.

Conclusions

The calculation of photoconversion cross-sections has implications for both photoreceptor and synthetic biologists. Our method allows, for the first time, direct comparisons of photoconversion cross-sections and response speeds of photoreceptors in different cellular environments and synthetic tools. Due to the generality of our procedure, as shown by the application to multiple test cases, the photoconversion cross-sections and quantum yields of any photoreceptor might now, in principle, be obtained.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0368-y) contains supplementary material, which is available to authorized users.

Gene Regulatory Network Inferences Using a Maximum-Relevance and Maximum-Significance Strategy

Recovering gene regulatory networks from expression data is a challenging problem in systems biology that provides valuable information on the regulatory mechanisms of cells. A number of algorithms based on computational models are currently used to recover network topology. However, most of these algorithms have limitations. For example, many models tend to be complicated because of the "large p, small n" problem. In this paper, we propose a novel regulatory network inference method called the maximum-relevance and maximum-significance network (MRMSn) method, which converts the problem of recovering networks into a problem of how to select the regulator genes for each gene. To solve the latter problem, we present an algorithm that is based on information theory and selects the regulator genes for a specific gene by maximizing the relevance and significance. A first-order incremental search algorithm is used to search for regulator genes. Eventually, a strict constraint is adopted to adjust all of the regulatory relationships according to the obtained regulator genes and thus obtain the complete network structure. We performed our method on five different datasets and compared our method to five state-of-the-art methods for network inference based on information theory. The results confirm the effectiveness of our method.