Refine
Year of publication
Document Type
- Article (113) (remove)
Has Fulltext
- no (113)
Is part of the Bibliography
- no (113)
Keywords
- MSM (2)
- Convex Optimization (1)
- GenPCCA (1)
- Girsanov (1)
- Markov State Models (1)
- Markov chain (1)
- Mixed-Integer Programming (1)
- NESS (1)
- Non-reversible (1)
- Non-reversible Markov Processes (1)
Institute
- Numerical Mathematics (91)
- Computational Molecular Design (81)
- Modeling and Simulation of Complex Processes (34)
- Computational Systems Biology (9)
- Mathematics for Life and Materials Science (6)
- Parallel and Distributed Computing (3)
- Mathematical Algorithmic Intelligence (2)
- AI in Society, Science, and Technology (1)
- Mathematical Optimization (1)
- Mathematical Optimization Methods (1)
A Coarse Graining Method for the Identification of Transition rates between Molecular Conformations
(2007)
Efficient Simulation of Ligand-Receptor Binding Processes Using the Conformation Dynamics Approach
(2009)
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 experiments, protein function and protein folding , DNA and RNA, and ligand-receptor interaction in drug design and more complicated multivalent 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 properties 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 energy landscapes are located differently and an alchemical transformation of butane to pentane where the dimension of the state space is changed.
A hands-off linear interaction energy approach to binding mode and affinity estimation of estrogens
(2013)
Upon ligand binding or during chemical reactions the state of a molecular system changes in time. Usually we consider a finite set of (macro-) states of the system (e.g., ’bound’ vs. ’unbound’), although the process itself takes place in a continuous space. In this context, the formula chi=XA connects the micro-dynamics of the molecular system to its macro-dynamics. Chi can be understood as a clustering of micro-states of a molecular system into a few macro-states. X is a basis of an invariant subspace of a transfer operator describing the micro-dynamics of the system. The formula claims that there is an unknown linear relation A between these two objects. With the aid of this formula we can understand rebinding effects, the electron flux in pericyclic reactions, and systematic changes of binding rates in kinetic ITC experiments. We can also analyze sequential spectroscopy experiments and rare event systems more easily. This article provides an explanation of the formula and an overview of some of its consequences.
Stimulation of renal collecting duct principal cells with antidiuretic hormone (arginine-vasopressin, AVP) results in inhibition of the small GTPase RhoA and the enrichment of the water channel aquaporin-2 (AQP2) in the plasma membrane. The membrane insertion facilitates water reabsorption from primary urine and fine-tuning of body water homeostasis. Rho guanine nucleotide exchange factors (GEFs) interact with RhoA, catalyze the exchange of GDP for GTP and thereby activate the GTPase. However, GEFs involved in the control of AQP2 in renal principal cells are unknown. The A-kinase anchoring protein, AKAP-Lbc, possesses GEF activity, specifically activates RhoA, and is expressed in primary renal inner medullary collecting duct principal (IMCD) cells. Through screening of 18,431 small molecules and synthesis of a focused library around one of the hits, we identified an inhibitor of the interaction of AKAP-Lbc and RhoA. This molecule, Scaff10-8, bound to RhoA, inhibited the AKAP-Lbc-mediated RhoA activation but did not interfere with RhoA activation through other GEFs or activities of other members of the Rho family of small GTPases, Rac1 and Cdc42. Scaff10-8 promoted the redistribution of AQP2 from intracellular vesicles to the periphery of IMCD cells. Thus, our data demonstrate an involvement of AKAP-Lbc-mediated RhoA activation in the control of AQP2 trafficking.
Fuzzy spectral clustering by PCCA+: application to Markov state models and data classification
(2013)
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.
In addition to the conventional Isothermal Titration Calorimetry (ITC), kinetic ITC (kinITC) not only gains thermodynamic information, but also kinetic data from a biochemical binding process. Moreover, kinITC gives insights into reactions consisting of two separate kinetic steps, such as protein folding or sequential binding processes. The ITC method alone cannot deliver kinetic parameters, especially not for multivalent bindings. This paper describes how to solve the problem using kinITC and an invariant subspace projection. The algorithm is tested for multivalent systems with different valencies.
Mussel glue‐proteins undergo structural transitions at material interfaces to optimize adhesive surface contacts. Those intriguing structure responses are mimicked by a mussel‐glue mimetic peptide (HSY*SGWSPY*RSG (Y* = l‐Dopa)) that was previously selected by phage‐display to adhere to Al2O3 after enzymatic activation. Molecular level insights into the full‐length adhesion domain at Al2O3 surfaces are provided by a divergent‐convergent analysis, combining nuclear Overhauser enhancement based 2D NOESY and saturation transfer difference NMR analysis of submotifs along with molecular dynamics simulations of the full‐length peptide. The peptide is divided into two submotifs, each containing one Dopa “anchor” (Motif‐1 and 2). The analysis proves Motif‐1 to constitute a dynamic Al2O3 binder and adopting an “M”‐structure with multiple surface contacts. Motif‐2 binds stronger by two surface contacts, forming a compact “C”‐structure. Taking these datasets as constraints enables to predict the structure and propose a binding process model of the full‐length peptide adhering to Al2O3.
The non-selective activation of central and peripheral opioid receptors is a major shortcoming of currently available opioids. Targeting peripheral opioid receptors is a promising strategy to preclude side effects. Recently, we showed that fentanyl-derived μ-opioid receptor (MOR) agonists with reduced acid dissociation constants (pKa) due to introducing single fluorine atoms produced injury-restricted antinociception in rat models of inflammatory, postoperative and neuropathic pain. Here, we report that a new double-fluorinated compound (FF6) and fentanyl show similar pKa, MOR affinity and [35S]-GTPγS binding at low and physiological pH values. In vivo, FF6 produced antinociception in injured and non-injured tissue, and induced sedation and constipation. The comparison of several fentanyl derivatives revealed a correlation between pKa values and pH-dependent MOR activation, antinociception and side effects. An opioid ligand's pKa value may be used as discriminating factor to design safer analgesics.
A new ion mobility (IM) spectrometer, enabling mobility measurements in the pressure range between 5 and 500 mbar and in the reduced field strength range E/N of 5–90 Td, was developed and characterized. Reduced mobility (K0) values were studied under low E/N (constant value) as well as high E/N (deviation from low field K0) for a series of molecular ions in nitrogen. Infrared matrix-assisted laser desorption ionization (IR-MALDI) was used in two configurations: a source working at atmospheric pressure (AP) and, for the first time, an IR-MALDI source working with a liquid (aqueous) matrix at sub-ambient/reduced pressure (RP). The influence of RP on IR-MALDI was examined and new insights into the dispersion process were gained. This enabled the optimization of the IM spectrometer for best analytical performance. While ion desolvation is less efficient at RP, the transport of ions is more efficient, leading to intensity enhancement and an increased number of oligomer ions. When deciding between AP and RP IR-MALDI, a trade-off between intensity and resolving power has to be considered. Here, the low field mobility of peptide ions was first measured and compared with reference values from ESI-IM spectrometry (at AP) as well as collision cross sections obtained from molecular dynamics simulations. The second application was the determination of the reduced mobility of various substituted ammonium ions as a function of E/N in nitrogen. The mobility is constant up to a threshold at high E/N. Beyond this threshold, mobility increases were observed. This behavior can be explained by the loss of hydrated water molecules.
The Augmented Jump Chain
(2021)
Modern methods of simulating molecular systems are based on the mathematical theory of Markov operators with a focus on autonomous equilibrated systems. However, non-autonomous physical systems or non-autonomous simulation processes are becoming more and more important. A representation of non-autonomous Markov jump processes is presented as autonomous Markov chains on space-time. Augmenting the spatial information of the embedded Markov chain by the temporal information of the associated jump times, the so-called augmented jump chain is derived. The augmented jump chain inherits the sparseness of the infinitesimal generator of the original process and therefore provides a useful tool for studying time-dependent dynamics even in high dimensions. Furthermore, possible generalizations and applications to the computation of committor functions and coherent sets in the non-autonomous setting are discussed. After deriving the theoretical foundations, the concepts with a proof-of-concept Galerkin discretization of the transfer operator of the augmented jump chain applied to simple examples are illustrated.
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.
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).
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.
Small-molecule oxoanions are often imprinted noncovalently as carboxylates into molecularly imprinted polymers (MIPs), requiring the use of an organic counterion. Popular species are either pentamethylpiperidine (PMP) as a protonatable cation or tetraalkylammonium (TXA) ions as permanent cations. The present work explores the influence of the TXA as a function of their alkyl chain length, from methyl to octyl, using UV/vis absorption, fluorescence titrations, and HPLC as well as MD simulations. Protected phenylalanines (Z-L/D-Phe) served as templates/analytes. While the influence of the counterion on the complex stability constants and anion-induced spectral changes shows a monotonous trend with increasing alkyl chain length at the prepolymerization stage, the cross-imprinting/rebinding studies showed a unique pattern that suggested the presence of adaptive cavities in the MIP matrix, related to the concept of induced fit of enzyme–substrate interaction. Larger cavities formed in the presence of larger counterions can take up pairs of Z-x-Phe and smaller TXA, eventually escaping spectroscopic detection. Correlation of the experimental data with the MD simulations revealed that counterion mobility, the relative distances between the three partners, and the hydrogen bond lifetimes are more decisive for the response features observed than actual distances between interacting atoms in a complex or the orientation of binding moieties. TBA has been found to yield the highest imprinting factor, also showing a unique dual behavior regarding the interaction with template and fluorescent monomer. Finally, interesting differences between both enantiomers have been observed in both theory and experiment, suggesting true control of enantioselectivity. The contribution concludes with suggestions for translating the findings into actual MIP development.
This work addresses the problem of determining the number of components from sequential spectroscopic data analyzed by non-negative matrix factorization without separability assumption (SepFree NMF). These data are stored in a matrix M of dimension “measured times” versus “measured wavenumbers” and can be decomposed to obtain the spectral fingerprints of the states and their evolution over time. SepFree NMF assumes a memoryless (Markovian) process to underline the dynamics and decomposes M so that M=WH, with W representing the components’ fingerprints and H their kinetics. However, the rank of this decomposition (i.e., the number of physical states in the process) has to be guessed from pre-existing knowledge on the observed process. We propose a measure for determining the number of components with the computation of the minimal memory effect resulting from the decomposition; by quantifying how much the obtained factorization is deviating from the Markovian property, we are able to score factorizations of a different number of components. In this way, we estimate the number of different entities which contribute to the observed system, and we can extract kinetic information without knowing the characteristic spectra of the single components. This manuscript provides the mathematical background as well as an analysis of computer generated and experimental sequentially measured Raman spectra.
Plasma Proteome Fingerprints Reveal Distinctiveness and Clinical Outcome of SARS-CoV-2 Infection
(2021)
We evaluated how plasma proteomic signatures in patients with suspected COVID-19 can unravel the pathophysiology, and determine kinetics and clinical outcome of the infection. We identified distinct plasma proteins linked to the presence and course of COVID-19. These plasma proteomic findings may translate to a protein fingerprint, helping to assist clinical management decisions.
Dynamical reweighting methods permit to estimate kinetic observables of a stochastic process governed by a target potential U(x) from trajectories that have been generated at a different potential V(x). In this article, we present Girsanov reweighting and Square Root Approximation (SqRA): the first method reweights path probabilities exploiting the Girsanov theorem and can be applied to Markov State Models (MSMs) to reweight transition probabilities; the second method was originally developed to discretize the Fokker-Planck operator into a transition rate matrix, but here we implement it into a reweighting scheme for transition rates. We begin by reviewing the theoretical background of the methods, then present two applications relevant to Molecular Dynamics (MD), highlighting their strengths and weaknesses.
Cyanine Dye Coupling Mediates Self-assembly of a pH Sensitive Peptide into Novel 3D Architectures
(2022)
A conjugated Cy5 dye-peptide system reveals the formation of two novel and structurally distinct supramolecular assemblies with photo-physical characteristics of H-type dimers or tetramers, respectively. The molecular ultrastructures are triggered by the complementary interplay of mutual chromophore coupling and pH induced changes in the peptide charge pattern.
We study consistency of cell-centered finite difference methods for elliptic equations with degenerate coefficients in any space dimension $d \geq 2$. This results in order of convergence estimates in the natural weighted energy norm and in the weighted discrete $L^2$-norm on admissible meshes. The cells of meshes under consideration may be very irregular in size. We particularly allow the size of certain cells to remain bounded from below even in the asymptotic limit. For uniform meshes we show that the order of convergence is at least 1 in the energy semi-norm, provided the discrete and continuous solutions exist and the continuous solution has $H^2$ regularity.
The choice of solvents influences crystalline solid formed during the crystallization of active pharmaceutical ingredients (API). The underlying effects are not always well understood because of the complexity of the systems. Theoretical models are often insufficient to describe this phenomenon. In this study, the crystallization behavior of the model drug paracetamol in different solvents was studied based on experimental and molecular dynamics data. The crystallization process was followed in situ using time-resolved Raman spectroscopy. Molecular dynamics with simulated annealing algorithm was used for an atomistic understanding of the underlying processes. The experimental and theoretical data indicate that paracetamol molecules adopt a particular geometry in a given solvent predefining the crystallization of certain polymorphs.
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.
Molecular dynamics (MD) are extremely complex, yet understanding the slow components of their dynamics is essential to understanding their macroscopic properties. To achieve this, one models the MD as a stochastic process and analyses the dominant eigenfunctions of the associated Fokker–Planck operator, or of closely related transfer operators. So far, the calculation of the discretized operators requires extensive MD simulations. The square-root approximation of the Fokker–Planck equation is a method to calculate transition rates as a ratio of the Boltzmann densities of neighboring grid cells times a flux, and can in principle be calculated without a simulation. In a previous work we still used MD simulations to determine the flux. Here, we propose several methods to calculate the exact or approximate flux for various grid types, and thus estimate the rate matrix without a simulation. Using model potentials we test computational efficiency of the methods, and the accuracy with which they reproduce the dominant eigenfunctions and eigenvalues. For these model potentials, rate matrices with up to $\mathcal{O}\left(1{0}^{6}\right)$ states can be obtained within seconds on a single high-performance compute server if regular grids are used.
We present a method to estimate the transition rates of molecular systems under different environmental conditions which cause the formation or the breaking of bonds and require the sampling of the Grand Canonical Ensemble. For this purpose, we model the molecular system in terms of probable "scenarios", governed by different potential energy functions, which are separately sampled by classical MD simulations. Reweighting the canonical distribution of each scenario according to specific environmental variables, we estimate the grand canonical distribution, then we use the Square Root Approximation (SqRA) method to discretize the Fokker-Planck operator into a rate matrix and the robust Perron Cluster Cluster Analysis (PCCA+) method to coarse-grain the kinetic model. This permits to efficiently estimate the transition rates of conformational states as functions of environmental variables, for example, the local pH at a cell membrane. In this work we formalize the theoretical framework of the procedure and we present a numerical experiment comparing the results with those provided by a constant-pH method based on non-equilibrium Molecular Dynamics Monte Carlo simulations. The method is relevant for the development of new drug design strategies which take into account how the cellular environment influences biochemical processes.
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.
An automatic adaptive importance sampling algorithm for molecular dynamics in reaction coordinates
(2018)
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.
The aim of this paper is to investigate the rebinding effect, a phenomenon describing a "short-time memory" which can occur when projecting a Markov process onto a smaller state space. For guaranteeing a correct mapping by the Markov State Model, we assume a fuzzy clustering in terms of membership functions, assigning degrees of membership to each state. The macro states are represented by the membership functions and may be overlapping. The magnitude of this overlap is a measure for the strength of the rebinding effect, caused by the projection and stabilizing the system. A minimal bound for the rebinding effect included in a given system is computed as the solution of an optimization problem. Based on membership functions chosen as a linear combination of Schur vectors, this generalized approach includes reversible as well as non-reversible processes.
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.
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.
Computation of temperature-dependent dissociation rates of metastable protein–ligand complexes
(2019)
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.
Raman spectroscopy is a well established tool for the analysis of vibration spectra, which then allow for the determination of individual substances in a chemical sample, or for their phase transitions. In the Time-Resolved-Raman-Sprectroscopy the vibration spectra of a chemical sample are recorded sequentially over a time interval, such that conclusions for intermediate products (transients) can be drawn within a chemical process. The observed data-matrix M from a Raman spectroscopy can be regarded as a matrix product of two unknown matrices W and H, where the first is representing the contribution of the spectra and the latter represents the chemical spectra. One approach for obtaining W and H is the non-negative matrix factorization. We propose a novel approach, which does not need the commonly used separability assumption. The performance of this approach is shown on a real world chemical example.
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.
The Wiseman fitting can be used to extract binding parameters from ITC data sets, such as heat of binding, number of binding sites, and the overall dissociation rate. The classical Wiseman fitting assumes a direct binding process and neglects the possibility of intermediate binding steps. In principle, it only provides thermodynamic information and not the kinetics of the process. In this article we show that a concentration dependent dissociation constant could possibly stem from intermediate binding steps. The mathematical form of this dependency can be exploited with the aid of the Robust Perron Cluster Cluster Analysis method. Our proposed extension of the Wiseman fitting rationalizes the concentration dependency, and can probably also be used to determine the kinetic parameters of intermediate binding steps of a multivalent binding process. The novelty of this paper is to assume that the binding rate varies per titration step due to the change of the ligand concentration and to use this information in the Wiseman fitting. We do not claim to produce the most accurate values of the binding parameters, we rather present a novel method of how to approach multivalent bindings from a different angle.
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.
We consider two disjoint sets of points with a distance metric, or a
proximity function, associated with each set. If each set can be separately
embedded into separate Euclidean spaces, then we provide sufficient conditions
for the two sets to be jointly embedded in one Euclidean space. In this joint
Euclidean embedding, the distances between the points are generated by a
specific relation-preserving function. Consequently, the mutual distances
between two points of the same set are specific qualitative transformations of
their mutual distances in their original space; the pairwise distances between
the points of different sets can be constructed from an arbitrary proximity
function (might require scaling).
Understanding the kinetics between the components of time-resolved spectra is a crucial step in the study of photo-activatedprocesses. However, modeling the kinetics requires usually some a priori knowledge about the system. In our approach, webuild a Markov State Model (MSM) from the spectral data, and obtain a Koopman transition matrix K(t). With genPCCA,an invariant subspace projection, we project the process into its metastable components. The result of the application of gen-PCCA is a transition matrix Kc(t), from which we can read the transition probability between the metastable components of the reaction. We discuss the application of this analysis method to the transient absorption spectrum of brominated Al-corrole
This article addresses the problem of estimating the Koopman generator of a Markov process. The direct computation of the infinitesimal generator is not easy because of the discretization of the state space, in particular because of the trade-off inherent in the choice of the best lag time to study the process. Short lag times implies a strong discretization of the state space and a consequent loss of Markovianity. Large lag times bypass events on fast timescales. We propose a method to approximate the generator with the computation of the Newton polynomial extrapolation. This technique is a multistep approach which uses as its input Koopman transfer operators evaluated for a series of lag times. Thus, the estimated infinitesimal generator combines information from different time resolutions and does not bias only fast- or slow-decaying dynamics. We show that the multi-scale Newton method can improve the estimation of the generator in comparison to the computation using finite difference or matrix logarithm methods.
Crystallization is a complex phenomenon with far-reaching implications for the production and formulation of active pharmaceutical ingredients. Understanding this process is critical for achieving control over key physicochemical properties that can affect, for example, the bioavailability and stability of a drug. In this study, we were able to reveal intricate and diverse dynamics of the formation of metastable intermediates of paracetamol crystallization varying with the choice of solvent. We demonstrate the efficacy of our novel approach utilizing an objective function-based non-negative matrix factorization technique for the analysis of time-resolved Raman spectroscopy data, in conjunction with time-lapse photography. Furthermore, we emphasize the crucial importance of integrating Raman spectroscopy with supplementary experimental instrumentation for the mathematical analysis of the obtained spectra.
Markov state models (MSMs) have received an unabated increase in popularity in recent years, as they are very
well suited for the identification and analysis of metastable states and related kinetics. However, the state-of-the-art Markov state modeling methods and tools enforce the fulfillment of a
detailed balance condition, restricting their applicability to equilibrium MSMs. To date, they are unsuitable to deal with
general dominant data structures including cyclic processes, which are essentially associated with nonequilibrium systems.
To overcome this limitation, we developed a generalization of the common robust Perron Cluster Cluster Analysis (PCCA+) method, termed generalized PCCA (G-PCCA). This method handles equilibrium and nonequilibrium simulation data, utilizing Schur vectors instead of eigenvectors. G-PCCA is not limited to the detection of metastable states but enables the identification of dominant structures in a general sense, unraveling cyclic processes. This is exemplified by application of G-PCCA on nonequilibrium molecular dynamics data of the Amyloid β (1−40) peptide, periodically driven by an oscillating electric field.
Ergopeptides, like ergocornine and a-ergocryptine, exist in an S- and in an R-configuration. Kinetic experiments imply that certain configurations are preferred depending on the solvent. The experimental methods are explained in this article. Furthermore, computational methods are used to understand this configurational preference. Standard quantum chemical methods can predict the favored configurations by using minimum energy calculations on the potential energy landscape. However, the explicit role of the solvent is not revealed by this type of methods. In order to better understand its influence, classical mechanical molecular simulations are applied. It appears from our research that “folding” the ergopeptide molecules into an intermediate state (between the S- and the R-configuration) is mechanically hindered for the preferred configurations.
Medizin aus dem Computer
(2013)
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.
Applications of the cross-entropy method to importance sampling and optimal control of diffusions
(2014)
In this paper we discuss the notion of research data for the field of mathematics and report on the status quo of research-data management and planning. A number of decentralized approaches are presented and compared to needs and challenges faced in three use cases from different mathematical subdisciplines. We highlight the importance of tailoring research-data management plans to mathematicians’ research processes and discuss their usage all along the data life cycle.