Cancer Discovery (2021), accepted
Cell (2019), vol. 178(4), 949-963
The target of many drugs are established early in their development and are not systematically assessed in different biological models. Thus, it is often unclear whether therapeutics with the same nominal targets but different chemical structures are functionally equivalent and induce similar biological effects. To address these questions, I analyzed data from different phenotypic and biochemical assays to compare therapeutics.
At Genentech, we characterized the role of degradation in estrogen therapies and identify chromatin open and transcription factor mobility to be the main determinants of response. In another project, we identified the role of protein degradation in the mechanism of action of the latest generation of PI3Ka inhibitor.
During my postdoc, we focused on cyclin dependent kinases 4/6. We find that transcriptional, proteomic, and phenotypic changes induced by palbociclib, ribociclib, and abemaciclib differ significantly; abemaciclib in particular has advantageous activities partially overlapping those of alvocidib, an older polyselective CDK inhibitor. I also analyzed the transcriptional response of six cell lines to more than 100 different kinase inhibitors. By developing a new clustering approach, I found that the transcriptional signature of RTK inhibitors clusters with either inhibitors of the PI3K/AKT or the MAPK pathway depending on the cell line. Based on these signatures, I inferred the changes in transcription factor activity and identified synergetic drug combinations that we validated.
Measurements of dose-response curves are fundamental to investigate the efficacy of therapeutic drugs and identify genes involved in sensitivity and resistance. However, one source of variability that strongly biases the traditional metrics of drug response, IC50 and Emax, and could potentially explain variations between studies is growth rate. To address this problem, I developed a new method to quantify drug response: the normalized growth rate inhibition (GR), from which I derived two metrics (GR50 andGRmax), that are independent of the duration of the assay and the division rate of the cell line in contrast to IC50 and Emax (see figure). As patient-derived cell lines generally grow slowly, the problem is exacerbated and thus the GR metrics are essential to properly compare results across established and patient-derived cell lines, which will be crucial for my research. Because the GR values are independent of the length of the assay, we can adjust the duration of the experiment to accommodate variations in growth rate between cell lines. We expanded on the GR framework to make it possible to distinguish cytotoxic and cytostatic drug effects based on Deep Dye Drop protocol.
I am now combining the GR with the Loewe additivity model to quantify and compare synergies across multiple cell lines. This new synergy metric will allow me to quantify differences in phenotype where, for example, two agents are cytostatic individually and cytotoxic when combined – something that other metrics fail to quantify.
Most studies on drug sensitivity focus on genomic features. Here we investigated how measurements of the basal and perturbed states of signal transduction networks relate to drug response in breast cancer cell lines. We found that growth factor signaling predicts sensitivity to RTK and PI3K inhibitors. In particular, response to heregulin (HRG) is more predictive of sensitivity to inhibitors of the PI3K/AKT pathway than AKT activity. Further analyses showed that HRG response depends on the levels of ErbB2 and ErbB3. In this work, I developed new statistical approaches to overcome the problems of multiple hypothesis testing which will be essential for my future research. My results led to a better understanding of the relation between RTK signaling and kinase inhibitor sensitivity.
The partial response of tumor cells to TRAIL and therapeutic antibody agonists is limiting the efficacy of this class of death-inducing agent as anti‐cancer drugs. Previous studies in the Sorger Lab raised the question of the role of non‐genetic cell‐to‐cell variability in incomplete killing by TRAIL. Based on single cell live microscopy data reporting caspase-8 activity and cell death in response to TRAIL treatment, I developed a model of initiator caspase dynamics. This model allowed me to identify a threshold in caspase activity that separates dying and surviving cells. Together with Dr. J. Roux, an experimental postdoctoral fellow, we used the model to predict how co-treatment with other therapeutic agents impacts on caspase-8 activity and cell death. We also identified treatment regimens that overcome overexpression of oncogenes FLIP-S/-L and BCL-2/-XL.
During my PhD, I developed novel computational methods to quantify the robustness of biological networks. Robustness is a fundamental property of many biological systems and accounting for it in biological modeling is a challenge because it requires quantitative metrics capable of comparing robustness across models. By formalizing the evaluation of robustness, I was able to use it to discriminate between models and parameters, and in the absence of other criteria, decide that a model, or a parameter set, would be judged superior if it is more robust than other models to some class of perturbations. I applied my methods to different biological systems that are notionally robust or designed to be robust. For example, I identified causes of robustness in models of the cyanobacterial circadian clock, found the best network architecture to obtain stable oscillations in a model of the Suprachiasmatic nucleus, and designed models robust to external perturbations for synthetic circuits for pancreatic tissue homeostasis.