TY - GEN A1 - Winkelmann, Stefanie A1 - Schütte, Christof A1 - von Kleist, Max T1 - Markov Control Processes with Rare State Observation: Theory and Application to Treatment Scheduling in HIV-1 N2 - Markov Decision Processes (MDP) or Partially Observable MDPs (POMDP) are used for modelling situations in which the evolution of a process is partly random and partly controllable. These MDP theories allow for computing the optimal control policy for processes that can continuously or frequently be observed, even if only partially. However, they cannot be applied if state observation is very costly and therefore rare (in time). We present a novel MDP theory for rare, costly observations and derive the corresponding Bellman equation. In the new theory, state information can be derived for a particular cost after certain, rather long time intervals. The resulting information costs enter into the total cost and thus into the optimization criterion. This approach applies to many real world problems, particularly in the medical context, where the medical condition is examined rather rarely because examination costs are high. At the same time, the approach allows for efficient numerical realization. We demonstrate the usefulness of the novel theory by determining, from the national economic perspective, optimal therapeutic policies for the treatment of the human immunodefficiency virus (HIV) in resource-rich and resource-poor settings. Based on the developed theory and models, we discover that available drugs may not be utilized efficiently in resource-poor settings due to exorbitant diagnostic costs. T3 - ZIB-Report - 13-34 KW - information costs KW - hidden state KW - bellmann equation KW - optimal therapeutic policies KW - diagnostic frequency KW - resource-poor KW - resource-rich Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-41955 SN - 1438-0064 ER - TY - JOUR A1 - Gelß, Patrick A1 - Klus, Stefan A1 - Eisert, Jens A1 - Schütte, Christof T1 - Multidimensional Approximation of Nonlinear Dynamical Systems JF - Journal of Computational and Nonlinear Dynamics N2 - A key task in the field of modeling and analyzing nonlinear dynamical systems is the recovery of unknown governing equations from measurement data only. There is a wide range of application areas for this important instance of system identification, ranging from industrial engineering and acoustic signal processing to stock market models. In order to find appropriate representations of underlying dynamical systems, various data-driven methods have been proposed by different communities. However, if the given data sets are high-dimensional, then these methods typically suffer from the curse of dimensionality. To significantly reduce the computational costs and storage consumption, we propose the method multidimensional approximation of nonlinear dynamical systems (MANDy) which combines data-driven methods with tensor network decompositions. The efficiency of the introduced approach will be illustrated with the aid of several high-dimensional nonlinear dynamical systems. Y1 - 2019 U6 - https://doi.org/10.1115/1.4043148 VL - 14 IS - 6 ER - TY - JOUR A1 - Koltai, Péter A1 - Schütte, Christof T1 - A multiscale perturbation expansion approach for Markov state modeling of non-stationary molecular dynamics JF - SIAM J. Multiscale Model. Simul. N2 - We investigate metastable dynamical systems subject to non-stationary forcing as they appear in molecular dynamics for systems driven by external fields. We show, that if the strength of the forcing is inversely proportional to the length of the slow metastable time scales of the unforced system, then the effective behavior of the forced system on slow time scales can be described by a low-dimensional reduced master equation. Our construction is explicit and uses the multiscale perturbation expansion method called two-timing, or method of multiple scales. The reduced master equation—a Markov state model—can be assembled by constructing two equilibrium Markov state models; one for the unforced system, and one for a slightly perturbed one. Y1 - 2018 U6 - https://doi.org/10.1137/17M1146403 VL - 16 IS - 4 SP - 1455 EP - 1485 PB - SIAM ER - TY - JOUR A1 - Klus, Stefan A1 - Bittracher, Andreas A1 - Schuster, Ingmar A1 - Schütte, Christof T1 - A kernel-based approach to molecular conformation analysis JF - Journal of Chemical Physics N2 - We present a novel machine learning approach to understanding conformation dynamics of biomolecules. The approach combines kernel-based techniques that are popular in the machine learning community with transfer operator theory for analyzing dynamical systems in order to identify conformation dynamics based on molecular dynamics simulation data. We show that many of the prominent methods like Markov State Models, EDMD, and TICA can be regarded as special cases of this approach and that new efficient algorithms can be constructed based on this derivation. The results of these new powerful methods will be illustrated with several examples, in particular the alanine dipeptide and the protein NTL9. Y1 - 2018 U6 - https://doi.org/10.1063/1.5063533 VL - 149 IS - 24 ER - TY - JOUR A1 - Shao, Borong A1 - Bjaanaes, Maria A1 - Helland, Aslaug A1 - Schütte, Christof A1 - Conrad, Tim T1 - EMT network-based feature selection improves prognosis prediction in lung adenocarcinoma JF - PLOS ONE N2 - Various feature selection algorithms have been proposed to identify cancer prognostic biomarkers. In recent years, however, their reproducibility is criticized. The performance of feature selection algorithms is shown to be affected by the datasets, underlying networks and evaluation metrics. One of the causes is the curse of dimensionality, which makes it hard to select the features that generalize well on independent data. Even the integration of biological networks does not mitigate this issue because the networks are large and many of their components are not relevant for the phenotype of interest. With the availability of multi-omics data, integrative approaches are being developed to build more robust predictive models. In this scenario, the higher data dimensions create greater challenges. We proposed a phenotype relevant network-based feature selection (PRNFS) framework and demonstrated its advantages in lung cancer prognosis prediction. We constructed cancer prognosis relevant networks based on epithelial mesenchymal transition (EMT) and integrated them with different types of omics data for feature selection. With less than 2.5% of the total dimensionality, we obtained EMT prognostic signatures that achieved remarkable prediction performance (average AUC values above 0.8), very significant sample stratifications, and meaningful biological interpretations. In addition to finding EMT signatures from different omics data levels, we combined these single-omics signatures into multi-omics signatures, which improved sample stratifications significantly. Both single- and multi-omics EMT signatures were tested on independent multi-omics lung cancer datasets and significant sample stratifications were obtained. Y1 - 2019 U6 - https://doi.org/10.1371/journal.pone.0204186 VL - 14 IS - 1 ER - TY - JOUR A1 - Gupta, Pooja A1 - Peter, Sarah A1 - Jung, Markus A1 - Lewin, Astrid A1 - Hemmrich-Stanisak, Georg A1 - Franke, Andre A1 - von Kleist, Max A1 - Schütte, Christof A1 - Einspanier, Ralf A1 - Sharbati, Soroush A1 - zur Bruegge, Jennifer T1 - Analysis of long non-coding RNA and mRNA expression in bovine macrophages brings up novel 2 aspects of Mycobacterium avium subspecies paratuberculosis infections JF - Scientific Reports in Nature N2 - Paratuberculosis is a major disease in cattle that severely affects animal welfare and causes huge economic losses worldwide. Development of alternative diagnostic methods is of urgent need to control the disease. Recent studies suggest that long non-coding RNAs (lncRNAs) play a crucial role in regulating immune function and may confer valuable information about the disease. However, their role has not yet been investigated in cattle with respect to infection towards Paratuberculosis. Therefore, we investigated the alteration in genomic expression profiles of mRNA and lncRNA in bovine macrophages in response to Paratuberculosis infection using RNA-Seq. We identified 397 potentially novel lncRNA candidates in macrophages of which 38 were differentially regulated by the infection. A total of 820 coding genes were also significantly altered by the infection. Co-expression analysis of lncRNAs and their neighbouring coding genes suggest regulatory functions of lncRNAs in pathways related to immune response. For example, this included protein coding genes such as TNIP3, TNFAIP3 and NF-κB2 that play a role in NF-κB2 signalling, a pathway associated with immune response. This study advances our understanding of lncRNA roles during Paratuberculosis infection. Y1 - 2019 U6 - https://doi.org/10.1038/s41598-018-38141-x VL - 9 ER - TY - JOUR A1 - Schütte, Christof A1 - Wulkow, Michael T1 - A hybrid Galerkin-Monte-Carlo approach to higher-dimensional population balances in polymerization kinetics JF - Macromol. React. Eng. Y1 - 2010 VL - 4 SP - 562 EP - 577 ER - TY - BOOK A1 - Schütte, Christof A1 - Sarich, Marco T1 - Metastability and Markov State Models in Molecular Dynamics: Modeling, Analysis, Algorithmic Approaches Y1 - 2013 N1 - Courant Lecture Notes No. 32, American Mathematical Society PB - American Mathematical Society ER - TY - JOUR A1 - Paulus, Florian A1 - Weiss, Maximilian A1 - Steinhilber, Dirk A1 - Nikitin, Anatoly A1 - Schütte, Christof A1 - Haag, Rainer T1 - Anionic Ring-Opening Polymerization Simulations for Hyperbranched Polyglycerols with Defined Molecular Weights JF - Macromolecules Y1 - 2013 U6 - https://doi.org/10.1021/ma401712w VL - 46 IS - 21 SP - 8458 EP - 8466 ER - TY - JOUR A1 - Hartmann, Carsten A1 - Banisch, Ralf A1 - Sarich, Marco A1 - Badowski, Thomas A1 - Schütte, Christof T1 - Characterization of Rare Events in Molecular Dynamics JF - Entropy (Special Issue) Y1 - 2013 U6 - https://doi.org/10.3390/e16010350 VL - 16 IS - 1 SP - 350 EP - 376 ER - TY - JOUR A1 - Sarich, Marco A1 - Banisch, Ralf A1 - Hartmann, Carsten A1 - Schütte, Christof T1 - Markov State Models for Rare Events in Molecular Dynamics JF - Entropy (Special Issue) Y1 - 2013 U6 - https://doi.org/10.3390/e16010258 VL - 16 IS - 1 SP - 258 EP - 286 ER - TY - JOUR A1 - Wang, Han A1 - Hartmann, Carsten A1 - Schütte, Christof T1 - Linear response theory and optimal control for a molecular system under nonequilibrium conditions JF - Molecular Physics Y1 - 2013 UR - http://publications.imp.fu-berlin.de/1297/ U6 - https://doi.org/10.1080/00268976.2013.844370 VL - 111 SP - 3555 EP - 3564 ER - TY - JOUR A1 - Nikitin, B. A. A1 - Wulkow, Michael A1 - Schütte, Christof T1 - Modeling of Free Radical Styrene/Divinylbenzene Copolymerization with the Numerical Fractionation Technique JF - Macromolecular Theory and Simulation Y1 - 2013 U6 - https://doi.org/10.1002/mats.201300125 VL - 22 IS - 9 SP - 475 EP - 489 ER - TY - JOUR A1 - Enciso, Marta A1 - Schütte, Christof A1 - Delle Site, Luigi T1 - pH-dependent Response of Coiled Coils: A Coarse-Grained Molecular Simulation Study JF - Molecular Physics Y1 - 2013 U6 - https://doi.org/10.1080/00268976.2013.827254 VL - 111 IS - 22-23 SP - 3363 EP - 3371 ER - TY - JOUR A1 - Wang, Han A1 - Schütte, Christof A1 - Ciccotti, Giovanni A1 - Delle Site, Luigi T1 - Exploring the conformational dynamics of alanine dipeptide in solution subjected to an external electric field: A nonequilibrium molecular dynamics simulation JF - Journal of Chemical Theory and Computation Y1 - 2014 U6 - https://doi.org/10.1186/s12859-017-1565-4 VL - 10 IS - 4 SP - 1376 EP - 1386 ER - TY - GEN A1 - Weber, Marcus A1 - Fackeldey, Konstantin A1 - Schütte, Christof T1 - Set-free Markov State Building N2 - Molecular dynamics (MD) simulations face challenging problems since the timescales of interest often are much longer than what is possible to simulate and even if sufficiently long simulation are possible the complex nature of the resulting simulation data makes interpretation difficult. Markov State Models (MSMs) help to overcome these problems by making experimentally relevant timescales accessible via coarse grained representations that also allows for convenient interpretation. However, standard set-based MSMs exhibit some caveats limiting their approximation quality and statistical significance. One of the main caveats results from the fact that typical MD trajectories repeatedly re-cross the boundary between the sets used to build the MSM which causes statistical bias in estimating the transition probabilities between these sets. In this article, we present a set-free approach to MSM building utilizing smooth overlapping ansatz functions instead of sets and an adaptive refinement approach. This kind of meshless discretization helps to overcome the recrossing problem and yields an adaptive refinement procedure that allows to improve the quality of the model while exploring state space and inserting new ansatz functions into the MSM. T3 - ZIB-Report - 17-10 Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-62167 SN - 1438-0064 ER - TY - GEN A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - The spatiotemporal master equation: approximation of reaction-diffusion dynamics via Markov state modeling N2 - Accurate modeling and numerical simulation of reaction kinetics is a topic of steady interest. We consider the spatiotemporal chemical master equation (ST-CME) as a model for stochastic reaction-diffusion systems that exhibit properties of metastability. The space of motion is decomposed into metastable compartments and diffusive motion is approximated by jumps between these compartments. Treating these jumps as first-order reactions, simulation of the resulting stochastic system is possible by the Gillespie method. We present the theory of Markov state models (MSM) as a theoretical foundation of this intuitive approach. By means of Markov state modeling, both the number and shape of compartments and the transition rates between them can be determined. We consider the ST-CME for two reaction-diffusion systems and compare it to more detailed models. Moreover, a rigorous formal justification of the ST-CME by Galerkin projection methods is presented. T3 - ZIB-Report - 16-60 KW - reaction-diffusion KW - stochastic chemical kinetics KW - chemical master equation Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-60999 SN - 1438-0064 ER - TY - GEN A1 - Klebanov, Ilja A1 - Sikorski, Alexander A1 - Schütte, Christof A1 - Röblitz, Susanna T1 - Empirical Bayes Methods, Reference Priors, Cross Entropy and the EM Algorithm N2 - When estimating a probability density within the empirical Bayes framework, the non-parametric maximum likelihood estimate (NPMLE) usually tends to overfit the data. This issue is usually taken care of by regularization - a penalization term is subtracted from the marginal log-likelihood before the maximization step, so that the estimate favors smooth solutions, resulting in the so-called maximum penalized likelihood estimation (MPLE). The majority of penalizations currently in use are rather arbitrary brute-force solutions, which lack invariance under transformation of the parameters(reparametrization) and measurements. This contradicts the principle that, if the underlying model has several equivalent formulations, the methods of inductive inference should lead to consistent results. Motivated by this principle and using an information-theoretic point of view, we suggest an entropy-based penalization term that guarantees this kind of invariance. The resulting density estimate can be seen as a generalization of reference priors. Using the reference prior as a hyperprior, on the other hand, is argued to be a poor choice for regularization. We also present an insightful connection between the NPMLE, the cross entropy and the principle of minimum discrimination information suggesting another method of inference that contains the doubly-smoothed maximum likelihood estimation as a special case. T3 - ZIB-Report - 16-56 KW - parameter estimation KW - Bayesian inference KW - Bayesian hierarchical modeling KW - hyperparameter KW - hyperprior KW - EM algorithm KW - NPMLE KW - MPLE KW - DS-MLE KW - principle of maximum entropy KW - cross entropy KW - minimum discrimination information KW - reference prior KW - Jeffreys prior Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-61230 SN - 1438-0064 ER - TY - GEN A1 - Klebanov, Ilja A1 - Sikorski, Alexander A1 - Schütte, Christof A1 - Röblitz, Susanna T1 - Empirical Bayes Methods for Prior Estimation in Systems Medicine N2 - One of the main goals of mathematical modelling in systems medicine related to medical applications is to obtain patient-specific parameterizations and model predictions. In clinical practice, however, the number of available measurements for single patients is usually limited due to time and cost restrictions. This hampers the process of making patient-specific predictions about the outcome of a treatment. On the other hand, data are often available for many patients, in particular if extensive clinical studies have been performed. Therefore, before applying Bayes’ rule separately to the data of each patient (which is typically performed using a non-informative prior), it is meaningful to use empirical Bayes methods in order to construct an informative prior from all available data. We compare the performance of four priors - a non-informative prior and priors chosen by nonparametric maximum likelihood estimation (NPMLE), by maximum penalized lilelihood estimation (MPLE) and by doubly-smoothed maximum likelihood estimation (DS-MLE) - by applying them to a low-dimensional parameter estimation problem in a toy model as well as to a high-dimensional ODE model of the human menstrual cycle, which represents a typical example from systems biology modelling. T3 - ZIB-Report - 16-57 KW - Parameter estimation KW - Bayesian inference KW - Bayesian hierarchical modelling KW - NPMLE KW - MPLE KW - DS-MLE KW - EM algorithm KW - Jeffreys prior KW - reference prior KW - hyperparameter KW - hyperprior KW - principle of maximum entropy Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-61307 SN - 1438-0064 ER - TY - JOUR A1 - Hartmann, Carsten A1 - Schütte, Christof A1 - Weber, Marcus A1 - Zhang, Wei T1 - Importance sampling in path space for diffusion processes with slow-fast variables JF - Probability Theory and Related Fields N2 - Importance sampling is a widely used technique to reduce the variance of a Monte Carlo estimator by an appropriate change of measure. In this work, we study importance sampling in the framework of diffusion process and consider the change of measure which is realized by adding a control force to the original dynamics. For certain exponential type expectation, the corresponding control force of the optimal change of measure leads to a zero-variance estimator and is related to the solution of a Hamilton–Jacobi–Bellmann equation. We focus on certain diffusions with both slow and fast variables, and the main result is that we obtain an upper bound of the relative error for the importance sampling estimators with control obtained from the limiting dynamics. We demonstrate our approximation strategy with an illustrative numerical example. Y1 - 2017 UR - http://rdcu.be/oA51 U6 - https://doi.org/10.1007/s00440-017-0755-3 N1 - If you are not subscribing to the journal please use the http link below as part of the Springer Nature SharedIt initiative SP - 1 EP - 52 ER - TY - JOUR A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - The Spatiotemporal Master Equation: Approximation of Reaction-Diffusion Dynamics via Markov State Modeling JF - Journal of Chemical Physics N2 - Accurate modeling and numerical simulation of reaction kinetics is a topic of steady interest.We consider the spatiotemporal chemical master equation (ST-CME) as a model for stochastic reaction-diffusion systems that exhibit properties of metastability. The space of motion is decomposed into metastable compartments and diffusive motion is approximated by jumps between these compartments. Treating these jumps as first-order reactions, simulation of the resulting stochastic system is possible by the Gillespie method. We present the theory of Markov state models (MSM) as a theoretical foundation of this intuitive approach. By means of Markov state modeling, both the number and shape of compartments and the transition rates between them can be determined. We consider the ST-CME for two reaction-diffusion systems and compare it to more detailed models. Moreover, a rigorous formal justification of the ST-CME by Galerkin projection methods is presented. Y1 - 2016 U6 - https://doi.org/10.1063/1.4971163 VL - 145 IS - 21 ER - TY - JOUR A1 - Duwal, Sulav A1 - Winkelmann, Stefanie A1 - Schütte, Christof A1 - von Kleist, Max T1 - Optimal Treatment Strategies in the Context of 'Treatment for Prevention' against HIV/1 in Resource-Poor Settings JF - PloS Computational Biology N2 - An estimated 2.7 million new HIV-1 infections occurred in 2010. `Treatment-for-prevention’ may strongly prevent HIV-1 transmission. The basic idea is that immediate treatment initiation rapidly decreases virus burden, which reduces the number of transmittable viruses and thereby the probability of infection. However, HIV inevitably develops drug resistance, which leads to virus rebound and nullifies the effect of `treatment-for-prevention’ for the time it remains unrecognized. While timely conducted treatment changes may avert periods of viral rebound, necessary treatment options and diagnostics may be lacking in resource-constrained settings. Within this work, we provide a mathematical platform for comparing different treatment paradigms that can be applied to many medical phenomena. We use this platform to optimize two distinct approaches for the treatment of HIV-1: (i) a diagnostic-guided treatment strategy, based on infrequent and patient-specific diagnostic schedules and (ii) a pro-active strategy that allows treatment adaptation prior to diagnostic ascertainment. Both strategies are compared to current clinical protocols (standard of care and the HPTN052 protocol) in terms of patient health, economic means and reduction in HIV-1 onward transmission exemplarily for South Africa. All therapeutic strategies are assessed using a coarse-grained stochastic model of within-host HIV dynamics and pseudo-codes for solving the respective optimal control problems are provided. Our mathematical model suggests that both optimal strategies (i)-(ii) perform better than the current clinical protocols and no treatment in terms of economic means, life prolongation and reduction of HIV-transmission. The optimal diagnostic-guided strategy suggests rare diagnostics and performs similar to the optimal pro-active strategy. Our results suggest that ‘treatment-for-prevention’ may be further improved using either of the two analyzed treatment paradigms. Y1 - 2015 U6 - https://doi.org/10.1371/journal.pcbi.1004200 VL - 11 IS - 4 ER - TY - JOUR A1 - Fischer, Alexander A1 - Cordes, Frank A1 - Schütte, Christof T1 - Hybrid Monte Carlo with Adaptive Temperature in Mixed–Canonical Ensemble: Efficient conformational analysis of RNA JF - J. Comp. Chem. Y1 - 1998 U6 - https://doi.org/10.1002/(SICI)1096-987X(19981130)19:15<1689::AID-JCC2>3.0.CO;2-J VL - 19 IS - 15 SP - 1689 EP - 1697 ER - TY - CHAP A1 - Fischer, Alexander A1 - Schütte, Christof A1 - Deuflhard, Peter A1 - Cordes, Frank ED - Schlick, T. ED - Gan, H. T1 - Hierarchical Uncoupling-Coupling of Metastable Conformations T2 - Computational Methods for Macromolecules Y1 - 2002 UR - http://www.zib.de/PaperWeb/abstracts/ZR-01-03 IS - 24 SP - 235 EP - 259 PB - Springer ER - TY - JOUR A1 - Metzner, Ph. A1 - Weber, Marcus A1 - Schütte, Christof T1 - Observation uncertainty in reversible Markov chains JF - Phys. Rev. E Y1 - 2010 U6 - https://doi.org/10.1103/PhysRevE.82.031114 VL - 82 IS - 3 SP - 031114 PB - American Physical Society ER - TY - JOUR A1 - Schütte, Christof A1 - Fischer, Alexander A1 - Huisinga, Wilhelm A1 - Deuflhard, Peter T1 - A Direct Approach to Conformational Dynamics Based on Hybrid Monte Carlo JF - J. Comput. Phys. Y1 - 1999 UR - http://www.zib.de/PaperWeb/abstracts/SC-98-45 VL - 151 SP - 146 EP - 168 ER - TY - CHAP A1 - Schütte, Christof A1 - Huisinga, Wilhelm A1 - Deuflhard, Peter ED - Fiedler, Bernold T1 - Transfer operator approach to conformational dynamics in biomolecular systems T2 - Ergodic theory, analysis, and efficient simulation of dynamical systems Y1 - 2001 SP - 191 EP - 223 PB - Berlin: Springer ER - TY - JOUR A1 - Wang, Han A1 - Hartmann, Carsten A1 - Schütte, Christof A1 - Site, Luigi Delle T1 - Grand-canonical-like molecular-dynamics simulations by using an adaptive-resolution technique JF - Phys. Rev. X Y1 - 2013 U6 - https://doi.org/10.1103/PhysRevX.3.011018 VL - 3 SP - 011018 ER - TY - JOUR A1 - Sarich, Marco A1 - Djurdjevac Conrad, Natasa A1 - Bruckner, Sharon A1 - Conrad, Tim A1 - Schütte, Christof T1 - Modularity revisited: A novel dynamics-based concept for decomposing complex networks JF - Journal of Computational Dynamics Y1 - 2014 UR - http://publications.mi.fu-berlin.de/1127/ U6 - https://doi.org/10.3934/jcd.2014.1.191 VL - 1 IS - 1 SP - 191 EP - 212 ER - TY - JOUR A1 - Aiche, Stephan A1 - Reinert, Knut A1 - Schütte, Christof A1 - Hildebrand, Diana A1 - Schlüter, Hartmut A1 - Conrad, Tim T1 - Inferring Proteolytic Processes from Mass Spectrometry Time Series Data Using Degradation Graphs JF - PLoS ONE Y1 - 2012 UR - http://publications.imp.fu-berlin.de/1143/ U6 - https://doi.org/10.1371/journal.pone.0040656 VL - 7 IS - 7 SP - e40656 PB - Public Library of Science ER - TY - JOUR A1 - Djurdjevac Conrad, Natasa A1 - Bruckner, Stefanie A1 - Conrad, Tim A1 - Schütte, Christof T1 - Random Walks on Complex Modular Networks JF - Journal of Numerical Analysis, Industrial and Applied Mathematics Y1 - 2012 UR - http://publications.imp.fu-berlin.de/1090/ VL - 6 IS - 1-2 SP - 29 EP - 50 PB - European Society of Computational Methods in Sciences and Engineering ER - TY - JOUR A1 - Rath, Barbara A1 - Yousef, Kaveh A1 - Katzenstein, D. A1 - Shafer, R. A1 - Schütte, Christof A1 - von Kleist, Max A1 - Merigan, T. T1 - HIV-1 Evolution in Response to Triple Reverse Transcriptase Inhibitor Induced Selective Pressure in vitro JF - PLoS ONE Y1 - 2013 UR - http://publications.imp.fu-berlin.de/1178/ U6 - https://doi.org/10.1371/journal.pone.0061102 VL - 8 IS - 4 PB - Public Library of Science ER - TY - GEN A1 - Zhang, Wei A1 - Hartmann, Carsten A1 - Schütte, Christof T1 - Effective Dynamics Along Given Reaction Coordinates, and Reaction Rate Theory N2 - In molecular dynamics and related fields one considers dynamical descriptions of complex systems in full (atomic) detail. In order to reduce the overwhelming complexity of realistic systems (high dimension, large timescale spread, limited computational resources) the projection of the full dynamics onto some reaction coordinates is examined in order to extract statistical information like free energies or reaction rates. In this context, the effective dynamics that is induced by the full dynamics on the reaction coordinate space has attracted considerable attention in the literature. In this article, we contribute to this discussion: We first show that if we start with an ergodic diffusion processes whose invariant measure is unique then these properties are inherited by the effective dynamics. Then, we give equations for the effective dynamics, discuss whether the dominant timescales and reaction rates inferred from the effective dynamics are accurate approximations of such quantities for the full dynamics, and compare our findings to results from approaches like Zwanzig-Mori, averaging, or homogenization. Finally, by discussing the algorithmic realization of the effective dynamics, we demonstrate that recent algorithmic techniques like the ”equation-free” approach and the ”heterogeneous multiscale method” can be seen as special cases of our approach. T3 - ZIB-Report - 16-35 KW - Ergodic diffusion KW - reaction coordinate KW - effective dynamics KW - model reduction KW - equation-free KW - heterogeneous multiscale method Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-59706 SN - 1438-0064 ER - TY - GEN A1 - Ciccotti, Giovanni A1 - Ferrario, Mauro A1 - Schütte, Christof T1 - Molecular Dynamics vs. Stochastic Processes: Are We Heading Anywhere? T2 - Special Issue: Understanding Molecular Dynamics via Stochastic Processes, Entropy Y1 - 2018 U6 - https://doi.org/10.3390/e20050348 VL - 20 IS - 5 ER - TY - JOUR A1 - Zhang, Wei A1 - Schütte, Christof T1 - Reliable approximation of long relaxation timescales in molecular dynamics JF - Entropy Y1 - 2017 U6 - https://doi.org/10.3390/e19070367 VL - 19 IS - 7 ER - TY - JOUR A1 - Zhang, Wei A1 - Hartmann, Carsten A1 - Schütte, Christof T1 - Effective dynamics along given reaction coordinates, and reaction rate theory JF - Faraday Discussions Y1 - 2016 U6 - https://doi.org/10.1039/C6FD00147E IS - 195 SP - 365 EP - 394 ER - TY - JOUR A1 - Klebanov, Ilja A1 - Sikorski, Alexander A1 - Schütte, Christof A1 - Röblitz, Susanna T1 - Objective priors in the empirical Bayes framework JF - Scandinavian Journal of Statistics N2 - When dealing with Bayesian inference the choice of the prior often remains a debatable question. Empirical Bayes methods offer a data-driven solution to this problem by estimating the prior itself from an ensemble of data. In the nonparametric case, the maximum likelihood estimate is known to overfit the data, an issue that is commonly tackled by regularization. However, the majority of regularizations are ad hoc choices which lack invariance under reparametrization of the model and result in inconsistent estimates for equivalent models. We introduce a nonparametric, transformation-invariant estimator for the prior distribution. Being defined in terms of the missing information similar to the reference prior, it can be seen as an extension of the latter to the data-driven setting. This implies a natural interpretation as a trade-off between choosing the least informative prior and incorporating the information provided by the data, a symbiosis between the objective and empirical Bayes methodologies. Y1 - 2021 U6 - https://doi.org/10.1111/sjos.12485 VL - 48 IS - 4 SP - 1212 EP - 1233 PB - Wiley Online Library ER - TY - JOUR A1 - Hartmann, Carsten A1 - Richter, Lorenz A1 - Schütte, Christof A1 - Zhang, Wei T1 - Variational characterization of free energy: theory and algorithms JF - Entropy Y1 - 2017 U6 - https://doi.org/10.3390/e19110626 VL - 19 IS - 11 ER - TY - JOUR A1 - Dibak, Manuel A1 - J. del Razo, Mauricio A1 - De Sancho, David A1 - Schütte, Christof A1 - Noé, Frank T1 - MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations JF - Journal of Chemical Physics N2 - Molecular dynamics (MD) simulations can model the interactions between macromolecules with high spatiotemporal resolution but at a high computational cost. By combining high-throughput MD with Markov state models (MSMs), it is now possible to obtain long-timescale behavior of small to intermediate biomolecules and complexes. To model the interactions of many molecules at large lengthscales, particle-based reaction-diffusion (RD) simulations are more suitable but lack molecular detail. Thus, coupling MSMs and RD simulations (MSM/RD) would be highly desirable, as they could efficiently produce simulations at large time- and lengthscales, while still conserving the characteristic features of the interactions observed at atomic detail. While such a coupling seems straightforward, fundamental questions are still open: Which definition of MSM states is suitable? Which protocol to merge and split RD particles in an association/dissociation reaction will conserve the correct bimolecular kinetics and thermodynamics? In this paper, we make the first step towards MSM/RD by laying out a general theory of coupling and proposing a first implementation for association/dissociation of a protein with a small ligand (A + B <--> C). Applications on a toy model and CO diffusion into the heme cavity of myoglobin are reported. Y1 - 2018 U6 - https://doi.org/10.1063/1.5020294 VL - 148 IS - 214107 ER - TY - JOUR A1 - Djurdjevac Conrad, Natasa A1 - Helfmann, Luzie A1 - Zonker, Johannes A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - Human mobility and innovation spreading in ancient times: a stochastic agent-based simulation approach JF - EPJ Data Science N2 - Human mobility always had a great influence on the spreading of cultural, social and technological ideas. Developing realistic models that allow for a better understanding, prediction and control of such coupled processes has gained a lot of attention in recent years. However, the modeling of spreading processes that happened in ancient times faces the additional challenge that available knowledge and data is often limited and sparse. In this paper, we present a new agent-based model for the spreading of innovations in the ancient world that is governed by human movements. Our model considers the diffusion of innovations on a spatial network that is changing in time, as the agents are changing their positions. Additionally, we propose a novel stochastic simulation approach to produce spatio-temporal realizations of the spreading process that are instructive for studying its dynamical properties and exploring how different influences affect its speed and spatial evolution. Y1 - 2018 U6 - https://doi.org/10.1140/epjds/s13688-018-0153-9 VL - 7 IS - 1 PB - EPJ Data Science ET - EPJ Data Science ER - TY - JOUR A1 - Möller, Jan A1 - Isbilir, Ali A1 - Sungkaworn, Titiwat A1 - Osberg, Brenda A1 - Karathanasis, Christos A1 - Sunkara, Vikram A1 - Grushevsky, Eugene O A1 - Bock, Andreas A1 - Annibale, Paolo A1 - Heilemann, Mike A1 - Schütte, Christof A1 - Lohse, Martin J. T1 - Single molecule mu-opioid receptor membrane-dynamics reveal agonist-specific dimer formation with super-resolved precision JF - Nature Chemical Biology Y1 - 2020 U6 - https://doi.org/10.1038/s41589-020-0566-1 VL - 16 SP - 946 EP - 954 ER - TY - JOUR A1 - Bittracher, Andreas A1 - Schütte, Christof T1 - A probabilistic algorithm for aggregating vastly undersampled large Markov chains JF - Physica D: Nonlinear Phenomena Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-75874 VL - 416 ER - TY - JOUR A1 - Bittracher, Andreas A1 - Klus, Stefan A1 - Hamzi, Boumediene A1 - Schütte, Christof T1 - Dimensionality Reduction of Complex Metastable Systems via Kernel Embeddings of Transition Manifolds JF - Journal of Nonlinear Science N2 - We present a novel kernel-based machine learning algorithm for identifying the low-dimensional geometry of the effective dynamics of high-dimensional multiscale stochastic systems. Recently, the authors developed a mathematical framework for the computation of optimal reaction coordinates of such systems that is based on learning a parameterization of a low-dimensional transition manifold in a certain function space. In this article, we enhance this approach by embedding and learning this transition manifold in a reproducing kernel Hilbert space, exploiting the favorable properties of kernel embeddings. Under mild assumptions on the kernel, the manifold structure is shown to be preserved under the embedding, and distortion bounds can be derived. This leads to a more robust and more efficient algorithm compared to the previous parameterization approaches. Y1 - 2021 U6 - https://doi.org/10.1007/s00332-020-09668-z VL - 31 ER - TY - JOUR A1 - Klus, Stefan A1 - Nüske, Feliks A1 - Peitz, Sebastian A1 - Niemann, Jan-Hendrik A1 - Clementi, Cecilia A1 - Schütte, Christof T1 - Data-driven approximation of the Koopman generator: Model reduction, system identification, and control JF - Physica D: Nonlinear Phenomena Y1 - 2020 U6 - https://doi.org/10.1016/j.physd.2020.132416 VL - 406 ER - TY - BOOK A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - Stochastic Dynamics in Computational Biology T3 - Frontiers in Applied Dynamical Systems: Reviews and Tutorials Y1 - 2020 SN - 978-3-030-62386-9 U6 - https://doi.org/10.1007/978-3-030-62387-6 VL - 8 PB - Springer International Publishing ER - TY - JOUR A1 - Schulz, Robert A1 - Yamamoto, Kenji A1 - Klossek, André A1 - Rancan, Fiorenza A1 - Vogt, Annika A1 - Schütte, Christof A1 - Rühl, Eckart A1 - Netz, Roland R. T1 - Modeling of Drug Diffusion Based on Concentration Profiles in Healthy and Damaged Human Skin JF - Biophysical Journal N2 - Based on experimental drug concentration profiles in healthy as well as tape-stripped ex vivo human skin, we model the penetration of the antiinflammatory drug dexamethasone into the skin layers by the one-dimensional generalized diffusion equation. We estimate the position-dependent free-energy and diffusivity profiles by solving the conjugated minimization problem, in which the only inputs are concentration profiles of dexamethasone in skin at three consecutive penetration times. The resulting free-energy profiles for damaged and healthy skin show only minor differences. In contrast, the drug diffusivity in the first 10 μm of the upper skin layer of damaged skin is 200-fold increased compared to healthy skin, which reflects the corrupted barrier function of tape-stripped skin. For the case of healthy skin, we examine the robustness of our method by analyzing the behavior of the extracted skin parameters when the number of input and output parameters are reduced. We also discuss techniques for the regularization of our parameter extraction method. Y1 - 2019 U6 - https://doi.org/10.1016/j.bpj.2019.07.027 VL - 117 IS - 5 SP - 998 EP - 1008 ER - TY - JOUR A1 - Becker, Fabian A1 - Djurdjevac Conrad, Natasa A1 - Eser, Raphael A. A1 - Helfmann, Luzie A1 - Schütt, Brigitta A1 - Schütte, Christof A1 - Zonker, Johannes T1 - The Furnace and the Goat—A spatio-temporal model of the fuelwood requirement for iron metallurgy on Elba Island, 4th century BCE to 2nd century CE JF - PLOS ONE Y1 - 2020 UR - https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0241133 U6 - https://doi.org/10.1371/journal.pone.0241133 VL - 15 SP - 1 EP - 37 PB - Public Library of Science ER - TY - GEN A1 - Bittracher, Andreas A1 - Schütte, Christof ED - Junge, Oliver ED - Schütze, O. ED - Froyland, Gary ED - Ober-Blobaum, S. ED - Padberg-Gehle, E. T1 - A weak characterization of slow variables in stochastic dynamical systems T2 - Advances in Dynamics, Optimization and Computation. Series: Studies in Systems, Decision and Control. A volume dedicated to Michael Dellnitz on the occasion of his 60th birthday Y1 - 2020 SN - 978-3-030-51264-4 U6 - https://doi.org/10.1007/978-3-030-51264-4_6 VL - 304 SP - 132 EP - 150 PB - Springer International ER - TY - GEN A1 - Mollenhauer, Mattes A1 - Schuster, Ingmar A1 - Klus, Stefan A1 - Schütte, Christof ED - Junge, Oliver ED - Schütze, O. ED - Froyland, Gary ED - Ober-Blobaum, S. ED - Padberg-Gehle, K. T1 - Singular Value Decomposition of Operators on Reproducing Kernel Hilbert Spaces T2 - Advances om Dynamics, Optimization and Computation. Series: Studies in Systems, Decision and Control. A volume dedicated to Michael Dellnitz on his 60th birthday Y1 - 2020 SN - 978-3-030-51264-4 U6 - https://doi.org/10.1007/978-3-030-51264-4_5 VL - 304 SP - 109 EP - 131 PB - Springer International ER - TY - JOUR A1 - Helfmann, Luzie A1 - Ribera Borrell, Enric A1 - Schütte, Christof A1 - Koltai, Peter T1 - Extending Transition Path Theory: Periodically Driven and Finite-Time Dynamics JF - Journal of Nonlinear Science Y1 - 2020 U6 - https://doi.org/https://doi.org/10.1007/s00332-020-09652-7 VL - 30 SP - 3321 EP - 3366 ER -