TY - CHAP A1 - Yousefian, Maryam A1 - Donati, Luca A1 - Sikorski, Alexander A1 - Weber, Marcus A1 - Röblitz, Susanna T1 - Exploring Metastable Dynamics of Gene Regulatory Networks with ISOKANN T2 - Computational Methods in Systems Biology. CMSB 2025 N2 - Stochastic dynamical systems like gene regulatory networks (GRNs) often exhibit behavior characterized by metastable sets (representing cellular phenotypes), in which trajectories remain for long times, whereas switches between these sets in the phase space are rare events. One way to capture these rare events is to infer the system’s long-term behavior from the spectral characteristics (eigenvalues and eigenvectors) of its Koopman operator. For GRNs, the Koopman operator is based on the chemical master equation (CME), which provides a precise mathematical modeling framework for stochastic GRNs. Since the CME is typically analytically intractable, methods based on discretizing the CME operator have been developed. However, determining the number and location of metastable sets in the phase space as well as the transition rates between them remains computationally challenging, especially for large GRNs with many genes and interactions. A promising alternative method, called ISOKANN (invariant subspaces of Koopman operators with artificial neural networks) has been developed in the context of molecular dynamics. ISOKANN uses a combination of the power iteration and neural networks to learn the basis functions of an invariant subspace of the Koopman operator. In this paper, we extend the application of ISOKANN to the CME operator and apply it to two small GRNs: a genetic toggle switch model and a model for macrophage polarization. Our work opens a new field of application for the ISOKANN algorithm and demonstrates the potential of this algorithm for studying large GRNs. Y1 - 2026 U6 - https://doi.org/10.1007/978-3-032-01436-8_8 VL - 15959 SP - 126 EP - 149 ER - TY - JOUR A1 - Donati, Luca A1 - Chewle, Surahit A1 - St. Pierre, Dominik A1 - Natarajan, Vijay A1 - Weber, Marcus T1 - Topological analysis reveals multiple pathways in molecular dynamics JF - Journal of Chemical Theory and Computation N2 - 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. Y1 - 2025 U6 - https://doi.org/10.1021/acs.jctc.5c00819 VL - 21 IS - 20 SP - 10385 EP - 10397 CY - J. Chem. Theory Comput. ER - TY - JOUR A1 - Sikorski, Alexander A1 - Niknejad, Amir A1 - Weber, Marcus A1 - Donati, Luca T1 - Tensor-SqRA: Modeling the transition rates of interacting molecular systems in terms of potential energies JF - Journal of Chemical Physics N2 - 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. Y1 - 2024 U6 - https://doi.org/10.1063/5.0187792 VL - 160 SP - 104112 ER - TY - JOUR A1 - Donati, Luca A1 - Schütte, Christof A1 - Weber, Marcus T1 - The Kramers turnover in terms of a macro-state projection on phase space JF - Molecular Physics N2 - 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. Y1 - 2024 U6 - https://doi.org/10.1080/00268976.2024.2356748 IS - Ciccotti Special Issue (by invitation only) SP - e2356748 PB - Taylor & Francis ER - TY - JOUR A1 - Donati, Luca A1 - Weber, Marcus T1 - Efficient Estimation of Transition Rates as Functions of pH JF - Proceedings in Applied Mathematics & Mechanics N2 - 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. Y1 - 2023 U6 - https://doi.org/10.1002/pamm.202300264 VL - 23 ER - TY - JOUR A1 - Donati, Luca A1 - Weber, Marcus T1 - Assessing transition rates as functions of environmental variables JF - The Journal of Chemical Physics N2 - We present a method to estimate the transition rates of molecular systems under different environmental conditions which cause the formation or the breaking of bonds and require the sampling of the Grand Canonical Ensemble. For this purpose, we model the molecular system in terms of probable "scenarios", governed by different potential energy functions, which are separately sampled by classical MD simulations. Reweighting the canonical distribution of each scenario according to specific environmental variables, we estimate the grand canonical distribution, then we use the Square Root Approximation (SqRA) method to discretize the Fokker-Planck operator into a rate matrix and the robust Perron Cluster Cluster Analysis (PCCA+) method to coarse-grain the kinetic model. This permits to efficiently estimate the transition rates of conformational states as functions of environmental variables, for example, the local pH at a cell membrane. In this work we formalize the theoretical framework of the procedure and we present a numerical experiment comparing the results with those provided by a constant-pH method based on non-equilibrium Molecular Dynamics Monte Carlo simulations. The method is relevant for the development of new drug design strategies which take into account how the cellular environment influences biochemical processes. Y1 - 2022 U6 - https://doi.org/10.1063/5.0109555 VL - 157 IS - 22 SP - 224103-1 EP - 224103-14 PB - AIP Publishing ER - TY - JOUR A1 - Donati, Luca A1 - Weber, Marcus A1 - Keller, Bettina G. T1 - A review of Girsanov Reweighting and of Square Root Approximation for building molecular Markov State Models JF - Journal of Mathematical Physics N2 - Dynamical reweighting methods permit to estimate kinetic observables of a stochastic process governed by a target potential U(x) from trajectories that have been generated at a different potential V(x). In this article, we present Girsanov reweighting and Square Root Approximation (SqRA): the first method reweights path probabilities exploiting the Girsanov theorem and can be applied to Markov State Models (MSMs) to reweight transition probabilities; the second method was originally developed to discretize the Fokker-Planck operator into a transition rate matrix, but here we implement it into a reweighting scheme for transition rates. We begin by reviewing the theoretical background of the methods, then present two applications relevant to Molecular Dynamics (MD), highlighting their strengths and weaknesses. Y1 - 2022 U6 - https://doi.org/10.1063/5.0127227 VL - 63 IS - 12 SP - 123306-1 EP - 123306-21 PB - AIP Publishing ER - TY - JOUR A1 - Donati, Luca A1 - Weber, Marcus A1 - Keller, Bettina G. T1 - Markov models from the square root approximation of the Fokker–Planck equation: calculating the grid-dependent flux JF - Journal of Physics: Condensed Matter N2 - 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. Y1 - 2021 U6 - https://doi.org/10.1088/1361-648X/abd5f7 VL - 33 IS - 11 SP - 115902 ER - TY - JOUR A1 - Donati, Luca A1 - Heida, Martin A1 - Keller, Bettina G. A1 - Weber, Marcus T1 - Estimation of the infinitesimal generator by square-root approximation JF - J. Phys.: Condens. Matter N2 - In recent years, for the analysis of molecular processes, the estimation of time-scales and transition rates has become fundamental. Estimating the transition rates between molecular conformations is—from a mathematical point of view—an invariant subspace projection problem. We present a method to project the infinitesimal generator acting on function space to a low-dimensional rate matrix. This projection can be performed in two steps. First, we discretize the conformational space in a Voronoi tessellation, then the transition rates between adjacent cells is approximated by the geometric average of the Boltzmann weights of the Voronoi cells. This method demonstrates that there is a direct relation between the potential energy surface of molecular structures and the transition rates of conformational changes. We will show also that this approximation is correct and converges to the generator of the Smoluchowski equation in the limit of infinitely small Voronoi cells. We present results for a two dimensional diffusion process and alanine dipeptide as a high-dimensional system. Y1 - 2018 U6 - https://doi.org/10.1088/1361-648X/aadfc8 VL - 30 IS - 42 SP - 425201 EP - 425201 ER - TY - JOUR A1 - Quer, Jannes A1 - Donati, Luca A1 - Keller, Bettina A1 - Weber, Marcus T1 - An automatic adaptive importance sampling algorithm for molecular dynamics in reaction coordinates JF - SIAM Journal on Scientific Computing N2 - In this article we propose an adaptive importance sampling scheme for dynamical quantities of high dimensional complex systems which are metastable. The main idea of this article is to combine a method coming from Molecular Dynamics Simulation, Metadynamics, with a theorem from stochastic analysis, Girsanov's theorem. The proposed algorithm has two advantages compared to a standard estimator of dynamic quantities: firstly, it is possible to produce estimators with a lower variance and, secondly, we can speed up the sampling. One of the main problems for building importance sampling schemes for metastable systems is to find the metastable region in order to manipulate the potential accordingly. Our method circumvents this problem by using an assimilated version of the Metadynamics algorithm and thus creates a non-equilibrium dynamics which is used to sample the equilibrium quantities. Y1 - 2018 U6 - https://doi.org/10.1137/17m1124772 VL - 40 IS - 2 SP - A653 EP - A670 ER -