Docking-based Virtual Screening Research Articles

Evaluation of AlphaFold2 Structures for Hit Identification across Multiple Scenarios.

The introduction of AlphaFold2 (AF2) has sparked significant enthusiasm and generated extensive discussion within the scientific community, particularly among drug discovery researchers. Although previous studies have addressed the performance of AF2 structures in virtual screening (VS), a more comprehensive investigation is still necessary considering the paramount importance of structural accuracy in drug design. In this study, we evaluate the performance of AF2 structures in VS across three common drug discovery scenarios: targets with holo, apo, and AF2 structures; targets with only apo and AF2 structures; and targets exclusively with AF2 structures. We utilized both the traditional physics-based Glide and the deep-learning-based scoring function RTMscore to rank the compounds in the DUD-E, DEKOIS 2.0, and DECOY data sets. The results demonstrate that, overall, the performance of VS on AF2 structures is comparable to that on apo structures but notably inferior to that on holo structures across diverse scenarios. Moreover, when a target has solely AF2 structure, selecting the holo structure of the target from different subtypes within the same protein family produces comparable results with the AF2 structure for VS on the data set of the AF2 structures, and significantly better results than the AF2 structures on its own data set. This indicates that utilizing AF2 structures for docking-based VS may not yield most satisfactory outcomes, even when solely AF2 structures are available. Moreover, we rule out the possibility that the variations in VS performance between the binding pockets of AF2 and holo structures arise from the differences in their biological assembly composition.

Design, synthesis and biological activity evaluation of novel inhibitors targeting TNKS2, based on docking, virtual screening and molecular dynamics simulation studies

In recent years, tankyrases (TNKS) have garnered substantial interest as a promising new drug target. An escalating amount of research has been dedicated to the development of TNKS-targeting drugs, highlighting the substantial potential in anti-cancer therapeutics. In this study, a potential TNKS2 lead compound, ZINC11726230, was identified through docking-based virtual screening, molecular dynamics simulation and bioassay. Structural modification of the lead compound, ZINC11726230, yielded two new TNKS2 inhibitors (compound 1-2), which exhibited remarkable inhibitory activity, with an IC50 value of 0.013 μM, outperforming the positive control compound G007-LK (IC50 = 0.035 μM). This study offers unique insights into the discovery of effective TNKS2 inhibitors.

Discovery of Novel Anti-Resistance AR Antagonists Guided by Funnel Metadynamics Simulation.

Androgen receptor (AR) antagonists are widely used for the treatment of prostate cancer (PCa), but their therapeutic efficacy is usually compromised by the rapid emergence of drug resistance. However, the lack of the detailed interaction between AR and its antagonists poses a major obstacle to the design of novel AR antagonists. Here, funnel metadynamics is employed to elucidate the inherent regulation mechanisms of three AR antagonists (hydroxyflutamide, enzalutamide, and darolutamide) on AR. For the first time it is observed that the binding of antagonists significantly disturbed the C-terminus of AR helix-11, thereby disrupting the specific internal hydrophobic contacts of AR-LBD and correspondingly the communication between AR ligand binding pocket (AR-LBP), activation function 2 (AF2), and binding function 3 (BF3). The subsequent bioassays verified the necessity of the hydrophobic contacts for AR function. Furthermore, it is found that darolutamide, a newly approved AR antagonist capable of fighting almost all reported drug resistant AR mutants, can induce antagonistic binding structure. Subsequently, docking-based virtual screening toward the dominant binding conformation of AR for darolutamide is conducted, and three novel AR antagonists with favorable binding affinity and strong capability to combat drug resistance are identified by in vitro bioassays. This work provides a novel rational strategy for the development of anti-resistant AR antagonists.

Identification of 3-aryl-5-methyl-isoxazole-4-carboxamide derivatives and analogs as novel HIF-2α agonists through docking-based virtual screening and structural modification

Hypoxia-inducible factor-2 (HIF-2) serves as the pivotal transcription factor in cellular responses to low oxygen levels, particularly concerning the regulation of erythropoietin (EPO) production. A docking-based virtual screening on crystal structures of HIF-2α inhibitors unexpectedly identified 3-phenyl-5-methyl-isoxazole-4-carboxamide derivative v19 as a hit of HIF-2α agonist. Further structural optimizations of compound v19 led to the discovery of a series of HIF-2α agonists with novel scaffolds. The most promising compounds 12g and 14d exhibited potent HIF-2α agonistic activities in vitro with EC50 values of 2.29 μM and 1.78 μM, respectively. Molecular dynamics simulations have revealed their capacity to allosterically enhance HIF-2 dimerization, which shed light on their mechanism of action. Moreover, compound 14d demonstrated a favorable pharmacokinetic (PK) profile, boasting an impressive oral bioavailability value of 68.71 %. These findings strongly suggest that compound 14d is an auspicious lead compound for the treatment of renal anemia.

Sterically hindered phenolic derivatives: effect on the production of Pseudomonas aeruginosa virulence factors, high-throughput virtual screening and ADME properties prediction.

Inhibition of quorum sensing is considered to be an effective strategy of control and treatment of a wide range of acute and persistent infections. Pseudomonas aeruginosa is an opportunistic bacterium with a high adaptation potential that contributes to healthcare-associated infections. In the present study, the effects of the synthesized hybrid structures bearing sterically hindered phenolic and heterocyclic moieties in a single scaffold on the production of virulence factors by P. aeruginosa were determined. It has been shown that the obtained compounds significantly reduce both pyocyanin and alginate production and stimulate the biosynthesis of siderophores in vitro, which may be attributed to their iron-chelating properties. The results of docking-based inverse high-throughput virtual screening indicate that transcription regulator LasR and Cu-transporter OPRC could be potential molecular targets for these compounds. Investigation of the impact small molecules exert on the molecular mechanisms of the production of bacterial virulence factors may pave the way for the design and development of novel antibacterial agents.

Identification of New EGFR Inhibitors by Structure-Based Virtual Screening and Biological Evaluation.

Epidermal growth factor receptor (EGFR) inhibitors have been used in clinical for the treatment of non-small-cell lung cancer for years. However, the emergence of drug resistance continues to be a major problem. To identify potential inhibitors, molecular docking-based virtual screening was conducted on ChemDiv and Enamine commercial databases using the Glide program. After multi-step VS and visual inspection, a total of 23 compounds with novel and varied structures were selected, and the predicted ADMET properties were within the satisfactory range. Further molecular dynamics simulations revealed that the reprehensive compound ZINC49691377 formed a stable complex with the allosteric pocket of EGFR and exhibited conserved hydrogen bond interactions with Lys 745 and Asp855 of EGFR over the course of simulation. All compounds were further tested in experiments. Among them, the most promising hit ZINC49691377 demonstrated excellent anti-proliferation activity against H1975 and PC-9 cells, while showing no significant anti-proliferation activity against A549 cells. Meanwhile, apoptosis analysis indicated that the compound ZINC49691377 can effectively induce apoptosis of H1975 and PC-9 cells in a dose-dependent manner, while having no significant effect on the apoptosis of A549 cells. The results indicate that ZINC49691377 exhibits good selectivity. Based on virtual screening and bioassays, ZINC4961377 can be considered as an excellent starting point for the development of new EGFR inhibitors.

Structure-based Virtual Screening, Molecular Docking, Molecular Dynamics Simulation, and Metabolic Reactivity studies of Quinazoline Derivatives for their Anti-EGFR Activity Against Tumor Angiogenesis

Epidermal growth factor receptor (EGFR/HER-1) and its role in tumor development and progression through the mechanism of tumor angiogenesis is prevalent in non-small lung cancer, head and neck cancer, cholangiocarcinoma & glioblastoma. Previous treatments targeting the oncogenic activity of EGFR's kinase domain have been hindered by acquired mutational resistance and side effects from existing drugs like erlotinib, highlighting the need for new EGFR inhibitors through structure- based drug designing The research aims to develop novel quinazoline derivatives through structure-based virtual screening, molecular docking, and molecular dynamics simulation to potentially interact with EGFR's kinase domain and impede tumor angiogenic phenomenon. Quinazoline derivatives were retrieved and filtered from the PubChem database using structure- based virtual screening and the Lipinski rule of five drug-likeness studies. Molecular docking-based virtual screening methods and molecular dynamics simulation were then carried out to identify top leads A total of 1000 quinazoline derivatives were retrieved, with 671 compounds possessing drug-- like properties after applying Lipinski filters. Further filtration using ADME and toxicity filters yielded 28 compounds with good pharmacokinetic profiles. Docking-based virtual screening identified seven compounds with better binding scores than the control drug, dacomitinib. After cross-checking binding scores, three top compounds QU524, QU571, and QU297 were selected for molecular dynamics simulation study of 100 ns interval using Desmond module of Schrodinger maestro to understand their conformational stability The research results showed that the selected quinazoline leads exhibited better binding affinity and conformational stability than the control drug, erlotinib. These compounds also had good pharmacokinetic and pharmacodynamic profiles and did not violate Lipinski’s rule of five limits. The findings suggest that these leads have the potential to target EGFR's kinase domain and inhibit the EGFR-associated phenomenon of tumor angiogenesis.

Targeting GluR3 in Depression and Alzheimer's Disease: Novel Compounds and Therapeutic Prospects.

The present study investigates the interrelated pathophysiology of depression and Alzheimer's disease (AD), with the objective of elucidating common underlying mechanisms. Our objective is to identify previously undiscovered biogenic compounds from the NuBBE database that specifically interact with GluR3. This study examines the bidirectional association between depression and AD, specifically focusing on the role of depression as a risk factor in the onset and progression of the disease. In this study, we utilize pharmacokinetics, homology modeling, and molecular docking-based virtual screening techniques to examine the GluR3 AMPA receptor subunit. The compounds, namely ZINC000002558953, ZINC000001228056, ZINC000000187911, ZINC000003954487, and ZINC000002040988, exhibited favorable pharmacokinetic profiles and drug-like characteristics, displaying high binding affinities to the GluR3 binding pocket. These findings suggest that targeting GluR3 could hold promise for the development of therapies for depression and AD. Further validation through in vitro, in vivo, and clinical studies is necessary to explore the potential of these compounds as lead candidates for potent and selective GluR3 inhibitors. The shared molecular mechanisms between depression and AD provide an opportunity for novel treatment approaches that address both conditions simultaneously.

An ensemble docking-based virtual screening and molecular dynamics simulation of phytochemical compounds from Malaysian Kelulut Honey (KH) against SARS-CoV-2 target enzyme, human angiotensin-converting enzyme 2 (ACE-2)

The human angiotensin-converting enzyme 2 (ACE-2) receptor is a metalloenzyme that plays an important role in regulating blood pressure by modulating angiotensin II. This receptor facilitates SARS-CoV-2 entry into human cells via receptor-mediated endocytosis, causing the global COVID-19 pandemic and a major health crisis. Kelulut honey (KH), one of Malaysian honey recently gained attention for its distinct flavour and taste while having many nutritional and medicinal properties. Recent study demonstrates the antiviral potential of KH against SARS-CoV-2 by inhibiting ACE-2 in vitro, but the bioactive compound pertaining to the ACE-2 inhibition is yet unknown. An ensemble docking-based virtual screening was employed to screen the phytochemical compounds from KH with high binding affinity against the 10 best representative structures of ACE-2 that mostly formed from MD simulation. From 110 phytochemicals previously identified in KH, 27 compounds passed the ADMET analysis and proceeded to docking. Among the docked compound, SDC and FMN consistently exhibited strong binding to ACE-2's active site (-9.719 and −9.473 kcal/mol) and allosteric site (-7.305 and −7.464 kcal/mol) as compared to potent ACE-2 inhibitor, MLN 4760. Detailed trajectory analysis of MD simulation showed stable binding interaction towards active and allosteric sites of ACE-2. KH's compounds show promise in inhibiting SARS-CoV-2 binding to ACE-2 receptors, indicating potential for preventive use or as a supplement to other COVID-19 treatments. Additional research is needed to confirm KH's antiviral effects and its role in SARS-CoV-2 therapy, including prophylaxis and adjuvant treatment with vaccination. Communicated by Ramaswamy H. Sarma

Identification of amentoflavone as a potent SARS-CoV-2 Mpro inhibitor: a combination of computational studies and in vitro biological evaluation

Small-molecule inhibitors of SARS-CoV-2 Mpro that block the active site pocket of the viral main protease have been considered potential therapeutics for the development of drugs against SARS-CoV-2. Here, we report the identification of amentoflavone (a biflavonoid) through docking-based virtual screening of a library comprised of 231 compounds consisting of flavonoids and isoflavonoids. The docking results were further substantiated through extensive analysis of the data obtained from all-atom 150 ns MD simulation. End-state effective free energy calculations using MM-PBSA calculations further suggested that (Ra )-amentoflavone (C3’-C8’’-atropisomer) may show a greater binding affinity towards the Mpro than (Sa )-amentoflavone. In vitro cytotoxicity assay established that amentoflavone showed a high CC50 value indicating much lower toxicity. Further, potent inhibition of the Mpro by amentoflavone was established by studying the effect on HEK293T cells treated with SARS-CoV-2 Mpro expressing plasmid. Communicated by Ramaswamy H. Sarma

Discovery of furopyridine-based compounds as novel inhibitors of Janus kinase 2: In silico and in vitro studies

Janus kinase 2 (JAK2), one of the JAK isoforms participating in a JAK/STAT signaling cascade, has been considered a potential clinical target owing to its critical role in physiological processes involved in cell growth, survival, development, and differentiation of various cell types, especially immune and hematopoietic cells. Substantial studies have proven that the inhibition of this target could disrupt the JAK/STAT pathway and provide therapeutic outcomes for cancer, immune disorders, inflammation, and COVID-19. Herein, we performed docking-based virtual screening of 63 in-house furopyridine-based compounds and verified the first-round screened compounds by in vitro enzyme- and cell-based assays. By shedding light on the integration of both in silico and in vitro methods, we could elucidate two promising compounds. PD19 showed cytotoxic effects on human erythroblast cell lines (TF-1 and HEL) with IC50 values of 57.27 and 27.28 μM, respectively, while PD12 exhibited a cytotoxic effect on TF-1 with an IC50 value of 83.47 μM by suppressing JAK2/STAT5 autophosphorylation. In addition, all screened compounds were predicted to meet drug-like criteria based on Lipinski's rule of five, and none of the extreme toxicity features were found. Molecular dynamic simulations revealed that PD12 and PD19 could form stable complexes with JAK2 in an aqueous environment, and the van der Waals interactions were the main force driving the complex formation. Besides, all compounds sufficiently interacted with surrounding amino acids in all crucial regions, including glycine, catalytic, and activation loops. Altogether, PD12 and PD19 identified here could potentially be developed as novel therapeutic inhibitors disrupting the JAK/STAT pathway.

Identification of the myxobacterial secondary metabolites Aurachin A and Soraphinol A as promising inhibitors of thymidylate kinase of the Monkeypox virus.

Thymidylate kinase (TMPK) of monkeypox virus (MPXV) has emerged as a promising target for potential therapeutics due to its significant role in pyrimidine metabolism. While smallpox drugs are advised for treating monkeypox, the European Medicine Agency has sanctioned Tecovirimat due to its potent nanomolar activity. Nonetheless, there is a need for monkeypox-specific therapeutic options. In this work, we employed docking-based virtual screening and molecular dynamics (MD) simulations to identify myxobacterial secondary metabolites as promising anti-viral natural compounds capable of inhibiting thymidylate kinase. The computational pharmacokinetics and manual curation of top-scoring compounds identified six lead compounds that were compared in terms of protein-ligand contacts and protein-essential dynamics. The study shows that among the six candidates, Aurachin A and the Soraphinol analogues such as Soraphinol A and Soraphinol C remain very stable compared to other compounds, enabling the active site integrity via a stable dynamics pattern. We also show that other compounds such as Phenoxan, Phenylnannolone C, and 8E-Aurafuron B remain unstable and have a negative impact on the active site integrity and may not be suitable binders for TMPK protein. Analyzing the Aurachin A and Soraphinol A binding, the established hydrogen bonds with Arg93 and the conserved hydrophobic interaction with Tyr101 are consistent with previous experimental interactions. Additionally, a deeper insight into the indole and the aromatic ring interaction through π-π stacking and π-cation interactions, as well as the background of Aurachin A and Soraphinol A as a bioactive compound, has significant implications not only for its potential as a promising drug but also for directing future drug discovery efforts targeting the TMPK protein.

Pareto optimization to accelerate multi-objective virtual screening

Pareto optimization is suited to multi-objective problems when the relative importance of objectives is not known a priori. We report an open source tool to accelerate docking-based virtual screening with strong empirical performance.

In-silico Investigations for the Identification of Novel Inhibitors Targeting Hepatitis C Virus RNA-dependent RNA Polymerase.

Hepatitis C is an inflammatory condition of the liver caused by the hepatitis C virus, exhibiting acute and chronic manifestations with severity ranging from mild to severe and lifelong illnesses leading to liver cirrhosis and cancer. According to the World Health Organization's global estimates, a population of about 58 million have chronic hepatitis C virus infection, with around 1.5 million new infections occurring every year. The present study aimed to identify novel molecules targeting the Hepatitis C viral RNA Dependent RNA polymerases, which play a crucial role in genome replication, mRNA synthesis, etc. Methods: Structure-based virtual screening of chemical libraries of small molecules was done using AutoDock/Vina. The top-ranking pose for every ligand was complexed with the protein and used for further protein-ligand interaction analysis using the Protein-ligand interaction Profiler. Molecules from virtual screening were further assessed using the pkCSM web server. The proteinligand interactions were further subjected to molecular dynamics simulation studies to establish dynamic stability. Molecular docking-based virtual screening of the database of small molecules, followed by screening based on pharmacokinetic and toxicity parameters, yielded eight probable RNA Dependent RNA polymerase inhibitors. The docking scores for the proposed candidates ranged from - 8.04 to -9.10 kcal/mol. The potential stability of the ligands bound to the target protein was demonstrated by molecular dynamics simulation studies. Data from exhaustive computational studies proposed eight molecules as potential anti-viral candidates, targeting Hepatitis C viral RNA Dependent RNA polymerases, which can be further evaluated for their biological potential.

GPU-optimized approaches to molecular docking-based virtual screening in drug discovery: A comparative analysis

Finding a novel drug is a very long and complex procedure. Using computer simulations, it is possible to accelerate the preliminary phases by performing a virtual screening that filters a large set of drug candidates to a manageable number. This paper presents the implementations and comparative analysis of two GPU-optimized implementations of a virtual screening algorithm targeting novel GPU architectures. This work focuses on the analysis of parallel computation patterns and their mapping onto the target architecture. The first method adopts a traditional approach that spreads the computation for a single molecule across the entire GPU. The second uses a novel batched approach that exploits the parallel architecture of the GPU to evaluate more molecules in parallel. Experimental results showed a different behavior depending on the size of the database to be screened, either reaching a performance plateau sooner or having a more extended initial transient period to achieve a higher throughput (up to 5x), which is more suitable for extreme-scale virtual screening campaigns.

Structure-guided screening of protein-protein interaction for the identification of Myc-Max heterodimer complex modulators

De-regulation of oncogenic myelocytomatosis (c-Myc or Myc) transcription factor is one of the most common molecular anomalies encountered in human cancers, and it is typically linked to many aggressive malignancies including breast, lung, cervix, colon glioblastomas, and other haematological organs. The Myc belongs to the basic helix-loop-helix zipper protein family (bHLH-ZIP), and its dimerization with another principal interactor protein partner Myc-associated factor X (Max) is essentially required for cellular transformation, cell growth and proliferation, and transcriptional activation. Intermolecular interactions have been evaluated between hetero-dimer Myc-Max protein, which identified protein-protein interaction (PPI) specific modulators using highly précised molecular docking study followed by long-range interaction stability analyzed through molecular dynamic (MD) simulation. Moreover, ADME profile analyses have been estimated for screened hit compounds. MM-GBSA-based binding free energy (ΔG) estimations have been performed for all screened hit compounds obtained from multi-step molecular docking-based virtual screening technique. According to the employed various rigorous multi-chemometric techniques, four identified inhibitors/modulators appear to have a considerable number of intermolecular contacts with hotspot residues in the hetero-dimer interface region of the Myc-Max PPI complex. However, identified hit compounds might need further structural optimization or extensive biophysical analyses for better understanding of the molecular mechanism for exhibiting the Myc-Max PPI interface binding stability. Communicated by Ramaswamy H. Sarma

Molecular docking studies of some benzoxazole and benzothiazole derivatives as VEGFR-2 target inhibitors: In silico design, MD simulation, pharmacokinetics and DFT studies

Breast cancer, a deadly disease among women, demands effective interventions due to its global impact, with over one million annual cases. Current anti-breast cancer drugs displays several side effects, also most patients resists to these drugs during early treatment stage. Hence, global search for better drugs with less side effects became necessary. The objective of this study is to identify potential inhibitors of vascular endothelial growth factor receptor-2 (VEGFR-2), a critical target in breast cancer treatment. Molecular docking-based virtual screening of 45 benzoxazole/thiazole derivatives was conducted, followed by molecular dynamics simulations to explore ligand-protein interactions. Seven ligands (compounds 7, 10, 12, 13, 14, 20, and 26) demonstrated superior binding affinities ranging from −157.85 to −173.88 ​kcal/mol (MolDock scores) and −109.96 to −129.23 ​kcal/mol (Re-rank scores) compared to Sorafenib (−156.35 ​kcal/mol MolDock score and −102.63 ​kcal/mol Re-rank score). Compound 7, identified as a potential hit, exhibited stability in a 100-ns dynamic simulation. It was chosen as a template for designing novel inhibitors, resulting in five compounds with improved binding affinities ranging from −177.84 to −184.69 ​kcal/mol (MolDock scores) and −137.34 to −143.44 ​kcal/mol (Re-rank scores). Pharmacological profiling confirmed the drug-like properties of both the potential hit molecules and the designed compounds, with 0–1 violations against Lipinski's rule of five and favorable pharmacokinetics status. Density functional theory (DFT) studies illustrated the reactive nature of the designed compounds. These findings suggest the potential of these molecules as novel VEGFR-2 inhibitors for breast cancer treatment, providing promising prospects for future drug development.

Discovery of new inhibitors of Cdc25B phosphatases by molecular docking-based virtual screening

Cell division cycle 25 (Cdc25) phosphatases play key roles in both normal and abnormal cell proliferation, which represent attractive drug targets for anticancer therapies. To discover subtype selective Cdc25B inhibitors with diverse scaffolds, a molecular docking-based virtual screening approach was performed by docking >3.8 million compounds into the Cdc25B catalytic site. An initial subset of 19 compounds was selected and assayed, and most compounds showed Cdc25B enzyme inhibition activity at 200 μM. Among these, 2 structurally diverse compounds Y11 (IC50 = 156.30 ± 7.19 μM) and Y19 (IC50 = 118.85 ± 12.07 μM) displayed micromolar inhibitory activity and subtype selectivity to Cdc25B phosphatases. Y19 was found to have anti-growth activity against the A549 cancer cell at micromolar profile. Additionally, molecular dynamics (MD) simulations and binding free energy calculations (MM/GBSA) were done to elucidate the binding mechanisms. The predicated ADMET properties suggested that Y11 and Y19 could serve as starting points for development of subtype selective Cdc25B phosphatases inhibitors.

Discovery of small molecule c-Maf inhibitors using molecular docking-based virtual screening, molecular dynamics simulation, and biological evaluation.

Multiple myeloma (MM) is a prevalent plasma cell malignancy in the blood system that remains incurable. Given the abnormally high expression of c-Maf in most MM patients, targeting c-Maf presents an attractive therapeutic approach for treating MM malignancies. In this study, we employed a combined strategy involving molecular docking-based virtual screening, molecular dynamics (MD) simulation, and molecular mechanics/generalized Born surface area (MM/GBSA) free energy calculation on existing FDA-approved drugs. Six compounds were selected for further experimental assay: vemurafenib, sorafenib, sildenafil, fluvastatin, erlotinib, and glimepiride. Among these compounds, sorafenib and glimepiride exhibited significant inhibition of myeloma cell proliferation in the RPMI-8226 cell line. Moreover, both compounds simultaneously downregulated c-Maf protein expression to induce G1 phase arrest and apoptosis in myeloma cells. Collectively, sorafenib and glimepiride may be considered promising candidates for developing more potent c-Maf inhibitors in the future.

Drug repurposing for targeting fibronectin in treatment of endometriosis and cancers

Increased concentrations of the fibronectin glycoprotein can cause ectopic tissue growth patients with endometriosis and the formation of various cancerous tumors. Furthermore, fibronectin binding to its receptors from the EDA (Extra Domain A) region contributes to promote tumorigenesis, metastasis and vasculogenesis. Thus, the EDA region can be considered a unique target for therapeutic intervention. Therefore, the present study used computational methods to identify the best fibronectin inhibitor(s) among FDA-approved drugs. First, docking-based virtual screening was performed using PyRx 0.8. Next, FDA-approved drugs that obtained favorable results in the docking phase were selected for further studies and analysis using molecular dynamics (MD) simulation. The preliminary findings of the virtual screening showed that 17 FDA-approved drugs (from 2471) had more favorable energy with their binding energy less than −9 kcal/mol. The MD simulation results of these 17 drugs showed that Avapritinib had a lower RMSD value and higher binding energy and hydrogen bonding than the other complexes in the EDA domain. Also, analyses related to the second structure changes displayed that Avapritinib in the EDA domain led to more changes in the second structure. According to the results, the anticancer drug Avapritinib forms a more stable complex with fibronectin than other FDA-approved drugs. Furthermore, this drug leads to more changes in the second EDA structure, which may have more serious potential for inhibiting EDA fibronectin. Communicated by Ramaswamy H. Sarma