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 - Montefusco, Alberto A1 - Schütte, Christof A1 - Winkelmann, Stefanie T1 - A route to the hydrodynamic limit of a reaction-diffusion master equation using gradient structures JF - SIAM Journal on Applied Mathematics N2 - The reaction-diffusion master equation (RDME) is a lattice-based stochastic model for spatially resolved cellular processes. It is often interpreted as an approximation to spatially continuous reaction-diffusion models, which, in the limit of an infinitely large population, may be described by means of reaction-diffusion partial differential equations. Analyzing and understanding the relation between different mathematical models for reaction-diffusion dynamics is a research topic of steady interest. In this work, we explore a route to the hydrodynamic limit of the RDME which uses gradient structures. Specifically, we elaborate on a method introduced in [J. Maas and A. Mielke, J. Stat. Phys., 181 (2020), pp. 2257–2303] in the context of well-mixed reaction networks by showing that, once it is complemented with an appropriate limit procedure, it can be applied to spatially extended systems with diffusion. Under the assumption of detailed balance, we write down a gradient structure for the RDME and use the method in order to produce a gradient structure for its hydrodynamic limit, namely, for the corresponding RDPDE. Y1 - 2023 U6 - https://doi.org/10.1137/22M1488831 VL - 83 IS - 2 SP - 837 EP - 861 ER - TY - JOUR A1 - Wehlitz, Nathalie A1 - Sadeghi, Mohsen A1 - Montefusco, Alberto A1 - Schütte, Christof A1 - Pavliotis, Grigorios A. A1 - Winkelmann, Stefanie T1 - Approximating particle-based clustering dynamics by stochastic PDEs JF - SIAM Journal on Applied Dynamical Systems N2 - This work proposes stochastic partial differential equations (SPDEs) as a practical tool to replicate clustering effects of more detailed particle-based dynamics. Inspired by membrane mediated receptor dynamics on cell surfaces, we formulate a stochastic particle-based model for diffusion and pairwise interaction of particles, leading to intriguing clustering phenomena. Employing numerical simulation and cluster detection methods, we explore the approximation of the particle-based clustering dynamics through mean-field approaches. We find that SPDEs successfully reproduce spatiotemporal clustering dynamics, not only in the initial cluster formation period, but also on longer time scales where the successive merging of clusters cannot be tracked by deterministic mean-field models. The computational efficiency of the SPDE approach allows us to generate extensive statistical data for parameter estimation in a simpler model that uses a Markov jump process to capture the temporal evolution of the cluster number. Y1 - 2025 U6 - https://doi.org/10.1137/24M1676661 VL - 24 IS - 2 SP - 1231 EP - 1250 ER - TY - JOUR A1 - Schütte, Christof A1 - Klus, Stefan A1 - Hartmann, Carsten T1 - Overcoming the Timescale Barrier in Molecular Dynamics: Transfer Operators, Variational Principles, and Machine Learning JF - Acta Numerica N2 - One of the main challenges in molecular dynamics is overcoming the ‘timescale barrier’: in many realistic molecular systems, biologically important rare transitions occur on timescales that are not accessible to direct numerical simulation, even on the largest or specifically dedicated supercomputers. This article discusses how to circumvent the timescale barrier by a collection of transfer operator-based techniques that have emerged from dynamical systems theory, numerical mathematics and machine learning over the last two decades. We will focus on how transfer operators can be used to approximate the dynamical behaviour on long timescales, review the introduction of this approach into molecular dynamics, and outline the respective theory, as well as the algorithmic development, from the early numerics-based methods, via variational reformulations, to modern data-based techniques utilizing and improving concepts from machine learning. Furthermore, its relation to rare event simulation techniques will be explained, revealing a broad equivalence of variational principles for long-time quantities in molecular dynamics. The article will mainly take a mathematical perspective and will leave the application to real-world molecular systems to the more than 1000 research articles already written on this subject. Y1 - 2023 U6 - https://doi.org/10.1017/S0962492923000016 VL - 32 SP - 517 EP - 673 ER - TY - JOUR A1 - Montefusco, Alberto A1 - Helfmann, Luzie A1 - Okunola, Toluwani A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - Partial mean-field model for neurotransmission dynamics JF - Mathematical Biosciences N2 - This article addresses reaction networks in which spatial and stochastic effects are of crucial importance. For such systems, particle-based models allow us to describe all microscopic details with high accuracy. However, they suffer from computational inefficiency if particle numbers and density get too large. Alternative coarse-grained-resolution models reduce computational effort tremendously, e.g., by replacing the particle distribution by a continuous concentration field governed by reaction-diffusion PDEs. We demonstrate how models on the different resolution levels can be combined into hybrid models that seamlessly combine the best of both worlds, describing molecular species with large copy numbers by macroscopic equations with spatial resolution while keeping the stochastic-spatial particle-based resolution level for the species with low copy numbers. To this end, we introduce a simple particle-based model for the binding dynamics of ions and vesicles at the heart of the neurotransmission process. Within this framework, we derive a novel hybrid model and present results from numerical experiments which demonstrate that the hybrid model allows for an accurate approximation of the full particle-based model in realistic scenarios. Y1 - 2024 U6 - https://doi.org/10.1016/j.mbs.2024.109143 VL - 369 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 - GEN A1 - Schütte, Christof A1 - Klus, Stefan A1 - Hartmann, Carsten T1 - Overcoming the Timescale Barrier in Molecular Dynamics: Transfer Operators, Variational Principles, and Machine Learning N2 - One of the main challenges in molecular dynamics is overcoming the “timescale barrier”, a phrase used to describe that in many realistic molecular systems, biologically important rare transitions occur on timescales that are not accessible to direct numerical simulation, not even on the largest or specifically dedicated supercomputers. This article discusses how to circumvent the timescale barrier by a collection of transfer operator-based techniques that have emerged from dynamical systems theory, numerical mathematics, and machine learning over the last two decades. We will focus on how transfer operators can be used to approximate the dynamical behavior on long timescales, review the introduction of this approach into molecular dynamics, and outline the respective theory as well as the algorithmic development from the early numerics-based methods, via variational reformulations, to modern data-based techniques utilizing and improving concepts from machine learning. Furthermore, its relation to rare event simulation techniques will be explained, revealing a broad equivalence of variational principles for long-time quantities in MD. The article will mainly take a mathematical perspective and will leave the application to real-world molecular systems to the more than 1000 research articles already written on this subject. T3 - ZIB-Report - 22-25 Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-88637 SN - 1438-0064 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 -