ABSTRACT
Single-cell (sc)RNA-seq, together with RNA velocity and metabolic labeling, reveals cellular states and transitions at unprecedented resolution. Fully exploiting these data, however, requires kinetic models capable of unveiling governing regulatory functions. Here, we introduce an analytical framework dynamo (https://github.com/aristoteleo/dynamo-release), which infers absolute RNA velocity, reconstructs continuous vector fields that predict cell fates, employs differential geometry to extract underlying regulations, and ultimately predicts optimal reprogramming paths and perturbation outcomes. We highlight dynamo's power to overcome fundamental limitations of conventional splicing-based RNA velocity analyses to enable accurate velocity estimations on a metabolically labeled human hematopoiesis scRNA-seq dataset. Furthermore, differential geometry analyses reveal mechanisms driving early megakaryocyte appearance and elucidate asymmetrical regulation within the PU.1-GATA1 circuit. Leveraging the least-action-path method, dynamo accurately predicts drivers of numerous hematopoietic transitions. Finally, in silico perturbations predict cell-fate diversions induced by gene perturbations. Dynamo, thus, represents an important step in advancing quantitative and predictive theories of cell-state transitions.
Subject(s)
Single-Cell Analysis , Transcriptome/genetics , Algorithms , Female , Gene Expression Regulation , HL-60 Cells , Hematopoiesis/genetics , Hematopoietic Stem Cells/metabolism , Humans , Kinetics , Models, Biological , RNA, Messenger/metabolism , Staining and LabelingABSTRACT
Ferroptosis is a form of programmed cell death that is pathogenic to several acute and chronic diseases and executed via oxygenation of polyunsaturated phosphatidylethanolamines (PE) by 15-lipoxygenases (15-LO) that normally use free polyunsaturated fatty acids as substrates. Mechanisms of the altered 15-LO substrate specificity are enigmatic. We sought a common ferroptosis regulator for 15LO. We discovered that PEBP1, a scaffold protein inhibitor of protein kinase cascades, complexes with two 15LO isoforms, 15LO1 and 15LO2, and changes their substrate competence to generate hydroperoxy-PE. Inadequate reduction of hydroperoxy-PE due to insufficiency or dysfunction of a selenoperoxidase, GPX4, leads to ferroptosis. We demonstrated the importance of PEBP1-dependent regulatory mechanisms of ferroptotic death in airway epithelial cells in asthma, kidney epithelial cells in renal failure, and cortical and hippocampal neurons in brain trauma. As master regulators of ferroptotic cell death with profound implications for human disease, PEBP1/15LO complexes represent a new target for drug discovery.
Subject(s)
Acute Kidney Injury/pathology , Asthma/pathology , Brain Injuries, Traumatic/pathology , Cell Death , Phosphatidylethanolamine Binding Protein/metabolism , Acute Kidney Injury/metabolism , Animals , Apoptosis , Asthma/metabolism , Brain Injuries, Traumatic/metabolism , Cell Death/drug effects , Cell Line , Humans , Isoenzymes/metabolism , Lipoxygenase/chemistry , Lipoxygenase/metabolism , Mice , Models, Molecular , Oxazolidinones/pharmacology , Oxidation-Reduction , Phosphatidylethanolamine Binding Protein/chemistryABSTRACT
Chemokine-like receptor 1 (CMKLR1), also known as chemerin receptor 23 (ChemR23) or chemerin receptor 1, is a chemoattractant G protein-coupled receptor (GPCR) that responds to the adipokine chemerin and is highly expressed in innate immune cells, including macrophages and neutrophils. The signaling pathways of CMKLR1 can lead to both pro- and anti-inflammatory effects depending on the ligands and physiological contexts. To understand the molecular mechanisms of CMKLR1 signaling, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of the CMKLR1-Gi signaling complex with chemerin9, a nanopeptide agonist derived from chemerin, which induced complex phenotypic changes of macrophages in our assays. The cryo-EM structure, together with molecular dynamics simulations and mutagenesis studies, revealed the molecular basis of CMKLR1 signaling by elucidating the interactions at the ligand-binding pocket and the agonist-induced conformational changes. Our results are expected to facilitate the development of small molecule CMKLR1 agonists that mimic the action of chemerin9 to promote the resolution of inflammation.
Subject(s)
Intercellular Signaling Peptides and Proteins , Signal Transduction , Cryoelectron Microscopy , Receptors, G-Protein-Coupled/physiology , Chemokines/physiologyABSTRACT
Programmed ferroptotic death eliminates cells in all major organs and tissues with imbalanced redox metabolism due to overwhelming iron-catalyzed lipid peroxidation under insufficient control by thiols (Glutathione (GSH)). Ferroptosis has been associated with the pathogenesis of major chronic degenerative diseases and acute injuries of the brain, cardiovascular system, liver, kidneys, and other organs, and its manipulation offers a promising new strategy for anticancer therapy. This explains the high interest in designing new small-molecule-specific inhibitors against ferroptosis. Given the role of 15-lipoxygenase (15LOX) association with phosphatidylethanolamine (PE)-binding protein 1 (PEBP1) in initiating ferroptosis-specific peroxidation of polyunsaturated PE, we propose a strategy of discovering antiferroptotic agents as inhibitors of the 15LOX/PEBP1 catalytic complex rather than 15LOX alone. Here we designed, synthesized, and tested a customized library of 26 compounds using biochemical, molecular, and cell biology models along with redox lipidomic and computational analyses. We selected two lead compounds, FerroLOXIN-1 and 2, which effectively suppressed ferroptosis in vitro and in vivo without affecting the biosynthesis of pro-/anti-inflammatory lipid mediators in vivo. The effectiveness of these lead compounds is not due to radical scavenging or iron-chelation but results from their specific mechanisms of interaction with the 15LOX-2/PEBP1 complex, which either alters the binding pose of the substrate [eicosatetraenoyl-PE (ETE-PE)] in a nonproductive way or blocks the predominant oxygen channel thus preventing the catalysis of ETE-PE peroxidation. Our successful strategy may be adapted to the design of additional chemical libraries to reveal new ferroptosis-targeting therapeutic modalities.
Subject(s)
Ferroptosis , Phosphatidylethanolamine Binding Protein , Glutathione/metabolism , Iron/metabolism , Lipid Peroxidation , Lipids , Oxidation-Reduction , Phosphatidylethanolamine Binding Protein/antagonists & inhibitorsABSTRACT
SUMMARY: We introduce WatFinder, a tool designed to identify and visualize protein-water interactions (water bridges, water-mediated associations, or water channels, fluxes, and clusters) relevant to protein stability, dynamics, and function. WatFinder is integrated into ProDy, a Python API broadly used for structure-based prediction of protein dynamics. WatFinder provides a suite of functions for generating raw data as well as outputs from statistical analyses. The ProDy framework facilitates comprehensive automation and efficient analysis of the ensembles of structures resolved for a given protein or the time-evolved conformations from simulations in explicit water, as illustrated in five case studies presented in the Supplementary Material. AVAILABILITY AND IMPLEMENTATION: ProDy is open-source and freely available under MIT License from https://github.com/ProDy/ProDy.
Subject(s)
Proteins , Software , Water , Proteins/chemistry , Proteins/metabolism , Water/chemistry , Protein Conformation , Molecular Dynamics SimulationABSTRACT
Class B G protein-coupled receptors (GPCRs) are notoriously difficult to target by small molecules because their large orthosteric peptide-binding pocket embedded deep within the transmembrane domain limits the identification and development of nonpeptide small molecule ligands. Using the parathyroid hormone type 1 receptor (PTHR) as a prototypic class B GPCR target, and a combination of molecular dynamics simulations and elastic network model-based methods, we demonstrate that PTHR druggability can be effectively addressed. Here we found a key mechanical site that modulates the collective dynamics of the receptor and used this ensemble of PTHR conformers to identify selective small molecules with strong negative allosteric and biased properties for PTHR signaling in cell and PTH actions in vivo. This study provides a computational pipeline to detect precise druggable sites and identify allosteric modulators of PTHR signaling that could be extended to GPCRs to expedite discoveries of small molecules as novel therapeutic candidates.
Subject(s)
Receptor, Parathyroid Hormone, Type 1 , Receptors, G-Protein-Coupled , Ligands , Molecular Dynamics Simulation , Signal TransductionABSTRACT
The balance between NLRP3 inflammasome activation and mitophagy is essential for homeostasis and cellular health, but this relationship remains poorly understood. Here we found that interleukin-1α (IL-1α)-deficient macrophages have reduced caspase-1 activity and diminished IL-1ß release, concurrent with reduced mitochondrial damage, suggesting a role for IL-1α in regulating this balance. LPS priming of macrophages induced pro-IL-1α translocation to mitochondria, where it directly interacted with mitochondrial cardiolipin (CL). Computational modeling revealed a likely CL binding motif in pro-IL-1α, similar to that found in LC3b. Thus, binding of pro-IL-1α to CL in activated macrophages may interrupt CL-LC3b-dependent mitophagy, leading to enhanced Nlrp3 inflammasome activation and more robust IL-1ß production. Mutation of pro-IL-1α residues predicted to be involved in CL binding resulted in reduced pro-IL-1α-CL interaction, a reduction in NLRP3 inflammasome activity, and increased mitophagy. These data identify a function for pro-IL-1α in regulating mitophagy and the potency of NLRP3 inflammasome activation.
Subject(s)
Cardiolipins/metabolism , Interleukin-1alpha/metabolism , NLR Family, Pyrin Domain-Containing 3 Protein/metabolism , Animals , Autophagy , Cardiolipins/physiology , Caspase 1/metabolism , Female , HEK293 Cells , Humans , Inflammasomes/metabolism , Interleukin-1alpha/physiology , Macrophages/metabolism , Male , Mice , Mice, Knockout , Microtubule-Associated Proteins/metabolism , Mitochondria/metabolism , Mitophagy/physiology , NLR Family, Pyrin Domain-Containing 3 Protein/physiology , Protein Binding/physiology , Protein Domains/physiology , Reactive Oxygen Species/metabolismABSTRACT
The vast majority of membrane phospholipids (PLs) include two asymmetrically positioned fatty acyls: oxidizable polyunsaturated fatty acids (PUFA) attached predominantly at the sn2 position, and non-oxidizable saturated/monounsaturated acids (SFA/MUFA) localized at the sn1 position. The peroxidation of PUFA-PLs, particularly sn2-arachidonoyl(AA)- and sn2-adrenoyl(AdA)-containing phosphatidylethanolamines (PE), has been associated with the execution of ferroptosis, a program of regulated cell death. There is a minor subpopulation (≈1-2â mol %) of doubly PUFA-acylated phospholipids (di-PUFA-PLs) whose role in ferroptosis remains enigmatic. Here we report that 15-lipoxygenase (15LOX) exhibits unexpectedly high pro-ferroptotic peroxidation activity towards di-PUFA-PEs. We revealed that peroxidation of several molecular species of di-PUFA-PEs occurred early in ferroptosis. Ferrostatin-1, a typical ferroptosis inhibitor, effectively prevented peroxidation of di-PUFA-PEs. Furthermore, co-incubation of cells with di-AA-PE and 15LOX produced PUFA-PE peroxidation and induced ferroptotic death. The decreased contents of di-PUFA-PEs in ACSL4 KO A375 cells was associated with lower levels of di-PUFA-PE peroxidation and enhanced resistance to ferroptosis. Thus, di-PUFA-PE species are newly identified phospholipid peroxidation substrates and regulators of ferroptosis, representing a promising therapeutic target for many diseases related to ferroptotic death.
Subject(s)
Arachidonate 15-Lipoxygenase , Phosphatidylethanolamines , Phosphatidylethanolamines/metabolism , Arachidonate 15-Lipoxygenase/metabolism , Cell Death , Phospholipids/metabolism , Fatty Acids, Unsaturated/metabolism , Lipid PeroxidationABSTRACT
Ferroptosis, triggered by discoordination of iron, thiols and lipids, leads to the accumulation of 15-hydroperoxy (Hp)-arachidonoyl-phosphatidylethanolamine (15-HpETE-PE), generated by complexes of 15-lipoxygenase (15-LOX) and a scaffold protein, phosphatidylethanolamine (PE)-binding protein (PEBP)1. As the Ca2+-independent phospholipase A2ß (iPLA2ß, PLA2G6 or PNPLA9 gene) can preferentially hydrolyze peroxidized phospholipids, it may eliminate the ferroptotic 15-HpETE-PE death signal. Here, we demonstrate that by hydrolyzing 15-HpETE-PE, iPLA2ß averts ferroptosis, whereas its genetic or pharmacological inactivation sensitizes cells to ferroptosis. Given that PLA2G6 mutations relate to neurodegeneration, we examined fibroblasts from a patient with a Parkinson's disease (PD)-associated mutation (fPDR747W) and found selectively decreased 15-HpETE-PE-hydrolyzing activity, 15-HpETE-PE accumulation and elevated sensitivity to ferroptosis. CRISPR-Cas9-engineered Pnpla9R748W/R748W mice exhibited progressive parkinsonian motor deficits and 15-HpETE-PE accumulation. Elevated 15-HpETE-PE levels were also detected in midbrains of rotenone-infused parkinsonian rats and α-synuclein-mutant SncaA53T mice, with decreased iPLA2ß expression and a PD-relevant phenotype. Thus, iPLA2ß is a new ferroptosis regulator, and its mutations may be implicated in PD pathogenesis.
Subject(s)
Ferroptosis/physiology , Group VI Phospholipases A2/metabolism , Animals , Arachidonate 15-Lipoxygenase/metabolism , Disease Models, Animal , Female , Group VI Phospholipases A2/physiology , Humans , Iron/metabolism , Leukotrienes/metabolism , Lipid Metabolism/physiology , Lipid Peroxides/metabolism , Lipids/physiology , Male , Mice , Mice, Inbred C57BL , Oxidation-Reduction , Parkinson Disease/metabolism , Phosphatidylethanolamine Binding Protein/metabolism , Phospholipases/metabolism , Phospholipids/metabolism , Rats , Rats, Inbred LewABSTRACT
The high-level organization of the cell is embedded in indirect relationships that connect distinct cellular processes. Existing computational approaches for detecting indirect relationships between genes typically consist of propagating abstract information through network representations of the cell. However, the selection of genes to serve as the source of propagation is inherently biased by prior knowledge. Here, we sought to derive an unbiased view of the high-level organization of the cell by identifying the genes that propagate and receive information most effectively in the cell, and the indirect relationships between these genes. To this aim, we adapted a perturbation-response scanning strategy initially developed for identifying allosteric interactions within proteins. We deployed this strategy onto an elastic network model of the yeast genetic interaction profile similarity network. This network revealed a superior propensity for information propagation relative to simulated networks with similar topology. Perturbation-response scanning identified the major distributors and receivers of information in the network, named effector and sensor genes, respectively. Effectors formed dense clusters centrally integrated into the network, whereas sensors formed loosely connected antenna-shaped clusters and contained genes with previously characterized involvement in signal transduction. We propose that indirect relationships between effector and sensor clusters represent major paths of information flow between distinct cellular processes. Genetic similarity networks for fission yeast and human displayed similarly strong propensities for information propagation and clusters of effector and sensor genes, suggesting that the global architecture enabling indirect relationships is evolutionarily conserved across species. Our results demonstrate that elastic network modeling of cellular networks constitutes a promising strategy to probe the high-level organization and cooperativity in the cell.
Subject(s)
Gene Regulatory Networks , Proteins , Gene Regulatory Networks/genetics , Humans , Proteins/metabolism , Saccharomyces cerevisiae/genetics , Saccharomyces cerevisiae/metabolism , Signal Transduction/geneticsABSTRACT
Multisystem Inflammatory Syndrome in Children (MIS-C) associated with COVID-19 is a newly recognized condition in children with recent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. These children and adult patients with severe hyperinflammation present with a constellation of symptoms that strongly resemble toxic shock syndrome, an escalation of the cytotoxic adaptive immune response triggered upon the binding of pathogenic superantigens to T cell receptors (TCRs) and/or major histocompatibility complex class II (MHCII) molecules. Here, using structure-based computational models, we demonstrate that the SARS-CoV-2 spike (S) glycoprotein exhibits a high-affinity motif for binding TCRs, and may form a ternary complex with MHCII. The binding epitope on S harbors a sequence motif unique to SARS-CoV-2 (not present in other SARS-related coronaviruses), which is highly similar in both sequence and structure to the bacterial superantigen staphylococcal enterotoxin B. This interaction between the virus and human T cells could be strengthened by a rare mutation (D839Y/N/E) from a European strain of SARS-CoV-2. Furthermore, the interfacial region includes selected residues from an intercellular adhesion molecule (ICAM)-like motif shared between the SARS viruses from the 2003 and 2019 pandemics. A neurotoxin-like sequence motif on the receptor-binding domain also exhibits a high tendency to bind TCRs. Analysis of the TCR repertoire in adult COVID-19 patients demonstrates that those with severe hyperinflammatory disease exhibit TCR skewing consistent with superantigen activation. These data suggest that SARS-CoV-2 S may act as a superantigen to trigger the development of MIS-C as well as cytokine storm in adult COVID-19 patients, with important implications for the development of therapeutic approaches.
Subject(s)
Betacoronavirus/immunology , Coronavirus Infections/immunology , Pneumonia, Viral/immunology , Receptors, Antigen, T-Cell/metabolism , Spike Glycoprotein, Coronavirus/metabolism , Superantigens/metabolism , Systemic Inflammatory Response Syndrome/immunology , Amino Acid Motifs , Betacoronavirus/chemistry , Betacoronavirus/genetics , Betacoronavirus/metabolism , COVID-19 , Coronavirus Infections/genetics , Coronavirus Infections/pathology , Enterotoxins/chemistry , Epitopes, T-Lymphocyte , Humans , Intercellular Adhesion Molecule-1/chemistry , Models, Molecular , Mutation , Neurotoxins/chemistry , Pandemics , Pneumonia, Viral/genetics , Pneumonia, Viral/pathology , Protein Binding , Receptors, Antigen, T-Cell/chemistry , Receptors, Antigen, T-Cell/genetics , SARS-CoV-2 , Spike Glycoprotein, Coronavirus/chemistry , Spike Glycoprotein, Coronavirus/genetics , Superantigens/chemistry , Superantigens/genetics , Systemic Inflammatory Response Syndrome/genetics , Systemic Inflammatory Response Syndrome/pathologyABSTRACT
Temporally harmonized elimination of damaged or unnecessary organelles and cells is a prerequisite of health. Under Type 2 inflammatory conditions, human airway epithelial cells (HAECs) generate proferroptotic hydroperoxy-arachidonoyl-phosphatidylethanolamines (HpETE-PEs) as proximate death signals. Production of 15-HpETE-PE depends on activation of 15-lipoxygenase-1 (15LO1) in complex with PE-binding protein-1 (PEBP1). We hypothesized that cellular membrane damage induced by these proferroptotic phospholipids triggers compensatory prosurvival pathways, and in particular autophagic pathways, to prevent cell elimination through programmed death. We discovered that PEBP1 is pivotal to driving dynamic interactions with both proferroptotic 15LO1 and the autophagic protein microtubule-associated light chain-3 (LC3). Further, the 15LO1-PEBP1-generated ferroptotic phospholipid, 15-HpETE-PE, promoted LC3-I lipidation to stimulate autophagy. This concurrent activation of autophagy protects cells from ferroptotic death and release of mitochondrial DNA. Similar findings are observed in Type 2 Hi asthma, where high levels of both 15LO1-PEBP1 and LC3-II are seen in HAECs, in association with low bronchoalveolar lavage fluid mitochondrial DNA and more severe disease. The concomitant activation of ferroptosis and autophagy by 15LO1-PEBP1 complexes and their hydroperoxy-phospholipids reveals a pathobiologic pathway relevant to asthma and amenable to therapeutic targeting.
Subject(s)
Arachidonate 15-Lipoxygenase/metabolism , Asthma/immunology , Autophagy/immunology , Epithelial Cells/pathology , Ferroptosis/immunology , Phosphatidylethanolamine Binding Protein/metabolism , Adult , Animals , Asthma/diagnosis , Asthma/pathology , Bronchoalveolar Lavage Fluid/cytology , Cell Line , Cell Survival/immunology , Epithelial Cells/immunology , Female , Gene Knockout Techniques , Humans , Hydroxyeicosatetraenoic Acids/immunology , Hydroxyeicosatetraenoic Acids/metabolism , Interleukin-13/immunology , Interleukin-13/metabolism , Male , Mice , Microtubule-Associated Proteins/metabolism , Molecular Dynamics Simulation , Phosphatidylethanolamine Binding Protein/genetics , Phosphatidylethanolamines/immunology , Phosphatidylethanolamines/metabolism , Primary Cell Culture , Protein Binding/immunology , Severity of Illness IndexABSTRACT
The insertion or deletion (indel) of amino acids has a variety of effects on protein function, ranging from disease-forming changes to gaining new functions. Despite their importance, indels have not been systematically characterized towards protein engineering or modification goals. In the present work, we focus on deletions composed of multiple contiguous amino acids (mAA-dels) and their effects on the protein (mutant) folding ability. Our analysis reveals that the mutant retains the native fold when the mAA-del obeys well-defined structural dynamics properties: localization in intrinsically flexible regions, showing low resistance to mechanical stress, and separation from allosteric signaling paths. Motivated by the possibility of distinguishing the features that underlie the adaptability of proteins to mAA-dels, and by the rapid evaluation of these features using elastic network models, we developed a positive-unlabeled learning-based classifier that can be adopted for protein design purposes. Trained on a consolidated set of features, including those reflecting the intrinsic dynamics of the regions where the mAA-dels occur, the new classifier yields a high recall of 84.3% for identifying mAA-dels that are stably tolerated by the protein. The comparative examination of the relative contribution of different features to the prediction reveals the dominant role of structural dynamics in enabling the adaptation of the mutant to mAA-del without disrupting the native fold.
Subject(s)
Amino Acids , Proteins , Amino Acids/genetics , Proteins/chemistry , INDEL Mutation , Protein EngineeringABSTRACT
Aberrant dopamine (DA) signaling is associated with several psychiatric disorders, such as autism, bipolar disorder, addiction, and Parkinson's disease, and several medications that target the DA transporter (DAT) can induce or treat these disorders. In addition, psychostimulants, such as cocaine and D-amphetamine (AMPH), rely on the competitive interactions with the transporter's substrate binding site to produce their rewarding effects. Agents that exhibit noncompetitive, allosteric modulation of DAT remain an important topic of investigation due to their potential therapeutic applications. We previously identified a novel allosteric modulator of human DAT, KM822, that can decrease the affinity of cocaine for DAT and attenuate cocaine-elicited behaviors; however, whether DAT is the sole mediator of KM822 actions in vivo is unproven given the large number of potential off-target sites. Here, we provide in silico and in vitro evidence that the allosteric site engaged by KM822 is conserved between human DAT and Caenorhabditis elegans DAT-1. KM822 binds to a similar pocket in DAT-1 as previously identified in human DAT. In functional dopamine uptake assays, KM822 affects the interaction between AMPH and DAT-1 by reducing the affinity of AMPH for DAT-1. Finally, through a combination of genetic and pharmacological in vivo approaches we provide evidence that KM822 diminishes the behavioral actions of AMPH on swimming-induced paralysis through a direct allosteric modulation of DAT-1. More broadly, our findings demonstrate allosteric modulation of DAT as a behavior modifying strategy and suggests that Caenorhabditis elegans can be operationalized to identify and investigate the interactions of DAT allosteric modulators. SIGNIFICANCE STATEMENT: We previously demonstrated that the dopamine transporter (DAT) allosteric modulator KM822 decreases cocaine affinity for human DAT. Here, using in silico and in vivo genetic approaches, we extend this finding to interactions with amphetamine, demonstrating evolutionary conservation of the DAT allosteric site. In Caenorhabditis elegans, we report that KM822 suppresses amphetamine behavioral effects via specific interactions with DAT-1. Our findings reveal Caenorhabditis elegans as a new tool to study allosteric modulation of DAT and its behavioral consequences.
Subject(s)
Amphetamine/metabolism , Caenorhabditis elegans Proteins/metabolism , Dopamine Agents/metabolism , Dopamine Plasma Membrane Transport Proteins/metabolism , Allosteric Regulation/drug effects , Allosteric Regulation/physiology , Amphetamine/pharmacology , Animals , COS Cells , Caenorhabditis elegans , Caenorhabditis elegans Proteins/chemistry , Chlorocebus aethiops , Dopamine Agents/pharmacology , Dopamine Plasma Membrane Transport Proteins/chemistry , Dose-Response Relationship, Drug , Drosophila melanogaster , Protein Binding/drug effects , Protein Binding/physiology , Protein Structure, SecondaryABSTRACT
Dopamine transporter (DAT) mediates the reuptake of synaptically released dopamine, and thus controls the duration and intensity of dopamine neurotransmission. Mammalian DAT has been observed to form oligomers, although the mechanisms of oligomerization and its role in DAT activity and trafficking remain largely unknown. We discovered a series of small molecule compounds that stabilize trimers and induce high-order oligomers of DAT and concomitantly promote its clathrin-independent endocytosis. Using a combination of chemical cross-linking, fluorescence resonance energy transfer microscopy, antibody-uptake endocytosis assay, live-cell lattice light sheet microscopy, ligand binding and substrate transport kinetics analyses, and molecular modeling and simulations, we investigated molecular basis of DAT oligomerization and endocytosis induced by these compounds. Our study showed that small molecule-induced DAT oligomerization and endocytosis are favored by the inward-facing DAT conformation and involve interactions of four hydrophobic residues at the interface between transmembrane (TM) helices TM4 and TM9. Surprisingly, a corresponding quadruple DAT mutant displays altered dopamine transport kinetics and increased cocaine-analog binding. The latter is shown to originate from an increased preference for outward-facing conformation and inward-to-outward transition. Taken together, our results demonstrate a direct coupling between conformational dynamics of DAT, functional activity of the transporter, and its oligomerization leading to endocytosis. The high specificity of such coupling for DAT makes the TM4-9 hub a new target for pharmacological modulation of DAT activity and subcellular localization.
Subject(s)
Dopamine Plasma Membrane Transport Proteins/drug effects , Dopamine Plasma Membrane Transport Proteins/metabolism , Dopamine/metabolism , Animals , Cell Line , Clathrin/metabolism , Dopamine Plasma Membrane Transport Proteins/physiology , Endocytosis/drug effects , Endocytosis/physiology , Endothelial Cells/metabolism , Fluorescence Resonance Energy Transfer/methods , Humans , Models, Molecular , Protein Binding , Protein Conformation , Small Molecule Libraries/pharmacology , SwineABSTRACT
SUMMARY: Efficient sampling of conformational space is essential for elucidating functional/allosteric mechanisms of proteins and generating ensembles of conformers for docking applications. However, unbiased sampling is still a challenge especially for highly flexible and/or large systems. To address this challenge, we describe a new implementation of our computationally efficient algorithm ClustENMD that is integrated with ProDy and OpenMM softwares. This hybrid method performs iterative cycles of conformer generation using elastic network model for deformations along global modes, followed by clustering and short molecular dynamics simulations. ProDy framework enables full automation and analysis of generated conformers and visualization of their distributions in the essential subspace. AVAILABILITY AND IMPLEMENTATION: ClustENMD is open-source and freely available under MIT License from https://github.com/prody/ProDy. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Proteins , Software , Proteins/metabolism , Protein Conformation , Algorithms , Molecular Dynamics SimulationABSTRACT
SUMMARY: ProDy, an integrated application programming interface developed for modelling and analysing protein dynamics, has significantly evolved in recent years in response to the growing data and needs of the computational biology community. We present major developments that led to ProDy 2.0: (i) improved interfacing with databases and parsing new file formats, (ii) SignDy for signature dynamics of protein families, (iii) CryoDy for collective dynamics of supramolecular systems using cryo-EM density maps and (iv) essential site scanning analysis for identifying sites essential to modulating global dynamics. AVAILABILITY AND IMPLEMENTATION: ProDy is open-source and freely available under MIT License from https://github.com/prody/ProDy. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
ABSTRACT
Understanding the mechanism of SARS-CoV-2 infection and identifying potential therapeutics are global imperatives. Using a quantitative systems pharmacology approach, we identified a set of repurposable and investigational drugs as potential therapeutics against COVID-19. These were deduced from the gene expression signature of SARS-CoV-2-infected A549 cells screened against Connectivity Map and prioritized by network proximity analysis with respect to disease modules in the viral-host interactome. We also identified immuno-modulating compounds aiming at suppressing hyperinflammatory responses in severe COVID-19 patients, based on the transcriptome of ACE2-overexpressing A549 cells. Experiments with Vero-E6 cells infected by SARS-CoV-2, as well as independent syncytia formation assays for probing ACE2/SARS-CoV-2 spike protein-mediated cell fusion using HEK293T and Calu-3 cells, showed that several predicted compounds had inhibitory activities. Among them, salmeterol, rottlerin, and mTOR inhibitors exhibited antiviral activities in Vero-E6 cells; imipramine, linsitinib, hexylresorcinol, ezetimibe, and brompheniramine impaired viral entry. These novel findings provide new paths for broadening the repertoire of compounds pursued as therapeutics against COVID-19.
Subject(s)
Antiviral Agents/pharmacology , COVID-19 Drug Treatment , Drug Evaluation, Preclinical/methods , Virus Internalization/drug effects , Angiotensin-Converting Enzyme 2/genetics , Angiotensin-Converting Enzyme 2/metabolism , Animals , Anti-Inflammatory Agents, Non-Steroidal/pharmacology , COVID-19/genetics , COVID-19/virology , Chlorocebus aethiops , Drug Repositioning , HEK293 Cells , Host-Pathogen Interactions/drug effects , Host-Pathogen Interactions/physiology , Humans , Imidazoles/pharmacology , Pyrazines/pharmacology , SARS-CoV-2/drug effects , SARS-CoV-2/pathogenicity , Salmeterol Xinafoate/pharmacology , Vero CellsABSTRACT
Peptide ligands of class B G-protein-coupled receptors act via a two-step binding process, but the essential mechanisms that link their extracellular binding to intracellular receptor-arrestin interactions are not fully understood. Using NMR, crosslinking coupled to mass spectrometry, signaling experiments and computational approaches on the parathyroid hormone (PTH) type 1 receptor (PTHR), we show that initial binding of the PTH C-terminal part constrains the conformation of the flexible PTH N-terminal signaling epitope before a second binding event occurs. A 'hot-spot' PTH residue, His9, that inserts into the PTHR transmembrane domain at this second step allosterically engages receptor-arrestin coupling. A conformational change in PTHR intracellular loop 3 permits favorable interactions with ß-arrestin's finger loop. These results unveil structural determinants for PTHR-arrestin complex formation and reveal that the two-step binding mechanism proceeds via cooperative fluctuations between ligand and receptor, which extend to other class B G-protein-coupled receptors.
Subject(s)
Arrestin/metabolism , Parathyroid Hormone/metabolism , Arrestin/chemistry , Calcium Phosphates , Cryoelectron Microscopy , Cyclic AMP , Escherichia coli , HEK293 Cells , Humans , Molecular Dynamics Simulation , Parathyroid Hormone/chemistry , Receptors, G-Protein-CoupledABSTRACT
Ferroptotic death is the penalty for losing control over three processes-iron metabolism, lipid peroxidation and thiol regulation-that are common in the pro-inflammatory environment where professional phagocytes fulfill their functions and yet survive. We hypothesized that redox reprogramming of 15-lipoxygenase (15-LOX) during the generation of pro-ferroptotic signal 15-hydroperoxy-eicosa-tetra-enoyl-phosphatidylethanolamine (15-HpETE-PE) modulates ferroptotic endurance. Here, we have discovered that inducible nitric oxide synthase (iNOS)/NOâ¢-enrichment of activated M1 (but not alternatively activated M2) macrophages/microglia modulates susceptibility to ferroptosis. Genetic or pharmacologic depletion/inactivation of iNOS confers sensitivity on M1 cells, whereas NO⢠donors empower resistance of M2 cells to ferroptosis. In vivo, M1 phagocytes, in comparison to M2 phagocytes, exert higher resistance to pharmacologically induced ferroptosis. This resistance is diminished in iNOS-deficient cells in the pro-inflammatory conditions of brain trauma or the tumour microenvironment. The nitroxygenation of eicosatetraenoyl (ETE)-PE intermediates and oxidatively truncated species by NO⢠donors and/or suppression of NO⢠production by iNOS inhibitors represent a novel redox mechanism of regulation of ferroptosis in pro-inflammatory conditions.