@misc{MeyerMollCordesetal.2006, author = {Meyer, Holger and Moll, Sebastian and Cordes, Frank and Weber, Marcus}, title = {ConFlow? A new space-based Application for complete Conformational Analysis of Molecules}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-9258}, number = {06-31}, year = {2006}, abstract = {Biochemical interactions are determined by the 3D-structure of the involved components - thus the identification of conformations is a key for many applications in rational drug design. {\sf ConFlow} is a new multilevel approach to conformational analysis with main focus on completeness in investigation of conformational space. In contrast to known conformational analysis, the starting point for design is a space-based description of conformational areas. A tight integration of sampling and analysis leads to an identification of conformational areas simultaneously during sampling. An incremental decomposition of high-dimensional conformational space is used to guide the analysis. A new concept for the description of conformations and their path connected components based on convex hulls and {\em Hypercubes}is developed. The first results of the {\sf ConFlow} application constitute a 'proof of concept' and are further more highly encouraging. In comparison to conventional industrial applications, {\sf ConFlow} achieves higher accuracy and a specified degree of completeness with comparable effort.}, language = {en} } @misc{KubeWeber2006, author = {Kube, Susanna and Weber, Marcus}, title = {Coarse Grained Molecular Kinetics}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-9299}, number = {06-35}, year = {2006}, abstract = {The dynamic behavior of molecules can often be described by Markov processes. From computational molecular simulations one can derive transition rates or transition probabilities between subsets of the discretized conformational space. On the basis of this dynamic information, the spatial subsets are combined into a small number of so-called metastable molecular conformations. This is done by clustering methods like the Robust Perron Cluster Analysis (PCCA+). Up to now it is an open question how this coarse graining in space can be transformed to a coarse graining of the Markov chain while preserving the essential dynamic information. In the following article we aim at a consistent coarse graining of transition probabilities or rates on the basis of metastable conformations such that important physical and mathematical relations are preserved. This approach is new because PCCA+ computes molecular conformations as linear combinations of the dominant eigenvectors of the transition matrix which does not hold for other clustering methods.}, language = {en} } @misc{KubeLasserWeber2007, author = {Kube, Susanna and Lasser, Caroline and Weber, Marcus}, title = {Monte Carlo sampling of Wigner functions and surface hopping quantum dynamics}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-9604}, number = {07-17}, year = {2007}, abstract = {Wigner transformation provides a one-to-one correspondence between functions on position space (wave functions) and functions on phase space (Wigner functions). Weighted integrals of Wigner functions yield quadratic quantities of wave functions like position and momentum densities or expectation values. For molecular quantum systems, suitably modified classical transport of Wigner functions provides an asymptotic approximation of the dynamics in the high energy regime. The article addresses the computation of Wigner functions by Monte Carlo quadrature. An ad aption of the Metropolis algorithm for the approximation of signed measures with disconnected support is systematically tested in combination with a surface hopping algorithm for non-adiabatic quantum dynamics. The numerical experiments give expectation values and level populations with an error of two to three percent, which agrees with the theoretically expected accuracy.}, language = {en} } @article{ThiesSunkaraRayetal.2023, author = {Thies, Arne and Sunkara, Vikram and Ray, Sourav and Wulkow, Hanna and Celik, M. {\"O}zg{\"u}r and Yerg{\"o}z, Fatih and Sch{\"u}tte, Christof and Stein, Christoph and Weber, Marcus and Winkelmann, Stefanie}, title = {Modelling altered signalling of G-protein coupled receptors in inflamed environment to advance drug design}, volume = {13}, journal = {Scientific Reports}, number = {607}, doi = {10.1038/s41598-023-27699-w}, year = {2023}, abstract = {We previously reported the successful design, synthesis and testing of the prototype opioid painkiller NFEPP that does not elicit adverse side effects. The design process of NFEPP was based on mathematical modelling of extracellular interactions between G-protein coupled receptors (GPCRs) and ligands, recognizing that GPCRs function differently under pathological versus healthy conditions. We now present an additional and novel stochastic model of GPCR function that includes intracellular dissociation of G-protein subunits and modulation of plasma membrane calcium channels and their dependence on parameters of inflamed and healthy tissue (pH, radicals). The model is validated against in vitro experimental data for the ligands NFEPP and fentanyl at different pH values and radical concentrations. We observe markedly reduced binding affinity and calcium channel inhibition for NFEPP at normal pH compared to lower pH, in contrast to the effect of fentanyl. For increasing radical concentrations, we find enhanced constitutive G-protein activation but reduced ligand binding affinity. Assessing the different effects, the results suggest that, compared to radicals, low pH is a more important determinant of overall GPCR function in an inflamed environment. Future drug design efforts should take this into account.}, language = {en} } @misc{RayThiesSunkaraetal.2021, author = {Ray, Sourav and Thies, Arne and Sunkara, Vikram and Wulkow, Hanna and Celik, {\"O}zg{\"u}r and Yerg{\"o}z, Fatih and Sch{\"u}tte, Christof and Stein, Christoph and Weber, Marcus and Winkelmann, Stefanie}, title = {Modelling altered signalling of G-protein coupled receptors in inflamed environment to advance drug design}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-82797}, year = {2021}, abstract = {Initiated by mathematical modelling of extracellular interactions between G-protein coupled receptors (GPCRs) and ligands in normal versus diseased (inflamed) environments, we previously reported the successful design, synthesis and testing of the prototype opioid painkiller NFEPP that does not elicit adverse side effects. Uniquely, this design recognised that GPCRs function differently under pathological versus healthy conditions. We now present a novel stochastic model of GPCR function that includes intracellular dissociation of G-protein subunits and modulation of plasma membrane calcium channels associated with parameters of inflamed tissue (pH, radicals). By means of molecular dynamics simulations, we also assessed qualitative changes of the reaction rates due to additional disulfide bridges inside the GPCR binding pocket and used these rates for stochastic simulations of the corresponding reaction jump process. The modelling results were validated with in vitro experiments measuring calcium currents and G-protein activation. We found markedly reduced G-protein dissociation and calcium channel inhibition induced by NFEPP at normal pH, and enhanced constitutive G-protein activation but lower probability of ligand binding with increasing radical concentrations. These results suggest that, compared to radicals, low pH is a more important determinant of overall GPCR function in an inflamed environment. Future drug design efforts should take this into account.}, language = {en} } @article{HeidaSikorskiWeber2022, author = {Heida, Martin and Sikorski, Alexander and Weber, Marcus}, title = {Consistency and order 1 convergence of cell-centered finite volume discretizations of degenerate elliptic problems in any space dimension}, journal = {SIAM Journal on Numerical Analysis}, doi = {10.20347/WIAS.PREPRINT.2913}, year = {2022}, abstract = {We study consistency of cell-centered finite difference methods for elliptic equations with degenerate coefficients in any space dimension \$d \geq 2\$. This results in order of convergence estimates in the natural weighted energy norm and in the weighted discrete \$L^2\$-norm on admissible meshes. The cells of meshes under consideration may be very irregular in size. We particularly allow the size of certain cells to remain bounded from below even in the asymptotic limit. For uniform meshes we show that the order of convergence is at least 1 in the energy semi-norm, provided the discrete and continuous solutions exist and the continuous solution has \$H^2\$ regularity.}, language = {en} } @article{FackeldeyGorgullaWeber2021, author = {Fackeldey, Konstantin and Gorgulla, Christoph and Weber, Marcus}, title = {Neue Medikamente dank Supercomputern}, journal = {Spektrum der Wissenschaft}, number = {11}, pages = {40 -- 46}, year = {2021}, abstract = {Die aktuelle Pandemie verdeutlicht, wie wichtig es ist, rasch geeignete Arzneimittel zu finden. In Computer­simulationen gelingt das erheblich schneller als im Labor. Gegen das Coronavirus ließen sich auf diese Weise bereits Wirkstoffkandidaten identifizieren.}, language = {de} } @article{BirkRaharinirinaFackeldeyetal.2021, author = {Birk, Ralph and Raharinirina, N. Alexia and Fackeldey, Konstantin and Richter, Tonio Sebastian and Weber, Marcus}, title = {Inferring cultural and social processes based on patterns of statistical relationships between Synodal texts}, year = {2021}, abstract = {In this paper, we explore the relationship patterns between Ancient Egyptian texts of the corpus ``Synodal decrees'', which are originating between 243 and 185 BCE, during the Ptolemaic period. Particularly, we are interested in analyzing the grammatical features of the different texts. Conventional data analysis methods such as correspondence Analysis are very useful to explore the patterns of statistical interdependence between categories of variables. However, it is based on a PCA-like dimension-reduction method and turned out to be unsuitable for our dataset due to the high dimensionality of our data representations. Additionally, the similarity between pairs of texts and pairs of grammatical features is observed through the distance between their representation, but the degree of association between a particular grammatical feature and a text is not. Here, we applied a qualitative Euclidean embedding method that provides a new Euclidean representation of the categories of variables. This new representation of the categories is constructed in such a way that all the patterns of statistical interdependence, similarity, and association, are seen through the Euclidean distance between them. Nevertheless, the PCA-like dimension-reduction method also performed poorly on our new representation. Therefore, we obtained a two-dimensional visualization using non-linear methods such UMAP or t-SNE. Although these dimension-reduction methods reduced the interpretability of interpoint distances, we were still able to identify important similarity patterns between the Synodal text as well as their association patterns with the grammatical features.}, language = {en} } @article{SechiWeberHeyne2021, author = {Sechi, Renata and Weber, Marcus and Heyne, Karsten}, title = {MSM building and projection for the analysis of time-resolved spectra}, journal = {Proceedings in Applied Mathematics and Mechanics}, number = {21}, doi = {https://doi.org/10.1002/pamm.202100102}, url = {http://nbn-resolving.de/https://onlinelibrary.wiley.com/doi/10.1002/pamm.202100102}, pages = {e202100102}, year = {2021}, abstract = {Understanding the kinetics between the components of time-resolved spectra is a crucial step in the study of photo-activatedprocesses. However, modeling the kinetics requires usually some a priori knowledge about the system. In our approach, webuild a Markov State Model (MSM) from the spectral data, and obtain a Koopman transition matrix K(t). With genPCCA,an invariant subspace projection, we project the process into its metastable components. The result of the application of gen-PCCA is a transition matrix Kc(t), from which we can read the transition probability between the metastable components of the reaction. We discuss the application of this analysis method to the transient absorption spectrum of brominated Al-corrole}, language = {en} } @article{ChewleWeberEmmerling2023, author = {Chewle, Surahit and Weber, Marcus and Emmerling, Franziska}, title = {Revealing Kinetics of Paracetamol Crystallization Using Time-Resolved Raman Spectroscopy, Orthogonal Time-Lapse Photography, and Non-Negative Matrix Factorization (OSANO)}, volume = {23}, journal = {Crystal Growth \& Design}, number = {9}, edition = {Crystal Growth \& Design}, publisher = {American Chemical Society}, doi = {10.1021/acs.cgd.3c00617}, pages = {6737 -- 6746}, year = {2023}, abstract = {Crystallization is a complex phenomenon with far-reaching implications for the production and formulation of active pharmaceutical ingredients. Understanding this process is critical for achieving control over key physicochemical properties that can affect, for example, the bioavailability and stability of a drug. In this study, we were able to reveal intricate and diverse dynamics of the formation of metastable intermediates of paracetamol crystallization varying with the choice of solvent. We demonstrate the efficacy of our novel approach utilizing an objective function-based non-negative matrix factorization technique for the analysis of time-resolved Raman spectroscopy data, in conjunction with time-lapse photography. Furthermore, we emphasize the crucial importance of integrating Raman spectroscopy with supplementary experimental instrumentation for the mathematical analysis of the obtained spectra.}, language = {en} } @article{MarcoMarcus2025, author = {Marco, Reidelbach and Marcus, Weber}, title = {From Documentation to Discovery - MaRDMO and the MaRDI Ecosystem}, journal = {2nd Conference on Research Data Infrastructure (CoRDI)}, doi = {10.5281/zenodo.16735889}, year = {2025}, language = {en} } @inproceedings{ReidelbachSchemberaWeber2024, author = {Reidelbach, Marco and Schembera, Bj{\"o}rn and Weber, Marcus}, title = {Towards a FAIR Documentation of Workflows and Models in Applied Mathematics}, volume = {14749}, booktitle = {Mathematical Software - ICMS 2024}, edition = {Lecture Notes in Computer Science}, publisher = {Springer Nature Switzerland}, address = {Cham}, arxiv = {http://arxiv.org/abs/2403.17778}, doi = {10.1007/978-3-031-64529-7_27}, pages = {254 -- 262}, year = {2024}, abstract = {Modeling-Simulation-Optimization workflows play a fundamental role in applied mathematics. The Mathematical Research Data Initiative, MaRDI, responded to this by developing a FAIR and machine-interpretable template for a comprehensive documentation of such workflows. MaRDMO, a Plugin for the Research Data Management Organiser, enables scientists from diverse fields to document and publish their workflows on the MaRDI Portal seamlessly using the MaRDI template. Central to these workflows are mathematical models. MaRDI addresses them with the MathModDB ontology, offering a structured formal model description. Here, we showcase the interaction between MaRDMO and the MathModDB Knowledge Graph through an algebraic modeling workflow from the Digital Humanities. This demonstration underscores the versatility of both services beyond their original numerical domain.}, language = {en} } @article{RaharinirinaSunkaravonKleistetal.2024, author = {Raharinirina, Nomenjanahary Alexia and Sunkara, Vikram and von Kleist, Max and Fackeldey, Konstantin and Weber, Marcus}, title = {Multi-Input data ASsembly for joint Analysis (MIASA): A framework for the joint analysis of disjoint sets of variables}, volume = {19}, journal = {PLOS ONE}, number = {5}, publisher = {Public Library of Science}, doi = {10.1371/journal.pone.0302425}, pages = {26}, year = {2024}, language = {en} } @article{ReidelbachFerrerWeber2023, author = {Reidelbach, Marco and Ferrer, Eloi and Weber, Marcus}, title = {MaRDMO Plugin}, volume = {1}, journal = {Proceedings of the Conference on Research Data Infrastructure}, publisher = {TIB Open Publishing}, issn = {2941-296X}, doi = {10.52825/cordi.v1i.254}, year = {2023}, abstract = {MaRDMO, a plugin for the Research Data Management Organiser, was developed in the Mathematical Research Data Initiative to document interdisciplinary workflows using a standardised scheme. Interdisciplinary workflows recorded this way are published directly on the MaRDI portal. In addition, central information is integrated into the MaRDI knowledge graph. Next to the documentation, MaRDMO offers the possibility to retrieve existing interdisciplinary workflows from the MaRDI Knowledge Graph to allow the reproduction of the initial work and to provide scientists with new researchimpulses. Thus, MaRDMO creates a community-driven knowledge loop that could help to overcome the replication crisis.}, language = {en} } @article{KlasseWeber2025, author = {Klasse, Sarah Michelle and Weber, Marcus}, title = {A Case Study of the Use of Logical Data Analysis in the Workmen's Village in Tell el-Amarna, Egypt}, volume = {66}, journal = {it - Information Technology}, number = {6}, publisher = {de Gruyter}, doi = {10.1515/itit-2024-0088}, pages = {256 -- 268}, year = {2025}, abstract = {The search for patterns in archaeological contexts is crucial to understanding ancient civilizations, including the organization of their societies and economies. This study focuses on the interdisciplinary application of mathematical logical data analysis to Egyptology, particularly examining food production in the Workmen's Village of Amarna. The introduction provides historical context about Amarna, founded by Pharaoh Akhenaten, and outlines the settlement structure and excavations of the Workmen's Village. The research utilizes archaeobotanical and archaeological data, along with digital humanities tools, to identify activity areas related to food production. The findings reveal spatial patterns in the distribution of tools such as mortars, pestles, and ovens, providing insights into the self-sufficiency and daily life of the villagers. Logical data analysis helps in finding rules and patterns in a dataset which are not visible otherwise, predicting the distribution of unexcavated areas and thus in this case enhancing our understanding of ancient Egyptian food production practices in state- planned settlements.}, language = {en} } @article{DonatiChewleStPierreetal.2025, author = {Donati, Luca and Chewle, Surahit and St. Pierre, Dominik and Natarajan, Vijay and Weber, Marcus}, title = {Topological analysis reveals multiple pathways in molecular dynamics}, volume = {21}, journal = {Journal of Chemical Theory and Computation}, number = {20}, address = {J. Chem. Theory Comput.}, arxiv = {http://arxiv.org/abs/2412.20580}, doi = {10.1021/acs.jctc.5c00819}, pages = {10385 -- 10397}, year = {2025}, abstract = {Molecular Dynamics simulations are indispensable tools for comprehending the dynamic behavior of biomolecules, yet extracting meaningful molecular pathways from these simulations remains challenging due to the vast amount of high dimensional data. In this work, we present Molecular Kinetics via Topology (MoKiTo), a novel approach that combines the ISOKANN algorithm to determine the membership function of a molecular system with a topological analysis tool inspired by the Mapper algorithm. Our strategy efficiently identifies and characterizes distinct molecular pathways, enabling the detection and visualization of critical conformational transitions and rare events. This method offers deeper insights into molecular mechanisms, facilitating the design of targeted interventions in drug discovery and protein engineering.}, language = {en} } @article{BauerWeberDiehlWieseneckeretal.2021, author = {Bauer, Wolfgang and Weber, Marcus and Diehl-Wiesenecker, Eva and Galtung, Noa and Prpic, Monika and Somasundaram, Rajan and Tauber, Rudolf and Schwenk, Jochen and Micke, Patrick and Kappert, Kai}, title = {Plasma Proteome Fingerprints Reveal Distinctiveness and Clinical Outcome of SARS-CoV-2 Infection}, volume = {13}, journal = {Viruses}, number = {12}, doi = {10.3390/v13122456}, pages = {2456}, year = {2021}, abstract = {We evaluated how plasma proteomic signatures in patients with suspected COVID-19 can unravel the pathophysiology, and determine kinetics and clinical outcome of the infection. We identified distinct plasma proteins linked to the presence and course of COVID-19. These plasma proteomic findings may translate to a protein fingerprint, helping to assist clinical management decisions.}, language = {en} } @article{KontouWalterAlkaetal.2023, author = {Kontou, Eftychia E. and Walter, Axel and Alka, Oliver and Pfeuffer, Julianus and Sachsenberg, Timo and Mohite, Omkar and Nuhamunanda, Matin and Kohlbacher, Oliver and Weber, Tilmann}, title = {UmetaFlow: An untargeted metabolomics workflow for high-throughput data processing and analysis}, volume = {15}, journal = {Journal of Cheminformatics}, doi = {10.1186/s13321-023-00724-w}, year = {2023}, abstract = {Metabolomics experiments generate highly complex datasets, which are time and work-intensive, sometimes even error-prone if inspected manually. Therefore, new methods for automated, fast, reproducible, and accurate data processing and dereplication are required. Here, we present UmetaFlow, a computational workflow for untargeted metabolomics that combines algorithms for data pre-processing, spectral matching, molecular formula and structural predictions, and an integration to the GNPS workflows Feature-Based Molecular Networking and Ion Identity Molecular Networking for downstream analysis. UmetaFlow is implemented as a Snakemake workflow, making it easy to use, scalable, and reproducible. For more interactive computing, visualization, as well as development, the workflow is also implemented in Jupyter notebooks using the Python programming language and a set of Python bindings to the OpenMS algorithms (pyOpenMS). Finally, UmetaFlow is also offered as a web-based Graphical User Interface for parameter optimization and processing of smaller-sized datasets. UmetaFlow was validated with in-house LC-MS/MS datasets of actinomycetes producing known secondary metabolites, as well as commercial standards, and it detected all expected features and accurately annotated 76\% of the molecular formulas and 65\% of the structures. As a more generic validation, the publicly available MTBLS733 and MTBLS736 datasets were used for benchmarking, and UmetaFlow detected more than 90\% of all ground truth features and performed exceptionally well in quantification and discriminating marker selection.}, language = {en} } @article{DonatiWeber2023, author = {Donati, Luca and Weber, Marcus}, title = {Efficient Estimation of Transition Rates as Functions of pH}, volume = {23}, journal = {Proceedings in Applied Mathematics \& Mechanics}, arxiv = {http://arxiv.org/abs/2307.04439}, doi = {10.1002/pamm.202300264}, year = {2023}, abstract = {Extracting the kinetic properties of a system whose dynamics depend on the pH of the environment with which it exchanges energy and atoms requires sampling the Grand Canonical Ensemble. As an alternative, we present a novel strategy that requires simulating only the most recurrent Canonical Ensembles that compose the Grand Canonical Ensemble. The simulations are used to estimate the Gran Canonical distribution for a specific pH value by reweighting and to construct the transition rate matrix by discretizing the Fokker-Planck equation by Square Root Approximation and robust Perron Cluster Cluster Analysis. As an application, we have studied the tripeptide Ala-Asp-Ala.}, language = {en} } @article{DonatiWeberKeller2021, author = {Donati, Luca and Weber, Marcus and Keller, Bettina G.}, title = {Markov models from the square root approximation of the Fokker-Planck equation: calculating the grid-dependent flux}, volume = {33}, journal = {Journal of Physics: Condensed Matter}, number = {11}, doi = {10.1088/1361-648X/abd5f7}, pages = {115902}, year = {2021}, abstract = {Molecular dynamics (MD) are extremely complex, yet understanding the slow components of their dynamics is essential to understanding their macroscopic properties. To achieve this, one models the MD as a stochastic process and analyses the dominant eigenfunctions of the associated Fokker-Planck operator, or of closely related transfer operators. So far, the calculation of the discretized operators requires extensive MD simulations. The square-root approximation of the Fokker-Planck equation is a method to calculate transition rates as a ratio of the Boltzmann densities of neighboring grid cells times a flux, and can in principle be calculated without a simulation. In a previous work we still used MD simulations to determine the flux. Here, we propose several methods to calculate the exact or approximate flux for various grid types, and thus estimate the rate matrix without a simulation. Using model potentials we test computational efficiency of the methods, and the accuracy with which they reproduce the dominant eigenfunctions and eigenvalues. For these model potentials, rate matrices with up to \$\mathcal{O}\left(1{0}^{6}\right)\$ states can be obtained within seconds on a single high-performance compute server if regular grids are used.}, language = {en} } @article{YousefianFrankWeberetal.2025, author = {Yousefian, Maryam and Frank, Anne-Simone and Weber, Marcus and R{\"o}blitz, Susanna}, title = {Efficient construction of Markov state models for stochastic gene regulatory networks by domain decomposition}, volume = {26}, journal = {BMC Bioinformatics}, number = {147}, doi = {10.1186/s12859-025-06174-5}, year = {2025}, abstract = {The dynamics of many gene regulatory networks (GRNs) is characterized by the occurrence of metastable phenotypes and stochastic phenotype switches. The chemical master equation (CME) is the most accurate description to model such stochastic dynamics, whereby the long-time dynamics of the system is encoded in the spectral properties of the CME operator. Markov State Models (MSMs) provide a general framework for analyzing and visualizing stochastic multistability and state transitions based on these spectral properties. Until now, however, this approach is either limited to low-dimensional systems or requires the use of high-performance computing facilities, thus limiting its usability.}, language = {en} }