Lorenz Adlung

Data2Dynamics: a modeling environment tailored to parameter estimation in dynamical systems

UNLABELLED Modeling of dynamical systems using ordinary differential equations is a popular approach in the field of systems biology. Two of the most critical steps in this approach are to construct dynamical models of biochemical reaction networks for large datasets and complex experimental conditions and to perform efficient and reliable parameter estimation for model fitting. We present a modeling environment for MATLAB that pioneers these challenges. The numerically expensive parts of the calculations such as the solving of the differential equations and of the associated sensitivity system are parallelized and automatically compiled into efficient C code. A variety of parameter estimation algorithms as well as frequentist and Bayesian methods for uncertainty analysis have been implemented and used on a range of applications that lead to publications. AVAILABILITY AND IMPLEMENTATION The Data2Dynamics modeling environment is MATLAB based, open source and freely available at http://www.data2dynamics.org. CONTACT andreas.raue@fdm.uni-freiburg.de SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Creating functional engineered variants of the single-module non-ribosomal peptide synthetase IndC by T domain exchange

Non-ribosomal peptide synthetases (NRPSs) are enzymes that catalyze ribosome-independent production of small peptides, most of which are bioactive. NRPSs act as peptide assembly lines where individual, often interconnected modules each incorporate a specific amino acid into the nascent chain. The modules themselves consist of several domains that function in the activation, modification and condensation of the substrate. NRPSs are evidently modular, yet experimental proof of the ability to engineer desired permutations of domains and modules is still sought. Here, we use a synthetic-biology approach to create a small library of engineered NRPSs, in which the domain responsible for carrying the activated amino acid (T domain) is exchanged with natural or synthetic T domains. As a model system, we employ the single-module NRPS IndC from Photorhabdus luminescens that produces the blue pigment indigoidine. As chassis we use Escherichia coli. We demonstrate that heterologous T domain exchange is possible, even for T domains derived from different organisms. Interestingly, substitution of the native T domain with a synthetic one enhanced indigoidine production. Moreover, we show that selection of appropriate inter-domain linker regions is critical for functionality. Taken together, our results extend the engineering avenues for NRPSs, as they point out the possibility of combining domain sequences coming from different pathways, organisms or from conservation criteria. Moreover, our data suggest that NRPSs can be rationally engineered to control the level of production of the corresponding peptides. This could have important implications for industrial and medical applications.

Identification of isoform-specific dynamics in phosphorylation-dependent STAT5 dimerization by quantitative mass spectrometry and mathematical modeling

STAT5A and STAT5B are important transcription factors that dimerize and transduce activation signals of cytokine receptors directly to the nucleus. A typical cytokine that mediates STAT5 activation is erythropoietin (Epo). Differential functions of STAT5A and STAT5B have been reported. However, the extent to which phosphorylated STAT5A and STAT5B (pSTAT5A, pSTAT5B) form homo- or heterodimers is not understood, nor is how this might influence the signal transmission to the nucleus. To study this, we designed a concept to investigate the isoform-specific dimerization behavior of pSTAT5A and pSTAT5B that comprises isoform-specific immunoprecipitation (IP), measurement of the degree of phosphorylation, and isoform ratio determination between STAT5A and STAT5B. For the main analytical method, we employed quantitative label-free and -based mass spectrometry. For the cellular model system, we used Epo receptor (EpoR)-expressing BaF3 cells (BaF3-EpoR) stimulated with Epo. Three hypotheses of dimer formation between pSTAT5A and pSTAT5B were used to explain the analytical results by a static mathematical model: formation of (i) homodimers only, (ii) heterodimers only, and (iii) random formation of homo- and heterodimers. The best agreement between experimental data and model simulations was found for the last case. Dynamics of cytoplasmic STAT5 dimerization could be explained by distinct nuclear import rates and individual nuclear retention for homo- and heterodimers of phosphorylated STAT5.

The quorum-sensing regulator ComA from Bacillus subtilis activates transcription using topologically distinct DNA motifs

ComA-like transcription factors regulate the quorum response in numerous Gram-positive bacteria. ComA proteins belong to the tetrahelical helix-turn-helix superfamily of transcriptional activators, which bind as homodimers to inverted sequence repeats in the DNA. Here, we report that ComA from Bacillus subtilis recognizes a topologically distinct motif, in which the binding elements form a direct repeat. We provide in vitro and in vivo evidence that the canonical and non-canonical site play an important role in facilitating type I and type II promoter activation, respectively, by interacting with different subunits of RNA polymerase. We furthermore show that there is a variety of contexts in which the non-canonical site can occur and identify new direct target genes that are located within the integrative and conjugative element ICEBs1. We therefore suggest that ComA acts as a multifunctional transcriptional activator and provides a striking example for complexity in protein–DNA interactions that evolved in the context of quorum sensing.

Protein abundance of AKT and ERK pathway components governs cell type-specific regulation of proliferation

Signaling through the AKT and ERK pathways controls cell proliferation. However, the integrated regulation of this multistep process, involving signal processing, cell growth and cell cycle progression, is poorly understood. Here, we study different hematopoietic cell types, in which AKT and ERK signaling is triggered by erythropoietin (Epo). Although these cell types share the molecular network topology for pro‐proliferative Epo signaling, they exhibit distinct proliferative responses. Iterating quantitative experiments and mathematical modeling, we identify two molecular sources for cell type‐specific proliferation. First, cell type‐specific protein abundance patterns cause differential signal flow along the AKT and ERK pathways. Second, downstream regulators of both pathways have differential effects on proliferation, suggesting that protein synthesis is rate‐limiting for faster cycling cells while slower cell cycles are controlled at the G1‐S progression. The integrated mathematical model of Epo‐driven proliferation explains cell type‐specific effects of targeted AKT and ERK inhibitors and faithfully predicts, based on the protein abundance, anti‐proliferative effects of inhibitors in primary human erythroid progenitor cells. Our findings suggest that the effectiveness of targeted cancer therapy might become predictable from protein abundance.