@article{ConradGenzelCvetkovicetal.2017, author = {Conrad, Tim and Genzel, Martin and Cvetkovic, Nada and Wulkow, Niklas and Vybiral, Jan and Kutyniok, Gitta and Sch{\"u}tte, Christof}, title = {Sparse Proteomics Analysis - a compressed sensing-based approach for feature selection and classification of high-dimensional proteomics mass spectrometry data}, volume = {18}, journal = {BMC Bioinformatics}, number = {160}, doi = {10.1186/s12859-017-1565-4}, pages = {1 -- 20}, year = {2017}, abstract = {Motivation: High-throughput proteomics techniques, such as mass spectrometry (MS)-based approaches, produce very high-dimensional data-sets. In a clinical setting one is often interested how MS spectra dier between patients of different classes, for example spectra from healthy patients vs. spectra from patients having a particular disease. Machine learning algorithms are needed to (a) identify these discriminating features and (b) classify unknown spectra based on this feature set. Since the acquired data is usually noisy, the algorithms should be robust to noise and outliers, and the identied feature set should be as small as possible. Results: We present a new algorithm, Sparse Proteomics Analysis (SPA), based on the theory of Compressed Sensing that allows to identify a minimal discriminating set of features from mass spectrometry data-sets. We show how our method performs on artificial and real-world data-sets.}, language = {en} } @article{GelssKlusEisertetal.2019, author = {Gelß, Patrick and Klus, Stefan and Eisert, Jens and Sch{\"u}tte, Christof}, title = {Multidimensional Approximation of Nonlinear Dynamical Systems}, volume = {14}, journal = {Journal of Computational and Nonlinear Dynamics}, number = {6}, doi = {10.1115/1.4043148}, year = {2019}, abstract = {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.}, language = {en} } @article{KoltaiSchuette2018, author = {Koltai, P{\´e}ter and Sch{\"u}tte, Christof}, title = {A multiscale perturbation expansion approach for Markov state modeling of non-stationary molecular dynamics}, volume = {16}, journal = {SIAM J. Multiscale Model. Simul.}, number = {4}, publisher = {SIAM}, doi = {10.1137/17M1146403}, pages = {1455 -- 1485}, year = {2018}, abstract = {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.}, language = {en} } @article{ShaoBjaanaesHellandetal.2019, author = {Shao, Borong and Bjaanaes, Maria and Helland, Aslaug and Sch{\"u}tte, Christof and Conrad, Tim}, title = {EMT network-based feature selection improves prognosis prediction in lung adenocarcinoma}, volume = {14}, journal = {PLOS ONE}, number = {1}, doi = {10.1371/journal.pone.0204186}, year = {2019}, abstract = {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.}, language = {en} } @article{GuptaPeterJungetal.2019, author = {Gupta, Pooja and Peter, Sarah and Jung, Markus and Lewin, Astrid and Hemmrich-Stanisak, Georg and Franke, Andre and von Kleist, Max and Sch{\"u}tte, Christof and Einspanier, Ralf and Sharbati, Soroush and zur Bruegge, Jennifer}, title = {Analysis of long non-coding RNA and mRNA expression in bovine macrophages brings up novel 2 aspects of Mycobacterium avium subspecies paratuberculosis infections}, volume = {9}, journal = {Scientific Reports in Nature}, doi = {10.1038/s41598-018-38141-x}, year = {2019}, abstract = {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.}, language = {en} } @article{WinkelmannSchuette2016, author = {Winkelmann, Stefanie and Sch{\"u}tte, Christof}, title = {The Spatiotemporal Master Equation: Approximation of Reaction-Diffusion Dynamics via Markov State Modeling}, volume = {145}, journal = {Journal of Chemical Physics}, number = {21}, doi = {10.1063/1.4971163}, year = {2016}, abstract = {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.}, language = {en} } @article{DuwalWinkelmannSchuetteetal.2015, author = {Duwal, Sulav and Winkelmann, Stefanie and Sch{\"u}tte, Christof and von Kleist, Max}, title = {Optimal Treatment Strategies in the Context of 'Treatment for Prevention' against HIV/1 in Resource-Poor Settings}, volume = {11}, journal = {PloS Computational Biology}, number = {4}, doi = {10.1371/journal.pcbi.1004200}, year = {2015}, abstract = {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.}, language = {en} } @article{SarichDjurdjevacConradBruckneretal.2014, author = {Sarich, Marco and Djurdjevac Conrad, Natasa and Bruckner, Sharon and Conrad, Tim and Sch{\"u}tte, Christof}, title = {Modularity revisited: A novel dynamics-based concept for decomposing complex networks}, volume = {1}, journal = {Journal of Computational Dynamics}, number = {1}, doi = {10.3934/jcd.2014.1.191}, pages = {191 -- 212}, year = {2014}, language = {en} } @article{DibakJdelRazoDeSanchoetal.2018, author = {Dibak, Manuel and J. del Razo, Mauricio and De Sancho, David and Sch{\"u}tte, Christof and No{\´e}, Frank}, title = {MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations}, volume = {148}, journal = {Journal of Chemical Physics}, number = {214107}, doi = {10.1063/1.5020294}, year = {2018}, abstract = {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.}, language = {en} } @article{KlusNueskePeitzetal.2020, author = {Klus, Stefan and N{\"u}ske, Feliks and Peitz, Sebastian and Niemann, Jan-Hendrik and Clementi, Cecilia and Sch{\"u}tte, Christof}, title = {Data-driven approximation of the Koopman generator: Model reduction, system identification, and control}, volume = {406}, journal = {Physica D: Nonlinear Phenomena}, doi = {10.1016/j.physd.2020.132416}, year = {2020}, language = {en} } @misc{MollenhauerSchusterKlusetal.2020, author = {Mollenhauer, Mattes and Schuster, Ingmar and Klus, Stefan and Sch{\"u}tte, Christof}, title = {Singular Value Decomposition of Operators on Reproducing Kernel Hilbert Spaces}, volume = {304}, journal = {Advances om Dynamics, Optimization and Computation. Series: Studies in Systems, Decision and Control. A volume dedicated to Michael Dellnitz on his 60th birthday}, editor = {Junge, Oliver and Sch{\"u}tze, O. and Froyland, Gary and Ober-Blobaum, S. and Padberg-Gehle, K.}, publisher = {Springer International}, isbn = {978-3-030-51264-4}, doi = {10.1007/978-3-030-51264-4_5}, pages = {109 -- 131}, year = {2020}, language = {en} } @article{RaySunkaraSchuetteetal.2020, author = {Ray, Sourav and Sunkara, Vikram and Sch{\"u}tte, Christof and Weber, Marcus}, title = {How to calculate pH-dependent binding rates for receptor-ligand systems based on thermodynamic simulations with different binding motifs}, volume = {46}, journal = {Molecular Simulation}, number = {18}, publisher = {Taylor and Francis}, doi = {10.1080/08927022.2020.1839660}, pages = {1443 -- 1452}, year = {2020}, abstract = {Molecular simulations of ligand-receptor interactions are a computational challenge, especially when their association- ('on'-rate) and dissociation- ('off'-rate) mechanisms are working on vastly differing timescales. One way of tackling this multiscale problem is to compute the free-energy landscapes, where molecular dynamics (MD) trajectories are used to only produce certain statistical ensembles. The approach allows for deriving the transition rates between energy states as a function of the height of the activation-energy barriers. In this article, we derive the association rates of the opioids fentanyl and N-(3-fluoro-1-phenethylpiperidin-4-yl)-N-phenyl propionamide (NFEPP) in a μ-opioid receptor by combining the free-energy landscape approach with the square-root-approximation method (SQRA), which is a particularly robust version of Markov modelling. The novelty of this work is that we derive the association rates as a function of the pH level using only an ensemble of MD simulations. We also verify our MD-derived insights by reproducing the in vitro study performed by the Stein Lab.}, language = {en} } @article{KoltaiWuNoeetal.2018, author = {Koltai, Peter and Wu, Hao and No{\´e}, Frank and Sch{\"u}tte, Christof}, title = {Optimal data-driven estimation of generalized Markov state models for non-equilibrium dynamics}, volume = {6}, journal = {Computation}, number = {1}, publisher = {MDPI}, address = {Basel, Switzerland}, doi = {10.3390/computation6010022}, pages = {22}, year = {2018}, language = {en} } @article{DibakdelRazodeSanchoetal.2018, author = {Dibak, Manuel and del Razo, Mauricio J. and de Sancho, David and Sch{\"u}tte, Christof and No{\´e}, Frank}, title = {MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations}, volume = {148}, journal = {Journal of Chemical Physics}, number = {21}, doi = {10.1063/1.5020294}, year = {2018}, abstract = {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 time-scale behavior of small to intermediate biomolecules and complexes. To model the interactions of many molecules at large length scales, 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 length scales, 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 toward 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.}, language = {en} } @article{WulkowConradDjurdjevacConradetal.2021, author = {Wulkow, Hanna and Conrad, Tim and Djurdjevac Conrad, Natasa and M{\"u}ller, Sebastian A. and Nagel, Kai and Sch{\"u}tte, Christof}, title = {Prediction of Covid-19 spreading and optimal coordination of counter-measures: From microscopic to macroscopic models to Pareto fronts}, volume = {16}, journal = {PLOS One}, number = {4}, publisher = {Public Library of Science}, doi = {10.1371/journal.pone.0249676}, year = {2021}, language = {en} } @article{BanischDjurdjevacConradSchuette2015, author = {Banisch, Ralf and Djurdjevac Conrad, Natasa and Sch{\"u}tte, Christof}, title = {Reactive flows and unproductive cycles for random walks on complex networks}, journal = {The European Physical Journal Special Topics, vol. 224, iss. 12 (2015) pp. 2369-2387}, doi = {10.1140/epjst/e2015-02417-8}, year = {2015}, language = {en} } @misc{BanischDjurdjevacConradSchuette2015, author = {Banisch, Ralf and Djurdjevac Conrad, Natasa and Sch{\"u}tte, Christof}, title = {Reactive flows and unproductive cycles for random walks on complex networks}, issn = {1438-0064}, doi = {10.1140/epjst/e2015-02417-8}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-54239}, year = {2015}, abstract = {We present a comprehensive theory for analysis and understanding of transition events between an initial set A and a target set B for general ergodic finite-state space Markov chains or jump processes, including random walks on networks as they occur, e.g., in Markov State Modelling in molecular dynamics. The theory allows us to decompose the probability flow generated by transition events between the sets A and B into the productive part that directly flows from A to B through reaction pathways and the unproductive part that runs in loops and is supported on cycles of the underlying network. It applies to random walks on directed networks and nonreversible Markov processes and can be seen as an extension of Transition Path Theory. Information on reaction pathways and unproductive cycles results from the stochastic cycle decomposition of the underlying network which also allows to compute their corresponding weight, thus characterizing completely which structure is used how often in transition events. The new theory is illustrated by an application to a Markov State Model resulting from weakly damped Langevin dynamics where the unproductive cycles are associated with periodic orbits of the underlying Hamiltonian dynamics.}, language = {en} } @article{WinkelmannSchuettevonKleist2014, author = {Winkelmann, Stefanie and Sch{\"u}tte, Christof and von Kleist, Max}, title = {Markov Control Processes with Rare State Observation: Theory and Application to Treatment Scheduling in HIV-1}, volume = {12}, journal = {Communications in Mathematical Sciences}, number = {5}, doi = {10.4310/CMS.2014.v12.n5.a4}, pages = {859 -- 877}, year = {2014}, abstract = {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 immunodeficiency 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.}, language = {en} } @article{DjurdjevacConradBanischSchuette2015, author = {Djurdjevac Conrad, Natasa and Banisch, Ralf and Sch{\"u}tte, Christof}, title = {Modularity of Directed Networks: Cycle Decomposition Approach}, journal = {Journal of Computational Dynamics 2 (2015) pp. 1-24}, doi = {10.3934/jcd.2015.2.1}, year = {2015}, abstract = {The problem of decomposing networks into modules (or clusters) has gained much attention in recent years, as it can account for a coarsegrained description of complex systems, often revealing functional subunits of these systems. A variety of module detection algorithms have been proposed, mostly oriented towards finding hard partitionings of undirected networks. Despite the increasing number of fuzzy clustering methods for directed networks, many of these approaches tend to neglect important directional information. In this paper, we present a novel random walk based approach for finding fuzzy partitions of directed, weighted networks, where edge directions play a crucial role in defining how well nodes in a module are interconnected. We will show that cycle decomposition of a random walk process connects the notion of network modules and information transport in a network, leading to a new, symmetric measure of node communication. Finally, we will use this measure to introduce a communication graph, for which we will show that although being undirected it inherits all necessary information about modular structures from the original network.}, language = {en} } @misc{WangSchuette2014, author = {Wang, Han and Sch{\"u}tte, Christof}, title = {Building Markov State Models for Periodically Driven Non-Equilibrium Systems}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-53167}, year = {2014}, abstract = {Recent years have seen an increased interest in non-equilibrium molecular dynamics (NEMD) simulations, especially for molecular systems with periodic forcing by external fields, e.g., in the context of studying effects of electromagnetic radiation on the human body tissue. Lately, an NEMD methods with local thermostating has been proposed that allows for studying non-equilibrium processes in a statistically reliable and thermodynamically consistent way. In this article, we demonstrate how to construct Markov State Models (MSMs) for such NEMD simulations. MSM building has been well-established for systems in equilibrium where MSMs with just a few (macro-)states allow for accurate reproduction of the essential kinetics of the molecular system under consideration. Non-equilibrium MSMs have been lacking so far. The article presents how to construct such MSMs and illustrates their validity and usefulness for the case of conformation dynamics of alanine dipeptide in an external electric field.}, language = {en} } @misc{GulSchuetteBernhard2015, author = {Gul, Raheem and Sch{\"u}tte, Christof and Bernhard, Stefan}, title = {Mathematical modeling and sensitivity analysis of arterial anastomosis in arm arteries}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-54339}, year = {2015}, abstract = {Cardiovascular diseases are one of the major problems in medicine today and the number of patients increases worldwide. To find the most efficient treatment, prior knowledge about function and dysfunction of the cardiovas- cular system is required and methods need to be developed that identify the disease in an early stage. Mathematical modeling is a powerful tool for prediction and investigation of cardiovascular diseases. It has been shown that the Windkessel model, being based on an analogy between electrical circuits and fluid flow, is a simple but effective method to model the human cardiovascular system. In this paper, we have applied parametric local sensitivity analysis (LSA) to a linear elastic model of the arm arteries, to find and rank sensitive param- eters that may be helpful in clinical diagnosis. A computational model for end-to-side anastomosis (superior ulnar collateral anastomosis with posterior ulnar recurrent, SUC-PUR) is carried out to study the effects of some clinically relevant haemodynamic parameters like blood flow resistance and terminal re- sistance on pressure and flow at different locations of the arm artery. In this context, we also discuss the spatio-temporal dependency of local sensitivities. The sensitivities with respect to cardiovascular parameters reveal the flow resistance and diameter of the vessels as most sensitive parameters. These parameters play a key role in diagnosis of severe stenosis and aneurysms. In contrast, wall thickness and elastic modulus are found to be less sensitive.}, language = {en} }