TY - JOUR A1 - Thies, Arne A1 - Sunkara, Vikram A1 - Ray, Sourav A1 - Wulkow, Hanna A1 - Celik, M. Özgür A1 - Yergöz, Fatih A1 - Schütte, Christof A1 - Stein, Christoph A1 - Weber, Marcus A1 - Winkelmann, Stefanie T1 - Modelling altered signalling of G-protein coupled receptors in inflamed environment to advance drug design JF - Scientific Reports N2 - We previously reported the successful design, synthesis and testing of the prototype opioid painkiller NFEPP that does not elicit adverse side effects. The design process of NFEPP was based on mathematical modelling of extracellular interactions between G-protein coupled receptors (GPCRs) and ligands, recognizing that GPCRs function differently under pathological versus healthy conditions. We now present an additional and novel stochastic model of GPCR function that includes intracellular dissociation of G-protein subunits and modulation of plasma membrane calcium channels and their dependence on parameters of inflamed and healthy tissue (pH, radicals). The model is validated against in vitro experimental data for the ligands NFEPP and fentanyl at different pH values and radical concentrations. We observe markedly reduced binding affinity and calcium channel inhibition for NFEPP at normal pH compared to lower pH, in contrast to the effect of fentanyl. For increasing radical concentrations, we find enhanced constitutive G-protein activation but reduced ligand binding affinity. Assessing the different effects, the results suggest that, compared to radicals, low pH is a more important determinant of overall GPCR function in an inflamed environment. Future drug design efforts should take this into account. Y1 - 2023 U6 - https://doi.org/10.1038/s41598-023-27699-w VL - 13 IS - 607 ER - TY - JOUR A1 - Donati, Luca A1 - Chewle, Surahit A1 - St. Pierre, Dominik A1 - Natarajan, Vijay A1 - Weber, Marcus T1 - Topological analysis reveals multiple pathways in molecular dynamics JF - Journal of Chemical Theory and Computation N2 - Molecular Dynamics simulations are indispensable tools for comprehending the dynamic behavior of biomolecules, yet extracting meaningful molecular pathways from these simulations remains challenging due to the vast amount of high dimensional data. In this work, we present Molecular Kinetics via Topology (MoKiTo), a novel approach that combines the ISOKANN algorithm to determine the membership function of a molecular system with a topological analysis tool inspired by the Mapper algorithm. Our strategy efficiently identifies and characterizes distinct molecular pathways, enabling the detection and visualization of critical conformational transitions and rare events. This method offers deeper insights into molecular mechanisms, facilitating the design of targeted interventions in drug discovery and protein engineering. Y1 - 2025 U6 - https://doi.org/10.1021/acs.jctc.5c00819 VL - 21 IS - 20 SP - 10385 EP - 10397 CY - J. Chem. Theory Comput. ER - TY - BOOK A1 - Ternes, Thomas A1 - Bauer, Karl-Heinz A1 - Brauer, Frank A1 - Drewes, Jörg A1 - Joss, Adriano A1 - Hiller, Georg A1 - Jewell, Kevin A1 - Oehlmann, Jörg A1 - Radke, Michael A1 - Schulte-Oehlmann, Ulrike A1 - Schwartz, Thomas A1 - Seel, Peter A1 - Völker, Jeanette A1 - Weber, Lilo A1 - Weber, Marcus ED - Wilhelm, Christian T1 - Handlungsempfehlung zur integrativen Bewertung der weitergehenden Abwasserbehandlung von kommunalen Kläranlagen N2 - Das vorliegende Statuspapier beschreibt ein Konzept zur weitergehenden Abwasserbehandlung für die Bewertung von Aufbereitungsverfahren, sowohl in einer Pilotphase zur Auswahl von Verfah- rensoptionen als auch für die Bewertung großtechnischer Anlagen. Y1 - 2023 SN - 978-3-96862-563-8 VL - T1/2023 PB - DWA / GDCh ER - TY - JOUR A1 - Donati, Luca A1 - Weber, Marcus T1 - Efficient Estimation of Transition Rates as Functions of pH JF - Proceedings in Applied Mathematics & Mechanics N2 - Extracting the kinetic properties of a system whose dynamics depend on the pH of the environment with which it exchanges energy and atoms requires sampling the Grand Canonical Ensemble. As an alternative, we present a novel strategy that requires simulating only the most recurrent Canonical Ensembles that compose the Grand Canonical Ensemble. The simulations are used to estimate the Gran Canonical distribution for a specific pH value by reweighting and to construct the transition rate matrix by discretizing the Fokker-Planck equation by Square Root Approximation and robust Perron Cluster Cluster Analysis. As an application, we have studied the tripeptide Ala-Asp-Ala. Y1 - 2023 U6 - https://doi.org/10.1002/pamm.202300264 VL - 23 ER - TY - JOUR A1 - Yousefian, Maryam A1 - Frank, Anne-Simone A1 - Weber, Marcus A1 - Röblitz, Susanna T1 - Efficient construction of Markov state models for stochastic gene regulatory networks by domain decomposition JF - BMC Bioinformatics N2 - The dynamics of many gene regulatory networks (GRNs) is characterized by the occurrence of metastable phenotypes and stochastic phenotype switches. The chemical master equation (CME) is the most accurate description to model such stochastic dynamics, whereby the long-time dynamics of the system is encoded in the spectral properties of the CME operator. Markov State Models (MSMs) provide a general framework for analyzing and visualizing stochastic multistability and state transitions based on these spectral properties. Until now, however, this approach is either limited to low-dimensional systems or requires the use of high-performance computing facilities, thus limiting its usability. Y1 - 2025 U6 - https://doi.org/10.1186/s12859-025-06174-5 VL - 26 IS - 147 ER - TY - JOUR A1 - Coomber, Celvic A1 - Chewle, Surahit A1 - Secker, Christopher A1 - Fackeldey, Konstantin A1 - Weber, Marcus A1 - Winkelmann, Stefanie A1 - Schütte, Christof A1 - Sunkara, Vikram T1 - Investigating Endogenous Opioids Unravels the Mechanisms Behind Opioid-Induced Constipation, a Mathematical Modeling Approach JF - International Journal of Molecular Sciences N2 - Endogenous opioids, such as Endomorphin-2, are not typically associated with severe constipation, unlike pharmaceutical opioids, which induce opioid-induced constipation (OIC) by activating μ-opioid receptors in the gastrointestinal tract. In this study, we present a mathematical model, which integrates the serotonergic and opioid pathways, simulating the interaction between serotonin and opioid signaling within the enteric nervous system (ENS). The model explores the mechanisms underlying OIC, with a focus on the change in adenylyl cyclase (AC) activity, cAMP accumulation, and the distinct functionalities of Endomorphin-2 compared to commonly used pharmaceutical opioids. We study the effects of Morphine, Fentanyl, and Methadone and contrast them with Endomorphin-2. Our findings reveal that opioids do not perturb the signaling of serotonin, but only the activity of AC, suggesting that serotonin levels have no influence on improving opioid-induced constipation. Furthermore, this study reveals that the primary difference between endogenous and pharmaceutical opioids is their degradation rates. This finding shows that modulating opioid degradation rates significantly improves cAMP recovery. In conclusion, our insights steer towards exploring opioid degrading enzymes, localized to the gut, as a strategy for mitigating OIC. Y1 - 2025 U6 - https://doi.org/10.3390/ijms26136207 VL - 26 IS - 13 ER - TY - JOUR A1 - Secker, Christopher A1 - Fackeldey, Konstantin A1 - Weber, Marcus A1 - Ray, Sourav A1 - Gorgulla, Christoph A1 - Schütte, Christof T1 - Novel multi-objective affinity approach allows to identify pH-specific μ-opioid receptor agonists JF - Journal of Cheminformatics N2 - Opioids are essential pharmaceuticals due to their analgesic properties, however, lethal side effects, addiction, and opioid tolerance are extremely challenging. The development of novel molecules targeting the μ-opioid receptor (MOR) in inflamed, but not in healthy tissue, could significantly reduce these unwanted effects. Finding such novel molecules can be achieved by maximizing the binding affinity to the MOR at acidic pH while minimizing it at neutral pH, thus combining two conflicting objectives. Here, this multi-objective optimal affinity approach is presented, together with a virtual drug discovery pipeline for its practical implementation. When applied to finding pH-specific drug candidates, it combines protonation state-dependent structure and ligand preparation with high-throughput virtual screening. We employ this pipeline to characterize a set of MOR agonists identifying a morphine-like opioid derivative with higher predicted binding affinities to the MOR at low pH compared to neutral pH. Our results also confirm existing experimental evidence that NFEPP, a previously described fentanyl derivative with reduced side effects, and recently reported β-fluorofentanyls and -morphines show an increased specificity for the MOR at acidic pH when compared to fentanyl and morphine. We further applied our approach to screen a >50K ligand library identifying novel molecules with pH-specific predicted binding affinities to the MOR. The presented differential docking pipeline can be applied to perform multi-objective affinity optimization to identify safer and more specific drug candidates at large scale. Y1 - 2023 U6 - https://doi.org/10.1186/s13321-023-00746-4 VL - 15 ER - TY - JOUR A1 - Donati, Luca A1 - Weber, Marcus T1 - Assessing transition rates as functions of environmental variables JF - The Journal of Chemical Physics N2 - 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. Y1 - 2022 U6 - https://doi.org/10.1063/5.0109555 VL - 157 IS - 22 SP - 224103-1 EP - 224103-14 PB - AIP Publishing ER - TY - JOUR A1 - Sechi, Renata A1 - Fackeldey, Konstantin A1 - Chewle, Surahit A1 - Weber, Marcus T1 - SepFree NMF: A Toolbox for Analyzing the Kinetics of Sequential Spectroscopic Data JF - Algorithms N2 - 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. Y1 - 2022 U6 - https://doi.org/10.3390/a15090297 VL - 15 IS - 9 SP - 297 ER - TY - JOUR A1 - Donati, Luca A1 - Weber, Marcus A1 - Keller, Bettina G. T1 - A review of Girsanov Reweighting and of Square Root Approximation for building molecular Markov State Models JF - Journal of Mathematical Physics N2 - 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. Y1 - 2022 U6 - https://doi.org/10.1063/5.0127227 VL - 63 IS - 12 SP - 123306-1 EP - 123306-21 PB - AIP Publishing ER - TY - JOUR A1 - Ray, Sourav A1 - Fackeldey, Konstantin A1 - Stein, Christoph A1 - Weber, Marcus T1 - Coarse Grained MD Simulations of Opioid interactions with the µ-opioid receptor and the surrounding lipid membrane JF - Biophysica N2 - In our previous studies, a new opioid (NFEPP) was developed to only selectively bind to the 𝜇-opoid receptor (MOR) in inflamed tissue and thus avoid the severe side effects of fentanyl. We know that NFEPP has a reduced binding affinity to MOR in healthy tissue. Inspired by the modelling and simulations performed by Sutcliffe et al., we present our own results of coarse-grained molecular dynamics simulations of fentanyl and NFEPP with regards to their interaction with the 𝜇-opioid receptor embedded within the lipid cell membrane. For technical reasons, we have slightly modified Sutcliffe’s parametrisation of opioids. The pH-dependent opioid simulations are of interest because while fentanyl is protonated at the physiological pH, NFEPP is deprotonated due to its lower pKa value than that of fentanyl. Here, we analyse for the first time whether pH changes have an effect on the dynamical behaviour of NFEPP when it is inside the cell membrane. Besides these changes, our analysis shows a possible alternative interaction of NFEPP at pH 7.4 outside the binding region of the MOR. The interaction potential of NFEPP with MOR is also depicted by analysing the provided statistical molecular dynamics simulations with the aid of an eigenvector analysis of a transition rate matrix. In our modelling, we see differences in the XY-diffusion profiles of NFEPP compared with fentanyl in the cell membrane. Y1 - 2023 U6 - https://doi.org/10.3390/biophysica3020017 VL - 3 IS - 2 SP - 263 EP - 275 ER - TY - JOUR A1 - Sikorski, Alexander A1 - Niknejad, Amir A1 - Weber, Marcus A1 - Donati, Luca T1 - Tensor-SqRA: Modeling the transition rates of interacting molecular systems in terms of potential energies JF - Journal of Chemical Physics N2 - Estimating the rate of rare conformational changes in molecular systems is one of the goals of molecular dynamics simulations. In the past few decades, a lot of progress has been done in data-based approaches toward this problem. In contrast, model-based methods, such as the Square Root Approximation (SqRA), directly derive these quantities from the potential energy functions. In this article, we demonstrate how the SqRA formalism naturally blends with the tensor structure obtained by coupling multiple systems, resulting in the tensor-based Square Root Approximation (tSqRA). It enables efficient treatment of high-dimensional systems using the SqRA and provides an algebraic expression of the impact of coupling energies between molecular subsystems. Based on the tSqRA, we also develop the projected rate estimation, a hybrid data-model-based algorithm that efficiently estimates the slowest rates for coupled systems. In addition, we investigate the possibility of integrating low-rank approximations within this framework to maximize the potential of the tSqRA. Y1 - 2024 U6 - https://doi.org/10.1063/5.0187792 VL - 160 SP - 104112 ER - TY - JOUR A1 - Donati, Luca A1 - Schütte, Christof A1 - Weber, Marcus T1 - The Kramers turnover in terms of a macro-state projection on phase space JF - Molecular Physics N2 - We have investigated how Langevin dynamics is affected by the friction coefficient using the novel algorithm ISOKANN, which combines the transfer operator approach with modern machine learning techniques. ISOKANN describes the dynamics in terms of an invariant subspace projection of the Koopman operator defined in the entire state space, avoiding approximations due to dimensionality reduction and discretization. Our results are consistent with the Kramers turnover and show that in the low and moderate friction regimes, metastable macro-states and transition rates are defined in phase space, not only in position space. Y1 - 2025 U6 - https://doi.org/10.1080/00268976.2024.2356748 VL - 123 IS - 7-8: Giovanni Ciccotti: A Renaissance Physicist SP - e2356748 PB - Taylor & Francis ER -