Epithelial-to-mesenchymal transition (EMT) is discussed to be centrally involved in invasion, stemness, and drug resistance. Experimental models to evaluate this process in its biological complexity are limited. To shed light on EMT impact and test drug response more reliably, we use a lung tumor test system based on a decellularized intestinal matrix showing more in vivo-like proliferation levels and enhanced expression of clinical markers and carcinogenesis-related genes. In our models, we found evidence for a correlation of EMT with drug resistance in primary and secondary resistant cells harboring KRAS^G12C or EGFR mutations, which was simulated in silico based on an optimized signaling network topology. Notably, drug resistance did not correlate with EMT status in KRAS-mutated patient-derived xenograft (PDX) cell lines, and drug efficacy was not affected by EMT induction via TGF-β. To investigate further determinants of drug response, we tested several drugs in combination with a KRAS^G12C inhibitor in KRAS^G12C mutant HCC44 models, which, besides EMT, display mutations in P53, LKB1, KEAP1, and high c-MYC expression. We identified an aurora-kinase A (AURKA) inhibitor as the most promising candidate. In our network, AURKA is a centrally linked hub to EMT, proliferation, apoptosis, LKB1, and c-MYC. This exemplifies our systemic analysis approach for clinical translation of biomarker signatures.

To improve and focus preclinical testing, we combine tumor models based on a decellularized tissue matrix with bioinformatics to stratify tumors according to stage-specific mutations that are linked to central cancer pathways. We generated tissue models with BRAF-mutant colorectal cancer (CRC) cells (HROC24 and HROC87) and compared treatment responses to two-dimensional (2D) cultures and xenografts. As the BRAF inhibitor vemurafenib is—in contrast to melanoma—not effective in CRC, we combined it with the EGFR inhibitor gefitinib. In general, our 3D models showed higher chemoresistance and in contrast to 2D a more active HGFR after gefitinib and combination-therapy. In xenograft models murine HGF could not activate the human HGFR, stressing the importance of the human microenvironment. In order to stratify patient groups for targeted treatment options in CRC, an in silico topology with different stages including mutations and changes in common signaling pathways was developed. We applied the established topology for in silico simulations to predict new therapeutic options for BRAF-mutated CRC patients in advanced stages. Our in silico tool connects genome information with a deeper understanding of tumor engines in clinically relevant signaling networks which goes beyond the consideration of single drivers to improve CRC patient stratification.

Precision medicine has changed thinking in cancer therapy, highlighting a better understanding of the individual clinical interventions. But what role do the drivers and pathways identified from pan-cancer genome analysis play in the tumor? In this letter, we will highlight the importance of in silico modeling in precision medicine. In the current era of big data, tumor engines and pathways derived from pan-cancer analysis should be integrated into in silico models to understand the mutational tumor status and individual molecular pathway mechanism at a deeper level. This allows to pre-evaluate the potential therapy response and develop optimal patient-tailored treatment strategies which pave the way to support precision medicine in the clinic of the future.

High attrition-rates entailed by drug testing in 2D cell culture and animal models stress the need for improved modeling of human tumor tissues. In previous studies our 3D models on a decellularized tissue matrix have shown better predictivity and higher chemoresistance. A single porcine intestine yields material for 150 3D models of breast, lung, colorectal cancer (CRC) or leukemia. The uniquely preserved structure of the basement membrane enables physiological anchorage of endothelial cells and epithelial-derived carcinoma cells. The matrix provides different niches for cell growth: on top as monolayer, in crypts as aggregates and within deeper layers. Dynamic culture in bioreactors enhances cell growth. Comparing gene expression between 2D and 3D cultures, we observed changes related to proliferation, apoptosis and stemness. For drug target predictions, we utilize tumor-specific sequencing data in our in silico model finding an additive effect of metformin and gefitinib treatment for lung cancer in silico, validated in vitro. To analyze mode-of-action, immune therapies such as trispecific T-cell engagers in leukemia, as well as toxicity on non-cancer cells, the model can be modularly enriched with human endothelial cells (hECs), immune cells and fibroblasts. Upon addition of hECs, transmigration of immune cells through the endothelial barrier can be investigated. In an allogenic CRC model we observe a lower basic apoptosis rate after applying PBMCs in 3D compared to 2D, which offers new options to mirror antigen-specific immunotherapies in vitro. In conclusion, we present modular human 3D tumor models with tissue-like features for preclinical testing to reduce animal experiments.