Refine
Year of publication
Document Type
- Article (112)
- ZIB-Report (40)
- In Proceedings (5)
- Book chapter (4)
- In Collection (3)
- Other (3)
- Research data (3)
- Doctoral Thesis (2)
- Poster (2)
- Book (1)
Is part of the Bibliography
- no (176)
Keywords
- Markov State Models (5)
- metastability (5)
- cluster analysis (4)
- molecular dynamics (4)
- rare events (4)
- Conformation Dynamics (3)
- MSM (3)
- Perron cluster analysis (3)
- Schur decomposition (3)
- dynamical systems (3)
Institute
- Numerical Mathematics (136)
- Computational Molecular Design (117)
- Modeling and Simulation of Complex Processes (39)
- ZIB Allgemein (15)
- Computational Systems Biology (13)
- Mathematics for Life and Materials Science (7)
- Parallel and Distributed Computing (3)
- Mathematical Algorithmic Intelligence (2)
- Mathematical Optimization (2)
- Mathematical Optimization Methods (2)
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.
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.
Molecular processes such as protein folding or ligand-receptor-binding can be understood by analyzing the free energy landscape. Those processes are often metastable, i.e. the molecular systems remain in basins around local minima of the free energy landscape, and in rare cases undergo gauche transitions between metastable states by passing saddle-points of this landscape. By discretizing the configuration space, this can be modeled as a discrete Markov process. One way to compute the transition rates between conformations of a molecular system is by utilizing Transition Path Theory and the concept of committor functions. A fundamental problem from the computational point of view is that many time-scales are involved, ranging from 10^(-14) sec for the fastest motion to 10^(-6) sec or more for conformation changes that cause biological effects.
The goal of our work is to provide a better understanding of such transitions in configuration space on various time-scales by analyzing characteristic scalar functions topologically and geometrically. We are developing suitable visualization and interaction techniques to support our analysis. For example, we are analyzing a transition rate indicator function by computing and visualizing its Reeb graph together with the sets of molecular states corresponding to maxima of the transition rate indicator function. A particular challenge is the high dimensionality of the domain which does not allow for a straightforward visualization of the function.
The computational topology approach to the analysis of the transition rate indicator functions for a molecular system allows to explore different time scales of the system by utilizing coarser or finer topological partitioning of the function. A specific goal is the development of tools for analyzing the hierarchy of these partitionings. This approach tackles the analysis of a complex and sparse dataset from a different angle than the well-known spectral analysis of Markov State Models.
We analyze how Langevin dynamics is affected by the friction coefficient using an invariant subspace projection of the associated Koopman operator. This provides the friction-dependent metastable macro-states of the dynamical system as well as the transition rates in the entire phase space. We used the algorithm ISOKANN for a wide range of friction coefficient values and reproduced results consistent with the Kramers turnover.