Next Article in Journal
Taurine and Polyphenol Complex Repaired Epidermal Keratinocyte Wounds by Regulating IL8 and TIMP2 Expression
Next Article in Special Issue
Marine Bioactive Molecules as Inhibitors of the Janus Kinases: A Comparative Molecular Docking and Molecular Dynamics Simulation Approach
Previous Article in Journal
TAAR8 Mediates Increased Migrasome Formation by Cadaverine in RPE Cells
Previous Article in Special Issue
Docking, MD Simulations, and DFT Calculations: Assessing W254’s Function and Sartan Binding in Furin
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
CIMB
Volume 46
Issue 8
10.3390/cimb46080511
cimb-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Potential Furin Protease Inhibitory Properties of Pioglitazone, Rosiglitazone, and Pirfenidone: An In Silico Docking and Molecular Dynamics Simulation Approach

by
Ahtziri Socorro Carranza-Aranda
1,
Carlos Daniel Diaz-Palomera
2,
Eduardo Lepe-Reynoso
2,
Anne Santerre
3,
José Francisco Muñoz-Valle
2 and
Oliver Viera-Segura
2,*
1
Doctorado en Ciencias Biomédicas, Centro Universitario de Ciencias de la Salud, Universidad de Guadalajara, Guadalajara 44340, Jalisco, Mexico
2
Instituto de Investigación en Ciencias Biomédicas, Centro Universitario de Ciencias de la Salud, Universidad de Guadalajara, Guadalajara 44340, Jalisco, Mexico
3
Departamento de Biología Celular y Molecular, Centro Universitario de Ciencias Biológicas y Agropecuarias, Universidad de Guadalajara, Zapopan 45221, Jalisco, Mexico
*
Author to whom correspondence should be addressed.
Curr. Issues Mol. Biol. 2024, 46(8), 8665-8684; https://doi.org/10.3390/cimb46080511
Submission received: 26 April 2024 / Revised: 23 July 2024 / Accepted: 29 July 2024 / Published: 8 August 2024
(This article belongs to the Special Issue Synthesis and Theoretical Study of Bioactive Molecules)
Browse Figures
Review Reports Versions Notes

Abstract

:
Furin (Fur) is a member of the protease convertase family; its expression is crucial for cleaving and maturing many proteins. Fur also represents a therapeutic target in cancer, autoimmune diseases, and viral infections. Pioglitazone (PGZ) and rosiglitazone (RGZ) are thiazolidinediones prescribed to type 2 diabetes patients and are structurally similar to the known Fur inhibitors naphthofluorescein (NPF) and pirfenidone (PFD). Thus, this study used molecular docking and molecular dynamics to assess and compare the affinities and the molecular interactions of these four ligands with the Fur active site (FurAct) and the recently described Fur allosteric site (FurAll). The 7QXZ Fur structure was used for molecular dockings, and for the best pose complexes, molecular dynamics were run for 100 ns. The best affinities of the ligand/FurAct and ligand/FurAll complexes were with NPF, PGZ, and RGZ, while PFD presented the lowest affinity. Asp154 was the central residue involved in FurAct complex formation, while Glu488 and Asn310 were the central residues involved in FurAll complex formation. This study shows the potential of RGZ, PGZ, and PFD as Fur competitive (FurAct) and non-competitive (FurAll) inhibitors. Therefore, they are candidates for repurposing in response to future emerging diseases through the modulation of Fur activity.

1. Introduction

Secretion protein processing is crucial in developing and functioning higher organism proteins [1,2]. Many of these proteins are synthesized as inactive precursors; therefore, proteolytic cleavage is a critical step for bioactive protein production [3]. Proprotein/prohormone convertase (PC) is one of the most important groups of endoproteolytic processing proteins in mammals and consists of enzymes that exhibit a catalytic subtilisin-like serine protease domain [4,5]. The PC family comprises nine members: PC1/3, PC2, PC4, PC5/6, PACE4, PC7, SKI-1/S1P, PCSK9, and furin. Only four enzymes (PC5/6, PACE4, PC7, and furin) are ubiquitously expressed in tissues and are responsible for the majority of protein processing events [6].
One of the most widely studied proteases is furin (Fur), also known as a dibasic-processing enzyme, paired basic amino acid residue-cleaving enzyme (PACE), and protein convertase subtilisin/kexin 3 (PCSK3). This protease comprises 794 amino acids divided into catalytic and control domains, the latter of which is known as the P domain, whose primary function is stabilizing the catalytic domain [5,6,7]. Due to its ubiquitous expression and localization in the trans-Golgi network, cell surface, and endosomes, Fur plays crucial roles in cleaving and maturing many receptor proteins, enzymes, ligands, and growth factors. Specifically, it is involved in processing Notch, Wnt, and fibroblast growth factor receptors (FGFRs), as well as several cytokines and chemokines. It is also important in regulating cell proliferation, differentiation, and migration.
Fur is also implicated in pathological processes such as cancer, autoimmune diseases, and viral infection, affecting mRNA and protein expression in several tissues [8]. The role of Fur in cancer has been extensively studied, showing that it promotes tumorigenesis and increases malignant tumor phenotype in human head and neck cancer cells [9], as well as in breast [10,11], colorectal [12], and ovarian [13] cancers. Fur also plays a role in the immune tolerance of cancer by contributing to the production of extracellular matrix components to the differentiation of cancer-associated fibroblasts, the epithelial–mesenchymal transition, and the Treg phenotype through TGF-β [14,15,16]. In this sense, Fur inhibition leads to decreased tumorigenicity, motility, and invasiveness of cancer cells [17].
Concerning autoimmune diseases, Fur activity is found to be increased in inflammatory conditions; therefore, it is strongly related to the progression of atherosclerosis [18], rheumatoid arthritis [19,20], lupus erythematosus [21], and Sjogren’s syndrome [22].
In addition, Fur is involved in bacterial and viral infections, increasing pathogenicity. During bacterial infections, this protein processes and modifies toxins, such as anthrax, diphtheria, Shiga toxin, pseudomonas exotoxin A, and Clostridium septicum alpha-toxin, enhancing their toxicity [23,24,25]. Fur also plays a critical role in regulating the viral life cycle during viral infection by processing and maturing glycoproteins on the surface of viral particles [8,26,27]. For example, the human papillomavirus (HPV) requires Fur-mediated cleavage of the L2 minor protein for de novo infections. However, this cleavage can also occur independently from Fur during viral assembly within the host cell [27].
Similarly, flaviviruses, including dengue, Zika, Japanese encephalitis, and West Nile viruses, also rely on Fur-mediated cleavage of the pre-membrane protein for maturation and transitioning from inert to infectious viral particles [8,27,28,29]. In the case of the hepatitis B virus (HBV), Fur serves a dual purpose, aiding in the viral life cycle and processing the E-antigen essential for immune tolerance, potentially leading to chronic infections [30,31]. Due to the above reasons, Fur has become a key target for treating multiple pathologies.
In terms of viruses affecting the respiratory system, influenza virus proteins, specifically those containing hemagglutinin types 5, 7, and 9, are susceptible to Fur cleavage, which significantly increases their pathogenicity and infectivity [8,27,32]. In particular, most coronaviruses use Fur to cleave the spike protein (S), one of the most promising targets for therapy against MERS-CoV, SARS-CoV, and SARS-CoV-2. The spike protein cleavage facilitates interaction with ACE2 for entry into the host cell [33,34,35,36]. Thus, recent efforts to find therapeutic drugs for SARS-CoV-2 have led to the evaluation and repurposing of potential Fur inhibitory molecules using in silico and in vitro assays. In this regard, naphthofluorescein (NPF) has been recently reported as a Fur inhibitor useful against SARS-CoV-2 replication [37], while pirfenidone (PFD), an FDA-approved molecule for treating idiopathic pulmonary fibrosis [38], has recently been recommended as a potential treatment for COVID-19, ameliorating the cytokine storm induced by SARS-CoV-2 infection [39]. Previous studies have shown that PFD also induces a modest Fur inhibition (~50%) at relatively high concentrations (5–8 mM) [40]. However, the mechanism of Fur inhibition by NPF or PFD is unknown.
The current emphasis in pharmacological repurposing is focused on molecules identified as effective inhibitors of Fur activity containing imidazole, such as the sulconazole [10], and thiazolidine rings [5,41], such as thiazolidinediones (TZDs), a group of molecules to which pioglitazone (PGZ) and rosiglitazone (RGZ) belong, which are primarily used in treating type 2 diabetes [42], acting as an α-lipoic acid and activating peroxisome proliferator-activated receptor γ (PPARγ), a nuclear transcription factor that controls glucose metabolism and lipid homeostasis [43]. Beyond this well-established role, TZDs exhibit antifibrotic, anti-inflammatory, and antitumoral effects in various tissues [44]. These effects encompass a reduction in TGF-β expression, extracellular matrix deposition, NFkβ activation, TNFα levels, and inhibition of p38 MAPK activation. TZDs have recently been reported for their antiviral potential against SARS-CoV-2 [45,46,47,48,49]. Although these ligands show great promise, their precise mechanisms of action remain incompletely understood. Considering the promising Fur inhibitor activity of the thiazolidine rings of the TZDs, PGZ, and RGZ and well-known Fur inhibitors such as NPF and PFD, this study focused on the in silico analysis of their molecular interactions and affinity with the active and allosteric sites of the Fur protein and assessed their potential for repositioning as Fur inhibitors for future medical interventions.

2. Materials and Methods

2.1. Molecular Structures

The three-dimensional (3D) structures of the Fur protein were obtained from the Protein Data Bank (PDB) under IDs 7QXZ, 7LCU, 5MIM, and 7O1Y [50,51,52,53]. The potential Fur ligands included four heterocyclic compounds: NPF (a heteropolycyclic molecule), PFD, PGZ, and RGZ (heteromonocyclic molecules). NPF was obtained from the PubChem CID database, while PFD, PGZ, and RGZ structures were retrieved from the ZINC15 database [54] (Table 1). Additionally, the 3-(3,5-dichlorophenyl)-pyridine derivative (3DPP-D) molecule was used as a reference/control during molecular re-dockings. 3DPP-D is a Fur inhibitor and has been previously reported to be co-crystallized with the Fur structure (Figure S1) [50].

2.2. Preparation of Protein and Ligand

In a preliminary step, different Fur structures (PDB IDs: 7QZX, 7LCU, 5MIM, and 7O1Y) were selected to determine the optimum structures for the dockings analysis. All structures were subjected to molecular adaptations, such as the removal of water molecules and extra ligands, addition of hydrogen atoms, and determination of AM1-BCC atomic charges with bond charge corrections (BCCs); energy minimization was also performed, and the structures were refined [55]. Likewise, the structures of the four ligands were prepared by adding hydrogen atoms (polar and non-polar) and Gasteiger atomic charges [56]. All these adjustments were performed using the UCSF Chimera v.1.1.6 (University of California, Oakland, CA, USA), while AutoDock Tools v.1.5.6 (Center for Computational Structural Biology, La Jolla, CA, USA) was used for the creation of the pdbqt files [57,58,59,60].

2.3. Molecular Docking Analysis

Using redocking, we identified and validated the coordinates of the active site of Fur (FurAct) based on the binding site location of the following co-crystallographic ligands previously reported as highly potent Fur inhibitors: 3DPP-D (PDB ID: 7QXZ), BOS-318 (PDB ID: 7LCU), 2,5-dideoxystreptamine-derived (PDB ID: 5MIM), and guanylhydrazone-based inhibitor 2 (mi307) (PDB ID: 7O1Y) [50,51,52,53]. This strategy allowed us to determine that the 7QXZ Fur structure was optimum and select it for subsequent molecular dockings (Table S1 and Figure S1).
Considering the resulting parameters, 100 run modes of localized molecular docking on FurAct were carried out for the four molecules of interest (NPF, PFD, PGZ, and RGZ) with a box size of 26 Å × 26 Å × 26 Å and a center x = 51.42, y = −32.10, and z = −6.66 using AutoDock Vina v.2.5.4 and AutoDock4 v.4.2.6 (Center for Computational Structural Biology, La Jolla, CA, USA) and a box size of 40 Å × 40 Å × 40 Å and a center x = 49.35, y = 33.95, and z = −4.297 using the Lamarckian genetic algorithm [61,62].
Furthermore, to identify possible interactions of the four evaluated molecules with important allosteric inhibition sites of Fur, we also ran 100 modes of blind molecular docking, i.e., considering the entire Fur structure, and identified a second interaction site that we named FurAll. This site corresponds to a recently described allosteric site of Fur, as discussed below. The calculations were carried out with the Blind Docking Server [63] and confirmed by AutoDock Vina v.2.5.4 (box size of 65 Å × 58 Å × 70.89 Å with a center x = 35.50, y = −38.54, and z = −1.63) and AutoDock4 v.4.2.6 (box size of 126 Å × 126 Å × 126 Å with a center x = 35.73, y = −38.043, and z = 0.156). The complexes with the highest affinity scores during redocking and docking assays were visualized and analyzed for covalent and non-covalent bonds using PyMOL2 v.2.5.4 [64] and Discovery Studio visualizer v.21.1.0.20298 (BIOVIA, San Diego, CA, USA). Contacts between molecules were considered “real” when bonds presented a length lower or equal to 3.3 Å [65]. The results of the blind docking were corroborated using AutoDock Vina v.2.5.4 and AutoDock4 v.4.2.6 (Figure S3, Table S3) [65]).

2.4. Molecular Dynamics Simulation

The best ligand–Fur complexes obtained with each of the four molecules under study and FurAct and FurAll were then prepared for molecular dynamics (MD) simulations using the Nanoscale Molecular Dynamics (NAMD) v.2.14 software (University of Illinois at Urbana–Champaign, Champaign, IL, USA) [66]. The force field parameters for Fur were generated with the CHARMM 27 package, and for the ligands, the topology files were obtained from the CHARMM-GUI server [67,68].
The complete process was performed in five steps according to the NAMD protocols found in https://www.ks.uiuc.edu/Training/Tutorials/ (accessed on 15 September 2023). First, all complexes were prepared by minimizing the structures of the proteins and ligands, eliminating water molecules in 2000 steps. Second, all complexes were solvated in a 3D TIP3P water cube (5 Å dimension), and the periodic boundary conditions presented the following dimensions: x = 75.33 Å, y = 96.17 Å, and z = 87.57 Å for the active site (FurAct) and x = 69.65 Å, y = 65.61 Å, and z = 75.20 Å for the allosteric site (FurAll). Third, all systems were neutralized and minimized by adding sodium, calcium, and chlorine ions. The temperature was set at 310 K and the pressure at 1 atm for two fs for 100 ps. Fourth, all systems were subjected to an equilibrium phase for 200 ps to subsequently run the production dynamics for 100 ns. During the 100 ns of MD generation, the pressure was held at 1 atm by employing the Nose–Hoover–Langevin piston barostat 1.0132 with a Langevin piston decay of 0.1 ps and a period of 0.2 ps. Finally, the MD simulation was analyzed with Carma v.2.01 (Alexandroupolis, Greece) and VMD v.1.9.4 (University of Illinois at Urbana–Champaign, IL, USA) software to examine the root mean square deviation (RMSD), root mean square fluctuation (RMSF), and radius of gyration of amino acids (Rg) through the fluctuations of the furin c-alpha [69,70]. RMSD, Rg, and RMSF are indicators of the element in the MD of protein/ligand complexes and provide us with valuable information to compare the stability of the two protein structures, measure the stability of electrons from their center of gravity, and determine the flexibility of each protein residue during a simulation period. Lower RMSD, Rg, and RMSF values indicate more stability and feasibility of the complex formation [71,72,73]. The graphs were elaborated with RStudio v.4.2.1, and the complexes obtained from MD were compared with those analyzed by molecular docking to corroborate the stability of the interactions; this was conducted using PyMOL2 v.2.5.4 (Schrödinger, New York, NY, USA) [63] and Discovery Studio visualizer v.21.1.0.20298 (BIOVIA, San Diego, CA, USA. Figure 1 shows the workflow diagram for the studied molecules.

3. Results

3.1. Molecular Docking

The Fur structure presents two important regulatory regions: the catalytic domain and the P domain. Within the catalytic domain, the active site (Asp154, His194, Asn295, and Ser368) plays a critical role in the peptidase activity; meanwhile, the P domain (Gly265, Ans310, Gln488, and Ala532) stabilizes the catalytic domain, modulating its activity (Figure 2).
In this study, we assessed the molecular interactions of NPF, PFD, PGZ, and RGZ at both sites and determined their binding affinities as well as the inhibitory potential of Fur’s enzymatic activity. The localized molecular docking assays showed that NPF presented the highest affinity to the Fur active site (FurAct) with a score of −9.9 kcal/mol, followed by PGZ (−8.6 kcal/mol), RGZ (−8.1 kcal/mol), and PFD (−6.9 kcal/mol). Asp154 was the central amino acid residue involved in the ligand/protein complexes with NPF, PGZ, and RGZ. The variation in the affinity score of these three molecules was influenced by the number of polar contacts (1, 3, and 1, respectively) and the distance of the polar bond between Asp154 and its ligands. This distance was longer for NPF (2.7 Å) compared to PGZ (2.6 Å) and RGZ (2.2 Å). On the other hand, the complex formed between FurAct and PFD (−6.9 kcal/mol) presented low affinity and was orchestrated by a single polar bond with Gly255 (Table 2 and Figure 3A–F).
Remarkably, the blind docking with FurAll showed that three of the four ligands presented higher affinity forces for the Fur allosteric than for the Fur active site. NPF showed the highest affinity score (−10.9 kcal/mol), followed by PGZ (−8.7 kcal/mol), RGZ (−7.9 kcal/mol), and PFD (−7.1 kcal/mol). Glu488 and Asn310 were the main amino acid residues involved in the complexes’ formation and stabilization (Table 3 and Figure 4A–F). In all cases, non-covalent interactions, such as electrostatic forces (van der Waals, Pi–cation, and Pi–anion), hydrophobic (Pi and alkyl hydrophobic), and sulfur bonds, favored the stability of the ligand/protein complexes (Table 3, Figure 3A–F and Figure 4A–F).
As for the redocking of the co-crystallized ligand derived from 3DPP-D, it only showed high affinity towards FurAct, where the molecular interactions involved Asp233 (conventional H-bond), Pro256, Asp264, Glu236, and Asp154 (electrostatic interactions) and Trp254, Ala252, Met226, Leu240, Trp291, Val231, and Gly255 (hydrophobic interactions) (Table S2 and Figure S2).
Thus, the molecular docking assays indicated that according to the number of H-bonds, NPF, RGZ, and PFD presented a higher affinity for FurAll than FurAct. In comparison, PGZ showed a higher affinity for the active site, although the binding affinity was more significant with FurAll. The non-polar bindings, such as hydrophobic interactions, probably promoted the increase in affinity to FurAct.
Based on the molecular docking assay using AutoDock4, we obtained the inhibition constant (Ki) of Fur for each of the study molecules. We observed that the theoretical inhibition concentration from the active site was very low for PGZ (1.28 μM), followed by RGZ (3.33 μM) and PFD (24.87 μM); however, for NPF, it could not be calculated. On the other hand, the Ki values towards the allosteric site calculated from blind molecular docking showed a decrease in inhibition concentration for NFP, PGZ, and RGZ. In contrast, for PFD, it could not be calculated due to their preference for FurAct (Table 4).

3.2. Molecular Dynamic

To further validate the stability of the ligand/protein complexes, especially those targeting the active and allosteric sites of Fur, MD simulations were conducted, starting with the most optimal conformations identified through molecular docking. The movements of atoms and molecules of Fur alone (ligand-free) in systems (water and ions) were also considered to compare the conformational changes in the presence of molecules.
The structural integrity of all complexes was assessed using the root mean square deviation of atomic position value (RMSD). The mean of RMSD of the ligand-free Fur protein (3.98 ± 0.07 Å), PFD/FurAct (4.21 ± 0.08 Å), and PGZ/FurAct (4.12 ± 0.06 Å) complexes were similar, unlike the NPF/FurAct (4.67 ± 0.13 Å) and RGZ/FurAct (4.77 ± 0.14 Å) complexes, which displayed higher RMSD values (Table 5 and Figure 5A). Notably, the interactions of NPF and RGZ with the active site of Fur (FurAct) exhibited greater RMSD fluctuations in the last 10% of the MD simulation.
Among the four analyzed ligands/FurAll complexes, NPF presented the lowest RMSD value (4.52 ± 0.05 Å), followed by PGZ (4.66 ± 0.17 Å), PFD (5.22 ± 0.17 Å), and RGZ (5.57 ± 0.17 Å) (Table 5 and Figure 5B). The RGZ/FurAll, PGZ/FurAll, and PFD/FurAll complexes showed higher RMSD values compared to their RGZ/FurAct, PGZ/FurAct, and PFD/FucAct counterparts, suggesting that they established more stable contacts at the Fur active site than at the Fur allosteric site, probably due to a stabilizing failure at the end of the MD simulations. Contrary to the above, NPF/FurAll established a more stable complex, as supported by the statistical similarity between the RMSD values of these complexes; also, the RMSD values in NPF/FurAll and PGZ/FurAll were similar to those in ligand-free Fur protein (Table 5 and Figure 5B).
Based on this information, we assessed the ligand/protein complexes’ general stability and structural compactness for both FurAct and FurAll using the radius of gyration (Rg) values. For FurAct, we observed a high consistency of Rg values in the union of NPF and PGZ throughout the entire trajectory, compared to those with PFD and RGZ. The Rg value of PFD/FurAct remained similar to those of the ligand-free Fur protein during the entire MD simulation.
The degree of compactness of the complexes was similar (lowest Rg values) for FurAll and FurAct, highlighting the NPF/FurAll complex as the most compact one, followed by PFD/FurAll and PGZ/FurAll (Table 6 and Figure 6A,B). Regarding the complex RGZ/FurAll, an increase in complex size was observed in the last 10 ns (Table 6 and Figure 6B).
Another critical parameter that contributes to understanding the dynamics of the ligand/protein complexes during MD simulation is the RMSF value (root mean square fluctuation). This value calculates the flexibility of specific amino acid residues and how they contribute to molecular motion within the complexes. High RMSF values indicate high amino acid mobility, as opposed to low RMSF values, which are associated with secondary protein structures with high stability.
The complexes between the active site of Fur and NPF, PFD, PGZ, and RGZ showed a more significant number of unique fluctuating regions in contrast to the ligand-free Fur. In particular, the NPF/FurAct complex showed seven highly fluctuating regions located in the catalytic domain, with Regions 2 (Arg185, Tyr186, Gln188, Met189, and Asn190) and 5 (Phe323, Gly324, Asn325, Val326, and Pro327) being the most prominent, with distances of approximately 5.3 and 4 Å, respectively (Figure 7A).
The PFD/FurAct and RGZ/FurAct complexes also showed nine regions with high fluctuation, located mainly in the catalytic domain and, to a lesser extent, in the P domain (Figure 7A). Finally, the PGZ/FurAct complex showed the fewest regions of the fluctuating protein as the residue modifications, mainly in the catalytic domain, were smaller than for the other complexes; however, Region 2 (Pro184, Arg185, Tyr186, Thr187, Gln188, Met189, and Asn190) of this complex showed the most significant movement in its structure, reaching a peak of ~5.6 Å (Figure 7A).
In contrast to previous observations of FurAct, the FurAll complexes showed lower RMSF values and fewer fluctuating regions, comparable to that of ligand-free Fur. Similar to what we observed for the NPF/FurAll complex, only four highly fluctuating regions were observed compared to the Fur-alone protein, highlighting a region located in the P domain (Leu475 and Gly476) with an RMSF value of 4.15 Å (Figure 7B).
Also, the NPF/FurAll, PFD/FurAll, and RGZ/FurAll complexes showed five regions with high fluctuation with respect to the ligand-free Fur observed with FurAct. In both cases, the region comprising amino acids Asp258 to Val263 was the most variable, with an RMSF peak of ~5 A. The other regions identified in both complexes showed fluctuations ranging from 4 to 1.7 Å (Figure 7B).
Interestingly, the PGZ/FurAll complex’s fluctuating regions increased to eight when the ligand interacted with the Fur allosteric site. The regions with the greatest movement correspond to amino acids Thr187 to Met189 and Asp258 to Thr262, located near the Fur active site (Figure 7B).
The data obtained from RMSD, Rg, and RMSF suggest that the complexes formed with FurAll are more stable and compact than those directed to the FurAct, which could confirm that FurAll contains an important binding and control site of Fur and is adequate for binding with our ligands. Furthermore, when evaluating the ligands/Fur complexes’ specific interactions towards the active site in the last 10% of the MD simulation trajectories, we observed that all the ligands, except for NPF, remained in the determined site of interaction. Similarly, the MD simulations showed that all ligands respected their interaction site towards FurAll (Figures S4–S7). However, the amino acid residues involved in establishing the complexes at the end of the MD simulation of FurAct the FurAll sites differed from those observed by molecular docking (Figures S4–S7).
Despite the difference in amino acid residues involved in the formation of the complexes, all models, except the ligand-free Fur protein model, showed a similar conformation to that observed by molecular docking (Figure S8).

4. Discussion

Most of the existing Fur inhibitors are based on polypeptide molecules that resemble the native substrates of this enzyme, occupying the space within the catalytic site and inhibiting its function, as reported for peptidyl chloromethyl ketone (CMK), α1-antitrypsin variants, and meta-guanidinomethyl-Phac-RVR-Amba, among others [25,53,74,75]. Besides this, multiple low-molecular-weight non-peptide molecules may have an inhibitory capacity of Fur proteolytic activity, for example, the guanidinylated aryl 2,5-dideoxystreptamine derivatives (GADDs). However, these are not often used as therapeutic molecules due to their potential toxicity and lack of research efforts. Thus, the repositioning of therapeutic molecules (used for treating known conditions) represents an interesting and relatively rapid strategy to establish new therapeutic alternatives [37,40,76,77].
The use of non-peptide molecules as Fur inhibitors presents several advantages compared to peptide inhibitors, such as improved solubility, bioavailability, cell permeability, and metabolic and proteolytic stability. Some small non-peptide molecules that have shown inhibitory effects in Fur are labdane diterpenes, terpyridine/quinoline complexes of Zn+2/Cu+2 ions, dideoxy streptamine analogs, dicumarine derivates, flavonoids, and NPF analogs. Most prior studies of these molecules have been carried out in vitro, providing information about their efficacy or toxicity [17,78,79], but lack information on their interaction with their target protein. For these reasons, the in silico study of the molecular interactions of non-peptide molecules currently approved to treat different diseases is crucial for their rapid repositioning to treat other existing or emerging ailments. Specifically, the non-peptide NPF and PFD molecules have been reported for their inhibitory function of Fur in in vitro models of several ailments, such as cancer, inflammatory and metabolic diseases, idiopathic pulmonary fibrosis, viral and bacterial infections, and neurological diseases such as Alzheimer’s disease [80,81]. However, the mechanisms and molecular interactions conferring their inhibitory property have not been studied so far. The in silico approach used here assessed the molecular interactions of NPF and, for the first time, PFD with active (FurAct) and allosteric (FurAll) sites of Fur and calculated the affinity forces of the ligand/protein complexes. Our results show that both molecules have a high affinity for FurAct and even higher affinity for FurAll.
The FurAct site consists of a highly twisted β-sheet composed of seven parallel and one antiparallel β-strand, flanked by five adjacent and two peripheral α-helices and two β-hairpin loops, which regulate Fur protease activity, specifically through the surrounding negatively charged aspartic acid and glutamic acid residues (Asp258, Asp306, Asp154, Asp19, Glu236, Glu264, Glu257, Glu230, and Glu233) as well as the catalytic amino acid triad (Asp153, His194, and Ser368) [41,82,83,84]. Furthermore, it has been reported that besides the key role of the active site in the activation of Fur, in the P domain, there may exist an allosteric site (interdomain), which plays a crucial role in the stabilization and functions of this protein [2,82]. This interdomain connection is located between the eight-stranded β-sandwich of the P domain (Pβ5 and Pβ6) and the jelly-roll β-barrels of the catalytic domain (loops Cβ5-Calpha4, Cβ6-Cβ7, Cβ7-Cβ8, and Cβ11-Cβ12). It is highly conserved between proprotein/prohormone convertases (PC) and is believed to stabilize the catalytic domain [82]. A hydrophobic core and several polar–salt bridges mediate the property of this region. Also, the observation that these interactions shield it from solvent surface patches in this area [82] may be an ideal characteristic for interactions with other inhibitory molecules [52,74].
In a recent study, Feng et al. (2023) [74], using molecular docking, predicted the existence of a new potential allosteric site and determined its localization between the S8 region of the catalytic domain (aa121-435) and the P domain (aa444-576) of Fur. This localization corresponds with the amino acid residues observed in the present study (Gly265, Ans310, Gln488, and Ala532) during MD with FurAll, predicting the interaction of our four ligands with the potential allosteric site within the interdomain region. In this regard, Feng and Cols (2023) [75] showed that the activity of Fur decreased upon permethrin (tested inhibitor) binding in a dose-dependent way even though Fur´s Km did not, suggesting that this protein presents an alternative control region, which the authors named the allosteric site. The authors proposed this site as a non-competitive binding site that mediates the action of Fur. Thus, the molecular dockings carried out in the present study suggest that NPF and PFD presented higher affinity and preferential binding for FurAll than FurAct, possibly due to their non-competitive binding and the ease of forming polar bonds in the complexes.
In the case of the NPF/FurAll complex, we observed three polar bonds involving Gly265, Ans310, and Gln488 residues of FurAll but only one polar bond involving Asp154 of FurAct. Peiter et al. (2022) [81] reported similar results, where the amino acids Gln488 and Ala532 were involved in stabilizing the NPF/Fur complex, which may be related to the high inhibitory capacity of this drug.
On the other hand, for the PFD/FurAll complex, we observed that the main amino acid involved through polar bonds was Ans310, while for the PFD/FurAct complex, the interactions were mainly due to non-polar forces through Pi–sulfur, Pi–sigma, and amide–Pi stacked bonds and only one polar contact with Gly255, which makes the complex more unstable. Dahms et al. (2016) [85] reported that Gly255 was a key substrate recognition element.
Similar to the observations for NPF and PFD, our results showed a higher affinity for the two thiazolidinediones tested here (PGZ and RGZ) for FurAll than for FurAct structures. The molecular dockings calculated affinity forces of PGZ and RGZ with FurAct of −8.6 and −8.1 kcal/mol, respectively, whereas with FurAll, the affinities were −8.7 and −7.9 kcal/mol, respectively. The decrease in affinity forces may be due to the number and type of interactions, the distance between interactions, and the amino acid residues involved in forming these complexes.
The PGZ/FurAct and PGZ/FurAll complexes presented three and two polar bonds, respectively, which stabilized the ligand/protein interactions. The most significant differences in the calculated affinities in both complexes were the distances of the polar bonds, which were smaller for the FurAct (2.5, 2.5, and 2.6 Å) than for FurAll (2.28 and 2.7 Å). Likewise, the PGZ/FurAct complex presented many hydrophobic forces. Compared to PGZ/FurAll, hydrophobic forces were the main forces involved in these ligand/protein interactions. They may contribute to a lipophilic environment of the binding cavity, which would facilitate ligand/protein interactions [86,87].
PGZ/FurAct and PGZ/FurAll complexes showed high affinity forces. The number of polar bonds and the hydrophobic forces were higher for PGZ/FurAll (Ser311 and Gln488) than for PGZ/FurAct (Asp154). However, the amino acids that participate in H-bond formation and hydrophobic forces of the PGZ/FurAll complex promoted high stability due to their localization of Pβ6, Calpha4, and Cβ7, which could increase the stability of the catalytic site and the activity of Fur [82].
NFP, PFD, PGZ, and RGZ are molecules that contain heterocyclic systems in their structure; these are relevant in pharmacology because they can modulate properties at both the pharmacokinetic and pharmacodynamic levels [88] and add important properties, such as lipophilicity, interaction capacity, high polarity, and toxic effect. Specifically, TZDs have a heterocycle sulfur functional group, which promotes various pharmacological activities [89]. In addition, the presence of oxygen-based functional groups seems to be an important factor for the stability and high affinity in all complexes.
An MD simulation was then run to validate the molecular docking observations. We observed that in the last 10% of the MD trajectory, the average RMSD was lower in the complexes directed towards the FurAct than towards the FurAll site. PGZ/FurAll (4.66 Å) and NPF/FurAll (4.52 Å) complexes were the most stable during the trajectory compared to the complexes towards the catalytic site (FurAct), which were similar to that of the ligand-free Fur trajectory (3.98 Å). Similarly, the Rg and RMSF values of the PGZ/FurAll complex were 22.93 Å, comparable to the ligand-free protein (22.46 Å), with few fluctuating regions during complex formation. However, PFD and RGZ were more stable. They presented more compactness towards FurAct (RMSD: 4.21 ± 0.08 Å, Rg: 2.57 ± 0.05 vs. RMSD: 4.12 ± 0.06 Å, Rg: 22.65 ± 0. 06) than the allosteric site (FurAll), contrary to what was observed with NPF (RMSD: 4.67 ± 0.13 Å, Rg: 22.62 ± 0.09) and RGZ (RMSD: 4.77 Å, Rg: 22.77 ± 0.10) regarding the allosteric site and Fur alone (RMSD: 3.98 ± 0.07 Å, Rg: 22.45 ± 0.06). Thus, PGZ could be a candidate for a non-competitive Fur inhibitor. At the same time, RGZ could be a candidate for a competitive and non-competitive Fur inhibitor for stably interacting with both the FurAct and the FurAll sites.
Feng et al. (2023) [75] highlighted the significance of the potential allosteric site of Fur for the design and repositioning of specific inhibitors of Fur, and other authors have also suggested the feasibility of inhibiting Fur through non-competitive means by allosteric binding of drugs [81,90]. Sheybani et al. (2021) [90] demonstrated that folic acid and folinic acid bind Fur through an H-bond with amino acid residues that belong to the allosteric site (interdomain region) (folic acid: Gly307, Glu271, Tyr313, Gln488, Ala532, Arg490, and Asp530; phenolic acid: Ser311, Glu271, Arg490, Lys449, and Gln488). On the other hand, Peiter et al. (2022) [90] showed that NPF might inhibit Fur by binding Ala532 and Gln488.
The molecular docking data presented here showed slight differences between the AutoDock4 v.4.2.6 and AutoDock Vina software v.2.5.4. This is most probably due to the properties of each program, as explained by Nguyen et al. (2020) [91], who mentioned that AutoDock Vina better illustrates the binding and poses of the ligands with the protein of interest, while AutoDock4 is more accurate in the calculation of the binding affinity between the ligand and the protein. However, molecular docking assays with both programs showed interaction patterns similar to those observed in MD.
According to the molecular interactions, affinity, and stability of complexes, the Ki values calculated by AutoDock4 were consistent with the calculated affinity forces and showed low Ki values for the four molecules studied here, which is related to their high-affinity forces towards their specific interaction sites. PFD was better against FurAct compared to NFP, PGZ, and RGZ, in which Ki values were very low against FurAll and coincided with a high-affinity force of each specific molecule for these sites.
Ki is the inhibitor concentration required to decrease the enzyme’s reaction by half. Therefore, a lower Ki is associated with low or poor dissociation of the inhibitors from the target, i.e., the inhibition is more stable [92].
Therefore, the TZD studied here might have the capacity to bind to Fur as well as NPF and PFD. However, their mechanism of action may not only involve competitive inhibition through binding to the active site (as NPF and PFD do) but also a non-competitive inhibition through destabilization of the interdomain region (catalytic domain and P domain) of Fur. This function of the TZD will be further explored as an antiviral, a function recently observed for NPF against SARS-CoV2 by Cheng et al. (2020) [37]. Further studies are needed to scrutinize this point, as Fur convertase activity plays a key role in the entry of the virus into its host cells during the infectious process [5,49].
Some limitations should be considered for this study. While molecular dynamics simulations provide a valuable starting point for understanding how molecules interact with proteins and assessing their stability, additional in silico analyses like Poisson–Boltzmann, generalized Born, and surface area continuum solvation give more information regarding the union stability. Moreover, enzymatic assays further refine this understanding by confirming both stability and affinity. Future studies should investigate whether these effects hold true under more challenging conditions, as such molecules could offer promising leads for drug development.

5. Conclusions

In conclusion, this study employed in silico methods to identify PGZ and RGZ, which belong to the TZD group, as potential inhibitors of furin. Notably, these molecules have not been previously explored for this purpose. Our analysis revealed favorable binding interactions between these novel inhibitors and furin, similar to those observed with established inhibitors like NPF and PFD. These findings provide a valuable starting point for further research. Experimental testing is now crucial to validate the inhibitory potential of TZD against furin. Additionally, the insights gained regarding NPF and PFD binding modes bring new knowledge for their optimization as furin inhibitors. The potential for repurposing existing drugs like PGZ and RGZ, already tested for safety and side effects in type 2 diabetes patients, could significantly accelerate the development of new therapies. This approach holds promise for diseases with highly unmet medical needs or emerging infectious diseases like SARS-CoV-2. Further research is now necessary to determine their efficacy and safety for specific applications.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/cimb46080511/s1: Figure S1: The 3D representation of redocking of crystallographic structures and co-crystalized ligands. (A) 7QXZ structure with -(3,5-dichlorophenyl) pyridine-derived inhibitor, (B) 5MIM structure with 2,5-dideoxystreptamine-derived inhibitor, (C) 7O1Y with guanylhydrazone-based inhibitor 2 (mi307), and (D) 7LCU structure with BOS-318 inhibitor. Figures obtained with PyMOL v.2.5.4. Blue-green color represents the crystallographic structures with its ligand inhibitors, and pink represents structure redocking with its ligand inhibitors; Figure S2: Representation of the 3D and 2D structures of the molecular interactions between the Fur active site (FurAct) and 3-(3,5-dichlorophenyl) pyridine-derived (3DPP-D) inhibitor. (A) Full complex, (B) specific interactions of 3DPP-D. Figures obtained with the PyMOL v.2.5.4 and Discovery Studio v.21.1.0.20298 software; Figure S3: The 3D and 2D representations of the Fur active site (FurAct) and its molecular interactions with the potential ligands. (A,B) Overview of overlay complexes, (C) NPF, (D) PFD, (E) PGZ, and (F) RGZ. Figures obtained with the PyMOL v.2.5.4 and Discovery Studio v.21.1.0.20298 software; Figure S4: The 3D representation of the models of the complexes between FurAct and its potential ligands during the last 10% of molecular dynamics simulation. (A) NPF (naphtofluorescein), (B) PFD (pirfenidone), (C) PGZ (pioglitazone), and (D) RGZ (rosiglitazone). Figures obtained with PyMOL2; Figure S5: The 2D representation of models of the complex formed between FurAct and its potential ligands in the last 10% of molecular dynamics simulation. (A) NPF (naphtofluorescein), (B) PFD (pirfenidone), (C) PGZ (pioglitazone), and (D) RGZ (rosiglitazone). Figures obtained with the PyMOL v.2.5.4 and LigPlot+ v.2.2.8 software; Figure S6: The 3D representation of models of the complex between FurAll and its potential ligands during the last 10% of molecular dynamics simulation. (A) NPF (naphtofluorescein), (B) PFD (pirfenidone), (C) PGZ (pioglitazone), and (D) RGZ (rosiglitazone). Figures obtained with PyMOL v.2.5.4; Figure S7: The 2D representation of models of the complex formed between FurAll and its potential ligands in the last 10% of molecular dynamics simulation. (A) NPF (naphtofluorescein), (B) PFD (pirfenidone), (C) PGZ (pioglitazone), and (D) RGZ (rosiglitazone). Figures obtained with the PyMOL v.2.5.4 and LigPlot+ v.2.2.8 software; Figure S8: Comparison of the 3D representation of the complex formed between FurAct and FurAll and the potential ligands obtained by molecular docking and molecular dynamics simulation. Figures obtained with Discovery Studio v.2021. Yellow boxes represent the interaction site determined by molecular dockings, and green boxes represent the interaction site determined by molecular dynamics; Table S1: Redocking of structures of furin with their co-crystalized ligands; Table S2: Summary of the molecular interactions of 3DPP-D with FurAct and FurAll; Table S3: Summary of the molecular interactions of NPF, PFD, PGZ, and RGZ with FurAll by general molecular docking.

Author Contributions

A.S.C.-A. and C.D.D.-P. drafted the manuscript and were involved in computational analysis and the interpretation of results. E.L.-R. was involved in computational analysis. A.S. was involved in drafting and revising the manuscript. E.L.-R. and A.S.C.-A. were involved in the data analysis of molecular dynamics. J.F.M.-V. was involved in manuscript revision and study design. O.V.-S. was involved in the study design, computational analysis, results interpretation, and manuscript revision. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and Supplementary Materials. Further inquiries can be directed to the corresponding authors.

Acknowledgments

We appreciate the support given by the staff of LeoAtrox from Universidad de Guadalajara for computational analysis and Itzel A. Ruvalcaba for editing the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Praznik, A.; Fink, T.; Franko, N.; Lonzarić, J.; Benčina, M.; Jerala, N.; Plaper, T.; Roškar, S.; Jerala, R. Regulation of protein secretion through chemical regulation of endoplasmic reticulum retention signal cleavage. Nat. Commun. 2022, 13, 1323. [Google Scholar] [CrossRef] [PubMed]
  2. Solovyeva, N.I.; Gureeva, T.A.; Timoshenko, O.S.; Moskvitina, T.A.; Kugaevskaya, E.V. Furin as proprotein convertase and its role in normal and pathological biological processes. Biochem. Mosc. Suppl. B Biomed. Chem. 2017, 11, 87–100. [Google Scholar] [CrossRef]
  3. Tian, S.; Huang, Q.; Fang, Y.; Wu, J. FurinDB: A Database of 20-residue furin cleavage site motifs, substrates and their associated drugs. Int. J. Mol. Sci. 2011, 12, 1060–1065. [Google Scholar] [CrossRef]
  4. Steiner, D.F. The proprotein convertases. Curr. Opin. Chem. Biol. 1998, 2, 31–39. [Google Scholar] [CrossRef]
  5. Osman, E.E.A.; Rehemtulla, A.; Neamati, N. Why All the Fury over Furin? J. Med. Chem. 2022, 65, 2747–2784. [Google Scholar] [CrossRef] [PubMed]
  6. Seidah, N.G.; Sadr, M.S.; Chrétien, M.; Mbikay, M. The multifaceted proprotein convertases: Their unique, redundant, complementary, and opposite functions. J. Biol. Chem. 2013, 288, 21473–21481. [Google Scholar] [CrossRef]
  7. Cameron, A.; Appel, J.; Houghten, R.A.; Lindberg, I. Polyarginines are potent furin inhibitors. J. Biol. Chem. 2000, 275, 36741–36749. [Google Scholar] [CrossRef]
  8. Braun, E.; Sauter, D. Furin-mediated protein processing in infectious diseases and cancer. Clin. Transl. Immunol. 2019, 8, e1073. [Google Scholar] [CrossRef] [PubMed]
  9. Bassi, D.E.; Mahloogi, H.; De Cicco, R.L.; Klein-Szanto, A. Increased furin activity enhances the malignant phenotype of human head and neck cancer cells. Am. J. Pathol. 2003, 162, 439–447. [Google Scholar] [CrossRef]
  10. He, Z.; Khatib, A.M.; Creemers, J.W.M. Loss of the proprotein convertase Furin in T cells represses mammary tumorigenesis in oncogene-driven triple negative breast cancer. Cancer Lett. 2020, 484, 40–49. [Google Scholar] [CrossRef]
  11. Farhat, D.; Léon, S.; Ghayad, S.E.; Gadot, N.; Icard, P.; Le Romancer, M.; Hussein, N.; Lincet, H. Lipoic acid decreases breast cancer cell proliferation by inhibiting IGF-1R via furin downregulation. Br. J. Cancer. 2020, 122, 885–894. [Google Scholar] [CrossRef] [PubMed]
  12. He, Z.; Thorrez, L.; Siegfried, G.; Meulemans, S.; Evrard, S.; Tejpar, S.; Khatib, A.-M.; Creemers, J.W.M. The proprotein convertase furin is a pro-oncogenic driver in KRAS and BRAF driven colorectal cancer. Oncogene 2020, 39, 3571–3587. [Google Scholar] [CrossRef] [PubMed]
  13. Page, R.E.; Klein-Szanto, A.J.P.; Litwin, S.; Nicolas, E.; Al-Jumaily, R.; Alexander, P.; Godwin, A.K.; Ross, E.A.; Schilder, R.J.; Bassi, D.E. Increased expression of the pro-protein convertase furin predicts decreased survival in ovarian cancer. Cell Oncol. 2007, 29, 289–299. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, Y.; Zhou, M.; Wei, H.; Zhou, H.; He, J.; Lu, Y.; Wang, D.; Chen, B.; Zeng, J.; Peng, W.; et al. Furin promotes epithelial-mesenchymal transition in pancreatic cancer cells via Hippo-YAP pathway. Int. J. Oncol. 2017, 50, 1352–1362. [Google Scholar] [CrossRef] [PubMed]
  15. Tauriello, D.V.F.; Sancho, E.; Batlle, E. Overcoming TGFβ-mediated immune evasion in cancer. Nat. Rev. Cancer. 2022, 22, 25–44. [Google Scholar] [CrossRef] [PubMed]
  16. David, C.J.; Massagué, J. Contextual determinants of TGFβ action in development, immunity and cancer. Nat. Rev. Mol. Cell Biol. 2018, 19, 419–435. [Google Scholar] [CrossRef] [PubMed]
  17. Coppola, J.M.; Hamilton, C.A.; Bhojani, M.S.; Larsen, M.J.; Ross, B.D.; Rehemtulla, A. Identification of inhibitors using a cell-based assay for monitoring Golgi-resident protease activity. Anal. Biochem. 2007, 364, 19–29. [Google Scholar] [CrossRef]
  18. Ren, K.; Jiang, T.; Zheng, X.L.; Zhao, G.J. Proprotein convertase furin/PCSK3 and atherosclerosis: New insights and potential therapeutic targets. Atherosclerosis 2017, 262, 163–170. [Google Scholar] [CrossRef] [PubMed]
  19. Lin, H.; Kioon, M.-D.A.; Lalou, C.; Larghero, J.; Launay, J.-M.; Khatib, A.-M.; Cohen-Solal, M. Protective role of systemic furin in immune response-induced arthritis. Arthritis Rheum. 2012, 64, 2878–2886. [Google Scholar] [CrossRef]
  20. Cao, R.; Zhang, Y.; Du, J.; Chen, S.; Wang, N.; Ying, H.; Shen, B. Increased FURIN expression in rheumatoid arthritis patients and its anti-inflammatory effect. J. Clin. Lab. Anal. 2020, 34, e23530. [Google Scholar] [CrossRef]
  21. Wu, T.; Ding, H.; Han, J.; Arriens, C.; Wei, C.; Han, W.; Pedroza, C.; Jiang, S.; Anolik, J.; Petri, M.; et al. Antibody-Array-Based Proteomic Screening of Serum Markers in Systemic Lupus Erythematosus: A Discovery Study. J. Proteome Res. 2016, 15, 2102–2114. [Google Scholar] [CrossRef] [PubMed]
  22. Ranta, N.; Valli, A.; Grönholm, A.; Silvennoinen, O.; Isomäki, P.; Pesu, M.; Pertovaara, M. Proprotein convertase enzyme FURIN is upregulated in primary Sjögren’s syndrome. Clin. Exp. Rheumatol. 2018, 36 (Suppl. S112), 47–50. Available online: https://europepmc.org/article/med/29465367 (accessed on 30 July 2024).
  23. Gordon, V.M.; Benz, R.; Fujii, K.; Leppla, S.H.; Tweten, R.K. Clostridium septicum alpha-toxin is proteolytically activated by furin. Infect. Immun. 1997, 65, 4130–4134. [Google Scholar] [CrossRef]
  24. Shiryaev, S.A.; Remacle, A.G.; Ratnikov, B.I.; Nelson, N.A.; Savinov, A.Y.; Wei, G.; Bottini, M.; Rega, M.F.; Parent, A.; Desjardins, R.; et al. Targeting host cell furin proprotein convertases as a therapeutic strategy against bacterial toxins and viral pathogens. J. Biol. Chem. 2007, 282, 20847–20853. [Google Scholar] [CrossRef]
  25. Couture, F.; Kwiatkowska, A.; Dory, Y.L.; Day, R. Therapeutic uses of furin and its inhibitors: A patent review. Expert Opin. Ther. Pat. 2015, 25, 379–396. [Google Scholar] [CrossRef]
  26. Logeat, F.; Bessia, C.; Brou, C.; LeBail, O.; Jarriault, S.; Seidah, N.G.; Israël, A. The Notch1 receptor is cleaved constitutively by a furin-like convertase. Proc. Natl. Acad. Sci. USA 1998, 95, 8108–8112. [Google Scholar] [CrossRef]
  27. Izaguirre, G. The Proteolytic Regulation of Virus Cell Entry by Furin and Other Proprotein Convertases. Viruses 2019, 11, 837. [Google Scholar] [CrossRef] [PubMed]
  28. Pierson, T.C.; Diamond, M.S. Degrees of maturity: The complex structure and biology of flaviviruses. Curr. Opin. Virol. 2012, 2, 168–175. [Google Scholar] [CrossRef] [PubMed]
  29. Renner, M.; Dejnirattisai, W.; Carrique, L.; Martin, I.S.; Karia, D.; Ilca, S.L.; Ho, S.F.; Kotecha, A.; Keown, J.R.; Mongkolsapaya, J.; et al. Flavivirus maturation leads to the formation of an occupied lipid pocket in the surface glycoproteins. Nat. Commun. 2021, 12, 1238. [Google Scholar] [CrossRef]
  30. Pang, Y.J.; Tan, X.J.; Li, D.M.; Zheng, Z.H.; Lei, R.X.; Peng, X.M. Therapeutic potential of furin inhibitors for the chronic infection of hepatitis B virus. Liver International. 2013, 33, 1230–1238. [Google Scholar] [CrossRef]
  31. Messageot, F.; Salhi, S.; Eon, P.; Rossignol, J.M. Proteolytic processing of the Hepatitis B virus e antigen precursor: Cleavage at two furin consensus sequences. J. Biol. Chem. 2003, 278, 891–895. [Google Scholar] [CrossRef]
  32. Tse, L.V.; Hamilton, A.M.; Friling, T.; Whittaker, G.R. A Novel Activation Mechanism of Avian Influenza Virus H9N2 by Furin. J. Virol. 2014, 88, 1673–1683. [Google Scholar] [CrossRef] [PubMed]
  33. Wu, Y.; Zhao, S. Furin cleavage sites naturally occur in coronaviruses. Stem Cell Res. 2021, 50, 102115. [Google Scholar] [CrossRef]
  34. Johnson, B.A.; Xie, X.; Kalveram, B.; Lokugamage, K.G.; Muruato, A.; Zou, J.; Zhang, X.; Juelich, T.; Smith, J.K.; Zhang, L.; et al. Furin Cleavage Site Is Key to SARS-CoV-2 Pathogenesis. bioRxiv. 2020, preprint. [Google Scholar] [CrossRef]
  35. Johnson, B.A.; Xie, X.; Bailey, A.L.; Kalveram, B.; Lokugamage, K.G.; Muruato, A.; Zou, J.; Zhang, X.; Juelich, T.; Smith, J.K.; et al. Loss of furin cleavage site attenuates SARS-CoV-2 pathogenesis. Nature 2021, 591, 293–299. [Google Scholar] [CrossRef] [PubMed]
  36. Chan, Y.A.; Zhan, S.H. The Emergence of the Spike Furin Cleavage Site in SARS-CoV-2. Mol. Biol. Evol. 2022, 39, msab327. [Google Scholar] [CrossRef]
  37. Cheng, Y.-W.; Chao, T.-L.; Li, C.-L.; Chiu, M.-F.; Kao, H.-C.; Wang, S.-H.; Pang, Y.-H.; Lin, C.-H.; Tsai, Y.-M.; Lee, W.-H.; et al. Furin Inhibitors Block SARS-CoV-2 Spike Protein Cleavage to Suppress Virus Production and Cytopathic Effects. Cell Rep. 2020, 33, 108254. [Google Scholar] [CrossRef] [PubMed]
  38. Margaritopoulos, G.A.; Kokosi, M.A.; Wells, A.U. Diagnosing complications and co-morbidities of fibrotic interstitial lung disease. Expert Rev. Respir. Med. 2019, 13, 645–658. [Google Scholar] [CrossRef]
  39. Seifirad, S. Pirfenidone: A novel hypothetical treatment for COVID-19. Med. Hypotheses 2020, 144, 110005. [Google Scholar] [CrossRef]
  40. Burghardt, I.; Tritschler, F.; Opitz, C.A.; Frank, B.; Weller, M.; Wick, W. Pirfenidone inhibits TGF-β expression in malignant glioma cells. Biochem. Biophys. Res. Commun. 2007, 354, 542–547. [Google Scholar] [CrossRef]
  41. Villoutreix, B.O.; Creemers, J.; Léger, Y.; Siegfried, G.; Decroly, E.; Evrard, S.; Khatib, A.M. Targeting furin activity through in silico and in vitro drug repurposing strategy for SARS-CoV-2 spike glycoprotein cleavage repression. Res. Sq. 2020, preprint. [Google Scholar] [CrossRef]
  42. Miller, B.W.; Willett, K.C.; Desilets, A.R. Rosiglitazone and Pioglitazone for the Treatment of Alzheimer’s Disease. Ann. Pharmacother. 2011, 45, 1416–1424. [Google Scholar] [CrossRef] [PubMed]
  43. Orasanu, G.; Ziouzenkova, O.; Devchand, P.R.; Nehra, V.; Hamdy, O.; Horton, E.S.; Plutzky, J. The Peroxisome Proliferator-Activated Receptor-γ Agonist Pioglitazone Represses Inflammation in a Peroxisome Proliferator-Activated Receptor-α-Dependent Manner In vitro and In vivo in Mice. J. Am. Coll. Cardiol. 2008, 52, 869–881. [Google Scholar] [CrossRef]
  44. Tilekar, K.; Shelke, O.; Upadhyay, N.; Lavecchia, A.; Ramaa, C.S. Current status and future prospects of molecular hybrids with thiazolidinedione (TZD) scaffold in anticancer drug discovery. J. Mol. Struct. 2022, 1250, 131767. [Google Scholar] [CrossRef]
  45. Adeshara, K.A.; Agrawal, S.B.; Gaikwad, S.M.; Tupe, R.S. Pioglitazone inhibits advanced glycation induced protein modifications and down-regulates expression of RAGE and NF-κB in renal cells. Int. J. Biol. Macromol. 2018, 119, 1154–1163. [Google Scholar] [CrossRef]
  46. Nhu, Q.M.; Hsieh, L.; Dohil, L.; Dohil, R.; Newbury, R.O.; Kurten, R.; Moawad, F.J.; Aceves, S.S. Antifibrotic Effects of the Thiazolidinediones in Eosinophilic Esophagitis Pathologic Remodeling: A Preclinical Evaluation. Clin. Transl. Gastroenterol. 2020, 11, e00164. [Google Scholar] [CrossRef]
  47. Wei, Z.; Zhao, D.; Zhang, Y.; Chen, Y.; Zhang, S.; Li, Q.; Zeng, P.; Li, X.; Zhang, W.; Duan, Y.; et al. Rosiglitazone ameliorates bile duct ligation-induced liver fibrosis by down-regulating NF-κB-TNF-α signaling pathway in a PPARγ-dependent manner. Biochem. Biophys. Res. Commun. 2019, 519, 854–860. [Google Scholar] [CrossRef] [PubMed]
  48. Németh, Á.; Mózes, M.M.; Calvier, L.; Hansmann, G.; Kökény, G. The PPARγagonist pioglitazone prevents TGF-β induced renal fibrosis by repressing EGR-1 and STAT3. BMC Nephrol. 2019, 20, 245. [Google Scholar] [CrossRef] [PubMed]
  49. Setz, C.; Große, M.; Auth, J.; Fröba, M.; Rauch, P.; Bausch, A.; Wright, M.; Schubert, U. Synergistic Antiviral Activity of Pamapimod and Pioglitazone against SARS-CoV-2 and Its Variants of Concern. Int. J. Mol. Sci. 2022, 23, 6830. [Google Scholar] [CrossRef]
  50. Dahms, S.O.; Schnapp, G.; Winter, M.; Büttner, F.H.; Schlepütz, M.; Gnamm, C.; Pautsch, A.; Brandstetter, H. Dichlorophenylpyridine-Based Molecules Inhibit Furin through an Induced-Fit Mechanism. ACS Chem. Biol. 2022, 17, 816–821. [Google Scholar] [CrossRef]
  51. Douglas, L.E.; Reihill, J.A.; Ho, M.W.; Axten, J.M.; Campobasso, N.; Schneck, J.L.; Rendina, A.R.; Wilcoxen, K.M.; Martin, S.L. A highly selective, cell-permeable furin inhibitor BOS-318 rescues key features of cystic fibrosis airway disease. Cell Chem. Biol. 2022, 29, 947–957. [Google Scholar] [CrossRef] [PubMed]
  52. Dahms, S.O.; Jiao, G.S.; Than, M.E. Structural Studies Revealed Active Site Distortions of Human Furin by a Small Molecule Inhibitor. ACS Chem. Biol. 2017, 12, 1211–1216. [Google Scholar] [CrossRef] [PubMed]
  53. Dahms, S.O.; Haider, T.; Klebe, G.; Steinmetzer, T.; Brandstetter, H. OFF-State-Specific inhibition of the proprotein convertase furin. ACS Chem. Biol. 2021, 16, 1692–1700. [Google Scholar] [CrossRef]
  54. Irwin, J.J.; Shoichet, B.K. ZINC—A free database of commercially available compounds for virtual screening. J. Chem. Inf. Model. 2005, 45, 177–182. [Google Scholar] [CrossRef] [PubMed]
  55. Jakalian, A.; Jack, D.B.; Bayly, C.I. Fast, efficient generation of high-quality atomic charges. AM1-BCC model: II. Parameterization and validation. J. Comput. Chem. 2002, 23, 1623–1641. [Google Scholar] [CrossRef] [PubMed]
  56. Gasteiger, J.; Marsili, M. Iterative partial equalization of orbital electronegativity-a rapid access to atomic charges. Tetrahedron 1980, 36, 3219–3228. [Google Scholar] [CrossRef]
  57. Morris, G.M.; Goodsell, D.S.; Halliday, R.S.; Huey, R.; Hart, W.E.; Belew, R.K.; Olson, A.J. Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. J. Comput. Chem. 1998, 19, 1639–1662. [Google Scholar] [CrossRef]
  58. Jaqueline, A.; Zúñiga, P.; Rojo Domínguez, A. SIMULACIÓN DEL RECONOCIMIENTO ENTRE PROTEÍNAS Y MOLÉCULAS ORGÁNICAS O DOCKING. APLICACIÓN AL DISEÑO DE FÁRMACOS. Mensaje Bioquímico 2002, 26, 129–145. [Google Scholar]
  59. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera—A visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef] [PubMed]
  60. Plewczynski, D.; Łaźniewski, M.; Augustyniak, R.; Ginalski, K. Can we trust docking results? Evaluation of seven commonly used programs on PDBbind database. J. Comput. Chem. 2011, 32, 742–755. [Google Scholar] [CrossRef]
  61. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef] [PubMed]
  62. Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 2009, 30, 2785–2791. [Google Scholar] [CrossRef] [PubMed]
  63. Sánchez-Linares, I.; Pérez-Sánchez, H.; Cecilia, J.M.; García, J.M. High-Throughput parallel blind Virtual Screening using BINDSURF. BMC Bioinform. 2012, 13 (Suppl. S14), S13. [Google Scholar] [CrossRef] [PubMed]
  64. Schrödinger, L.; DeLano, W. PyMOL. 2020. Available online: http://www.pymol.org/pymol (accessed on 30 July 2024).
  65. McRee, D.E. 3—Computational techniques. In Practical Protein Cristalographyc, 2nd ed.; Academic Press: Cambridge, MA, USA, 1999. [Google Scholar]
  66. Kalé, L.; Skeel, R.; Bhandarkar, M.; Brunner, R.; Gursoy, A.; Krawetz, N.; Phillips, J.; Shinozaki, A.; Varadarajan, K.; Schulten, K. NAMD2: Greater Scalability for Parallel Molecular Dynamics. J. Comput. Phys. 1999, 151, 283–312. [Google Scholar] [CrossRef]
  67. Jo, S.; Kim, T.; Iyer, V.G.; Im, W. CHARMM-GUI: A web-based graphical user interface for CHARMM. J. Comput. Chem. 2008, 29, 1859–1865. [Google Scholar] [CrossRef]
  68. Lee, J.; Cheng, X.; Swails, J.M.; Yeom, M.S.; Eastman, P.K.; Lemkul, J.A.; Wei, S.; Buckner, J.; Jeong, J.C.; Qi, Y.; et al. CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field. J. Chem. Theory Comput. 2016, 12, 405–413. [Google Scholar] [CrossRef] [PubMed]
  69. Humphrey, W.; Andrew, D.; Klaus, S. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef] [PubMed]
  70. Glykos, N.M. Software news and updates carma: A molecular dynamics analysis program. J. Comput. Chem. 2006, 27, 1765–1768. [Google Scholar] [CrossRef] [PubMed]
  71. Reva, B.A.; Finkelstein, A.V.; Sanner, M.F.; Olson, A.J. Residue-residue mean-force potentials for protein structure recognition. Protein Eng. 1997, 10, 865–876. [Google Scholar] [CrossRef]
  72. Abdalla, M.; Eltayb, W.A.; El-Arabey, A.A.; Singh, K.; Jiang, X. Molecular dynamic study of SARS-CoV-2 with various S protein mutations and their effect on thermodynamic properties. Comput. Biol. Med. 2022, 141, 105025. [Google Scholar] [CrossRef]
  73. Anwar, M.A.; Choi, S. Structure-activity relationship in tlr4 mutations: Atomistic molecular dynamics simulations and residue interaction network analysis. Sci. Rep. 2017, 7, 43807. [Google Scholar] [CrossRef] [PubMed]
  74. Thomas, G.; Couture, F.; Kwiatkowska, A. The Path to Therapeutic Furin Inhibitors: From Yeast Pheromones to SARS-CoV-2. Int. J. Mol. Sci. 2022, 23, 3435. [Google Scholar] [CrossRef] [PubMed]
  75. Feng, D.; Ren, L.; Wu, J.; Guo, L.; Han, Z.; Yang, J.; Xie, W.; Wang, Y.; Xu, F.; Su, X.; et al. Permethrin as a Potential Furin Inhibitor through a Novel Non-Competitive Allosteric Inhibition. Molecules 2023, 28, 1883. [Google Scholar] [CrossRef] [PubMed]
  76. Devi, K.P.; Pourkarim, M.R.; Thijssen, M.; Sureda, A.; Khayatkashani, M.; Cismaru, C.A.; Neagoe, I.B.; Habtemariam, S.; Razmjouei, S.; Kashani, H.R.K. A perspective on the applications of furin inhibitors for the treatment of SARS-CoV-2. Pharmacol. Rep. 2022, 74, 425–430. [Google Scholar] [CrossRef] [PubMed]
  77. Sakamoto, T.; Fukui, Y.; Kondoh, Y.; Honda, K.; Shimizu, T.; Hara, T.; Hayashi, T.; Saitoh, Y.; Murakami, Y.; Inoue, J.-I.; et al. Pharmacological inhibition of Mint3 attenuates tumour growth, metastasis, and endotoxic shock. Commun. Biol. 2021, 4, 1165. [Google Scholar] [CrossRef]
  78. Ramos-Molina, B.; Lick, A.N.; Blanco, E.H.; Posada-Salgado, J.A.; Martinez-Mayorga, K.; Johnson, A.T.; Jiao, G.-S.; Lindberg, I. Identification of potent and compartment-selective small molecule furin inhibitors using cell-based assays. Biochem. Pharmacol. 2015, 96, 107–118. [Google Scholar] [CrossRef] [PubMed]
  79. Majumdar, S.; Mohanta, B.C.; Chowdhury, D.R.; Banik, R.; Dinda, B.; Basak, A. Proprotein Convertase Inhibitory Activities of Flavonoids Isolated from Oroxylum Indicum. Curr. Med. Chem. 2010, 17, 2049–2058. [Google Scholar] [CrossRef]
  80. Tanaka, N.; Sakamoto, T. Mint3 as a Potential Target for Cooling Down HIF-1α-Mediated Inflammation and Cancer Aggressiveness. Biomedicines 2023, 11, 549. [Google Scholar] [CrossRef] [PubMed]
  81. Peiter, G.C.; de Bem Torquato de Souza, C.; de Oliveira, L.M.; Pagliarin, L.G.; Anjos, V.N.F.D.; da Silva, F.A.F.; de Melo, F.F.; Teixeira, K.N. COVID-19 liver and gastroenterology findings: An in silico analysis of SARS-CoV-2 interactions with liver molecules. World J. Hepatol. 2022, 14, 1131–1141. [Google Scholar] [CrossRef]
  82. Henrich, S.; Cameron, A.; Bourenkov, G.P.; Kiefersauer, R.; Huber, R.; Lindberg, I.; Bode, W.; Than, M.E. The crystal structure of the proprotein processing proteinase furin explains its stringent specificity. Nat. Struct. Biol. 2003, 10, 520–526. [Google Scholar] [CrossRef]
  83. Pearce, K.H.; Overton, L.K.; Gampe, R.T.; Barrett, G.B.; Taylor, J.D.; McKee, D.D.; Campobasso, N.; Nolte, R.T.; Reid, R.A. BacMam production and crystal structure of nonglycosylated apo human furin at 1.89 Å resolution. Acta Crystallogr. F Struct. Biol. Commun. 2019, 75, 239–245. [Google Scholar] [CrossRef] [PubMed]
  84. Jiao, J.; Friedman, S.L.; Aloman, C. Hepatic fibrosis. Curr. Opin. Gastroenterol. 2009, 25, 223–229. [Google Scholar] [CrossRef] [PubMed]
  85. Dahms, S.O.; Arciniega, M.; Steinmetzer, T.; Huber, R.; Than, M.E. Structure of the unliganded form of the proprotein convertase furin suggests activation by a substrate-induced mechanism. Proc. Natl. Acad. Sci. USA 2016, 113, 11196–11201. [Google Scholar] [CrossRef] [PubMed]
  86. Ferreira De Freitas, R.; Schapira, M. A systematic analysis of atomic protein-ligand interactions in the PDB. Medchemcomm 2017, 8, 1970–1981. [Google Scholar] [CrossRef] [PubMed]
  87. Klebe, G. Protein–Ligand Interactions as the Basis for Drug Action. In Multifaceted Roles of Crystallography in Modern Drug Discovery; Springer: Dordrecht, The Netherlands, 2013; pp. 83–92. [Google Scholar]
  88. Jampilek, J. Heterocycles in Medicinal Chemistry. Molecules 2019, 24, 3839. [Google Scholar] [CrossRef] [PubMed]
  89. Long, N.; Le Gresley, A.; Wren, S.P. Thiazolidinediones: An In–Depth Study of Their Synthesis and Application to Medicinal Chemistry in the Treatment of Diabetes Mellitus. ChemMedChem 2021, 16, 1717–1736. [Google Scholar] [CrossRef] [PubMed]
  90. Sheybani, Z.; Dokoohaki, M.H.; Negahdaripour, M.; Dehdashti, M.; Zolghadr, H.; Moghadami, M.; Masoompour, S.M.; Zolghadr, A.R. The interactions of folate with the enzyme furin: A computational study. RSC Adv. 2021, 11, 23815. [Google Scholar] [CrossRef] [PubMed]
  91. Nguyen, N.T.; Nguyen, T.H.; Pham, T.N.H.; Huy, N.T.; Van Bay, M.; Pham, M.Q.; Nam, P.C.; Vu, V.V.; Ngo, S.T. Autodock Vina Adopts More Accurate Binding Poses but Autodock4 Forms Better Binding Affinity. J. Chem. Inf. Model. 2020, 60, 204–211. [Google Scholar] [CrossRef]
  92. Pandey, A.K.; Verma, S. An in-silico evaluation of dietary components for structural inhibition of SARS-Cov-2 main protease. J. Biomol. Struct. Dyn. 2022, 40, 136–142. [Google Scholar] [CrossRef]
Cimb 46 00511 g001
Figure 1. Workflow diagram.
Figure 1. Workflow diagram.
Cimb 46 00511 g001
Cimb 46 00511 g002
Figure 2. Representation of furin’s 3D general structure and domains (PDB ID: 7XQZ). Figure obtained with PyMOL v.2.5.4.
Figure 2. Representation of furin’s 3D general structure and domains (PDB ID: 7XQZ). Figure obtained with PyMOL v.2.5.4.
Cimb 46 00511 g002
Cimb 46 00511 g003
Figure 3. The 3D and 2D representations of the Fur active site (FurAct) and its molecular interactions with the potential ligands. (A,B) Overview of overlay complexes, (C) NPF, (D) PFD, (E) PGZ, and (F) RGZ. Figures obtained with the PyMOL v.2.5.4 and Discovery Studio v.21.1.0.20298 software.
Figure 3. The 3D and 2D representations of the Fur active site (FurAct) and its molecular interactions with the potential ligands. (A,B) Overview of overlay complexes, (C) NPF, (D) PFD, (E) PGZ, and (F) RGZ. Figures obtained with the PyMOL v.2.5.4 and Discovery Studio v.21.1.0.20298 software.
Cimb 46 00511 g003
Cimb 46 00511 g004
Figure 4. The 3D and 2D representations of the Fur allosteric site (FurAll) and its molecular interactions with the potential ligands. (A,B) Overview of overlay complexes, (C) NPF, (D) PFD, (E) PGZ, and (F) RGZ. Figures obtained with the PyMOL v.2.5.4 and Discovery Studio v.21.1.0.20298 software.
Figure 4. The 3D and 2D representations of the Fur allosteric site (FurAll) and its molecular interactions with the potential ligands. (A,B) Overview of overlay complexes, (C) NPF, (D) PFD, (E) PGZ, and (F) RGZ. Figures obtained with the PyMOL v.2.5.4 and Discovery Studio v.21.1.0.20298 software.
Cimb 46 00511 g004
Cimb 46 00511 g005
Figure 5. RMSD of ligands/Fur complexes during MD simulation with (A) active site and (B) allosteric site of Fur. RMSD of complexes of NPF (naphtofluorescein), PFD (pirfenidone), PGZ (pioglitazone), and RGZ (rosiglitazone).
Figure 5. RMSD of ligands/Fur complexes during MD simulation with (A) active site and (B) allosteric site of Fur. RMSD of complexes of NPF (naphtofluorescein), PFD (pirfenidone), PGZ (pioglitazone), and RGZ (rosiglitazone).
Cimb 46 00511 g005
Cimb 46 00511 g006
Figure 6. Rg of ligands/Fur complexes during MD simulation with (A) active site and (B) allosteric site of Fur. Rg of complexes of NPF (naphtofluorescein), PFD (pirfenidone), PGZ (pioglitazone), and RGZ (rosiglitazone).
Figure 6. Rg of ligands/Fur complexes during MD simulation with (A) active site and (B) allosteric site of Fur. Rg of complexes of NPF (naphtofluorescein), PFD (pirfenidone), PGZ (pioglitazone), and RGZ (rosiglitazone).
Cimb 46 00511 g006
Cimb 46 00511 g007
Figure 7. Root mean square fluctuation (RMSF) values of ligand/Fur complexes within the Fur active site (FurAct; (A)) and Fur allosteric site (FurAll; (B)) during MD simulation for 100 ns. NPF: naphthofluorescein; PFD: pirferidone, PGZ: pioglitazona, RGZ: rosiglitazone. Figures obtained with Carma v.2.01 and Excel software v. 2406. Arrows indicate the critical residues on the Fur active site. The purple box indicates the catalytic domain. The blue box indicates the P domain.
Figure 7. Root mean square fluctuation (RMSF) values of ligand/Fur complexes within the Fur active site (FurAct; (A)) and Fur allosteric site (FurAll; (B)) during MD simulation for 100 ns. NPF: naphthofluorescein; PFD: pirferidone, PGZ: pioglitazona, RGZ: rosiglitazone. Figures obtained with Carma v.2.01 and Excel software v. 2406. Arrows indicate the critical residues on the Fur active site. The purple box indicates the catalytic domain. The blue box indicates the P domain.
Cimb 46 00511 g007
Table 1. Zinc15 database IDs of the 3D structures of four potential Fur ligands: NPF, PFD, PGZ, and RGZ.
Table 1. Zinc15 database IDs of the 3D structures of four potential Fur ligands: NPF, PFD, PGZ, and RGZ.
IDName and AbbreviationStructure
PubChem CID 3124834Naphthofluorescein
NPF		Cimb 46 00511 i001
ZINC1958Pirfenidone
PFD		Cimb 46 00511 i002
ZINC968326Pioglitazone
PGZ		Cimb 46 00511 i003
ZINC968328Rosiglitazone
RGZ		Cimb 46 00511 i004
NPF: naphtofluorescein; PFD: pirfenidone; PGZ: pioglitazone; RGZ: rosiglitazone.
Table 2. Summary of NPF, PFD, PGZ, and RGZ molecular interactions with FurAct by localized molecular docking.
Table 2. Summary of NPF, PFD, PGZ, and RGZ molecular interactions with FurAct by localized molecular docking.
LigandsSite
Interaction		Free Energy of Binding (kcal/mol)
AutoDock4		Free Energy of Binding (kcal/mol)
Vina		H-Bond Classical and Non-ClassicalDistance H-Bond (Å)ElectrostaticHydrophobic
Naphthofluorescein
(NPF)		ActiveNA
(NA)		−9.9
(0.00)		Asp154—H
Asp191—H		2.7
3.1		Asp258Val231
Leu227
Pirfenidone
(PFD)		Active−6.28
(0.00)		−6.9
(1.709)		Gly255—O3.4Glu236Gly255
Trp254
Trp291
Leu227
Tyr308
Val231
Pro256
Pioglitazone
(PGZ)		Active−8.04
(0.00)		−8.6
(0.00)		Asp154—H
Asp191—H
Asn192—O		2.6
2.5
2.6		NALeu227
Val231
Leu152
Met226
Trp254
Leu240
Tyr308
Trp291
Rosiglitazone
(RGZ)		Active−7.47
(0.00)		−8.1
(0.00)		Asp154—H2.2NALeu227
Val231
Met226
Trp291
Ala252
Tyr308
NA: not available; Å = angstrom; NPF: naphthofluorescein; PFD = pirfenidone; PGZ: pioglitazone; RGZ: rosiglitazone.
Table 3. Summary of molecular interactions of NPF, PFD, PGZ, and RGZ with FurAll by blind docking.
Table 3. Summary of molecular interactions of NPF, PFD, PGZ, and RGZ with FurAll by blind docking.
LigandsSite
Interaction		Free Energy of Binding (kcal/mol)
Blind Docking Server		H-Bond Classical and Non-ClassicalDistance H-Bond (Å)ElectrostaticHydrophobic
Naphthofluorescein
(NPF)		Allosteric−10.9Gln488—H
Asn310—O
Gly265—O		2.3
2.4
2.3		NAPro266
Trp531
Ser311
Ala532
Glu271
Pirfenidone
(PFD)		Allosteric−7.1Asn310—O2.6Glu271Ala532
Trp531
Tyr313
Ile312
Pioglitazone
(PGZ)		Allosteric−8.7Ser311—H
Gln488—O		2.2
2.7		NATrp531
Val263
Ile312
Rosiglitazone
(RGZ)		Allosteric−7.9Asn310—H
Ser311—O
Val263—N
Glu271—O		1.9
2.2
3.3
3.4		NATrp531
Val263
Asp264
Gly307
NA: not available; Å = angstrom; NPF: naphthofluorescein; PFD = pirfenidone; PGZ: pioglitazone; RGZ: rosiglitazone.
Table 4. Inhibition constant (Ki) of Fur by NFP, PFD, PGZ, and RGZ.
Table 4. Inhibition constant (Ki) of Fur by NFP, PFD, PGZ, and RGZ.
Inhibition Constant (Ki)
LigandsLocalized Docking
(FurAct)		Blind Docking
(FurAll)
NFPNA201.04 nM
PFD24.87 μMNA
PGZ1.28 μM970.36 nM
RGZ3.33 μM1.49 μM
NA: not available. NFP: naphtofluorescein; PFD: pirfenidone; PGZ: pioglitazone; RGZ: rosiglitazone.
Table 5. Average RMSD of the ligands/FurAct and ligands/FurAll complexes during the last 10% of the molecular dynamics.
Table 5. Average RMSD of the ligands/FurAct and ligands/FurAll complexes during the last 10% of the molecular dynamics.
ComplexesRMSD (Å) of Active Site (100 ns)RMSD (Å) of Allosteric Site (100 ns)
Ligand-free Fur3.98 ± 0.073.98 ± 0.07
NPF/Fur4.67 ± 0.134.52 ± 0.05
PFD/Fur4.21 ± 0.085.22 ± 0.17
PGZ/Fur4.12 ± 0.064.66 ± 0.17
RGZ/Fur4.77 ± 0.145.57 ± 0.17
NPF: naphtofluorescein; PFD: pirfenidone; PGZ: pioglitazone; RGZ: rosiglitazone.
Table 6. Average Rg of ligand/Fur complexes between the Fur active and allosteric sites during the last 10% of the molecular dynamics.
Table 6. Average Rg of ligand/Fur complexes between the Fur active and allosteric sites during the last 10% of the molecular dynamics.
ComplexesRg (Å) of Active Site
(100 ns)		Rg (Å) of Allosteric Site
(100 ns)
Ligand-free Fur22.45 ± 0.0622.4 ± 0.06
NPF/Fur22.62 ± 0.0922.56 ± 0.06
PFD/Fur22.57 ± 0.0522.90 ± 0.06
PGZ/Fur22.65 ± 0.0622.93 ± 0.06
RGZ/Fur22.77 ± 0.1023.22 ± 0.06
NPF: naphtofluorescein; PFD: pirfenidone; PGZ: pioglitazone; RGZ: rosiglitazone.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Carranza-Aranda, A.S.; Diaz-Palomera, C.D.; Lepe-Reynoso, E.; Santerre, A.; Muñoz-Valle, J.F.; Viera-Segura, O. Evaluation of Potential Furin Protease Inhibitory Properties of Pioglitazone, Rosiglitazone, and Pirfenidone: An In Silico Docking and Molecular Dynamics Simulation Approach. Curr. Issues Mol. Biol. 2024, 46, 8665-8684. https://doi.org/10.3390/cimb46080511

AMA Style

Carranza-Aranda AS, Diaz-Palomera CD, Lepe-Reynoso E, Santerre A, Muñoz-Valle JF, Viera-Segura O. Evaluation of Potential Furin Protease Inhibitory Properties of Pioglitazone, Rosiglitazone, and Pirfenidone: An In Silico Docking and Molecular Dynamics Simulation Approach. Current Issues in Molecular Biology. 2024; 46(8):8665-8684. https://doi.org/10.3390/cimb46080511

Chicago/Turabian Style

Carranza-Aranda, Ahtziri Socorro, Carlos Daniel Diaz-Palomera, Eduardo Lepe-Reynoso, Anne Santerre, José Francisco Muñoz-Valle, and Oliver Viera-Segura. 2024. "Evaluation of Potential Furin Protease Inhibitory Properties of Pioglitazone, Rosiglitazone, and Pirfenidone: An In Silico Docking and Molecular Dynamics Simulation Approach" Current Issues in Molecular Biology 46, no. 8: 8665-8684. https://doi.org/10.3390/cimb46080511

APA Style

Carranza-Aranda, A. S., Diaz-Palomera, C. D., Lepe-Reynoso, E., Santerre, A., Muñoz-Valle, J. F., & Viera-Segura, O. (2024). Evaluation of Potential Furin Protease Inhibitory Properties of Pioglitazone, Rosiglitazone, and Pirfenidone: An In Silico Docking and Molecular Dynamics Simulation Approach. Current Issues in Molecular Biology, 46(8), 8665-8684. https://doi.org/10.3390/cimb46080511

Article Metrics

Back to TopTop