1
|
Thirunavukkarasu MK, Ramesh P, Karuppasamy R, Veerappapillai S. Transcriptome profiling and metabolic pathway analysis towards reliable biomarker discovery in early-stage lung cancer. J Appl Genet 2024:10.1007/s13353-024-00847-2. [PMID: 38443694 DOI: 10.1007/s13353-024-00847-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2023] [Revised: 02/21/2024] [Accepted: 02/22/2024] [Indexed: 03/07/2024]
Abstract
Earlier diagnosis of lung cancer is crucial for reducing mortality and morbidity in high-risk patients. Liquid biopsy is a critical technique for detecting the cancer earlier and tracking the treatment outcomes. However, noninvasive biomarkers are desperately needed due to the lack of therapeutic sensitivity and early-stage diagnosis. Therefore, we have utilized transcriptomic profiling of early-stage lung cancer patients to discover promising biomarkers and their associated metabolic functions. Initially, PCA highlights the diversity level of gene expression in three stages of lung cancer samples. We have identified two major clusters consisting of highly variant genes among the three stages. Further, a total of 7742, 6611, and 643 genes were identified as DGE for stages I-III respectively. Topological analysis of the protein-protein interaction network resulted in seven candidate biomarkers such as JUN, LYN, PTK2, UBC, HSP90AA1, TP53, and UBB cumulatively for the three stages of lung cancers. Gene enrichment and KEGG pathway analyses aid in the comprehension of pathway mechanisms and regulation of identified hub genes in lung cancer. Importantly, the medial survival rates up to ~ 70 months were identified for hub genes during the Kaplan-Meier survival analysis. Moreover, the hub genes displayed the significance of risk factors during gene expression analysis using TIMER2.0 analysis. Therefore, we have reason that these biomarkers may serve as a prospective targeting candidate with higher treatment efficacy in early-stage lung cancer patients.
Collapse
Affiliation(s)
| | - Priyanka Ramesh
- Bioinformatics Core, College of Agriculture, Agriculture Research and Graduate Education, Purdue University, West Lafayette, IN, 47907, USA
| | - Ramanathan Karuppasamy
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, 632014, Tamil Nadu, India
| | - Shanthi Veerappapillai
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, 632014, Tamil Nadu, India.
| |
Collapse
|
2
|
Thirunavukkarasu MK, Veerappapillai S, Karuppasamy R. Sequential virtual screening collaborated with machine-learning strategies for the discovery of precise medicine against non-small cell lung cancer. J Biomol Struct Dyn 2024; 42:615-628. [PMID: 36995235 DOI: 10.1080/07391102.2023.2194994] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2022] [Accepted: 03/17/2023] [Indexed: 03/31/2023]
Abstract
Dysregulation of MAPK pathway receptors are crucial in causing uncontrolled cell proliferation in many cancer types including non-small cell lung cancer. Due to the complications in targeting the upstream components, MEK is an appealing target to diminish this pathway activity. Hence, we have aimed to discover potent MEK inhibitors by integrating virtual screening and machine learning-based strategies. Preliminary screening was conducted on 11,808 compounds using the cavity-based pharmacophore model AADDRRR. Further, seven ML models were accessed to predict the MEK active compounds using six molecular representations. The LGB model with morgan2 fingerprints surpasses other models ensuing 0.92 accuracy and 0.83 MCC value versus test set and 0.85 accuracy and 0.70 MCC value with external set. Further, the binding ability of screened hits were examined using glide XP docking and prime-MM/GBSA calculations. Note that we have utilized three ML-based scoring functions to predict the various biological properties of the compounds. The two hit compounds such as DB06920 and DB08010 resulted excellent binding mechanism with acceptable toxicity properties against MEK. Further, 200 ns of MD simulation combined with MM-GBSA/PBSA calculations confirms that DB06920 may have stable binding conformations with MEK thus step forwarded to the experimental studies in the near future.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Muthu Kumar Thirunavukkarasu
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - Shanthi Veerappapillai
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - Ramanathan Karuppasamy
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| |
Collapse
|
3
|
Thirunavukkarasu MK, Veerappapillai S, Karuppasamy R. Computational biophysics approach towards the discovery of multi-kinase blockers for the management of MAPK pathway dysregulation. Mol Divers 2023; 27:2093-2110. [PMID: 36260173 DOI: 10.1007/s11030-022-10545-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2022] [Accepted: 10/06/2022] [Indexed: 10/24/2022]
Abstract
The MAPK pathway is important in human lung cancer and is improperly activated in a substantial proportion through number of ways. Strategies on dual-targeting RAF and MEK are an alternative option to diminish the limitations in this pathway inhibition. Hence, we implemented parallel pharmacophore screening of 11,808 DrugBank compounds against RAF and MEK. ADHRR and DHHRR were modeled as a pharmacophore hypothesis for RAF and MEK respectively. Importantly, these hypotheses resulted an AUC value of > 0.90 with the external data set. As a result of phase screening, glide docking, and prime-MM/GBSA scoring, it is determined that DB08424 and DB08907 have the best chances of acting as multi-kinase inhibitors. The pi-cation interaction with key amino acid residues of both target receptors may responsible for the stronger binding with these kinases. Cumulative 600 ns MD simulation studies validate the binding ability of these compounds. Significantly, the hit compounds resulted higher number of stable conformational state with less atomic movements than the reference compound against both targets. The anti-cancer efficacy of the lead compounds was validated through machine learning-based approaches. These findings suggest that DB08424 and DB08907 might be novel molecules to be explored further experimentally to block the MAPK signaling in lung cancer patients.
Collapse
Affiliation(s)
- Muthu Kumar Thirunavukkarasu
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, 632014, India
| | - Shanthi Veerappapillai
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, 632014, India
| | - Ramanathan Karuppasamy
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, 632014, India.
| |
Collapse
|
4
|
Thirunavukkarasu MK, Karuppasamy R. Drug repurposing combined with MM/PBSA based validation strategies towards MEK inhibitors screening. J Biomol Struct Dyn 2022; 40:12392-12403. [PMID: 34459701 DOI: 10.1080/07391102.2021.1970629] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Emergence of oncogenic mutations in the MAPK pathway gaining more impact in the recent years. Importantly, MEK is a core element of this pathway as it is easy to inhibit and is a gatekeeper of multiple malignancies. Therefore, we performed in-silico strategy to screen repurposed candidate for MEK protein using a library of 11,808 compounds from different clusters in the DrugBank database. Glide docking, Prime-MM/GBSA and QikProp analysis were implemented to retrieve the hits with high precision. The stability of the binding mode and binding affinity of the resultant hit were explored using molecular dynamic simulations and MM/PBSA approach. The results highlight that Nebivolol (DB04861) not only achieved a stable conformation in the MEK binding pocket but also displayed highest binding affinity than the other molecules investigated in our study. Taken together, we hypothesized that Nebivolol is an excellent candidate for the inhibition of MEK in NSCLC patients in future.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Muthu Kumar Thirunavukkarasu
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - Ramanathan Karuppasamy
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| |
Collapse
|
5
|
Thirunavukkarasu MK, Karuppasamy R. Forecasting determinants of recurrence in lung cancer patients exploiting various machine learning models. J Biopharm Stat 2022; 33:257-271. [PMID: 36397284 DOI: 10.1080/10543406.2022.2148162] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
Abstract
Lung cancer recurrence seems to be the most leading cause of death as well as deterioration of lifespan. Proper assessment of the probability of recurrence in early-stage lung cancer is necessary to push up the treatment progress. We therefore employed machine-learning technologies to forecast post-operative recurrence risks using 174 lung cancer patient records. Six classification algorithms logistic regression, SVM, decision tree classification, random forest classification, XGBoost and lightGBM were used to predict the cancer recurrence. The patient samples were divided into training and test group with the split ratio of 3:1 for model generation and the accuracy were validated using k-fold cross-validation method. It is worth noting that the logistic regression model outperformed all the models in both training (Accuracy = 0.82) and test set (Accuracy = 0.79) on k-fold validation. Further, the optimal features (n = 7) identified using the RFE method is certainly helpful to improve the model in a high precision. The imperative risk factors associated with recurrence were identified using three feature selection methods. Importantly, our research showed that age is an important prognostic factor to be considered during the recurrence prediction. Indeed, severe concern on the identified risk factors combined with predictive models assists the physician to reduce the cancer recurrence rate in patients with lung cancer.
Collapse
Affiliation(s)
- Muthu Kumar Thirunavukkarasu
- Department of Biotechnology, School of Bio Sciences and Technology Vellore Institute of Technology, Vellore Tamil Nadu, India
| | - Ramanathan Karuppasamy
- Department of Biotechnology, School of Bio Sciences and Technology Vellore Institute of Technology, Vellore Tamil Nadu, India
| |
Collapse
|
6
|
Thirunavukkarasu MK, Suriya U, Rungrotmongkol T, Karuppasamy R. In Silico Screening of Available Drugs Targeting Non-Small Cell Lung Cancer Targets: A Drug Repurposing Approach. Pharmaceutics 2021; 14:59. [PMID: 35056955 PMCID: PMC8778223 DOI: 10.3390/pharmaceutics14010059] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2021] [Revised: 12/16/2021] [Accepted: 12/20/2021] [Indexed: 01/03/2023] Open
Abstract
The RAS-RAF-MEK-ERK pathway plays a key role in malevolent cell progression in many tumors. The high structural complexity in the upstream kinases limits the treatment progress. Thus, MEK inhibition is a promising strategy since it is easy to inhibit and is a gatekeeper for the many malignant effects of its downstream effector. Even though MEK inhibitors are under investigation in many cancers, drug resistance continues to be the principal limiting factor to achieving cures in patients with cancer. Hence, we accomplished a high-throughput virtual screening to overcome this bottleneck by the discovery of dual-targeting therapy in cancer treatment. Here, a total of 11,808 DrugBank molecules were assessed through high-throughput virtual screening for their activity against MEK. Further, the Glide docking, MLSF and prime-MM/GBSA methods were implemented to extract the potential lead compounds from the database. Two compounds, DB012661 and DB07642, were outperformed in all the screening analyses. Further, the study results reveal that the lead compounds also have a significant binding capability with the co-target PIM1. Finally, the SIE-based free energy calculation reveals that the binding of compounds was majorly affected by the van der Waals interactions with MEK receptor. Overall, the in silico binding efficacy of these lead compounds against both MEK and PIM1 could be of significant therapeutic interest to overcome drug resistance in the near future.
Collapse
Affiliation(s)
- Muthu Kumar Thirunavukkarasu
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore 632014, India;
| | - Utid Suriya
- Program in Biotechnology, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand;
| | - Thanyada Rungrotmongkol
- Biocatalyst and Environmental Biotechnology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok 10330, Thailand
| | - Ramanathan Karuppasamy
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore 632014, India;
| |
Collapse
|
7
|
Thirunavukkarasu MK, Shin WH, Karuppasamy R. Exploring safe and potent bioactives for the treatment of non-small cell lung cancer. 3 Biotech 2021; 11:241. [PMID: 33968584 DOI: 10.1007/s13205-021-02797-6] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2020] [Accepted: 04/15/2021] [Indexed: 11/28/2022] Open
Abstract
Activating and suppressing mutations in the MAPK pathway receptors are the primary causes of NSCLC. Of note, MEK inhibition is considered a promising strategy because of the diverse structures and harmful effects of upstream receptors in MAPK pathway. Thus, we explore a total of 1574 plant-based bioactive compounds activity against MEK using an energy-based virtual screening strategy. Molecular docking, binding free energy, and drug-likeness analysis were performed through GLIDE, Prime MM-GBSA, and QikProp module, respectively. The findings indicate that 5-O-caffeoylshikimic acid has an increased binding affinity to MEK protein. Further, molecular dynamic simulations and MM-PBSA analysis were performed to explore the ligand activity in real-life situations. In essence, compounds inhibitory activity was validated across 77 lung cancer cell lines using multimodal attention-based neural network algorithm. Eventually, our analysis highlight that 5-O-caffeoylshikimic acid obtained from the bark of Rhizoma smilacis glabrae would be developed as a potential compound for treating NSCLC.
Collapse
Affiliation(s)
- Muthu Kumar Thirunavukkarasu
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Tamil Nadu, Vellore, 632014 India
| | - Woong-Hee Shin
- Department of Chemical Science Education, College of Education, Sunchon National University, Suncheon, Republic of Korea
| | - Ramanathan Karuppasamy
- Department of Biotechnology, School of Bio Sciences and Technology, Vellore Institute of Technology, Tamil Nadu, Vellore, 632014 India
| |
Collapse
|