The transcriptomic response of cells to a drug combination is more than the sum of the responses to the monotherapies

  1. Jennifer EL Diaz
  2. Mehmet Eren Ahsen
  3. Thomas Schaffter
  4. Xintong Chen
  5. Ronald B Realubit
  6. Charles Karan
  7. Andrea Califano
  8. Bojan Losic
  9. Gustavo Stolovitzky  Is a corresponding author
  1. Icahn School of Medicine at Mount Sinai, United States
  2. IBM Research, United States
  3. Columbia University, United States

Abstract

Our ability to discover effective drug combinations is limited, in part by insufficient understanding of how the transcriptional response of two monotherapies results in that of their combination. We analyzed matched time course RNAseq profiling of cells treated with single drugs and their combinations and found that the transcriptional signature of the synergistic combination was unique relative to that of either constituent monotherapy. The sequential activation of transcription factors in time in the gene regulatory network was implicated. The nature of this transcriptional cascade suggests that drug synergy may ensue when the transcriptional responses elicited by two unrelated individual drugs are correlated. We used these results as the basis of a simple prediction algorithm attaining an AUROC of 0.77 in the prediction of synergistic drug combinations in an independent dataset.

Data availability

Raw RNAseq data have been deposited in GEO under accession code GSE149428. Code is available at github.com/jennifereldiaz/drug-synergy.

The following previously published data sets were used

Article and author information

Author details

  1. Jennifer EL Diaz

    Department of Genetics and Genomics Sciences, Icahn School of Medicine at Mount Sinai, New York, United States
    Competing interests
    No competing interests declared.
  2. Mehmet Eren Ahsen

    Department of Genetics and Genomics Sciences, Icahn School of Medicine at Mount Sinai, New York, United States
    Competing interests
    No competing interests declared.
  3. Thomas Schaffter

    IBM Computational Biology Center, IBM Research, Yorktown Heights, United States
    Competing interests
    No competing interests declared.
  4. Xintong Chen

    Department of Genetics and Genomics Sciences, Icahn School of Medicine at Mount Sinai, New York, United States
    Competing interests
    No competing interests declared.
  5. Ronald B Realubit

    Department of Systems Biology, Sulzberger Columbia Genome Center, High Throughput Screening Facility, Columbia University, New York, United States
    Competing interests
    No competing interests declared.
  6. Charles Karan

    Department of Systems Biology, Sulzberger Columbia Genome Center, High Throughput Screening Facility, Columbia University, New York, United States
    Competing interests
    No competing interests declared.
  7. Andrea Califano

    Department of Systems Biology, Department of Biomedical Informatics, Department of Biochemistry and Molecular Biophysics, J.P. Sulzberger Columbia Genome Center, Herbert Irving Comprehensive Cancer, Columbia University, New York, United States
    Competing interests
    Andrea Califano, Dr. Califano is founder, equity holder, consultant, and director of DarwinHealth Inc., a company that has licensed some of the algorithms used in this manuscript from Columbia University. Columbia University is also an equity holder in DarwinHealth Inc..
  8. Bojan Losic

    Icahn Institute for Genomics and Multiscale Biology, Icahn School of Medicine at Mount Sinai, New York, United States
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3643-4612
  9. Gustavo Stolovitzky

    IBM Computational Biology Center, IBM Research, Yorktown Heights, United States
    For correspondence
    gustavo@us.ibm.com
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9618-2819

Funding

Internal funding from IBM Research and the Icahn School of Medicine at Mount Sinai to GS

  • Jennifer EL Diaz
  • Mehmet Eren Ahsen
  • Thomas Schaffter
  • Xintong Chen
  • Bojan Losic
  • Gustavo Stolovitzky

National Institutes of Health (NIH T32 GM007280)

  • Jennifer EL Diaz

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Reviewing Editor

  1. Alfonso Valencia, Barcelona Supercomputing Center - BSC, Spain

Version history

  1. Received: October 13, 2019
  2. Accepted: August 17, 2020
  3. Accepted Manuscript published: September 18, 2020 (version 1)
  4. Version of Record published: October 9, 2020 (version 2)

Copyright

© 2020, Diaz et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 7,940
    views
  • 810
    downloads
  • 19
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Jennifer EL Diaz
  2. Mehmet Eren Ahsen
  3. Thomas Schaffter
  4. Xintong Chen
  5. Ronald B Realubit
  6. Charles Karan
  7. Andrea Califano
  8. Bojan Losic
  9. Gustavo Stolovitzky
(2020)
The transcriptomic response of cells to a drug combination is more than the sum of the responses to the monotherapies
eLife 9:e52707.
https://doi.org/10.7554/eLife.52707

Share this article

https://doi.org/10.7554/eLife.52707

Further reading

    1. Biochemistry and Chemical Biology
    2. Cell Biology
    Ya-Juan Wang, Xiao-Jing Di ... Ting-Wei Mu
    Research Article Updated

    Protein homeostasis (proteostasis) deficiency is an important contributing factor to neurological and metabolic diseases. However, how the proteostasis network orchestrates the folding and assembly of multi-subunit membrane proteins is poorly understood. Previous proteomics studies identified Hsp47 (Gene: SERPINH1), a heat shock protein in the endoplasmic reticulum lumen, as the most enriched interacting chaperone for gamma-aminobutyric acid type A (GABAA) receptors. Here, we show that Hsp47 enhances the functional surface expression of GABAA receptors in rat neurons and human HEK293T cells. Furthermore, molecular mechanism study demonstrates that Hsp47 acts after BiP (Gene: HSPA5) and preferentially binds the folded conformation of GABAA receptors without inducing the unfolded protein response in HEK293T cells. Therefore, Hsp47 promotes the subunit-subunit interaction, the receptor assembly process, and the anterograde trafficking of GABAA receptors. Overexpressing Hsp47 is sufficient to correct the surface expression and function of epilepsy-associated GABAA receptor variants in HEK293T cells. Hsp47 also promotes the surface trafficking of other Cys-loop receptors, including nicotinic acetylcholine receptors and serotonin type 3 receptors in HEK293T cells. Therefore, in addition to its known function as a collagen chaperone, this work establishes that Hsp47 plays a critical and general role in the maturation of multi-subunit Cys-loop neuroreceptors.

    1. Cell Biology
    2. Immunology and Inflammation
    Kevin Portmann, Aline Linder, Klaus Eyer
    Research Article

    Cytokine polyfunctionality is a well-established concept in immune cells, especially T cells, and their ability to concurrently produce multiple cytokines has been associated with better immunological disease control and subsequent effectiveness during infection and disease. To date, only little is known about the secretion dynamics of those cells, masked by the widespread deployment of mainly time-integrated endpoint measurement techniques that do not easily differentiate between concurrent and sequential secretion. Here, we employed a single-cell microfluidic platform capable of resolving the secretion dynamics of individual PBMCs. To study the dynamics of poly-cytokine secretion, as well as the dynamics of concurrent and sequential polyfunctionality, we analyzed the response at different time points after ex vivo activation. First, we observed the simultaneous secretion of cytokines over the measurement time for most stimulants in a subpopulation of cells only. Second, polyfunctionality generally decreased with prolonged stimulation times and revealed no correlation with the concentration of secreted cytokines in response to stimulation. However, we observed a general trend towards higher cytokine secretion in polyfunctional cells, with their secretion dynamics being distinctly different from mono-cytokine-secreting cells. This study provided insights into the distinct secretion behavior of heterogenous cell populations after stimulation with well-described agents and such a system could provide a better understanding of various immune dynamics in therapy and disease.