Biosimulation Using the Cellworks Computational Omics Biology Model (CBM) Identifies Novel Biomarkers to Inform Mitoxantrone, Etoposide, and Cytarabine (MEC)-Based Combination Therapy in Refractory & Relapsed Acute Myeloid Leukemia (AML) Patients

Background: The optimal treatment strategy for managing Acute Myeloid Leukemia (AML) and the use of reliable and predictive biomarkers to guide selection of cytotoxic chemotherapy regimens among patients with diverse genomic profiles remain unmet needs in the clinic. The combination of MEC [mitoxantrone (MIT), etoposide (VP16), and cytarabine (ARA-C)] is a commonly used regimen for relapsed or refractory AML patients. Unfortunately, many patients do not respond to MEC, and which of the three drug agents matters most for each individual patient is not known. Predictors of response are needed urgently. Methods: The Computational Omics Biology Model (CBM) is a computational multi-omic biology software model created using artificial intelligence heuristics and literature sourced from PubMed to generate a patient-specific protein network map. The aberration and copy number variations from individual cases served as input into the CBM. Disease-biomarkers unique to each patient were identified within patient-specific protein network maps. Biosimulations were conducted on the Cellworks Biosimulation Platform by measuring the effect of chemotherapy on a cell growth score comprised of a composite of cell proliferation, viability, apoptosis, metastasis, and other cancer hallmarks. Biosimulation of drugs was conducted by mapping the interaction of various drug combinations with the patient's genomic and pathway alterations based on signaling pathway mechanisms and their phenotypic consequences. The Cellworks Biosimulation Platform identified unique chromosomal signatures that permit a stratification of patients that are most likely to respond to MIT, VP16, or ARA-C as well as their combinations. 65 AML patients were selected for this study largely based on genomic data published in TCGA and PubMed: •ARA-C [N=12, 7 responders (R) & 5 non-responders (NR)]•ARA-C + MIT [N=30, 29 R & 1 NR]•ARA-C + MIT + VP16 [N=23, 12 R & 11 NR] Results: Of the12 patients treated with ARA-C alone, 5 were predicted to be NR and 7 were predicted to be R. Of the 5 NR, 4 had 5q del which resulted in loss of APC, CSNK1A1 and SLC22A4 (nucleoside carrier) forming the non-response biomarkers for ARA-C. Notably, the biosimulation predicted lenalidomide to be beneficial for these patients. Out of 7 R, 4 patients also had 5q del, but were predicted to be R because of co-occurring aberrations involving CLSPN del, DHODH del, MSH2 del, EP300 del, CREBBP del, MSH6 Del, and RRM2 del. These genes were excl