Refine
Year of publication
Document Type
- Article (114)
- ZIB-Report (40)
- In Proceedings (6)
- Book chapter (4)
- In Collection (3)
- Other (3)
- Research data (3)
- Doctoral Thesis (2)
- Poster (2)
- Book (1)
Is part of the Bibliography
- no (179)
Keywords
- Markov State Models (5)
- metastability (5)
- cluster analysis (4)
- molecular dynamics (4)
- rare events (4)
- Conformation Dynamics (3)
- MSM (3)
- Perron cluster analysis (3)
- Schur decomposition (3)
- dynamical systems (3)
Institute
- Numerical Mathematics (136)
- Computational Molecular Design (117)
- Modeling and Simulation of Complex Processes (41)
- ZIB Allgemein (15)
- Computational Systems Biology (13)
- Mathematics for Life and Materials Science (7)
- Mathematical Algorithmic Intelligence (3)
- Parallel and Distributed Computing (3)
- Mathematical Optimization (2)
- Mathematical Optimization Methods (2)
Docking is a fundamental problem in computational biology and drug discovery that seeks to predict a ligand’s binding mode and affinity to a target protein. However, the large search space size and the complexity of the underlying physical interactions make docking a challenging task. Here, we review a docking method, based on the ant colony optimization algorithm, that ranks a set of candidate ligands by solving a minimization problem for each ligand individually. In addition, we propose an augmented version that takes into account all energy functions collectively, allowing only one minimization problem to be solved. The results show that our modification outperforms in accuracy and efficiency.
Estimating the rate of rare conformational changes in molecular systems is one of the goals of molecular dynamics simulations. In the past few decades, a lot of progress has been done in data-based approaches toward this problem. In contrast, model-based methods, such as the Square Root Approximation (SqRA), directly derive these quantities from the potential energy functions. In this article, we demonstrate how the SqRA formalism naturally blends with the tensor structure obtained by coupling multiple systems, resulting in the tensor-based Square Root Approximation (tSqRA). It enables efficient treatment of high-dimensional systems using the SqRA and provides an algebraic expression of the impact of coupling energies between molecular subsystems. Based on the tSqRA, we also develop the projected rate estimation, a hybrid data-model-based algorithm that efficiently estimates the slowest rates for coupled systems. In addition, we investigate the possibility of integrating low-rank approximations within this framework to maximize the potential of the tSqRA.
The dominant eigenfunctions of the Koopman operator characterize the metastabilities and slow-timescale dynamics of stochastic diffusion processes. In the context of molecular dynamics and Markov state modeling, they allow for a description of the location and frequencies of rare transitions, which are hard to obtain by direct simulation alone. In this article, we reformulate the eigenproblem in terms of the ISOKANN framework, an iterative algorithm that learns the eigenfunctions by alternating between short burst simulations and a mixture of machine learning and classical numerics, which naturally leads to a proof of convergence. We furthermore show how the intermediate iterates can be used to reduce the sampling variance by importance sampling and optimal control (enhanced sampling), as well as to select locations for further training (adaptive sampling). We demonstrate the usage of our proposed method in experiments, increasing the approximation accuracy by several orders of magnitude.
We have investigated how Langevin dynamics is affected by the friction
coefficient using the novel algorithm ISOKANN, which combines the transfer
operator approach with modern machine learning techniques. ISOKANN describes
the dynamics in terms of an invariant subspace projection of the Koopman
operator defined in the entire state space, avoiding approximations due to
dimensionality reduction and discretization. Our results are consistent with
the Kramers turnover and show that in the low and moderate friction regimes,
metastable macro-states and transition rates are defined in phase space, not
only in position space.
Markov processes serve as foundational models in many scientific disciplines,
such as molecular dynamics, and their simulation forms a common basis for
analysis. While simulations produce useful trajectories, obtaining macroscopic
information directly from microstate data presents significant challenges. This
paper addresses this gap by introducing the concept of membership functions
being the macrostates themselves. We derive equations for the holding times of
these macrostates and demonstrate their consistency with the classical definition.
Furthermore, we discuss the application of the ISOKANN method for learning
these quantities from simulation data. In addition, we present a novel method
for extracting transition paths based on the ISOKANN results and demonstrate
its efficacy by applying it to simulations of the 𝜇-opioid receptor. With this
approach we provide a new perspective on analyzing the macroscopic behaviour
of Markov systems.
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.
Der DWA-Themenband beschreibt ein Konzept zur weitergehenden Abwasserbehandlung für die Bewertung von Aufbereitungsverfahren, sowohl in einer Pilotphase zur Auswahl von Verfahrensoptionen als auch für die Bewertung großtechnischer Anlagen. Emissionsseitig basiert das Konzept auf bereits regulatorisch definierten Parametern wie anorganischen Stickstoff-Verbindungen oder Phosphat sowie auf neuen noch nicht in der Abwasserverordnung regulierten Parametern. Die immissionsseitige Betrachtung erfolgt auf Basis der rechtlich durch die Europäische Wasserrahmenrichtlinie und andere Anforderungen bindenden Instrumente. Hierfür werden spezifische Vorgehensweisen vorgeschlagen. Anhand zweier ausgewählter Praxisbeispiele wird deutlich, dass es zur Bewertung der Verfahrensoptionen an einem Standort dienlich ist, ausgewählte Reduktionen bzw. Entfernungen von Stoffen, Organismen und Effekten zu bestimmen.
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.