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,729
    views
  • 783
    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
    Christopher TA Lewis, Elise G Melhedegaard ... Julien Ochala
    Research Article

    Hibernation is a period of metabolic suppression utilized by many small and large mammal species to survive during winter periods. As the underlying cellular and molecular mechanisms remain incompletely understood, our study aimed to determine whether skeletal muscle myosin and its metabolic efficiency undergo alterations during hibernation to optimize energy utilization. We isolated muscle fibers from small hibernators, Ictidomys tridecemlineatus and Eliomys quercinus and larger hibernators, Ursus arctos and Ursus americanus. We then conducted loaded Mant-ATP chase experiments alongside X-ray diffraction to measure resting myosin dynamics and its ATP demand. In parallel, we performed multiple proteomics analyses. Our results showed a preservation of myosin structure in U. arctos and U. americanus during hibernation, whilst in I. tridecemlineatus and E. quercinus, changes in myosin metabolic states during torpor unexpectedly led to higher levels in energy expenditure of type II, fast-twitch muscle fibers at ambient lab temperatures (20 °C). Upon repeating loaded Mant-ATP chase experiments at 8 °C (near the body temperature of torpid animals), we found that myosin ATP consumption in type II muscle fibers was reduced by 77–107% during torpor compared to active periods. Additionally, we observed Myh2 hyper-phosphorylation during torpor in I. tridecemilineatus, which was predicted to stabilize the myosin molecule. This may act as a potential molecular mechanism mitigating myosin-associated increases in skeletal muscle energy expenditure during periods of torpor in response to cold exposure. Altogether, we demonstrate that resting myosin is altered in hibernating mammals, contributing to significant changes to the ATP consumption of skeletal muscle. Additionally, we observe that it is further altered in response to cold exposure and highlight myosin as a potentially contributor to skeletal muscle non-shivering thermogenesis.

    1. Cell Biology
    Jun Yang, Shitian Zou ... Xiaochun Bai
    Research Article

    Quiescence (G0) maintenance and exit are crucial for tissue homeostasis and regeneration in mammals. Here, we show that methyl-CpG binding protein 2 (Mecp2) expression is cell cycle-dependent and negatively regulates quiescence exit in cultured cells and in an injury-induced liver regeneration mouse model. Specifically, acute reduction of Mecp2 is required for efficient quiescence exit as deletion of Mecp2 accelerates, while overexpression of Mecp2 delays quiescence exit, and forced expression of Mecp2 after Mecp2 conditional knockout rescues cell cycle reentry. The E3 ligase Nedd4 mediates the ubiquitination and degradation of Mecp2, and thus facilitates quiescence exit. A genome-wide study uncovered the dual role of Mecp2 in preventing quiescence exit by transcriptionally activating metabolic genes while repressing proliferation-associated genes. Particularly disruption of two nuclear receptors, Rara or Nr1h3, accelerates quiescence exit, mimicking the Mecp2 depletion phenotype. Our studies unravel a previously unrecognized role for Mecp2 as an essential regulator of quiescence exit and tissue regeneration.