@article{SchuetteWinkelmannHartmann2012, author = {Sch{\"u}tte, Christof and Winkelmann, Stefanie and Hartmann, Carsten}, title = {Optimal control of molecular dynamics using Markov state models}, volume = {134}, journal = {Math. Program. (Series B)}, number = {1}, doi = {10.1007/s10107-012-0547-6}, pages = {259 -- 282}, year = {2012}, language = {en} } @article{SchuetteWulkow2010, author = {Sch{\"u}tte, Christof and Wulkow, Michael}, title = {A hybrid Galerkin-Monte-Carlo approach to higher-dimensional population balances in polymerization kinetics}, volume = {4}, journal = {Macromol. React. Eng.}, pages = {562 -- 577}, year = {2010}, language = {en} } @article{SchuetteWulkow2010, author = {Sch{\"u}tte, Christof and Wulkow, Michael}, title = {A Hybrid Galerkin-Monte-Carlo Approach to Higher-Dimensional Population Balances in Polymerization Kinetics}, volume = {4}, journal = {Macromol. React. Eng.}, number = {9-10}, doi = {10.1002/mren.200900073}, pages = {562 -- 577}, year = {2010}, language = {en} } @article{SchuetteWulkow1992, author = {Sch{\"u}tte, Christof and Wulkow, Michael}, title = {Quantum Theory with Discrete Spectra and Countable Systems of Differential Equations - A Numerical Treatment of Raman Spectroscopy.}, journal = {preprint}, year = {1992}, language = {en} } @misc{SchuetteWulkow1992, author = {Sch{\"u}tte, Christof and Wulkow, Michael}, title = {Quantum Theory with Discrete Spectra and Countable Systems of Differential Equations - A Numerical Treatment of RamanSpectroscopy.}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-778}, number = {SC-92-07}, year = {1992}, abstract = {Models for occupation dynamics in discrete quantum systems lead to large or even infinite systems of ordinary differential equations. Some new mathematical techniques, developed for the simulation of chemical processes, make a numerical solution of countable systems of ordinary differential equations possible. Both, a basic physical concept for the construction of such systems and the structure of the numerical tools for solving them are presented. These conceptual aspects are illustrated by a simulation of an occupation process from spectroscopy. In this example the structures of rotation spectra observed in infrared spectroscopy are explained and some possibilities for an extension of the model are shown.}, language = {en} } @article{SchuetteZumbuschBrinkmann1995, author = {Sch{\"u}tte, Christof and Zumbusch, Gerhard and Brinkmann, Ralf}, title = {Dynamics of Erbium-doped Waveguide Lasers}, journal = {preprint}, year = {1995}, language = {en} } @article{SeckerFackeldeyWeberetal.2023, author = {Secker, Christopher and Fackeldey, Konstantin and Weber, Marcus and Ray, Sourav and Gorgulla, Christoph and Sch{\"u}tte, Christof}, title = {Novel multi-objective affinity approach allows to identify pH-specific μ-opioid receptor agonists}, volume = {15}, journal = {Journal of Cheminformatics}, doi = {10.1186/s13321-023-00746-4}, year = {2023}, abstract = {Opioids are essential pharmaceuticals due to their analgesic properties, however, lethal side effects, addiction, and opioid tolerance are extremely challenging. The development of novel molecules targeting the μ-opioid receptor (MOR) in inflamed, but not in healthy tissue, could significantly reduce these unwanted effects. Finding such novel molecules can be achieved by maximizing the binding affinity to the MOR at acidic pH while minimizing it at neutral pH, thus combining two conflicting objectives. Here, this multi-objective optimal affinity approach is presented, together with a virtual drug discovery pipeline for its practical implementation. When applied to finding pH-specific drug candidates, it combines protonation state-dependent structure and ligand preparation with high-throughput virtual screening. We employ this pipeline to characterize a set of MOR agonists identifying a morphine-like opioid derivative with higher predicted binding affinities to the MOR at low pH compared to neutral pH. Our results also confirm existing experimental evidence that NFEPP, a previously described fentanyl derivative with reduced side effects, and recently reported β-fluorofentanyls and -morphines show an increased specificity for the MOR at acidic pH when compared to fentanyl and morphine. We further applied our approach to screen a >50K ligand library identifying novel molecules with pH-specific predicted binding affinities to the MOR. The presented differential docking pipeline can be applied to perform multi-objective affinity optimization to identify safer and more specific drug candidates at large scale.}, language = {en} } @article{SenneTrendelkampSchroerMeyetal.2012, author = {Senne, M. and Trendelkamp-Schroer, B. and Mey, A. and Sch{\"u}tte, Christof and No{\´e}, Frank}, title = {EMMA - A software package for Markov model building and analysis}, volume = {8}, journal = {Journal of Chemical Theory and Computation}, doi = {10.1021/ct300274u}, pages = {2223 -- 2238}, year = {2012}, 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{SherrattSrivastavaAinslieetal.2024, author = {Sherratt, Katharine and Srivastava, Ajitesh and Ainslie, Kylie and Singh, David E. and Cublier, Aymar and Marinescu, Maria Cristina and Carretero, Jesus and Garcia, Alberto Cascajo and Franco, Nicolas and Willem, Lander and Abrams, Steven and Faes, Christel and Beutels, Philippe and Hens, Niel and M{\"u}ller, Sebastian and Charlton, Billy and Ewert, Ricardo and Paltra, Sydney and Rakow, Christian and Rehmann, Jakob and Conrad, Tim and Sch{\"u}tte, Christof and Nagel, Kai and Abbott, Sam and Grah, Rok and Niehus, Rene and Prasse, Bastian and Sandmann, Frank and Funk, Sebastian}, title = {Characterising information gains and losses when collecting multiple epidemic model outputs}, volume = {47}, journal = {Epidemics}, publisher = {Elsevier BV}, issn = {1755-4365}, doi = {10.1016/j.epidem.2024.100765}, year = {2024}, abstract = {Collaborative comparisons and combinations of epidemic models are used as policy-relevant evidence during epidemic outbreaks. In the process of collecting multiple model projections, such collaborations may gain or lose relevant information. Typically, modellers contribute a probabilistic summary at each time-step. We compared this to directly collecting simulated trajectories. We aimed to explore information on key epidemic quantities; ensemble uncertainty; and performance against data, investigating potential to continuously gain information from a single cross-sectional collection of model results. Methods We compared July 2022 projections from the European COVID-19 Scenario Modelling Hub. Five modelling teams projected incidence in Belgium, the Netherlands, and Spain. We compared projections by incidence, peaks, and cumulative totals. We created a probabilistic ensemble drawn from all trajectories, and compared to ensembles from a median across each model's quantiles, or a linear opinion pool. We measured the predictive accuracy of individual trajectories against observations, using this in a weighted ensemble. We repeated this sequentially against increasing weeks of observed data. We evaluated these ensembles to reflect performance with varying observed data. Results. By collecting modelled trajectories, we showed policy-relevant epidemic characteristics. Trajectories contained a right-skewed distribution well represented by an ensemble of trajectories or a linear opinion pool, but not models' quantile intervals. Ensembles weighted by performance typically retained the range of plausible incidence over time, and in some cases narrowed this by excluding some epidemic shapes. Conclusions. We observed several information gains from collecting modelled trajectories rather than quantile distributions, including potential for continuously updated information from a single model collection. The value of information gains and losses may vary with each collaborative effort's aims, depending on the needs of projection users. Understanding the differing information potential of methods to collect model projections can support the accuracy, sustainability, and communication of collaborative infectious disease modelling efforts. Data availability All code and data available on Github: https://github.com/covid19-forecast-hub-europe/aggregation-info-loss}, language = {en} } @article{SikorskiWeberSchuette2021, author = {Sikorski, Alexander and Weber, Marcus and Sch{\"u}tte, Christof}, title = {The Augmented Jump Chain}, volume = {4}, journal = {Advanced Theory and Simulations}, number = {4}, publisher = {Wiley-VCH}, arxiv = {http://arxiv.org/abs/2008.04624}, doi = {10.1002/adts.202000274}, year = {2021}, abstract = {Modern methods of simulating molecular systems are based on the mathematical theory of Markov operators with a focus on autonomous equilibrated systems. However, non-autonomous physical systems or non-autonomous simulation processes are becoming more and more important. A representation of non-autonomous Markov jump processes is presented as autonomous Markov chains on space-time. Augmenting the spatial information of the embedded Markov chain by the temporal information of the associated jump times, the so-called augmented jump chain is derived. The augmented jump chain inherits the sparseness of the infinitesimal generator of the original process and therefore provides a useful tool for studying time-dependent dynamics even in high dimensions. Furthermore, possible generalizations and applications to the computation of committor functions and coherent sets in the non-autonomous setting are discussed. After deriving the theoretical foundations, the concepts with a proof-of-concept Galerkin discretization of the transfer operator of the augmented jump chain applied to simple examples are illustrated.}, language = {en} } @article{StraubeWinkelmannSchuetteetal.2021, author = {Straube, Arthur and Winkelmann, Stefanie and Sch{\"u}tte, Christof and H{\"o}fling, Felix}, title = {Stochastic pH oscillations in a model of the urea-urease reaction confined to lipid vesicles}, volume = {12}, journal = {J. Phys. Chem. Lett.}, doi = {10.1021/acs.jpclett.1c03016}, pages = {9888 -- 9893}, year = {2021}, abstract = {The urea-urease clock reaction is a pH switch from acid to basic that can turn into a pH oscillator if it occurs inside a suitable open reactor. We numerically study the confinement of the reaction to lipid vesicles, which permit the exchange with an external reservoir by differential transport, enabling the recovery of the pH level and yielding a constant supply of urea molecules. For microscopically small vesicles, the discreteness of the number of molecules requires a stochastic treatment of the reaction dynamics. Our analysis shows that intrinsic noise induces a significant statistical variation of the oscillation period, which increases as the vesicles become smaller. The mean period, however, is found to be remarkably robust for vesicle sizes down to approximately 200 nm, but the periodicity of the rhythm is gradually destroyed for smaller vesicles. The observed oscillations are explained as a canard-like limit cycle that differs from the wide class of conventional feedback oscillators.}, language = {en} } @misc{SunkaraRaharinirinaPeppertetal.2020, author = {Sunkara, Vikram and Raharinirina, N. Alexia and Peppert, Felix and von Kleist, Max and Sch{\"u}tte, Christof}, title = {Inferring Gene Regulatory Networks from Single Cell RNA-seq Temporal Snapshot Data Requires Higher Order Moments}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-79664}, year = {2020}, abstract = {Due to the increase in accessibility and robustness of sequencing technology, single cell RNA-seq (scRNA-seq) data has become abundant. The technology has made significant contributions to discovering novel phenotypes and heterogeneities of cells. Recently, there has been a push for using single-- or multiple scRNA-seq snapshots to infer the underlying gene regulatory networks (GRNs) steering the cells' biological functions. To date, this aspiration remains unrealised. In this paper, we took a bottom-up approach and curated a stochastic two gene interaction model capturing the dynamics of a complete system of genes, mRNAs, and proteins. In the model, the regulation was placed upstream from the mRNA on the gene level. We then inferred the underlying regulatory interactions from only the observation of the mRNA population through~time. We could detect signatures of the regulation by combining information of the mean, covariance, and the skewness of the mRNA counts through time. We also saw that reordering the observations using pseudo-time did not conserve the covariance and skewness of the true time course. The underlying GRN could be captured consistently when we fitted the moments up to degree three; however, this required a computationally expensive non-linear least squares minimisation solver. There are still major numerical challenges to overcome for inference of GRNs from scRNA-seq data. These challenges entail finding informative summary statistics of the data which capture the critical regulatory information. Furthermore, the statistics have to evolve linearly or piece-wise linearly through time to achieve computational feasibility and scalability.}, language = {en} } @article{ThiesSunkaraRayetal.2023, author = {Thies, Arne and Sunkara, Vikram and Ray, Sourav and Wulkow, Hanna and Celik, M. {\"O}zg{\"u}r and Yerg{\"o}z, Fatih and Sch{\"u}tte, Christof and Stein, Christoph and Weber, Marcus and Winkelmann, Stefanie}, title = {Modelling altered signalling of G-protein coupled receptors in inflamed environment to advance drug design}, volume = {13}, journal = {Scientific Reports}, number = {607}, doi = {10.1038/s41598-023-27699-w}, year = {2023}, abstract = {We previously reported the successful design, synthesis and testing of the prototype opioid painkiller NFEPP that does not elicit adverse side effects. The design process of NFEPP was based on mathematical modelling of extracellular interactions between G-protein coupled receptors (GPCRs) and ligands, recognizing that GPCRs function differently under pathological versus healthy conditions. We now present an additional and novel stochastic model of GPCR function that includes intracellular dissociation of G-protein subunits and modulation of plasma membrane calcium channels and their dependence on parameters of inflamed and healthy tissue (pH, radicals). The model is validated against in vitro experimental data for the ligands NFEPP and fentanyl at different pH values and radical concentrations. We observe markedly reduced binding affinity and calcium channel inhibition for NFEPP at normal pH compared to lower pH, in contrast to the effect of fentanyl. For increasing radical concentrations, we find enhanced constitutive G-protein activation but reduced ligand binding affinity. Assessing the different effects, the results suggest that, compared to radicals, low pH is a more important determinant of overall GPCR function in an inflamed environment. Future drug design efforts should take this into account.}, language = {en} } @article{VegaSchuetteConrad2016, author = {Vega, Iliusi and Sch{\"u}tte, Christof and Conrad, Tim}, title = {Finding metastable states in real-world time series with recurrence networks}, volume = {445}, journal = {Physica A: Statistical Mechanics and its Applications}, doi = {10.1016/j.physa.2015.10.041}, pages = {1 -- 17}, year = {2016}, abstract = {In the framework of time series analysis with recurrence networks, we introduce a self-adaptive method that determines the elusive recurrence threshold and identifies metastable states in complex real-world time series. As initial step, we introduce a way to set the embedding parameters used to reconstruct the state space from the time series. We set them as the ones giving the maximum Shannon entropy of the diagonal line length distribution for the first simultaneous minima of recurrence rate and Shannon entropy. To identify metastable states, as well as the transitions between them, we use a soft partitioning algorithm for module finding which is specifically developed for the case in which a system shows metastability. We illustrate our method with a complex time series example. Finally, we show the robustness of our method for identifying metastable states. Our results suggest that our method is robust for identifying metastable states in complex time series, even when introducing considerable levels of noise and missing data points.}, language = {en} } @misc{VegaSchuetteConrad2014, author = {Vega, Iliusi and Sch{\"u}tte, Christof and Conrad, Tim}, title = {SAIMeR: Self-adapted method for the identification of metastable states in real-world time series}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-50130}, year = {2014}, abstract = {In the framework of time series analysis with recurrence networks, we introduce SAIMeR, a heuristic self-adapted method that determines the elusive recurrence threshold and identifies metastable states in complex time series. To identify metastable states as well as the transitions between them, we use graph theory concepts and a fuzzy partitioning clustering algorithm. We illustrate SAIMeR by applying it to three real-world time series and show that it is able to identify metastable states in real-world data with noise and missing data points. Finally, we suggest a way to choose the embedding parameters used to construct the state space in which this method is performed, based on the analysis of how the values of these parameters affect two recurrence quantitative measurements: recurrence rate and entropy.}, language = {en} } @article{vonKleistMenzStockeretal.2011, author = {von Kleist, Max and Menz, Stephan and Stocker, Hartmut and Arasteh, Keikawus and Huisinga, Wilhelm and Sch{\"u}tte, Christof}, title = {HIV Quasispecies Dynamics during Pro-active Treatment Switching}, volume = {6}, journal = {Plos One}, number = {3}, publisher = {Public Library of Science}, doi = {10.1371/journal.pone.0018204}, pages = {e18204}, year = {2011}, language = {en} } @article{vonKleistMetznerMarquetetal.2012, author = {von Kleist, Max and Metzner, Ph. and Marquet, R. and Sch{\"u}tte, Christof}, title = {Polymerase Inhibition by Nucleoside Analogs}, volume = {8}, journal = {Plos Computational Biology}, number = {1}, doi = {10.1371/journal.pcbi.1002359}, pages = {e1002359}, year = {2012}, 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} } @inproceedings{WalterSchuette2006, author = {Walter, J. and Sch{\"u}tte, Christof}, title = {Conditional Averaging for Diffusive Fast-Slow Systems}, booktitle = {Analysis, Modeling and Simulation of Multiscale Problems}, editor = {Mielke, A.}, publisher = {Springer}, doi = {10.1007/3-540-35657-6_24}, pages = {647 -- 682}, year = {2006}, language = {en} }