0% READ
FULL TEXT
INTRODUCTION
Cancer is one of the second leading causes of death worldwide because the high expenditure on anticancer drugs becomes a burden on public health (Luo et al., 2023). Cancer can occur due to gene mutation and carcinogenic agents that influence the metabolism and change the normal growth pattern of cells (Kurniawan et al., 2021). In many developed countries, the number of deaths is increasing day by day due to a lack of effective treatment (Wang et al., 2021). According to the WHO report, global death due to cancer recorded in 2020 reached 9.90 million; out of this, 5.53 million were men and 4.43 million were women. The number of new cases of breast cancer has reached around 2.26 million (Arruebo et al., 2011; Chawla and Vaidya, 2025). Therefore, breast cancer has replaced lung cancer as the most common in type of cancer in the world (Bai et al., 2023). The research of anticancer medicine is needed because the suffering and number of death cases are unstoppable till now. Cancer treatment is one of the most challenging research areas in the medicinal and pharmacological fields (Brunner, 2016). Currently; there are more than 200 drugs under clinical use. All anticancer drugs are classified into two categories: target-based and toxicity-based agents. Recently, pyridopyrimidine derivatives are under clinical trial and approved for treatment of different kinds cancer, such as chronic myelogenous leukemia, CDK4/CDK6 inhibitors, MEK allosteric site inhibitors and amino pyridopyrimidine as inhibitors of cellular proliferation (YadavandShah, 2022). The resistance is a major concern for cancer patients and a challenge for cancer research scientists (Abourehab et al., 2022). Recently, targeting the cell cycle and mitosis process is an alternative therapeutic strategy for cancer treatment (Dominguez-Brauer et al., 2015). Monopolar spindle-1 is one of the most important targets of protein kinase responsible for regulating the various steps in mitosis, centrosomes, and several different signaling pathways (Yang et al., 2019; Sflakidou et al., 2024; Kang et al., 2007). MPS-1 is overexpressed in numerous types of cancer, such as breast cancer (Daniel et al., 2011; Lee et al., 2013), thyroid cancer (He et al., 2015), pancreatic cancer (Tannous et al., 2013; Skoulidis et al., 2022), glioblastoma (Kitajima et al., 2022), lung cancer (Maciejowski et al., 2010), and its role in neuroblastoma is still unknown. MPS1 is involved in the regulation of chromosome alignment during metaphase along with checkpoint function (Fabritius et al., 2024). Therefore, so many MPS1 inhibitors have been developed and shown anticancer activity, but safety has not reported that is why no MPS1 inhibitors have been entered into clinical testing in cancer patients. (Sugimoto et al., 2017) We are investigating the development of more effective anticancer agents by using concepts of computational chemistry, such as computer-aided drug design that could save time and the cost of synthesis for discovering new drug candidates (Pugh et al., 2022). Molecular docking is one of the tools that predict the binding affinity and best binding poses of a molecule with the target (MPS1) (Xing et al., 2021). We explored the design of pyridopyrimidine derivatives based on in silico virtual screening for an efficient method for drug development research (Hua et al., 2024). This study examined the design of pyridopyrimidine derivatives, ADMET prediction, and binding potential against MPS1.
MPS1 is an important regulator of facilitate massisive chromosome mis segregation during mitosis with subsequent reentry to G1 phase (Shakya et al., 2024; Gurjar, 2024) MPS1 have properties of activating cGAS-STING signaling in proliferating KL cell and facilitate tumor suppressive function.
The present work aims to explore the potential drug candidates against MPS1 to control the growth of abnormal cancer cells. The molecular docking analysis with the PyRx virtual screening tool was performed against MPS1. Moreover, ADMET analysis was performed for pyridopyrimidine derivatives to evaluate physicochemical, drug-likeness, lipophilicity, toxicity, and pharmacokinetic attributes inside the biological system using Protox-III and the Swiss ADME web server. PROTOX-III is a computational technique for predicting the toxicity of small molecules by analyzing their similarity to toxic substances and structural characteristics. PkCSM helps to predict pharmacokinetic features and drug metabolism, whereas Swiss ADME predicts ADME-T properties like lipophilicity and drug-likeness.
MATERIALS AND METHODS
Molecular docking studies
Legands preparation
A total 150 molecules designed and their structure was drawn through CS Chem Draw Ultra and saved in the mol files format. (YadavandShah, 2021) The structure of standard drug (palbociclib) was taken from Pub Chem database (https://pubchem.ncbi.nlm.nih.gov). All compounds were transferred to PyRx virtual screening tool by using Open Babel one by one then minimized energy and converted into auto-dock ligand pdbqt files format.
Gride preparation
The receptor gride was used to define the specific region of interaction of protein and ligand molecules. This tool was used to determine the region surrounding the active center based on x, y and z coordinate measurements. The resolution of the grid box was centered on the coordinates 5.83, 43.43, and 23.22, which corresponded to the x, y and z axis respectively.
Protein preparation
The crystal structure of target protein that is MPS1 (Monopolar spindle1) (PDB ID-6H3K) was downloaded from the Research Collaboratory for structure bioinformatics Protein Data Bank (RCSB PDB). The protein structure was prepared for docking by removal of water molecules and co-crystallized ligand; the polar hydrogen atom was added and file saved in the pdbqt format.
Molecular docking
For molecular docking, PyRx virtual screening tool (autodock vina) version 8.0 was used. A maximum of ten confirmations were set for each compound using (genetic logarithm) lamarckian, 3D prediction of obtained result. Discovery studio client 2021 was used for visualization of 2D and 3D interaction of docked molecules and protein target. All of the ligands were docked independently, lip dock score was used for analyze various poses of all ligands.
Docking calculation
For docking calculation, discovery studio client 2021 software was used. The target protein MSP1 followed by (PDB ID 6H3K) and resolution of this protein was 2.8Aº. On the basis of lowest score of binding energy, the top ranked conformations under gone for visualization and atomic level interaction investigation.
ADMET prediction
The ADMET analysis of the 150 compounds was screened by using three different web servers for analysis of important pharmacokinetics parameters. The web servers are as follows, https://lmmd.ecust.edu.cn/admetsar2, https://biosig.lab.uq.edu.au/pkcsm/prediction_single/adme http://www.swissadme.ch/, and https://tox.charite.de/protox3/. The descriptors such as molecular weight, hydrogen bond donors, hydrogen bond acceptors, QclogPo/w, QlogS, HIA, BBB permeability, GI absorption and lipinski’s rule five (Daina et al., 2017) Based on this result, we selected all compounds for the molecular docking study.
RESULTS
Consequently, these compounds might have a favorable pharmacokinetics profile based on prediction of the ADMET characteristics (Tables 1 and 2). All 150 molecules were screened for using docking studies to identify the most effective pyridopyrimidine bioactive anticancer agents. The top 18 molecules were identified from all of them based on binding energy and hydrogen bonding. The detail of all compounds, the evaluation of the molecular level of interaction showed (Table 3) that the designed scaffold can bind with amino acid residues present in the receptor active domain. The binding pose view of the selected molecules is depicted in Figure 1. but all of them, the docking results showed that P7 had the highest binding score (-9.7kcal/mol).
| Compound code | Molecular structure | Compound code | Molecular structure |
|---|---|---|---|
| P1 | P5 | ||
| P2 | P6 | ||
| P3 | P7 | ||
| P4 | P8 | ||
| P9 | P13 | ||
| P10 | P14 | ||
| P11 | P15 | ||
| P12 | P16 | ||
| P17 | P18 | ||
| Control (Palbociclib) |
| Sl. No. | Code No | Binding score(kca/mol) 6H3K | Rotatable bond | GI absorption | BBB permeability | P-gp substrate | CYP3A4 inhibitors | Log kp (cm/s) | Log S solubility | B Sa | S Ab | Toxicitypridiction (Pro-tox-3.0 software) | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Carcinigenecit | Mutagenicity | ||||||||||||
| 1 | P1 | -6.7 | 4 | High | No | No | Yes | -7.82 | -4.27 | 0.55 | 4.04 | No | No |
| 2 | P2 | -8 | 4 | High | No | Yes | Yes | -7.87 | -4.8 | 0.55 | 3.89 | No | No |
| 3 | P3 | -7.1 | 4 | High | No | Yes | Yes | -7.6 | -3.95 | 0.55 | 4.02 | No | No |
| 4 | P4 | -9.7 | 4 | High | No | Yes | Yes | -7.64 | -3.76 | 0.55 | 3.87 | No | No |
| 5 | P5 | -9.1 | 4 | High | No | No | Yes | -7.66 | -3.66 | 0.55 | 4.13 | No | No |
| 6 | P6 | -9.7 | 4 | High | No | Yes | Yes | -7.7 | -3.47 | 0.55 | 3.98 | No | No |
| 7 | P7 | -8.7 | 4 | High | No | Yes | Yes | -7.36 | -4.55 | 0.55 | 4.14 | No | No |
| 8 | P8 | -6.8 | 4 | High | No | Yes | Yes | -7.4 | -4.36 | 0.55 | 3.99 | No | No |
| 9 | P9 | -8.1 | 5 | Low | No | Yes | No | -8.23 | -3.42 | 0.55 | 4.13 | No | No |
| 10 | P10 | -7 | 5 | Low | No | Yes | No | -8.27 | -3.23 | 0.55 | 3.99 | No | No |
| 11 | P11 | -6.9 | 5 | High | No | Yes | Yes | -7.29 | -4.22 | 0.55 | 4.39 | No | No |
| 12 | P12 | -7.2 | 5 | High | No | Yes | Yes | -7.33 | -4.03 | 0.55 | 4.2 | No | No |
| 13 | P13 | -9.2 | 6 | Low | No | Yes | Yes | -8.21 | -3.53 | 0.55 | 4.34 | No | No |
| 14 | P14 | -6.6 | 6 | Low | No | Yes | Yes | -8.26 | -3.34 | 0.55 | 4.19 | No | No |
| 15 | P15 | -6.9 | 6 | Low | No | Yes | Yes | -8.59 | -3.37 | 0.55 | 4.39 | No | No |
| 16 | P16 | -6.6 | 6 | High | No | Yes | No | -8.63 | -3.18 | 0.55 | 4.24 | No | No |
| 17 | P17 | -6.6 | 5 | High | No | Yes | No | -8.39 | -3.29 | 0.55 | 4.22 | No | No |
| 18 | P18 | -8.6 | 5 | Low | No | Yes | No | -8.43 | -3.1 | 0.55 | 4.07 | No | No |
| Control | palbociclib | -7.9 | 5 | High | No | Yes | Yes | -7.74 | -3.78 | 0.55 | 3.57 | No | No |
| Code of compounds | PDB ID | H Bond interaction | Hydrophobic interaction |
|---|---|---|---|
| P4 | 6H3K | MET A 671, ASP A 664 | ILE A531, LEU A 654, ALA A551, MET A 602, VAL A 539, ASP A 664, ILE A 663 |
| P7 | GLY A 605, CYS A 604, MET A 671 | PRO A 673, ALA A 651, ILE A 531, ILE A 664, ILE A 663, ASP A 664, ILE A 531 | |
| P8 | CYS A 604, ASP A 664 | ALA A 651, MET A 673, ILE A 663, LYS A 553, MET A 602, LEU A 654, ILE A 532 | |
| P13 | ASP A 664, CYS A 604 | LEU A 654, ILE A 531, VEL A 531, LYS A 553, ILE A 663, MET A 671, ASP A 664 | |
| Control (palbociclib) | ASN A 606, GLY A 605 | LEU A 654, ILE A 531, CYS A 604 |
Against MPS1. Docking analysis of pyridopyrimidine with MPS1 showed that strong molecular interaction such as hydrogen and hydrophobic bonding with amino acids GLY A 605, CYS A 604, MET A 671, PRO A 673, ALA A 651, ILE A 531, ILE A 664, ILE A 663, ASP A 664, ILE A 531 (Figure 1). Still, detailed studies are required to understand the binding activity further. Based on this work, more research could be done to develop new drug candidate that greatly help to treat cancer.
DISCUSSION
All the designed pyrido-pymidines derivatives have wide range of medicinal and pharmacological importance in the area of research. Drug discovery and development process is an expensive and time taken phenomenon. The molecular docking and ADME-T analysis will help to identity potential drug candidate before going experimental procedure, this techniques will reduce time and man power required for synthesis of millions of molecules and wastage of chemicals. The ADME-T and drug-likeness prediction of all our designed pyridopyrimidines had investigated by using rotox-III, pkCSM, and the SwissADME web server. Based on the determined value, the maximum of compounds showed all descriptors of Lipinski rule of five, mention detailed in Table 4. The primary factors of significance are the molecular weight below 500 Daltons and log p value ranged from (1.93 to 3.02) indicates that passive intestinal absorption. The ADMET parameter is crucial for lipophilicity and capacity to cross critical biological barrier and membranes. Remember, if molecules with high lipophilicity show poor solubility, affecting bioavailability or accumulate in adipose tissues cause slow excretion and influences systemic toxicity. An important determinant oral bioavailability depends on total number of hydrogen donors and hydrogen acceptors. The calculated all parameters for all tested compounds came under standards limits.
| Lipinski rule of five and physicochemical properties | |||||||
|---|---|---|---|---|---|---|---|
| Descriptors | MW | HBD | HBA | cLogPo/w | No of rule violation | druglikeness | |
| Sl. No. | Code No | ≤500 Dalton | ≤5 | ≤10 | ≤5 | ≤ 2 violations | |
| 1 | P1 | 483.31 | 2 | 4 | 2.91 | 0 | Yes |
| 2 | P2 | 469.29 | 3 | 4 | 2.36 | 0 | Yes |
| 3 | P3 | 438.86 | 2 | 4 | 2.84 | 0 | Yes |
| 4 | P4 | 424.84 | 3 | 4 | 2.19 | 0 | Yes |
| 5 | P5 | 418.45 | 2 | 4 | 2.78 | 0 | Yes |
| 6 | P6 | 404.42 | 3 | 4 | 2.05 | 0 | Yes |
| 7 | P7 | 473.31 | 2 | 4 | 3.02 | 0 | Yes |
| 8 | P8 | 459.28 | 3 | 4 | 2.23 | 0 | Yes |
| 9 | P9 | 449.42 | 2 | 6 | 2.36 | 1 | Yes |
| 10 | P10 | 435.39 | 3 | 6 | 1.93 | 1 | Yes |
| 11 | P11 | 446.5 | 2 | 4 | 2.81 | 0 | Yes |
| 12 | P12 | 432.47 | 3 | 4 | 2.7 | 0 | Yes |
| 13 | P13 | 464.42 | 3 | 6 | 2.81 | 0 | Yes |
| 14 | P14 | 450.44 | 4 | 6 | 2.62 | 0 | Yes |
| 15 | P15 | 480.47 | 3 | 7 | 2.95 | 1 | Yes |
| 16 | P16 | 466.44 | 4 | 7 | 2.36 | 1 | Yes |
| 17 | P17 | 450.44 | 3 | 6 | 2.54 | 0 | Yes |
| 18 | P18 | 436.42 | 4 | 6 | 2.01 | 0 | Yes |
| Control | Palbociclib | 447.53 | 2 | 6 | 2.23 | 1 | Yes |
CONCLUSION
In this research work, we designed all 150 different derivatives of Pyrido-pyrimidines motifs but we are showing best docked score of 18 compounds and evaluated their molecular docking studies with receptor for anti-cancer activity through in silico studies. All of the compounds in this series showed favorable pharmacokinetic properties such as Absorption, Distribution, Metabolism, Elimination, and Toxicity (ADME-T). According to docking studies, substances P4, P7, P8, and P13 exhibited the highest binding affinities when compared to palbociclib (with mono-polar spindle-1 (PDB ID-6H3K)). Docking and ADMET studies demonstrated that these ligands interacted mostly with the MPS1 receptor for anticancer reports by forming hydrogen bonds and interacting hydrophobically with the domain. The above study signifies if the suggested compounds screened. They may have promising anticancer activity.
