- Department of Mathematics, University of Gujrat, Gujrat, Pakistan
The exploration of non-cancer medications with potential anti-cancer activity offers a promising avenue for drug repurposing, accelerating the development of new oncological therapies. This study employs Quantitative Structure-Property Relationship (QSPR) modeling to identify and predict the anti-cancer efficacy of various non-cancer drugs, utilizing topological indices as key descriptors. Topological indices, which capture the molecular structure’s geometric and topological characteristics, provide critical insights into the pharmacological interactions relevant to anti-cancer activity. By analyzing a comprehensive dataset of non-cancer medications, this research establishes robust QSPR models that correlate topological indices with anti-cancer activity. The models demonstrate significant predictive power, highlighting several non-cancer drugs with potential anti-cancer properties. Further, we will use linear, quadratic and logarithmic regression to understand the structures of anti-cancer drugs and strengthen our ability to manipulate the molecular structures. The findings underscore the utility of topological indices in drug repurposing strategies and pave the way for further experimental validation and clinical trials. This integrative approach enhances our understanding of drug action mechanisms and offers a cost-effective strategy for expanding the repertoire of anti-cancer agents.
1 Introduction
Cancer is a disease characterized by the uncontrolled growth and spread of certain body cells to other regions. It can originate in nearly any part of the human body, which consists of trillions of cells. Under normal circumstances, human cells grow and divide through cell division to create new cells as needed. When cells become old or damaged, they die and are replaced by new cells (Trichopoulos et al., 1996). Several symptoms and indications of this disease include weight loss, lumps, irregular bleeding, and prolonged coughing. Chewing tobacco, obesity, poor food, laziness, and increased alcohol consumption (Cleeland, 2000) are the main causes of this cancerous illness. To know more about this disease explore (Bailar and Gornik, 1997).
The rapid advancement of cancer research continually emphasizes the urgent need for new and effective oncological therapies. Traditional drug development processes are often time-consuming and expensive, prompting the scientific community to explore alternative strategies, such as drug repurposing. Drug repurposing involves identifying new therapeutic uses for existing medications, offering a cost-effective and expedited pathway to discovering anti-cancer agents (Figuerola and Avila, 2019; Gao et al., 2016; Kumar et al., 2015). This approach leverages the established safety profiles and pharmacokinetic properties of non-cancer drugs, significantly reducing the time and resources required for drug development.
In this context, Quantitative Structure-Property Relationship (QSPR) modeling emerges as a powerful tool. QSPR models utilize mathematical relationships to correlate chemical structure with biological activity, aiding in the prediction of a compound’s therapeutic potential. Specifically, topological indices, which are numerical representations capturing the geometric and topological characteristics of molecular structures, play a crucial role in these models. Topological indices have five main types: matching, mixed, eigenvalue, degree, and distance (Shanmukha et al., 2020). Degree-based topological indices on anti-cancer drugs are presented in this paper. Topological indices provide valuable insights into the molecular features that influence drug activity, making them indispensable in the identification of potential anti-cancer properties of non-cancer medications. For more details of QSPR analysis, see (Ren, 1999; Havare, 2021; Zhong et al., 2021; Adnan et al., 2022; Arockiaraj et al., 2023a; Arockiaraj et al., 2023b; Arockiaraj et al., 2023c; Zaman et al., 2023).
In topological analysis, stereochemistry refers to the study of the spatial arrangements of atoms within a molecule without considering their specific three dimensional orientation. Topological analysis techniques allow chemists to understand the connectivity of atoms and the overall shape of molecules which is crucial for predicting their chemical behavior and reactivity. Topological analysis focuses on understanding how atoms are connected within a molecule. Different types of bonds (single, double and triple) and their connectivity provide information about the stereochemistry of the molecules. For example, the presence of double bonds can lead to geometric isomerism, where the spatial arrangement of substituents around the double bond affects the molecule’s properties. By analyzing the connectivity of atom and their realative positions, topological analysis can distinguish between different stereoisomeric forms of a molecule. Overall, topological analysis provides valuable insight into the stereochemistry of molecules by focusing on their connectivity, symmetric, and spatial arrangement. By considering these factors, chmeists can predict the behavior of molecules in various chemical reactions and design new compounds with desired properties. Thus, Figure 1 plays a pivotal role in elucidating the stereochemical aspects of the molecules under consideration.
This study focuses on the application of QSPR modeling, with an emphasis on topological indices, to predict and identify non-cancer drugs that may exhibit anti-cancer activity. By analyzing a comprehensive dataset of non-cancer medications, we aim to establish robust QSPR models that accurately correlate topological indices with anti-cancer efficacy. The predictive power of these models is expected to highlight several non-cancer drugs as promising candidates for further experimental validation and clinical trials. Some of the common anti-cancer drugs approved by the Food and Drug Administration (FDA) and their molecular targets are shown in Table 1. Explore (Gupta et al., 2013; D Amato et al., 1994; Gasic et al., 1972; Michaelis et al., 2007; Jendrossek, 2013; Jiang et al., 2018; Wang et al., 2010; Zekri et al., 2014; Rezaei et al., 2023) to learn about these anti-cancer drugs. Drugs were selected based on a meticulous assessment encompassing structural characteristics, predicted anti-cancer activity, and suitability for repurposing. Leveraging topological indices and quantitative structure-property relationship (QSPR) modeling, we prioritized molecules with favorable molecular connectivity patterns indicative of potential anti-cancer properties. Each selected drug underwent scrutiny for safety profiles, availability, known mechanisms of action, and previous clinical data, ensuring a comprehensive approach towards identifying promising candidates for repurposing. This methodical selection process aimed to maximize the likelihood of uncovering non-cancer medications with significant potential for anti-cancer therapeutic applications (Singhal et al., 2022). The implications of this research are significant. Not only does it enhance our understanding of the molecular mechanisms underlying drug action, but it also paves the way for more efficient and cost-effective drug discovery processes. By integrating computational modeling with drug repurposing strategies, this study contributes to expanding the repertoire of available anti-cancer agents, ultimately improving therapeutic outcomes for cancer patients.
Consider G (V,E) a molecular graph with vertex and edge sets denoted by V and E, respectively. The number of vertices adjacent to a vertex is known as the degree of a vertex which is denoted as deg(v), while d (u,v) is the shortest distance between two vertices. In this work, we will use following topological indices.
1.1 Difference between ABC and R index ABC-R(G)
In present times, there is a growing interest in studying the correlation or comparison between topological indices, see (Das et al., 2016). Recently, Ali and Du (Ali and Du, 2017) explored extremal binary and chemical trees, specifically focusing on the difference between ABC and R indices. The ABC-R index provides additional information about the molecular structure beyond just the topological complexity measured by the ABC index alone. It accounts for the size of the molecule and provides insights into its topological properties relative to its size. The ABC-R index is defined as
1.2 Geometric arithmetic index GA(G)
Another recently conceived vertex-degree-based topological index utilizes the difference between the geometric and arithmetic means and is defined as
Where, of course
1.3 Multiplicative first exponential Zagreb index (G)
The multiplicative first exponential Zagreb index (Akgunes and Aydin, 2021) is a molecular descriptor in chemical graph theory. This index provides valuable structural information about the molecular graph, with a focus on the significance of highly connected vertices. This index is used in quantitative structure-property relationship (QSPR) studies to predict various chemical and physical properties of molecules. These properties may include boiling points, solubility, stability, and other characteristics relevant in chemistry and pharmacology.
1.4 Multiplicative second exponential Zagreb index (G)
The multiplicative second exponential Zagreb index (Akgunes and Aydin, 2021) provides valuable structural information about the molecular graph, focusing on the significance of pairs of highly connected vertices. It is used to predict various chemical and physical properties of molecules, particularly emphasizing the role of vertex degrees and their pairwise interactions in determining molecular behavior.
1.5 Multiplicative geometric arithmetic index GA (G)
Multiplicative GA index (Kulli, 2016) is a molecular descriptor that is calculated by takin the product of edge multiplicities and the geometric mean of vertex degrees in a molecular graph. This index provides information about the atom-bond connectivity in the molecular graph, considering both the geometric and arithmetic means of the degrees of the connected vertices.
1.6 Symmetric division degree index SDD(G)
Several years ago, D. Vuki
2 Methodology
The methodology for this study involves a systematic approach to identify and predict the anti-cancer potential of non-cancer medications using Quantitative Structure-Property Relationship (QSPR) modeling with topological indices. The process begins with data collection, where a comprehensive dataset of non-cancer medications is compiled from publicly available databases such as PubChem and ChemSpider. This dataset includes detailed molecular structures and known pharmacological properties of each drug. Next, molecular descriptor calculation is performed, where various topological indices are calculated for each drug in the dataset. These indices include the symmetric division degree index, geometric arithmetic index, multiplicative first exponential Zagreb index, difference between ABC and R index, multiplicative second exponential Zagreb index, and multiplicative geometric arithmetic index.
In the model construction phase, multiple QSPR models are developed using linear, quadratic, and logarithmic regression algorithms with the aid of SPSS software. Each type of regression model is built to explore different potential relationships between the topological indices and anti-cancer activity. Linear regression models assume a direct proportional relationship, quadratic regression models consider a parabolic relationship, and logarithmic regression models account for situations where the change in anti-cancer activity diminishes as the value of the topological indices increases. These models are trained on the selected topological indices and known anti-cancer activity data.
Finally, data analysis and interpretation are performed. The experimental data are analyzed using statistical methods to compare the predicted and observed anti-cancer activities. Correlation coefficients and significance tests are used to assess the accuracy and reliability of the QSPR models. Based on the experimental validation results, the QSPR models are refined and retrained to improve their predictive power and reliability. By following this structured methodology, the study aims to establish a reliable framework for repurposing non-cancer medications as potential anti-cancer agents, leveraging the power of QSPR modeling and topological indices. Figure 2 shows the basic flowchart of this work.
3 Analyzing molecular structures and computing topological indices for various anti-cancer medications
Topological indices are used to analyze the different structures. In this article, six TIs are used to describe the structures of fifteen anti-cancer drugs. The structures of these drugs are shown in Figure 3. The six physio-chemical properties are obtained from ChemSpider which are given in Table 2 and the value of each topological index corresponding to each structure is computed with the help of their formulas, as shown in Table 3.
Table 2. The physio-chemical properties of non-cancer medication to be used in the treatment of cancer.
Table 3. The values of topological indices of molecular structures of non-cancer medication to be used in the treatment of cancer.
The main purpose of this section is to give the overview of data. Derived from glutamic acid, thalidomide was initially formulated in the 1950s as a sedative-hypnotic to address pregnancy-related nausea. Recent research spanning the past decade suggests that thalidomide, initially explored for its potential anti-angiogenic properties, has the ability to influence various cell signaling pathways associated with cancer. Suppose G be the molecular structure of Thalidomide which is shown in Figure 3. The molecular graph of G has 19 vertices and 21 edges. It has four edge partition which is shown in Table 4 By using the same methodology we can determine the other topological indices. Table 1 shows the original indication and new anti-cancer indications of these drugs.
4 Regression models
The following equations are used to correlate the various physical properties of various non-cancer medications used for the treatment of cancer with some topological indices. We have used the following linear, quadratic and logarithmic regression models:
In the above Equations 1–3, p is dependent variable and TI is the independent variable. Where P is physical property of drug, A is constant, b and c are the regression coefficients. The six physiochemical properties: boiling point (BP) in
4.1 Difference between ABC and R index ABC-R(G)
4.2 Multiplicative first exponential Zagreb index
4.3 Multiplicative second exponential Zagreb index
4.4 Geometric arithmetic index GA(G)
4.5 Multiplicative geometric arithmetic index
4.6 Symmetric division degree index SDD(G)
5 Results and discussions
In this section, we delve into the statistical analysis of our regression models, which aimed to predict the anti-cancer properties of non-cancer medications based on their molecular characteristics represented by topological indices. The regression parameters for linear, quadratic, and logarithmic models were computed, providing insights into the relationships between the independent variables (topological indices) and the dependent variables (anti-cancer properties) The key parameters analyzed include the sample size (N), constant or Y-intercept (A), coefficients of the independent variables (b and c), correlation coefficient (r), and the percentage of variation explained by the linear model
The correlation coefficient (r) indicates the strength and direction of the relationship between variables, with values ranging from −1 to +1. A positive coefficient signifies a direct relationship, while a negative coefficient suggests an inverse relationship. The high correlation coefficients observed in our analysis indicate strong associations between the topological indices and anti-cancer properties, supporting the validity of our models. Moreover, the p-values associated with each term in the regression models were examined to test the null hypothesis that the coefficient is zero, implying no effect. A smaller p-value suggests that changes in the predictor variables are associated with changes in the response variable, indicating the significance of the predictors. In our analysis, all p-values were found to be zero, indicating the statistical significance of the regression models. Furthermore, we conducted additional tests to assess the overall predictive capability of our models. The F-value resulting from these tests helps to determine whether the model has predictive power beyond chance. The significant F-values obtained in our analysis confirm the predictive capability of our regression models in accurately predicting anti-cancer properties based on topological indices. These statistical characteristics for logarithmic, linear, and quadratic QSPR models for various topological indices are presented in Tables 5–10. Overall, the statistical analysis presented in this section underscores the robustness and reliability of our regression models in identifying non-cancer medications with potential anti-cancer properties. These findings provide a solid foundation for further exploration and validation of drug repurposing strategies in cancer therapy.
The study emphasizes the importance of topological descriptors in predicting the molecular structures of anti-cancer drugs, crucial for estimating the molecular weight of non-cancer medications. Quantitative Structure-Property Relationship (QSPR) methodology stands out as a robust approach utilized across diverse fields such as drug design, material science, environmental chemistry, cheminformatics, and computational chemistry. QSPR methodology is focused on establishing mathematical relationships between the chemical structure of compounds and their properties, enabling the prediction of properties based on molecular characteristics. Unlike QSAR, which typically addresses biological activities, QSPR specifically targets the physical and chemical properties of compounds. Molecular modeling techniques often entail complex simulations and calculations to forecast molecular behavior accurately. Ultimately, QSPR methodology provides a systematic and quantitative means to predict physicochemical properties of chemical compounds based on their molecular structure, distinguishing it from other methodologies like QSAR, MD simulation, DFT, machine learning models, and hybrid QSAR/QSPR models, each offering its unique advantages and limitations depending on the specific research objectives and applications.
6 Conclusion
This study demonstrates the efficacy of utilizing Quantitative Structure-Property Relationship (QSPR) modeling, specifically leveraging topological indices, to identify non-cancer medications with potential anti-cancer properties. The systematic approach outlined ranging from data collection and molecular descriptor calculation to model construction and experimental validation has proven to be a robust framework for drug repurposing. By focusing on topological indices, we have highlighted the significant role that molecular structure plays in determining pharmacological interactions relevant to anti-cancer activity.
Our findings indicate that several non-cancer drugs, identified through our QSPR models, exhibit promising anti-cancer properties, warranting further experimental validation and clinical trials. The use of topological indices in the modeling process has provided critical insights into the structural attributes that contribute to anti-cancer efficacy, thus enhancing our understanding of drug action mechanisms.
The integration of computational modeling with experimental validation offers a cost-effective and accelerated pathway for expanding the repertoire of anti-cancer agents. This research not only underscores the potential of drug repurposing strategies in oncology but also establishes a foundation for future studies to explore and refine the application of QSPR models. Ultimately, this approach holds the promise of improving therapeutic outcomes for cancer patients by identifying new uses for existing medications, thereby bridging the gap between drug discovery and clinical application.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.
Author contributions
SY: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing–original draft, Writing–review and editing. KS: Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Adnan, M., Bokhary, S. A. U. H., Abbas, G., and Iqbal, T. (2022). Degree-based topological indices and QSPR analysis of antituberculosis drugs. J. Chem. 2022, 1–17. doi:10.1155/2022/5748626
Akgunes, N., and Aydin, B. (2021). Introducing new exponential Zagreb indices for graphs. J. Math. 2021, 1–13. doi:10.1155/2021/6675321
Ali, A., and Du, Z. (2017). On the difference between atom-bond connectivity index and Randic index of binary and chemical trees. Int. J. Quantum Chem. 117 (23), e25446. doi:10.1002/qua.25446
Arockiaraj, M., Greeni, A. B., and Kalaam, A. A. (2023b). Comparative analysis of reverse degree and entropy topological indices for drug molecules in blood cancer treatment through QSPR regression models. Polycycl. Aromat. Compd., 1–18. doi:10.1080/10406638.2023.2271648
Arockiaraj, M., Greeni, A. B., and Kalaam, A. A. (2023c). Linear versus cubic regression models for analyzing generalized reverse degree based topological indices of certain latest corona treatment drug molecules. Int. J. Quantum Chem. 123 (16), e27136. doi:10.1002/qua.27136
Arockiaraj, M., Paul, D., Clement, J., Tigga, S., Jacob, K., and Balasubramanian, K. (2023a). Novel molecular hybrid geometric-harmonic-Zagreb degree based descriptors and their efficacy in QSPR studies of polycyclic aromatic hydrocarbons. SAR QSAR Environ. Res. 34 (7), 569–589. doi:10.1080/1062936x.2023.2239149
Bailar, J. C., and Gornik, H. L. (1997). Cancer undefeated. N. Engl. J. Med. 336 (22), 1569–1574. doi:10.1056/nejm199705293362206
Cleeland, C. S. (2000). Cancer-related symptoms. Seminars Radiat. Oncol. 10 (3), 175–190. doi:10.1053/srao.2000.6590
D Amato, R. J., Loughnan, M. S., Flynn, E., and Folkman, J. (1994). Thalidomide is an inhibitor of angiogenesis. Proc. Natl. Acad. Sci. 91 (9), 4082–4085. doi:10.1073/pnas.91.9.4082
Das, K. C., Das, S., and Zhou, B. (2016). Sum-connectivity index of a graph. Front. Math. china 11, 47–54. doi:10.1007/s11464-015-0470-2
Du, Z., Zhou, B., and Trinajstić, N. (2011). On geometric-arithmetic indices of (molecular) trees, unicyclic graphs and bicyclic graphs. MATCH Commun. Math. Comput. Chem. 66 (2), 681–697.
Figuerola, B., and Avila, C. (2019). The phylum bryozoa as a promising source of anticancer drugs. Mar. drugs 17 (8), 477. doi:10.3390/md17080477
Gao, W., Wang, W., and Farahani, M. R. (2016). Topological indices study of molecular structure in anticancer drugs. J. Chem. 2016, 1–8. doi:10.1155/2016/3216327
Gasic, G., Gasic, T., and Murphy, S. (1972). Anti-metastatic effect of aspirin. Lancet 300 (7783), 932–933. doi:10.1016/s0140-6736(72)92581-0
Gupta, S. C., Sung, B., Prasad, S., Webb, L. J., and Aggarwal, B. B. (2013). Cancer drug discovery by repurposing teaching new tricks to old dogs. Trends Pharmacol. Sci. 34 (9), 508–517. doi:10.1016/j.tips.2013.06.005
Havare, Ö. Ç. (2021). Topological indices and QSPR modeling of some novel drugs used in the cancer treatment. Int. J. Quantum Chem. 121 (24), e26813. doi:10.1002/qua.26813
Jendrossek, V. (2013). Targeting apoptosis pathways by Celecoxib in cancer. Cancer Lett. 332 (2), 313–324. doi:10.1016/j.canlet.2011.01.012
Jiang, L., Zhang, W., Li, W., Ling, C., and Jiang, M. (2018). Anti-inflammatory drug, leflunomide and its metabolite teriflunomide inhibit NSCLC proliferation in vivo and in vitro. Toxicol. Lett. 282, 154–165. doi:10.1016/j.toxlet.2017.10.013
Kulli, V. R. (2016). Multiplicative connectivity indices of certain nanotubes. Ann. Pure Appl. Math. 12 (2), 169–176. doi:10.22457/apam.v12n2a8
Kumar, S., Ahmad, M. K., Waseem, M., and Pandey, A. K. (2015). Drug targets for cancer treatment: an overview. Med. Chem. 5 (3), 115–123. doi:10.4172/2161-0444.1000252
Michaelis, M., Doerr, H. W., and Cinatl, J. (2007). Valproic acid as anti-cancer drug. Curr. Pharm. Des. 13 (33), 3378–3393. doi:10.2174/138161207782360528
Ren, B. (1999). A new topological index for QSPR of alkanes. J. Chem. Inf. Comput. Sci. 39 (1), 139–143. doi:10.1021/ci980098p
Rezaei, A., Moqadami, A., and Khalaj-Kondori, M. (2023). Minocycline as a prospective therapeutic agent for cancer and non-cancer diseases: a scoping review. Naunyn-Schmiedeberg’s Archives Pharmacol. 397, 2835–2848. doi:10.1007/s00210-023-02839-1
Shanmukha, M. C., Basavarajappa, N. S., Shilpa, K. C., and Usha, A. (2020). Degree-based topological indices on anticancer drugs with QSPR analysis. Heliyon 6 (6), e04235. doi:10.1016/j.heliyon.2020.e04235
Singhal, S., Maheshwari, P., Krishnamurthy, P. T., and Patil, V. M. (2022). Drug repurposing strategies for non-cancer to cancer therapeutics. Anti-Cancer Agents Med. Chem. Former.( Curr. Med. Chemistry-Anti-Cancer Agents) 22 (15), 2726–2756. doi:10.2174/1871520622666220317140557
Trichopoulos, D., Li, F. P., and Hunter, D. J. (1996). What causes cancer? Sci. Am. 275 (3), 80–87. doi:10.1038/scientificamerican0996-80
Vukicevic, D., and Furtula, B. (2009). Topological index based on the ratios of geometrical and arithmetical means of end-vertex degrees of edges. J. Math. Chem. 46, 1369–1376. doi:10.1007/s10910-009-9520-x
Vukičević, D., and Gašperov, M. (2010). Bond additive modeling 1. Adriatic indices. Croat. Chem. acta 83 (3), 243.
Wang, X., Wu, Q., Zhang, L., Wu, Y., and Shu, Y. (2010). Wortmannin induced apoptosis of leukemia cells by reducing PI3K/Akt. Chinese-German J. Clin. Oncol. 9, 734–738. doi:10.1007/s10330-010-0715-1
Zaman, S., Yaqoob, H. S. A., Ullah, A., and Sheikh, M. (2023). QSPR analysis of some novel drugs used in blood cancer treatment via degree based topological indices and regression models. Polycycl. Aromat. Compd. 44, 2458–2474. doi:10.1080/10406638.2023.2217990
Zekri, J., Mansour, M., and Karim, S. M. (2014). The anti-tumour effects of zoledronic acid. J. bone Oncol. 3 (1), 25–35. doi:10.1016/j.jbo.2013.12.001
Keywords: topological indices, QSPR models, regression model, anti-cancer activity, non-cancer medications
Citation: Yousaf S and Shahzadi K (2024) Utilizing topological indices in QSPR modeling to identify non-cancer medications with potential anti-cancer properties: a promising strategy for drug repurposing. Front. Chem. 12:1410882. doi: 10.3389/fchem.2024.1410882
Received: 02 April 2024; Accepted: 11 July 2024;
Published: 08 August 2024.
Edited by:
Sam P. De Visser, The University of Manchester, United KingdomReviewed by:
María Guadalupe Hernández-Linares, Benemérita Universidad Autónoma de Puebla, MexicoSevki Adem, Cankiri Karatekin University, Türkiye
Copyright © 2024 Yousaf and Shahzadi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shamaila Yousaf, shumaila.yousaf@uog.edu.pk