TY - GEN A1 - Weber, Marcus T1 - A Subspace Approach to Molecular Markov State Models via an Infinitesimal Generator N2 - Supercomputers can simulate complex molecular systems. However, there is a very large gap between the fastest oscillations of covalent bonds of a molecule and the time-scale of the dominant processes. In order to extract the dominant time-scales and to identify the dominant processes, a clustering of information is needed. This thesis shows that only the subspace-based Robust Perron Cluster Analysis (PCCA+) can solve this problem correctly by the construction of a Markov State Model. PCCA+ allows for time-extrapolation in molecular kinetics. This thesis shows the difference between molecular dynamics and molecular kinetics. Only in the molecular kinetics framework a definition of transition rates is possible. In this context, the existence of an infinitesimal generator of the dynamical processes is discussed. If the existence is assumed, the Theorem of Gauß can be applied in order to compute transition rates efficiently. Molecular dynamics, however, is not able to provide a suitable statistical basis for the determination of the transition pattern. T3 - ZIB-Report - 09-27 KW - Robuste Perron Cluster Analyse KW - Molekülkinetik KW - Übergangsraten KW - Robust Perron cluster analysis KW - molecular kinetics KW - transition rates Y1 - 2009 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-11432 SN - 1438-0064 ER - TY - GEN A1 - Weber, Marcus T1 - An efficient analysis of rare events in canonical ensemble dynamics N2 - For an analysis of a molecular system from a computational statistical thermodynamics point of view, extensive molecular dynamics simulations are very inefficient. During this procedure, at lot of redundant data is generated. Whereas the algorithms spend most of the computing time for a sampling of configurations within the basins of the potential energy landscape of the molecular system, the important information about the long-time behaviour of the molecules is given by transition regions and barriers between the basins, which are sampled rarely only. Thinking of molecular dynamics trajectories, researchers try to figure out which kind of dynamical model is suitable for an efficient simulation. This article suggests to change the point of view from extensive simulation of molecular dynamics trajectories to more efficient sampling strategies of the conformation dynamics approach. T3 - ZIB-Report - 08-36 KW - Moleküldynamik KW - Kanonische Gesamtheit KW - Metastabilität KW - molecular dynamics KW - canonical ensemble KW - metastability Y1 - 2008 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-10871 SN - 1438-0064 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 - THES A1 - Weber, Marcus T1 - A Subspace Approach to Molecular Markov State Models via a New Infinitesimal Generator N2 - Supercomputers can simulate complex molecular systems. However, there is a very large gap between the fastest oscillations of covalent bonds of a molecule and the time-scale of the dominant processes. In order to extract the dominant time-scales and to identify the dominant processes, a clustering of information is needed. This thesis shows that only the subspace-based Robust Perron Cluster Analysis (PCCA+) can solve this problem correctly by the construction of a Markov State Model. PCCA+ allows for time-extrapolation in molecular kinetics. This thesis shows the difference between molecular dynamics and molecular kinetics. Only in the molecular kinetics framework a definition of transition rates is possible. In this context, the existence of an infinitesimal generator of the dynamical processes is discussed. If the existence is assumed, the Theorem of Gauß can be applied in order to compute transition rates efficiently. Molecular dynamics, however, is not able to provide a suitable statistical basis for the determination of the transition pattern. KW - Conformation Dynamics KW - Molecular Kinetics KW - Transition Rates KW - Markov State Models Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-14025 ER - TY - GEN A1 - Weber, Marcus T1 - The funnel trap paradox N2 - In this article, an illustrative example is given for the coarse-graining of a Markov process which leads to a shift in the statistical weights of a two-states-system. The example is based on a 2D-funnel trap. The funnel trap is constructed in such a way, that the area inside and outside of the trap is identical. However, observing the flight of the insect as a Markov process, the probability for being “in the trap” is higher. This example can be transferred to several kinds of processes (like receptor-ligandbinding processes in chemistry) and describes the influence of “re-entering events”. T3 - ZIB-Report - 12-12 KW - Markov process KW - Robust Perron Cluster Analysis KW - Conformation Dynamics Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-14765 SN - 1438-0064 ER - TY - GEN A1 - Weber, Marcus A1 - Walter, Lionel A1 - Kube, Susanna A1 - Deuflhard, Peter T1 - Stable computation of probability densities for metastable dynamical systems N2 - Whenever the invariant stationary density of metastable dynamical systems decomposes into almost invariant partial densities, its computation as eigenvector of some transition probability matrix is an ill-conditioned problem. In order to avoid this computational difficulty, we suggest to apply an aggregation/disaggregation method which only addresses wellconditioned sub-problems and thus results in a stable algorithm. In contrast to existing methods, the aggregation step is done via a sampling algorithm which covers only small patches of the sampling space. Finally, the theoretical analysis is illustrated by two biomolecular examples. T3 - ZIB-Report - 06-39 KW - dynamical systems KW - metastability KW - molecular conformations KW - cluster analysis KW - sampling KW - aggregation/disaggregation KW - domain decomposition Y1 - 2006 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-9331 ER - TY - GEN A1 - Weber, Marcus T1 - Clustering by using a simplex structure N2 - In this paper we interpret clustering as a mapping of data into a simplex. If the data itself has simplicial struture this mapping becomes linear. Spectral analysis is an often used tool for clustering data. We will show that corresponding singular vectors or eigenvectors comprise simplicial structure. Therefore they lead to a cluster algorithm, which consists of a simple linear mapping. An example for this kind of algorithms is the Perron cluster analysis (PCCA). We have applied it in practice to identify metastable sets of molecular dynamical systems. In contrast to other algorithms, this kind of approach provides an a priori criterion to determine the number of clusters. In this paper we extend the ideas to more general problems like clustering of bipartite graphs. T3 - ZIB-Report - 04-03 KW - cluster algorithms KW - Perron cluster analysis KW - stochastic matrices KW - bipartite graphs Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7782 ER - TY - GEN A1 - Weber, Marcus A1 - Rungsarityotin, Wasinee A1 - Schliep, Alexander T1 - Perron Cluster Analysis and Its Connection to Graph Partitioning for Noisy Data N2 - The problem of clustering data can be formulated as a graph partitioning problem. Spectral methods for obtaining optimal solutions have reveceived a lot of attention recently. We describe Perron Cluster Cluster Analysis (PCCA) and, for the first time, establish a connection to spectral graph partitioning. We show that in our approach a clustering can be efficiently computed using a simple linear map of the eigenvector data. To deal with the prevalent problem of noisy and possibly overlapping data we introduce the min Chi indicator which helps in selecting the number of clusters and confirming the existence of a partition of the data. This gives a non-probabilistic alternative to statistical mixture-models. We close with showing favorable results on the analysis of gene expressi on data for two different cancer types. T3 - ZIB-Report - 04-39 KW - Perron cluster analysis KW - spectral graph theory KW - clustering KW - gene expression Y1 - 2004 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-8140 ER - TY - GEN A1 - Weber, Marcus A1 - Kube, Susanna A1 - Riemer, Alexander A1 - Bujotzek, Alexander T1 - Efficient Sampling of the Stationary Distribution of Metastable Dynamical Systems N2 - In this article we aim at an efficient sampling of the stationary distribution of dynamical systems in the presence of metastabilities. In the past decade many sophisticated algorithms have been inven ted in this field. We do not want to simply add a further one. We address the problem that one has applied a sampling algorithm for a dynamical system many times. This leads to different samplings which more or less represent the stationary distribution partially very well, but which are still far away from ergodicity or from the global stationary distribution. We will show how these samplings can be joined together in order to get one global sampling of the stationary distribution. T3 - ZIB-Report - 07-03 KW - dynamical systems KW - stationary distribution KW - rare events KW - metastability KW - cluster analysis Y1 - 2006 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-9467 ER - TY - GEN A1 - Kube, Susanna A1 - Weber, Marcus T1 - Identification of Metastabilities in Monomolecular Conformation Kinetics N2 - The identification of metastable conformations of molecules plays an important role in computational drug design. One main difficulty is the fact that the underlying dynamic processes take place in high dimensional spaces. Although the restriction of degrees of freedom to a few dihedral angles significantly reduces the complexity of the problem, the existing algorithms are time-consuming. They are based on the approximation of transition probabilities by an extensive sampling of states according to the Boltzmann distribution. We present a method which can identify metastable conformations without sampling the complete distribution. Our algorithm is based on local transition rates and uses only pointwise information about the potential energy surface. In order to apply the cluster algorithm PCCA+, we compute a few eigenvectors of the rate matrix by the Jacobi-Davidson method. Interpolation techniques are applied to approximate the thermodynamical weights of the clusters. The concluding example illustrates our approach for epigallocatechine, a molecule which can be described by seven dihedral angles. T3 - ZIB-Report - 06-01 KW - metastable conformations KW - potential energy surface KW - transition rates KW - eigenvectors KW - clustering Y1 - 2005 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-8956 ER - TY - GEN A1 - Walter, Lionel A1 - Weber, Marcus T1 - ConfJump : a fast biomolecular sampling method which drills tunnels through high mountains N2 - In order to compute the thermodynamic weights of the different metastable conformations of a molecule, we want to approximate the molecule's Boltzmann distribution in a reasonable time. This is an essential issue in computational drug design. The energy landscape of active biomolecules is generally very rough with a lot of high barriers and low regions. Many of the algorithms that perform such samplings (e.g. the hybrid Monte Carlo method) have difficulties with such landscapes. They are trapped in low-energy regions for a very long time and cannot overcome high barriers. Moving from one low-energy region to another is a very rare event. For these reasons, the distribution of the generated sampling points converges very slowly against the thermodynamically correct distribution of the molecule. The idea of ConfJump is to use $a~priori$ knowledge of the localization of low-energy regions to enhance the sampling with artificial jumps between these low-energy regions. The artificial jumps are combined with the hybrid Monte Carlo method. This allows the computation of some dynamical properties of the molecule. In ConfJump, the detailed balance condition is satisfied and the mathematically correct molecular distribution is sampled. T3 - ZIB-Report - 06-26 KW - Monte Carlo simulation KW - rare events KW - rough potential energy function KW - molecular dynamics Y1 - 2006 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-9204 ER - TY - THES A1 - Weber, Marcus T1 - Meshless Methods in Confirmation Dynamics Y1 - 2006 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-10232 ER - TY - GEN A1 - Weber, Marcus A1 - Becker, Roland A1 - Köppen, Robert A1 - Durmaz, Vedat T1 - Classical hybrid Monte-Carlo simulations of the interconversion of hexabromocyclododecane N2 - In this paper, we investigate the interconversion processes of the major flame retardant -- 1,2,5,6,9,10-hexabromocyclododecane (HBCD) -- by the means of statistical thermodynamics based on classical force-fields. Three ideas will be presented. First, the application of classical hybrid Monte-Carlo simulations for quantum mechanical processes will be justified. Second, the problem of insufficient convergence properties of hybrid Monte-Carlo methods for the generation of low temperature canonical ensembles will be solved by an interpolation approach. Furthermore, it will be shown how free energy differences can be used for a rate matrix computation. The results of our numerical simulations will be compared to experimental results. T3 - ZIB-Report - 07-31 KW - Markov process KW - molecular dynamics KW - rate matrix Y1 - 2007 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-10308 SN - 1438-0064 ER - TY - GEN A1 - Tunga, Burcu A1 - Weber, Marcus T1 - Free Energy Calculation Using Mayer Cluster Expansion and Fluctuation Free Integration N2 - This work aims to develop a new algorithm to calculate the free energy of water molecules by using a deterministic way. For this purpose, we assume a closed system confined to a physical volume, having water molecules in gas phase. To calculate the free energy of this sytem we utilized Mayer cluster expansion and the fluctuation free integration method. T3 - ZIB-Report - 12-35 Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-16031 SN - 1438-0064 ER - TY - GEN A1 - Weber, Marcus A1 - Fackeldey, Konstantin T1 - Computing the Minimal Rebinding Effect Included in a Given Kinetics N2 - The rebinding effect is a phenomenon which occurs when observing a ligand-receptor binding process. On the macro scale this process comprises the Markov property. This Makovian view is spoiled when switching to the atomistic scale of a binding process. We therefore suggest a model which accurately describes the rebinding effect on the atomistic scale by allowing ''intermediate'' bound states. This allows us to define an indicator for the magnitude of rebinding and to formulate an optimization problem. The results form our examples show good agreement with data form laboratory. T3 - ZIB-Report - 13-12 KW - Rebinding KW - Molecular Kinetics KW - Conformation Dynamics Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-17796 SN - 1438-0064 ER - TY - GEN A1 - Lie, Han Cheng A1 - Fackeldey, Konstantin A1 - Weber, Marcus T1 - A square root approximation of transition rates for a Markov State Model N2 - Trajectory- or mesh-based methods for analyzing the dynamical behavior of large molecules tend to be impractical due to the curse of dimensionality - their computational cost increases exponentially with the size of the molecule. We propose a method to break the curse by a novel square root approximation of transition rates, Monte Carlo quadrature and a discretization approach based on solving linear programs. With randomly sampled points on the molecular energy landscape and randomly generated discretizations of the molecular configuration space as our initial data, we construct a matrix describing the transition rates between adjacent discretization regions. This transition rate matrix yields a Markov State Model of the molecular dynamics. We use Perron cluster analysis and coarse-graining techniques in order to identify metastable sets in configuration space and approximate the transition rates between the metastable sets. Application of our method to a simple energy landscape on a two-dimensional configuration space provides proof of concept and an example for which we compare the performance of different discretizations. We show that the computational cost of our method grows only polynomially with the size of the molecule. However, finding discretizations of higher-dimensional configuration spaces in which metastable sets can be identified remains a challenge. T3 - ZIB-Report - 13-43 KW - Markov State Models KW - Markov chains KW - meshfree methods KW - metastability KW - Voronoi KW - linear programming Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-42195 SN - 1438-0064 ER - TY - GEN A1 - Nielsen, Adam A1 - Fackeldey, Konstantin A1 - Weber, Marcus T1 - On a Generalized Transfer Operator N2 - We introduce a generalized operator for arbitrary stochastic processes by using a pre-kernel, which is a generalization of the Markov kernel. For deterministic processes, such an operator is already known as the Frobenius-Perron operator, which is defined for a large class of measures. For Markov processes, there exists transfer operators being only well defined for stationary measures in $L^2$. Our novel generalized transfer operator is well defined for arbitrary stochastic processes, in particular also for deterministic ones. We can show that this operator is acting on $L^1$. For stationary measures, this operator is also an endomorphism of $L^2$ and, therefore, allows for a mathematical analysis in Hilbert spaces. T3 - ZIB-Report - 13-74 KW - Transfer Operator KW - Pre Kernel KW - Perron Frobenius Generalization Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-43162 SN - 1438-0064 ER - TY - JOUR A1 - Fackeldey, Konstantin A1 - Koltai, Peter A1 - Nevir, Peter A1 - Rust, Henning A1 - Schild, Axel A1 - Weber, Marcus T1 - From metastable to coherent sets - Time-discretization schemes JF - Chaos: An Interdisciplinary Journal of Nonlinear Science N2 - In this article, we show that these well-established spectral algorithms (like PCCA+, Perron Cluster Cluster Analysis) 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-time-discretization scheme. The article gives an overview about different time-discretization schemes and shows their applicability in two different fields of application. Y1 - 2019 U6 - https://doi.org/10.1063/1.5058128 VL - 29 SP - 012101 EP - 012101 ER - TY - JOUR A1 - Abbas, Aennes A1 - Schneider, Ilona A1 - Bollmann, Anna A1 - Funke, Jan A1 - Oehlmann, Jörg A1 - Prasse, Carsten A1 - Schulte-Oehlmann, Ulrike A1 - Seitz, Wolfram A1 - Ternes, Thomas A1 - Weber, Marcus A1 - Wesely, Henning A1 - Wagner, Martin T1 - What you extract is what you see: Optimising the preparation of water and wastewater samples for in vitro bioassays JF - Water Research N2 - The assessment of water quality is crucial for safeguarding drinking water resources and ecosystem integrity. To this end, sample preparation and extraction is critically important, especially when investigating emerging contaminants and the toxicity of water samples. As extraction methods are rarely optimised for bioassays but rather adopted from chemical analysis, this may result in a misrepresentation of the actual toxicity. In this study, surface water, groundwater, hospital and municipal wastewater were used to characterise the impacts of common sample preparation techniques (acidification, filtration and solid phase extraction (SPE)) on the outcomes of eleven in vitro bioassays. The latter covered endocrine activity (reporter gene assays for estrogen, androgen, aryl-hydrocarbon, retinoic acid, retinoid X, vitamin D, thyroid receptor), mutagenicity (Ames fluctuation test), genotoxicity (umu test) and cytotoxicity. Water samples extracted using different SPE sorbents (Oasis HLB, Supelco ENVI-Carb+, Telos C18/ENV) at acidic and neutral pH were compared for their performance in recovering biological effects. Acidification, commonly used for stabilisation, significantly altered the endocrine activity and toxicity of most (waste)water samples. Sample filtration did not affect the majority of endpoints but in certain cases affected the (anti-)estrogenic and dioxin-like activities. SPE extracts (10.4 × final concentration), including WWTP effluents, induced significant endocrine effects that were not detected in aqueous samples (0.63 × final concentration), such as estrogenic, (anti-)androgenic and dioxin-like activities. When ranking the SPE methods using multivariate Pareto optimisation an extraction with Telos C18/ENV at pH 7 was most effective in recovering toxicity. At the same time, these extracts were highly cytotoxic masking the endpoint under investigation. Compared to that, extraction at pH 2.5 enriched less cytotoxicity. In summary, our study demonstrates that sample preparation and extraction critically affect the outcome of bioassays when assessing the toxicity of water samples. Depending on the water matrix and the bioassay, these methods need to be optimised to accurately assess water quality. Y1 - 2019 U6 - https://doi.org/10.1016/j.watres.2018.12.049 VL - 152 SP - 47 EP - 60 ER - TY - JOUR A1 - Reuter, Bernhard A1 - Fackeldey, Konstantin A1 - Weber, Marcus T1 - Generalized Markov modeling of nonreversible molecular kinetics JF - The Journal of Chemical Physics N2 - Markov state models are to date the gold standard for modeling molecular kinetics since they enable the identification and analysis of metastable states and related kinetics in a very instructive manner. The state-of-the-art Markov state modeling methods and tools are very well developed for the modeling of reversible processes in closed equilibrium systems. On the contrary, they are largely not well suited to deal with nonreversible or even nonautonomous processes of nonequilibrium systems. Thus, we generalized the common Robust Perron Cluster Cluster Analysis (PCCA+) method to enable straightforward modeling of nonequilibrium systems as well. The resulting Generalized PCCA (G-PCCA) method readily handles equilibrium as well as nonequilibrium data by utilizing real Schur vectors instead of eigenvectors. This is implemented in the G-PCCA algorithm that enables the semiautomatic coarse graining of molecular kinetics. G-PCCA is not limited to the detection of metastable states but also enables the identification and modeling of cyclic processes. This is demonstrated by three typical examples of nonreversible systems. Y1 - 2019 U6 - https://doi.org/10.1063/1.5064530 VL - 17 IS - 150 SP - 174103 ER -