In silico targeting of colony-stimulating factor-1 receptor: delineating immunotherapy in cancer

Explor Target Antitumor Ther. 2023;4(4):727-742. doi: 10.37349/etat.2023.00164. Epub 2023 Aug 31.

Abstract

Aim: Delineate structure-based inhibition of colony-stimulating factor-1 receptor (CSF1R) by small molecule CSF1R inhibitors in clinical development for target identification and potential lead optimization in cancer therapeutics since CSF1R is a novel predictive biomarker for immunotherapy in cancer.

Methods: Compounds were in silico modelled by induced fit docking protocol in a molecular operating environment (MOE, MOE.v.2015). The 3-dimensional (3D) X-ray crystallized structure of CSF1R kinase (Protein Databank, ID 4R7H) was obtained from Research Collaboratory for Structural Bioinformatics (RSCB) Protein Databank. The 3D conformers of edicotinib, DCC-3014, ARRY-382, BLZ-945, chiauranib, dovitinib, and sorafenib were obtained from PubChem Database. These structures were modelled in Amber10:EHT molecular force field, and quick prep application was used to correct and optimize the structures for missing residues, H-counts, termini capping, and alternates. The binding site was defined within the vicinity of the co-crystallized ligand of CSF1R kinase. The compounds were docked by the triangular matcher placement method and ranked by the London dG scoring function. The docked poses were further refined by the induced fit method. The pose with the lowest binding score (ΔG) was used to model the ligand interaction profile in Discovery Studio Visualizer v17.2. The co-crystallized ligand was docked in its apo conformation, and root-mean-square deviation was computed to validate the docking protocol.

Results: All 7 CSF1R inhibitors interact with residue Met637 exhibiting selectivity except for edicotinib. The inhibitors maintain CSF1R in an auto-inhibitory conformation by interacting with Asp797 of the Asp-Phe-Gly (DFG) motif and/or hindering the conserved salt bridge formed between Glu633 and Lys616 thus stabilizing the activation loop, or interacting with tryptophan residue (Trp550) in the juxtamembrane domain. DCC-3014, ARRY-382, BLZ-945, and sorafenib bind with the lowest binding energy with CSF1R kinase.

Conclusions: Pyrimidines are potent inhibitors that interact with CSF1R residues. DCC-3014 and ARRY-382 exhibit exceptional pharmaceutical potential exhibiting great structural stability and affinity.

Keywords: Cancer drug development; computational modelling; drug resistance; lead optimization; molecular docking; target identification; targeted cancer therapy; tyrosine kinase inhibitors.