@misc{DeuflhardWeber2003, author = {Deuflhard, Peter and Weber, Marcus}, title = {Robust Perron Cluster Analysis in Conformation Dynamics}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7415}, number = {03-19}, year = {2003}, abstract = {The key to molecular conformation dynamics is the direct identification of metastable conformations, which are almost invariant sets of molecular dynamical systems. Once some reversible Markov operator has been discretized, a generalized symmetric stochastic matrix arises. This matrix can be treated by Perron cluster analysis, a rather recent method involving a Perron cluster eigenproblem. The paper presents an improved Perron cluster analysis algorithm, which is more robust than earlier suggestions. Numerical examples are included.}, language = {en} } @misc{WeberMeyer2005, author = {Weber, Marcus and Meyer, Holger}, title = {ZIBgridfree - Adaptive Conformation Analysis with qualified Support of Transition States and Thermodynamic Weights}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-8514}, number = {05-17}, year = {2005}, abstract = {This paper introduces a new algorithm of conformational analysis based on mesh-free methods as described in [M. Weber. Mehless methods in Conformation Dynamics.(2005)]. The adaptive decomposition of the conformational space by softly limiting functions avoids trapping effects and allows adaptive refinement strategies. These properties of the algorithm makes ZIBgridfree particularly suitable for the complete exploration of high-dimensional conformational space. The adaptive control of the algorithm benefits from the tight integration of molecular simulation and conformational analysis. An emphasized part of the analysis is the Robust Perron Cluster Analysis (PCCA+) based on the work of Peter Deuflhard and Marcus Weber. PCCA+ supports an almost-characteristic cluster definition with an outstanding mapping of transition states. The outcome is expressed by the metastable sets of conformations, their thermodynamic weights and flexibility.}, language = {en} } @misc{WeberFackeldey2015, author = {Weber, Marcus and Fackeldey, Konstantin}, title = {G-PCCA: Spectral Clustering for Non-reversible Markov Chains}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-55505}, year = {2015}, abstract = {Spectral clustering methods are based on solving eigenvalue problems for the identification of clusters, e.g., the identification of metastable subsets of a Markov chain. Usually, real-valued eigenvectors are mandatory for this type of algorithms. The Perron Cluster Analysis (PCCA+) is a well-known spectral clustering method of Markov chains. It is applicable for reversible Markov chains, because reversibility implies a real-valued spectrum. We extend this spectral clustering method also to non-reversible Markov chains and give some illustrative examples. The main idea is to replace the eigenvalue problem by a real-valued Schur decomposition. By this extension, non-reversible Markov chains can be analyzed. Furthermore, the chains need not have a positive stationary distribution. And additionally to metastabilities, dominant cycles and sinks can be identified, too.}, language = {en} } @misc{DjurdjevacConradWeberSchuette2015, author = {Djurdjevac Conrad, Natasa and Weber, Marcus and Sch{\"u}tte, Christof}, title = {Finding dominant structures of nonreversible Markov processes}, issn = {1438-0064}, doi = {10.1137/15M1032272}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-55739}, year = {2015}, abstract = {Finding metastable sets as dominant structures of Markov processes has been shown to be especially useful in modeling interesting slow dynamics of various real world complex processes. Furthermore, coarse graining of such processes based on their dominant structures leads to better understanding and dimension reduction of observed systems. However, in many cases, e.g. for nonreversible Markov processes, dominant structures are often not formed by metastable sets but by important cycles or mixture of both. This paper aims at understanding and identifying these different types of dominant structures for reversible as well as nonreversible ergodic Markov processes. Our algorithmic approach generalizes spectral based methods for reversible process by using Schur decomposition techniques which can tackle also nonreversible cases. We illustrate the mathematical construction of our new approach by numerical experiments.}, language = {en} } @misc{NielsenWeber2014, author = {Nielsen, Adam and Weber, Marcus}, title = {Computing the nearest reversible Markov chain}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-53292}, year = {2014}, abstract = {Reversible Markov chains are the basis of many applications. However, computing transition probabilities by a finite sampling of a Markov chain can lead to truncation errors. Even if the original Markov chain is reversible, the approximated Markov chain might be non-reversible and will lose important properties, like the real valued spectrum. In this paper, we show how to find the closest reversible Markov chain to a given transition matrix. It turns out that this matrix can be computed by solving a convex minimization problem.}, language = {en} } @misc{ZhangWangHartmannetal.2014, author = {Zhang, Wei and Wang, Han and Hartmann, Carsten and Weber, Marcus and Sch{\"u}tte, Christof}, title = {Applications of the cross-entropy method to importance sampling and optimal control of diffusions}, issn = {1438-0064}, doi = {10.1137/14096493X}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-49720}, year = {2014}, abstract = {We study the cross-entropy method for diffusions. One of the results is a versatile cross-entropy algorithm that can be used to design efficient importance sampling strategies for rare events or to solve optimal control problems. The approach is based on the minimization of a suitable cross-entropy functional, with a parametric family of exponentially tilted probability distributions. We illustrate the new algorithm with several numerical examples and discuss algorithmic issues and possible extensions of the method.}, language = {en} } @misc{WeberQuer2015, author = {Weber, Marcus and Quer, Jannes}, title = {Estimating exit rates in rare event dynamical systems via extrapolation}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-56622}, year = {2015}, abstract = {In this article we present a new idea for approximating exit rates for diffusion processes living in a craggy landscape. We are especially interested in the exit rates of a process living in a metastable regions. Due to the fact that Monte Carlo simulations perform quite poor and are very computational expensive in this setting we create several similar situations with a smoothed potential. For this we introduce a new parameter \$\lambda \in [0,1]\$ (\$\lambda = 1\$ very smoothed potential, \$\lambda=0\$ original potential) into the potential which controls the influence the smoothing. We then sample the exit rate for different parameters \$\lambda\$ the exit rate from a given region. Due to the fact that \$\lambda\$ is connected to the exit rate we can use this dependency to approximate the real exit rate. The method can be seen as something between hyperdynamics and temperature accelerated MC.}, language = {en} } @misc{SchuetteNielsenWeber2014, author = {Sch{\"u}tte, Christof and Nielsen, Adam and Weber, Marcus}, title = {Markov State Models and Molecular Alchemy}, issn = {1438-0064}, doi = {10.1080/00268976.2014.944597}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-46718}, year = {2014}, abstract = {In recent years Markov State Models (MSMs) have attracted a consid- erable amount of attention with regard to modelling conformation changes and associated function of biomolecular systems. They have been used successfully, e.g., for peptides including time-resolved spectroscopic ex- periments, protein function and protein folding , DNA and RNA, and ligand-receptor interaction in drug design and more complicated multi- valent scenarios. In this article a novel reweighting scheme is introduced that allows to construct an MSM for certain molecular system out of an MSM for a similar system. This permits studying how molecular proper- ties on long timescales differ between similar molecular systems without performing full molecular dynamics simulations for each system under con- sideration. The performance of the reweighting scheme is illustrated for simple test cases including one where the main wells of the respective en- ergy landscapes are located differently and an alchemical transformation of butane to pentane where the dimension of the state space is changed.}, language = {en} } @misc{FackeldeyKoltaiNeviretal.2017, author = {Fackeldey, Konstantin and Koltai, P{\´e}ter and N{\´e}vir, Peter and Rust, Henning and Schild, Axel and Weber, Marcus}, title = {From Metastable to Coherent Sets - time-discretization schemes}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-66074}, year = {2017}, abstract = {Given a time-dependent stochastic process with trajectories x(t) in a space \$\Omega\$, there may be sets such that the corresponding trajectories only very rarely cross the boundaries of these sets. We can analyze such a process in terms of metastability or coherence. Metastable sets M are defined in space \$M\subset\Omega\$, coherent sets \$M(t)\subset\Omega\$ are defined in space and time. Hence, if we extend the space by the time-variable t, coherent sets are metastable sets in \$\Omega\times[0,\infty]\$. This relation can be exploited, because there already exist spectral algorithms for the identification of metastable sets. In this article we show that these well-established spectral algorithms (like PCCA+) also identify coherent sets of non-autonomous dynamical systems. For the identification of coherent sets, one has to compute a discretization (a matrix T) of the transfer operator of the process using a space-timediscretization scheme. The article gives an overview about different time-discretization schemes and shows their applicability in two different fields of application.}, language = {en} } @misc{WeberGalliat2002, author = {Weber, Marcus and Galliat, Tobias}, title = {Characterization of Transition States in Conformational Dynamics using Fuzzy Sets}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6791}, number = {02-12}, year = {2002}, abstract = {Recently, a novel approach for the analysis of molecular dynamics on the basis of a transfer operator has been introduced. Therein conformations are considered to be disjoint metastable clusters within position space of a molecule. These clusters are defined by almost invariant characteristic functions that can be computed via {\em Perron Cluster} analysis. The present paper suggests to replace crisp clusters with {\em fuzzy} clusters, i.e. to replace characteristic functions with membership functions. This allows a more sufficient characterization of transiton states between different confor conformations and therefore leads to a better understanding of molecular dynamics. Fur thermore, an indicator for the uniqueness of metastable fuzzy clusters and a fast algorithm for the computation of these clusters are described. Numerical examples are included.}, language = {en} } @misc{Weber2007, author = {Weber, Marcus}, title = {Conformation-based transition state theory}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-9612}, number = {07-18}, year = {2007}, abstract = {For the treatment of equilibrated molecular systems in a heat bath we propose a transition state theory that is based on conformation dynamics. In general, a set-based discretization of a Markov operator \${\cal P}^\tau\$ does not preserve the Markov property. In this article, we propose a discretization method which is based on a Galerkin approach. This discretization method preserves the Markov property of the operator and can be interpreted as a decomposition of the state space into (fuzzy) sets. The conformation-based transition state theory presented here can be seen as a first step in conformation dynamics towards the computation of essential dynamical properties of molecular systems without time-consuming molecular dynamics simulations.}, language = {en} } @misc{KubeWeber2005, author = {Kube, Susanna and Weber, Marcus}, title = {Conformation Kinetics as a Reduced Model for Transition Pathways}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-8763}, number = {05-43}, year = {2005}, abstract = {The complexity of molecular kinetics can be reduced significantly by a restriction to metastable conformations which are almost invariant sets of molecular dynamical systems. With the Robust Perron Cl uster Analysis PCCA+, developed by Weber and Deuflhard, we have a tool available which can be used to identify these conformations from a transition probability matrix. This method can also be applied to the corresponding transition rate matrix which provides important information concerning transition pathways of single molecules. In the present paper, we explain the relationship between these tw o concepts and the extraction of conformation kinetics from transition rates. Moreover, we show how transition rates can be approximated and conclude with numerical examples.}, language = {en} } @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} } @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} } @misc{WitzigBeckenbachEifleretal.2016, author = {Witzig, Jakob and Beckenbach, Isabel and Eifler, Leon and Fackeldey, Konstantin and Gleixner, Ambros and Grever, Andreas and Weber, Marcus}, title = {Mixed-Integer Programming for Cycle Detection in Non-reversible Markov Processes}, issn = {1438-0064}, doi = {10.1137/16M1091162}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-60353}, year = {2016}, abstract = {In this paper, we present a new, optimization-based method to exhibit cyclic behavior in non-reversible stochastic processes. While our method is general, it is strongly motivated by discrete simulations of ordinary differential equations representing non-reversible biological processes, in particular molecular simulations. Here, the discrete time steps of the simulation are often very small compared to the time scale of interest, i.e., of the whole process. In this setting, the detection of a global cyclic behavior of the process becomes difficult because transitions between individual states may appear almost reversible on the small time scale of the simulation. We address this difficulty using a mixed-integer programming model that allows us to compute a cycle of clusters with maximum net flow, i.e., large forward and small backward probability. For a synthetic genetic regulatory network consisting of a ring-oscillator with three genes, we show that this approach can detect the most productive overall cycle, outperforming classical spectral analysis methods. Our method applies to general non-equilibrium steady state systems such as catalytic reactions, for which the objective value computes the effectiveness of the catalyst.}, language = {en} }