@misc{WeberWalterKubeetal.2006, author = {Weber, Marcus and Walter, Lionel and Kube, Susanna and Deuflhard, Peter}, title = {Stable computation of probability densities for metastable dynamical systems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-9331}, number = {06-39}, year = {2006}, abstract = {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.}, language = {en} } @misc{Weber2003, author = {Weber, Marcus}, title = {Clustering by using a simplex structure}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7782}, number = {04-03}, year = {2003}, abstract = {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.}, language = {en} } @misc{WeberRungsarityotinSchliep2004, author = {Weber, Marcus and Rungsarityotin, Wasinee and Schliep, Alexander}, title = {Perron Cluster Analysis and Its Connection to Graph Partitioning for Noisy Data}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-8140}, number = {04-39}, year = {2004}, abstract = {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.}, language = {en} } @misc{WeberKubeRiemeretal.2006, author = {Weber, Marcus and Kube, Susanna and Riemer, Alexander and Bujotzek, Alexander}, title = {Efficient Sampling of the Stationary Distribution of Metastable Dynamical Systems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-9467}, number = {07-03}, year = {2006}, abstract = {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.}, language = {en} } @misc{KubeWeber2005, author = {Kube, Susanna and Weber, Marcus}, title = {Identification of Metastabilities in Monomolecular Conformation Kinetics}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-8956}, number = {06-01}, year = {2005}, abstract = {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.}, language = {en} } @misc{WalterWeber2006, author = {Walter, Lionel and Weber, Marcus}, title = {ConfJump : a fast biomolecular sampling method which drills tunnels through high mountains}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-9204}, number = {06-26}, year = {2006}, abstract = {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.}, language = {en} } @phdthesis{Weber2006, author = {Weber, Marcus}, title = {Meshless Methods in Confirmation Dynamics}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-10232}, year = {2006}, language = {en} } @misc{WeberBeckerKoeppenetal.2007, author = {Weber, Marcus and Becker, Roland and K{\"o}ppen, Robert and Durmaz, Vedat}, title = {Classical hybrid Monte-Carlo simulations of the interconversion of hexabromocyclododecane}, organization = {Zuse-Institut Berlin}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-10308}, number = {07-31}, year = {2007}, abstract = {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.}, language = {en} } @misc{TungaWeber2012, author = {Tunga, Burcu and Weber, Marcus}, title = {Free Energy Calculation Using Mayer Cluster Expansion and Fluctuation Free Integration}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-16031}, year = {2012}, abstract = {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.}, language = {en} } @misc{WeberFackeldey2013, author = {Weber, Marcus and Fackeldey, Konstantin}, title = {Computing the Minimal Rebinding Effect Included in a Given Kinetics}, issn = {1438-0064}, doi = {10.1137/13091124X}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-17796}, year = {2013}, abstract = {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.}, language = {en} } @misc{LieFackeldeyWeber2013, author = {Lie, Han Cheng and Fackeldey, Konstantin and Weber, Marcus}, title = {A square root approximation of transition rates for a Markov State Model}, issn = {1438-0064}, doi = {10.1137/120899959}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-42195}, year = {2013}, abstract = {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.}, language = {en} } @misc{NielsenFackeldeyWeber2013, author = {Nielsen, Adam and Fackeldey, Konstantin and Weber, Marcus}, title = {On a Generalized Transfer Operator}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-43162}, year = {2013}, abstract = {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.}, language = {en} } @article{FackeldeyKoltaiNeviretal.2019, author = {Fackeldey, Konstantin and Koltai, Peter and Nevir, Peter and Rust, Henning and Schild, Axel and Weber, Marcus}, title = {From metastable to coherent sets - Time-discretization schemes}, volume = {29}, journal = {Chaos: An Interdisciplinary Journal of Nonlinear Science}, doi = {10.1063/1.5058128}, pages = {012101 -- 012101}, year = {2019}, abstract = {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.}, language = {en} } @article{AbbasSchneiderBollmannetal.2019, author = {Abbas, Aennes and Schneider, Ilona and Bollmann, Anna and Funke, Jan and Oehlmann, J{\"o}rg and Prasse, Carsten and Schulte-Oehlmann, Ulrike and Seitz, Wolfram and Ternes, Thomas and Weber, Marcus and Wesely, Henning and Wagner, Martin}, title = {What you extract is what you see: Optimising the preparation of water and wastewater samples for in vitro bioassays}, volume = {152}, journal = {Water Research}, doi = {10.1016/j.watres.2018.12.049}, pages = {47 -- 60}, year = {2019}, abstract = {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.}, language = {en} } @article{ReuterFackeldeyWeber2019, author = {Reuter, Bernhard and Fackeldey, Konstantin and Weber, Marcus}, title = {Generalized Markov modeling of nonreversible molecular kinetics}, volume = {17}, journal = {The Journal of Chemical Physics}, number = {150}, doi = {10.1063/1.5064530}, pages = {174103}, year = {2019}, abstract = {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.}, language = {en} } @article{ErnstFackeldeyVolkameretal.2019, author = {Ernst, Natalia and Fackeldey, Konstantin and Volkamer, Andrea and Opatz, Oliver and Weber, Marcus}, title = {Computation of temperature-dependent dissociation rates of metastable protein-ligand complexes}, volume = {45}, journal = {Molecular Simulation}, number = {11}, doi = {10.1080/08927022.2019.1610949}, pages = {904 -- 911}, year = {2019}, abstract = {Molecular simulations are often used to analyse the stability of protein-ligand complexes. The stability can be characterised by exit rates or using the exit time approach, i.e. by computing the expected holding time of the complex before its dissociation. However determining exit rates by straightforward molecular dynamics methods can be challenging for stochastic processes in which the exit event occurs very rarely. Finding a low variance procedure for collecting rare event statistics is still an open problem. In this work we discuss a novel method for computing exit rates which uses results of Robust Perron Cluster Analysis (PCCA+). This clustering method gives the possibility to define a fuzzy set by a membership function, which provides additional information of the kind 'the process is being about to leave the set'. Thus, the derived approach is not based on the exit event occurrence and, therefore, is also applicable in case of rare events. The novel method can be used to analyse the temperature effect of protein-ligand systems through the differences in exit rates, and, thus, open up new drug design strategies and therapeutic applications.}, language = {en} } @article{VillatoroWeberZuehlkeetal.2019, author = {Villatoro, Jos{\´e} and Weber, Marcus and Z{\"u}hlke, Martin and Lehmann, Andreas and Zechiowski, Karl and Riebe, Daniel and Beitz, Toralf and L{\"o}hmannsr{\"o}ben, Hans-Gerd and Kreuzer, Oliver}, title = {Structural characterization of synthetic peptides using electronspray ion mobility spectrometry and molecular dynamics simulations}, volume = {436}, journal = {International Journal of Mass Spectrometry}, doi = {10.1016/j.ijms.2018.10.036}, pages = {108 -- 117}, year = {2019}, abstract = {Electrospray ionization-ion mobility spectrometry was employed for the determination of collision cross sections (CCS) of 25 synthetically produced peptides in the mass range between 540-3310 Da. The experimental measurement of the CCS is complemented by their calculation applying two different methods. One prediction method is the intrinsic size parameter (ISP) method developed by the Clemmer group. The second new method is based on the evaluation of molecular dynamics (MD) simulation trajectories as a whole, resulting in a single, averaged collision cross-section value for a given peptide in the gas phase. A high temperature MD simulation is run in order to scan through the whole conformational space. The lower temperature conformational distribution is obtained through thermodynamic reweighting. In the first part, various correlations, e.g. CCS vs. mass and inverse mobility vs. m/z correlations, are presented. Differences in CCS between peptides are also discussed in terms of their respective mass and m/z differences, as well as their respective structures. In the second part, measured and calculated CCS are compared. The agreement between the prediction results and the experimental values is in the same range for both calculation methods. While the calculation effort of the ISP method is much lower, the MD method comprises several tools providing deeper insights into the conformations of peptides. Advantages and limitations of both methods are discussed. Based on the separation of two pairs of linear and cyclic peptides of virtually the same mass, the influence of the structure on the cross sections is discussed. The shift in cross section differences and peak shape after transition from the linear to the cyclic peptide can be well understood by applying different MD tools, e.g. the root-mean-square deviation (RMSD) and the root mean square fluctuation (RMSF).}, language = {en} } @article{DonatiHeidaKelleretal.2018, author = {Donati, Luca and Heida, Martin and Keller, Bettina G. and Weber, Marcus}, title = {Estimation of the infinitesimal generator by square-root approximation}, volume = {30}, journal = {J. Phys.: Condens. Matter}, number = {42}, doi = {10.1088/1361-648X/aadfc8}, pages = {425201 -- 425201}, year = {2018}, abstract = {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.}, language = {en} } @article{ReidelbachWeberImhof2018, author = {Reidelbach, Marco and Weber, Marcus and Imhof, Petra}, title = {Prediction of perturbed proton transfer networks}, volume = {13}, journal = {PLoS ONE}, number = {12}, doi = {https://doi.org/10.1371/journal.pone.0207718}, pages = {e0207718 -- e0207718}, year = {2018}, abstract = {The transfer of protons through proton translocating channels is a complex process, for which direct samplings of different protonation states and side chain conformations in a transition network calculation provide an efficient, bias-free description. In principle, a new transition network calculation is required for every unsampled change in the system of interest, e.g. an unsampled protonation state change, which is associated with significant computational costs. Transition networks void of or including an unsampled change are termed unperturbed or perturbed, respectively. Here, we present a prediction method, which is based on an extensive coarse-graining of the underlying transition networks to speed up the calculations. It uses the minimum spanning tree and a corresponding sensitivity analysis of an unperturbed transition network as initial guess and refinement parameter for the determination of an unknown, perturbed transition network. Thereby, the minimum spanning tree defines a sub-network connecting all nodes without cycles and minimal edge weight sum, while the sensitivity analysis analyzes the stability of the minimum spanning tree towards individual edge weight reductions. Using the prediction method, we are able to reduce the calculation costs in a model system by up to 80\%, while important network properties are maintained in most predictions.}, language = {en} } @misc{FackeldeySikorskiWeber2018, author = {Fackeldey, Konstantin and Sikorski, Alexander and Weber, Marcus}, title = {Spectral Clustering for Non-reversible Markov Chains}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-70218}, year = {2018}, 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 also extend this spectral clustering method 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 do not need to have a positive stationary distribution. In addition to metastabilities, dominant cycles and sinks can also be identified. This novel method is called GenPCCA (i.e. Generalized PCCA), since it includes the case of non reversible processes. We also apply the method to real world eye tracking data.}, language = {en} } @article{FackeldeySikorskiWeber2018, author = {Fackeldey, Konstantin and Sikorski, Alexander and Weber, Marcus}, title = {Spectral Clustering for Non-Reversible Markov Chains}, volume = {37}, journal = {Computational and Applied Mathematics}, number = {5}, doi = {https://doi.org/10.1007/s40314-018-0697-0}, pages = {6376 -- 6391}, year = {2018}, 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 also extend this spectral clustering method 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 do not need to have a positive stationary distribution. In addition to metastabilities, dominant cycles and sinks can also be identified. This novel method is called GenPCCA (i.e., generalized PCCA), since it includes the case of non-reversible processes. We also apply the method to real-world eye-tracking data.}, language = {en} } @article{LieFackeldeyWeber2013, author = {Lie, Han Cheng and Fackeldey, Konstantin and Weber, Marcus}, title = {A Square Root Approximation of Transition Rates for a Markov State Model}, volume = {34}, journal = {SIAM. J. Matrix Anal. Appl.}, number = {2}, doi = {10.1137/120899959}, pages = {738 -- 756}, year = {2013}, language = {en} } @article{ScharkoiFackeldeyMerkulowetal.2013, author = {Scharkoi, Olga and Fackeldey, Konstantin and Merkulow, Igor and Andrae, Karsten and Weber, Marcus and Nehls, Irene}, title = {Conformational Analysis of Alternariol on the Quantum Level}, volume = {19}, journal = {J. Mol. Model.}, number = {6}, doi = {10.1007/s00894-013-1803-2}, pages = {2567 -- 2572}, year = {2013}, language = {en} } @misc{BujotzekSchuettNielsenetal.2013, author = {Bujotzek, Alexander and Sch{\"u}tt, Ole and Nielsen, Adam and Fackeldey, Konstantin and Weber, Marcus}, title = {Efficient Conformational Analysis by Partition-of-Unity Coupling}, journal = {Math Chem}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-42570}, year = {2013}, abstract = {Obtaining a sufficient sampling of conformational space is a common problem in molecular simulation. We present the implementation of an umbrella-like adaptive sampling approach based on function-based meshless discretization of conformational space that is compatible with state of the art molecular dynamics code and that integrates an eigenvector-based clustering approach for conformational analysis and the computation of inter-conformational transition rates. The approach is applied to three example systems, namely n-pentane, alanine dipeptide, and a small synthetic host-guest system, the latter two including explicitly modeled solvent.}, language = {en} } @misc{QuerDonatiKelleretal.2017, author = {Quer, Jannes and Donati, Luca and Keller, Bettina and Weber, Marcus}, title = {An automatic adaptive importance sampling algorithm for molecular dynamics in reaction coordinates}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-62075}, year = {2017}, abstract = {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.}, language = {en} } @article{SpahnDelVecchioLabuzetal.2017, author = {Spahn, Viola and Del Vecchio, Giovanna and Labuz, Dominika and Rodriguez-Gaztelumendi, Antonio and Massaly, N. and Temp, Julia and Durmaz, Vedat and Sabri, P. and Reidelbach, Marco and Machelska, Halina and Weber, Marcus and Stein, Christoph}, title = {A nontoxic pain killer designed by modeling of pathological receptor conformations}, volume = {355}, journal = {Science}, number = {6328}, doi = {10.1126/science.aai8636}, pages = {966 -- 969}, year = {2017}, language = {en} } @misc{WeberFackeldeySchuette2017, author = {Weber, Marcus and Fackeldey, Konstantin and Sch{\"u}tte, Christof}, title = {Set-free Markov State Building}, issn = {1438-0064}, doi = {10.1063/1.4978501}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-62167}, year = {2017}, abstract = {Molecular dynamics (MD) simulations face challenging problems since the timescales of interest often are much longer than what is possible to simulate and even if sufficiently long simulation are possible the complex nature of the resulting simulation data makes interpretation difficult. Markov State Models (MSMs) help to overcome these problems by making experimentally relevant timescales accessible via coarse grained representations that also allows for convenient interpretation. However, standard set-based MSMs exhibit some caveats limiting their approximation quality and statistical significance. One of the main caveats results from the fact that typical MD trajectories repeatedly re-cross the boundary between the sets used to build the MSM which causes statistical bias in estimating the transition probabilities between these sets. In this article, we present a set-free approach to MSM building utilizing smooth overlapping ansatz functions instead of sets and an adaptive refinement approach. This kind of meshless discretization helps to overcome the recrossing problem and yields an adaptive refinement procedure that allows to improve the quality of the model while exploring state space and inserting new ansatz functions into the MSM.}, language = {en} } @misc{Weber2017, author = {Weber, Marcus}, title = {Eigenvalues of non-reversible Markov chains - A case study}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-62191}, year = {2017}, abstract = {Finite reversible Markov chains are characterized by a transition matrix P that has real eigenvalues and pi-orthogonal eigenvectors, where pi is the stationary distribution of P. This means, that a transition matrix with complex eigenvalues corresponds to a non-reversible Markov chain. This observation leads to the question, whether the imaginary part of that eigendecomposition corresponds to or indicates the "pattern" of the nonreversibility. This article shows that the direct relation between imaginary parts of eigendecompositions and the non-reversibility of a transition matrix is not given. It is proposed to apply the Schur decomposition of P instead of the eigendecomposition in order to characterize its nonreversibility.}, language = {en} } @article{HartmannSchuetteWeberetal.2017, author = {Hartmann, Carsten and Sch{\"u}tte, Christof and Weber, Marcus and Zhang, Wei}, title = {Importance sampling in path space for diffusion processes with slow-fast variables}, journal = {Probability Theory and Related Fields}, doi = {10.1007/s00440-017-0755-3}, pages = {1 -- 52}, year = {2017}, abstract = {Importance sampling is a widely used technique to reduce the variance of a Monte Carlo estimator by an appropriate change of measure. In this work, we study importance sampling in the framework of diffusion process and consider the change of measure which is realized by adding a control force to the original dynamics. For certain exponential type expectation, the corresponding control force of the optimal change of measure leads to a zero-variance estimator and is related to the solution of a Hamilton-Jacobi-Bellmann equation. We focus on certain diffusions with both slow and fast variables, and the main result is that we obtain an upper bound of the relative error for the importance sampling estimators with control obtained from the limiting dynamics. We demonstrate our approximation strategy with an illustrative numerical example.}, language = {en} } @article{VillatoroZuehlkeRiebeetal.2016, author = {Villatoro, Jos{\´e} and Z{\"u}hlke, Martin and Riebe, Daniel and Beitz, Toralf and Weber, Marcus and Riedel, Jens and L{\"o}hmannsr{\"o}ben, Hans-Gerd}, title = {IR-MALDI ion mobility spectrometry: physical source characterization and application as HPLC detector}, journal = {International Journal for Ion Mobility Spectrometry}, number = {19/4}, doi = {10.1007/s12127-016-0208-1}, pages = {197 -- 297}, year = {2016}, language = {en} } @article{GuerlerMollWeberetal.2008, author = {Guerler, A. and Moll, Sebastian and Weber, Marcus and Meyer, Holger and Cordes, Frank}, title = {Selection and flexible optimization of binding modes from conformation ensembles}, volume = {92}, journal = {Biosystems}, number = {1}, doi = {DOI: 10.1016/j.biosystems.2007.11.004}, pages = {42 -- 48}, year = {2008}, language = {en} } @article{HaackRoeblitzScharkoietal.2010, author = {Haack, Fiete and R{\"o}blitz, Susanna and Scharkoi, Olga and Schmidt, Burkhard and Weber, Marcus}, title = {Adaptive Spectral Clustering for Conformation Analysis}, volume = {1281}, journal = {AIP Conference Proceedings}, number = {1}, publisher = {AIP}, doi = {10.1063/1.3498116}, pages = {1585 -- 1588}, year = {2010}, language = {en} } @incollection{DurmazFackeldeyWeber2011, author = {Durmaz, Vedat and Fackeldey, Konstantin and Weber, Marcus}, title = {A rapidly Mixing Monte Carlo Method for the Simulation of Slow Molecular Processes}, booktitle = {Applications of Monte Carlo Methods in Biology, Medicine and Other Fields of Science}, editor = {Mode, Ch.}, publisher = {InTech}, year = {2011}, language = {en} } @article{GuerlerMollWeberetal.2007, author = {G{\"u}rler, A. and Moll, Sebastian and Weber, Marcus and Meyer, Holger and Cordes, Frank}, title = {Selection and flexible optimization of binding modes from conformation ensembles}, journal = {Biosystems}, year = {2007}, language = {en} } @article{KubeWeber2007, author = {Kube, Susanna and Weber, Marcus}, title = {A Coarse Graining Method for the Identification of Transition rates between Molecular Conformations}, volume = {126}, journal = {Journal of Chemical Physics}, number = {2}, doi = {10.1063/1.2404953}, year = {2007}, language = {en} } @inproceedings{WeberKube2005, author = {Weber, Marcus and Kube, Susanna}, title = {Robust Perron Cluster Analysis for Various Applications in Computational Life Science}, booktitle = {Computational Life Sciences}, pages = {57 -- 66}, year = {2005}, language = {en} } @article{MetznerWeberSchuette2010, author = {Metzner, Ph. and Weber, Marcus and Sch{\"u}tte, Christof}, title = {Observation uncertainty in reversible Markov chains}, volume = {82}, journal = {Phys. Rev. E}, number = {3}, publisher = {American Physical Society}, doi = {10.1103/PhysRevE.82.031114}, pages = {031114}, year = {2010}, language = {en} } @article{KlimmBujotzekWeber2011, author = {Klimm, Martina and Bujotzek, Alexander and Weber, Marcus}, title = {Direct Reweighting Strategies in Conformation Dynamics}, volume = {65(2)}, journal = {MATCH Commun. Math. Comp. Chem.}, pages = {333 -- 346}, year = {2011}, language = {en} } @inproceedings{KubeWeber2008, author = {Kube, Susanna and Weber, Marcus}, title = {Computation of equilibrium densities in metastable dynamical systems by domain decomposition}, volume = {1048}, booktitle = {Numerical Analysis and Applied Mathematics, International Conference on Numerical Analysis and Applied Mathematics 2008}, publisher = {AIP Conference Proceedings}, pages = {339 -- 342}, year = {2008}, language = {en} } @article{ScheibeBujotzekDerneddeetal.2011, author = {Scheibe, Ch. and Bujotzek, Alexander and Dernedde, Jens and Weber, Marcus and Seitz, O.}, title = {DNA-programmed spatial screening of carbohydrate-lectin interactions}, volume = {2}, journal = {Chem. Sci.}, pages = {770 -- 775}, year = {2011}, language = {en} } @article{SiegelAndraeProskeetal.2010, author = {Siegel, D. and Andrae, Karsten and Proske, Matthias and Kochan, C. and Koch, Matthias and Weber, Marcus and Nehls, Irene}, title = {Dynamic covalent hydrazine chemistry as a specific extraction and cleanup technique for the quantification of the Fusarium mycotoxin zearalenone in edible oils}, volume = {1217(15)}, journal = {Journal of Chromatography A}, pages = {2206 -- 15}, year = {2010}, language = {en} } @phdthesis{Weber2006, author = {Weber, Marcus}, title = {Meshless Methods in Conformation Dynamics}, year = {2006}, language = {en} } @incollection{WeberRungsarityotinSchliep2006, author = {Weber, Marcus and Rungsarityotin, Wasinee and Schliep, Alexander}, title = {An Indicator for the Number of Clusters}, booktitle = {From Data and Information Analysis to Knowledge Engineering}, editor = {Spiliopoulou, Myra and Kruse, Rudolf and Borgelt, Christian and N{\"u}rnberger, Andreas and Gaul, Wolfgang}, publisher = {Springer Berlin Heidelberg}, pages = {103 -- 110}, year = {2006}, language = {en} } @article{WeberAndrae2010, author = {Weber, Marcus and Andrae, Karsten}, title = {A simple method for the estimation of entropy differences}, volume = {63(2)}, journal = {MATCH Commun. Math. Comp. Chem. 2010}, pages = {319 -- 332}, year = {2010}, language = {en} } @article{WeberBeckerKoeppenetal.2008, author = {Weber, Marcus and Becker, Roland and K{\"o}ppen, Robert and Durmaz, Vedat}, title = {Classical hybrid Monte-Carlo simulations of the interconversion of hexabromocyclododecane}, volume = {34}, journal = {Journal of Molecular Simulation}, number = {7}, pages = {727 -- 736}, year = {2008}, language = {en} } @article{WeberBujotzekAndraeetal.2011, author = {Weber, Marcus and Bujotzek, Alexander and Andrae, Karsten and Weinhart, M. and Haag, Rainer}, title = {Computational entropy estimation of linear polyether modified surfaces and correlation with protein resistant properties of such surfaces}, journal = {J. Mol. Sim.}, year = {2011}, language = {en} } @article{FoersterBrauerFuersteetal.2007, author = {F{\"o}rster, C. and Brauer, Arnd B. E. and F{\"u}rste, J. and Betzel, C. and Weber, Marcus and Cordes, Frank and Erdmann, V.}, title = {Visualization of the tRNA(Ser) acceptor step binding site in the seryl-tRNA synthetase}, volume = {362}, journal = {BBRC}, number = {2}, pages = {415 -- 418}, year = {2007}, language = {en} } @inproceedings{KubeWeber2008, author = {Kube, Susanna and Weber, Marcus}, title = {Preserving the Markov Property of Reduced Reversible Markov Chains}, volume = {1048}, booktitle = {Numerical Analysis and Applied Mathematics, International Conference on Numerical Analysis and Applied Mathematics 2008}, pages = {593 -- 596}, year = {2008}, language = {en} } @article{WeberKubeWalteretal.2007, author = {Weber, Marcus and Kube, Susanna and Walter, Lionel and Deuflhard, Peter}, title = {Stable Computation of Probability Densities of Metastable Dynamical Systems}, volume = {6}, journal = {SIAM J. Multiscale Model. Simul.}, number = {2}, pages = {396 -- 416}, year = {2007}, language = {en} } @article{KoeppenRiedelProskeetal.2012, author = {K{\"o}ppen, Robert and Riedel, Juliane and Proske, Matthias and Drzymala, Sarah and Rasenko, Tatjana and Durmaz, Vedat and Weber, Marcus and Koch, Matthias}, title = {Photochemical trans-/cis-isomerization and quantification of zearalenone in edible oils}, volume = {60}, journal = {J. Agric. Food Chem.}, doi = {10.1021/jf3037775}, pages = {11733 -- 11740}, year = {2012}, language = {en} }