In-silico probing of AML related RUNX1 cancer-associated missense mutations: Predicted relationships to DNA binding and drug interactions

The molecular consequences of cancer associated mutations in Acute myeloid leukemia (AML) linked factors are not very well understood. Here, we interrogated the COSMIC database for missense mutations associated with the RUNX1 protein, that is frequently mis-regulated in AML, where we sought to identify recurrently mutated positions at the DNA-interacting interface. Indeed, six of the mutated residues, out of a total 417 residues examined within the DNA binding domain, evidenced reduced DNA association in in silico predictions. Further, given the prominence of RUNX1's compromised function in AML, we asked the question if the mutations themselves might alter RUNX1's interaction (off-target) with known FDA-approved drug molecules, including three currently used in treating AML. We identified several AML-associated mutations in RUNX1 that were calculated to enhance RUNX1's interaction with specific drugs. Specifically, we retrieved data from the COSMIC database for cancer-associated mutations of RUNX1 by using R package "data.table " and "ggplot2 " modules. In the presence of DNA and/or drug, we used docking scores and energetics of the complexes as tools to evaluate predicted interaction strengths with RUNX1. For example, we performed predictions of drug binding pockets involving Enasidenib, Giltertinib, and Midostaurin (AML associated), as well as ten different published cancer associated drug compounds. Docking of wild type RUNX1 with these 13 different cancer-associated drugs indicates that wild-type RUNX1 has a lower efficiency of binding while RUNX1 mutants R142K, D171N, R174Q, P176H, and R177Q suggested higher affinity of drug association. Literature evidence support our prediction and suggests the mutation R174Q affects RUNX1 DNA binding and could lead to compromised function. We conclude that specific RUNX1 mutations that lessen DNA binding facilitate the binding of a number of tested drug molecules. Further, we propose that molecular modeling and docking studies for RUNX1 in the presence of DNA and/or drugs enables evaluation of the potential impact of RUNX1 cancer associated mutations in AML.