@article{SeckerFackeldeyWeberetal., author = {Secker, Christopher and Fackeldey, Konstantin and Weber, Marcus and Ray, Sourav and Gorgulla, Christoph and Sch{\"u}tte, Christof}, title = {Novel multi-objective affinity approach allows to identify pH-specific μ-opioid receptor agonists}, series = {Journal of Cheminformatics}, volume = {15}, journal = {Journal of Cheminformatics}, doi = {10.1186/s13321-023-00746-4}, abstract = {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.}, language = {en} } @article{ThiesSunkaraRayetal., author = {Thies, Arne and Sunkara, Vikram and Ray, Sourav and Wulkow, Hanna and Celik, M. {\"O}zg{\"u}r and Yerg{\"o}z, Fatih and Sch{\"u}tte, Christof and Stein, Christoph and Weber, Marcus and Winkelmann, Stefanie}, title = {Modelling altered signalling of G-protein coupled receptors in inflamed environment to advance drug design}, series = {Scientific Reports}, volume = {13}, journal = {Scientific Reports}, number = {607}, doi = {10.1038/s41598-023-27699-w}, abstract = {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.}, language = {en} } @article{RaySunkaraSchuetteetal.2020, author = {Ray, Sourav and Sunkara, Vikram and Sch{\"u}tte, Christof and Weber, Marcus}, title = {How to calculate pH-dependent binding rates for receptor-ligand systems based on thermodynamic simulations with different binding motifs}, series = {Molecular Simulation}, volume = {46}, journal = {Molecular Simulation}, number = {18}, publisher = {Taylor and Francis}, doi = {10.1080/08927022.2020.1839660}, pages = {1443 -- 1452}, year = {2020}, abstract = {Molecular simulations of ligand-receptor interactions are a computational challenge, especially when their association- ('on'-rate) and dissociation- ('off'-rate) mechanisms are working on vastly differing timescales. One way of tackling this multiscale problem is to compute the free-energy landscapes, where molecular dynamics (MD) trajectories are used to only produce certain statistical ensembles. The approach allows for deriving the transition rates between energy states as a function of the height of the activation-energy barriers. In this article, we derive the association rates of the opioids fentanyl and N-(3-fluoro-1-phenethylpiperidin-4-yl)-N-phenyl propionamide (NFEPP) in a μ-opioid receptor by combining the free-energy landscape approach with the square-root-approximation method (SQRA), which is a particularly robust version of Markov modelling. The novelty of this work is that we derive the association rates as a function of the pH level using only an ensemble of MD simulations. We also verify our MD-derived insights by reproducing the in vitro study performed by the Stein Lab.}, language = {en} } @article{RayFackeldeySteinetal., author = {Ray, Sourav and Fackeldey, Konstantin and Stein, Christoph and Weber, Marcus}, title = {Coarse Grained MD Simulations of Opioid interactions with the µ-opioid receptor and the surrounding lipid membrane}, series = {Biophysica}, volume = {3}, journal = {Biophysica}, number = {2}, doi = {10.3390/biophysica3020017}, pages = {263 -- 275}, abstract = {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.}, language = {en} } @article{RabbenRayWeber, author = {Rabben, Robert Julian and Ray, Sourav and Weber, Marcus}, title = {ISOKANN: Invariant subspaces of Koopman operators learned by a neural network}, series = {The Journal of Chemical Physics}, volume = {153}, journal = {The Journal of Chemical Physics}, number = {11}, doi = {10.1063/5.0015132}, pages = {114109}, abstract = {The problem of determining the rate of rare events in dynamical systems is quite well-known but still difficult to solve. Recent attempts to overcome this problem exploit the fact that dynamic systems can be represented by a linear operator, such as the Koopman operator. Mathematically, the rare event problem comes down to the difficulty in finding invariant subspaces of these Koopman operators K. In this article, we describe a method to learn basis functions of invariant subspaces using an artificial neural Network.}, language = {en} } @article{PrasadTiwariRayetal., author = {Prasad, Anup K. and Tiwari, Chandni and Ray, Sourav and Holden, Stephanie and Armstrong, David A. and Rosengren, K. Johan and Rodger, Alison and Panwar, Ajay S. and Martin, Lisandra L.}, title = {Secondary Structure Transitions for a Family of Amyloidogenic, Antimircobial Uperin 3 Peptides in Contact with Sodium Dodecyl Sulfate}, series = {ChemPlusChem}, volume = {87}, journal = {ChemPlusChem}, doi = {10.1002/cplu.202100408}, pages = {e202100408}, language = {en} }