Computation-guided analysis of paroxetine binding to hSERT reveals functionally important structural elements and dynamics

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Dynamic extracellular vestibule of human SERT: Unveiling druggable potential with high-affinity allosteric inhibitors

The serotonin transporter (SERT) tightly regulates synaptic serotonin levels and has been the primary target of antidepressants. Binding of inhibitors to the allosteric site of human SERT (hSERT) impedes the dissociation of antidepressants bound at the central site and may enhance the efficacy of such antidepressants to potentially reduce their dosage and side effects. Here, we report the identification of a series of high-affinity allosteric inhibitors of hSERT in a unique scaffold, with the lead compound, Lu AF88273 (3-(1-(2-(1H-indol-3-yl)ethyl)piperidin-4-yl)-6-chloro-1H-indole), having 2.1 nM allosteric potency in inhibiting imipramine dissociation. In addition, we find that Lu AF88273 also inhibits serotonin transport in a noncompetitive manner. The binding pose of Lu AF88273 in the allosteric site of hSERT is determined with extensive molecular dynamics simulations and rigorous absolute binding free energy perturbation (FEP) calculations, which show that a part of the compound occupies a dynamically formed small cavity. The predicted binding location and pose are validated by site-directed mutagenesis and can explain much of the structure-activity relationship of these inhibitors using the relative binding FEP calculations. Together, our findings provide a promising lead compound and the structural basis for the development of allosteric drugs targeting hSERT. Further, they demonstrate that the divergent allosteric sites of neurotransmitter transporters can be selectively targeted.

Design, synthesis and cytotoxic evaluation of a selective serotonin reuptake inhibitor (SSRI) by virtual screening

Depression is one of the most common mental illnesses, affecting almost 300 million people. According to the WHO, depression is one of the world's leading causes of disability and morbidity. People with this illness require both psychological and pharmaceutical treatment because severe depressive episodes often result in suicide. Selective serotonin reuptake inhibitors (SSRI) are widely used antidepressants that target the human serotonin transporter (hSERT). The crystallization of hSERT and the experimental data available allows cost and time-efficient computational tools like virtual screening (VS) to be utilized in the development of therapeutic agents. Here, we synthesized, characterized, and evaluated the biological activity of a novel SSRI analog of paroxetine, rationally designed by applying an artificial neural network-based QSAR model and a molecular docking analysis on hSERT. The analog N-substituted 18a showed higher affinity for the transporter (-10.2 kcal/mol), lower Ki value (1.19 nM) and a safer toxicological profile than paroxetine and was synthesized with a 71% yield. The in vitro cytotoxicity of the analog was evaluated using human glioblastoma (U87 MG), human neuroblastoma (SH SY5Y) and murine fibroblast (L929) cell lines. Also, the hemolytic ability of the compound was assessed on human erythrocytes. Results showed that analog 18a did not exhibit cytotoxic activity on the cell lines used and has no hemolytic activity at any of the concentrations tested, whereas with paroxetine, hemolysis was observed at 2.3, 1.29 y 0.67 mM. Based on these results, it is possible to suggest that analog 18a could be a promising new SSRI candidate for the treatment of this illness.

Unravelling the Turn-On Fluorescence Mechanism of a Fluorescein-Based Probe in GABAA Receptors

GABAA (γ-aminobutyric acid type A) receptors are ligand-gated ion channels mediating fast inhibitory transmission in the mammalian brain. Here we report the molecular and electronic mechanism governing the turn-on emission of a fluorescein-based imaging probe able to target the human GABAA receptor. Multiscale calculations evidence a drastic conformational change of the probe from folded in solution to extended upon binding to the receptor. Intramolecular ππ-stacking interactions present in the folded probe are responsible for quenching fluorescence in solution. In contrast, unfolding within the GABAA receptor changes the nature of the bright excited state triggering emission. Remarkably, this turn-on effect only manifests for the dianionic prototropic form of the imaging probe, which is found to be the strongest binder to the GABAA receptor. This study is expected to assist the design of new photoactivatable screening tools for allosteric modulators of the GABAA receptor.

Unravelling the Turn‐On Fluorescence Mechanism of a Fluorescein‐Based Probe in GABA A Receptors

GABA_A (γ‐aminobutyric acid type A) receptors are ligand‐gated ion channels mediating fast inhibitory transmission in the mammalian brain. Here we report the molecular and electronic mechanism governing the turn‐on emission of a fluorescein‐based imaging probe able to target the human GABA_A receptor. Multiscale calculations evidence a drastic conformational change of the probe from folded in solution to extended upon binding to the receptor. Intramolecular ππ‐stacking interactions present in the folded probe are responsible for quenching fluorescence in solution. In contrast, unfolding within the GABA_A receptor changes the nature of the bright excited state triggering emission. Remarkably, this turn‐on effect only manifests for the dianionic prototropic form of the imaging probe, which is found to be the strongest binder to the GABA_A receptor. This study is expected to assist the design of new photoactivatable screening tools for allosteric modulators of the GABA_A receptor.

Recent developments in horizontal gene transfer with the adaptive innovation of fermented foods

Horizontal gene transfer (HGT) has contributed significantly to the adaptability of bacteria, yeast and mold in fermented foods, whose evidence has been found in several fermented foods. Although not every HGT has biological significance, it plays an important role in improving the quality of fermented foods. In this review, how HGT facilitated microbial domestication and adaptive evolution in fermented foods was discussed. HGT can assist in the industrial innovation of fermented foods, and this adaptive evolution strategy can improve the quality of fermented foods. Additionally, the mechanism underlying HGT in fermented foods were analyzed. Furthermore, the critical bottlenecks involved in optimizing HGT during the production of fermented foods and strategies for optimizing HGT were proposed. Finally, the prospect of HGT for promoting the industrial innovation of fermented foods was highlighted. The comprehensive report on HGT in fermented foods provides a new trend for domesticating preferable starters for food fermentation, thus optimizing the quality and improving the industrial production of fermented foods.

A comprehensive evaluation of the potential binding poses of fentanyl and its analogs at the µ-opioid receptor

Fentanyl and its analogs are selective agonists of the µ-opioid receptor (MOR). Among novel synthetic opioids (NSOs), they dominate the recreational drug market and are the main culprits for the opioid crisis, which has been exacerbated by the COVID-19 pandemic. By taking advantage of the crystal structures of the MOR, several groups have investigated the binding mechanism of fentanyl, but have not reached a consensus, in terms of both the binding orientation and the fentanyl conformation. Thus, the binding mechanism of fentanyl at the MOR remains an unsolved and challenging question. Here, we carried out a systematic computational study to investigate the preferred fentanyl conformations, and how these conformations are being accommodated in the MOR binding pocket. We characterized the free energy landscape of fentanyl conformations with metadynamics simulations, and compared and evaluated several possible fentanyl binding conditions in the MOR with long-timescale molecular dynamics simulations. Our results indicate that the most preferred binding pose in the MOR binding pocket corresponds well with the global minimum on the energy landscape of fentanyl in the absence of the receptor, while the energy landscape can be reconfigured by modifying the fentanyl scaffold. The interactions with the receptor may stabilize a slightly unfavored fentanyl conformation in an alternative binding pose. By extending similar investigations to fentanyl analogs, our findings establish a structure–activity relationship of fentanyl binding at the MOR. In addition to providing a structural basis to understand the potential toxicity of the emerging NSOs, such insights will contribute to developing new, safer analgesics.

Molecular Mechanism for the Allosteric Inhibition of the Human Serotonin Transporter by Antidepressant Escitalopram

Human serotine transporter (hSERT) is one of the most influential drug targets, and its allosteric modulators (e.g., escitalopram) have emerged to be the next-generation medication for psychiatric disorders. However, the molecular mechanism underlying the allosteric modulation of hSERT is still elusive. Here, the simulation strategies of conventional (cMD) and steered (SMD) molecular dynamics were applied to investigate this molecular mechanism from distinct perspectives. First, cMD simulations revealed that escitalopram's binding to hSERT's allosteric site simultaneously enhanced its binding to the orthosteric site. Then, SMD simulation identified that the occupation of hSERT's allosteric site by escitalopram could also block its dissociation from the orthosteric site. Finally, by comparing the simulated structures of two hSERT-escitalopram complexes with and without allosteric modulation, a new conformational coupling between an extracellular (Arg104-Glu494) and an intracellular (Lys490-Glu494) salt bridge was identified. In summary, this study explored the mechanism underlying the allosteric modulation of hSERT by collectively applying two MD simulation strategies, which could facilitate our understanding of the allosteric modulations of not only hSERT but also other clinically important therapeutic targets.

Toward Reducing hERG Affinities for DAT Inhibitors with a Combined Machine Learning and Molecular Modeling Approach

Psychostimulant drugs, such as cocaine, inhibit dopamine reuptake via blockading the dopamine transporter (DAT), which is the primary mechanism underpinning their abuse. Atypical DAT inhibitors are dissimilar to cocaine and can block cocaine- or methamphetamine-induced behaviors, supporting their development as part of a treatment regimen for psychostimulant use disorders. When developing these atypical DAT inhibitors as medications, it is necessary to avoid off-target binding that can produce unwanted side effects or toxicities. In particular, the blockade of a potassium channel, human ether-a-go-go (hERG), can lead to potentially lethal ventricular tachycardia. In this study, we established a counter screening platform for DAT and against hERG binding by combining machine learning-based quantitative structure-activity relationship (QSAR) modeling, experimental validation, and molecular modeling and simulations. Our results show that the available data are adequate to establish robust QSAR models, as validated by chemical synthesis and pharmacological evaluation of a validation set of DAT inhibitors. Furthermore, the QSAR models based on subsets of the data according to experimental approaches used have predictive power as well, which opens the door to target specific functional states of a protein. Complementarily, our molecular modeling and simulations identified the structural elements responsible for a pair of DAT inhibitors having opposite binding affinity trends at DAT and hERG, which can be leveraged for rational optimization of lead atypical DAT inhibitors with desired pharmacological properties.

Assessment of Paroxetine Molecular Interactions with Selected Monoamine and γ-Aminobutyric Acid Transporters

Thus far, many hypotheses have been proposed explaining the cause of depression. Among the most popular of these are: monoamine, neurogenesis, neurobiology, inflammation and stress hypotheses. Many studies have proven that neurogenesis in the brains of adult mammals occurs throughout life. The generation of new neurons persists throughout adulthood in the mammalian brain due to the proliferation and differentiation of adult neural stem cells. For this reason, the search for drugs acting in this mechanism seems to be a priority for modern pharmacotherapy. Paroxetine is one of the most commonly used antidepressants. However, the exact mechanism of its action is not fully understood. The fact that the therapeutic effect after the administration of paroxetine occurs after a few weeks, even if the levels of monoamine are rapidly increased (within a few minutes), allows us to assume a neurogenic mechanism of action. Due to the confirmed dependence of depression on serotonin, norepinephrine, dopamine and γ-aminobutyric acid levels, studies have been undertaken into paroxetine interactions with these primary neurotransmitters using in silico and in vitro methods. We confirmed that paroxetine interacts most strongly with monoamine transporters and shows some interaction with γ-aminobutyric acid transporters. However, studies of the potency inhibitors and binding affinity values indicate that the neurogenic mechanism of paroxetine's action may be determined mainly by its interactions with serotonin transporters.

Tuning the activity of known drugs via the introduction of halogen atoms, a case study of SERT ligands - Fluoxetine and fluvoxamine

The selective serotonin reuptake inhibitors (SSRIs), acting at the serotonin transporter (SERT), are one of the most widely prescribed antidepressant medications. All five approved SSRIs possess either fluorine or chlorine atoms, and there is a limited number of reports describing their analogs with heavier halogens, i.e., bromine and iodine. To elucidate the role of halogen atoms in the binding of SSRIs to SERT, we designed a series of 22 fluoxetine and fluvoxamine analogs substituted with fluorine, chlorine, bromine, and iodine atoms, differently arranged on the phenyl ring. The obtained biological activity data, supported by a thorough in silico binding mode analysis, allowed the identification of two partners for halogen bond interactions: the backbone carbonyl oxygen atoms of E493 and T497. Additionally, compounds with heavier halogen atoms were found to bind with the SERT via a distinctly different binding mode, a result not presented elsewhere. The subsequent analysis of the prepared XSAR sets showed that E493 and T497 participated in the largest number of formed halogen bonds. The XSAR library analysis led to the synthesis of two of the most active compounds (3,4-diCl-fluoxetine 42, SERT K_i = 5 nM and 3,4-diCl-fluvoxamine 46, SERT K_i = 9 nM, fluoxetine SERT K_i = 31 nM, fluvoxamine SERT K_i = 458 nM). We present an example of the successful use of a rational methodology to analyze binding and design more active compounds by halogen atom introduction. 'XSAR library analysis', a new tool in medicinal chemistry, was instrumental in identifying optimal halogen atom substitution.

Paroxetine-Overview of the Molecular Mechanisms of Action

In the 21st century and especially during a pandemic, the diagnosis and treatment of depression is an essential part of the daily practice of many family doctors. It mainly affects patients in the age category 15–44 years, regardless of gender. Anxiety disorders are often diagnosed in children and adolescents. Social phobias can account for up to 13% of these diagnoses. Social anxiety manifests itself in fear of negative social assessment and humiliation, which disrupts the quality of social functioning. Treatment of the above-mentioned disorders is based on psychotherapy and pharmacotherapy. Serious side effects or mortality from antidepressant drug overdose are currently rare. Recent studies indicate that paroxetine (ATC code: N06AB), belonging to the selective serotonin reuptake inhibitors, has promising therapeutic effects and is used off-label in children and adolescents. The purpose of this review is to describe the interaction of paroxetine with several molecular targets in various points of view including the basic chemical and pharmaceutical properties. The central point of the review is focused on the pharmacodynamic analysis based on the molecular mechanism of binding paroxetine to various therapeutic targets.

Chemical and structural investigation of the paroxetine-human serotonin transporter complex

Antidepressants target the serotonin transporter (SERT) by inhibiting serotonin reuptake. Structural and biochemical studies aiming to understand binding of small-molecules to conformationally dynamic transporters like SERT often require thermostabilizing mutations and antibodies to stabilize a specific conformation, leading to questions about relationships of these structures to the bonafide conformation and inhibitor binding poses of wild-type transporter. To address these concerns, we determined the structures of ∆N72/∆C13 and ts2-inactive SERT bound to paroxetine analogues using single-particle cryo-EM and x-ray crystallography, respectively. We synthesized enantiopure analogues of paroxetine containing either bromine or iodine instead of fluorine. We exploited the anomalous scattering of bromine and iodine to define the pose of these inhibitors and investigated inhibitor binding to Asn177 mutants of ts2-active SERT. These studies provide mutually consistent insights into how paroxetine and its analogues bind to the central substrate-binding site of SERT, stabilize the outward-open conformation, and inhibit serotonin transport.

The mechanism of a high-affinity allosteric inhibitor of the serotonin transporter

The serotonin transporter (SERT) terminates serotonin signaling by rapid presynaptic reuptake. SERT activity is modulated by antidepressants, e.g., S-citalopram and imipramine, to alleviate symptoms of depression and anxiety. SERT crystal structures reveal two S-citalopram binding pockets in the central binding (S1) site and the extracellular vestibule (S2 site). In this study, our combined in vitro and in silico analysis indicates that the bound S-citalopram or imipramine in S1 is allosterically coupled to the ligand binding to S2 through altering protein conformations. Remarkably, SERT inhibitor Lu AF60097, the first high-affinity S2-ligand reported and characterized here, allosterically couples the ligand binding to S1 through a similar mechanism. The SERT inhibition by Lu AF60097 is demonstrated by the potentiated imipramine binding and increased hippocampal serotonin level in rats. Together, we reveal a S1-S2 coupling mechanism that will facilitate rational design of high-affinity SERT allosteric inhibitors.

A Novel Bromine-Containing Paroxetine Analogue Provides Mechanistic Clues for Binding Ambiguity at the Central Primary Binding Site of the Serotonin Transporter

The serotonin transporter (SERT) is the primary target for the selective serotonin reuptake inhibitors (SSRIs). However, the structural basis for the extraordinarily high binding affinity of the widely prescribed SSRI, paroxetine, to human SERT (hSERT) has not yet been fully elucidated. Our previous findings unveiled a plausible ambiguity in paroxetine's binding orientations that may constitute an integral component of this SSRI's high affinity for hSERT. Herein, we investigate factors contributing to paroxetine's high affinity by modifying both the ligand and the protein. We generated a series of bromine (Br)-containing derivatives and found that the one in which the 4-F of paroxetine had been replaced with the chemically similar but more electron-rich Br atom (13) had the highest affinity. By comparatively characterizing the binding of paroxetine and 13 to both wild type (WT) and a construct harboring a paroxetine-sensitive mutation in the binding cavity, we identified a mechanistic determinant responsible for the pose ambiguity of paroxetine, which can guide future drug design.

The allosteric mechanism of substrate-specific transport in SLC6 is mediated by a volumetric sensor

Neurotransmitter:sodium symporters (NSSs) in the SLC6 family terminate neurotransmission by coupling the thermodynamically favorable transport of ions to the thermodynamically unfavorable transport of neurotransmitter back into presynaptic neurons. Results from many structural, functional, and computational studies on LeuT, a bacterial NSS homolog, have provided critical insight into the mechanism of sodium-coupled transport, but the mechanism underlying substrate-specific transport rates is still not understood. We present a combination of molecular dynamics simulations, single-molecule fluorescence resonance energy transfer (smFRET) imaging, and measurements of Na+ binding and substrate transport that reveals an allosteric substrate specificity mechanism. In this mechanism, residues F259 and I359 in the substrate binding pocket couple the binding of substrate to Na+ release from the Na2 site by allosterically modulating the stability of a partially open, inward-facing state. We propose a model for transport selectivity in which residues F259 and I359 act as a volumetric sensor that inhibits the transport of bulky amino acids.

Monoamine transporters: structure, intrinsic dynamics and allosteric regulation

Monoamine transporters (MATs) regulate neurotransmission via the reuptake of dopamine, serotonin and norepinephrine from extra-neuronal regions and thus maintain neurotransmitter homeostasis. As targets of a wide range of compounds, including antidepressants, substances of abuse and drugs for neuropsychiatric and neurodegenerative disorders, their mechanism of action and their modulation by small molecules have long been of broad interest. Recent advances in the structural characterization of dopamine and serotonin transporters have opened the way for structure-based modeling and simulations, which, together with experimental data, now provide mechanistic understanding of their transport function and interactions. Here we review recent progress in the elucidation of the structural dynamics of MATs and their conformational landscape and transitions, as well as allosteric regulation mechanisms.

Insights into the mechanism and pharmacology of neurotransmitter sodium symporters

Neurotransmitter sodium symporters (NSS) belong to the SLC6 family of solute carriers and play an essential role in neurotransmitter homeostasis throughout the body. In the past decade, structural studies employing bacterial orthologs of NSSs have provided insight into the mechanism of neurotransmitter transport. While the overall architecture of SLC6 transporters is conserved among species, in comparison to the bacterial homologs, the eukaryotic SLC6 family members harbor differences in amino acid sequence and molecular structure, which underpins their functional and pharmacological diversity, as well as their ligand specificity. Here, we review the structures and mechanisms of eukaryotic NSSs, focusing on the molecular basis for ligand recognition and on transport mechanism.