0% READ
FULL TEXT
INTRODUCTION
In 2020 alone, an estimated 19.3 million new cases and 10 million deaths from cancer were reported, making it one of the most serious health issues in the globe (Sung et al., 2021). Despite advancements in early detection, diagnostic techniques and treatment modalities, the high mortality rate and the significant burden of disease necessitate continuous exploration for novel therapeutic targets and strategies (Liu et al., 2024). One of the promising avenues in cancer therapy is the targeting of specific molecular pathways that are dysregulated in cancer cells, principal to uncontrolled proliferation and survival. Among these, protein kinases have emerged as critical regulators of cell signalling and are thus attractive targets for cancer therapy (Bhullar et al., 2018).
Aurora kinases, a family of serine/threonine kinases, show a essential role in the regulation of mitosis (Fu et al., 2007). They ensure proper chromosome alignment, segregation and cytokinesis, thus maintaining genomic stability (Roshan et al., 2023). This family comprises three members: Aurora A, Aurora B and Aurora C, each with distinct yet overlapping functions during cell division (Willems et al., 2018). Aurora Several crucial processes in mitosis, such as spindle formation, chromosomal alignment, and centrosome maturation and separation, depend on a kinase in particular (Magnaghi-Jaulin et al., 2019).
Overexpression and hyperactivation of Aurora A have been frequently observed in various cancers, including breast (Ingebriktsen et al., 2024), hepatocellular (Jeng et al., 2004), esophageal squamous cell carcinoma (Wang et al., 2006), head and neck (Mehra et al., 2013), and pancreatic cancers (Zhang et al., 2022). This aberrant expression correlates with poor prognosis, increased tumor violence, and resistance to conventional therapies (Zhang et al., 2004; Zheng et al., 2023). As a result, Aurora A kinase has gained significant attention as a potential therapeutic aim in oncology.
Aurora A kinase is located on chromosome 20q13, a region generally amplified in cancers (Keen and Taylor, 2004; Kitzen et al., 2010; Sen et al., 1997). Its overexpression is associated with aneuploidy and chromosomal instability, hallmark features of cancer cells (Takeshita et al., 2013). Aurora A contributes to the malignant phenotype by disrupting normal mitotic processes, leading to abnormal centrosome function, spindle assembly defects and improper chromosome segregation (Goldenson et al., 2015; Sasai et al., 2008). These aberrations result in genomic instability, a driving force behind cancer progression and metastasis (Torchia et al., 2009). Moreover, Aurora A interacts with several key oncogenic and tumor suppressor pathways (Lin et al., 2020). Aurora A also inactivates the tumor suppressor protein p53 by phosphorylation, leading to impaired DNA damage response and apoptosis (Katayama et al., 2012). These interactions determine the role of Aurora A in cancer development and highlight its possible as a therapeutic target.
Targeting the Aurora A kinase presents a promising therapeutic strategy due to its central role in mitosis and cancer pathogenesis. Inhibitors of Aurora A aim to disrupt its kinase activity, thereby inducing mitotic defects, cell cycle arrest, and apoptosis in cancer cells. The development of Aurora A inhibitors, along with combination strategies and novel inhibitors with improved selectivity holds promise for enhancing therapeutic outcomes and overcoming resistance.
MATERIALS AND METHODS
Retrieval and preparation of the target molecule
The crystal structure of human Aurora A protein kinase was downloaded from RCSB PDB database, with PDB ID: 1MQ4 (Nowakowski et al., 2002). The structure was prepared using a standard receptor protocol. Swiss-PDB Viewer tools were used to check and correct any structural issues in the protein molecule; such include erroneous bond arrangements or missing atoms, or unresolved residues. Further, energy minimization was performed to repair distorted geometries via moving atoms and releasing internal constraints on the protein structure. Hydrogen atoms and Kollman United Atom Charges were then additional to the structure. Finally, the structure was saved in. pdbqt format using MGLTools.
Retrieval and preparation of the ligand molecules
A small isoflavonoid library of 39 plant-derived natural compounds was curated and docked in contradiction of the target protein. The structures of these ligand molecules were downloaded in mol format and converted to the PDB format using Open Babel tools. Hydrogen atoms were then added to the ligand molecules using PyMOL. Swiss-PDB viewer was employed to carry out energy minimization and relax the ligand structure to its lowest energy conformation. The structure of each ligand molecule was verified for any missing atoms, incorrect bond orders, or other structural issues using molecular modeling software. Appropriate atomic charges were subsequently assigned to the ligand molecules. Finally, the optimized and corrected ligand structures were saved in PDBQT format using MGLTools, which is essential for docking using AutoDock Vina.
Molecular docking
We utilized AutoDock Vina for the molecular docking process, performing blind docking of ligand molecules with the Aurora A protein kinase (Trott et al., 2010). The grid dimensions for the X, Y and Z coordinates were set to 67, 83, and 88 Å, respectively, centered at -42.34, 16.21, and 71.13. The grid spacing was fixed at 1.00 Å with an exhaustiveness parameter of 8. Binding affinity was evaluated and the resulting docked complexes were visualized by means of PyMOL and Discovery Studio Visualizer to analyze the interaction patterns between the ligand molecules and the Aurora A protein kinase.
ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) characteristics of Duartin
To assess the ADMET characteristics of a ligand molecule, we utilized online tools such as SwissADME (http://www.swissadme.ch/index.php) and pkCSM (https://biosig.lab.uq.edu.au/pkcsm/). Integrating the results from both SwissADME and pkCSM provides a thorough understanding of the ADMET properties of the top hit.
Swiss Target Prediction tool and PASS online are used to predict Duartin's activity spectrum
This study utilized PASS Online (http://www.pharmaexpert.ru/passonline/) and the Swiss Target Prediction tool (http://www.swisstargetprediction.ch/) to identify the activity spectrum of Duartin. PASS Online (Prediction of Activity Spectra for Substances) is a web-based platform that predicts the biological and pharmacological activities of chemical compounds, helping researchers identify potential activities and properties of small molecules. Swiss Target Prediction is an online program that assists in identifying proteins that might interact with a particular substance by predicting possible targets for small molecules or medications.
Test compound and MDA-MB-231 breast cancer cells
Duartin (>97.0% GC(T)) was obtained from MedKoo Biosciences, Inc. and liquefied in DMSO (SRL) to prepare a stock solution with a concentration of 500 µg/mL. MDA-MB-231 cells were sourced from NCCS, Pune and cultured in DMEM supplemented by means of 10% FBS and 1% antibiotics.
Cell culture and cytotoxicity assay
MDA-MB-231 cells were cultured in DMEM by 10% FBS and 1% penicillin-streptomycin at 37ºC with 5% CO2 until confluency. The cells were then seeded into 96-well plates at 5,000 cells/well and grown for 48 hr. They were treated by Duartin (0-1000 µM) for 72 hr, with control wells left untreated. Subsequently treatment, 100 µL of 5 mg/mL MTT solution was additional to each well and incubated at 37ºC for 3-4 hr. The MTT solution was substituted with 100 µL of DMSO to solubilize the formazan crystals, and absorbance was measured at 570 nm (reference at 630 nm). Cell viability (%) was considered as [Absorbance of treated cells/Absorbance of control cells] × 100. Additionally, a time-dependent viability assay was performed at 24, 48, and 72 hr using the IC50 concentration of Duartin.
mRNA expression of Aurora A protein kinase in Duartin treated MDA-MB-231 breast cancer cells
Breast cancer cells (MDA-MB-231) were treated with the IC50 concentration of Duartin for 48 hr. RNA was taken out using TRIzol Reagent (Thermo Fisher). cDNA was made by using cDNA Synthesis Kit (Almanac, Cat# ALS-MB72). The quantitative qRT-PCR was conducted consuming SYBR Green (Almanac, Cat# ALS-MB76) on an Applied Biosystems RT-PCR System. The 2-ΔΔCT method was working to calculate the relative fold change in gene expression. The primer sequence used for the amplifications was as: Aurora A protein kinase, forward: 5′-CCACCTTCGGCATCCTAATA-3′ and reverse: 5′-TCCAAGTGGTGCATATTCCA-3′ and GAPDH, forward: 5′-GAAGGTGAAGGTCGGAGTC-3′ and reverse: 5′-GAAGATGGTGATGGGATTTC-3′ (Wang et al., 2019).
RESULTS
Molecular docking and ADMET properties
The highest binding affinity was observed between Duartin and the Aurora A protein kinase. The binding affinity results are represented in Table 1 and Figures 1A-C. The interaction patterns of the top five hits are depicted in Figure 1C. Duartin formed several conventional hydrogen bonds with the LYS41 and ALA213 residues of Aurora A protein kinase. Additionally, Duartin formed various other interactions, as well as van der Waals, carbon hydrogen, Pi-sigma, alkyl, and Pi-alkyl bonds, with the LYS143, GLY142, GLY140, GLU260, THR217, GLY216, LEU139, PRO214, TYR212, VAL147, LEU263, LEU210, ALA160, LYS162 and ASP274 residues of Aurora A.
| Name of the ligand | Target protein | Binding Energy (kcal/mol) | pKi | Target protein's amino acids interacting with ligand molecules |
|---|---|---|---|---|
| Duartin | Aurora-A protein kinase. | Aurora-A protein kinase. | 5.94 | Conventional hydrogen bonds: LYS41, ALA213 van der Waals and other interactions: LYS143, GLY142, GLY140, GLU260, THR217, GLY216, LEU139, PRO214, TYR212, VAL147, LEU263, LEU210, ALA160, LYS162, ASP274. |
ADMET predictions, PASS analysis and Swiss Target Prediction
ADMET estimates, PASS examination, and SwissTargetPrediction collectively encourage the drug-like characteristic of Duartin. The ADMET properties of Duartin, presented in Table 2, demonstrate favourable characteristics in terms of absorption, distribution, metabolism, excretion, and lack of toxicity. Furthermore, the outcomes of PASS Analysis and SwissTargetPrediction analysis, represented in Table 3 and Figure 2 respectively, not only confirm but also elaborate on the diverse biological properties associated with Duartin.
| Properties | |||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Absorption | Distribution | Metabolism | Excretion | Toxicity | |||||||||
| Models | Intestinal absorption (human) | VDss (human) | BBB permeability | CNS permeability | CYP | Total clearance | AMES toxicity/Hepatoxicity | ||||||
| Substrate | Inhibitor | ||||||||||||
| 2D6 | 3A4 | 1A2 | 2C19 | 2C9 | 2D6 | 3A4 | |||||||
| Unity | Numeric (% absorbed) | Numeric (log L/kg) | Numeric (Log BB) | Numeric (Log PS) | Categorical (yes/no) | Numeric (log mL/ min / kg) | Categorical (yes/no) | ||||||
| Predicted values | |||||||||||||
| Duartin | 91.013 | 0.059 | -0.525 | -2.987 | NO | Yes | Yes | Yes | Yes | No | Yes | 0.084 | No/No |
| Ligand | Structure | Pa | Pi | Activity |
|---|---|---|---|---|
| Duartin | 0,756 | 0,012 | JAK2 expression inhibitor. | |
| 0,743 | 0,011 | Apoptosis agonist. | ||
| 0,551 | 0,015 | Anticarcinogenic. | ||
| 0,531 | 0,062 | Antineoplastic. | ||
| 0,375 | 0,027 | Antineoplastic (lung cancer). | ||
| 0,392 | 0,070 | Kinase inhibitor. | ||
| 0,316 | 0,201 | Antineoplastic (non-Hodgkin’s lymphoma). | ||
| 0,289 | 0,030 | Antineoplastic (brain cancer) | ||
| 0,284 | 0,027 | Antineoplastic (lymphoma). | ||
| 0,289 | 0,052 | MAP kinase kinase 4 inhibitor. | ||
| 0,223 | 0,025 | Antineoplastic (carcinoma). | ||
| 0,202 | 0,061 | Antineoplastic (glioblastoma multiforme). | ||
| 0,188 | 0,040 | Antineoplastic (glioma). | ||
| 0,180 | 0,088 | Antileukemic. | ||
| 0,175 | 0,071 | Antineoplastic (melanoma). | ||
| 0,142 | 0,002 | Antimitotic, Podophyllotoxin-like. | ||
| 0,108 | 0,106 | MAP-kinase-activated kinase inhibitor. |
Concentration-dependent cytotoxic effect of Duartin on breast cancer cells (MDA-MB-231)
The outcomes of the concentration-dependent cell toxicity assay are presented in Figure 3. This assay demonstrated a significant impact of Duartin on the viability of MDA-MB-231 cells. Additionally, the IC50 value for Duartin was determined to be 13.42 µM in MDA-MB-231 cancer cells.
Time dependent cytotoxic effect of Duartin on breast cancer cells
The outcomes of the time-dependent cell toxicity assay are depicted in Figure 4. Increased exposure durations of breast cancer cells to the IC50 concentration of Duartin significantly impacted their viability for up to 48 hr.
Expression of Aurora A protein kinase in Duartin treated breast cancer cells
The outcomes represented in Figure 5 illustration a statistically significant effect of the IC50 concentration of Duartin on the mRNA expression of Aurora A protein kinase in breast cancer cells. Specifically, mRNA expression was significantly reduced (p<0.001) in Duartin-treated cells compared to untreated cells.
DISCUSSION
Initially, the study aimed to investigate the binding affinity of various compounds with the Aurora A protein kinase, a crucial enzyme involved in cell cycle regulation and mitotic spindle assembly (Carmena et al., 2009). The results demonstrated that Duartin exhibited the highest binding affinity with Aurora A protein kinase, suggesting its potential as a strong inhibitor of this kinase. The binding affinity indicated that Duartin forms multiple stable interactions with key residues of Aurora A protein kinase. Specifically, Duartin establishes conventional hydrogen bonds with the LYS41 and ALA213 residues, which are critical for the stabilization of the ligand-protein complex (Panigrahi and Desiraju, 2007). The presence of these hydrogen bonds likely contributes significantly to the high binding affinity observed (Chen et al., 2016). Moreover, Duartin also engaged in a diversity of other non-covalent interactions, as well as van der Waals forces, carbon hydrogen bonds, Pi-sigma interactions, alkyl bonds, and Pi-alkyl bonds. These interactions further stabilize the Duartin-Aurora A complex, enhancing the overall binding affinity. The interaction pattern of Duartin with Aurora A protein kinase, emphasize the unique and robust interaction network formed by Duartin. The diverse array of bonds and interactions not only highlights the strong binding capability of Duartin but also highlights its specificity towards Aurora A protein kinase (Pitsawong et al., 2018; Zorba et al., 2014). This specificity is crucial for the development of targeted inhibitors that can effectively disrupt Aurora A protein kinase activity without affecting other kinases, thereby minimizing potential off-target effects (Zorba et al., 2014).
The findings of this study provide a comprehensive evaluation of the drug-like properties of Duartin, supported by ADMET calculations, PASS examination, and Swiss Target Prediction. The ADMET properties direct that Duartin keeps promising features in terms of absorption, distribution, metabolism, and excretion, while also exhibiting a lack of toxicity. These properties are crucial for any potential therapeutic agent, as they ensure that the compound can be efficiently absorbed, distributed to the target tissues, metabolized effectively and excreted without causing adverse effects (Ononamadu et al., 2021). The PASS analysis further corroborates the potential of Duartin as a therapeutic agent by predicting a range of biological activities. PASS analysis employs a robust algorithm to envisage the activity (biological) spectrum of a compound and the results for Duartin indicate its potential efficacy in various pharmacological domains (Poroikov et al., 2000). This predictive capability is essential for early-stage drug discovery, as it helps in identifying promising candidates for further investigation. SwissTargetPrediction analysis supports the findings of the PASS analysis by providing insights into the potential molecular targets of Duartin (Daina et al., 2019). Together; these computational predictions affirm the drug-like properties of Duartin and highlights its potential as a good therapeutic agent, targeting Aurora A protein kinase in breast cancer. The convergence of favourable ADMET characteristics, a broad spectrum of predicted biological activities and multiple potential molecular targets positions Duartin as a promising candidate for further preclinical and clinical studies.
The concentration-dependent cell toxicity assay revealed that Duartin exerts a pronounced effect on the viability of breast cancer cells. Duartin, 13.42 µM, indicates that Duartin has a potent cytotoxic effect on these cells. This finding is crucial as it proposes that Duartin could potentially be developed as a therapeutic agent for breast cancer, specifically targeting MDA-MB-231 cells, which are known for their aggressive behaviour and resistance to certain treatments. Furthermore, the time-dependent cell toxicity assay provided additional insights into the cytotoxic dynamics of Duartin. The outcomes presented that prolonged exposure of MDA-MB-231 cells to the Duartin significantly reduced cell viability over a 48 hr period. This time-dependent decrease in viability highlights the sustained action of Duartin, suggesting that continuous exposure could be necessary for achieving optimal therapeutic effects. The significant reduction in cell viability with both concentration and time indicates that Duartin acts effectively at low micromolar concentrations and continues to impact cell survival over extended periods (Hashemnia et al., 2016; Pancholi et al., 2022). These findings are promising for the possible expansion of Duartin as a therapeutic agent.
Further, the results of this study also highlight the critical influence of Duartin on the regulation of Aurora A protein kinase mRNA expression in breast cancer cells (MDA-MB-231). The marked decrease in Aurora A mRNA levels following treatment with Duartin, suggests a potent inhibitory effect of Duartin. This finding is particularly significant with the well-established role of Aurora A kinase in the regulation of mitotic entry and progression, as well as its association with poor prognosis in various cancers, including breast cancer (Du et al., 2021; Nikhil and Shah, 2024). Furthermore, the impact of Duartin on Aurora A protein kinase mRNA expression may provide insights into its broader antitumor mechanisms. Aurora A kinase is a key player in the regulation of mitosis and its overexpression has been linked to chromosomal instability and tumorigenesis (Katayama et al., 2023). By significantly reducing the expression of this kinase, Duartin may help to restore normal cell cycle control and reduce the propagation of cancer cells.
CONCLUSION
The study aimed to examine the binding affinity of various compounds with Aurora A protein kinase, a crucial enzyme in cell cycle regulation. Duartin exhibited the highest binding affinity, forming stable interactions with Aurora A protein kinase. ADMET properties, PASS and SwissTargetPrediction analyses affirm its drug-like properties. Duartin significantly reduced cell viability in breast cancer cells, portentous it’s probable as a effective therapeutic agent. Its ability to decrease Aurora A mRNA levels indicates a strong inhibitory effect, potentially restoring normal cell cycle control and reducing cancer cell proliferation. Overall, the findings from this study indicate that the Duartin may emerge as a good therapeutic agent for further preclinical and clinical studies targeting breast cancer.
