key: cord-0125912-hrpzz7li authors: Dayer, Mohammad Reza title: New Candidates for Furin Inhibition as Probable Treat for COVID-19: Docking Output date: 2021-07-27 journal: nan DOI: nan sha: a503c88fe317f401de945fe53a0f0e8ec2340d88 doc_id: 125912 cord_uid: hrpzz7li Furin is a serine protease that takes part in the processing and activation of the host cell pre-proteins. The enzyme also plays an important role in the activation of several viruses like the newly emerging SARS-CoV-2 virus that causes COVID-19 disease with a high rate of virulence and mortality. Unlike viral enzymes, furin owns a constant sequence and active site characteristics and seems to be a better target for drug design for COVID-19 treatment. Considering furin active site as receptor and some approved drugs from different classes including antiviral, antibiotics, and anti protozoa/anti parasites with suspected beneficial effects on COVID-19, as ligands we have carried out docking experiments in HEX software to pickup those capable to bind furin active site with high affinity and suggest them as probable candidates for clinical trials assessments. Our docking experiments show that saquinavir, nelfinavir, and atazanavir with cumulative inhibitory effects of 2.52, 2.16, and 2.13 respectively seem to be the best candidates for furin inhibition even in severe cases of COVID-19 as adjuvant therapy, while clarithromycin, niclosamide, and erythromycin with cumulative inhibitory indices of 1.97, 1.90, and 1.84 respectively with lower side effects than antiviral drugs could be suggested as prophylaxes for the first stage of COVID-19 as a promising treat. represents the primary structure of the inactive form of furin on synthesis. As it is depicted furin comprises a short signal peptide at the N-terminal followed by a prodomain ending at a cleavage site of R107 upon activation. The prodomain helps correct the folding of the next catalytic domain initially. This domain is also known as inhibitory propeptide since it inhibits the proteolytic activity of furin until its removal. This domain is cleaved off by autoproteolytic of activated forms of furin or by other subtilisin-like proteases. The next domain is called the catalytic domain which contains an essential triad of Asp, His, and Ser that take part in proteolytic activity. P-domain is the following domain playing role in regulating enzyme activity at different pH and calcium ions concentrations [4] [5] . Cysteine-rich domain (CRD) is the next functional domain ending with the next cleavage site SR683 residues. At certain conditions, this site is also cleaved off by active furin or by other subtilisin-like enzymes. Under this condition furin leeks to the extracellular compartment with retained activity. The ultimate domain at the C-terminal site of furin is the transmembrane (TM) domain that tight the enzyme to membranes of the endoplasmic reticulum [6] [7] . Besides the physiological role, furin is also recruiting to process some pathogens proteins such as envelope proteins of HIV, influenza, and several filoviruses of Ebola and Marburg virus and also spike protein of SARS-CoV-2 and therefore fully activates the pathogens [8] [9] . The tertiary structure of furin, Figure 2 contains active site cavity lined by negative charge residues, where the enzyme substrates or inhibitors can bind the enzyme and come to contact with the catalytic triad of Asp153, His194, and Ser368 which is important for catalysis. The presence of negatively charged residues in the catalytic site explains the requirement for substrates/inhibitors to carry positive charges for effective binding and also explains the cleavage point of proteins to be at positive residues of Lys and/or Arg [10] [11] . As it is indicated in figure 2 there are two or three calcium ions depend on the origin of furin that is not involved in the catalytic process but are essential for enzyme native conformation [12] [13] . It is well documented that furin is up-regulated in several conditions, like diabetes mellitus, cancer, and viral infections that is suspected to play role in diseases deterioration and so comprises a potential target for drug development and inhibition [14] [15] . COVID-19 among the newly emerging viral disease that is caused by SARS-CoV-2 comprises a serious pandemic worldwide with more than 194 million cases and 4.16 million deaths by Jul 2021. The diseases encourage investigators to search for a vaccine for healthy individuals or drugs for patients. It is shown that extracellular furin plays a critical role in SARS-CoV-2 infectivity [16] [17] . The virus uses its surface spike or S protein to binds the host cell receptor called angiotensin-converting enzyme 2 (ACE2). This protein contains N-terminal or S1 domain, which is responsible for receptor binding, and a C-terminal or S2 domain, for host cell fusion and entrance. Furin by cleaving the spike at S1/S2 cleavage site accelerates virus entry and pathogenesis [18] [19] [20] . One strategy against viral infections is the application of inhibitors against viral enzymes to inhibit viral amplification. The main obstacle in this way is the drug resistance posed by a high rate of replication errors. Inhibition of furin with no drug resistance seems to be a good approach to prevent viral infections especially in the case of COVID-19 [21] [22] . According to this introduction and during this work we tried to search for drugs capable to bind furin active site and inhibition among known anti-viral and approved antibiotics based on their structural similarities through molecular docking calculations with the hope to find and suggest new candidates to fight the SARS-CoV-2 virus. Furin coordinate structure: The coordinate structure of furin with PDB ID of 4OMD was retrieved from the protein data bank (https://www.rcsb.org/). The structure was obtained by the X-ray diffraction methods and refined at the resolution of 2.70 Å. Since the protein structure was refined from dried crystal, so its conformation will be far from its native structure in physiological conditions. Accordingly, we optimized and equilibrate the protein structure via minimization to under 200kJ/mol in pH 7.5, 37 degrees centigrade, and 1atomsphere of pressure using in GROMACS 4.5.5 software (http://www.gromacs.org) and GROMOS force field and steepest descent algorithm. The structure was placed in a rectangular box with a dimension of 6.11×7.48×7.27nm filled with SPCE water [23] [24] . The chemical structures of candidate inhibitors, including antiviral drugs of amprenavir, atazanavir, baloxavir, darunavir, disoproxil, emtricitabine, indinavir, lamivudine, lopinavir, nelfinavir, nevirapine, oseltamivir, remdesivir, ritonavir, saquinavir, tenofovir, tipranavir, zidovudine; antibiotics of azithromycin, cefaclor, cefazolin, cefdinir, cefditoren, cefixime, cefotaxime, cefpodoxime, cefprozil, ceftizoxime, ceftriaxone, cefuroxime, ciprofloxacin, clarithromycin, doxycycline, erythromycin, fidaxomicin, gemifloxacin, imipenem, moxifloxacin, ofloxacin, sulfamethoxazole, tetracycline; anti protozoa/anti parasite of diminazene, i, and niclosamide, in SDF format were obtained from PubChem database (https://pubchem.ncbi.nlm.nih.gov/) and converted to PDB format in Open Babel software (http://openbabel.org/) and energy minimized in ArgusLab software (http://www.arguslab.com/) [25] . Furin Active Sites: Furin's active site was extracted from the PDB structure of 4OMD i.e., the binding site of competitive inhibitor using ArgusLab software (http://www.arguslab.com/) [25] . To carry molecular docking experiments the optimized structure of furin with PDBID:4OMD was used as a receptor and the coordinate structures of inhibitor in PDB formats were used as a ligand for blind docking experiments in Hex 8.0.0 (http://www.loria.fr/~ritchied/hex/) software [26] . The sahpe+electrostatic and macro sampling modes of docking were used for docking and the best 100 poses and their binding energies were saved for further analysis. The logP or partition coefficient for drugs is a known index for hydrophobicity, the more and positive values reflecting the more hydrophobic nature for a chemical and vice versa [27] . The Virtual Computational Chemistry Laboratory (http://www.vcclab.org/) server was used to calculated logP for each drug [28] . All the numerical data were exported to Excel and SPSS software for analysis. A P-value lower than 0.05 was considered as the significance level. It is well known that interruption of S1/S2 cleavage of S protein of SARS-CoV-2 by host protease and especially furin using inhibitors can prevent virus entry and pathogenesis [28] [29] [30] . Accordingly, furin becomes a good target for drug design against viral infections especially for the management of the COVID-19 pandemic threat currently [31] [32] . Considering the resolved crystal structure of furin, the in silico molecular docking experiments were conducted to predict the ability of small drugs approved for clinical application molecule as candidates for furin inhibitions and COVID-19 treatment. Furin inhibitors may bind as competitive inhibitors to the furin active site or to the interface cleft between P-domain and catalytic site as non-competitive inhibitors. Since the binding site seems to be a better target to followed precisely in this work our goal was to evaluate the binding capacity of the tested drugs to bind the enzyme active site and their potential was judged based on their binding energy (-ΔG). To perform docking experiments with furin as a receptor, we have selected the crystal structure of furin with PDB ID of 4OMD from the protein data bank. To assure that the crystal structure matches the wild-type protein and to survey their uniformity we performed pair-wise alignment on emboss sever (www.ebi.ac.uk/Tools/psa/emboss_needle/). As it is depicted in figure 3 , 4OMD structure completely has a wild-type structure without any mutations. Table 1 represents the docking results including the percent of binding site occupation and the binding energies in kJ/mol for our tested drugs. It should be mentioned that baloxavir, cefaclor, cefdinir, cefotaxime, cefuroxime, cefpodoxime, ciprofloxacin, emtricitabine, gemifloxacin, imipenem, moxifloxacin, ofloxacin, sulfamethoxazole, and tipranavir could not bind to enzyme active site (0%) and were deleted from further evaluation, however, they may bind to an allosteric site and inhibit the enzyme non-competitively. Table 1 also depicts the logP for drugs which reflects the hydrophobicity or bioavailability of drugs to reach the furin as a target for inhibition. A more hydrophobic drug (higher logP) with higher binding energy in accordance with more percent of binding site occupation will be a more effective inhibitor. In order to make a reasonable comparison we normalized the calculated percent of occupation, logP, and binding energies values and summate them in a cumulative index. Therefore the more cumulative index reveals the more effective drug for furin inhibition. As it is evident saquinavir, atazanavir, and nelfinavir, obtain a higher score in our series with 2.52, 2.16 and 2.13 cumulative indices respectively. In vitro assays show that nelfinavir exerts more effects on viral infections in contrast to ritonavir and lopinavir, while there is no experimental evidence regarding the clinical application of atazanavir and saquinavir in viral disease [33] [34] [35] . Among the rest anti viral drugs used in this work, ritonavir, indinavir, amprenavir, remdesivir, lopinavir, nevirapine, darunavir, lamivudine, oseltamivir, zidovudine and tenofovir with cumulative indices varies between 1.89 to 0.79 retain the next order of effectiveness in furin inhibition. It is important to notice that darunavir as example has least affinity to attacks furin active site while it can inhibits furin by binding to its allosteric site with high enough, this confirms our finding that darunavir only in 10 percent of probability binds furin active site [36] . Niclosamide among anti protozoa/parasite drugs used in this work include diminazene, niclosamide and ivermectin with cumulative index of 1.90 seems to be good candidate for clinical trail, however there are reports implying the positive effects of ivermectin in COVID-19 [40] [41] [42] . Table 3 summarize the interactions of drugs with their counterpart residues in furin active site analyzed and extracted on Protein-Ligand Interaction Profiler (PLIP) server (https://plip-tool.biotec.tu-dresden.de/plip-web/plip/index). This server gives all non-covalent interactions including, hydrogen bonds, salt bridges, hydrophobic interactions, and other non covalent interaction for docked complexes. As it is evident all the selected drugs make hydrophobic and hydrogen bond interactions with catalytic residues of the enzyme and their neighbors residues that make them effective binders for furin inhibition. Conclusion:Considering the maximum allowed dose for drugs and their cumulative indices as well as their biocompatibility and side effects we suggest erythromycin with more than 1gr/day and clarithromycin with 0.4gr/day among macrolides at the first step of disease as prophylaxis over niclosamide (with more than 1gr/day), and antiviral drugs of saquinavir, nelfinavir, indinavir and atazanavir with less than 1gr/day in sever state of the disease with the higher most cumulative indices. The biology and therapeutic targeting of the proprotein convertases Proprotein convertases in health and disease Therapeutic uses of furin and its inhibitors: a patent review A truncation allele in vascular endothelial growth factor c reveals distinct modes of signaling during lymphatic and vascular development Engineering protein processing of the mammary gland to produce abundant hemophilia B therapy in milk. Sci Rep Strategies for recombinant Furin employment in a biotechnological process: Complete target protein precursor cleavage Shed' furin: Mapping of the cleavage determinants and identification of its C-terminus The spike glycoprotein of the new coronavirus 2019-nCoV contains a furin-like cleavage site absent in CoV of the same clade A Multibasic Cleavage Site in the Spike Protein of SARS-CoV-2 Is Essential for Infection of Human Lung Cells Handbook of Proteolytic Enzymes X-ray Structures of Human Furin in Complex with Competitive Inhibitors The role of furin cleavage site in SARS-CoV-2 spike proteinmediated membrane fusion in the presence or absence of trypsin Structure of Furin Protease Binding to SARS-CoV-2 Spike Glycoprotein and Implications for Potential Targets and Virulence Proteolytic processing of Middle East respiratory syndrome coronavirus spikes expands virus tropism A unique protease cleavage site predicted in the spike protein of the novel pneumonia coronavirus (2019-nCoV) potentially related to viral transmissibility Furin Inhibitors Block SARS-CoV-2 Spike Protein Cleavage to Suppress Virus Production and Cytopathic Effects SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor Involvement of Spike Protein, Furin, and ACE2 in SARS-CoV-2 The Potential Role of Osteopontin and Furin in Worsening Disease Outcomes in COVID-19 Patients with Pre-Existing Diabetes. Cells Structure of the unliganded form of the proprotein convertase furin suggests activation by a substrate-induced mechanism The spike glycoprotein of the new coronavirus 2019-nCoV contains a furin-like cleavage site absent in CoV of the same clade The biology and therapeutic targeting of the proprotein convertases The spike glycoprotein of the new coronavirus 2019-nCoV contains a furin-like cleavage site absent in CoV of the same clade Homology modeling and molecular dynamics simulation of N-myristoyltransferase from protozoan parasites: active site characterization and insights into rational inhibitor design HexServer: an FFTbased protein docking server powered by graphics processors Docking Efficiency Comparison of Surflex Commercial Package and Arguslab, a Licensable Freeware Accelerating Protein-Protein Docking Correlations Using A Six-Dimensional Analytic FFT Generating Virtual computational chemistry laboratory -design and description The spike glycoprotein of the new coronavirus 2019-nCoV contains a furin-like cleavage site absent in CoV of the same clade Proteolytic processing of Middle East respiratory syndrome coronavirus spikes expands virus tropism A unique protease cleavage site predicted in the spike protein of the novel pneumonia coronavirus (2019-nCoV) potentially related to viral transmissibility A multibasic cleavage site in the spike protein of SARS-CoV-2 is essential for infection of human lung cells Furin: A Potential Therapeutic Target for COVID-19 iScience 23 Nelfinavir was predicted to be a potential inhibitor of 2019-nCov main protease by an integrative approach combining homology modelling, molecular docking and binding free energy calculation Erratum in: Lancet A Potent Drug against Coronavirus Infection: Insight from Molecular Docking Study, Arch Clin Infect Dis Discovery of an Allosteric Site on Furin, Contributing to Potent Inhibition: A Promising Therapeutic for the Anemia of Chronic Inflammation Are clarithromycin, azithromycin and their analogues effective in the treatment of COVID19? Repurposing of antibiotics for clinical management of COVID-19: a narrative review Old Drugs for a Newly Emerging Viral Disease, COVID-19: Bioinformatic Prospective Ivermectin: an award-winning drug with expected antiviral activity against COVID-19 Ivermectin to prevent hospitalizations in patients with COVID-19 (IVERCOR-COVID19): a structured summary of a study protocol for a randomized controlled trial A COVID-19 prophylaxis? Lower incidence associated with prophylactic administration of ivermectin The author would like to express his thanks to the vice-chancellor of research and technology of the Shahid Chamran University of Ahvaz for providing financial support for this study under Research Grant No: SCU.SB99.477.