key: cord-0789897-jzu20vzd
authors: Zaragoza-Huesca, David; Martínez-Cortés, Carlos; Banegas-Luna, Antonio Jesús; Pérez-Garrido, Alfonso; Vegara-Meseguer, Josefina María; Peñas-Martínez, Julia; Rodenas, Maria Carmen; Espín, Salvador; Pérez-Sánchez, Horacio; Martínez-Martínez, Irene
title: Identification of Kukoamine A, Zeaxanthin, and Clexane as New Furin Inhibitors
date: 2022-03-03
journal: Int J Mol Sci
DOI: 10.3390/ijms23052796
sha: 1a93cb18bf1a16abc9f320593ded270ef471c4a7
doc_id: 789897
cord_uid: jzu20vzd

The endogenous protease furin is a key protein in many different diseases, such as cancer and infections. For this reason, a wide range of studies has focused on targeting furin from a therapeutic point of view. Our main objective consisted of identifying new compounds that could enlarge the furin inhibitor arsenal; secondarily, we assayed their adjuvant effect in combination with a known furin inhibitor, CMK, which avoids the SARS-CoV-2 S protein cleavage by means of that inhibition. Virtual screening was carried out to identify potential furin inhibitors. The inhibition of physiological and purified recombinant furin by screening selected compounds, Clexane, and these drugs in combination with CMK was assayed in fluorogenic tests by using a specific furin substrate. The effects of the selected inhibitors from virtual screening on cell viability (293T HEK cell line) were assayed by means of flow cytometry. Through virtual screening, Zeaxanthin and Kukoamine A were selected as the main potential furin inhibitors. In fluorogenic assays, these two compounds and Clexane inhibited both physiological and recombinant furin in a dose-dependent way. In addition, these compounds increased physiological furin inhibition by CMK, showing an adjuvant effect. In conclusion, we identified Kukoamine A, Zeaxanthin, and Clexane as new furin inhibitors. In addition, these drugs were able to increase furin inhibition by CMK, so they could also increase its efficiency when avoiding S protein proteolysis, which is essential for SARS-CoV-2 cell infection.

Furin is a proprotein convertase whose main function consists of cleaving zymogens that, after cleavage, acquire biological active function [1] . The furin cleavage site in substrates presents arginine at the first and fourth amino acids from the N-terminal region [2] . Among these substrates, we highlight growth factors, cytokines, coagulation proteins, albumin, hormones, metalloproteinases, and many types of receptors [3] . Since furin acts on a wide variety of substrates, it regulates many different processes. For example, furin regulates pancreatic granular acidification and, consequently, it also affects insulin secretion [4] . This enzyme is also implicated in proliferation and apoptosis processes, as in granulosa cells from ovaries in rats [5] . Even during embryonic development, furin plays an important role, since it increases the migration and expansion of human trophoblast cells [6] . With these obtained QSAR models, consensus model virtual screening was carried out against the compounds of natural origin from the FoodBank database (FDB, https: //foodb.ca/, accessed on 1 July 2020) to extract possible candidates as furin inhibitors. Table 3 shows the seven most active compounds. 

Among the compounds highlighted by virtual screening, just Kukoamine A (Cat#CFN93215; ChemFaces, Wuhan, Hubei, China) and Zeaxanthin (Cat#Q444; AK Scientific, Inc., Union City, CA, USA) inhibited both recombinant and physiological furin, and they did so in a dosedependent manner. In the case of Kukoamine A, the IC50 for recombinant furin was 1.07 mM (95% CI; 0.817-4.577), whereas the IC50 for physiological furin was 193.2 µM (95% CI; 165.5-229) ( Figure 1A ). With respect to Zeaxanthin, the IC50 for recombinant furin was 90.55 µM (95% CI; 64.87-753), whereas the IC50 for physiological furin was 49.55 µM (95% CI; 38.85-64.52) ( Figure 1B ). In addition to the compounds selected by virtual screening, furin inhibition by Clexane (LMWH) was also tested. The IC50 for recombinant furin was 31.648 µM (95% CI; 10.08-53.22), whereas the IC50 for physiological furin was 6.45 µM (95% CI; 6.142-6.784) ( Figure 1C ). 

Among the compounds highlighted by virtual screening, just Kukoamine A (Cat#CFN93215; ChemFaces, Wuhan, Hubei, China) and Zeaxanthin (Cat#Q444; AK Scientific, Inc., Union City, CA, USA) inhibited both recombinant and physiological furin, and they did so in a dose-dependent manner. In the case of Kukoamine A, the IC50 for recombinant furin was 1.07 mM (95% CI; 0.817-4.577), whereas the IC50 for physiological furin was 193.2 µM (95% CI; 165.5-229) ( Figure 1A ). With respect to Zeaxanthin, the IC50 for recombinant furin was 90.55 µM (95% CI; 64.87-753), whereas the IC50 for physiological furin was 49.55 µM (95% CI; 38.85-64.52) ( Figure 1B ). In addition to the compounds selected by virtual screening, furin inhibition by Clexane (LMWH) was also tested. The IC50 for recombinant furin was 31.648 µM (95% CI; 10.08-53.22), whereas the IC50 for physiological furin was 6.45 µM (95% CI; 6.142-6.784) ( Figure 1C ). 

In addition to the effect of the single compounds on the inhibition of furin, we decided to test their coadjuvants with CMK, a known furin inhibitor that also avoids SARS-CoV-2 S protein cleavage [27] . The addition of Zeaxanthin, Kukoamine A, or Clexane IC50 to CMK increased physiological furin inhibition ( Figure 2 ). The IC50 for CMK alone was 4.205 µM (95% CI; 2.991-6.002), whereas the addition of Zexanthin, Kukoamine A or Clexane reduced the IC50 to 1.488 µM (95% CI; 0.885-2.501), 0.749 µM (95% CI; 0.552-0.998) or 0.096 µM (95% CI; 0.037-0.160), respectively. This showed an adjuvant effect on physiological furin inhibition using CMK in combination with Zeaxanthin, Kukoamine A, or Clexane. 

In addition to the effect of the single compounds on the inhibition of furin, we decided to test their coadjuvants with CMK, a known furin inhibitor that also avoids SARS-CoV-2 S protein cleavage [27] . The addition of Zeaxanthin, Kukoamine A, or Clexane IC50 to CMK increased physiological furin inhibition ( Figure 2 ). The IC50 for CMK alone was 4.205 µM (95% CI; 2.991-6.002), whereas the addition of Zexanthin, Kukoamine A or Clexane reduced the IC50 to 1.488 µM (95% CI; 0.885-2.501), 0.749 µM (95% CI; 0.552-0.998) or 0.096 µM (95% CI; 0.037-0.160), respectively. This showed an adjuvant effect on physiological furin inhibition using CMK in combination with Zeaxanthin, Kukoamine A, or Clexane. 

As the effects of CMK and heparin have already been tested on human cell lines [27, 30] , we decided to check if Zeaxanthin and Kukoamine A could alter cell viability at tested concentrations for physiological furin inhibition. Kukoamine A significantly reduced 293T cells' viability at double IC50 and full IC50 concentrations ( Figure 3 ). However, a half IC50 concentration did not affect cell viability, as the percentage of live cells was similar to cells in the absence of any compound (87.55% ± 1.48% (half of the IC50) vs. 94.25% ± 3.18% (no Kukoamine A)). In the case of Zeaxanthin, a half IC50 concentration caused a slight reduction in cell viability in comparison to cells in the absence of bioactive compounds (81.1% ± 4.1% (half of the IC50) vs 94.25% ± 3.18% (no Zeaxanthin)) ( Figure  3 ). The Zeaxanthin IC50 concentration reduced cell viability by 61%, whereas double the IC50 concentration reduced it to nearly total death. 

As the effects of CMK and heparin have already been tested on human cell lines [27, 30] , we decided to check if Zeaxanthin and Kukoamine A could alter cell viability at tested concentrations for physiological furin inhibition. Kukoamine A significantly reduced 293T cells' viability at double IC50 and full IC50 concentrations ( Figure 3 ). However, a half IC50 concentration did not affect cell viability, as the percentage of live cells was similar to cells in the absence of any compound (87.55% ± 1.48% (half of the IC50) vs. 94.25% ± 3.18% (no Kukoamine A)). In the case of Zeaxanthin, a half IC50 concentration caused a slight reduction in cell viability in comparison to cells in the absence of bioactive compounds (81.1% ± 4.1% (half of the IC50) vs 94.25% ± 3.18% (no Zeaxanthin)) ( Figure 3 ). The Zeaxanthin IC50 concentration reduced cell viability by 61%, whereas double the IC50 concentration reduced it to nearly total death. 

Furin inhibitors that were predicted by virtual screening and that were experimentally confirmed were submitted to molecular docking calculations, and the resulting poses can be observed in Figures 4 and 5 for Zeaxanthin and Kukoamine A, respectively. It can be appreciated that both compounds extensively blocked access to the active site by means of hydrogen bonds and hydrophobic interactions, thus inhibiting furin's catalytic action. 

Furin inhibitors that were predicted by virtual screening and that were experimentally confirmed were submitted to molecular docking calculations, and the resulting poses can be observed in Figures 4 and 5 for Zeaxanthin and Kukoamine A, respectively. It can be appreciated that both compounds extensively blocked access to the active site by means of hydrogen bonds and hydrophobic interactions, thus inhibiting furin's catalytic action. 

Proprotein convertase furin is a very versatile protease. It exhibits different functions that become essential in both health and disease. Among the pathophysiological processes in which furin has a crucial role, there are many types of cancer where furin cleaves and activates a wide range of proteins that promote tumor phenotypes [31] . 

Proprotein convertase furin is a very versatile protease. It exhibits different functions that become essential in both health and disease. Among the pathophysiological processes in which furin has a crucial role, there are many types of cancer where furin cleaves and activates a wide range of proteins that promote tumor phenotypes [31] . However, due to its heterogeneous functions, furin represses the tumorigenic properties of some cancer cells and its inhibition can lead to aggressive phenotypes in other tumors [31] . Furin is also involved in some infectious diseases. Viral substrates have gained importance among furin targets as they become activated as infectious proteins when they are cleaved [1] . This relates to viruses from different families, involving infections provoked by Ebola [32] , Influenza A [33] , or Metapneumovirus [34] .

In an endeavor to propose furin inhibition as a possible therapeutic target, there have been projects aimed at studying the biological effects of this inhibition in different pathological models. For example, furin autoinhibitory propeptide has been used to reduce metalloproteinase-9 activity in vitro in breast cancer [35, 36] . In clinical trials concerning patients with Ewing's sarcoma cancer, one of the strategies consisted of silencing furin by short hairpin RNAs [37] . In a murine model infected with Pseudomonas aeruginosa, nona-D-arginine-mediated furin inhibition diminished corneal adverse effects [38] . In the recent pandemic context, furin inhibition has also been tested as a possible preventive strategy against infection with SARS-CoV-2, by the usage of drugs such as CMK [27] .

In this study, we aimed to identify new furin inhibitors that could be added to the arsenal of compounds that target furin from a therapeutic point of view. From a virtual screening, we selected a series of candidate compounds, which were tested for the inhibition of purified recombinant and physiological secreted furin. Among these compounds, we selected Zeaxanthin and Kukoamine A, since they were the most effective inhibitors of both physiological and recombinant furin. In addition to these compounds, Clexane, a LMWH, directly inhibited purified recombinant furin but also furin substrate proteolysis in the U-251 MG secretome. These results support that Clexane can directly act against furin. Its IC50 was lower than those of Kukoamine A or Zeaxanthin. Zeaxanthin, as with other carotenoids, is highlighted for its antioxidant and preventive properties against cardiovascular or ocular diseases as well as cancer [39] . In the context of viral infections, lower levels of Zeaxanthin have been reported in the serum of patients infected by the human immunodeficiency virus (HIV), in comparison with noninfected subjects [40] . Although its IC50 value reduced the viability of healthy cells tested in our assay, Zeaxanthin is an FDA-approved drug for which concentrations of up to 100 µM have been administered daily in clinical trials in patients with age-related macular degeneration, with no adverse effects reported [41, 42] . In the case of Kukoamine A, some studies support its antitumor properties in glioblastoma [43] and its neuroprotective effects [44] . Similarly, Zeaxanthin's and Kukoamine A's IC50 values considerably reduced 293T HEK cell viability. However, some assays that administered Kukoamine A to rats or mice in doses between 5 and 20 mg/kg highlighted only its protective effects, with no significant reported adverse effects [45, 46] . Clexane is a wellknown venous thromboembolism prophylactic drug [47] , but despite its reported abilities to inhibit some furin-mediated pathways, a direct effect of LMWH on this enzyme has not been described [29] . However, our results with purified recombinant furin and Clexane could tip the scale towards a direct interaction. It is important to note that, although we have only tested Clexane, similar effects would be expected with other LMWHs.

Secondarily, we showed that Kukoamine A, Zeaxanthin, and Clexane increased CMK efficiency over furin, and CMK avoided SARS-CoV-2 S protein proteolysis through furin inhibition [27] . These results suggest that these new inhibitors could be considered, individually or in combination, as potential drugs against COVID-19. In this context, these compounds would bypass the immune escape carried out by SARS-CoV-2 variants, since they act on furin and not on the virus. However, future and extensive studies should be undertaken in order to examine these potential applications. In any case, drug coadjuvants are of vital importance, especially when dealing with infectious diseases, since excessive use of single-compound-based therapies is potentially responsible for the drug-resistant generation of many types of pathogens [48] [49] [50] .

The main limitations of this study are related to the high IC50s for the furin inhibitors proposed by virtual screening, taking into account that the QSAR models used to select the compounds were built considering active compounds with IC50/Ki values lower than 0.1 uM. However, there are clinical trials of our selected compounds using concentrations similar to those tested in our experiments. On the other hand, although the selected compounds could potentially be used against COVID-19, it would be necessary to test in clinical trials whether they could prevent contagion and whether they would also be effective in reducing viral replication once contagion has occurred. In addition, we do not know the adverse effects of these drugs in this context. In the specific case of heparin, there are numerous adverse effects described in addition to its contraindication with the use of other drugs [51] .

In conclusion, we have identified Kukoamine A, a natural origin compound; Zeaxanthin, an FDA-approved drug; and Clexane, a known antithrombotic compound, as new furin direct inhibitors. In addition, Clexane, Zeaxanthin, and Kukoamine A are able to increase furin inhibition by CMK, so they could also increase its efficiency when avoiding S protein proteolysis and could be tested in clinical trials for COVID-19 prevention.

The chemical substances, as well as the activity values (Ki, IC50, and % inhibition), were extracted from the CHEMBL database (https://www.ebi.ac.uk/chembl, accessed on 1 July 2020), for a total of 148 chemicals. All data, both chemicals and activities, were curated for the three activity sets. The data curation was carried out to reduce misannotation and rounding unit errors [52] . The chemical dataset was curated following the protocols proposed by Fourches et al. [53, 54] . These protocols included the structural normalization of specific chemotypes, such as aromatic and nitro groups, and the removal of inorganic salts, organometallic compounds, etc. Standardizer was used for chemical curation (Standardizer 17.21.0, ChemAxon https://www.chemaxon.com, Budapest, Hungary, accessed on 1 July 2020). To select the endpoint, the following method was carried out: we used data with values of Ki as the main dataset, using a threshold of pKi ≥ 7 to select actives (45 actives and 66 inactives). Next, we included chemicals with pIC50 values that were not in the main dataset, choosing a threshold of pIC50 ≥ 7 to select active compounds (3 actives and 10 inactives were included). Then, we included chemicals that were not present in the main dataset using % inhibition ≥ 50% as the active selection threshold (0 actives and 7 inactives were included).

More than 5000 different descriptors were calculated for the 131 compounds using DRAGON software [55] . The descriptors with low variance or that were highly correlated with each other (r 2 > 0.95) were removed. As the QSAR models were based on the principle of similarity (similar substances present similar activities), and to avoid heterogeneities in the data, the so-called activity cliffs (similar substances with very different activity) and outliers (substances very different from the rest of the compounds) were removed. The activity cliff elimination of inactive chemicals was conducted using the Castillo-González et al. [56] method with a threshold value of 101.21. The outlier selection was carried out with the applicability domain (AD) protocol defined by Melagraki et al. [57, 58] with Z = 0.5 (see applicability domain section). Once the analysis was completed, 18 compounds were eliminated (15 activity cliffs and 3 outliers). To ensure that the data will result in effective QSAR modeling, they must present a MODI value [59] greater than 0.65; for our data, the value obtained from MODI was 0.76.

The 113 chemicals selected were randomly divided into three subsets: training (80%), test (10%), and external (10%). With these three subsets, we used linear discriminant analysis for the QSAR model development and a genetic algorithm for the descriptor selection technique using the ROCED [60] parameter as a fitness function. We computed models using between 2 and 11 descriptors. Moreover, selectivity, specificity, correct classification rate (CCR), and area under the ROC curve (AUC) were calculated for training, test, and external sets.

In addition, the leave-one-out (LOO) approach was conducted on the training set to assess the internal predictivity. In this approach, one compound of the training set was omitted, and the statistical parameters were recalculated with the remaining substances. This process was conducted for all training sets, and the values of the sensitivity, specificity, AUC, and ROCED were reported. High values of sensitivity and specificity and an AUC with a low value for ROCED are indicative of a model's robustness.

To avoid models producing good classification due to chance, the Y-randomization technique was carried out. The activity values of the training set were randomized (Yrnd). This process was conducted 300 times and the average values of the sensitivity (Sensrnd), specificity (Sprnd), and ROCEDrnd for both training and test sets, as well as the AUC values, were reported. A high average in ROCED and a low average in sensitivity, specificity, and AUC are indicative of a good QSAR model.

To obtain the most statistically robust and predictive models, we employed the combinatorial QSAR strategy. For this purpose, the most different models based on the canonical measure of distance (CMD) [61] were chosen to preserve the most information and diversity and to construct a consensus model.

In this study, we defined AD as a distance threshold ∆T between a compound of interest and its nearest neighbors of the set considered. It was calculated as

where y and σ are the mean and the standard deviation of the distances that are below the mean of the distance matrix, and Z is defined at 0.5. For outlier exclusion, we defined the AD in the entire descriptor space, and a chemical that presented a Euclidean distance with its nearest neighbor that was higher than this value was considered as an outlier. For chemical prediction in the virtual screening test, we defined the AD in the selected descriptor space of the model, and for chemicals that presented a Euclidean distance with their nearest neighbors in the training set higher than this value, no prediction was made.

Compound class (i.e., active or inactive) assignment in the virtual screening was based on the majority vote across the independent models developed with one condition: the chemical must be inside the applicability domain defined above. To prioritize compounds with the same activity vote, we sorted them by the mean value of probability given by each model.

The human glioblastoma cell line U-251 MG expresses human physiological furin, which can be secreted or retained intracellularly [62] . The U-251 MG cells (European Collection for Authenticated Cell Cultures) were first cultured in Corning ® T-75 flasks (catalog #430641; ThermoFisher Scientific, Waltham, MA, USA) with DMEM (Dulbecco's Modified Eagle's Medium) containing 4.5 g/L glucose (Gibco Thermo Fisher, Madrid, Spain), supplemented with 10% fetal bovine serum (Gibco Thermo Fisher, Madrid, Spain), 1% GlutaMax (Gibco Thermo Fisher, Madrid, Spain), 1% nonessential amino acids (Gibco Thermo Fisher, Madrid, Spain) and 0.1% gentamicin (Gibco Thermo Fisher, Madrid, Spain). When 100% confluence was reached, the medium was recovered and centrifuged for 5 at 280 g to discard cell debris, and finally, supernatant (secretome) was collected. The activity of secreted furin was confirmed by using a fluorogenic specific substrate (pERTKR-AMC fluorogenic peptide substrate; ES013; LOT#PYO02; R&D Systems, Minneapolis, MN, USA) ( Figure A1 ). These results were achieved by incubating U-251 MG secretome concentrated 10X with 50 µM furin substrate in a 96-well black plate (10030581, ThermoFisher Scientific, Waltham, MA, USA), and then, measuring fluorescence in a Multiskan Go (ThermoFisher Scientific, Waltham, MA, USA) plate reader at 380 and 460 nm for emission and excitation wavelengths, respectively, at 37 • C for 30 min.

Clexane (LMWH, enoxaparin sodium, Clexane, Sanofi Aventis S.A., Barcelona, Spain) and each candidate compound selected by virtual screening were assayed as inhibitors of human recombinant furin (1503-SE-010; Lot#INK2320031; R&D Systems) and physiological furin secreted from U-251 MG cells. For both types of furin, the fluorogenic specific substrate (mentioned above) was used. In the inhibition assays, activity buffer (Tris Base 25 mM (BP152-1; LOT#165920; Fisher BioReagents), CaCl2 1 mM and Brij- 35 (Cat. No. 20150 ; ThermoFisher Scientific, Waltham, MA, USA) 0,5% (w/v), pH 9) was the medium for recombinant furin, whereas distilled water was used for U-251 MG secretome since it showed less furin activity in the presence of activity buffer. Compounds were reconstituted in their appropriate solvents (distilled water, ethanol, PBS, etc.), according to the supplier's instructions. In a 96-well black plate, U-251 MG secretome concentrated 10× or 3.6 nM human recombinant furin were incubated with different compound concentrations for 10 min at room temperature, and then, 50 µM furin specific fluorogenic substrate was added to every well, except one blank. Fluorescence was recorded in the Multiskan Go plate reader at 380 and 460 nm for emission and excitation wavelengths, respectively, at 37 • C for 30 min.

The IC50 was calculated only for those compounds that inhibited recombinant and physiological furin. For each compound concentration, furin activity's maximum velocity (mRFU/min) was recorded. Then, by using GraphPad software (v. 8), the correlation between the decimal logarithm (compound concentration) and furin activity's maximum velocity was calculated by nonlinear regression (curve fit) for XY analyses. From this statistical method, the IC50 was calculated with a 95% confidence interval (CI).

Since furin inhibition by CMK had been tested only in physiological conditions [27] , adjuvant inhibition assays were performed with furin from the U-251 MG secretome. Inhibition was assessed by using different concentrations of CMK (Cat.No. B5437; APExBIO) and the IC50 concentrations of selected virtual screening compounds or Clexane. In a 96-well black plate, various CMK (0-50 µM) concentrations were incubated with the IC50 concentrations of the rest of the compounds and with 10× U-251 MG secretome as described above. After 10 min incubation at room temperature, 50 µM furin specific fluorogenic substrate was added and the plate was read in a Multiskan Go plate reader at 380 and 460 nm for emission and excitation wavelengths, respectively, at 37 • C for 30 min.

New furin inhibitors identified by virtual screening were tested in cell viability assays. The 293T HEK cell line (human embryonic kidney cell line purchased from American Type Culture Collection) was used as healthy cell model. They were cultured in DMEM containing 4,5 g/L glucose, supplemented with 10% fetal bovine serum, 1% GlutaMax, and 0,1% gentamicin. They were maintained in Corning ® T-75 flasks (catalog #430641; ThermoFisher Scientific, Waltham, MA, USA) at 37 • C and 5% carbon dioxide. At 100% cellular confluence, cells were subcultured in 12-well plates (flat bottom, sterile, NUNC brand, Biolab, Madrid, Spain) at 100,000 cells/mL for 24 h. Then, different concentrations (double, half, and full IC50 concentrations for physiological furin) of furin inhibitors were added to cells for 24 h, and cells were collected afterward. The experiment was run in duplicates. Then, the cell pellet was resuspended in running buffer (autoMACS ® Running Buffer-MACS Separation Buffer; order no.:130-091-221; Miltenyi Biotec, Madrid, Spain) and 7-AAD (Cat: 51-68981E; BD Biosciences, Madrid, Spain). Finally, this solution was evaluated in a flow cytometer (BD Accuri™ C6 Plus Flow Cytometer; BD Biosciences, Madrid, Spain) to calculate mortality rates. A blank solution with cell pellet resuspended in running buffer was used to discard cell autofluorescence. The results were analyzed in FlowJo software (v. 10).

To provide mechanistic insights into the structure and main interactions established between furin and characterized inhibitors derived from virtual screening, molecular docking calculations were carried out.

The selected inhibitors were set up for docking simulations using AmberTools (AMBER 2017, University of California, San Francisco, CA, USA) [63] . Molecular parameters were calculated by computing partial charges by the MMFF94 force field, by adding hydrogen atoms, and by minimizing energies (default parameters) [64] .

The crystal structure of furin (Protein Data Bank code 5MIM) was used to build the protein model system. At an early stage, bond orders were assigned, hydrogens were added, and cap termini were included with the Protein Preparation Wizard module as implemented in Maestro (Schrödinger Release 2021-2: Maestro, Schrödinger, LLC, New York, NY, USA) [65] . Protonation states of all side chains were subsequently defined using PROPKA3.1. Partial charges over all atoms were finally assigned within the AMBER99 force field scheme as implemented in AmberTools. Docking simulations were performed with Lead Finder software v1.1.20 [66] via MetaScreener (https://github.com/bio-hpc/ metascreener, accessed on 20 December 2021) in the coordinates of the active site of the protein. All docking parameters were set to default for the calculations. The best-ranked docking score pose for every compound was retained for further analysis.

Two patents in connection with these findings were filed in Spain, 

Furin-mediated protein processing in infectious diseases and cancer

Furin: A Potential Therapeutic Target for COVID-19

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

Role of furin in granular acidification in the endocrine pancreas: Identification of the V-ATPase subunit Ac45 as a candidate substrate

Proprotein convertase furin regulates apoptosis and proliferation of granulosa cells in the rat ovary

The proprotein convertase furin in human trophoblast: Possible role in promoting trophoblast cell migration and invasion

The proprotein convertase furin in tumour progression: The PC furin in tumour progression

Elevated furin expression in aggressive human head and neck tumors and tumor cell lines

Furin inhibition results in absent or decreased invasiveness and tumorigenicity of human cancer cells

Increased furin activity enhances the malignant phenotype of human head and neck cancer cells

Human carcinoma cell growth and invasiveness is impaired by the propeptide of the ubiquitous proprotein convertase furin

Inhibition of proprotein convertases is associated with loss of growth and tumorigenicity of HT-29 human colon carcinoma cells: Importance of insulin-like growth factor-1 (IGF-1) receptor processing in IGF-1-mediated functions

Furin-mediated protein processing in infectious diseases and cancer

Furin: A Potential Therapeutic Target for COVID-19

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

Role of furin in granular acidification in the endocrine pancreas: Identification of the V-ATPase subunit Ac45 as a candidate substrate

Proprotein convertase furin regulates apoptosis and proliferation of granulosa cells in the rat ovary

The proprotein convertase furin in human trophoblast: Possible role in promoting trophoblast cell migration and invasion

The proprotein convertase furin in tumour progression: The PC furin in tumour progression

Elevated furin expression in aggressive human head and neck tumors and tumor cell lines

Furin inhibition results in absent or decreased invasiveness and tumorigenicity of human cancer cells

Increased furin activity enhances the malignant phenotype of human head and neck cancer cells

Human carcinoma cell growth and invasiveness is impaired by the propeptide of the ubiquitous proprotein convertase furin

Inhibition of proprotein convertases is associated with loss of growth and tumorigenicity of HT-29 human colon carcinoma cells: Importance of insulin-like growth factor-1 (IGF-1) receptor processing in IGF-1-mediated functions

Selective inhibition of proprotein convertases represses the metastatic potential of human colorectal tumor cells

Proteolytic activation of bacterial toxins: Role of bacterial and host cell proteases

Evidence for involvement of furin in cleavage and activation of diphtheria toxin

Cell-mediated cleavage of Pseudomonas exotoxin between Arg279 and Gly280 generates the enzymatically active fragment which translocates to the cytosol

Human cytomegalovirus strain Towne glycoprotein B is processed by proteolytic cleavage

Partial amino terminal sequencing of spike polypeptide S2 identifies the sequence Arg-Arg-Phe-Arg-Arg at the cleavage site of the spike precursor propolypeptide of IBV strains Beaudette and M41

Mechanism of Borna disease virus entry into cells

Biochemical features and mutations of key proteins in SARS-CoV-2 and their impacts on RNA therapeutics

Fragmented health systems in COVID-19: Rectifying the misalignment between global health security and universal health coverage

Gut Microbiota and COVID-19

Deploying Biomolecules as Anti-COVID-19 Agents

Weekly Epidemiological Update on COVID-19-25

Waves and variants of SARS-CoV-2: Understanding the causes and effect of the COVID-19 catastrophe

The old but new: Can unfractioned heparin and low molecular weight heparins inhibit proteolytic activation and cellular internalization of SARS-CoV2 by inhibition of host cell proteases?

Furin Inhibitors Block SARS-CoV-2 Spike Protein Cleavage to Suppress Virus Production and Cytopathic Effects

Alpha1-Antitrypsin Portland, a bioengineered serpin highly selective for furin: Application as an antipathogenic agent

Heparin inhibits activation of latent transforming growth factor-β1

Heparin prevents programmed cell death in human trophoblast

The proprotein convertase furin in cancer: More than an oncogene

Processing of the Ebola virus glycoprotein by the proprotein convertase furin

Molecular analyses of the hemagglutinin genes of H5 influenza viruses: Origin of a virulent turkey strain

Modulation of protease activated receptor 1 influences human metapneumovirus disease severity in a mouse model

The prosegments of furin and PC7 as potent inhibitors of proprotein convertases. In vitro and ex vivo assessment of their efficacy and selectivity

Opposing function of the proprotein convertases furin and PACE4 on breast cancer cells malignant phenotypes: Role of tissue inhibitors of metalloproteinase-1

Vigil + Irinotecan and Temozolomide in Ewing's Sarcoma (VITA)

Nona-D-arginine therapy for Pseudomonas aeruginosa keratitis

Carotenoids: Biochemistry, pharmacology and treatment

Improved antioxidant status among HIV-infected injecting drug users on potent antiretroviral therapy

Age-related macular degeneration: Effects of a short-term intervention with an oleaginous kale extract-a pilot study

Macular xanthophylls and ω-3 long-chain polyunsaturated fatty acids in age-related macular degeneration: A randomized trial

Kukoamine A inhibits human glioblastoma cell growth and migration through apoptosis induction and epithelial-mesenchymal transition attenuation

Kukoamine A Protects against NMDA-Induced Neurotoxicity Accompanied with Down-Regulation of GluN2B-Containing NMDA Receptors and Phosphorylation of PI3K

Signaling Pathway in Cultured Primary Cortical Neurons

Neuroprotective Effects of Kukoamine a against Radiation-induced Rat Brain Injury through Inhibition of Oxidative Stress and Neuronal Apoptosis

Kukoamine A attenuates insulin resistance and fatty liver through downregulation of Srebp-1c

Safety and effectiveness of aspirin and enoxaparin for venous thromboembolism prophylaxis after total hip and knee arthroplasty: A systematic review

HIV-1 Resistance Dynamics in Patients with Virologic Failure to Dolutegravir Maintenance Monotherapy

Antiviral drug resistance testing in patients with chronic hepatitis B

A5230 team. Drug susceptibility and resistance mutations after first-line failure in resource limited settings

Stat Pearls

Low Molecular Weight Heparin (LMWH)

The Experimental Uncertainty of Heterogeneous Public Ki Data

On the Importance of Chemical Structure Curation in Cheminformatics and QSAR Modeling Research

Curation of chemogenomics data

Molecular Descriptors for Chemoinformatics

Harmonization of QSAR Best Practices and Molecular Docking Provides an Efficient Virtual Screening Tool for Discovering New G-Quadruplex Ligands

In Silico Exploration for Identifying Structure-Activity Relationship of MEK Inhibition and Oral Bioavailability for Isothiazole Derivatives

Ligand-based virtual screening procedure for the prediction and the identification of novel β -amyloid aggregation inhibitors using Kohonen Maps and Counterpropagation Artificial Neural Networks

Data Set Modelability by QSAR

Two new parameters based on distances in a receiver operating characteristic chart for the selection of classification models

Canonical measure of correlation (CMC) and canonical measure of distance (CMD) between sets of data. Part 1. Theory and simple chemometric applications

Processing of immunosuppressive pro-TGF-beta 1,2 by human glioblastoma cells involves cytoplasmic and secreted furin-like proteases

Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments

Lead finder: An approach to improve accuracy of protein-ligand docking, binding energy estimation, and virtual screening

Acknowledgments: Supercomputing resources in this work were supported by the Poznan Supercomputing Center's infrastructures, the e-infrastructure program of the Research Council of Norway,