Characterizing the pharmacological interaction of the antimalarial combination artefenomel-piperaquine in healthy volunteers with induced blood-stage Plasmodium falciparum to predict efficacy in patients with malaria.

Background: The combination antimalarial artefenomel-piperaquine failed to achieve target efficacy in a phase 2b study in Africa and Vietnam. We retrospectively evaluated whether characterizing the pharmacological interaction of this antimalarial combination in a volunteer infection study (VIS) would have enabled prediction of the phase 2b study results.

Methods: Twenty-four healthy adults enrolled over three consecutive cohorts were inoculated with Plasmodium falciparum-infected erythrocytes on day 0.

BlotIt-Optimal alignment of Western blot and qPCR experiments.

Biological systems are frequently analyzed by means of mechanistic mathematical models. In order to infer model parameters and provide a useful model that can be employed for systems understanding and hypothesis testing, the model is often calibrated on quantitative, time-resolved data. To do so, it is typically important to compare experimental measurements over broad time ranges and various experimental conditions, e.

Dynamic modeling of Nrf2 pathway activation in liver cells after toxicant exposure.

Cells are exposed to oxidative stress and reactive metabolites every day. The Nrf2 signaling pathway responds to oxidative stress by upregulation of antioxidants like glutathione (GSH) to compensate the stress insult and re-establish homeostasis. Although mechanisms describing the interaction between the key pathway constituents Nrf2, Keap1 and p62 are widely reviewed and discussed in literature, quantitative dynamic models bringing together these mechanisms with time-resolved data are limited.

Efficient simulation of clinical target response surfaces.

Simulation of combination therapies is challenging due to computational complexity. Either a simple model is used to simulate the response for many combinations of concentration to generate a response surface but parameter variability and uncertainty are neglected and the concentrations are constant-the link to the doses to be administered is difficult to make-or a population pharmacokinetic/pharmacodynamic model is used to predict the response to combination therapy in a clinical trial taking into account the time-varying concentration profile, interindividual variability (IIV), and parameter uncertainty but simulations are limited to only a few selected doses. We devised new algorithms to efficiently search for the combination doses that achieve a predefined efficacy target while taking into account the IIV and parameter uncertainty.

Mathematical modeling of drug-induced receptor internalization in the HER2-positive SKBR3 breast cancer cell-line.

About 20% of breast cancer tumors over-express the HER2 receptor. Trastuzumab, an approved drug to treat this type of breast cancer, is a monoclonal antibody directly binding at the HER2 receptor and ultimately inhibiting cancer cell growth. The goal of our study was to understand the early impact of trastuzumab on HER2 internalization and recycling in the HER2-overexpressing breast cancer cell line SKBR3.

Local Riemannian geometry of model manifolds and its implications for practical parameter identifiability.

When non-linear models are fitted to experimental data, parameter estimates can be poorly constrained albeit being identifiable in principle. This means that along certain paths in parameter space, the log-likelihood does not exceed a given statistical threshold but remains bounded. This situation, denoted as practical non-identifiability, can be detected by Monte Carlo sampling or by systematic scanning using the profile likelihood method.

Model-based identification of TNFα-induced IKKβ-mediated and IκBα-mediated regulation of NFκB signal transduction as a tool to quantify the impact of drug-induced liver injury compounds.

Drug-induced liver injury (DILI) has become a major problem for patients and for clinicians, academics and the pharmaceutical industry. To date, existing hepatotoxicity test systems are only poorly predictive and the underlying mechanisms are still unclear. One of the factors known to amplify hepatotoxicity is the tumor necrosis factor alpha (TNFα), especially due to its synergy with commonly used drugs such as diclofenac.

A Dynamic Mathematical Model of Bile Acid Clearance in HepaRG Cells.

A dynamic model based on ordinary differential equations that describes uptake, basolateral and canalicular export of taurocholic acid (TCA) in human HepaRG cells is presented. The highly reproducible inter-assay experimental data were used to reliably estimate model parameters. Primary human hepatocytes were similarly evaluated to establish a mathematical model, but with notably higher inter-assay differences in TCA clearance and bile canaliculi dynamics.

Driving the Model to Its Limit: Profile Likelihood Based Model Reduction.

In systems biology, one of the major tasks is to tailor model complexity to information content of the data. A useful model should describe the data and produce well-determined parameter estimates and predictions. Too small of a model will not be able to describe the data whereas a model which is too large tends to overfit measurement errors and does not provide precise predictions.

Customized Steady-State Constraints for Parameter Estimation in Non-Linear Ordinary Differential Equation Models.

Ordinary differential equation models have become a wide-spread approach to analyze dynamical systems and understand underlying mechanisms. Model parameters are often unknown and have to be estimated from experimental data, e.g.

Higher-order Lie symmetries in identifiability and predictability analysis of dynamic models.

Parameter estimation in ordinary differential equations (ODEs) has manifold applications not only in physics but also in the life sciences. When estimating the ODE parameters from experimentally observed data, the modeler is frequently concerned with the question of parameter identifiability. The source of parameter nonidentifiability is tightly related to Lie group symmetries.

RPPanalyzer toolbox: an improved R package for analysis of reverse phase protein array data.

Analysis of large-scale proteomic data sets requires specialized software tools, tailored toward the requirements of individual approaches. Here we introduce an extension of an open-source software solution for analyzing reverse phase protein array (RPPA) data. The R package RPPanalyzer was designed for data preprocessing followed by basic statistical analyses and proteomic data visualization.

Pre-clustering of the B cell antigen receptor demonstrated by mathematically extended electron microscopy.

The B cell antigen receptor (BCR) plays a crucial role in adaptive immunity, since antigen-induced signaling by the BCR leads to the activation of the B cell and production of antibodies during an immune response. However, the spatial nano-scale organization of the BCR on the cell surface prior to antigen encounter is still controversial. Here, we fixed murine B cells, stained the BCRs on the cell surface with immuno-gold and visualized the distribution of the gold particles by transmission electron microscopy.

Lessons learned from quantitative dynamical modeling in systems biology.

Due to the high complexity of biological data it is difficult to disentangle cellular processes relying only on intuitive interpretation of measurements. A Systems Biology approach that combines quantitative experimental data with dynamic mathematical modeling promises to yield deeper insights into these processes. Nevertheless, with growing complexity and increasing amount of quantitative experimental data, building realistic and reliable mathematical models can become a challenging task: the quality of experimental data has to be assessed objectively, unknown model parameters need to be estimated from the experimental data, and numerical calculations need to be precise and efficient.

Profile likelihood in systems biology.

Inferring knowledge about biological processes by a mathematical description is a major characteristic of Systems Biology. To understand and predict system's behavior the available experimental information is translated into a mathematical model. Since the availability of experimental data is often limited and measurements contain noise, it is essential to appropriately translate experimental uncertainty to model parameters as well as to model predictions.

Semi-automatic determination of cell surface areas used in systems biology.

Quantitative biology requires high precision measurement of cellular parameters such as surface areas or volumes. Here, we have developed an integrated approach in which the data from 3D confocal microscopy and 2D high-resolution transmission electron microscopy were combined. The volumes and diameters of the cells within one population were automatically measured from the confocal data sets.

Heterogeneous kinetics of AKT signaling in individual cells are accounted for by variable protein concentration.

In most solid cancers, cells harboring oncogenic mutations represent only a sub-fraction of the entire population. Within this sub-fraction the expression level of mutated proteins can vary significantly due to cellular variability limiting the efficiency of targeted therapy. To address the causes of the heterogeneity, we performed a systematic analysis of one of the most frequently mutated pathways in cancer cells, the phosphatidylinositol 3 kinase (PI3K) signaling pathway.

A variational approach to parameter estimation in ordinary differential equations.

Background: Ordinary differential equations are widely-used in the field of systems biology and chemical engineering to model chemical reaction networks. Numerous techniques have been developed to estimate parameters like rate constants, initial conditions or steady state concentrations from time-resolved data. In contrast to this countable set of parameters, the estimation of entire courses of network components corresponds to an innumerable set of parameters.

Division of labor by dual feedback regulators controls JAK2/STAT5 signaling over broad ligand range.

Cellular signal transduction is governed by multiple feedback mechanisms to elicit robust cellular decisions. The specific contributions of individual feedback regulators, however, remain unclear. Based on extensive time-resolved data sets in primary erythroid progenitor cells, we established a dynamic pathway model to dissect the roles of the two transcriptional negative feedback regulators of the suppressor of cytokine signaling (SOCS) family, CIS and SOCS3, in JAK2/STAT5 signaling.

Model-based extension of high-throughput to high-content data.

Background: High-quality quantitative data is a major limitation in systems biology. The experimental data used in systems biology can be assigned to one of the following categories: assays yielding average data of a cell population, high-content single cell measurements and high-throughput techniques generating single cell data for large cell populations. For modeling purposes, a combination of data from different categories is highly desirable in order to increase the number of observable species and processes and thereby maximize the identifiability of parameters.