Inhibition of prohormone convertases PC1/3 and PC2 by 2,5-dideoxystreptamine derivatives

The Role of Proprotein Convertases in Upper Airway Remodeling

Chronic rhinosinusitis (CRS) is a multifactorial, heterogeneous disease characterized by persistent inflammation of the sinonasal mucosa and tissue remodeling, which can include basal/progenitor cell hyperplasia, goblet cell hyperplasia, squamous cell metaplasia, loss or dysfunction of ciliated cells, and increased matrix deposition. Repeated injuries can stimulate airway epithelial cells to produce inflammatory mediators that activate epithelial cells, immune cells, or the epithelial-mesenchymal trophic unit. This persistent inflammation can consequently induce aberrant tissue remodeling. However, the molecular mechanisms driving disease within the different molecular CRS subtypes remain inadequately characterized. Numerous secreted and cell surface proteins relevant to airway inflammation and remodeling are initially synthesized as inactive precursor proteins, including growth/differentiation factors and their associated receptors, enzymes, adhesion molecules, neuropeptides, and peptide hormones. Therefore, these precursor proteins require post-translational cleavage by proprotein convertases (PCs) to become fully functional. In this review, we summarize the roles of PCs in CRS-associated tissue remodeling and discuss the therapeutic potential of targeting PCs for CRS treatment.

Altered islet prohormone processing: A cause or consequence of diabetes?

Peptide hormones are first produced as larger precursor prohormones that require endoproteolytic cleavage to liberate the mature hormones. A structurally conserved but functionally distinct family of nine prohormone convertase enzymes (PCs) are responsible for cleavage of protein precursors of which PC1/3 and PC2 are known to be exclusive to neuroendocrine cells and responsible for prohormone cleavage. Differential expression of PCs within tissues define prohormone processing; whereas glucagon is the major product liberated from proglucagon via PC2 in pancreatic α-cells, proglucagon is preferentially processed by PC1/3 in intestinal L cells to produce glucagon-like peptides 1 and 2 (GLP-1, GLP-2). Beyond our understanding of processing of islet prohormones in healthy islets, there is convincing evidence that proinsulin, proIAPP, and proglucagon processing is altered during prediabetes and diabetes. There is predictive value of elevated circulating proinsulin or proinsulin : C-peptide ratio for progression to type 2 diabetes and elevated proinsulin or proinsulin : C-peptide is predictive for development of type 1 diabetes in at risk groups. After onset of diabetes, patients have elevated circulating proinsulin and proIAPP and proinsulin may be an autoantigen in type 1 diabetes. Further, preclinical studies reveal that α-cells have altered proglucagon processing during diabetes leading to increased GLP-1 production. We conclude that despite strong associative data, current evidence is inconclusive on the potential causal role of impaired prohormone processing in diabetes, and suggest that future work should focus on resolving the question of whether altered prohormone processing is a causal driver or merely a consequence of diabetes pathology.

Mouse Models of Human Proprotein Convertase Insufficiency

The kexin-like proprotein convertases perform the initial proteolytic cleavages that ultimately generate a variety of different mature peptide and proteins, ranging from brain neuropeptides to endocrine peptide hormones, to structural proteins, among others. In this review, we present a general introduction to proprotein convertase structure and biochemistry, followed by a comprehensive discussion of each member of the kexin-like subfamily of proprotein convertases. We summarize current knowledge of human proprotein convertase insufficiency syndromes, including genome-wide analyses of convertase polymorphisms, and compare these to convertase null and mutant mouse models. These mouse models have illuminated our understanding of the roles specific convertases play in human disease and have led to the identification of convertase-specific substrates; for example, the identification of procorin as a specific PACE4 substrate in the heart. We also discuss the limitations of mouse null models in interpreting human disease, such as differential precursor cleavage due to species-specific sequence differences, and the challenges presented by functional redundancy among convertases in attempting to assign specific cleavages and/or physiological roles. However, in most cases, knockout mouse models have added substantively both to our knowledge of diseases caused by human proprotein convertase insufficiency and to our appreciation of their normal physiological roles, as clearly seen in the case of the furin, proprotein convertase 1/3, and proprotein convertase 5/6 mouse models. The creation of more sophisticated mouse models with tissue- or temporally-restricted expression of specific convertases will improve our understanding of human proprotein convertase insufficiency and potentially provide support for the emerging concept of therapeutic inhibition of convertases.

Attacking COVID-19 Progression Using Multi-Drug Therapy for Synergetic Target Engagement

COVID-19 is a devastating respiratory and inflammatory illness caused by a new coronavirus that is rapidly spreading throughout the human population. Over the past 12 months, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus responsible for COVID-19, has already infected over 160 million (>20% located in United States) and killed more than 3.3 million people around the world (>20% deaths in USA). As we face one of the most challenging times in our recent history, there is an urgent need to identify drug candidates that can attack SARS-CoV-2 on multiple fronts. We have therefore initiated a computational dynamics drug pipeline using molecular modeling, structure simulation, docking and machine learning models to predict the inhibitory activity of several million compounds against two essential SARS-CoV-2 viral proteins and their host protein interactors-S/Ace2, Tmprss2, Cathepsins L and K, and Mpro-to prevent binding, membrane fusion and replication of the virus, respectively. All together, we generated an ensemble of structural conformations that increase high-quality docking outcomes to screen over >6 million compounds including all FDA-approved drugs, drugs under clinical trial (>3000) and an additional >30 million selected chemotypes from fragment libraries. Our results yielded an initial set of 350 high-value compounds from both new and FDA-approved compounds that can now be tested experimentally in appropriate biological model systems. We anticipate that our results will initiate screening campaigns and accelerate the discovery of COVID-19 treatments.

A Virtual Screening Platform Identifies Chloroethylagelastatin A as a Potential Ribosomal Inhibitor

Chloroethylagelastatin A (CEAA) is an analogue of agelastatin A (AA), a natural alkaloid derived from a marine sponge. It is under development for therapeutic use against brain tumors as it has excellent central nervous system (CNS) penetration and pre-clinical therapeutic activity against brain tumors. Recently, AA was shown to inhibit protein synthesis by binding to the ribosomal A-site. In this study, we developed a novel virtual screening platform to perform a comprehensive screening of various AA analogues showing that AA analogues with proven therapeutic activity including CEAA have significant ribosomal binding capacity whereas therapeutically inactive analogues show poor ribosomal binding and revealing structural fingerprint features essential for drug-ribosome interactions. In particular, CEAA was found to have greater ribosomal binding capacity than AA. Biological tests showed that CEAA binds the ribosome and contributes to protein synthesis inhibition. Our findings suggest that CEAA may possess ribosomal inhibitor activity and that our virtual screening platform may be a useful tool in discovery and development of novel ribosomal inhibitors.

Revisiting Proinsulin Processing: Evidence That Human β-Cells Process Proinsulin With Prohormone Convertase (PC) 1/3 but Not PC2

Insulin is first produced in pancreatic β-cells as the precursor prohormone proinsulin. Defective proinsulin processing has been implicated in the pathogenesis of both type 1 and type 2 diabetes. Though there is substantial evidence that mouse β-cells process proinsulin using prohormone convertase 1/3 (PC1/3) and then prohormone convertase 2 (PC2), this finding has not been verified in human β-cells. Immunofluorescence with validated antibodies revealed that there was no detectable PC2 immunoreactivity in human β-cells and little PCSK2 mRNA by in situ hybridization. Similarly, rat β-cells were not immunoreactive for PC2. In all histological experiments, PC2 immunoreactivity in neighboring α-cells acted as a positive control. In donors with type 2 diabetes, β-cells had elevated PC2 immunoreactivity, suggesting that aberrant PC2 expression may contribute to impaired proinsulin processing in β-cells of patients with diabetes. To support histological findings using a biochemical approach, human islets were used for pulse-chase experiments. Despite inhibition of PC2 function by temperature blockade, brefeldin A, chloroquine, and multiple inhibitors that blocked production of mature glucagon from proglucagon, β-cells retained the ability to produce mature insulin. Conversely, suppression of PC1/3 blocked processing of proinsulin but not proglucagon. By demonstrating that healthy human β-cells process proinsulin by PC1/3 but not PC2, we suggest that there is a need to revise the long-standing theory of proinsulin processing.

Structural And Computational Perspectives of Selectively Targeting Mutant Proteins

Diseases are often caused by mutant proteins. Many drugs have limited effectiveness and/or toxic side effects because of a failure to selectively target the disease-causing mutant variant, rather than the functional wild type protein. Otherwise, the drugs may even target different proteins with similar structural features. Designing drugs that successfully target mutant proteins selectively represents a major challenge. Decades of cancer research have led to an abundance of potential therapeutic targets, often touted to be "master regulators". For many of these proteins, there are no FDA-approved drugs available; for others, off-target effects result in dose-limiting toxicity. Cancer-related proteins are an excellent medium to carry the story of mutant-specific targeting, as the disease is both initiated and sustained by mutant proteins; furthermore, current chemotherapies generally fail at adequate selective distinction. This review discusses some of the challenges associated with selective targeting from a structural biology perspective, as well as some of the developments in algorithm approach and computational workflow that can be applied to address those issues. One of the most widely researched proteins in cancer biology is p53, a tumor suppressor. Here, p53 is discussed as a specific example of a challenging target, with contemporary drugs and methodologies used as examples of burgeoning successes. The oncogene KRAS, which has been described as "undruggable", is another extensively investigated protein in cancer biology. This review also examines KRAS to exemplify progress made towards selective targeting of diseasecausing mutant proteins. Finally, possible future directions relevant to the topic are discussed.

Molecular Inhibitor of QSOX1 Suppresses Tumor Growth In Vivo

Quiescin sulfhydryl oxidase 1 (QSOX1) is an enzyme overexpressed by many different tumor types. QSOX1 catalyzes the formation of disulfide bonds in proteins. Because short hairpin knockdowns (KD) of QSOX1 have been shown to suppress tumor growth and invasion in vitro and in vivo, we hypothesized that chemical compounds inhibiting QSOX1 enzymatic activity would also suppress tumor growth, invasion, and metastasis. High throughput screening using a QSOX1-based enzymatic assay revealed multiple potential QSOX1 inhibitors. One of the inhibitors, known as "SBI-183," suppresses tumor cell growth in a Matrigel-based spheroid assay and inhibits invasion in a modified Boyden chamber, but does not affect viability of nonmalignant cells. Oral administration of SBI-183 inhibits tumor growth in 2 independent human xenograft mouse models of renal cell carcinoma. We conclude that SBI-183 warrants further exploration as a useful tool for understanding QSOX1 biology and as a potential novel anticancer agent in tumors that overexpress QSOX1.

Pax4 synergistically acts with Pdx1, Ngn3 and MafA to induce HuMSCs to differentiate into functional pancreatic β-cells

It has been indicated that the combination of pancreatic and duodenal homeobox 1 (Pdx1), MAF bZIP transcription factor A (MafA) and neurogenin 3 (Ngn3) was able to reprogram various cell types towards pancreatic β-like cells (pβLCs). Paired box 4 (Pax4), a transcription factor, has a key role in regulating the maturation of pancreatic β-cells (pβCs). In the present study, it was investigated whether Pax4 is able to synergistically act with Pdx1, Ngn3 and MafA to induce human umbilical cord mesenchymal stem cells (HuMSCs) to differentiate into functional pβCs in vitro. HuMSCs were isolated, cultured and separately transfected with adenovirus (Ad) expressing enhanced green fluorescence protein, Pax4 (Ad-Pax4), Pdx1+MafA+Ngn3 (Ad-3F) or Ad-Pxa4 + Ad-3F. The expression of C-peptide, insulin and glucagon was detected by immunofluorescence. The transcription of a panel of genes was determined by reverse transcription-quantitative PCR, including glucagon (GCG), insulin (INS), NK6 homeobox 1 (NKX6-1), solute carrier family 2 member 2 (SLC2A2), glucokinase (GCK), proprotein convertase subtilisin/kexin type 1 (PCSK1), neuronal differentiation 1 (NEUROD1), ISL LIM homeobox 1 (ISL 1), Pax6 and PCSK type 2 (PCSK2). Insulin secretion stimulated by glucose was determined using ELISA. The results suggested that, compared with Ad-3F alone, cells co-transfected with Ad-Pax4 and Ad-3F expressed higher levels of INS and C-peptide, as well as genes expressed in pancreatic β precursor cells, and secreted more insulin in response to high glucose. Furthermore, the expression of GCG in cells transfected with Ad-3F was depressed by Ad-Pax4. The present study demonstrated that Pax4 was able to synergistically act with the transcription factors Pdx1, Ngn3 and MafA to convert HuMSCs to functional pβLCs. HuMSCs may be potential seed cells for generating functional pβLCs in the therapy of diabetes.

Functional analysis of PCSK2 coding variants: A founder effect in the Old Order Amish population

AIMS

In humans, noncoding variants of PCSK2, the gene encoding prohormone convertase 2 (PC2), have been previously associated with risk for and age of onset of type 2 diabetes (T2D). The aims of this study were to identify coding variants in PCSK2; to determine their possible association with glucose handling; and to determine functional outcomes for coding variants in biochemical studies.

METHODS

Exome-wide genotyping was performed on 1725 Old Order Amish (OOA) subjects. PCSK2 coding variants were tested for association with diabetes-related phenotypes. In vitro analyses using transfected human PC2-encoding constructs were performed to determine the impact of each mutation on PC2 activity.

RESULTS

We identified 10 rare missense coding variants in PCSK2 in various genomic databases. R430W (rs200711626) is greatly enriched in the OOA population (MAF 4.3%). This variant is almost twice as common (MAF 7.4%) in OOA individuals with T2D as in OOA individuals with normal or with normal/impaired glucose tolerance (MAF 3.9% and 2.9%, respectively; p=0.25 and p=0.10). In vitro experiments revealed a broadening of the pH optimum for the R430W variant, which may result in increased activity against PCSK2 substrates.

CONCLUSIONS

Although the association of the R430W variation with T2D in the OOA population did not reach significance, based upon the broadened pH profile of R430W PC2, we speculate that the presence of this substitution may result in altered processing of PCSK2 substrates, ultimately leading to increased conversion to diabetes.

Selective Targeting of Extracellular Insulin-Degrading Enzyme by Quasi-Irreversible Thiol-Modifying Inhibitors

Many therapeutically important enzymes are present in multiple cellular compartments, where they can carry out markedly different functions; thus, there is a need for pharmacological strategies to selectively manipulate distinct pools of target enzymes. Insulin-degrading enzyme (IDE) is a thiol-sensitive zinc-metallopeptidase that hydrolyzes diverse peptide substrates in both the cytosol and the extracellular space, but current genetic and pharmacological approaches are incapable of selectively inhibiting the protease in specific subcellular compartments. Here, we describe the discovery, characterization, and kinetics-based optimization of potent benzoisothiazolone-based inhibitors that, by virtue of a unique quasi-irreversible mode of inhibition, exclusively inhibit extracellular IDE. The mechanism of inhibition involves nucleophilic attack by a specific active-site thiol of the enzyme on the inhibitors, which bear an isothiazolone ring that undergoes irreversible ring opening with the formation of a disulfide bond. Notably, binding of the inhibitors is reversible under reducing conditions, thus restricting inhibition to IDE present in the extracellular space. The identified inhibitors are highly potent (IC50(app) = 63 nM), nontoxic at concentrations up to 100 μM, and appear to preferentially target a specific cysteine residue within IDE. These novel inhibitors represent powerful new tools for clarifying the physiological and pathophysiological roles of this poorly understood protease, and their unusual mechanism of action should be applicable to other therapeutic targets.

Liquid chromatography-electrospray linear ion trap mass spectrometry analysis of targeted neuropeptides in Tac1(-/-) mouse spinal cords reveals significant lower concentration of opioid peptides

Tachykinin and opioid peptides play a central role in pain transmission, modulation and inhibition. The treatment of pain is very important in medicine and many studies using NK1 receptor antagonists failed to show significant analgesic effects in humans. Recent investigations suggest that both pronociceptive tachykinins and the analgesic opioid systems are important for normal pain sensation. The analysis of opioid peptides in Tac1(-/-) spinal cord tissues offers a great opportunity to verify the influence of the tachykinin system on specific opioid peptides. The objectives of this study were to develop an HPLC-MS/MRM assay to quantify targeted peptides in spinal cord tissues. Secondly, we wanted to verify if the Tac1(-/-) mouse endogenous opioid system is hampered and therefore affects significantly the pain modulatory pathways. Targeted neuropeptides were analyzed by high performance liquid chromatography linear ion trap mass spectrometry. Our results reveal that EM-2, Leu-Enk and Dyn A were down-regulated in Tac1(-/-) spinal cord tissues. Interestingly, Dyn A was almost 3 fold down-regulated (p<0.0001). No significant concentration differences were observed in mouse Tac1(-/-) spinal cords for Met-Enk and CGRP. The analysis of Tac1(-/-) mouse spinal cords revealed noteworthy decreases of EM-2, Leu-Enk and Dyn A concentrations which strongly suggest a significant impact on the endogenous pain-relieving mechanisms. These observations may have insightful impact on future analgesic drug developments and therapeutic strategies.

Cationic Cell-Penetrating Peptides Are Potent Furin Inhibitors

Cationic cell-penetrating peptides have been widely used to enhance the intracellular delivery of various types of cargoes, such as drugs and proteins. These reagents are chemically similar to the multi-basic peptides that are known to be potent proprotein convertase inhibitors. Here, we report that both HIV-1 TAT_47-57 peptide and the Chariot reagent are micromolar inhibitors of furin activity in vitro. In agreement, HIV-1 TAT_47-57 reduced HT1080 cell migration, thought to be mediated by proprotein convertases, by 25%. In addition, cyclic polyarginine peptides containing hydrophobic moieties which have been previously used as transfection reagents also exhibited potent furin inhibition in vitro and also inhibited intracellular convertases. Our finding that cationic cell-penetrating peptides exert potent effects on cellular convertase activity should be taken into account when biological effects are assessed.

Synthesis and evaluation of derivatives of the proteasome deubiquitinase inhibitor b-AP15

The ubiquitin-proteasome system (UPS) is increasingly recognized as a therapeutic target for the development of anticancer therapies. The success of the 20S proteasome core particle (20S CP) inhibitor bortezomib in the clinical management of multiple myeloma has raised the possibility of identifying other UPS components for therapeutic intervention. We previously identified the small molecule b-AP15 as an inhibitor of 19S proteasome deubiquitinase (DUB) activity. Building upon our previous data, we performed a structure-activity relationship (SAR) study on b-AP15 and identified VLX1570 as an analog with promising properties, including enhanced potency and improved solubility in aqueous solution. In silico modeling was consistent with interaction of VLX1570 with key cysteine residues located at the active sites of the proteasome DUBs USP14 and UCHL5. VLX1570 was found to inhibit proteasome deubiquitinase activity in vitro in a manner consistent with competitive inhibition. Furthermore, using active-site-directed probes, VLX1570 also inhibited proteasome DUB activity in exposed cells. Importantly, VLX1570 did not show inhibitory activity on a panel of recombinant non-proteasome DUBs, on recombinant kinases, or on caspase-3 activity, suggesting that VLX1570 is not an overtly reactive general enzyme inhibitor. Taken together, our data shows the chemical and biological properties of VLX1570 as an optimized proteasome DUB inhibitor.

Identification of potent and compartment-selective small molecule furin inhibitors using cell-based assays

The proprotein convertase furin is implicated in a variety of pathogenic processes such as bacterial toxin activation, viral propagation, and cancer. Several groups have identified non-peptide compounds with high inhibitory potency against furin in vitro, although their efficacy in various cell-based assays is largely unknown. In this study we show that certain guanidinylated 2,5-dideoxystreptamine derivatives exhibit interesting ex vivo properties. Compound 1b (1,1'-(4-((2,4-diguanidino-5-(4-guanidinophenoxy)cyclohexyl)oxy)-1,3-phenylene)diguanidine) is a potent and cell-permeable inhibitor of cellular furin, since it was able to retard tumor cell migration, block release of a Golgi reporter, and protect cells against Bacillus anthracis (anthrax) and Pseudomonas aeruginosa intoxication, with no evident cell toxicity. Other compounds based on the 2,5-dideoxystreptamine scaffold, such as compound 1g (1,1'-(4,6-bis(4-guanidinophenoxy)cyclohexane-1,3-diyl)diguanidine) also efficiently protected cells against anthrax, but displayed only moderate protection against Pseudomonas exotoxin A and did not inhibit cell migration, suggesting poor cell permeability. Certain bis-guanidinophenyl ether derivatives such as 2f (1,3-bis(2,4-diguanidinophenoxy) benzene) exhibited micromolar potency against furin in vitro, low cell toxicity, and highly efficient protection against anthrax toxin; this compound only slightly inhibited intracellular furin. Thus, compounds 1g and 2f both represent potent furin inhibitors at the cell surface with low intracellular inhibitory action, and these particular compounds might therefore be of preferred therapeutic interest in the treatment of certain bacterial and viral infections.

Nimbolide targets BCL2 and induces apoptosis in preclinical models of Waldenströms macroglobulinemia

Neem leaf extract (NLE) has medicinal properties, which have been attributed to its limonoid content. We identified the NLE tetranorterpenoid, nimbolide, as being the key limonoid responsible for the cytotoxicity of NLE in various preclinical models of human B-lymphocyte cancer. Of the models tested, Waldenströms macroglobulinemia (WM) cells were most sensitive to nimbolide, undergoing significant mitochondrial mediated apoptosis. Notably, nimbolide toxicity was also observed in drug-resistant (bortezomib or ibrutinib) WM cells. To identify putative targets of nimbolide, relevant in WM, we used chemoinformatics-based approaches comprised of virtual in silico screening, molecular modeling and target–ligand reverse docking. In silico analysis revealed the antiapoptotic protein BCL2 was the preferential binding partner of nimbolide. The significance of this finding was further tested in vitro in RS4;11 (BCL2-dependent) tumor cells, in which nimbolide induced significantly more apoptosis compared with BCL2 mutated (Jurkat BCL2^Ser70-Ala) cells. Lastly, intraperitoneal administration of nimbolide in WM tumor xenografted mice, significantly reduced tumor growth and IgM secretion in vivo, while modulating the expression of several proteins as seen on immunohistochemistry. Overall, our data demonstrate that nimbolide is highly active in WM cells, as well as other B-cell cancers, and engages BCL2 to exert its cytotoxic activity.

Identification of a small molecule that selectively inhibits mouse PC2 over mouse PC1/3: a computational and experimental study

The calcium-dependent serine endoproteases prohormone convertase 1/3 (PC1/3) and prohormone convertase 2 (PC2) play important roles in the homeostatic regulation of blood glucose levels, hence implicated in diabetes mellitus. Specifically, the absence of PC2 has been associated with chronic hypoglycemia. Since there is a reasonably good conservation of the catalytic domain between species translation of inhibitory effects is likely. In fact, similar results have been found using both mouse and human recombinant enzymes. Here, we employed computational structure-based approaches to screen 14,400 compounds from the Maybridge small molecule library towards mouse PC2. Our most remarkable finding was the identification of a potent and selective PC2 inhibitor. Kinetic data showed the compound to be an allosteric inhibitor. The compound identified is one of the few reported selective, small-molecule inhibitors of PC2. In addition, this new PC2 inhibitor is structurally different and of smaller size than those reported previously. This is advantageous for future studies where structural analogues can be built upon.

The biology and therapeutic targeting of the proprotein convertases

The mammalian proprotein convertases constitute a family of nine secretory serine proteases that are related to bacterial subtilisin and yeast kexin. Seven of these (proprotein convertase 1 (PC1), PC2, furin, PC4, PC5, paired basic amino acid cleaving enzyme 4 (PACE4) and PC7) activate cellular and pathogenic precursor proteins by cleavage at single or paired basic residues, whereas subtilisin kexin isozyme 1 (SKI-1) and proprotein convertase subtilisin kexin 9 (PCSK9) regulate cholesterol and/or lipid homeostasis via cleavage at non-basic residues or through induced degradation of receptors. Proprotein convertases are now considered to be attractive targets for the development of powerful novel therapeutics. In this Review, we summarize the physiological functions and pathological implications of the proprotein convertases, and discuss proposed strategies to control some of their activities, including their therapeutic application and validation in selected disease states.

Furin inhibitors: importance of the positive formal charge and beyond

Furin is the prototype member of the proprotein convertases superfamily. Proprotein convertases are associated with hormonal response, neural degeneration, viral and bacterial activation, and cancer. Several studies over the last decade have examined small molecules, natural products, peptides and peptide derivatives as furin inhibitors. Currently, subnanomolar inhibition of furin is possible. Herein, we report the analysis of 115 furin inhibitors reported in the literature. Analysis of the physicochemical properties of these compounds highlights the dependence of the inhibitory potency with the total formal charge and also shows how the most potent (peptide-based) furin inhibitors have physicochemical properties similar to drugs. In addition, we report docking studies of 26 furin inhibitors using Glide XP. Inspection of binding interactions shows that the two putative binding modes derived from our study are reasonable. Analysis of the binding modes and protein-ligand interaction fingerprints, used here as postdocking procedure, shows that electrostatic interactions predominate on S1, S2 and S4 subsites but are seldom in S3. Our models also show that the benzimidamide group, present in the most active inhibitors, can be accommodated in the S1 subsite. These results are valuable for the design of new furin inhibitors.