A mechanistic pan-cancer pathway model informed by multi-omics data interprets stochastic cell fate responses to drugs and mitogens.

Most cancer cells harbor multiple drivers whose epistasis and interactions with expression context clouds drug and drug combination sensitivity prediction. We constructed a mechanistic computational model that is context-tailored by omics data to capture regulation of stochastic proliferation and death by pan-cancer driver pathways. Simulations and experiments explore how the coordinated dynamics of RAF/MEK/ERK and PI-3K/AKT kinase activities in response to synergistic mitogen or drug combinations control cell fate in a specific cellular context. In this MCF10A cell context, simulations suggest that synergistic ERK and AKT inhibitor-induced death is likely mediated by BIM rather than BAD, which is supported by prior experimental studies. AKT dynamics explain S-phase entry synergy between EGF and insulin, but simulations suggest that stochastic ERK, and not AKT, dynamics seem to drive cell-to-cell proliferation variability, which in simulations is predictable from pre-stimulus fluctuations in C-Raf/B-Raf levels. Simulations suggest MEK alteration negligibly influences transformation, consistent with clinical data. Tailoring the model to an alternate cell expression and mutation context, a glioma cell line, allows prediction of increased sensitivity of cell death to AKT inhibition. Our model mechanistically interprets context-specific landscapes between driver pathways and cell fates, providing a framework for designing more rational cancer combination therapy.

Antiinflammatory Effects of Budesonide in Human Fetal Lung.

Lung inflammation in premature infants contributes to the development of bronchopulmonary dysplasia (BPD), a chronic lung disease with long-term sequelae. Pilot studies administering budesonide suspended in surfactant have found reduced BPD without the apparent adverse effects that occur with systemic dexamethasone therapy. Our objective was to determine budesonide potency, stability, and antiinflammatory effects in human fetal lung. We cultured explants of second-trimester fetal lung with budesonide or dexamethasone and used microscopy, immunoassays, RNA sequencing, liquid chromatography/tandem mass spectrometry, and pulsating bubble surfactometry. Budesonide suppressed secreted chemokines IL-8 and CCL2 (MCP-1) within 4 hours, reaching a 90% decrease at 12 hours, which was fully reversed 72 hours after removal of the steroid. Half-maximal effects occurred at 0.04-0.05 nM, representing a fivefold greater potency than for dexamethasone. Budesonide significantly induced 3.6% and repressed 2.8% of 14,500 sequenced mRNAs by 1.6- to 95-fold, including 119 genes that contribute to the glucocorticoid inflammatory transcriptome; some are known targets of nuclear factor-κB. By global proteomics, 22 secreted inflammatory proteins were hormonally regulated. Two glucocorticoid-regulated genes of interest because of their association with lung disease are CHI3L1 and IL1RL1. Budesonide retained activity in the presence of surfactant and did not alter its surface properties. There was some formation of palmitate-budesonide in lung tissue but no detectable metabolism to inactive 16α-hydroxy prednisolone. We concluded that budesonide is a potent and stable antiinflammatory glucocorticoid in human fetal lung in vitro, supporting a beneficial antiinflammatory response to lung-targeted budesonide:surfactant treatment of infants for the prevention of BPD.

Effects of Advancing Gestation and Non-Caucasian Race on Ductus Arteriosus Gene Expression.

OBJECTIVE: To identify genes affected by advancing gestation and racial/ethnic origin in human ductus arteriosus (DA). STUDY DESIGN: We collected 3 sets of DA tissue (n = 93, n = 89, n = 91; total = 273 fetuses) from second trimester pregnancies. We examined four genes, with DNA polymorphisms that distribute along racial lines, to identify "Caucasian" and "non-Caucasian" DA. We used real time polymerase chain reaction to measure RNA expression of 48 candidate genes involved in functional closure of the DA, and used multivariable regression analyses to examine the relationships between advancing gestation, "non-Caucasian" race, and gene expression. RESULTS: Mature gestation and non-Caucasian race are significant predictors for identifying infants who will close their patent DA when treated with indomethacin. Advancing gestation consistently altered gene expression in pathways involved with oxygen-induced constriction (eg, calcium-channels, potassium-channels, and endothelin signaling), contractile protein maturation, tissue remodeling, and prostaglandin and nitric oxide signaling in all 3 tissue sets. None of the pathways involved with oxygen-induced constriction appeared to be altered in "non-Caucasian" DA. Two genes, SLCO2A1 and NOS3, (involved with prostaglandin reuptake/metabolism and nitric oxide production, respectively) were consistently decreased in "non-Caucasian" DA. CONCLUSIONS: Prostaglandins and nitric oxide are the most important vasodilators opposing DA closure. Indomethacin inhibits prostaglandin production, but not nitric oxide production. Because decreased SLCO2A1 and NOS3 expression can lead to increased prostaglandin and decreased nitric oxide concentrations, we speculate that prostaglandin-mediated vasodilation may play a more dominant role in maintaining the "non-Caucasian" patent DA, making it more likely to close when inhibited by indomethacin.

Anti-invasive efficacy and survival benefit of the YAP-TEAD inhibitor verteporfin in preclinical glioblastoma models.

BACKGROUND: Glioblastoma (GBM) remains a largely incurable disease as current therapy fails to target the invasive nature of glioma growth in disease progression and recurrence. Here, we use the FDA-approved drug and small molecule Hippo inhibitor Verteporfin (VP) to target YAP-TEAD activity, known to mediate convergent aspects of tumor invasion/metastasis, and assess the drug's efficacy and survival benefit in GBM models. METHODS: Up to 8 low-passage patient-derived GBM cell lines with distinct genomic drivers, including 3 primary/recurrent pairs, were treated with VP or vehicle (VEH) to assess in vitro effects on proliferation, migration, invasion, YAP-TEAD activity, and transcriptomics. Patient-derived orthotopic xenograft (PDX) models were used to assess VP's brain penetrance and effects on tumor burden and survival. RESULTS: VP treatment disturbed YAP/TAZ-TEAD activity; disrupted transcriptome signatures related to invasion, epithelial-to-mesenchymal, and proneural-to-mesenchymal transition, phenocopying TEAD1-knockout effects; and impaired tumor migration/invasion dynamics across primary and recurrent GBM lines. In an aggressive orthotopic PDX GBM model, short-term VP treatment consistently diminished core and infiltrative tumor burden, which was associated with decreased tumor expression of Ki67, nuclear YAP, TEAD1, and TEAD-associated targets EGFR, CDH2, and ITGB1. Finally, long-term VP treatment appeared nontoxic and conferred survival benefit compared to VEH in 2 PDX models: as monotherapy in primary (de novo) GBM and in combination with Temozolomide chemoradiation in recurrent GBM, where VP treatment associated with increased MGMT methylation. CONCLUSIONS: We demonstrate combined anti-invasive and anti-proliferative efficacy for VP with survival benefit in preclinical GBM models, indicating potential therapeutic value of this already FDA-approved drug if repurposed for GBM patients.

A Multi-center Study on the Reproducibility of Drug-Response Assays in Mammalian Cell Lines.

Evidence that some high-impact biomedical results cannot be repeated has stimulated interest in practices that generate findable, accessible, interoperable, and reusable (FAIR) data. Multiple papers have identified specific examples of irreproducibility, but practical ways to make data more reproducible have not been widely studied. Here, five research centers in the NIH LINCS Program Consortium investigate the reproducibility of a prototypical perturbational assay: quantifying the responsiveness of cultured cells to anti-cancer drugs. Such assays are important for drug development, studying cellular networks, and patient stratification. While many experimental and computational factors impact intra- and inter-center reproducibility, the factors most difficult to identify and control are those with a strong dependency on biological context. These factors often vary in magnitude with the drug being analyzed and with growth conditions. We provide ways to identify such context-sensitive factors, thereby improving both the theory and practice of reproducible cell-based assays.

Integrating Transcriptomic Data with Mechanistic Systems Pharmacology Models for Virtual Drug Combination Trials.

Monotherapy clinical trials with mutation-targeted kinase inhibitors, despite some success in other cancers, have yet to impact glioblastoma (GBM). Besides insufficient blood-brain barrier penetration, combinations are key to overcoming obstacles such as intratumoral heterogeneity, adaptive resistance, and the epistatic nature of tumor genomics that cause mutation-targeted therapies to fail. With now hundreds of potential drugs, exploring the combination space clinically and preclinically is daunting. We are building a simulation-based approach that integrates patient-specific data with a mechanistic computational model of pan-cancer driver pathways (receptor tyrosine kinases, RAS/RAF/ERK, PI3K/AKT/mTOR, cell cycle, apoptosis, and DNA damage) to prioritize drug combinations by their simulated effects on tumor cell proliferation and death. Here we illustrate a first step, tailoring the model to 14 GBM patients from The Cancer Genome Atlas defined by an mRNA-seq transcriptome, and then simulating responses to three promiscuous FDA-approved kinase inhibitors (bosutinib, ibrutinib, and cabozantinib) with evidence for blood-brain barrier penetration. The model captures binding of the drug to primary targets and off-targets based on published affinity data and simulates responses of 100 heterogeneous tumor cells within a patient. Single drugs are marginally effective or even counterproductive. Common copy number alterations (PTEN loss, EGFR amplification, and NF1 loss) have a negligible correlation with single-drug or combination efficacy, reinforcing the importance of postgenetic approaches that account for kinase inhibitor promiscuity to match drugs to patients. Drug combinations tend to be either cytostatic or cytotoxic, but seldom both, highlighting the need for considering targeted and nontargeted therapy. Although we focus on GBM, the approach is generally applicable.

Validating Antibodies for Quantitative Western Blot Measurements with Microwestern Array.

Fluorescence-based western blots are quantitative in principal, but require determining linear range for each antibody. Here, we use microwestern array to rapidly evaluate suitable conditions for quantitative western blotting, with up to 192 antibody/dilution/replicate combinations on a single standard size gel with a seven-point, two-fold lysate dilution series (~100-fold range). Pilot experiments demonstrate a high proportion of investigated antibodies (17/24) are suitable for quantitative use; however this sample of antibodies is not yet comprehensive across companies, molecular weights, and other important antibody properties, so the ubiquity of this property cannot yet be determined. In some cases microwestern struggled with higher molecular weight membrane proteins, so the technique may not be uniformly applicable to all validation tasks. Linear range for all validated antibodies is at least 8-fold, and up to two orders of magnitude. Phospho-specific and total antibodies do not have discernable trend differences in linear range or limit of detection. Total antibodies generally required higher working concentrations, but more comprehensive antibody panels are required to better establish whether this trend is general or not. Importantly, we demonstrate that results from microwestern analyses scale to normal "macro" western for a subset of antibodies.

Expression of Carcinoembryonic Cell Adhesion Molecule 6 and Alveolar Epithelial Cell Markers in Lungs of Human Infants with Chronic Lung Disease.

The membrane protein carcinoembryonic antigen cell adhesion molecule (CEACAM6) is expressed in the epithelium of various tissues, participating in innate immune defense, cell proliferation and differentiation, with overexpression in gastrointestinal tract, pancreatic and lung tumors. It is developmentally and hormonally regulated in fetal human lung, with an apparent increased production in preterm infants with respiratory failure. To further examine the expression and cell localization of CEACAM6, we performed immunohistochemical and biochemical studies in lung specimens from infants with and without chronic lung disease. CEACAM6 protein and mRNA were increased ~4-fold in lungs from infants with chronic lung disease as compared with controls. By immunostaining, CEACAM6 expression was markedly increased in the lung parenchyma of infants and children with a variety of chronic lung disorders, localizing to hyperplastic epithelial cells with a ~7-fold elevated proliferative rate by PCNA staining. Some of these cells also co-expressed membrane markers of both type I and type II cells, which is not observed in normal postnatal lung, suggesting they are transitional epithelial cells. We suggest that CEACAM6 is both a marker of lung epithelial progenitor cells and a contributor to the proliferative response after injury due to its anti-apoptotic and cell adhesive properties.

Expression of human carcinoembryonic antigen-related cell adhesion molecule 6 and alveolar progenitor cells in normal and injured lungs of transgenic mice.

Carcinoembryonic antigen-related cell adhesion molecule 6 (CEACAM6) is expressed in the epithelium of various primate tissues, including lung airway and alveoli. In human lung, CEACAM6 is developmentally and hormonally regulated, protects surfactant function, has anti-apoptotic activity and is dysregulated in cancers. We hypothesized that alveolar CEACAM6 expression increases in lung injury and promotes cell proliferation during repair. Studies were performed in CEABAC transgenic mice-containing human CEACAM genes. The level of CEACAM6 in adult CEABAC lung was comparable to that in human infants; expression occurred in epithelium of airways and of some alveoli but rarely co-localized with markers of type I or type II cells. Ten days after bleomycin instillation, both the number of CEACAM6(+) cells and immunostaining intensity were elevated in injured lung areas, and there was increased co-localization with type I and II cell markers. To specifically address type II cells, we crossed CEABAC mice with animals expressing EGFP driven by the SP-C promoter. After bleomycin injury, partially flattened, elongated epithelial cells were observed that expressed type I cell markers and were primarily either EGFP(+) or CEACAM6(+). In cell cycle studies, mitosis was greater in CEACAM6(+) non-type II cells versus CEACAM6(+)/EGFP(+) cells. CEACAM6 epithelial expression was also increased after hyperoxic exposure and LPS instillation, suggesting a generalized response to acute lung injuries. We conclude that CEACAM6 expression is comparable in human lung and the CEABAC mouse. CEACAM6 in this model appears to be a marker of a progenitor cell population that contributes to alveolar epithelial cell replenishment after lung injury.