Potential inhibitors of VEGFR1, VEGFR2, and VEGFR3 developed through Deep Learning for the treatment of Cervical Cancer

Cervical cancer stands as a prevalent gynaecologic malignancy affecting women globally, often linked to persistent human papillomavirus infection. Biomarkers associated with cervical cancer, including VEGF-A, VEGF-B, VEGF-C, VEGF-D, and VEGF-E, show upregulation and are linked to angiogenesis and lymphangiogenesis. This research aims to employ in-silico methods to target tyrosine kinase receptor proteins-VEGFR-1, VEGFR-2, and VEGFR-3, and identify novel inhibitors for Vascular Endothelial Growth Factors receptors (VEGFRs). A comprehensive literary study was conducted which identified 26 established inhibitors for VEGFR-1, VEGFR-2, and VEGFR-3 receptor proteins. Compounds with high-affinity scores, including PubChem ID-25102847, 369976, and 208908 were chosen from pre-existing compounds for creating Deep Learning-based models. RD-Kit, a Deep learning algorithm, was used to generate 43 million compounds for VEGFR-1, VEGFR-2, and VEGFR-3 targets. Molecular docking studies were conducted on the top 10 molecules for each target to validate the receptor-ligand binding affinity. The results of Molecular Docking indicated that PubChem IDs-71465,645 and 11152946 exhibited strong affinity, designating them as the most efficient molecules. To further investigate their potential, a Molecular Dynamics Simulation was performed to assess conformational stability, and a pharmacophore analysis was also conducted for indoctrinating interactions.

Unveiling the ESR1 Conformational Stability and Screening Potent Inhibitors for Breast Cancer Treatment

BACKGROUND

The current study recognizes the significance of estrogen receptor alpha (ERα) as a member of the nuclear receptor protein family, which holds a central role in the pathophysiology of breast cancer. ERα serves as a valuable prognostic marker, with its established relevance in predicting disease outcomes and treatment responses.

METHODS

In this study, computational methods are utilized to search for suitable drug-like compounds that demonstrate analogous ligand binding kinetics to ERα.

RESULTS

Docking-based simulation screened out the top 5 compounds - ZINC13377936, NCI35753, ZINC35465238, ZINC14726791, and NCI663569 against the targeted protein. Further, their dynamics studies reveal that the compounds ZINC13377936 and NCI35753 exhibit the highest binding stability and affinity.

CONCLUSION

Anticipating the competitive inhibition of ERα protein expression in breast cancer, we envision that both ZINC13377936 and NCI35753 compounds hold substantial promise as potential therapeutic agents. These candidates warrant thorough consideration for rigorous In vitro and In vivo evaluations within the context of clinical trials. The findings from this current investigation carry significant implications for the advancement of future diagnostic and therapeutic approaches for breast cancer.

Experimental and Computational Approaches to Improve Binding Affinity in Chemical Biology and Drug Discovery

Drug discovery is one of the most complicated processes and establishment of a single drug may require multidisciplinary attempts to design efficient and commercially viable drugs. The main purpose of drug design is to identify a chemical compound or inhibitor that can bind to an active site of a specific cavity on a target protein. The traditional drug design methods involved various experimental based approaches including random screening of chemicals found in nature or can be synthesized directly in chemical laboratories. Except for the long cycle design and time, high cost is also the major issue of concern. Modernized computer-based algorithm including structure-based drug design has accelerated the drug design and discovery process adequately. Surprisingly from the past decade remarkable progress has been made concerned with all area of drug design and discovery. CADD (Computer Aided Drug Designing) based tools shorten the conventional cycle size and also generate chemically more stable and worthy compounds and hence reduce the drug discovery cost. This special edition of editorial comprises the combination of seven research and review articles set emphasis especially on the computational approaches along with the experimental approaches using a chemical synthesizing for the binding affinity in chemical biology and discovery as a salient used in de-novo drug designing. This set of articles exfoliates the role that systems biology and the evaluation of ligand affinity in drug design and discovery for the future.

Artificial Intelligence, Big Data and Machine Learning Approaches in Precision Medicine & Drug Discovery

Artificial Intelligence revolutionizes the drug development process that can quickly identify potential biologically active compounds from millions of candidate within a short period. The present review is an overview based on some applications of Machine Learning based tools, such as GOLD, Deep PVP, LIB SVM, etc. and the algorithms involved such as support vector machine (SVM), random forest (RF), decision tree and Artificial Neural Network (ANN), etc. at various stages of drug designing and development. These techniques can be employed in SNP discoveries, drug repurposing, ligand-based drug design (LBDD), Ligand-based Virtual Screening (LBVS) and Structure- based Virtual Screening (SBVS), Lead identification, quantitative structure-activity relationship (QSAR) modeling, and ADMET analysis. It is demonstrated that SVM exhibited better performance in indicating that the classification model will have great applications on human intestinal absorption (HIA) predictions. Successful cases have been reported which demonstrate the efficiency of SVM and RF models in identifying JFD00950 as a novel compound targeting against a colon cancer cell line, DLD-1, by inhibition of FEN1 cytotoxic and cleavage activity. Furthermore, a QSAR model was also used to predict flavonoid inhibitory effects on AR activity as a potent treatment for diabetes mellitus (DM), using ANN. Hence, in the era of big data, ML approaches have been evolved as a powerful and efficient way to deal with the huge amounts of generated data from modern drug discovery to model small-molecule drugs, gene biomarkers and identifying the novel drug targets for various diseases.

Most Promising Compounds for Treating COVID-19 and Recent Trends in Antimicrobial & Antifungal Agents

Multidrug resistance in microbes poses a major health crisis and demands for the discovery of novel antimicrobial agents. The recent pandemic of SARS-CoV-2 has raised a public health emergency in almost all the countries of the world. Unlike viruses, a bacterium plays a significant role in various environmental issues such as bioremediation. Furthermore, biosurfactants produced by various bacterial species have an edge over traditionally produced chemical surfactants for its biodegradability, low toxicity and better interfacial activity with various applications in agriculture and industry. This special issue focuses on the global perspective of drug discovery for various antimicrobial, antiviral, and antifungal agents for infectious diseases. The issue also emphasizes the ongoing developments and the role of microbes in environmental remediation. We wish the articles published in this issue will enhance the current understanding in microbiology among the readers, and serve as the "seed of an idea" for drug development for ongoing COVID-19 pandemic.

Shape-based Machine Learning Models for the Potential Novel COVID-19 Protease Inhibitors Assisted by Molecular Dynamics Simulation

BACKGROUND

The vast geographical expansion of novel coronavirus and an increasing number of COVID-19 affected cases have overwhelmed health and public health services. Artificial Intelligence (AI) and Machine Learning (ML) algorithms have extended their major role in tracking disease patterns, and in identifying possible treatments.

OBJECTIVE

This study aims to identify potential COVID-19 protease inhibitors through shape-based Machine Learning assisted by Molecular Docking and Molecular Dynamics simulations.

METHODS

31 Repurposed compounds have been selected targeting the main coronavirus protease (6LU7) and a machine learning approach was employed to generate shape-based molecules starting from the 3D shape to the pharmacophoric features of their seed compound. Ligand-Receptor Docking was performed with Optimized Potential for Liquid Simulations (OPLS) algorithms to identify highaffinity compounds from the list of selected candidates for 6LU7, which were subjected to Molecular Dynamic Simulations followed by ADMET studies and other analyses.

RESULTS

Shape-based Machine learning reported remdesivir, valrubicin, aprepitant, and fulvestrant as the best therapeutic agents with the highest affinity for the target protein. Among the best shape-based compounds, a novel compound identified was not indexed in any chemical databases (PubChem, Zinc, or ChEMBL). Hence, the novel compound was named 'nCorv-EMBS'. Further, toxicity analysis showed nCorv-EMBS to be suitable for further consideration as the main protease inhibitor in COVID-19.

CONCLUSION

Effective ACE-II, GAK, AAK1, and protease 3C blockers can serve as a novel therapeutic approach to block the binding and attachment of the main COVID-19 protease (PDB ID: 6LU7) to the host cell and thus inhibit the infection at AT2 receptors in the lung. The novel compound nCorv- EMBS herein proposed stands as a promising inhibitor to be evaluated further for COVID-19 treatment.

In silico Insights on IL-6: A Potential Target for Multicentric Castleman Disease

Background:

Multicentric Castleman Disease (MCD) is a confrontational lymphoproliferative disorder described by symptoms such as lymph node proliferation, unwarranted secretion of inflammatory cytokines, hyperactive immune system, and in severe cases, multiple organ dysfunction. Interleukin-6 (IL-6) is a pleiotropic cytokine which is involved in a large range of physiological processes in our body such as pro-inflammation, anti-inflammation, differentiation of T-cells and is reported to be a key pathological factor in MCD. In the case of MCD, it was observed that IL-6 is overproduced from T-cells and macrophages which disturb Hepcidin, a vital regulator of iron trafficking in macrophage. The present study endeavour to expound the inhibitor which binds to IL-6 protein receptor with high affinity.

Methods:

MolegroVirtual Docker software was employed to find the best-established drug from the list of selected inhibitors of IL-6. This compound was subjected to virtual screening against PubChem database to get inhibitors with a very similar structure. These inhibitors were docked to obtain a compound binding with high affinity to the target protein. The established compound and the virtual screened compound were subjected to relative analysis of interactivity energy variables and ADMET profile studies.

Results:

Among all the selected inhibitors, the virtual screened compound PubChem CID: 101119084 is seen to possess the highest affinity with the target protein. Comparative studies and ADMET analysis further implicate this compound as a better inhibitor of the IL-6 protein.

Conclusion:

Hence, this compound recognized in the study possesses high potential as an IL-6 inhibitor which might assist in the treatment of Multicentric Castleman Disease and should be examined for its efficiency by in vivo studies.

Computer-Aided Structure Prediction of Bluetongue Virus Coat Protein VP2 Assisted by Optimized Potential for Liquid Simulations (OPLS)

Background:

The capsid coated protein of Bluetongue virus (BTV) VP2 is responsible for BTV transmission by the Culicoides vector to vertebrate hosts. Besides, VP2 is responsible for BTV entry into permissive cells and hence plays a major role in disease progression. However, its mechanism of action is still unknown.

Objective:

The present investigation aimed to predict the 3D structure of Viral Protein 2 of the bluetongue virus assisted by Optimized Potential for Liquid Simulations (OPLS), structure validation, and an active site prediction.

Methods:

The 3D structure of the VP2 protein was built using a Python-based Computational algorithm. The templates were identified using Smith waterman’s Local alignment. The VP2 protein structure validated using PROCHECK. Molecular Dynamics Simulation (MDS) studies were performed using an academic software Desmond, Schrodinger dynamics, for determining the stability of a model protein. The Ligand-Binding site was predicted by structure comparison using homology search and proteinprotein network analysis to reveal their stability and inhibition mechanism, followed by the active site identification.

Results:

The secondary structure of the VP2 reveals that the protein contains 220 alpha helix atoms, 40 310 helix, 151 beta sheets, 134 coils and 424 turns, whereas the 3D structure of Viral Protein 2 of BTV has been found to have 15774 total atoms in the structure. However, 961 amino acids were found in the final model. The dynamical cross-correlation matrix (DCCM) analysis tool identifies putative protein domains and also confirms the stability of the predicted model and their dynamical behavior difference with the correlative fluctuations in motion.

Conclusion:

The biological interpretation of the Viral Protein 2 was carried out. DCCM maps were calculated, using a different coordinate reference frame, through which, protein domain boundaries and protein domain residue constituents were identified. The obtained model shows good reliability. Moreover, we anticipated that this research should play a promising role in the identification of novel candidates with the target protein to inhibit their functional significance.

In Silico Insights on GD2 : A Potential Target for Pediatric Neuroblastoma

Background:

Originating from the abnormal growth of neuroblasts, pediatric neuroblastoma affects the age group below 15 years. It is an aggressive heterogenous cancer with a high morbidity rate. Biological marker GD2 synthesised by the GD2 gene acts as a powerful predictor of neuroblastoma cells. GD2 gangliosides are sialic acid-containing glycosphingolipids. Differential expression during brain development governs the function of the GD2. The present study explains the interaction of the GD2 with its established inhibitors and discovers the compound having a high binding affinity against the target protein. Technically, during the development of new compounds through docking studies, the best drug among all pre-exist inhibitors was filtered. Hence in reference to the best docked compound, the study proceeded further.

Methodology:

The In silico approach provides a platform to determine and establish potential inhibitor against GD2 in Pediatric neuroblastoma. The 3D structure of GD2 protein was modelled by homology base fold methods using Smith-Watermans’ Local alignment. A total of 18 established potent compounds were subjected to molecular docking and Etoposide (CID: 36462) manifested the highest affinity. The similarity search presented 336 compounds similar to Etoposide.

Results:

Through virtual screening, the compound having PubChem ID 10254934 showed a better affinity towards GD2 than the established inhibitor. The comparative profiling of the two compounds based on various interactions such as H-bond interaction, aromatic interactions, electrostatic interactions and ADMET profiling and toxicity studies were performed using various computational tools.

Conclusion:

The docking separated the virtual screened drug (PubChemID: 10254934) from the established inhibitor with a better re-rank score of -136.33. The toxicity profile of the virtual screened drug was also lesser (less lethal) than the established drug. The virtual screened drug was observed to be bioavailable as it does not cross the blood-brain barrier. Conclusively, the virtual screened compound obtained in the present investigation is better than the established inhibitor and can be further augmented by In vitro analysis, pharmacodynamics and pharmacokinetic studies.

Identification of Potent VEGF Inhibitors for the Clinical Treatment of Glioblastoma, A Virtual Screening Approach

Vascular endothelial growth factor (VEGF) expression could be found in all glioblastomas. VEGF takes part in numerous changes including the endothelial cell proliferation, the vasculature of solid tumor: its survival invasion, and migration, chemotaxis of bone marrow-derived progenitor cells, vasodilation and vascular permeability. VEGF inhibition can be a smart therapeutic strategy because it is extremely specific and less toxic than cytotoxic therapy. To establish better inhibition of VEGF than the current inhibitors, present study approach is by molecular docking, virtual screening to illustrate the inhibitor with superior affinity against VEGF to have a cautious pharma profile. To retrieve the best established and high-affinity high affinity molecule, Molegro Virtual Docker software was executed. The high-affinity scoring compounds were subjected to further similarity search to retrieve the drugs with similar properties from pubchem database. The completion of virtual screening reveals that PubChem compound SCHEMBL1250485 (PubChem CID: 66965667) has the highest affinity. The study of the drug-likeness was verified using OSIRIS Property Explorer software which supported the virtual screened result. Further ADMET study and drug comparative study strongly prove the superiority of the new established inhibitor with lesser rerank score and toxicity. Overall, the new inhibitor has higher potential to stop the expression of VEGF in glioblastoma and positively can be further analysed through In vitro studies.

Identification of High-affinity Small Molecules Targeting Gamma Secretase for the Treatment of Alzheimer’s Disease

Background:

Alzheimers Disease (AD) is a neurodegenerative disease which is characterized by the deposition of amyloid plaques in the brain- a concept supported by most of the researchers worldwide. The main component of the plaques being amyloid-beta (Aβ42) results from the sequential cleavage of Amyloid precursor protein (APP) by beta and gamma secretase. This present study intends to inhibit the formation of amyloid plaques by blocking the action of gamma secretase protein with Inhibitors (GSI).

Methods:

A number of Gamma Secretase Inhibitors (GSI) were targeted to the protein by molecular docking. The inhibitor having the best affinity was used as a subject for further virtual screening methods to obtain similar compounds. The generated compounds were docked again at the same docking site on the protein to find a compound with higher affinity to inhibit the protein. The highlights of virtually screened compound consisted of Pharmacophore Mapping of the docking site. These steps were followed by comparative assessments for both the compounds, obtained from the two aforesaid docking studies, which included interaction energy descriptors, ADMET profiling and PreADMET evaluations.

Results:

111 GSI classified as azepines, sulfonamides and peptide isosteres were used in the study. By molecular docking an amorpholino-amide, compound (22), was identified to be the high affinity compound GSI along with its better interaction profiles.The virtually screened pubchem compound AKOS001083915 (CID:24462213) shows the best affinity with gamma secretase. Collective Pharmacophore mapping (H bonds, electrostatic profile, binding pattern and solvent accesibility) shows a stable interaction. The resulting ADMETand Descriptor values were nearly equivalent.

Conclusion:

These compounds identified herein hold a potential as Gamma Secretase inhibitors.According to PreADMET values the compound AKOS001083915 is effective and specific to the target protein. Its BOILED-egg plot analysis infers the compound permeable to blood brain barrier.Comparative study for both the compounds resulted in having nearly equivalent properties. These compounds have the capacity to inhibit the protein which is indirectly responsible for the formation of amyloid plaques and can be further put to in vitro pharmacokinetic and dynamic studies.

A Computer - Aided Drug Designing for Pharmacological Inhibition of Mutant ALK for the Treatment of Non-small Cell Lung Cancer

Background:

Lung cancer is the most common among all the types of cancer worldwide with 1.8 million people diagnosed every year, leading to 1.6 million deaths every year according to the American cancer society. The involvement of mutated Anaplasic Lymphoma Kinase (ALK) positive fusion protein in the progression of NSCLC has made a propitious target to inhibit and treat NSCLC. In the present study, the main motif is to screen the most effective inhibitor against ALK protein with the potential pharmacological profile. The ligands selected were docked with Molegro Virtual Docker (MVD) and CEP-37440 (PubChem CID- 71721648) was the best docked pre-established compound with a permissible pharmacological profile.

Methods:

The selected ligands were docked with Molegro Virtual Docker (MVD). With reference to the obtained compound with the lowest re-rank score, PubChem database was virtually screened to retrieve a large set of similar compounds which were docked to find the compound with higher affinity. Further comparative studies and in silico prediction included pharmacophore studies, proximity energy parameters, ADMET and BOILED-egg plot analysis.

Results:

CEP-37440 (PubChem CID- 71721648) was the best docked pre-established compound with preferable pharmacological profile and PubChem compound CID-123449015 came out as the most efficient virtually screened inhibitor. Interestingly, the contours of the virtual screened compound PubChem CID- 123449015 fall within our desired high volume cavity of protein having admirable property to control the ALK regulation to prevent carcinogenesis in NSCLC. BOILED-Egg plot analysis depicts that both the compounds have analogous characteristics in the divergent aspects. Moreover, in the evaluations of Blood Brain Barrier, Human Intestinal Absorption, AMES toxicity, and LD50, the virtually screened compound (PubChem CID-123449015) was found within high optimization.

Conclusion:

These investigations denote that the virtually screened compound (PubChem CID- 123449015) is more efficient to be a better prospective candidate for NSCLC treatment having good pharmacological profile than the pre-established compound CEP-37440 (PubChem CID- 71721648) with low re-rank score. The identified virtually screened compound has high potential to act as an ALK inhibitor and can show promising results in the research of non-small cell lung cancer (NSCLC).

Identification of High-Affinity Small Molecule Targeting IDH2 for the Clinical Treatment of Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is symbolized by an increase in the number of myeloid cells in the bone marrow and an arrest in their maturation, frequently resulting in hematopoietic insufficiency (granulocytopenia, thrombocytopenia, or anemia) with or without leukocytosis either by a predominance of immature forms or a loss of normal hematopoiesis. IDH2 gene encodes for isocitrate dehydrogenase enzyme which is involved in the TCA cycle domino effect and converts isocitrate to alpha-ketoglutarate. In the U.S, the annual incidence of AML progressively increases with age to a peak of 12.6 per 100,000 adults of 65 years or older. Mutations in isocitrate dehydrogenase 2 (arginine 132) have been demonstrated to be recurrent gene alterations in acute myeloid leukemia (AML) by forming 2-Hydroxy alpha ketoglutarate which, instead of participating in TCA cycle, accumulates to form AML. The current study approaches by molecular docking and virtual screening to elucidate inhibitor with superior affinity against IDH2 and achieve a pharmacological profile. To obtain the best established drug Molegro Virtual Docker algorithm was executed. The compound AG-221 (Pub CID 71299339) having the high affinity score was subjected to similarity search to retrieve the drugs with similar properties. The virtual screened compound SCHEMBL16391748 (PubChem CID-117816179)shows high affinity for the protein. Comparative study and ADMET study for both the above compounds resulted in equivalent chemical properties. Virtual screened compound SCHEMBL16391748 (PubChem CID-117816179) shows the lowest re-rank score. These drugs are identified as high potential IDH2 inhibitors and can halt AML when validated through further In vitro screening.

An In silico Approach to Identify High Affinity Small Molecule Targeting m-TOR Inhibitors for the Clinical Treatment of Breast Cancer

Breast cancer is the most frequent malignancy among women. It is a heterogeneous disease with different subtypes defined by its hormone receptor. A hormone receptor is mainly concerned with the progression of the PI3K/AKT/mTOR pathway which is often dysregulated in breast cancer. This is a major signaling pathway that controls the activities such as cell growth, cell division, and cell proliferation. The present study aims to suppress mTOR protein by its various inhibitors and to select one with the highest binding affinity to the receptor protein. Out of 40 inhibitors of mTOR against breast cancer, SF1126 was identified to have the best docking score of -8.705, using Schrodinger Suite which was further subjected for high throughput screening to obtain best similar compound using Lipinski’s filters. The compound obtained after virtual screening, ID: ZINC85569445 is seen to have the highest affinity with the target protein mTOR. The same result based on the binding free energy analysis using MM-GBSA showed that the compound ZINC85569445 to have the the highest binding free energy. The next study of interaction between the ligand and receptor protein with the pharmacophore mapping showed the best conjugates, and the ZINC85569445 can be further studied for future benefits of treatment of breast cancer.

Design of PD-L1 inhibitors for lung cancer

The progression of lung cancer is associated with inactivation of programmed cell death protein 1, abbreviated as PD- 1 which regulates the suppression of the body's immune system by suppressing T- cell inflammatory activity and is responsible for preventing cancer cell growth. It is of interest to identify inhibitors for PD-L1 dimeric structure through molecular docking and virtual screening. The virtual screened compound XGIQBUNWFCCMAS-UHFFFAOYSA-N (PubChem CID: 127263272) displays a high affinity with the target protein. ADMET analysis and cytotoxicity studies further add weight to this compound as a potential inhibitor of PD-L1. The established compound BMS-202 still shows the high re-rank score, but the virtual screened drug possesses a better ADMET profile with a higher intestinal absorption value and lower toxicity.

Design of novel JAK3 Inhibitors towards Rheumatoid Arthritis using molecular docking analysis

Multiple cytokines play a pivotal role in the pathogenesis of Rheumatoid Arthritis by inducing intracellular signaling and it is known that the members of the Janus kinase (JAK) family are essential for such signal transduction. Janus kinase 3 is a tyrosine kinase that belongs to the Janus family of kinases. Drugs targeting JAK3 in the treatment of Rheumatoid arthritis is relevant. Therefore, it is of interest to design suitable inhibitors for JAK3 dimer using molecular docking with Molegro Virtual Docker. The compound possessing the highest affinity score is subjected to virtual screening to retrieve inhibitors. The compound SCHEMBL19100243 (PubChem CID- 76749591) displays a high affinity with the target protein. The affinity scores of this compound are more than known drugs. ADMET analysis and BOILED Egg plot provide insights into this compound as a potent inhibitor of JAK3.

FLT3 inhibitor design using molecular docking based virtual screening for acute myeloid leukemia

Acute Myeloid Leukaemia (AML) is a blood cancer, which affects the red blood cells in the bone marrow. Of the possible proteins that are affected in AML, fms-like tyrosine kinase 3 (FLT3) has long been recognized as a potential therapeutic target as it affects the other signaling pathways and leads to a cascade of events. First-generation inhibitors sorafenib and midostaurin, as well as secondgeneration agents such as quizartinib and crenolanib are known. It is of interest to identify new compounds against FLT3 with improved activity using molecular docking and virtual screening. Molecular docking of existing inhibitors selected a top scoring bestestablished candidate Quizartinib having PubChem CID: 24889392. Similarity searching resulted in compound XGIQBUNWFCCMASUHFFFAOYSA-NPubChemCID: 44598530 which shows higher affinity scores. A comparative study of both the compounds using a drug-drug comparison, ADMET studies, boiled egg plot and pharmacophore parameters and properties confirmed the result and predicted the ligand to be an efficient inhibitor of FLT3.

Virtual Screening of IL-6 Inhibitors for Idiopathic Arthritis

Juvenile idiopathic arthritis (JIA) is a heterogeneous disease characterized by the arthritis of unknown origin and IL6 is a known target for JIA. 20 known inhibitors towards IL-6 were screened and Methotrexate (MTX) having PubChem ID: 126941 showed high binding capacity with the receptor IL-6. The similarity searching with this compound gave 269 virtual screened compounds. The said screening presented 269 possible drugs having structural similarity to Methotrexate. The docking studies of the screened drugs separated the compound having PubChem CID: 122677576 (re-rank value of -140.262). Toxicity and interaction profile validated this compound for having a better affinity with the target protein. Conclusively, this study shows that according to ADMET profile and BOILED-Egg plot, the compound (PubChem CID: 122677576) obtained from Virtual Screen could be the best drug in future during the prevention of juvenile idiopathic arthritis. In the current study, the drug CID: 122677576 is a potent candidate for treating JIA. The pharmacophore study revealed that the drug CID: 122677576 is a non-inhibitor of CYP450 microsomal enzymes and was found to be non-toxic, similar to the established drug Methotrexate (CID: 126941). It has a lower LD50 value of 2.6698mol/kg as compared to the established compound having LD50 value as 23.4955mol/kg. Moreover, the compound was found to be non-carcinogenic.