TY - THES A1 - Röblitz, Susanna T1 - Statistical Error Estimation and Grid-free Hierarchical Refinement in Conformation Dynamics N2 - The understanding of geometric structures and dynamical properties of molecular conformations gives insight into molecular long-term behavior. The identification of metastable conformations together with their life times and transition patterns is the intention of conformation dynamics. Conformation dynamics is a multi-scale approach that leads to a reduced description of the dynamical system in terms of a stochastic transition probability matrix. The present thesis deals with the error analysis of computed matrices and the resulting matrix functions. Since conformational membership vectors, as they are computed by the Robust Perron Cluster Analysis (PCCA+), form an invariant subspace of the transition matrix, subspace-based error estimators are of particular interest. The decomposition of the state space into basis functions and the approximation of integrals by Monte-Carlo quadrature give rise to row-wise correlated random matrices, for which stochastic norms are computed. Together with an appropriate statistical model for the distribution of matrix rows, this allows for the calculation of error bounds and error distributions of the invariant subspace and other variables of interest. Equilibration of errors among the basis functions can be achieved by enhanced sampling in regions where the trajectories are mixing slowly. Hierarchical refinement of such basis functions systematically improves the clustering into metastable conformations by reducing the error in the corresponding invariant subspace. These techniques allow for an evaluation of simulation results and pave the way for the analysis of larger molecules. Moreover, the extension of PCCA+ to non-reversible Markov chains, verified by the corresponding perturbation theory, and the modification of the objective function for the case of soft membership vectors represent a further generalization of the clustering method, thus continuing the development from PCCA over PCCA+ to PCCA++. The methods developed in this thesis are useful for but not limited to conformation dynamics. In fact, they are applicable to a broader class of problems which combine domain decomposition with Monte-Carlo quadrature. Possible application areas may include the chemical master equation or quantum dynamical systems. N2 - Das Verständnis von geometrischen Strukturen und dynamischen Eigenschaften molekularer Konformationen ist essentiell für die Vorhersage des Langzeitverhaltens von Molekülen. Die Identifikation metastabiler Konformationen sowie die Bestimmung von Übergangswahrscheinlichkeiten und Haltezeiten sind Bestandteil der Konformationdynamik. Dabei handelt es sich um eine Mehrskalenmethode, die auf eine reduzierte Beschreibung des Systems mittels einer stochastischen Übergangsmatrix führt. In der vorliegenden Dissertation wurde untersucht, wie man die Genauigkeit der Matrizen sowie der daraus berechneten Größen quantifizieren kann. Im Mittelpunkt stehen dabei Fehlerschätzer für den invarianten Unterraum, da die rechten Eigenvektoren als Grundlage der Robusten Perron Cluster Analyse (PCCA+) zur Identifizierung der metastabilen Konformationen dienen. Die Zerlegung des Zustandsraumes in Basisfunktionen sowie die Approximation der Matrixeinträge mittels Monte-Carlo-Quadratur führen zu zeilenweise korrelierten Zufallsmatrizen. Mit Hilfe einer stochastischen Norm sowie einem geeigneten statistischen Modell für die Verteilung der Matrixzeilen können u.a. Fehlerschranken und -verteilungen für den invarianten Unterraum brechnet werden. Eine Equilibrierung des Fehlers zwischen den Basisfunktionen kann durch erweitertes Sampling in solchen Regionen erreicht werden, in denen die Trajektorien nur langsam mischen.Eine hierarchische Zerlegung dieser Basisfunktionen verbessert systematisch die Zerlegung in metastabile Konformationen, indem sie den Fehler im invarianten Unterraum reduziert. Diese Techniken gestatten eine Evaluierung der Simulationsergebnisse und ebnen den Weg zur Behandlung komplexerer Moleküle. Desweiteren wurden Verallgemeinerungen der PCCA+ untersucht. Die Erweiterung der PCCA+ auf nicht-reversible Markov-Ketten sowie die Modifizierung der Zielfunktion für den Fall der weichen Clusterung setzen die Entwicklung von der PCCA über PCCA+ zu PCCA++ fort. Somit können neue Anwendungsfelder für dieses Cluster-Verfahren erschlossen werden. Die Methoden wurden zwar in Rahmen der Konformationsdynamik entwickelt, jedoch lassen sie sich auf eine weite Problemklasse anwenden, in der Gebietszerlegungsverfahren mit Monte-Carlo-Quadratur kombiniert werden. Mögliche Anwendungsgebiete umfassen die chemische Master-Gleichung oder quantenchemische Systeme. KW - Konformationsdynamik KW - Metastabilität KW - gitterfreie Methoden KW - Monte-Carlo-Quadratur KW - Perron-Cluster-Analyse KW - conformation dynamics KW - metastability KW - mesh-free methods KW - Monte Carlo quadrature KW - Perron cluster analysis Y1 - 2009 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:188-fudissthesis000000008079-9 ER - TY - GEN A1 - Lasser, Caroline A1 - Röblitz, Susanna T1 - Computing expectation values for molecular quantum dynamics N2 - We compute expectation values for the solution of the nuclear Schrödinger equation. The proposed particle method consists of three steps: sampling of the initial Wigner function, classical transport of the sampling points, weighted phase space summation for the final computation of the expectation values. The Egorov theorem guarantees that the algorithm is second order accurate with respect to the semiclassical parameter. We present numerical experiments for a two-dimensional torsional potential with three different sets of initial data and for a six-dimensional Henon-Heiles potential. By construction, the computing times scale linearly with the number of initial sampling points and range between three seconds and one hour. T3 - ZIB-Report - 09-30 KW - time dependent Schrödinger equation KW - molecular quantum dynamics KW - expectation values KW - Wigner functions Y1 - 2009 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-11461 SN - 1438-0064 ER - TY - GEN A1 - Boer, H. Marike T. A1 - Stötzel, Claudia A1 - Röblitz, Susanna A1 - Deuflhard, Peter A1 - Veerkamp, Roel F. A1 - Woelders, Henri T1 - A simple mathematical model of the bovine estrous cycle: follicle development and endocrine interactions N2 - Bovine fertility is the subject of extensive research in animal sciences, especially because fertility of dairy cows has declined during the last decades. The regulation of estrus is controlled by the complex interplay of various organs and hormones. Mathematical modeling of the bovine estrous cycle could help in understanding the dynamics of this complex biological system. In this paper we present a mathematical model of the bovine estrous cycle that includes the processes of follicle and corpus luteum development and the key hormones that interact to control these processes. Focus in this paper is on development of the model, but also some simulation results are presented, showing that a set of equations and parameters is obtained that describes the system consistent with empirical knowledge. Even though the majority of the mechanisms that are included are only known qualitatively as stimulatory or inhibitory effects, the model surprisingly well features quantitative observations made in reality. This model of the bovine estrous cycle could be used as a basis for more elaborate models with the ability to study effects of external manipulations and genetic differences. T3 - ZIB-Report - 10-06 KW - cow KW - reproduction KW - hormone patterns KW - ordinary differential equations KW - systems biology Y1 - 2010 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-11682 SN - 1438-0064 ER - TY - GEN A1 - Röblitz, Susanna A1 - Stötzel, Claudia A1 - Deuflhard, Peter A1 - Jones, Hannah M. A1 - Azulay, David-Olivier A1 - van der Graaf, Piet A1 - Martin, Steven T1 - A mathematical model of the human menstrual cycle for the administration of GnRH analogues N2 - This study presents a differential equation model for the feedback mechanisms between Gonadotropin-releasing Hormone (GnRH), Follicle-Stimulating Hormone (FSH), Luteinizing Hormone (LH), development of follicles and corpus luteum, and the production of estradiol (E2), progesterone (P4), inhibin A (IhA), and inhibin B (IhB) during the female menstrual cycle. In contrast to other models, this model does not involve delay differential equations and is based on deterministic modelling of the GnRH pulse pattern, which allows for faster simulation times and efficient parameter identification. These steps were essential to tackle the task of developing a mathematical model for the administration of GnRH analogues. The focus of this paper is on model development for GnRH receptor binding and the integration of a pharmacokinetic/pharmacodynamic model for the GnRH agonist Nafarelin and the GnRH antagonist Cetrorelix into the menstrual cycle model. The final mathematical model describes the hormone profiles (LH, FSH, P4, E2) throughout the menstrual cycle in 12 healthy women. Moreover, it correctly predicts the changes in the cycle following single and multiple dose administration of Nafarelin or Cetrorelix at different stages in the cycle. T3 - ZIB-Report - 11-16 KW - human menstrual cycle KW - systems biology KW - mathematical modelling KW - GnRH analogues KW - differential equations Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-12732 ER - TY - GEN A1 - Stötzel, Claudia A1 - Plöntzke, Julia A1 - Röblitz, Susanna T1 - Advances in modelling of the bovine estrous cycle: Administration of PGF2alpha N2 - Our model of the bovine estrous cycle is a set of ordinary differential equations which generates hormone profiles of successive estrous cycles with several follicular waves per cycle. It describes the growth and decay of the follicles and the corpus luteum, as well as the change of the key substances over time. In this work we describe recent improvements of this model, including the introduction of new components, and elimination of time delays. We validate our model by showing that the simulations agree with observations from synchronization studies and with measured progesterone data after a single dose administration of synthetic prostaglandin F2alpha. T3 - ZIB-Report - 11-17 KW - cow KW - reproduction KW - hormone patterns KW - differential equations KW - systems biology Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-12740 ER - TY - GEN A1 - Berg, Mascha A1 - Ritter, Julia A1 - Röblitz, Susanna T1 - A mathematical model for the Zinc response of young pigs N2 - In a number of studies the feeding of pharmacological levels of Zinc to weanling pigs was examined. The publications agree in the fact that the achieved performance increase of the pigs depends on the amount of Zinc as well as on the type of mineral complex used as feed additive. The results are based on measurements of the animals' body weight and the Zinc concentration in the blood within the first two weeks of the experiments. From the measurement values, conclusions are drawn about adequate Zinc doses, but a deeper understanding of the physiological processes that influence the Zinc distribution in the organism is still missing. In this report, we want to close this gap by using a systems biology approach. We present a small mathematical model that describes the organic feedback mechanisms with the help of differential equations. The simulation results not only match with measurement values, but also give insight into non-measurable processes, which will hopefully contribute to animal experiments being organized in a more selective and reduced way. T3 - ZIB-Report - 11-05 KW - systems biology KW - ordinary differential equations KW - feeding study KW - weanling pigs Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-12191 ER - TY - JOUR A1 - Boer, H. Marike T. A1 - Stötzel, Claudia A1 - Röblitz, Susanna A1 - Deuflhard, Peter A1 - Veerkamp, Roel F. A1 - Woelders, Henri T1 - A simple mathematical model of the bovine estrous cycle JF - Journal of Theoretical Biology Y1 - 2011 U6 - https://doi.org/10.1016/j.jtbi.2011.02.021 VL - 278 IS - 1 SP - 20 EP - 31 ER - TY - JOUR A1 - Boer, H. Marike T. A1 - Röblitz, Susanna A1 - Stötzel, Claudia A1 - Veerkamp, Roel F. A1 - Kemp, Bas A1 - Woelders, Henri T1 - Mechanisms regulating follicle wave patterns in the bovine estrous cycle investigated with a mathematical model JF - Journal of Dairy Science Y1 - 2011 U6 - https://doi.org/10.3168/jds.2011-4400 VL - 94 IS - 12 SP - 5987 EP - 6000 ER - TY - CHAP A1 - Boer, H. Marike T. A1 - Stötzel, Claudia A1 - Röblitz, Susanna A1 - Woelders, Henri T1 - A differential equation model to investigate the dynamics of the bovine estrous cycle T2 - Advances in Systems Biology Y1 - 2012 U6 - https://doi.org/10.1007/978-1-4419-7210-1_35 VL - 736 SP - 597 EP - 606 ER - TY - GEN A1 - Plöntzke, Julia A1 - Berg, Mascha A1 - Olany, Amiya A1 - Leonhard-Marek, Sabine A1 - Müller, Kerstin-Elisabeth A1 - Röblitz, Susanna T1 - Modeling potassium balance in dairy cows N2 - Potassium is fundamental for cell functioning including signal transduction, acid-base- and water-metabolism. Since diet of dairy cows is generally rich in potassium, hypokalemia was not in the focus of research for long time. Furthermore, hypokalemia was not frequently diagnosed because blood potassium content is difficult to measure. In recent years, measurement methods have been improved. Nowadays hypokalemia is increasingly diagnosed in cows with disorders such as abomasal displacement, ketosis or down cow syndrome, calling for intensified research on this topic. In this report we describe the development of a basic mechanistic, dynamic model of potassium balance based on ordinary differential and algebraic equations. Parameter values are obtained from data of a clinical trial in which potassium balance and the influence of therapeutic intervention in glucose and electrolyte metabolism on potassium balance in non-lactating dairy cows were studied. The model is formulated at a high abstraction level and includes information and hypotheses from literature. This work represents a first step towards the understanding and design of effective prophylactic feed additives and treatment strategies. T3 - ZIB-Report - 13-09 KW - systems biology KW - differential equations KW - metabolism Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-17749 SN - 1438-0064 ER - TY - JOUR A1 - Haack, Fiete A1 - Röblitz, Susanna A1 - Scharkoi, Olga A1 - Schmidt, Burkhard A1 - Weber, Marcus T1 - Adaptive Spectral Clustering for Conformation Analysis JF - AIP Conference Proceedings Y1 - 2010 UR - http://link.aip.org/link/?APC/1281/1585/1 U6 - https://doi.org/10.1063/1.3498116 VL - 1281 IS - 1 SP - 1585 EP - 1588 PB - AIP ER - TY - JOUR A1 - Lasser, Caroline A1 - Röblitz, Susanna T1 - Computing expectation values for molecular quantum dynamics JF - SIAM Journal on Scientific Computing Y1 - 2010 U6 - https://doi.org//10.1137/090770461 VL - 32 IS - 3 SP - 1465 EP - 1483 ER - TY - JOUR A1 - Stötzel, Claudia A1 - Plöntzke, Julia A1 - Heuwieser, Wolfgang A1 - Röblitz, Susanna T1 - Advances in modelling of the bovine estrous cycle: Synchronization with PGF2alpha JF - Theriogenology Y1 - 2012 U6 - https://doi.org/10.1016/j.theriogenology.2012.04.017 VL - 78 SP - 1415 EP - 1428 ER - TY - JOUR A1 - Stötzel, Claudia A1 - Ehrig, Rainald A1 - Boer, H. Marike T. A1 - Plöntzke, Julia A1 - Röblitz, Susanna T1 - Exploration of different wave patterns in a model of the bovine estrous cycle by Fourier analysis JF - BIOMAT - Proceedings of the 14th International Symposium on Mathematical and Computational Biology, Bedlewo, Poland N2 - Cows typically have different numbers of follicular waves during their hormonal cycle. Understanding the underlying regulations leads to insights into the reasons for declined fertility, a phenomenon that has been observed during the last decades. We present a systematic approach based on Fourier analysis to examine how parameter changes in a model of the bovine estrous cycle lead to different wave patterns. Even without any biological considerations, this allows to detect the responsible model parameters that control the type of periodicity of the solution, thus supporting experimental planning of animal scientists. Y1 - 2014 ER - TY - GEN A1 - Stötzel, Claudia A1 - Apri, Mochamad A1 - Röblitz, Susanna T1 - A reduced ODE model of the bovine estrous cycle N2 - This work deals with the reduction of a previously developed ODE model for the bovine estrous cycle. After applying a method for exploring the parameter space by Apri et al. (2012), we perform structure-based reduction steps and several system specific adaptations. Overall, the original model consisting of 15 ODEs and 60 parameters is reduced to a model of 10 odes and 38 parameters. The reduced model qualitatively reproduces the state trajectories of the original model. T3 - ZIB-Report - 14-33 Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-51594 SN - 1438-0064 ER - TY - JOUR A1 - Berg, Mascha A1 - Plöntzke, Julia A1 - Leonhard-Marek, Sabine A1 - Müller, Kerstin-Elisabeth A1 - Röblitz, Susanna T1 - A dynamic model to simulate potassium balance in dairy cows JF - Journal of Dairy Science N2 - High-performing dairy cows require a particular composition of nutritional ingredients, adapted to their individual requirements and depending on their production status. The optimal dimensioning of minerals in the diet, one being potassium, is indispensable for the prevention of imbalances. Potassium balance in cows is the result of potassium intake, distribution in the organism, and excretion, and it is closely related to glucose and electrolyte metabolism. In this paper, we present a dynamical model for potassium balance in lactating and nonlactating dairy cows based on ordinary differential equations. Parameter values were obtained from clinical trial data and from the literature. To verify the consistency of the model, we present simulation outcomes for 3 different scenarios: potassium balance in (1) nonlactating cows with varying feed intake, (2) nonlactating cows with varying potassium fraction in the diet, and (3) lactating cows with varying milk production levels. The results give insights into the short- and long-term potassium metabolism, providing an important step toward the understanding of the potassium network, the design of prophylactic feed additives, and possible treatment strategies. Y1 - 2017 U6 - https://doi.org/10.3168/jds.2016-12443 VL - 100 IS - 12 SP - 9799 EP - 9814 ER - TY - JOUR A1 - Klebanov, Ilja A1 - Sikorski, Alexander A1 - Schütte, Christof A1 - Röblitz, Susanna T1 - Objective priors in the empirical Bayes framework JF - Scandinavian Journal of Statistics N2 - When dealing with Bayesian inference the choice of the prior often remains a debatable question. Empirical Bayes methods offer a data-driven solution to this problem by estimating the prior itself from an ensemble of data. In the nonparametric case, the maximum likelihood estimate is known to overfit the data, an issue that is commonly tackled by regularization. However, the majority of regularizations are ad hoc choices which lack invariance under reparametrization of the model and result in inconsistent estimates for equivalent models. We introduce a nonparametric, transformation-invariant estimator for the prior distribution. Being defined in terms of the missing information similar to the reference prior, it can be seen as an extension of the latter to the data-driven setting. This implies a natural interpretation as a trade-off between choosing the least informative prior and incorporating the information provided by the data, a symbiosis between the objective and empirical Bayes methodologies. Y1 - 2021 U6 - https://doi.org/10.1111/sjos.12485 VL - 48 IS - 4 SP - 1212 EP - 1233 PB - Wiley Online Library ER - TY - JOUR A1 - Leeners, Brigitte A1 - Krüger, Tillmann H. C. A1 - Geraedts, Kirsten A1 - Tronci, Enrico A1 - Mancini, Toni A1 - Ille, Fabian A1 - Egli, Marcel A1 - Röblitz, Susanna A1 - Saleh, Lanja A1 - Spanhaus, K. A1 - Schippert, Cordula A1 - Zhang, Y. A1 - Hengartner, M. P. T1 - Lack of association between female hormone levels and visuospatial working memory, divided attention and cognitive bias across two consecutive menstrual cycles JF - Front. Behav. Neurosci. Y1 - 2017 U6 - https://doi.org/10.3389/fnbeh.2017.00120 ER - TY - GEN A1 - Klebanov, Ilja A1 - Sikorski, Alexander A1 - Schütte, Christof A1 - Röblitz, Susanna T1 - Prior estimation and Bayesian inference from large cohort data sets N2 - One of the main goals of mathematical modelling in systems biology related to medical applications is to obtain patient-specific parameterisations and model predictions. In clinical practice, however, the number of available measurements for single patients is usually limited due to time and cost restrictions. This hampers the process of making patient-specific predictions about the outcome of a treatment. On the other hand, data are often available for many patients, in particular if extensive clinical studies have been performed. Using these population data, we propose an iterative algorithm for contructing an informative prior distribution, which then serves as the basis for computing patient-specific posteriors and obtaining individual predictions. We demonsrate the performance of our method by applying it to a low-dimensional parameter estimation problem in a toy model as well as to a high-dimensional ODE model of the human menstrual cycle, which represents a typical example from systems biology modelling. T3 - ZIB-Report - 16-09 Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-57475 SN - 1438-0064 ER - TY - GEN A1 - Gupta, Pooja A1 - Röblitz, Susanna A1 - Krause, Carola A1 - Knaus, Petra A1 - Schütte, Christof T1 - Mathematical modeling of the Smad and Non-Smad BMP signaling pathways in context of cell density T2 - Computational Models in biology and medicine, 2013, Dresden, Germany Y1 - 2013 ER - TY - CHAP A1 - Gupta, Pooja A1 - Krause, Carola A1 - Rikeit, Paul A1 - Röblitz, Susanna A1 - Knaus, Petra A1 - Schütte, Christof T1 - Modeling of the BMP mediated co-regulation of the Smad and Non-Smad pathways in the context of cell density T2 - 10th International BMP conference, 2014, Berlin, Germany Y1 - 2014 ER - TY - CHAP A1 - Schäfer, Stefan A1 - Plöntzke, Julia A1 - Röblitz, Susanna T1 - Mathematical Modelling of Follicular Maturation in Cows and Women T2 - 49. Jahrestagung der Physiologie und Pathologie der Fortpflanzung und gleichzeitig 41. Veterinär-Humanmedizinische Gemeinschaftstagung 2016, Leipzig. Y1 - 2016 ER - TY - CHAP A1 - Igde, Sinaida A1 - Wölk, Hendrik A1 - Röblitz, Susanna A1 - Reidelbach, Marco A1 - Weber, Marcus A1 - Hartmann, Laura T1 - Identifying Multivalent Binding Kinetics of Precision Glycomacromolecules: A Kinetic Study Using kinITC T2 - Münster Symposium on Cooperative Effects 2015 - SFB 858, at Westfälische Wilhelms-Universität Münster, 2015 N2 - Multivalent sugar/protein interactions are well-known to proceed through different binding modes 1-5 which in turn can be described by their binding kinetics 3-5. This study provides additional insight into the association and dissociation reaction rates of complex multivalent sugar/protein interactions. Binding kinetics of recently introduced multivalent precision glycomacromolecules 6-8 to Concanavalin A (Con A) were studied by " kinetic Isothermal Titration Calorimetry " (kinITC) 9-11. The effect of multivalency is evaluated by comparing rate constants of glycomacromolecules obtaining the same and different valency of mannose ligands and by variation of the overall backbone properties, such as hydrophilic/ hydrophoboc. In addition, binding kinetics were studied using different conformations of Con A (homodimer vs.-tetramer) and thus a different protein valency. Our results show that precision glycomacromolecule/Con A binding proceeds non-cooperatively. Further, association and dissociation rates are mainly described by intermolecular complex formation. Together with the so-called functional valency, we can discriminate between " bound " and " unbound " states for macroscopic on-and off-rates, even for such complex glycooligomer/protein systems. By comparing e.g. a mono-to a divalent glycomacromolecule for their binding to dimeric Con A, we see a lower dissociation rate for the latter. As both bind monovalently to Con A, this is a strong indication for a statistical rebinding event. Further, there is a strong dependence of multivalent binding kinetics on the ligand density of glycomacromolecules as well as the Con A conformation and thus the overall on-and off-rates. Y1 - 2015 ER - TY - GEN A1 - Klebanov, Ilja A1 - Sikorski, Alexander A1 - Schütte, Christof A1 - Röblitz, Susanna T1 - Empirical Bayes Methods for Prior Estimation in Systems Medicine N2 - One of the main goals of mathematical modelling in systems medicine related to medical applications is to obtain patient-specific parameterizations and model predictions. In clinical practice, however, the number of available measurements for single patients is usually limited due to time and cost restrictions. This hampers the process of making patient-specific predictions about the outcome of a treatment. On the other hand, data are often available for many patients, in particular if extensive clinical studies have been performed. Therefore, before applying Bayes’ rule separately to the data of each patient (which is typically performed using a non-informative prior), it is meaningful to use empirical Bayes methods in order to construct an informative prior from all available data. We compare the performance of four priors - a non-informative prior and priors chosen by nonparametric maximum likelihood estimation (NPMLE), by maximum penalized lilelihood estimation (MPLE) and by doubly-smoothed maximum likelihood estimation (DS-MLE) - by applying them to a low-dimensional parameter estimation problem in a toy model as well as to a high-dimensional ODE model of the human menstrual cycle, which represents a typical example from systems biology modelling. T3 - ZIB-Report - 16-57 KW - Parameter estimation KW - Bayesian inference KW - Bayesian hierarchical modelling KW - NPMLE KW - MPLE KW - DS-MLE KW - EM algorithm KW - Jeffreys prior KW - reference prior KW - hyperparameter KW - hyperprior KW - principle of maximum entropy Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-61307 SN - 1438-0064 ER - TY - CHAP A1 - Lange, Alexander A1 - Plöntzke, Julia A1 - Schäfer, Stefan A1 - Röblitz, Susanna T1 - Follicular maturation in cows: mathematical models and data T2 - 10. European Conference on Mathematical and Theoretical Biology Y1 - 2016 ER - TY - GEN A1 - Omari, Mohamed A1 - Plöntzke, Julia A1 - Berg, Mascha A1 - Röblitz, Susanna T1 - Modelling of glucose-insulin metabolism and its effect on the estrous cycle in bovine T2 - Annual Meeting of the Population Approach Group in Europe Y1 - 2016 ER - TY - JOUR A1 - Hengartner, Michael A1 - Krüger, Tillmann H. C. A1 - Geraedts, Kirsten A1 - Tronci, Enrico A1 - Mancini, Toni A1 - Ille, Fabian A1 - Egli, Marcel A1 - Röblitz, Susanna A1 - Ehrig, Rainald A1 - Saleh, Lanja A1 - Spanaus, Katharina A1 - Schippert, Cordula A1 - Zhang, Yuanyuan A1 - Leeners, Brigitte T1 - Negative affect is unrelated to fluctuations in hormone levels across the menstrual cycle: evidence from a multisite observational study across two successive cycles JF - Journal of Psychosomatic Research Y1 - 2017 U6 - https://doi.org/10.1016/j.jpsychores.2017.05.018 VL - 99 SP - 21 EP - 27 ER - TY - GEN A1 - Lange, Alexander A1 - Schwieger, Robert A1 - Plöntzke, Julia A1 - Schäfer, Stefan A1 - Röblitz, Susanna T1 - Follicular competition: the selection of dominant follicles as a synergistic effect N2 - The estrous cycle of mono-ovulatory species such as cows or humans, is known to show two or more waves of follicular growth and decline between two successive ovulations. Within each wave, there is one dominant follicle escorted by subordinate follicles of varying number. Under the surge of the luteinizing hormone a growing dominant follicle ovulates. Rarely the number of ovulating follicles exceeds one. In the biological literature, the change of hormonal concentrations and individually varying numbers of follicular receptors are made responsible for the selection of exactly one dominant follicle, yet a clear cause has not been identified. In this paper, we suggest a synergistic explanation based on competition, formulated by a parsimoniously defined system of ordinary differential equations (ODEs) that quantifies the time evolution of multiple follicles and their competitive interaction during one wave. Not discriminating between follicles, growth and decline are given by fixed rates. Competition is introduced via a growth-suppressing term, equally supported by all follicles. We prove that the number of dominant follicles is determined exclusively by the ratio of follicular growth and competition. This number turns out to be independent of the number of subordinate follicles. The asymptotic behavior of the corresponding dynamical system is investigated rigorously, where we demonstrate that the omega-limit set only contains fixed points. When also including follicular decline, our ODEs perfectly resemble ultrasound data of bovine follicles. Implications for the involved but not explicitly modeled hormones are discussed. T3 - ZIB-Report - 17-21 Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-63863 SN - 1438-0064 ER - TY - GEN A1 - Bockmayr, Alexander A1 - Siebert, Heike A1 - Röblitz, Susanna A1 - Schütte, Christof A1 - Deuflhard, Peter ED - Deuflhard, Peter ED - Grötschel, Martin ED - Hömberg, Dietmar ED - Kramer, Jürg ED - Mehrmann, Volker ED - Polthier, Konrad ED - Schmidt, Frank ED - Schütte, Christof ED - Skutela, Martin ED - Sprekels, Jürgen T1 - Advanced mathematical modeling in systems biology T2 - MATHEON-Mathematics for Key Technologies Y1 - 2014 VL - 1 SP - 29 EP - 44 PB - European Mathematical Society ER - TY - JOUR A1 - Lang, Annemarie A1 - Volkamer, Andrea A1 - Behm, Laura A1 - Röblitz, Susanna A1 - Ehrig, Rainald A1 - Schneider, Marlon A1 - Geris, Lisbet A1 - Wichard, Joerg A1 - Buttgereit, Frank T1 - In silico Methods - Computational Alternatives to Animal Testing JF - ALTEX N2 - A seminar and interactive workshop on “In silico Methods – Computational Alternatives to Animal Testing” was held in Berlin, Germany, organized by Annemarie Lang, Frank Butt- gereit and Andrea Volkamer at the Charité-Universitätsmedizin Berlin, on August 17-18, 2017. During the half-day seminar, the variety and applications of in silico methods as alternatives to animal testing were presented with room for scientific discus- sions with experts from academia, industry and the German fed- eral ministry (Fig. 1). Talks on computational systems biology were followed by detailed information on predictive toxicology in order to display the diversity of in silico methods and the potential to embrace them in current approaches (Hartung and Hoffmann, 2009; Luechtefeld and Hartung, 2017). The follow- ing interactive one-day Design Thinking Workshop was aimed at experts, interested researchers and PhD-students interested in the use of in silico as alternative methods to promote the 3Rs (Fig. 2). Forty participants took part in the seminar while the workshop was restricted to sixteen participants. Y1 - 2018 U6 - https://doi.org/10.14573/altex.1712031 VL - 35 IS - 1 SP - 126 EP - 128 ER - TY - JOUR A1 - Igde, Sinaida A1 - Röblitz, Susanna A1 - Müller, Anne A1 - Kolbe, Katharina A1 - Boden, Sophia A1 - Fessele, Claudia A1 - Lindhorst, Thisbe A1 - Weber, Marcus A1 - Hartmann, Laura T1 - Linear Precision Glycomacromolecules with Varying Interligand Spacing and Linker Functionalities Binding to Concanavalin A and the Bacterial Lectin FimH JF - Marcomolecular Bioscience Y1 - 2017 U6 - https://doi.org/10.1002/mabi.201700198 VL - 17 IS - 12 SP - 1700198 ER - TY - GEN A1 - Schwieger, Robert A1 - Siebert, Heike A1 - Röblitz, Susanna T1 - Correspondence of Trap Spaces in Different Models of Bioregulatory Networks N2 - Mathematical models for bioregulatory networks can be based on different formalisms, depending on the quality of available data and the research question to be answered. Discrete boolean models can be constructed based on qualitative data, which are frequently available. On the other hand, continuous models in terms of ordinary differential equations (ODEs) can incorporate time-series data and give more detailed insight into the dynamics of the underlying system. A few years ago, a method based on multivariate polynomial interpolation and Hill functions has been developed for an automatic conversion of boolean models to systems of ordinary differential equations. This method is frequently used by modellers in systems biology today, but there are only a few results available about the conservation of mathematical structures and properties across the formalisms. Here, we consider subsets of the phase space where some components stay fixed, called trap spaces, and demonstrate how boolean trap spaces can be linked to invariant sets in the continuous state space. This knowledge is of practical relevance since finding trap spaces in the boolean setting, which is relatively easy, allows for the construction of reduced ODE models. T3 - ZIB-Report - 18-10 KW - dynamical systems KW - trap spaces KW - attractors KW - invariant sets KW - ordinary differential equations KW - Boolean networks Y1 - 2018 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-67121 SN - 1438-0064 ER - TY - JOUR A1 - Röblitz, Susanna A1 - Stötzel, Claudia A1 - Deuflhard, Peter A1 - Jones, Hannah M. A1 - Azulay, David-Olivier A1 - van der Graaf, Piet A1 - Martin, Steven T1 - A mathematical model of the human menstrual cycle for the administration of GnRH analogues JF - Journal of Theoretical Biology Y1 - 2013 U6 - https://doi.org/10.1016/j.jtbi.2012.11.020 VL - 321 SP - 8 EP - 27 ER - TY - GEN A1 - Bondouy, Manon A1 - Röblitz, Susanna T1 - Mathematical modeling of follicular development in bovine estrous cycles N2 - Normal follicular growth is a prerequisite for successful fertilization of dairy cows. However, submission rates have been decreasing during the last decades because animals are not ovulating at the planned time of mating or insemination. In this study we use a mathematical model to investigate mechanisms that lead to anestrus. This model is derived by coupling two previously published models: a small model for the development of multiple follicles (Smith et al., 2004), and a large estrous cycle model (St\"otzel et al., 2012). We first investigate the influence of synchronization protocols on the time-shift of ovulation. In a second scenario we simulate an extended period of anestrus as it typically occurs after calving. T3 - ZIB-Report - 12-26 KW - ordinary differential equations KW - systems biology KW - reproduction KW - dairy cows Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-15585 SN - 1438-0064 ER - TY - GEN A1 - Meisel, Kristin A1 - Röblitz, Susanna T1 - Mathematische Modellierung der Knochenremodellierung N2 - Knochenremodellierung beinhaltet den Auf- und Abbau der Knochenmasse durch die verschiedenen Knochenzellen und findet fast überall, auch am Skelett Erwachsener, statt. Für die Erneuerung der Knochensubstanz sind die Osteoblasten zuständig. Sie ersetzen exakt die Menge der Knochenmasse, welche zuvor durch Osteoklasten abgebaut wurde. Störungen dieses, durch viele Faktoren beeinflussten Prozesses führen zu pathologischen Veränderungen, beispielsweise zu Osteoporose oder Arthritis. Auf der Grundlage analysierter publizierter Modelle wurde ein Modell entwickelt, welches den Einfluss einiger dieser Faktoren realistisch abbildet. Auf diese Weise kann die Wirkung von basalem PTH, des RANKL-OPG-Systems ( RANKL: Receptor activator of nuclear factor kappa-B ligand, OPG: Osteoprotegerin) und des Estradiols auf den Knochenstoffwechsel am vorgestellten System untersucht werden. Außerdem wurde durch Estradiolmangel hervorgerufene Osteoporose und der kurative Effekt von synthetischen Medikamenten wie Estradiol oder intermittierend verabreichtem PTH modelliert. Mit der Parameterschätzung anhand des Gauß-Newton-Verfahrens wird des Weiteren eine Methode vorgestellt, die es ermöglicht, mathematische Modelle bestmöglich durch Variation der Parameterwerte an experimentelle Daten anzupassen. Das am Zuse-Institut Berlin entwickelte Softwarepaket POEM wendet diesen Algorithmus an und wird ebenfalls erläutert. T3 - ZIB-Report - 13-41 Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-42010 SN - 1438-0064 ER - TY - GEN A1 - Stötzel, Claudia A1 - Röblitz, Susanna A1 - Siebert, Heike T1 - Complementing ODE-based System Analysis Using Boolean Networks Derived from an Euler-like Transformation N2 - In this paper, we present a systematic transition scheme for a large class of ordinary differential equations (ODEs) into Boolean networks. Our transition scheme can be applied to any system of ODEs whose right hand sides can be written as sums and products of monotone functions. It performs an Euler-like step which uses the signs of the right hand sides to obtain the Boolean update functions for every variable of the corresponding discrete model. The discrete model can, on one hand, be considered as another representation of the biological system or, alternatively, it can be used to further the analysis of the original ODE model. Since the generic transformation method does not guarantee any property conservation, a subsequent validation step is required. Depending on the purpose of the model this step can be based on experimental data or ODE simulations and characteristics. Analysis of the resulting Boolean model, both on its own and in comparison with the ODE model, then allows to investigate system properties not accessible in a purely continuous setting. The method is exemplarily applied to a previously published model of the bovine estrous cycle, which leads to new insights regarding the regulation among the components, and also indicates strongly that the system is tailored to generate stable oscillations. T3 - ZIB-Report - 15-46 Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-56164 SN - 1438-0064 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 - GEN A1 - Stötzel, Claudia A1 - Ehrig, Rainald A1 - Boer, H. Marike T. A1 - Plöntzke, Julia A1 - Röblitz, Susanna T1 - Exploration of different wave patterns in a model of the bovine estrous cycle by Fourier analysis N2 - Cows typically have different numbers of follicular waves during their hormonal cycle. Understanding the underlying regulations leads to insights into the reasons for declined fertility, a phenomenon that has been observed during the last decades. We present a systematic approach based on Fourier analysis to examine how parameter changes in a model of the bovine estrous cycle lead to different wave patterns. Even without any biological considerations, this allows to detect the responsible model parameters that control the type of periodicity of the solution, thus supporting experimental planning of animal scientists. T3 - ZIB-Report - 16-02 Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-57039 SN - 1438-0064 ER - TY - GEN A1 - Plöntzke, Julia A1 - Berg, Mascha A1 - Stötzel, Claudia A1 - Röblitz, Susanna T1 - A systems biology approach to bovine fertility and metabolism: Development of a glucose insulin model. N2 - To counteract the antagonistic relationship between milk yield and fertility in dairy cow, a deeper understanding of the underlying biological mechanisms is required. For this purpose, we study physiological networks related to reproduction and metabolism in dairy cows. We interactively develop dynamic, mechanistic models by fitting the models to experimental data and mechanistic knowledge. We have already developed models for potassium balance and hormonal regulation of fertility in the dairy cow, which will briefly be reviewed here. The main focus of this article is a glucose-insulin model currently developed by us. This model links the bovine hormonal cycle and the potassium balance to glucose and thus to energy metabolism. The models can be applied in scientific research, education, experimental planning, drug development and production on farms. T3 - ZIB-Report - 15-51 Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-56453 SN - 1438-0064 ER - TY - JOUR A1 - Stötzel, Claudia A1 - Röblitz, Susanna A1 - Siebert, Heike T1 - Complementing ODE-based System Analysis Using Boolean Networks Derived from an Euler-like Transformation JF - PLOS ONE Y1 - 2015 U6 - https://doi.org/10.1371/journal.pone.0140954 ER - TY - CHAP A1 - Plöntzke, Julia A1 - Berg, Mascha A1 - Stötzel, Claudia A1 - Röblitz, Susanna T1 - A systems biology approach to bovine fertility and metabolism: Introduction of a glucose insulin model T2 - 15th International Symposium on Mathematical and Computational Biology, Rorkee, India Y1 - 2015 ER - TY - BOOK A1 - Deuflhard, Peter A1 - Röblitz, Susanna T1 - A guide to numerical modelling in systems biology Y1 - 2015 SN - 978-3-319-20058-3 U6 - https://doi.org/10.1007/978-3-319-20059-0 VL - 12 PB - Springer ER - TY - JOUR A1 - Deuflhard, Peter A1 - Röblitz, Susanna A1 - Stötzel, Claudia T1 - Mathematical Systems Biology at Zuse Institute Berlin (ZIB) JF - BioTOPics - Journal of Biotechnology in Berlin-Brandenburg Y1 - 2010 VL - 39 ER - TY - JOUR A1 - Röblitz, Susanna A1 - Weber, Marcus T1 - Fuzzy spectral clustering by PCCA+: application to Markov state models and data classification JF - Advances in Data Analysis and Classification Y1 - 2013 U6 - https://doi.org/10.1007/s11634-013-0134-6 VL - 7 IS - 2 SP - 147 EP - 179 ER - TY - JOUR A1 - Röblitz, Susanna A1 - Weber, Marcus T1 - Fuzzy Spectral Clustering by PCCA+ JF - Classification and Clustering: Models, Software and Applications Y1 - 2009 UR - http://www.wias-berlin.de/publications/wias-publ/run.jsp?template=abstract&type=Report&year=2009&number=26 IS - WIAS Report No. 26 SP - 73 EP - 79 ER - TY - GEN A1 - Moualeu-Ngangue, Dany Pascal A1 - Röblitz, Susanna A1 - Ehrig, Rainald A1 - Deuflhard, Peter T1 - Parameter Identification for a Tuberculosis Model in Cameroon N2 - A deterministic model of tuberculosis in sub-Saharan Africa in general and Cameroon in particular including lack of access to the treatment and weak diagnose capacity is designed and analyzed with respect to its transmission dynamics. The model includes both frequency- and density-dependent transmissions. It is shown that the model is mathematically well-posed and epidemiologically reasonable. Solutions are non-negative and bounded whenever the initial values are non-negative. A sensitivity analysis of model parameters is performed and most sensitive parameters of the model are identified using a state-of-the-art Gauss-Newton Method. In particular, parameters representing the proportion of individuals having access to medical facilities have a large impact on the dynamics of the disease. It has been shown that an increase of these parameter values over the time can significantly reduce the disease burden in the population within the next 15 years. T3 - ZIB-Report - 13-72 KW - Tuberculosis, parameter identification, sensitivity analysis, Gauss-Newton method, Nonlinear dynamical systems. Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-43087 SN - 1438-0064 ER - TY - JOUR A1 - Moualeu-Ngangue, Dany Pascal A1 - Röblitz, Susanna A1 - Ehrig, Rainald A1 - Deuflhard, Peter T1 - Parameter Identification in a Tuberculosis Model for Cameroon JF - PLOS ONE N2 - A deterministic model of tuberculosis in Cameroon is designed and analyzed with respect to its transmission dynamics. The model includes lack of access to treatment and weak diagnosis capacity as well as both frequency- and density-dependent transmissions. It is shown that the model is mathematically well-posed and epidemiologically reasonable. Solutions are non-negative and bounded whenever the initial values are non-negative. A sensitivity analysis of model parameters is performed and the most sensitive ones are identified by means of a state-of-the-art Gauss-Newton method. In particular, parameters representing the proportion of individuals having access to medical facilities are seen to have a large impact on the dynamics of the disease. The model predicts that a gradual increase of these parameters could significantly reduce the disease burden on the population within the next 15 years. Y1 - 2015 U6 - https://doi.org/10.1371/journal.pone.0120607 ER - TY - JOUR A1 - Weber, Marie-Christin A1 - Fischer, Lisa A1 - Damerau, Alexandra A1 - Ponomarev, Igor A1 - Pfeiffenberger, Moritz A1 - Gaber, Timo A1 - Götschel, Sebastian A1 - Lang, Jens A1 - Röblitz, Susanna A1 - Buttgereit, Frank A1 - Ehrig, Rainald A1 - Lang, Annemarie T1 - In vitro and in silico modeling of cellular and matrix-related changes during the early phase of osteoarthritis JF - BioRxiv N2 - Understanding the pathophysiological processes of osteoarthritis (OA) require adequate model systems. Although different in vitro or in vivo models have been described, further comprehensive approaches are needed to study specific parts of the disease. This study aimed to combine in vitro and in silico modeling to describe cellular and matrix-related changes during the early phase of OA. We developed an in vitro OA model based on scaffold-free cartilage-like constructs (SFCCs), which was mathematically modeled using a partial differential equation (PDE) system to resemble the processes during the onset of OA. SFCCs were produced from mesenchymal stromal cells and analyzed weekly by histology and qPCR to characterize the cellular and matrix-related composition. To simulate the early phase of OA, SFCCs were treated with interleukin-1β (IL-1β), tumor necrosis factor α (TNFα) and examined after 3 weeks or cultivated another 3 weeks without inflammatory cytokines to validate the regeneration potential. Mathematical modeling was performed in parallel to the in vitro experiments. SFCCs expressed cartilage-specific markers, and after stimulation an increased expression of inflammatory markers, matrix degrading enzymes, a loss of collagen II (Col-2) and a reduced cell density was observed which could be partially reversed by retraction of stimulation. Based on the PDEs, the distribution processes within the SFCCs, including those of IL-1β, Col-2 degradation and cell number reduction was simulated. By combining in vitro and in silico methods, we aimed to develop a valid, efficient alternative approach to examine and predict disease progression and new therapeutic strategies. Y1 - 2019 U6 - https://doi.org/10.1101/725317 ER - TY - JOUR A1 - Lang, Annemarie A1 - Fischer, Lisa A1 - Weber, Marie-Christin A1 - Gaber, Timo A1 - Ehrig, Rainald A1 - Röblitz, Susanna A1 - Buttgereit, Frank T1 - Combining in vitro simulation and in silico modelling towards a sophisticated human osteoarthritis model JF - Osteoarthritis and Cartilage N2 - Our project aimed at building an in silico model based on our recently developed in vitro osteoarthritis (OA) model seeking for refinement of the model to enhance validity and translatability towards the more sophisticated simulation of OA. In detail, the previously 3D in vitro model is based on 3D chondrogenic constructs generated solely from human bone marrow derived mesenchymal stromal cells (hMSCs). Besides studying the normal state of the model over 3 weeks, the in vitro model was treated with interleukin-1β (IL-1β) and tumor necrosis factor alpha (TNFα) to mimic an OA-like environment. Y1 - 2019 U6 - https://doi.org/10.1016/j.joca.2019.02.277 VL - 27 SP - S183 ER - TY - GEN A1 - Omari, Mohamed A1 - Plöntzke, Julia A1 - Röblitz, Susanna T1 - A pharmacokinetic-pharmacodynamic model for single dose administration of Dexamethasone in dairy cows N2 - We present a mechanistic pharmacokinetic-pharmacodynamic model to simulate the effect of dexamethasone on the glucose metabolism in dairy cows. The coupling of the pharmacokinetic model to the pharmacodynamic model is based on mechanisms underlying homeostasis regulation by dexamethasone. In particular, the coupling takes into account the predominant role of dexamethasone in stimulating glucagon secretion, glycogenolysis and lipolysis and in impairing the sensitivity of cells to insulin. Simulating the effect of a single dose of dexamethasone on the physiological behaviour of the system shows that the adopted mechanisms are able to induce a temporary hyperglycemia and hyperinsulinemia, which captures the observed data in non-lactating cows. In lactating cows, the model simulations show that a single dose of dexamethasone reduces the lipolytic effect, owing to the reduction of glucose uptake by the mammary gland. T3 - ZIB-Report - 19-53 KW - systems biology; modelling; differential equations; metabolism; nutrition; bovine; fertility; reproduction; hormones; follicles; Dexamethasone; PK-PD Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-75043 SN - 1438-0064 ER - TY - JOUR A1 - Weber, Marie-Christin A1 - Fischer, Lisa A1 - Damerau, Alexandra A1 - Ponomarev, Igor A1 - Pfeiffenberger, Moritz A1 - Gaber, Timo A1 - Götschel, Sebastian A1 - Lang, Jens A1 - Röblitz, Susanna A1 - Buttgereit, Frank A1 - Ehrig, Rainald A1 - Lang, Annemarie T1 - Macroscale mesenchymal condensation to study cytokine-driven cellular and matrix-related changes during cartilage degradation JF - Biofabrication N2 - Understanding the pathophysiological processes of cartilage degradation requires adequate model systems to develop therapeutic strategies towards osteoarthritis (OA). Although different in vitro or in vivo models have been described, further comprehensive approaches are needed to study specific disease aspects. This study aimed to combine in vitro and in silico modeling based on a tissue-engineering approach using mesenchymal condensation to mimic cytokine-induced cellular and matrix-related changes during cartilage degradation. Thus, scaffold-free cartilage-like constructs (SFCCs) were produced based on self-organization of mesenchymal stromal cells (mesenchymal condensation) and i) characterized regarding their cellular and matrix composition or secondly ii) treated with interleukin-1β (IL-1β) and tumor necrosis factor α (TNFα) for 3 weeks to simulate OA-related matrix degradation. In addition, an existing mathematical model based on partial differential equations was optimized and transferred to the underlying settings to simulate distribution of IL-1β, type II collagen degradation and cell number reduction. By combining in vitro and in silico methods, we aim to develop a valid, efficient alternative approach to examine and predict disease progression and effects of new therapeutics. Y1 - 2020 U6 - https://doi.org/10.1088/1758-5090/aba08f VL - 12 IS - 4 ER - TY - JOUR A1 - Erlekam, Franziska A1 - Igde, Sinaida A1 - Röblitz, Susanna A1 - Hartmann, Laura A1 - Weber, Marcus T1 - Modeling of Multivalent Ligand-Receptor Binding Measured by kinITC JF - Computation N2 - In addition to the conventional Isothermal Titration Calorimetry (ITC), kinetic ITC (kinITC) not only gains thermodynamic information, but also kinetic data from a biochemical binding process. Moreover, kinITC gives insights into reactions consisting of two separate kinetic steps, such as protein folding or sequential binding processes. The ITC method alone cannot deliver kinetic parameters, especially not for multivalent bindings. This paper describes how to solve the problem using kinITC and an invariant subspace projection. The algorithm is tested for multivalent systems with different valencies. Y1 - 2019 U6 - https://doi.org/10.3390/computation7030046 VL - 7 IS - 3 SP - 46 ER - TY - JOUR A1 - Omari, Mohamed A1 - Lange, Alexander A1 - Plöntzke, Julia A1 - Röblitz, Susanna T1 - Model-based exploration of the impact of glucose metabolism on the estrous cycle dynamics in dairy cows JF - Biology Direct Y1 - 2019 U6 - https://doi.org/10.1186/s13062-019-0256-7 VL - 15 ER - TY - JOUR A1 - Fischer, Sophie A1 - Ehrig, Rainald A1 - Schäfer, Stefan A1 - Tronci, Enrico A1 - Mancini, Toni A1 - Egli, Marcel A1 - Ille, Fabian A1 - Krüger, Tillmann H. C. A1 - Leeners, Brigitte A1 - Röblitz, Susanna T1 - Mathematical Modelling and Simulation Provides Evidence for New Strategies of Ovarian Stimulation JF - Frontiers in Endocrinology N2 - New approaches to ovarian stimulation protocols, such as luteal start, random start or double stimulation, allow for flexibility in ovarian stimulation at different phases of the menstrual cycle which is especially useful when time for assisted reproductive technology is limited, e.g. for emergency fertility preservation in cancer patients. It has been proposed that the success of these methods is based on the continuous growth of multiple cohorts ("waves") of follicles throughout the menstrual cycle which leads to the availability of ovarian follicles for ovarian controlled stimulation at several time points. Though several preliminary studies have been published, their scientific evidence has not been considered as being strong enough to integrate these results into routine clinical practice. This work aims at adding further scientific evidence about the efficiency of variable-start protocols and underpinning the theory of follicular waves by using mathematical modelling and numerical simulations. For this purpose, we have modified and coupled two previously published models, one describing the time course of hormones and one describing competitive follicular growth in a normal menstrual cycle. The coupled model is used to test stimulation protocols in silico. Simulation results show the occurrence of follicles in a wave-like manner during a normal menstrual cycle and qualitatively predict the outcome of ovarian stimulation initiated at different time points of the menstrual cycle. Y1 - 2021 U6 - https://doi.org/10.3389/fendo.2021.613048 VL - 12 ER - TY - GEN A1 - Omari, Mohamed A1 - Lange, Alexander A1 - Plöntzke, Julia A1 - Röblitz, Susanna T1 - A Mathematical Model for the Influence of Glucose-Insulin Dynamics on the Estrous Cycle in Dairy Cows N2 - Nutrition plays a crucial role in regulating reproductive hormones and follicular development in cattle. This is visible particularly during the time of negative energy balance at the onset of milk production after calving. Here, elongated periods of anovulation have been observed, resulting from alterations in luteiniz- ing hormone concentrations, likely caused by lower glucose and insulin concen- trations in the blood. The mechanisms that result in a reduced fertility are not completely understood, although a close relationship to the glucose-insulin metabolism is widely supported. Following this idea, a mathematical model of the hormonal network combining reproductive hormones and hormones that are coupled to the glucose compartments within the body of the cow was developed. The model is built on ordinary differential equations and relies on previously introduced models on the bovine estrous cycle and the glucose-insulin dynam- ics. Necessary modifications and coupling mechanisms are thoroughly discussed. Depending on the composition and the amount of food, in particular the glu- cose content in the dry matter, the model quantifies reproductive hormones and follicular development over time. Simulation results for different nutritional regimes in lactating and non-lactating dairy cows are examined and compared with experimental studies. Regarding its applicability, this work is an early attempt towards developing in silico feeding strategies and may eventually help refining and reducing animal experiments. T3 - ZIB-Report - 19-22 Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-73475 SN - 1438-0064 ER - TY - GEN A1 - Plöntzke, Julia A1 - Berg, Mascha A1 - Röblitz, Susanna T1 - A mathematical modelling approach to the insight of dynamic networks: Potassium homeostasis and glucose-insulin in dairy cows T2 - ADSA 2018 Annual Meeting N2 - Lactating dairy cows require a particular composition of nutritional ingredients depending on their production status. The optimal supply of energy and minerals in diet, one of them potassium, is indispensable for the prevention of disbalances like hypokalemia or hypoglycaemia. Potassium balance in cows is the result of potassium intake, distribution in the organism, and excretion, and closely interacts with glucose and electrolyte metabolism, in which postpartum veterinary treatments frequently intervene. We present a mechanistic, dynamic model for potassium balance together with a glucose insulin model in non-lactating and lactating dairy cows based on ordinary differential equations. Parameter values were obtained from data of a clinical trial as well as from literature. To verify the mechanistic functioning of the model, we validate the model by comparing simulation outcomes with clinical study findings. Furthermore we perform numerical experiments and compare them with expected behaviour according to mechanistic knowledge. The results give insight into the dynamic behaviour of the network and open the way for further open questions and hypothesis to be tested. Y1 - 2018 UR - https://www.adsa.org/Portals/0/SiteContent/Docs/Meetings/PastMeetings/Annual/2018/223.pdf VL - 101 Suppl. 2 SP - 226 EP - 227 PB - American Dairy Science Association ET - Journal of Dairy Science ER - TY - JOUR A1 - Schwieger, Robert A1 - Siebert, Heike A1 - Röblitz, Susanna T1 - Correspondence of Trap Spaces in Different Models of Bioregulatory Networks JF - SIAM Journal on Applied Dynamical Systems N2 - Mathematical models for bioregulatory networks can be based on different formalisms, depending on the quality of available data and the research question to be answered. Discrete boolean models can be constructed based on qualitative data, which are frequently available. On the other hand, continuous models in terms of ordinary differential equations (ODEs) can incorporate time-series data and give more detailed insight into the dynamics of the underlying system. A few years ago, a method based on multivariate polynomial interpolation and Hill functions has been developed for an automatic conversion of boolean models to systems of ordinary differential equations. This method is frequently used by modellers in systems biology today, but there are only a few results available about the conservation of mathematical structures and properties across the formalisms. Here, we consider subsets of the phase space where some components stay fixed, called trap spaces, and demonstrate how boolean trap spaces can be linked to invariant sets in the continuous state space. This knowledge is of practical relevance since finding trap spaces in the boolean setting, which is relatively easy, allows for the construction of reduced ODE models. Y1 - 2018 U6 - https://doi.org/https://doi.org/10.1137/17M1155302 VL - 17 IS - 2 SP - 1742 EP - 1765 ER - TY - JOUR A1 - Lange, Alexander A1 - Schwieger, Robert A1 - Plöntzke, Julia A1 - Schäfer, Stefan A1 - Röblitz, Susanna T1 - Follicular competition in cows: the selection of dominant follicles as a synergistic effect JF - Journal of Mathematical Biology N2 - The reproductive cycle of mono-ovulatory species such as cows or humans is known to show two or more waves of follicular growth and decline between two successive ovulations. Within each wave, there is one dominant follicle escorted by subordinate follicles of varying number. Under the surge of the luteinizing hormone a growing dominant follicle ovulates. Rarely the number of ovulating follicles exceeds one. In the biological literature, the change of hormonal concentrations and individually varying numbers of follicular receptors are made responsible for the selection of exactly one dominant follicle, yet a clear cause has not been identified. In this paper, we suggest a synergistic explanation based on competition, formulated by a parsimoniously defined system of ordinary differential equations (ODEs) that quantifies the time evolution of multiple follicles and their competitive interaction during one wave. Not discriminating between follicles, growth and decline are given by fixed rates. Competition is introduced via a growth-suppressing term, equally supported by all follicles. We prove that the number of dominant follicles is determined exclusively by the ratio of follicular growth and competition. This number turns out to be independent of the number of subordinate follicles. The asymptotic behavior of the corresponding dynamical system is investigated rigorously, where we demonstrate that the ω-limit set only contains fixed points. When also including follicular decline, our ODEs perfectly resemble ultrasound data of bovine follicles. Implications for the involved but not explicitly modeled hormones are discussed. Y1 - 2018 U6 - https://doi.org/https://doi.org/10.1007/s00285-018-1284-0 ER - TY - CHAP A1 - Mancini, Toni A1 - Mari, Federico A1 - Massini, Annalisa A1 - Melatti, Igor A1 - Salvo, Ivano A1 - Sinisi, Stefano A1 - Tronci, Enrico A1 - Ehrig, Rainald A1 - Röblitz, Susanna A1 - Leeners, Brigitte T1 - Computing Personalised Treatments through In Silico Clinical Trials. A Case Study on Downregulation in Assisted Reproduction T2 - Experimental Evaluation of Algorithms for Solving Problems with Combinatorial Explosion, Workshop of the Federated Logic Conference 2018 (FLoC 2018), Oxford N2 - In Silico Clinical Trials (ISCT), i.e., clinical experimental campaigns carried out by means of computer simulations, hold the promise to decrease time and cost for the safety and efficacy assessment of pharmacological treatments, reduce the need for animal and human testing, and enable precision medicine. In this paper we present a case study aiming at quantifying, by means of a multi-arm ISCT supervised by intelligent search, the potential impact of precision medicine approaches on a real pharmacological treatment, namely the downregulation phase of a complex clinical protocol for assisted reproduction. Y1 - 2018 U6 - https://doi.org/10.29007/g864 SP - 1 EP - 16 ER - TY - GEN A1 - Berg, Mascha A1 - Plöntzke, Julia A1 - Leonhard-Marek, Sabine A1 - Müller, Kerstin-Elisabeth A1 - Röblitz, Susanna T1 - A dynamic model to simulate potassium balance in dairy cows. N2 - High performing dairy cows require a particular composition of nutritional ingredients, adapted to their individual requirements and depending on their production status. The optimal dimensioning of minerals in the diet, one of them being potassium, is indispensable for the prevention of imbalances. The potassium balance in cows is the result of potassium intake, distribution in the organism, and excretion, it is closely related with the glucose and electrolyte metabolism. In this paper, we present a dynamical model for the potassium balance in lactating and non-lactating dairy cows based on ordinary differential equations. Parameter values are obtained from clinical trial data and from the literature. To verify the consistency of the model, we present simulation outcomes for three different scenarios: potassium balance in (i) non-lactating cows with varying feed intake, (ii) non-lactating cows with varying potassium fraction in the diet, and (iii) lactating cows with varying milk production levels. The results give insights into the short and long term potassium metabolism, providing an important step towards the understanding of the potassium network, the design of prophylactic feed additives, and possible treatment strategies. T3 - ZIB-Report - 17-47 KW - cow, potassium, model, ode Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-64756 SN - 1438-0064 ER - TY - JOUR A1 - Berg, Mascha A1 - Plöntzke, Julia A1 - Siebert, Heike A1 - Röblitz, Susanna T1 - Modelling Oscillatory Patterns in the Bovine Estrous Cycle with Boolean Delay Equations JF - Bulletin of Mathematical Biology N2 - Boolean delay equations (BDEs), with their relatively simple and intuitive mode of modelling, have been used in many research areas including, for example, climate dynamics and earthquake propagation. Their application to biological systems has been scarce and limited to the molecular level. Here, we derive and present two BDE models. One is directly derived from a previously published ordinary differential equation (ODE) model for the bovine estrous cycle, whereas the second model includes a modification of a particular biological mechanism. We not only compare the simulation results from the BDE models with the trajectories of the ODE model, but also validate the BDE models with two additional numerical experiments. One experiment induces a switch in the oscillatory pattern upon changes in the model parameters, and the other simulates the administration of a hormone that is known to shift the estrous cycle in time. The models presented here are the first BDE models for hormonal oscillators, and the first BDE models for drug administration. Even though automatic parameter estimation still remains challenging, our results support the role of BDEs as a framework for the systematic modelling of complex biological oscillators. Y1 - 2021 U6 - https://doi.org/10.1007/s11538-021-00942-z VL - 83 IS - 121 SP - 1 EP - 25 ER - TY - GEN A1 - Klebanov, Ilja A1 - Sikorski, Alexander A1 - Schütte, Christof A1 - Röblitz, Susanna T1 - Empirical Bayes Methods, Reference Priors, Cross Entropy and the EM Algorithm N2 - When estimating a probability density within the empirical Bayes framework, the non-parametric maximum likelihood estimate (NPMLE) usually tends to overfit the data. This issue is usually taken care of by regularization - a penalization term is subtracted from the marginal log-likelihood before the maximization step, so that the estimate favors smooth solutions, resulting in the so-called maximum penalized likelihood estimation (MPLE). The majority of penalizations currently in use are rather arbitrary brute-force solutions, which lack invariance under transformation of the parameters(reparametrization) and measurements. This contradicts the principle that, if the underlying model has several equivalent formulations, the methods of inductive inference should lead to consistent results. Motivated by this principle and using an information-theoretic point of view, we suggest an entropy-based penalization term that guarantees this kind of invariance. The resulting density estimate can be seen as a generalization of reference priors. Using the reference prior as a hyperprior, on the other hand, is argued to be a poor choice for regularization. We also present an insightful connection between the NPMLE, the cross entropy and the principle of minimum discrimination information suggesting another method of inference that contains the doubly-smoothed maximum likelihood estimation as a special case. T3 - ZIB-Report - 16-56 KW - parameter estimation KW - Bayesian inference KW - Bayesian hierarchical modeling KW - hyperparameter KW - hyperprior KW - EM algorithm KW - NPMLE KW - MPLE KW - DS-MLE KW - principle of maximum entropy KW - cross entropy KW - minimum discrimination information KW - reference prior KW - Jeffreys prior Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-61230 SN - 1438-0064 ER - TY - JOUR A1 - Fischer-Holzhausen, Sophie A1 - Röblitz, Susanna T1 - Hormonal regulation of ovarian follicle growth in humans: Model-based exploration of cycle variability and parameter sensitivities JF - Journal of Theoretical Biology N2 - We present a modelling and simulation framework for the dynamics of ovarian follicles and key hormones along the hypothalamic-pituitary-gonadal axis throughout consecutive human menstrual cycles. All simulation results (hormone concentrations and ovarian follicle sizes) are in biological units and can easily be compared to clinical data. The model takes into account variability in follicles' response to stimulating hormones, which introduces variability between cycles. The growth of ovarian follicles in waves is an emergent property in our model simulations and further supports the hypothesis that follicular waves are also present in humans. We use Approximate Bayesian Computation and cluster analysis to construct a population of virtual subjects and to study parameter distributions and sensitivities. The model can be used to compare and optimize treatment protocols for ovarian hyperstimulation, thus potentially forming the integral part of a clinical decision support system in reproductive endocrinology. Y1 - 2022 U6 - https://doi.org/https://doi.org/10.1016/j.jtbi.2022.111150 ER - TY - JOUR A1 - Krüger, Tillmann A1 - Leeners, Brigitte A1 - Tronci, Enrico A1 - Mancini, Toni A1 - Ille, Fabian A1 - Egli, Marcel A1 - Engler, Harald A1 - Röblitz, Susanna A1 - Frieling, Helge A1 - Sinke, Christopher A1 - Jahn, Kirsten T1 - The androgen system across the menstrual cycle: Hormonal, (epi-)genetic and psychometric alterations JF - Physiology and Behavior Y1 - 2022 U6 - https://doi.org/https://doi.org/10.1016/j.physbeh.2022.114034 ER - TY - JOUR A1 - Fischer-Holzhausen, Sophie A1 - Röblitz, Susanna T1 - Mathematical modelling of follicular growth and ovarian stimulation JF - Current Opinion in Endocrine and Metabolic Research Y1 - 2022 U6 - https://doi.org/https://doi.org/10.1016/j.coemr.2022.100385 ER - TY - GEN A1 - Ehrig, Rainald A1 - Dierkes, Thomas A1 - Schäfer, Stefan A1 - Röblitz, Susanna A1 - Tronci, Enrico A1 - Mancini, Toni A1 - Salvo, Ivano A1 - Alimguzhin, Vadim A1 - Mari, Federico A1 - Melatti, Igor A1 - Massini, Annalisa A1 - Krüger, Tillmann H. C. A1 - Egli, Marcel A1 - Ille, Fabian A1 - Leeners, Brigitte T1 - An Integrative Approach for Model Driven Computation of Treatments in Reproductive Medicine N2 - We present an overview of the current status of the European collaborative project PAEON. The challenge of PAEON is to provide specialists in reproduc- tive medicine with a computerised model of the menstrual cycle under normal and various pathological conditions, which will allow them to get further in- sight in fertility dynamics. This model also enables the simulation of treatment protocols, which were used within in vitro fertilization. By the definition of virtual patients through biologically admissible parametrizations our approach allows not only the evaluation of a given treatment strategy in silico, but also the design and optimization of such protocols. Once a protocol is formalized in the virtual hospital, the success can be controlled by a treatment execution monitor, which works then as a clinical decision support system. All these tools will be combined in a virtual hospital environment, enabling the access to the PAEON services through the web. T3 - ZIB-Report - 16-04 Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-57103 SN - 1438-0064 ER - TY - CHAP A1 - Ehrig, Rainald A1 - Dierkes, Thomas A1 - Schäfer, Stefan A1 - Röblitz, Susanna A1 - Tronci, Enrico A1 - Mancini, Toni A1 - Salvo, Ivano A1 - Alimguzhin, Vadim A1 - Mari, Federico A1 - Melatti, Igor A1 - Massini, Annalisa A1 - Krüger, Tillmann H. C. A1 - Egli, Marcel A1 - Ille, Fabian A1 - Leeners, Brigitte T1 - An integrative approach for model driven computation of treatments in reproductive medicine T2 - BIOMAT - Proceedings of the 15th International Symposium on Mathematical and Computational Biology, Rorkee, India Y1 - 2016 ER - TY - GEN A1 - Dierkes, Thomas A1 - Wade, Moritz A1 - Nowak, Ulrich A1 - Röblitz, Susanna T1 - BioPARKIN - Biology-related Parameter Identification in Large Kinetic Networks N2 - Modelling, parameter identification, and simulation play an important rôle in Systems Biology. In recent years, various software packages have been established for scientific use in both licencing types, open source as well as commercial. Many of these codes are based on inefficient and mathematically outdated algorithms. By introducing the package BioPARKIN recently developed at ZIB, we want to improve this situation significantly. The development of the software BioPARKIN involves long standing mathematical ideas that, however, have not yet entered the field of Systems Biology, as well as new ideas and tools that are particularly important for the analysis of the dynamics of biological networks. BioPARKIN originates from the package PARKIN, written by P.Deuflhard and U.Nowak, that has been applied successfully for parameter identification in physical chemistry for many years. T3 - ZIB-Report - 11-15 KW - systems biology KW - ordinary differential equations KW - sensitivity analysis KW - parameter identification KW - affine invariant Gauss-Newton algorithm KW - numerical library KW - graphical user interface Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-12706 ER - TY - CHAP A1 - Mancini, Toni A1 - Salvo, Ivano A1 - Mari, Federico A1 - Melatti, Igor A1 - Massini, Annalisa A1 - Sinisi, Stefano A1 - Tronci, Enrico A1 - Davi, Francesco A1 - Dierkes, Thomas A1 - Ehrig, Rainald A1 - Röblitz, Susanna A1 - Leeners, Brigitte A1 - Krüger, Tillmann H. C. A1 - Egli, Marcel A1 - Ille, Fabian T1 - Patient-Specific Models from Inter-Patient Biological Models and Clinical Records T2 - Proceedings of the 14th Conference on Formal Methods in Computer-Aided Design (FMCAD 2014) Y1 - 2014 SP - 207 EP - 214 ER - TY - JOUR A1 - Dierkes, Thomas A1 - Röblitz, Susanna A1 - Wade, Moritz A1 - Deuflhard, Peter T1 - Parameter identification in large kinetic networks with BioPARKIN JF - CoRR - Computing Research Repository - arXiv.org Y1 - 2013 UR - arxiv.org/abs/1303.4928v2 VL - abs/1303.4928 ER - TY - JOUR A1 - Plöntzke, Julia A1 - Berg, Mascha A1 - Ehrig, Rainald A1 - Leonhard-Marek, Sabine A1 - Müller, Kerstin-Elisabeth A1 - Röblitz, Susanna T1 - Model Based Exploration of Hypokalemia in Dairy Cows JF - Scientific Reports N2 - Hypokalemia, characterized by too low serum potassium levels, is a severe mineral disorder which can be life threatening. It is increasingly diagnosed in veterinarian healthcare and a topic of ongoing research. In this paper, we explore the different originating conditions of hypokalemia: reduced potassium intake, increased excretion, acid-base disturbances, or increased insulin, by using a dynamic mathematical model for potassium balance in non-lactating and lactating cows. Simulations are compared with literature. The results give insights into the network dynamics and point to scenarios on which experimental effort should be focused. Application of mathematical models can assist in experimental planning as well as the reduction, refinement and replacement of animal experiments. Y1 - 2022 U6 - https://doi.org/10.1038/s41598-022-22596-0 VL - 12, 19781 ER - TY - GEN A1 - Fackeldey, Konstantin A1 - Röblitz, Susanna A1 - Scharkoi, Olga A1 - Weber, Marcus T1 - Soft Versus Hard Metastable Conformations in Molecular Simulations N2 - Particle methods have become indispensible in conformation dynamics to compute transition rates in protein folding, binding processes and molecular design, to mention a few. Conformation dynamics requires at a decomposition of a molecule's position space into metastable conformations. In this paper, we show how this decomposition can be obtained via the design of either ``soft'' or ``hard'' molecular conformations. We show, that the soft approach results in a larger metastabilitiy of the decomposition and is thus more advantegous. This is illustrated by a simulation of Alanine Dipeptide. T3 - ZIB-Report - 11-27 KW - Proteins, Conformation Space, Meshfree Methods Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-13189 ER - TY - JOUR A1 - Haack, Fiete A1 - Fackeldey, Konstantin A1 - Röblitz, Susanna A1 - Scharkoi, Olga A1 - Weber, Marcus A1 - Schmidt, Burkhard T1 - Adaptive spectral clustering with application to tripeptide conformation analysis JF - The Journal of Chemical Physics Y1 - 2013 U6 - https://doi.org/10.1063/1.4830409 VL - 139 SP - 110 EP - 194 ER - TY - JOUR A1 - Reuter, Bernhard A1 - Weber, Marcus A1 - Fackeldey, Konstantin A1 - Röblitz, Susanna A1 - Garcia, Martin E. T1 - Generalized Markov State Modeling Method for Nonequilibrium Biomolecular Dynamics: Exemplified on Amyloid β Conformational Dynamics Driven by an Oscillating Electric Field JF - Journal of Chemical Theory and Computation N2 - Markov state models (MSMs) have received an unabated increase in popularity in recent years, as they are very well suited for the identification and analysis of metastable states and related kinetics. However, the state-of-the-art Markov state modeling methods and tools enforce the fulfillment of a detailed balance condition, restricting their applicability to equilibrium MSMs. To date, they are unsuitable to deal with general dominant data structures including cyclic processes, which are essentially associated with nonequilibrium systems. To overcome this limitation, we developed a generalization of the common robust Perron Cluster Cluster Analysis (PCCA+) method, termed generalized PCCA (G-PCCA). This method handles equilibrium and nonequilibrium simulation data, utilizing Schur vectors instead of eigenvectors. G-PCCA is not limited to the detection of metastable states but enables the identification of dominant structures in a general sense, unraveling cyclic processes. This is exemplified by application of G-PCCA on nonequilibrium molecular dynamics data of the Amyloid β (1−40) peptide, periodically driven by an oscillating electric field. Y1 - 2018 U6 - https://doi.org/10.1021/acs.jctc.8b00079 VL - 14 IS - 7 SP - 3579 EP - 3594 ER - TY - JOUR A1 - Fackeldey, Konstantin A1 - Röblitz, Susanna A1 - Scharkoi, O. A1 - Weber, Marcus T1 - Soft Versus Hard Metastable Conformations in Molecular Simulations JF - Particle Methods II, Fundamentals and Applications, Barcelona, Spain 26-28 Oct. 2011, E. Onate and D.R.J. Owen (eds.) Y1 - 2011 SP - 899 EP - 909 ER - TY - JOUR A1 - Frank, Anna-Simone A1 - Sikorski, Alexander A1 - Röblitz, Susanna T1 - Spectral clustering of Markov chain transition matrices with complex eigenvalues JF - Journal of Computational and Applied Mathematics N2 - The Robust Perron Cluster Analysis (PCCA+) has become a popular spectral clustering algorithm for coarse-graining transition matrices of nearly decomposable Markov chains with transition states. Originally developed for reversible Markov chains, the algorithm only worked for transition matrices with real eigenvalues. In this paper, we therefore extend the theoretical framework of PCCA+ to Markov chains with a complex eigen-decomposition. We show that by replacing a complex conjugate pair of eigenvectors by their real and imaginary components, a real representation of the same subspace is obtained, which is suitable for the cluster analysis. We show that our approach leads to the same results as the generalized PCCA+ (GPCCA), which replaces the complex eigen-decomposition by a conceptually more difficult real Schur decomposition. We apply the method on non-reversible Markov chains, including circular chains, and demonstrate its efficiency compared to GPCCA. The experiments are performed in the Matlab programming language and codes are provided. Y1 - 2024 U6 - https://doi.org/10.1016/j.cam.2024.115791 VL - 444 SP - 115791 ER -