TY - JOUR A1 - Koltai, Peter A1 - Wu, Hao A1 - Noé, Frank A1 - Schütte, Christof T1 - Optimal data-driven estimation of generalized Markov state models for non-equilibrium dynamics JF - Computation Y1 - 2018 U6 - https://doi.org/10.3390/computation6010022 VL - 6 IS - 1 PB - MDPI CY - Basel, Switzerland ER - TY - JOUR A1 - Dibak, Manuel A1 - del Razo, Mauricio J. 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 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. Y1 - 2018 U6 - https://doi.org/10.1063/1.5020294 VL - 148 IS - 21 ER - TY - CHAP A1 - Estacio, Laura A1 - Ehlke, Moritz A1 - Tack, Alexander A1 - Castro-Gutierrez, Eveling A1 - Lamecker, Hans A1 - Mora, Rensso A1 - Zachow, Stefan T1 - Unsupervised Detection of Disturbances in 2D Radiographs T2 - 2021 IEEE 18th International Symposium on Biomedical Imaging (ISBI) N2 - We present a method based on a generative model for detection of disturbances such as prosthesis, screws, zippers, and metals in 2D radiographs. The generative model is trained in an unsupervised fashion using clinical radiographs as well as simulated data, none of which contain disturbances. Our approach employs a latent space consistency loss which has the benefit of identifying similarities, and is enforced to reconstruct X-rays without disturbances. In order to detect images with disturbances, an anomaly score is computed also employing the Frechet distance between the input X-ray and the reconstructed one using our generative model. Validation was performed using clinical pelvis radiographs. We achieved an AUC of 0.77 and 0.83 with clinical and synthetic data, respectively. The results demonstrated a good accuracy of our method for detecting outliers as well as the advantage of utilizing synthetic data. Y1 - 2021 U6 - https://doi.org/10.1109/ISBI48211.2021.9434091 SP - 367 EP - 370 ER - TY - JOUR A1 - Sekuboyina, Anjany A1 - Bayat, Amirhossein A1 - Husseini, Malek E. A1 - Löffler, Maximilian A1 - Li, Hongwei A1 - Tetteh, Giles A1 - Kukačka, Jan A1 - Payer, Christian A1 - Štern, Darko A1 - Urschler, Martin A1 - Chen, Maodong A1 - Cheng, Dalong A1 - Lessmann, Nikolas A1 - Hu, Yujin A1 - Wang, Tianfu A1 - Yang, Dong A1 - Xu, Daguang A1 - Ambellan, Felix A1 - Amiranashvili, Tamaz A1 - Ehlke, Moritz A1 - Lamecker, Hans A1 - Lehnert, Sebastian A1 - Lirio, Marilia A1 - de Olaguer, Nicolás Pérez A1 - Ramm, Heiko A1 - Sahu, Manish A1 - Tack, Alexander A1 - Zachow, Stefan A1 - Jiang, Tao A1 - Ma, Xinjun A1 - Angerman, Christoph A1 - Wang, Xin A1 - Wei, Qingyue A1 - Brown, Kevin A1 - Wolf, Matthias A1 - Kirszenberg, Alexandre A1 - Puybareau, Élodie A1 - Valentinitsch, Alexander A1 - Rempfler, Markus A1 - Menze, Björn H. A1 - Kirschke, Jan S. T1 - VerSe: A Vertebrae Labelling and Segmentation Benchmark for Multi-detector CT Images JF - arXiv Y1 - 2020 ER - TY - JOUR A1 - Wulkow, Hanna A1 - Conrad, Tim A1 - Djurdjevac Conrad, Natasa A1 - Müller, Sebastian A. A1 - Nagel, Kai A1 - Schütte, Christof T1 - Prediction of Covid-19 spreading and optimal coordination of counter-measures: From microscopic to macroscopic models to Pareto fronts JF - PLOS One Y1 - 2021 U6 - https://doi.org/10.1371/journal.pone.0249676 VL - 16 IS - 4 PB - Public Library of Science ER - TY - JOUR A1 - Banisch, Ralf A1 - Djurdjevac Conrad, Natasa A1 - Schütte, Christof T1 - Reactive flows and unproductive cycles for random walks on complex networks JF - The European Physical Journal Special Topics, vol. 224, iss. 12 (2015) pp. 2369-2387 Y1 - 2015 U6 - https://doi.org/10.1140/epjst/e2015-02417-8 ER - TY - GEN A1 - Banisch, Ralf A1 - Djurdjevac Conrad, Natasa A1 - Schütte, Christof T1 - Reactive flows and unproductive cycles for random walks on complex networks N2 - 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. T3 - ZIB-Report - 15-19 KW - Complex networks KW - molecular transition networks KW - transition path theory KW - cycle decomposition KW - reactive trajectories KW - Markow State Methods Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-54239 SN - 1438-0064 ER - TY - JOUR 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 JF - Communications in Mathematical Sciences 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 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. Y1 - 2014 U6 - https://doi.org/10.4310/CMS.2014.v12.n5.a4 VL - 12 IS - 5 SP - 859 EP - 877 ER - TY - JOUR A1 - Djurdjevac Conrad, Natasa A1 - Banisch, Ralf A1 - Schütte, Christof T1 - Modularity of Directed Networks: Cycle Decomposition Approach JF - Journal of Computational Dynamics 2 (2015) pp. 1-24 N2 - 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. Y1 - 2015 U6 - https://doi.org/10.3934/jcd.2015.2.1 ER - TY - GEN A1 - Wang, Han A1 - Schütte, Christof T1 - Building Markov State Models for Periodically Driven Non-Equilibrium Systems N2 - 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. T3 - ZIB-Report - 14-46 KW - Non-equilibrium KW - Markov states model KW - Alanine dipeptide KW - Electric field KW - Floquet theory Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-53167 SN - 1438-0064 N1 - To appear in: Journal of Chemical Theory and Computation ER - TY - GEN A1 - Gul, Raheem A1 - Schütte, Christof A1 - Bernhard, Stefan T1 - Mathematical modeling and sensitivity analysis of arterial anastomosis in arm arteries N2 - 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. T3 - ZIB-Report - 15-22 KW - computational cardiovascular model KW - cardiovascular parameters KW - sensitivity analysis KW - anastomosis KW - Windkessel model Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-54339 SN - 1438-0064 ER - TY - CHAP A1 - Mukhopadhyay, Anirban T1 - Total Variation Random Forest: Fully automatic MRI segmentation in congenital heart disease T2 - RAMBO 2016, HVSMR 2016: Reconstruction, Segmentation, and Analysis of Medical Images Y1 - 2016 U6 - https://doi.org/10.1007/978-3-319-52280-7_17 VL - LNCS 10129 SP - 165 EP - 171 ER - TY - CHAP A1 - Mukhopadhyay, Anirban A1 - Morillo, Oscar A1 - Zachow, Stefan A1 - Lamecker, Hans T1 - Robust and Accurate Appearance Models Based on Joint Dictionary Learning Data from the Osteoarthritis Initiative T2 - Lecture Notes in Computer Science, Patch-Based Techniques in Medical Imaging. Patch-MI 2016 N2 - Deformable model-based approaches to 3D image segmentation have been shown to be highly successful. Such methodology requires an appearance model that drives the deformation of a geometric model to the image data. Appearance models are usually either created heuristically or through supervised learning. Heuristic methods have been shown to work effectively in many applications but are hard to transfer from one application (imaging modality/anatomical structure) to another. On the contrary, supervised learning approaches can learn patterns from a collection of annotated training data. In this work, we show that the supervised joint dictionary learning technique is capable of overcoming the traditional drawbacks of the heuristic approaches. Our evaluation based on two different applications (liver/CT and knee/MR) reveals that our approach generates appearance models, which can be used effectively and efficiently in a deformable model-based segmentation framework. Y1 - 2016 U6 - https://doi.org/10.1007/978-3-319-47118-1_4 VL - 9993 SP - 25 EP - 33 ER - TY - JOUR A1 - Sahu, Manish A1 - Mukhopadhyay, Anirban A1 - Szengel, Angelika A1 - Zachow, Stefan T1 - Addressing multi-label imbalance problem of Surgical Tool Detection using CNN JF - International Journal of Computer Assisted Radiology and Surgery N2 - Purpose: A fully automated surgical tool detection framework is proposed for endoscopic video streams. State-of-the-art surgical tool detection methods rely on supervised one-vs-all or multi-class classification techniques, completely ignoring the co-occurrence relationship of the tools and the associated class imbalance. Methods: In this paper, we formulate tool detection as a multi-label classification task where tool co-occurrences are treated as separate classes. In addition, imbalance on tool co-occurrences is analyzed and stratification techniques are employed to address the imbalance during Convolutional Neural Network (CNN) training. Moreover, temporal smoothing is introduced as an online post-processing step to enhance run time prediction. Results: Quantitative analysis is performed on the M2CAI16 tool detection dataset to highlight the importance of stratification, temporal smoothing and the overall framework for tool detection. Conclusion: The analysis on tool imbalance, backed by the empirical results indicates the need and superiority of the proposed framework over state-of-the-art techniques. Y1 - 2017 UR - https://link.springer.com/article/10.1007/s11548-017-1565-x U6 - https://doi.org/10.1007/s11548-017-1565-x N1 - Selected for final oral presentation VL - 12 IS - 6 SP - 1013 EP - 1020 PB - Springer ER - TY - JOUR A1 - Klus, Stefan A1 - Schütte, Christof T1 - Towards tensor-based methods for the numerical approximation of the Perron-Frobenius and Koopman operator JF - Journal of Computational Dynamics N2 - The global behavior of dynamical systems can be studied by analyzing the eigenvalues and corresponding eigenfunctions of linear operators associated with the system. Two important operators which are frequently used to gain insight into the system's behavior are the Perron-Frobenius operator and the Koopman operator. Due to the curse of dimensionality, computing the eigenfunctions of high-dimensional systems is in general infeasible. We will propose a tensor-based reformulation of two numerical methods for computing finite-dimensional approximations of the aforementioned infinite-dimensional operators, namely Ulam's method and Extended Dynamic Mode Decomposition (EDMD). The aim of the tensor formulation is to approximate the eigenfunctions by low-rank tensors, potentially resulting in a significant reduction of the time and memory required to solve the resulting eigenvalue problems, provided that such a low-rank tensor decomposition exists. Typically, not all variables of a high-dimensional dynamical system contribute equally to the system's behavior, often the dynamics can be decomposed into slow and fast processes, which is also reflected in the eigenfunctions. Thus, the weak coupling between different variables might be approximated by low-rank tensor cores. We will illustrate the efficiency of the tensor-based formulation of Ulam's method and EDMD using simple stochastic differential equations. Y1 - 2016 U6 - https://doi.org/10.3934/jcd.2016007 VL - 3 IS - 2 SP - 139 EP - 161 ER - TY - JOUR A1 - Klus, Stefan A1 - Koltai, Peter A1 - Schütte, Christof T1 - On the numerical approximation of the Perron-Frobenius and Koopman operator JF - Journal of Computational Dynamics N2 - Information about the behavior of dynamical systems can often be obtained by analyzing the eigenvalues and corresponding eigenfunctions of linear operators associated with a dynamical system. Examples of such operators are the Perron-Frobenius and the Koopman operator. In this paper, we will review di� fferent methods that have been developed over the last decades to compute � infinite-dimensional approximations of these in� finite-dimensional operators - in particular Ulam's method and Extended Dynamic Mode Decomposition (EDMD) - and highlight the similarities and di� fferences between these approaches. The results will be illustrated using simple stochastic di� fferential equations and molecular dynamics examples. Y1 - 2016 U6 - https://doi.org/10.3934/jcd.2016003 VL - 3 IS - 1 SP - 51 EP - 77 ER - TY - JOUR A1 - Huttary, Rudolf A1 - Goubergrits, Leonid A1 - Schütte, Christof A1 - Bernhard, Stefan T1 - Simulation, Identification and Statistical Variation in Cardiovascular Analysis (SISCA) - a Software Framework for Multi-compartment Lumped Modeling JF - Computers in Biology and Medicine Y1 - 2017 U6 - https://doi.org/10.1016/j.compbiomed.2017.05.021 VL - 87 SP - 104 EP - 123 ER - TY - GEN A1 - Gupta, Pooja A1 - Gramatke, Annika A1 - Einspanier, Ralf A1 - Schütte, Christof A1 - von Kleist, Max A1 - Sharbati, Jutta T1 - In silicio cytotoxicity assessment on cultured rat intestinal cells deduced from cellular impedance measurements N2 - 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. T3 - ZIB-Report - 17-08 KW - Real-time cell analyzer KW - Toxicity KW - Mathematical modeling KW - IC_{50} Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-62666 SN - 1438-0064 ER - TY - JOUR A1 - Klus, Stefan A1 - Nüske, Feliks A1 - Koltai, Peter A1 - Wu, Hao A1 - Kevrekidis, Ioannis A1 - Schütte, Christof A1 - Noé, Frank T1 - Data-driven model reduction and transfer operator approximation JF - Journal of Nonlinear Science Y1 - 2018 UR - https://link.springer.com/article/10.1007/s00332-017-9437-7 U6 - https://doi.org/10.1007/s00332-017-9437-7 VL - 28 IS - 3 SP - 985 EP - 1010 ER - TY - JOUR A1 - Koltai, Peter A1 - Ciccotti, Giovanni A1 - Schütte, Christof T1 - On Markov state models for non-equilibrium molecular dynamics JF - The Journal of Chemical Physics Y1 - 2016 U6 - https://doi.org/10.1063/1.4966157 N1 - 2016 Editor's Choice of The Journal of Chemical Physics VL - 145 IS - 174103 ER - TY - JOUR A1 - Conrad, Tim A1 - Karsch, K. A1 - Obermeier, Patrick A1 - Seeber, L. A1 - Chen, X. A1 - Tief, Franziska A1 - Muehlhans, S. A1 - Hoppe, Christian A1 - Boettcher, Sindy A1 - Diedrich, S. A1 - Rath, Barbara T1 - Human Parechovirus Infections Associated with Seizures and Rash: A Syndromic Surveillance Study in Children JF - The Pediatric Infectious Disease Journal Y1 - 2015 VL - 34 IS - 10 ER - TY - JOUR A1 - Conrad, Tim A1 - Mireles, Victor T1 - Minimum-overlap clusterings and the sparsity of overcomplete decompositions of binary matrics JF - Procedia Computer Science N2 - Given a set of n binary data points, a widely used technique is to group its features into k clusters. In the case where n {\ensuremath{<}} k, the question of how overlapping are the clusters becomes of interest. In this paper we approach the question through matrix decomposition, and relate the degree of overlap with the sparsity of one of the resulting matrices. We present analytical results regarding bounds on this sparsity, and a heuristic to estimate the minimum amount of overlap that an exact grouping of features into k clusters must have. As shown below, adding new data will not alter this minimum amount of overlap. Y1 - 2015 U6 - https://doi.org/10.1016/j.procs.2015.05.500 VL - 51 SP - 2967 EP - 2971 ER - TY - JOUR A1 - Conrad, Tim A1 - Shao, Borong T1 - Are NoSQL data stores useful for bioinformatics researchers? JF - International Journal on Recent and Innovation Trends in Computing and Communication (IJRITCC) Y1 - 2015 U6 - https://doi.org/10.17762/ijritcc2321-8169.1503176 VL - 3 IS - 3 SP - 1704 EP - 1708 ER - TY - JOUR A1 - Conrad, Tim A1 - Leichtle, Alexander Benedikt A1 - Ceglarek, Uta A1 - Weinert, P. A1 - Nakas, C.T. A1 - Nuoffer, Jean-Marc A1 - Kase, Julia A1 - Witzigmann, Helmut A1 - Thiery, Joachim A1 - Fiedler, Georg Martin T1 - Pancreatic carcinoma, pancreatitis, and healthy controls - metabolite models in a three-class diagnostic dilemma JF - Metabolomics N2 - Background: Metabolomics as one of the most rapidly growing technologies in the ?-omics?field denotes the comprehensive analysis of low molecular-weight compounds and their pathways. Cancer-specific alterations of the metabolome can be detected by high-throughput massspectrometric metabolite profiling and serve as a considerable source of new markers for the early differentiation of malignant diseases as well as their distinction from benign states. However, a comprehensive framework for the statistical evaluation of marker panels in a multi-class setting has not yet been established. Methods: We collected serum samples of 40 pancreatic carcinoma patients, 40 controls, and 23 pancreatitis patients according to standard protocols and generated amino acid profiles by routine mass-spectrometry. In an intrinsic three-class bioinformatic approach we compared these profiles, evaluated their selectivity and computed multi-marker panels combined with the conventional tumor marker CA 19-9. Additionally, we tested for non-inferiority and superiority to determine the diagnostic surplus value of our multi-metabolite marker panels.  Results: Compared to CA 19-9 alone, the combined amino acid-based metabolite panel had a superior selectivity for the discrimination of healthy controls, pancreatitis, and pancreatic carcinoma patients [Volume under ROC surface (VUS) = 0.891 (95\% CI 0.794 - 0.968)]. Conclusions: We combined highly standardized samples, a three-class study design, a highthroughput mass-spectrometric technique, and a comprehensive bioinformatic framework to identify metabolite panels selective for all three groups in a single approach. Our results suggest that metabolomic profiling necessitates appropriate evaluation strategies and ?despite all its current limitations? can deliver marker panels with high selectivity even in multi-class settings. Y1 - 2013 U6 - https://doi.org/10.1007/s11306-012-0476-7 VL - 9 IS - 3 SP - 677 EP - 687 ER - TY - JOUR A1 - Gupta, Pooja A1 - Reinsch, Norbert A1 - Spötter, Andreas A1 - Conrad, Tim A1 - Bienefeld, Kaspar T1 - Accuracy of the unified approach in maternally influenced traits - illustrated by a simulation study in the honey bee (Apis mellifera) JF - BMC Genetics Y1 - 2013 U6 - https://doi.org/10.1186/1471-2156-14-36 VL - 14 IS - 36 ER - TY - JOUR A1 - Conrad, Tim A1 - You, Xintian T1 - Acfs: accurate circRNA identification and quantification from NGS data JF - Nature Scientific Reports N2 - Circular RNAs (circRNAs) are a group of single-stranded RNAs in closed circular form. They are splicing-generated, widely expressed in various tissues and have functional implications in development and diseases. To facilitate genome-wide characterization of circRNAs using RNA-Seq data, we present a freely available software package named acfs. Acfs allows de novo, accurate and fast identification and abundance quantification of circRNAs from single- and paired-ended RNA-Seq data. On simulated datasets, acfs achieved the highest F1 accuracy and lowest false discovery rate among current state-of-the-art tools. On real-world datasets, acfs efficiently identified more bona fide circRNAs. Furthermore, we demonstrated the power of circRNA analysis on two leukemia datasets. We identified a set of circRNAs that are differentially expressed between AML and APL samples, which might shed light on the potential molecular classification of complex diseases using circRNA profiles. Moreover, chromosomal translocation, as manifested in numerous diseases, could produce not only fusion transcripts but also fusion circRNAs of clinical relevance. Featured with high accuracy, low FDR and the ability to identify fusion circRNAs, we believe that acfs is well suited for a wide spectrum of applications in characterizing the landscape of circRNAs from non-model organisms to cancer biology. Y1 - 2016 U6 - https://doi.org/10.1038/srep38820 VL - 6 ER - TY - JOUR A1 - Hoppe, Christian A1 - Obermeier, Patrick A1 - Mehlhans, S. A1 - Alchikh, Maren A1 - Seeber, L. A1 - Tief, Franziska A1 - Karsch, K. A1 - Chen, X. A1 - Boettcher, Sindy A1 - Diedrich, S. A1 - Conrad, Tim T1 - Innovative Digital Tools and Surveillance Systems for the Timely Detection of Adverse Events at the Point of Care: A Proof-of-Concept Study JF - Drug Safety N2 - Regulatory authorities often receive poorly structured safety reports requiring considerable effort to investigate potential adverse events post hoc. Automated question-and-answer systems may help to improve the overall quality of safety information transmitted to pharmacovigilance agencies. This paper explores the use of the VACC-Tool (ViVI Automated Case Classification Tool) 2.0, a mobile application enabling physicians to classify clinical cases according to 14 pre-defined case definitions for neuroinflammatory adverse events (NIAE) and in full compliance with data standards issued by the Clinical Data Interchange Standards Consortium. METHODS: The validation of the VACC-Tool 2.0 (beta-version) was conducted in the context of a unique quality management program for children with suspected NIAE in collaboration with the Robert Koch Institute in Berlin, Germany. The VACC-Tool was used for instant case classification and for longitudinal follow-up throughout the course of hospitalization. Results were compared to International Classification of Diseases , Tenth Revision (ICD-10) codes assigned in the emergency department (ED). RESULTS: From 07/2013 to 10/2014, a total of 34,368 patients were seen in the ED, and 5243 patients were hospitalized; 243 of these were admitted for suspected NIAE (mean age: 8.5 years), thus participating in the quality management program. Using the VACC-Tool in the ED, 209 cases were classified successfully, 69 \% of which had been missed or miscoded in the ED reports. Longitudinal follow-up with the VACC-Tool identified additional NIAE. CONCLUSION: Mobile applications are taking data standards to the point of care, enabling clinicians to ascertain potential adverse events in the ED setting and during inpatient follow-up. Compliance with Clinical Data Interchange Standards Consortium (CDISC) data standards facilitates data interoperability according to regulatory requirements. Y1 - 2016 U6 - https://doi.org/10.1007/s40264-016-0437-6 VL - 39 IS - 10 SP - 977 EP - 988 ER - TY - CHAP A1 - Shao, Borong A1 - Conrad, Tim T1 - Epithelial Mesenchymal Transition Regulatory Network-based Feature Selection in Lung Cancer Prognosis Prediction T2 - Lecture Notes in Computer Science (LNCS) N2 - Feature selection technique is often applied in identifying cancer prognosis biomarkers. However, many feature selection methods are prone to over-fitting or poor biological interpretation when applied on biological high-dimensional data. Network-based feature selection and data integration approaches are proposed to identify more robust biomarkers. We conducted experiments to investigate the advantages of the two approaches using epithelial mesenchymal transition regulatory network, which is demonstrated as highly relevant to cancer prognosis. We obtained data from The Cancer Genome Atlas. Prognosis prediction was made using Support Vector Machine. Under our experimental settings, the results showed that network-based features gave significantly more accurate predictions than individual molecular features, and features selected from integrated data (RNA-Seq and micro-RNA data) gave significantly more accurate predictions than features selected from single source data (RNA-Seq data). Our study indicated that biological network-based feature transformation and data integration are two useful approaches to identify robust cancer biomarkers. Y1 - 2016 U6 - https://doi.org/10.1007/978-3-319-31744-1_13 VL - 9656 SP - 1235 EP - 146 ER - TY - JOUR A1 - Tief, Franziska A1 - Hoppe, Christian A1 - Seeber, L. A1 - Obermeier, Patrick A1 - Chen, X. A1 - Karsch, K. A1 - Muehlhans, S. A1 - Adamou, E. A1 - Conrad, Tim A1 - Schweiger, Brunhilde A1 - Adam, T. A1 - Rath, Barbara T1 - An inception cohort study assessing the role of bacterial co-infections in children with influenza and ILI and a clinical decision model for stringent antibiotic use JF - Antiviral Therapy N2 - BACKGROUND: Influenza-like illness (ILI) is a common reason for paediatric consultations. Viral causes predominate, but antibiotics are used frequently. With regard to influenza, pneumococcal coinfections are considered major contributors to morbidity/mortality. METHODS: In the context of a perennial quality management (QM) programme at the Charit{\'e} Departments of Paediatrics and Microbiology in collaboration with the Robert Koch Institute, children aged 0-18 years presenting with signs and symptoms of ILI were followed from the time of initial presentation until hospital discharge (Charit{\'e} Influenza-Like Disease = ChILD Cohort). An independent QM team performed highly standardized clinical assessments using a disease severity score based on World Health Organization criteria for uncomplicated and complicated/progressive disease. Nasopharyngeal and pharyngeal samples were collected for viral reverse transcription polymerase chain reaction and bacterial culture/sensitivity and MaldiTOF analyses. The term 'detection' was used to denote any evidence of viral or bacterial pathogens in the (naso)pharyngeal cavity. With the ChILD Cohort data collected, a standard operating procedure (SOP) was created as a model system to reduce the inappropriate use of antibiotics in children with ILI. Monte Carlo simulations were performed to assess cost-effectiveness. RESULTS: Among 2,569 ChILD Cohort patients enrolled from 12/2010 to 04/2013 (55\% male, mean age 3.2 years, range 0-18, 19\% {\ensuremath{>}}5 years), 411 patients showed laboratory-confirmed influenza, with bacterial co-detection in 35\%. Influenza and pneumococcus were detected simultaneously in 12/2,569 patients, with disease severity clearly below average. Pneumococcal vaccination rates were close to 90\%. Nonetheless, every fifth patient was already on antibiotics upon presentation; new antibiotic prescriptions were issued in an additional 20\%. Simulation of the model SOP in the same dataset revealed that the proposed decision model could have reduced the inappropriate use of antibiotics significantly (P{\ensuremath{<}}0.01) with an incremental cost-effectiveness ratio of -99.55?. CONCLUSIONS: Physicians should be made aware that in times of pneumococcal vaccination the prevalence and severity of influenza infections complicated by pneumococci may decline. Microbiological testing in combination with standardized disease severity assessments and review of vaccination records could be cost-effective, as well as promoting stringent use of antibiotics and a personalized approach to managing children with ILI. Y1 - 2016 U6 - https://doi.org/10.3851/IMP3034 VL - 21 SP - 413 EP - 424 ER - TY - JOUR A1 - Obermeier, Patrick A1 - Muehlhans, S. A1 - Hoppe, Christian A1 - Karsch, K. A1 - Tief, Franziska A1 - Seeber, L. A1 - Chen, X. A1 - Conrad, Tim A1 - Boettcher, Sindy A1 - Diedrich, S. A1 - Rath, Barbara T1 - Enabling Precision Medicine With Digital Case Classification at the Point-of-Care JF - EBioMedicine N2 - Infectious and inflammatory diseases of the central nervous system are difficult to identify early. Case definitions for aseptic meningitis, encephalitis, myelitis, and acute disseminated encephalomyelitis (ADEM) are available, but rarely put to use. The VACC-Tool (Vienna Vaccine Safety Initiative Automated Case Classification-Tool) is a mobile application enabling immediate case ascertainment based on consensus criteria at the point-of-care. The VACC-Tool was validated in a quality management program in collaboration with the Robert-Koch-Institute. Results were compared to ICD-10 coding and retrospective analysis of electronic health records using the same case criteria. Of 68,921 patients attending the emergency room in 10/2010-06/2013, 11,575 were hospitalized, with 521 eligible patients (mean age: 7.6 years) entering the quality management program. Using the VACC-Tool at the point-of-care, 180/521 cases were classified successfully and 194/521 ruled out with certainty. Of the 180 confirmed cases, 116 had been missed by ICD-10 coding, 38 misclassified. By retrospective application of the same case criteria, 33 cases were missed. Encephalitis and ADEM cases were most likely missed or misclassified. The VACC-Tool enables physicians to ask the right questions at the right time, thereby classifying cases consistently and accurately, facilitating translational research. Future applications will alert physicians when additional diagnostic procedures are required. Y1 - 2016 U6 - https://doi.org/10.1016/j.ebiom.2016.01.008 VL - 4 SP - 191 EP - 196 ER - TY - JOUR A1 - Conrad, Tim A1 - Bruckner, Sharon A1 - Kayser, Bastian T1 - Finding Modules in Networks with Non-modular Regions JF - Lecture Notes in Computer Science (Proceedings of SEA 2013) N2 - Most network clustering methods share the assumption that the network can be completely decomposed into modules, that is, every node belongs to (usually exactly one) module. Forcing this constraint can lead to misidentification of modules where none exist, while the true modules are drowned out in the noise, as has been observed e.g. for protein interaction networks. We thus propose a clustering model where networks contain both a modular region consisting of nodes that can be partitioned into modules, and a transition region containing nodes that lie between or outside modules. We propose two scores based on spectral properties to determine how well a network fits this model. We then evaluate three (partially adapted) clustering algorithms from the literature on random networks that fit our model, based on the scores and comparison to the ground truth. This allows to pinpoint the types of networks for which the different algorithms perform well. Y1 - 2013 U6 - https://doi.org/10.1007/978-3-642-38527-8_18 VL - 7933 SP - 188 EP - 199 ER - TY - JOUR A1 - Conrad, Tim A1 - Rath, Barbara A1 - Tief, Franziska A1 - Karsch, K. A1 - Muehlhans, S. A1 - Obermeier, Patrick A1 - Adamou, E. A1 - Chen, X. A1 - Seeber, L. A1 - Peiser, Ch. A1 - Hoppe, Christian A1 - von Kleist, Max A1 - Schweiger, Brunhilde T1 - Towards a personalized approach to managing of influenza infections in infants and children - food for thought and a note on oseltamivir JF - Infectious Disorders - Drug Targets Y1 - 2013 VL - 13 IS - 1 SP - 25 EP - 33 ER - TY - JOUR A1 - Conrad, Tim A1 - Leichtle, Alexander Benedikt A1 - Nuoffer, Jean-Marc A1 - Ceglarek, Uta A1 - Kase, Julia A1 - Witzigmann, Helmut A1 - Thiery, Joachim A1 - Fiedler, Georg Martin T1 - Serum amino acid profiles and their alterations in colorectal cancer JF - Metabolomics N2 - Mass spectrometry-based serum metabolic profiling is a promising tool to analyse complex cancer associated metabolic alterations, which may broaden our pathophysiological understanding of the disease and may function as a source of new cancer-associated biomarkers. Highly standardized serum samples of patients suffering from colon cancer (n = 59) and controls (n = 58) were collected at the University Hospital Leipzig. We based our investigations on amino acid screening profiles using electrospray tandem-mass spectrometry. Metabolic profiles were evaluated using the Analyst 1.4.2 software. General, comparative and equivalence statistics were performed by R 2.12.2. 11 out of 26 serum amino acid concentrations were significantly different between colorectal cancer patients and healthy controls. We found a model including CEA, glycine, and tyrosine as best discriminating and superior to CEA alone with an AUROC of 0.878 (95\% CI 0.815?0.941). Our serum metabolic profiling in colon cancer revealed multiple significant disease-associated alterations in the amino acid profile with promising diagnostic power. Further large-scale studies are necessary to elucidate the potential of our model also to discriminate between cancer and potential differential diagnoses. In conclusion, serum glycine and tyrosine in combination with CEA are superior to CEA for the discrimination between colorectal cancer patients and controls. Y1 - 2012 U6 - https://doi.org/10.1007/s11306-011-0357-5 ER - TY - JOUR A1 - Gupta, Pooja A1 - Conrad, Tim A1 - Spötter, Andreas A1 - Reinsch, Norbert A1 - Bienefeld, Kaspar T1 - Simulating a base population in honey bee for molecular genetic studies JF - Genetics Selection Evolution N2 - Over the past years, reports have indicated that honey bee populations are declining and that infestation by an ecto-parasitic mite (Varroa destructor) is one of the main causes. Selective breeding of resistant bees can help to prevent losses due to the parasite, but it requires that a robust breeding program and genetic evaluation are implemented. Genomic selection has emerged as an important tool in animal breeding programs and simulation studies have shown that it yields more accurate breeding values estimates, higher genetic gain and low rates of inbreeding. Since genomic selection relies on marker data, simulations conducted on a genomic dataset are a pre-requisite before selection can be implemented. Although genomic datasets have been simulated in other species undergoing genetic evaluation, simulation of a genomic dataset specific to the honey bee is required since this species has distinct genetic and reproductive biology characteristics. Our software program was aimed at constructing a base population by simulating a random mating honey bee population. A forward-time population simulation approach was applied since it allows modeling of genetic characteristics and reproductive behavior specific to the honey bee.  Results: Our software program yielded a genomic dataset for a base population in linkage disequilibrium. In addition, information was obtained on (1) the position of markers on each chromosome, (2) allele frequency, (3) ?2 statistics for Hardy- Weinberg equilibrium, (4) a sorted list of markers with a minor allele frequency less than or equal to the input value, (5) average r2 values of linkage disequilibrium between all simulated marker loci pair for all generations and (6) average r2 value of linkage disequilibrium in the last generation for selected markers with the highest minor allele frequency. Conclusion: We developed a software program that takes into account the genetic and reproductive biology characteristics specific to the honey bee and that can be used to constitute a genomic dataset compatible with the simulation studies necessary to optimize breeding programs. The source code together with an instruction file is freely accessible at http://msproteomics.org/Research/Misc/honeybeepopulationsimulator.html Y1 - 2012 U6 - https://doi.org/10.1186/1297-9686-44-14 ER - TY - JOUR A1 - Conrad, Tim T1 - New Appraches for Visualizing and Analyzing Metabolic Pathways JF - Proceedings of the Second Australian Undergraduate Students? Computing Conference N2 - Visualizing of metabolic pathways (or networks) has been done by many differentapproaches. In this work, we implemented and tested existing graph layout algorithms, and present a new approach to lay-out medium size metabolic pathways (500-20,000 vertices) by implementing and combining three well known graph lay-out algorithms (high dimension embedding, spring-embedder preprocessing, spring-embedder), through 3D space density analysis facilitated by the Octree technique. For the analysis of the results of metabolic pathways simulations we present two new techniques: rstly, a powerful technique to visualize pathways simulation data was created to unveil and understand concentration ows through metabolic pathways. This was achieved by mapping the color encoded concentration value of every substance from each time step of the simulation to its graphical representation in the layout. By combining all resulting images (from each time step) and displaying them as a movie, many characteristics such as subnetworks, alternative routes through the network, and differences between a modied pathway and its unmodied version can be revealed. Secondly, a new method to detect co-regulated substances in metabolic pathways and to recognize differences between two versions of a pathway, was established. To do this, we transformed the simulation data into a row-based representation, color-coded these rows, and reordered them with respect to similarity by using a Genetic Algorithm variant. From the arising discrete 2-dimensional matrix consisting of concentration values, a continuous 2-dimensional fourier row function was computed. This function can be used to measure properties, such as similarities in a pathway between time steps, or substances, or to detect and evaluate differences between modied versions of the same pathway. Y1 - 2004 ER - TY - JOUR A1 - Gelß, Patrick A1 - Matera, Sebastian A1 - Schütte, Christof T1 - Solving the master equation without kinetic Monte Carlo: Tensor train approximations for a CO oxidation model JF - Journal of Computational Physics N2 - In multiscale modeling of heterogeneous catalytic processes, one crucial point is the solution of a Markovian master equation describing the stochastic reaction kinetics. Usually, this is too high-dimensional to be solved with standard numerical techniques and one has to rely on sampling approaches based on the kinetic Monte Carlo method. In this study we break the curse of dimensionality for the direct solution of the Markovian master equation by exploiting the Tensor Train Format for this purpose. The performance of the approach is demonstrated on a first principles based, reduced model for the CO oxidation on the RuO2(110) surface. We investigate the complexity for increasing system size and for various reaction conditions. The advantage over the stochastic simulation approach is illustrated by a problem with increased Y1 - 2016 U6 - https://doi.org/10.1016/j.jcp.2016.03.025 VL - 314 SP - 489 EP - 502 ER - TY - JOUR A1 - Vega, Iliusi A1 - Schütte, Christof A1 - Conrad, Tim T1 - Finding metastable states in real-world time series with recurrence networks JF - Physica A: Statistical Mechanics and its Applications N2 - 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. Y1 - 2016 U6 - https://doi.org/10.1016/j.physa.2015.10.041 VL - 445 SP - 1 EP - 17 ER - TY - JOUR A1 - Conrad, Tim A1 - Genzel, Martin A1 - Cvetkovic, Nada A1 - Wulkow, Niklas A1 - Leichtle, Alexander Benedikt A1 - Vybiral, Jan A1 - Kytyniok, Gitta A1 - Schütte, Christof T1 - Sparse Proteomics Analysis – a compressed sensing-based approach for feature selection and classification of high-dimensional proteomics mass spectrometry data JF - BMC Bioinfomatics N2 - Background: 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 in how mass spectra differ 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 against noise and outliers, while the identified feature set should be as small as possible. Results: We present a new algorithm, Sparse Proteomics Analysis (SPA),based on thet heory of compressed sensing that allows us to identify a minimal discriminating set of features from mass spectrometry data-sets. We show (1) how our method performs on artificial and real-world data-sets, (2) that its performance is competitive with standard (and widely used) algorithms for analyzing proteomics data, and (3) that it is robust against random and systematic noise. We further demonstrate the applicability of our algorithm to two previously published clinical data-sets. Y1 - 2017 U6 - https://doi.org/10.1186/s12859-017-1565-4 VL - 18 IS - 160 ER - TY - JOUR A1 - Gupta, Pooja A1 - Gramatke, Annika A1 - Einspanier, Ralf A1 - Schütte, Christof A1 - von Kleist, Max A1 - Sharbati, Jutta T1 - In silico cytotoxicity assessment on cultured rat intestinal cells deduced from cellular impedance measurements JF - Toxicology in Vitro N2 - 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 IC50 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. IC50 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. KW - Real-time cell analyzer KW - Toxicity KW - Mathematical modeling Y1 - 2017 SN - 1438-0064 VL - 41 SP - 179 EP - 188 ER - TY - JOUR A1 - Wilson, David A1 - Anglin, Carolyn A1 - Ambellan, Felix A1 - Grewe, Carl Martin A1 - Tack, Alexander A1 - Lamecker, Hans A1 - Dunbar, Michael A1 - Zachow, Stefan T1 - Validation of three-dimensional models of the distal femur created from surgical navigation point cloud data for intraoperative and postoperative analysis of total knee arthroplasty JF - International Journal of Computer Assisted Radiology and Surgery N2 - Purpose: Despite the success of total knee arthroplasty there continues to be a significant proportion of patients who are dissatisfied. One explanation may be a shape mismatch between pre and post-operative distal femurs. The purpose of this study was to investigate a method to match a statistical shape model (SSM) to intra-operatively acquired point cloud data from a surgical navigation system, and to validate it against the pre-operative magnetic resonance imaging (MRI) data from the same patients. Methods: A total of 10 patients who underwent navigated total knee arthroplasty also had an MRI scan less than 2 months pre-operatively. The standard surgical protocol was followed which included partial digitization of the distal femur. Two different methods were employed to fit the SSM to the digitized point cloud data, based on (1) Iterative Closest Points (ICP) and (2) Gaussian Mixture Models (GMM). The available MRI data were manually segmented and the reconstructed three-dimensional surfaces used as ground truth against which the statistical shape model fit was compared. Results: For both approaches, the difference between the statistical shape model-generated femur and the surface generated from MRI segmentation averaged less than 1.7 mm, with maximum errors occurring in less clinically important areas. Conclusion: The results demonstrated good correspondence with the distal femoral morphology even in cases of sparse data sets. Application of this technique will allow for measurement of mismatch between pre and post-operative femurs retrospectively on any case done using the surgical navigation system and could be integrated into the surgical navigation unit to provide real-time feedback. Y1 - 2017 UR - https://link.springer.com/content/pdf/10.1007%2Fs11548-017-1630-5.pdf U6 - https://doi.org/10.1007/s11548-017-1630-5 VL - 12 IS - 12 SP - 2097 EP - 2105 PB - Springer ER - TY - JOUR A1 - von Tycowicz, Christoph A1 - Ambellan, Felix A1 - Mukhopadhyay, Anirban A1 - Zachow, Stefan T1 - An Efficient Riemannian Statistical Shape Model using Differential Coordinates JF - Medical Image Analysis N2 - We propose a novel Riemannian framework for statistical analysis of shapes that is able to account for the nonlinearity in shape variation. By adopting a physical perspective, we introduce a differential representation that puts the local geometric variability into focus. We model these differential coordinates as elements of a Lie group thereby endowing our shape space with a non-Euclidean structure. A key advantage of our framework is that statistics in a manifold shape space becomes numerically tractable improving performance by several orders of magnitude over state-of-the-art. We show that our Riemannian model is well suited for the identification of intra-population variability as well as inter-population differences. In particular, we demonstrate the superiority of the proposed model in experiments on specificity and generalization ability. We further derive a statistical shape descriptor that outperforms the standard Euclidean approach in terms of shape-based classification of morphological disorders. Y1 - 2018 U6 - https://doi.org/10.1016/j.media.2017.09.004 VL - 43 IS - 1 SP - 1 EP - 9 ER - TY - JOUR A1 - Benn, Andreas A1 - Hiepen, Christian A1 - Osterland, Marc A1 - Schütte, Christof A1 - Zwijsen, An A1 - Knaus, Petra T1 - Role of bone morphogenetic proteins in sprouting angiogenesis: differential BMP receptor-dependent signaling pathways balance stalk vs. tip cell competence JF - FASEB Journal N2 - 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. Y1 - 2017 U6 - https://doi.org/10.1096/fj.201700193RR VL - 31 IS - 11 SP - 4720 EP - 4733 ER - TY - JOUR A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - Hybrid models for chemical reaction networks: Multiscale theory and application to gene regulatory systems JF - The Journal of Chemical Physics N2 - Well-mixed stochastic chemical kinetics are properly modeled 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 expressing 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. Especially, we reveal the cause of error in the case of small volume approximations. Y1 - 2017 U6 - https://doi.org/10.1063/1.4986560 VL - 147 IS - 11 SP - 114115-1 EP - 114115-18 ER - TY - GEN A1 - Schütte, Christof A1 - Conrad, Tim ED - Deuflhard, Peter ED - Grötschel, Martin ED - Hömberg, Dietmar ED - Horst, Ulrich ED - Kramer, Jürg ED - Mehrmann, Volker ED - Polthier, Konrad ED - Schmidt, Frank ED - Skutella, Martin ED - Sprekels, Jürgen T1 - Showcase 3: Information-based medicine T2 - MATHEON-Mathematics for Key Technologies Y1 - 2014 VL - 1 SP - 66 EP - 67 PB - European Mathematical Society ER - TY - JOUR A1 - Kryven, Ivan A1 - Röblitz, Susanna A1 - Schütte, Christof T1 - Solution of the chemical master equation by radial basis functions approximation with interface tracking JF - BMC Systems Biology N2 - Background. The chemical master equation is the fundamental equation of stochastic chemical kinetics. This differential-difference equation describes temporal evolution of the probability density function for states of a chemical system. A state of the system, usually encoded as a vector, represents the number of entities or copy numbers of interacting species, which are changing according to a list of possible reactions. It is often the case, especially when the state vector is high-dimensional, that the number of possible states the system may occupy is too large to be handled computationally. One way to get around this problem is to consider only those states that are associated with probabilities that are greater than a certain threshold level. Results. We introduce an algorithm that significantly reduces computational resources and is especially powerful when dealing with multi-modal distributions. The algorithm is built according to two key principles. Firstly, when performing time integration, the algorithm keeps track of the subset of states with significant probabilities (essential support). Secondly, the probability distribution that solves the equation is parametrised with a small number of coefficients using collocation on Gaussian radial basis functions. The system of basis functions is chosen in such a way that the solution is approximated only on the essential support instead of the whole state space. Discussion. In order to demonstrate the effectiveness of the method, we consider four application examples: a) the self-regulating gene model, b) the 2-dimensional bistable toggle switch, c) a generalisation of the bistable switch to a 3-dimensional tristable problem, and d) a 3-dimensional cell differentiation model that, depending on parameter values, may operate in bistable or tristable modes. In all multidimensional examples the manifold containing the system states with significant probabilities undergoes drastic transformations over time. This fact makes the examples especially challenging for numerical methods. Conclusions. The proposed method is a new numerical approach permitting to approximately solve a wide range of problems that have been hard to tackle until now. A full representation of multi-dimensional distributions is recovered. The method is especially attractive when dealing with models that yield solutions of a complex structure, for instance, featuring multi-stability. Electronic version: http://www.biomedcentral.com/1752-0509/9/67 Y1 - 2015 U6 - https://doi.org/10.1186/s12918-015-0210-y VL - 9 IS - 67 SP - 1 EP - 12 ER - TY - CHAP A1 - Lamecker, Hans A1 - Zachow, Stefan T1 - Statistical Shape Modeling of Musculoskeletal Structures and Its Applications T2 - Computational Radiology for Orthopaedic Interventions N2 - Statistical shape models (SSM) describe the shape variability contained in a given population. They are able to describe large populations of complex shapes with few degrees of freedom. This makes them a useful tool for a variety of tasks that arise in computer-aided madicine. In this chapter we are going to explain the basic methodology of SSMs and present a variety of examples, where SSMs have been successfully applied. Y1 - 2016 SN - 978-3-319-23481-6 U6 - https://doi.org/10.1007/978-3-319-23482-3 VL - 23 SP - 1 EP - 23 PB - Springer ER - TY - CHAP A1 - Mukhopadhyay, Anirban A1 - Oksuz, Ilkay A1 - Bevilacqua, Marco A1 - Dharmakumar, Rohan A1 - Tsaftaris, Sotirios T1 - Data-Driven Feature Learning for Myocardial Segmentation of CP-BOLD MRI T2 - Functional Imaging and Modeling of the Heart N2 - Cardiac Phase-resolved Blood Oxygen-Level-Dependent (CP- BOLD) MR is capable of diagnosing an ongoing ischemia by detecting changes in myocardial intensity patterns at rest without any contrast and stress agents. Visualizing and detecting these changes require significant post-processing, including myocardial segmentation for isolating the myocardium. But, changes in myocardial intensity pattern and myocardial shape due to the heart’s motion challenge automated standard CINE MR myocardial segmentation techniques resulting in a significant drop of segmentation accuracy. We hypothesize that the main reason behind this phenomenon is the lack of discernible features. In this paper, a multi scale discriminative dictionary learning approach is proposed for supervised learning and sparse representation of the myocardium, to improve the myocardial feature selection. The technique is validated on a challenging dataset of CP-BOLD MR and standard CINE MR acquired in baseline and ischemic condition across 10 canine subjects. The proposed method significantly outperforms standard cardiac segmentation techniques, including segmentation via registration, level sets and supervised methods for myocardial segmentation. Y1 - 2015 U6 - https://doi.org/10.1007/978-3-319-20309-6_22 VL - 9126 SP - 189 EP - 197 PB - Springer ER - TY - CHAP A1 - Mukhopadhyay, Anirban A1 - Oksuz, Ilkay A1 - Bevilacqua, Marco A1 - Dharmakumar, Rohan A1 - Tsaftaris, Sotirios T1 - Unsupervised myocardial segmentation for cardiac MRI T2 - Medical Image Computing and Computer-Assisted Intervention -- MICCAI 2015 N2 - Though unsupervised segmentation was a de-facto standard for cardiac MRI segmentation early on, recently cardiac MRI segmentation literature has favored fully supervised techniques such as Dictionary Learning and Atlas-based techniques. But, the benefits of unsupervised techniques e.g., no need for large amount of training data and better potential of handling variability in anatomy and image contrast, is more evident with emerging cardiac MR modalities. For example, CP-BOLD is a new MRI technique that has been shown to detect ischemia without any contrast at stress but also at rest conditions. Although CP-BOLD looks similar to standard CINE, changes in myocardial intensity patterns and shape across cardiac phases, due to the heart’s motion, BOLD effect and artifacts affect the underlying mechanisms of fully supervised segmentation techniques resulting in a significant drop in segmentation accuracy. In this paper, we present a fully unsupervised technique for segmenting myocardium from the background in both standard CINE MR and CP-BOLD MR. We combine appearance with motion information (obtained via Optical Flow) in a dictionary learning framework to sparsely represent important features in a low dimensional space and separate myocardium from background accordingly. Our fully automated method learns background-only models and one class classifier provides myocardial segmentation. The advantages of the proposed technique are demonstrated on a dataset containing CP-BOLD MR and standard CINE MR image sequences acquired in baseline and ischemic condition across 10 canine subjects, where our method outperforms state-of-the-art supervised segmentation techniques in CP-BOLD MR and performs at-par for standard CINE MR. Y1 - 2015 U6 - https://doi.org/10.1007/978-3-319-24574-4_2 VL - LNCS 9351 SP - 12 EP - 20 ER - TY - CHAP A1 - Oksuz, Ilkay A1 - Mukhopadhyay, Anirban A1 - Bevilacqua, Marco A1 - Dharmakumar, Rohan A1 - Tsaftaris, Sotirios T1 - Dictionary Learning Based Image Descriptor for Myocardial Registration of CP-BOLD MR T2 - Medical Image Computing and Computer-Assisted Intervention -- MICCAI 2015 N2 - Cardiac Phase-resolved Blood Oxygen-Level-Dependent (CP- BOLD) MRI is a new contrast agent- and stress-free imaging technique for the assessment of myocardial ischemia at rest. The precise registration among the cardiac phases in this cine type acquisition is essential for automating the analysis of images of this technique, since it can potentially lead to better specificity of ischemia detection. However, inconsistency in myocardial intensity patterns and the changes in myocardial shape due to the heart’s motion lead to low registration performance for state- of-the-art methods. This low accuracy can be explained by the lack of distinguishable features in CP-BOLD and inappropriate metric defini- tions in current intensity-based registration frameworks. In this paper, the sparse representations, which are defined by a discriminative dictionary learning approach for source and target images, are used to improve myocardial registration. This method combines appearance with Gabor and HOG features in a dictionary learning framework to sparsely represent features in a low dimensional space. The sum of squared differences of these distinctive sparse representations are used to define a similarity term in the registration framework. The proposed descriptor is validated on a challenging dataset of CP-BOLD MR and standard CINE MR acquired in baseline and ischemic condition across 10 canines. Y1 - 2015 U6 - https://doi.org/10.1007/978-3-319-24571-3_25 VL - 9350 SP - 205 EP - 213 PB - Springer ER - TY - GEN A1 - Osterland, Marc A1 - Benn, Andreas A1 - Prohaska, Steffen A1 - Schütte, Christof T1 - Single Cell Tracking in Phase-Contrast Microscopy T2 - EMBL Symposium 2015 - Seeing is Believing - Imaging the Processes of Life N2 - In this work, we developed an automatic algorithm to analyze cell migration in chemotaxis assays, based on phase-contrast time-lapse microscopy. While manual approaches are still widely used in recent publications, our algorithm is able to track hundreds of single cells per frame. The extracted paths are analysed with traditional geometrical approaches as well as diffusion-driven Markov state models (MSM). Based on these models, a detailed view on spatial and temporal effects is possible. Using our new approach on experimental data, we are able to distinguish between directed migration (e.g. towards a VEGF gradient) and random migration without favored direction. A calculation of the committor probabilities reveals that cells of the whole image area are more likely to migrate directly towards the VEGF than away from it during the first four hours. However, in absence of a chemoattractant, cells migrate more likely to their nearest image border. These conclusions are supported by the spatial mean directions. In a next step, the cell-cell interaction during migration and the migration of cell clusters will be analyzed. Furthermore, we want to observe phenotypical changes during migration based on fluorescence microscopy and machine learning. The algorithm is part of a collaborative platform which brings the experimental expertise of scientists from life sciences and the analytical knowledge of computer scientists together. This platform is built using web-based technologies with a responsive real-time user interface. All data, including raw and metadata as well as the accompanying results, will be stored in a secure and scalable compute cluster. The compute cluster provides sufficient space and computational power for modern image-based experiments and their analyses. Specific versions of data and results can be tagged to keep immutable records for archival. Y1 - 2015 ER - TY - GEN A1 - Koltai, Peter A1 - Ciccotti, Giovanni A1 - Schütte, Christof T1 - On metastability and Markov state models for non-stationary molecular dynamics BT - 2016 Editor's Choice of The Journal of Chemical Physics T2 - The Journal of Chemical Physics N2 - We utilize the theory of coherent sets to build Markov state models for non- equilibrium molecular dynamical systems. Unlike for systems in equilibrium, “meta- stable” sets in the non-equilibrium case may move as time evolves. We formalize this concept by relying on the theory of coherent sets, based on this we derive finite-time non-stationary Markov state models, and illustrate the concept and its main differences to equilibrium Markov state modeling on simple, one-dimensional examples. T3 - ZIB-Report - 16-11 KW - coherent set, KW - Markov state model KW - non-equilibrium molecular dynamics KW - metastability Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-57869 SN - 1438-0064 VL - 174103 ET - 145 ER - TY - JOUR A1 - Djurdjevac Conrad, Natasa A1 - Weber, Marcus A1 - Schütte, Christof T1 - Finding dominant structures of nonreversible Markov processes JF - Multiscale Modeling and Simulation Y1 - 2016 U6 - https://doi.org/10.1137/15M1032272 VL - 14 IS - 4 SP - 1319 EP - 1340 ER -