key: cord-0745393-6mkd5yx9
authors: Dhanalakshmi, M.; Das, Kajari; Pandya, Medha; Shah, Sejal; Gadnayak, Ayushman; Dave, Sushma; Das, Jayashankar
title: Artificial Neural Network-Based Study Predicts GS-441524 as a Potential Inhibitor of SARS-CoV-2 Activator Protein Furin: a Polypharmacology Approach
date: 2022-05-04
journal: Appl Biochem Biotechnol
DOI: 10.1007/s12010-022-03928-2
sha: c9b46faef0c17e124f1e4a07f24049a4cb6ab1f4
doc_id: 745393
cord_uid: 6mkd5yx9

Furin, a pro-protein convertase, plays a significant role as a biological scissor in bacterial, viral, and even mammalian substrates which in turn decides the fate of many viral and bacterial infections along with the numerous ailments caused by cancer, diabetes, inflammations, and neurological disorders. In the wake of the current pandemic caused by the virus SARS-CoV-2, furin has become the center of attraction for researchers as the spike protein contains a polybasic furin cleavage site. In the present work, we have searched for novel inhibitors against this interesting human target from FDA-approved antiviral. To enhance the selection of new inhibitors, we employed Kohonen’s artificial neural network-based self-organizing maps for ligand-based virtual screening. Promising results were obtained which can help in drug repurposing and network pharmacology studies can address the errors generated due to promiscuity/polypharmacology. We found 15 existing FDA antiviral drugs having the potential to inhibit furin. Among these, six compounds have targets on important human proteins (LDLR, FCGR1A, PCK1, TLR7, DNA, and PNP). The role of these 15 drugs inhibiting furin can be established by studying further on patients infected with number of viruses including SARS-CoV-2. Here we propose two promising candidate FDA drugs GS-441524 and Grazoprevir (MK-5172) for repurposing as inhibitors of furin. The best results were observed with GS-441524. GRAPHICAL ABSTRACT: [Image: see text] SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s12010-022-03928-2.

In the recent past even after a huge development in pharma industries, new drugs or treatments approved by USFDA have remained constant [1] [2] [3] . Repurposing of the drug is always promising as not only it can shorten the time for designing but also be costeffective with randomized clinical trials [4] . Different computational platforms can be used with novel testable hypotheses for systematic drug repositioning [5, 6] . Due to the unavailability of 3D-based protein structures in some cases, traditional structure-based methods are not always followed [7] . Furthermore, drug resistance is a big challenge due to the single virus protein target as the viral genome is rapidly evolving [8] . Our group recently reported for the first time in India in 2021 about two lineages B.1.112 and B.1.99 and also revealed 1,143 unique single nucleotide variants along with 73 novel variants using the COVIDSeq approach for SARS-CoV-2 [9] . Most antiviral therapies do not completely eradicate the viral infection. Instead, they limit the immune recovery, and their long-term use results in resistance and toxicity, making space for the search of new antiviral classes with new targets. Specific antiviral agents should essentially target either the viral life cycle or the host immune clearance. There are two categories of targets required for the completion of virus replication and infection: viral molecules and host molecules. In our recent work on molecular docking-based screening, direct antiviral activity by remarkable binding of vitamins as inhibitors to various important proteins of SARS-CoV-2 is suggested. Moreover, this work revealed vitamin B12 to be an efficient nutraceutical for almost all the nCoV-2 drug targets including the host protease furin which is a human host molecule and one of the promising targets for therapeutic intervention [10] . It cleaves more than 150 mammalian, viral, and bacterial proteins including the disastrous virus coronaviruses (CoVs). Researchers in their recent works have found furin inhibitors as a new avenue for the treatment of COVID-19 [11] . The mechanism for furin-mediated protein processing for health and disease is well-known. Furin is a member of the pro-protein convertase family, PCSK 3 representing furin. The viral and bacterial infection is caused by the production of various toxins by proteolytic cleavage such as Anthrax toxin, Diphtheria toxin, Pseudomonas exotoxin A, and Shiga toxin. Furin activation of viral proteins is essential for their entry to the host cell. Furin has basic amino acid motifs as a cleavage site and known as "PACE" paired basic amino acid cleaving enzymes. It is located at chromosome 15 q26.1 and expressed by the FUR (FES upstream region) [12] . The hot spot cleavage site of furin is known as R-X-K/R-R↓ [13] . Bioinformatic analyses have revealed more than 100 furin cleavage sites in humans such as growth factors, cytokines, collagens, adhesion molecules, MMP (metalloproteinases), and receptors [14] . Many recent reports have mentioned peptidic and nonpeptidic inhibitors for therapeutic inhibition of furin resulting in its reduced proteolytic activity. Numerous pharmaceutical companies have come up with solutions for the current global pandemic caused by SARS-CoV-2 with different vaccines and existing medicines. However, vaccination is going on in spite of which many cases are being observed and no effective medications are there due to the appearance of double mutant and triple mutant strains. Hence, the need of the hour is a better prognosis of the Covid patients.

Cheminformatics incorporating machine learning are breakthrough techniques imperative in the pharmaceutical field today as we live in an era of infectious outbreaks, rising lifestyle diseases, and the age-old time-consuming procedures in drug discovery [15] . Artificial intelligence (AI) and machine learning (ML) tools are applied for furin targeting while screening of molecules including higher-dimensional similarity/differences and networking observation has been performed to establish a relationship (interactome) with selected drug compounds and their respective important host factors along with furin.

Kohonen's self-organizing maps (SOM) are extensively used today to establish the relationship of structure to biological activity and drug discovery in many pharmaceutical research. This methodology is widely applied in drug repurposing and scaffold hopping. Self-organizing maps (SOM) are artificial neural network (ANN) that implements virtual screening for assembling thematic molecular libraries. The concept of the SOM technique, case studies, and potential future applications are presented according to the report of Schneider et al. [16] . Here the ANN is employed to identify the nonlinear connections of arbitrary dimensions into processable real numbers of one or two-dimensional maps at less computational costs. This mapping technique permits high-dimensional input space to be transformed into lower-dimensional output space [17] . Studies with SOM supervised methods were proven more accurate than standard linear QSAR methods [18] . ANN is not a single algorithm; instead, it comprises a framework of machine learning algorithms to process complex data. Kohonen's self-organizing map is employed here to screen and predict the active FDA antivirals on furin, basically to transform a pattern of arbitrary dimension into a 2D map adaptively.

Managing large amounts of heterogeneous data in a short period is becoming possible through AI technology. AI/ML, systems biology, including computer-aided drug repurposing are speedily becoming an integral part of the global workflow, to prioritize novel drug targets, new chemical entities, and to evaluate off-label or drug repositioning proposals [19] . All these above-mentioned techniques are utilized in our present work where we focused to present an integrative antiviral drug repurposing methodology for the first time with the combination of systems pharmacology-based network medicine platform that quantifies the interplay between the virus-host interaction and drug targets in the human PPI network. An ANN employed polypharmacology-based computational approach has been reported in the present paper with multiple antiviral drugs for furin, which uncovers the therapeutic interventions.

Human furin protein had 29 bioassays containing 221066 tested compounds in PubChem Repository on 07-07-2020 (https:// pubch em. ncbi. nlm. nih. gov/ prote in/ EAX02 111). All furin actives and inactives were taken for analysis except inactive of AID 463115 (~2.2 lakh) to minimize the bias due to imbalance of diversity of actives and inactives [20] . The 386 actives and 93 inactives in structure download format (SDF) were downloaded from the PubChem repository which was tested on human furin. These compounds were used for descriptor generation and subsequent model building.

Among the 21 PDB structures available, 5 MIM along with the bound ligand C 28 H 37 N 13 O 2 (CID -44593169, IUPAC Name:2-(4-((1R,2R,4S,5S)-2-(bis (azanyl) methylideneamino)-5-(4-(bis(azanyl)methylideneamino)phenoxy)-4-((4-(bis(azanyl)methylideneamino)phenyl) amino)cyclohexyl)oxyphenyl) guanidine), which is considered reference molecule in this study, was downloaded from the protein data bank as it has a single chain and zero mutations with a good X-ray Structure Validation Report (https:// www. rcsb. org/ struc ture/ 5MIM). Antiviral drugs which come under the approved/investigational category in drugBank (https:// go. drugb ank. com/) were selected as screening dataset and their SDFs (total 180) were downloaded from the PubChem repository.

The SDF files of furin actives and inactives downloaded from PubChem bioassays were used to build a model using self-organizing map, the unsupervised machine learning clustering algorithm in the Canvas module of the Schrodinger suite (Schrodinger, Inc., LLC, New York, USA), and to screen the activity of the FDA approved/investigational antiviral drugs. The physicochemical, topological, LigFilter, and the QikProp descriptors generated were fed as input into the neural network which in turn had undergone artificial neural network-based competitive learning to develop the model cluster. The SOM is generated by the property-based clustering using all available descriptors of categories mentioned above generated by the software.

The protein was pre-processed in the protein preparation wizard of Maestro 9.3 (Schrodinger, Inc., LLC, New York, USA) and the states were generated at a pH 7.0 ± 2.0. The metals and the water molecules (Beyond 5A o from the hit groups) were removed and minimized in the OPLS-3e force field. The ligand is also prepared using the LigPrep wizard of maestro, taking all stable conformations up to 32. The grid for docking is prepared by taking the centroid of the PDB inbuilt inhibitor in Glide. Ligand docking is done with the prepared protein on the generated grid with the prepared ligand (the ligands from the best clusters). The glide score along with the glide e-model is tabulated and analyzed for binding affinity.

The interactome was built for (i) FDA drug and human target network and (ii) FDA drug and viral target network in connection with human target furin. All the target proteins examined in this study are of antiviral FDA drugs (supplementary table S7 and S8) which are considered active on furin and obtained from the drug bank. The information of targets for each drug is collected from the UniProt ID and the interactions downloaded from biogrid interactions. The network is built and visualized in Cytoscape 3.8.2 (https:// cytos cape. org/).

The Desmond module of Schrodinger was used to undertake a molecular dynamic (MD) simulation of the Furin-GS441524 complex. MD simulations of these compounds were run for 100 ns. Initially, the protein-ligand complex was prepared using the Schrodinger suite's Protein Preparation Wizard with default settings that included hydrogen and bond reassignment, addition of missing side chain atoms of amino acid residues, optimization of loops residues, and water orientation sampling at pH 7.4. Using the System Builder module, a periodic simulation box was created, the system was solvated using the transferable intermolecular potential with 3 points (TIP3P) water model and neutralized by adding counter ions, and energy minimization was achieved using the OPLS with 1000 iterations of the steepest descent technique (optimized potentials for liquid simulations) allatom force field. After reaching equilibrium, an unrestrained production phase with NPT (constant atoms, pressure, and temperature) ensemble was progressed for 100 ns under the supervision of a Nose-Hoover thermostat (300 K, relaxation time = 1 ps) and an isotropic Martyna-Tobias-Klein barostat (1.01325 bar, relaxation time = 2 ps). The smooth particle mesh Ewald (PME) approach with the RESPA integrator was used to evaluate shortrange interactions (cut-off = 9 Å) and long-range Coulomb interactions. At a 5 ps period, the frames capturing the system's dynamic motions were exported. Plotting histograms for RMSD, root mean square fluctuations (RMSF), hydrogen bond analysis, radius of gyration (Rg), and torsional bonds was used to investigate the system's stability [21] .

The furin actives and inactives are manually curated from the PubChem database (supplementary table S1and S2) for model building for the ligand-based virtual screening (LB-VS). ANN-based self-organizing maps were generated using unsupervised machine learning clustering algorithms in the canvas Schrodinger suite, and using the above-mentioned model, the activity of the FDA-approved/investigational antiviral drugs was screened (supplementary table S3). The descriptors generated in the canvas were applied for the ANN modeling inputs and the SOM was generated by the property-based clustering. However, SOM is considered a supervised training method in this study and we employed the actives as input data without labeling while submitting the job; Most of the time, the output-data analysis became more straightforward with known labels [22] . The number of molecules from each inadequate dataset is presented in Table 1 .

The SOM results are presented in distinct rows and columns in specific color codes that help us to scrutinize the high dimensional data of complex datasets into lower-dimensional . 1a ). Fig. 1b illustrates the Euclidean distances in the hyperspace between the molecules. The difference between Euclidean distances of data points was large and was depicted as red, whereas small distances were in blue. The variations among the population and the distance were interpreted by various colors between blue and red. The bioassay actives are presented in blue. The slightly variant molecules are depicted in green for a very small distance in this analysis and they resemble property space as well. The greener section owned higher numbers of FDA drugs. Each cluster was manually examined to confirm the most desirable molecules comprising only actives along with FDA drugs and restrict inactives molecules. The clusters in the combination of all three active/inactive/FDA were discarded. Out of 15 clusters, the clusters holding more than 10 molecules were selected for further analysis. Cluster number 2 holding 33 molecules (Supplementary Table S4 ) has the maximum number of FDA antivirals and actives together (20 actives + 13 FDA). Cluster number 4 carrying 27 molecules (25 actives + 2 FDA) is the second-largest membered cluster (Supplementary Table S5 ) yet possessing just two FDA drugs (GS-441524 and Grazoprevir). Cluster 4 possessing active ligand C 28 H 37 N 13 O 2 (PubChem CID: 44593169) which was in interaction with furin from the protein structure database (PDB Id 5MIM) is acknowledged as the reference active molecule in this study [23] . Cluster numbers 5, 6, and 8 also own the active and FDA drugs screened out and granted for additional analysis.

The physicochemical descriptors are generated in the canvas, and the important ones are tabulated in Table 2 . The higher values of molecular weight, AlogP, rotatable bonds, 

Molecular docking analysis is a vital tool for structure-based virtual screening (SB-VS) to predict the promising structural conformations among ligand and active sites of the receptor. The SOM examination revealed five clusters 2, 4, 5, 6, and 8 to be significant for additional knowledge; thus, all active ligands as well as FDA drugs from these clusters were docked against furin. The docking result revealed the conformation, and the compound with the lowest docking score was chosen as the most efficient antiviral. The top 15 FDA antiviral drugs with an immeasurable docking score (Fig. 3a) with furin are from clusters 2 and 4. The findings indicated that 60 compounds from SOM clusters 2 (33 molecules) and 4 (27 molecules) had docking scores that are quite similar to each other (Fig. 2a) and with the majority of their energy falling in the −40 to −90 range (Fig. 2b) . Only a few drugs with high docking scores are not comparable to reference ligands, such as Porfimer sodium, best docked FDA drug which has too many contacts even beyond the specified grid box and in the model. The top-scored FDA drug was considered an outlier in the context of pharmacokinetic and topological descriptors ( Table 2 ). The docked position and orientation of C 28 H 37 N 13 O 2 (considered active/ (Fig. 3a, b) as well as key receptor contacts ( Supplementary Fig S2-S5 ) better than co-crystal ligand (C 28 H 37 N 13 O 2 ). Considerable amino acid residues of high interaction frequency were common among the top molecules obtained with a variety of interaction types like hydrogen bonds, π-π stacks, and hydrophobic contacts depicted in supplementary Table S6 . All FDA drugs were screened from docking analysis and their former targets were further considered for the network analysis. 

The pharmacokinetic descriptors are high for most of the drugs listed in Table 2 and can have promiscuity [24] . It is necessary to analyze all possible polypharmacology before we practice FDA-approved drug candidates on any other targets. For this purpose, we utilized the network-based interaction study of the drugs selected to repurpose. This analysis aims to check the interaction network of present targets before we move further, whether it is in host or pathogen and how costly it becomes while repurposing it on furin. Concerning FDA drugs with multiple targets to be active against furin, its target proteins should be within or near the corresponding subnetwork in the host-pathogen interactome. The structure-based virtual screening outcomes uncovered 15 promising antiviral FDA drugs for repurposing against the furin even though one was an outlier from the reference ligand values. The molecular targets of FDA antivirals are arranged in Table 3 .

The drug-target network exploration represents that nodes (red with label) at the center are human targets and white with label is their interactions. The plotted networks manifest six molecules in the interaction network shown in red as Porfimer sodium (LDLR, FCGR1A), Isatoribine (TLR7), GppCp (PCK1), Didanosine (PNP), Valganciclovir, and Entecavir (Fig. 4) . A drug with multiple targets can be efficient for a disease or may cause off-target side effects. Hereabouts, the targets are associated with multiple significant biological pathways; therefore, FDA drugs associated with these targets are eliminated in the further analysis. The proteins that act as a drug target for one disease may also be appropriate drug targets for possible antiviral infections due to common protein-protein interactions (PPI) and functional pathways explained by the human interactome [25] . From network analysis, we found that human protein furin is directly associated with viral drug targets in the PPIs. The most abundant group related to furin is SARS-CoV-2-related interactome. We identified that furin is exclusively found in the GS-441524 drug's subnetwork. Conversely, Valaciclovir, Penciclovir, Acyclovir, and Ganciclovir presently acting on human herpesvirus (HHV) protein UL30 (DNA polymerase), and no interaction is observed with human proteins. Grazoprevir is associated with the NS3/4A protein of the hepatitis C virus but no interactions to analyze. The network-based drug-disease vicinity sheds light on the connection between drug targets and molecular factors in disease pathobiology modules within the PPIs and can serve as a useful tool for efficient screening of potentially new suggestions for FDA-approved drugs. In brief, the up-to-date interactome examination disclosed that GS-441524 can be a repurposed drug on furin (Fig. 5 ).

The molecular dynamics simulation approach has been systematically employed to inspect features of conformational changes leading to modifications in protein-ligand interactions, molecular dynamics, and protein folding [26] . In this study, according to SOM clustering and network-based polypharmacology analysis, it was found that GS-441524 can be repurposed on human furin without many complications as per the available pieces of information today. The MD simulations of the furin-GS-441524 complex were accomplished with the Schrodinger Desmond module. The resulting MD trajectory was utilized to examine the dynamics equilibrium as a function of time. The root mean square deviations (RMSD) metric of the initial structure and the average simulated structure of all MD trajectory frames can be used (Fig. 6a) the stability of the docked complex which attained only after 20 ns which is in good agreement with the simulation of ligand. A thorough analysis of RMSD of these fluctuating residues was found to be enriched with loop structure that showed peaks about ~1.8 Å in respective plots. The GS-441524 complex scale of fluctuation and the divergence between the average RMSD values acclaimed that simulation produced stable trajectories. The secondary structural elements corresponding to pocket residues were intact in the β-sheets and loop regions of furin (Supplementary Fig S6) . Therefore, these outcomes established that FDA drug GS-441524 is more enduring complex with furin. The backbone RMS fluctuations were also calculated as they can postulate the clue about backbone mobility. RMSF protein fluctuations were highest in the residue index window at 75-80, 120-130, 145-160, and 345-370 positions as some of the amino acid residues of this window were pocket residues promoting ligand binding (Fig. 7) . In Figure 6d , the radius of gyration (Rg) graph illustrates that the circulation of atoms of a protein around its axis in all the complexes is stable. The Rg of furin-GS-441524 complex is 3.1 to 3.3 Å which qualifies one to assess the compactness changes of a complex. Intramolecular hydrogen bonding was detected (Fig. 6c) which accounted for stronger contacts for prolonged binding residence within the ligand-binding site. The molecule surface area (MSA) of ligand was around 450 Å 2 in most of the frames in the trajectory (Fig. 6d) . The solvent-accessible surface area (SASA) and polar surface area (PSA) for ligand are more than 250-270 Å 2 in most of the intervals (Fig. 6 e and f) . Overall analysis, spectacles furin in complex with GS-441524 was Fig. 8a . The 2D interaction plot of re-docked ligand exhibiting contact retention along the simulation course is shown in Fig. 8b .The interaction fraction is shown on the Y-axis, indicating the percentage of that particular contact, while numerous connections of the residue with ligand are shown in two or three different colors in a bar column. GS-441524 in its docking profile developed four hydrogen contacts (ASH 153, SER 253, PRO 256, and ASP 258) and three hydrophobic contacts (LEU 227, TRP 254, PRO 256) and three polar contacts (SER 253, SER 368, ASN 295) which were observed in the intermolecular interaction profiles. After the MD simulation, the protein-ligand complex is proximal to the in vitro environment; hence, further wet-lab experiments may qualify the GS-441524 as a furin inhibitor. 

In the present article, we focused on our particular interest in drug repurposing, since already approved drugs can go through shortened clinical trials. As our knowledge is exponentially expanding with constantly emerging "omics" (e.g., transcriptomic, proteomic, metabolomic), data on drug/disease-associated mechanisms could easily be incorporated and analyzed. Our group working in this area has reported through immunoinformatics approach for immunogenic protein from the proteome of Leptospira and a codon-based analytical study for susceptible hosts of Nipah Virus [27, 28] . Our work will add to a broad arena towards the analysis of quantitative structure-activity relationships, pharmacophores, homology models, and other molecular modeling approaches like deciding conformational stability [29] as well as network analysis of biological functions, machine learning, and any other analysis tools that include using a computer. In this study, each molecule is considered a data point, no labeling was done during submission of data for analysis, but generated descriptors propose the input into the neuron vectors which in turn make the input nodes. The data processing can be categorized into three broad interconnected networks in ANN, the input layer, the hidden layer, and the output layer. The connections in the network were assigned with weighted values which help the learning process and create the nodes in the hidden layers and finally the optimized signal after all the iterations are passed on to the output layer. According to the input given through descriptors of each data point, the molecules, within the different neurons in the network, there will be competitive and cooperative learning iterations by the end of which the winning neuron among the neuron vectors was determined and will be the closest to the molecule or the data point. After the determination Fig. 8 a Various intermolecular interactions with color scheme (H-bonds: green; hydrophobic: purple; ionic: pink; water bridges: blue) made by furin pocket residues with ligand GS-441524. b Preserved contacts of ligand with furin target, captured during molecular dynamics simulations of the winning neuron, update it with the weight of the neighboring neurons to place each one in possible closest positions, this process is continued until the best optimization is achieved and the clustered molecules have fallen into the cells of the map [30] (Fig. 1a) . Each cluster was manually analyzed to confirm the best ones containing actives together with FDA drugs and have no inactives. Out of 180 FDA antivirals, 15 predicted by clusters 2 and 4 demonstrated good results in molecular docking studies. However, it is necessary to analyze polypharmacology before we take a possible drug onto any other targets and not that binding affinity and matching properties alone. In recent years, researchers are becoming dependent on a computational approach for their study on the prediction of targets, off-targets, and polypharmacology of hit compounds. Thus, the network pharmacology comprises components functioning in stepwise order "multi drugs to multi-targets and multi-diseases [31] ."

There are many important host proteins like FCGRIA, LDLR, PNP, PCK1, and TLR 7 that have well-established interactions with the selected drugs. All of them are having more or less connections with various diseases (viral, tumor, etc.) positively or negatively. Purine-nucleoside phosphorylase (PNP) is a very important enzyme that catalyzes nucleotide synthesis in the salvage pathway. Inherited deficiency of PNP causes T cell deficiency and abnormal B cell function in some patients as adenosine deaminase activity is affected. The enzyme deficiency was observed in HIV patients [32] . FCGRIA is a leukocyte membrane receptor protein capable of binding the Fc domain of IgG with high affinity and signaling for a variety of inflammatory and proinflammatory responses. In human immunodeficiency virus (HIV) infection, it has a role in sepsis and rheumatoid arthritis. FCGRIA has been proposed as an attractive target for immunotherapy by various workers [33] . Phosphoenolpyruvate carboxykinase (PEPCK also known as PCK, EC number 4.1.1.32), a cytoplasmic protein, catalyzes gluconeogenesis implicating its important role in cancer biology and virology. Anti-tumor treatment based on gluconeogenesis has become an advanced concept nowadays. According to a recent report, gluconeogenic pathway inhibition of PCK1 enhanced HBV replication [34] . Several repeated amino acid sequences are present in the ligand-binding domain for LDL in LDLR (low-density lipoprotein receptor). Lipoprotein receptor-mediated HCV entry into the host cell is also well documented.

We target furin as part of the treatment of infectious diseases by making sure that the other targets plotted here will not get disturbed. But if we take the drugs presently acting on the 5 targets of human (FCGRIA, LDLR, PNP, PCK1, and TLR 7) and repurpose them on furin, their existing network of interactions will be disrupted; All the four drugs which are known to interact with these targets, namely Porfimer sodium, Isatoribine, 5′-Guanylyl methylene bisphosphonate, and Didanosine, are not to be carried further to repurpose without detailed studies. We strongly recommend collecting detailed data regarding the effect on furin from the patients who take these medicines and the repercussions due to this while infected by viral diseases like SARS-CoV-2. The drugs Entecavir and Valganciclovir also need proper studies before attempting to repurpose as their known current target is human DNA. Although no interactions are available, no inference can be made from this study at present. We collected information from drug bank followed by taking details about all the drugs and in this discussion the two best-docked molecules Porfimer Sodium and Isatoribine are used, former is for esophageal cancer treatment and the latter is a compound that elevates levels of interferon-alpha.

After the poly-pharmacology analysis followed by molecular dynamics simulation, we found GS-441524 to be the most promising molecules to repurpose against the furin. GS-441524 is the nucleoside, a metabolite of the phosphoramidate prodrug Remdesivir. GS-441524 was initially developed by Gilead Sciences, Inc. for the treatment of filovirus infections. It is a 1′-cyano substituted C-nucleoside derived from 4-aza-7, 9-dideazaadenosine. Remdesivir (GS-5734) is a prodrug of GS-441524 used for the treatment of Ebola infections [35] . Current investigations demonstrate that this molecule has been designed to target the RNA-dependent RNA polymerase (RdRp) of several viruses, including the Ebola virus, MERS-CoV, SARS-CoV, and potentially SARS-CoV-2. The latest study has further confirmed that the direct use of GS-441524 can potently inhibit the replication of SARS-CoV-2 in a mouse model [36] . In the present study, another prominent antiviral molecule Grazoprevir was screened out as a furin inhibitor. Grazoprevir is a hepatitis C protease inhibitor given in combination with Elbasvir (under the brand name Zepatier) for the treatment of chronic hepatitis C virus (HCV) infection in adults. Up-to-date examinations unveiled that it has striking anti-COVID activity for the essential proteins involved in viral entry into host cells ACE2 and TMPRSS2 [37] . The HCV drugs Simeprevir, Vaniprevir, Paritaprevir, and Grazoprevir in combination with Remdesivir gained the effectiveness of Remdesivir alone by 10-fold [38] . Nguyenla et al. (2020) published that Remdesivir has a synergistic effect with Velpatasvir and Elbasvir the HCV drugs [39] . Additionally in combination with HCV partner drugs, Sofosbuvir and Grazoprevir promote efficacy, synergistically with Remdesivir in a human lung epithelial cell line infected with a clinical isolate of SARS-CoV-2. Herewith we strongly recommend that Remdesivir analog GS-441524 alone or in combination with Grazoprevir may have a possible inhibitory effect on furin. Nonetheless, further wet-lab experiments may shed light on this investigation.

A number of very recent researches are evidencing the importance of GS-441524 as a potent antiviral drug due to its strong binding properties with NSP3 macrodomain and NSP12-NSP8-NSP7 of SARS-CoV-2 [40, 41] . It has also been proven that this molecule has very less retention value inside the host body as most of the residue is eliminated with urine suspecting no side effect [42] . Amazingly, out of the majority of Remdesivir in clinical trial experiments, recovered in urine was found to be in the form of GS-441524, whereas only 10% recovered as Remdesivir [43, 44] . Moreover, no report has shown any host protein interaction except some metabolizing enzymes like CPY450. Our work has achieved the interesting interaction of GS-441524 with furin, though indirectly. Hence, results of our work could be a very powerful support to the scientific community to take a final decision towards repurposing GS-441524 as SARS-CoV-2 antiviral drug.

The present study evidenced that all the 15 FDA-approved antiviral drugs have the potential to inhibit the pro-protein convertase furin. These drugs are presently being successfully used to prevent many diseases (viral, tumor, etc.). However, furin inhibiting properties of these drugs in these patients can be investigated to support our results. With proper studies digging deep into the biological/interactome/pharmacokinetics studies, all these 15 FDA-approved drugs can be repurposed according to the pros and cons of each disease, but we can propose with existing data that GS-441524 be repurposed on human furin without developing many complications as per the available information today. These results can open up more options to reduce death tolls in comorbidity conditions during SARS-CoV-2 infections.

The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s12010-022-03928-2.

Putting the patient back together -Social medicine, network medicine, and the limits of reductionism

Macrolactin A as a novel inhibitory agent for SARS-CoV-2 Mpro

Molecular landscape and computational screening of the natural inhibitors against HPV16 E6 oncoprotein

A comprehensive map of molecular drug targets

Individualized network-based drug repositioning infrastructure for precision oncology in the panomics era

In silico oncology drug repositioning and polypharmacology

Nature potential for COVID-19: Targeting SARS-CoV-2 Mpro Inhibitor with bioactive compound

Coronaviruses-drug discovery and therapeutic options

High throughput detection and genetic epidemiology of SARS-CoV-2 using COVIDSeq next-generation sequencing

Unravelling Vitamin B12 as a potential inhibitor against SARS-CoV-2: A computational approach

Furin: A potential therapeutic target for COVID-19. iScience

Tissue-based map of the human proteome. Science (80-)

20 residues motif delineates the furin cleavage site and its physical properties may influence viral fusion

FurinDB: A database of 20-residue furin cleavage site motifs, substrates and their associated drugs

Applications of self-organizing neural networks in virtual screening and diversity selection

Self-organizing maps in drug discovery: Compound library design, scaffold-hopping, repurposing

Virtual screening applications: A study of ligandbased methods and different structure representations in four different scenarios

Supervised selforganizing maps in drug discovery. 1. Robust behavior with overdetermined data sets

Artificial intelligence, drug repurposing and peer review

Virtual screening and repositioning of inconclusive molecules of beta-lactamase bioassays-A data mining approach

Making optimal use of empirical energy functions: Force-field parameterization in crystal space

Feature optimization in high dimensional chemical space: Statistical and data mining solutions

Structural studies revealed active site distortions of human furin by a small molecule inhibitor

Improving drug candidates by design: A focus on physicochemical properties as a means of improving compound disposition and safety

Network-based approach to prediction and population-based validation of in silico drug repurposing

Nanoinformatics: An emerging trend in cancer therapeutics

Analysis of Nipah virus codon usage and adaptation to hosts

Exploring Leptospiral proteomes to identify potential candidates for vaccine design against Leptospirosis using an immunoinformatics approach

A decision making components in cyclisation of mannose derivatives-A computational approach

Artificial neural networks for computer-based molecular design

Polypharmacology: Challenges and opportunities in drug discovery: miniperspective

Fcγ receptor 1 (CD64), a target beyond cancer

Elevated adenosine deaminase and purine nucleoside phosphorylase activity in peripheral blood null lymphocytes from patients with acquired immune deficiency syndrome

Increased expression of phosphoenolpyruvate carboxykinase cytoplasmic isoform by hepatitis B virus X protein affects hepatitis B virus replication

Therapeutic efficacy of the small molecule GS-5734 against Ebola virus in rhesus monkeys

Remdesivir metabolite GS-441524 effectively inhibits SARS-CoV-2 infection in mouse models

Computational drug repurposing study elucidating simultaneous inhibition of entry and replication of novel corona virus by Grazoprevir

Hepatitis C virus drugs that inhibit SARS-CoV-2 papain-like protease synergize with remdesivir to suppress viral replication in cell culture

Stroumza J, (2020) et al. Discovery of SARS-CoV-2 antiviral synergy between remdesivir and approved drugs in human lung cells. bioRxiv

Structural insights into plasticity and discovery of remdesivir metabolite GS-441524 binding in SARS-CoV-2 macrodomain

Comparison of binding site of remdesivir and its metabolites with NSP12-NSP7-NSP8, and NSP3 of SARS CoV-2 virus and alternative potential drugs for COVID-19 treatment

Pharmacokinetic, pharmacodynamic, and drug-interaction profile of remdesivir, a SARS-CoV-2 replication inhibitor

What Do we know about remdesivir drug interactions?

Reviews on biological activity, clinical trial and synthesis progress of small molecules for the treatment of COVID-19

The authors are grateful for the support given by the ministry of education and innova-

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

M. Dhanalakshmi 1 · Kajari Das 2 · Medha Pandya 3 · Sejal Shah 4 · Ayushman Gadnayak 5 · Sushma Dave 6 · Jayashankar Das 7 1 Research and Development Centre, Bharathiar University, Marudhamalai Rd, Coimbatore, Tamil Nadu, India