@article{BennHiepenOsterlandetal.2017, author = {Benn, Andreas and Hiepen, Christian and Osterland, Marc and Sch{\"u}tte, Christof and Zwijsen, An and Knaus, Petra}, title = {Role of bone morphogenetic proteins in sprouting angiogenesis: differential BMP receptor-dependent signaling pathways balance stalk vs. tip cell competence}, volume = {31}, journal = {FASEB Journal}, number = {11}, doi = {10.1096/fj.201700193RR}, pages = {4720 -- 4733}, year = {2017}, abstract = {Before the onset of sprouting angiogenesis, the endothelium is prepatterned for the positioning of tip and stalk cells. Both cell identities are not static, as endothelial cells (ECs) constantly compete for the tip cell position in a dynamic fashion. Here, we show that both bone morphogenetic protein (BMP) 2 and BMP6 are proangiogenic in vitro and ex vivo and that the BMP type I receptors, activin receptor-like kinase (ALK)3 and ALK2, play crucial and distinct roles in this process. BMP2 activates the expression of tip cell-associated genes, such as DLL4 (delta-like ligand 4) and KDR (kinase insert domain receptor), and p38-heat shock protein 27 (HSP27)-dependent cell migration, thereby generating tip cell competence. Whereas BMP6 also triggers collective cell migration via the p38-HSP27 signaling axis, BMP6 induces in addition SMAD1/5 signaling, thereby promoting the expression of stalk cell-associated genes, such as HES1 (hairy and enhancer of split 1) and FLT1 (fms-like tyrosine kinase 1). Specifically, ALK3 is required for sprouting from HUVEC spheroids, whereas ALK2 represses sprout formation. We demonstrate that expression levels and respective complex formation of BMP type I receptors in ECs determine stalk vs. tip cell identity, thus contributing to endothelial plasticity during sprouting angiogenesis. As antiangiogenic monotherapies that target the VEGF or ALK1 pathways have not fulfilled efficacy objectives in clinical trials, the selective targeting of the ALK2/3 pathways may be an attractive new approach.}, language = {en} } @misc{WinkelmannSchuette2017, author = {Winkelmann, Stefanie and Sch{\"u}tte, Christof}, title = {Hybrid Models for Chemical Reaction Networks: Multiscale Theory and Application to Gene Regulatory Systems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-64264}, year = {2017}, abstract = {Well-mixed stochastic chemical kinetics are properly modelled by the chemical master equation (CME) and associated Markov jump processes in molecule number space. If the reactants are present in large amounts, however, corresponding simulations of the stochastic dynamics become computationally expensive and model reductions are demanded. The classical model reduction approach uniformly rescales the overall dynamics to obtain deterministic systems characterized by ordinary differential equations, the well-known mass action reaction rate equations. For systems with multiple scales there exist hybrid approaches that keep parts of the system discrete while another part is approximated either using Langevin dynamics or deterministically. This paper aims at giving a coherent overview of the different hybrid approaches, focusing on their basic concepts and the relation between them. We derive a novel general description of such hybrid models that allows to express various forms by one type of equation. We also check in how far the approaches apply to model extensions of the CME for dynamics which do not comply with the central well-mixed condition and require some spatial resolution. A simple but meaningful gene expression system with negative self-regulation is analysed to illustrate the different approximation qualities of some of the hybrid approaches discussed.}, language = {en} } @misc{BittracherKoltaiKlusetal.2017, author = {Bittracher, Andreas and Koltai, P{\´e}ter and Klus, Stefan and Banisch, Ralf and Dellnitz, Michael and Sch{\"u}tte, Christof}, title = {Transition manifolds of complex metastable systems: Theory and data-driven computation of effective dynamics}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-63822}, year = {2017}, abstract = {We consider complex dynamical systems showing metastable behavior but no local separation of fast and slow time scales. The article raises the question of whether such systems exhibit a low-dimensional manifold supporting its effective dynamics. For answering this question, we aim at finding nonlinear coordinates, called reaction coordinates, such that the projection of the dynamics onto these coordinates preserves the dominant time scales of the dynamics. We show that, based on a specific reducibility property, the existence of good low-dimensional reaction coordinates preserving the dominant time scales is guaranteed. Based on this theoretical framework, we develop and test a novel numerical approach for computing good reaction coordinates. The proposed algorithmic approach is fully local and thus not prone to the curse of dimension with respect to the state space of the dynamics. Hence, it is a promising method for data-based model reduction of complex dynamical systems such as molecular dynamics.}, language = {en} } @misc{ZhangSchuette2017, author = {Zhang, Wei and Sch{\"u}tte, Christof}, title = {Reliable approximation of long relaxation timescales in molecular dynamics}, issn = {1438-0064}, doi = {10.3390/e19070367}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-63718}, year = {2017}, abstract = {Many interesting rare events in molecular systems like ligand association, protein folding or con- formational changes happen on timescales that often are not accessible by direct numerical simulation. Therefore rare event approximation approaches like interface sampling, Markov state model building or advanced reaction coordinate based free energy estimation have attracted huge attention recently. In this article we analyze the reliability of such approaches: How precise is an estimate of long relaxation timescales of molecular systems resulting from various forms of rare event approximation methods? Our results give a theoretical answer to this question by relating it with the transfer operator approach to molecular dynamics. By doing so they also allow for understanding deep connections between the different approaches.}, language = {en} } @article{KoltaiCiccottiSchuette2016, author = {Koltai, Peter and Ciccotti, Giovanni and Sch{\"u}tte, Christof}, title = {On Markov state models for non-equilibrium molecular dynamics}, volume = {145}, journal = {The Journal of Chemical Physics}, number = {174103}, doi = {10.1063/1.4966157}, year = {2016}, language = {en} } @article{KlusNueskeKoltaietal.2018, author = {Klus, Stefan and N{\"u}ske, Feliks and Koltai, Peter and Wu, Hao and Kevrekidis, Ioannis and Sch{\"u}tte, Christof and No{\´e}, Frank}, title = {Data-driven model reduction and transfer operator approximation}, volume = {28}, journal = {Journal of Nonlinear Science}, number = {3}, doi = {10.1007/s00332-017-9437-7}, pages = {985 -- 1010}, year = {2018}, language = {en} } @misc{GuptaGramatkeEinspanieretal.2017, author = {Gupta, Pooja and Gramatke, Annika and Einspanier, Ralf and Sch{\"u}tte, Christof and von Kleist, Max and Sharbati, Jutta}, title = {In silicio cytotoxicity assessment on cultured rat intestinal cells deduced from cellular impedance measurements}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-62666}, year = {2017}, abstract = {Early and reliable identification of chemical toxicity is of utmost importance. At the same time, reduction of animal testing is paramount. Therefore, methods that improve the interpretability and usability of in vitro assays are essential. xCELLigence's real-time cell analyzer (RTCA) provides a novel, fast and cost effective in vitro method to probe compound toxicity. We developed a simple mathematical framework for the qualitative and quantitative assessment of toxicity for RTCA measurements. Compound toxicity, in terms of its 50\% inhibitory concentration IC_{50} on cell growth, and parameters related to cell turnover were estimated on cultured IEC-6 cells exposed to 10 chemicals at varying concentrations. Our method estimated IC50 values of 113.05, 7.16, 28.69 and 725.15 μM for the apparently toxic compounds 2-acetylamino-fluorene, aflatoxin B1, benzo-[a]-pyrene and chloramphenicol in the tested cell line, in agreement with literature knowledge. IC_{50} values of all apparent in vivo non-toxic compounds were estimated to be non-toxic by our method. Corresponding estimates from RTCA's in-built model gave false positive (toxicity) predictions in 5/10 cases. Taken together, our proposed method reduces false positive predictions and reliably identifies chemical toxicity based on impedance measurements. The source code for the developed method including instructions is available at https://git.zib.de/bzfgupta/toxfit/tree/master.}, language = {en} } @article{HuttaryGoubergritsSchuetteetal.2017, author = {Huttary, Rudolf and Goubergrits, Leonid and Sch{\"u}tte, Christof and Bernhard, Stefan}, title = {Simulation, Identification and Statistical Variation in Cardiovascular Analysis (SISCA) - a Software Framework for Multi-compartment Lumped Modeling}, volume = {87}, journal = {Computers in Biology and Medicine}, doi = {10.1016/j.compbiomed.2017.05.021}, pages = {104 -- 123}, year = {2017}, language = {en} } @article{vonKleistSchuetteZhang2018, author = {von Kleist, Max and Sch{\"u}tte, Christof and Zhang, Wei}, title = {Statistical analysis of the first passage path ensemble of jump processes}, volume = {170}, journal = {Journal of Statistical Physics}, doi = {10.1007/s10955-017-1949-x}, pages = {809 -- 843}, year = {2018}, abstract = {The transition mechanism of jump processes between two different subsets in state space reveals important dynamical information of the processes and therefore has attracted considerable attention in the past years. In this paper, we study the first passage path ensemble of both discrete-time and continuous-time jump processes on a finite state space. The main approach is to divide each first passage path into nonreactive and reactive segments and to study them separately. The analysis can be applied to jump processes which are non-ergodic, as well as continuous-time jump processes where the waiting time distributions are non-exponential. In the particular case that the jump processes are both Markovian and ergodic, our analysis elucidates the relations between the study of the first passage paths and the study of the transition paths in transition path theory. We provide algorithms to numerically compute statistics of the first passage path ensemble. The computational complexity of these algorithms scales with the complexity of solving a linear system, for which efficient methods are available. Several examples demonstrate the wide applicability of the derived results across research areas.}, language = {en} } @article{WeberFackeldeySchuette2017, author = {Weber, Marcus and Fackeldey, Konstantin and Sch{\"u}tte, Christof}, title = {Set-Free Markov State Model Building}, volume = {146}, journal = {Journal of Chemical Physics}, number = {12}, doi = {10.1063/1.4978501}, year = {2017}, language = {en} }