Computational investigation of adenosine 5′-(α,β-methylene)-diphosphate (AMPCP) derivatives as ecto-5′-nucleotidase (CD73) inhibitors by using 3D-QSAR, molecular docking, and molecular dynamics simulations

CD73, as a surface enzyme anchored on the outside of the cell membrane via glycosylphosphatidylinositol (GPI), can convert the AMP in the tumor cell microenvironment into adenosine to promote the growth of tumor cells. It has been overexpressed in many different types of human tumors, such as gastric cancer, pancreatic cancer, liver cancer, and other tumor cells. Therefore, targeted inhibitors of CD73 are considered potential tumor treatment methods. Due to the low bioavailability of nucleoside CD73 inhibitors, it is necessary to develop new inhibitors. In this study, through molecular docking, three-dimensional quantitative structure–activity relationship (3D-QSAR) and molecular dynamics (MD) simulations, a series of CD73 inhibitors were calculated and studied to reveal their structure–activity relationships. Through molecular docking studies, the possible mode of interaction between inhibitors and protein is explored. Subsequently, a 3D-QSAR model was established by comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). For the best CoMFA model, the Q 2 and R 2 values are 0.708 and 0.983, respectively, while for the best CoMSIA model, the Q 2 and R 2 values are 0.809 and 0.992, respectively. Based on the contour maps, we designed ten new CD73 inhibitors and predicted their activity by the model, all of them are better than molecules in the dataset. In addition, in order to select potential drug candidates, ADMET prediction was performed on template molecules and designed compounds. Moreover, the stability of the complex formed by the two inhibitors and CD73 was evaluated by molecular dynamics simulation, and the results are consistent with the results of molecular docking and 3D-QSAR research. Finally, the binding free energy was calculated by the surface area method (MM-GBSA), and the results are consistent with the activities that van der Waals and Coulomb contribute the most during the binding process of the molecule to the CD73 protein. In conclusion, our research provides valuable information for the further development of CD73 inhibitors.