Skip to main content

ORIGINAL RESEARCH article

Front. Chem., 12 October 2021
Sec. Medicinal and Pharmaceutical Chemistry
This article is part of the Research Topic Computer-Aided Drug Design: Drug Discovery, Computational Modelling, and Artificial Intelligence View all 15 articles

Virtual Screening on Marine Natural Products for Discovering TMPRSS2 Inhibitors

  • 1The Persian Gulf Marine Biotechnology Research Center, The Persian Gulf Biomedical Sciences Research Institute, Bushehr University of Medical Sciences, Bushehr, Iran
  • 2The Persian Gulf Tropical Medicine Research Center, The Persian Gulf Biomedical Sciences Research Institute, Bushehr University of Medical Sciences, Bushehr, Iran

Although SARS-CoV-2 entry to cells strictly depends on angiotensin-converting enzyme 2 (ACE2), the virus also needs transmembrane serine protease 2 (TMPRSS2) for its spike protein priming. It has been shown that the entrance of SARS-CoV-2 through ACE2 can be blocked by cellular TMPRSS2 blockers. The main aim of this study was to find potential inhibitor(s) of TMPRSS2 through virtual screening against a homology model of TMPRSS2 using the library of marine natural products (MNPs). The homology modeling technique for generating a three-dimensional structure of TMPRSS2 was applied. Molecular docking, MM-GBSA and absorption, distribution, metabolism, excretion (ADME) evaluations were performed to investigate the inhibitory activity of marine natural products (MNPs) against TMPRSS2 and their pharmacokinetic properties. Camostat and nafamostat mesylate were used as the standard inhibitory molecules. Seven MNPs were able to inhibit TMPRSS2 better than the standard compounds. MNP 10 with CAS number 107503-09-3, called Watasenia β-D- Preluciferyl glucopyrasoiuronic acid, was found to be the best inhibitor of TMPRSS2 with acceptable pharmacokinetic properties. Herein, for the first time, a new marine natural product was introduced with potent inhibitory effects against TMPRSS2. MNP 10 exhibited favorable drug-like pharmacokinetic properties and it promises a novel TMPRSS2 blocker to combat SARS-CoV-2.

Highlights

1) Marine natural products (MNPs) are a valuable source for anti-SARS-CoV-2 drugs.

2) MNP 10 is a potent TMPRSS2 inhibitor to combat SARS-CoV-2.

3) MNP 10 has favorable drug-like pharmacokinetic characteristics.

Introduction

The devastating pandemic caused by SARS-CoV-2 (Astuti, 2020; Guan et al., 2020; Wang et al., 2020; Zhu et al., 2020) that first broke out in Wuhan, China in late 2019, has become the most important global health and socioeconomic issue. Although there is a worldwide effort to develop an effective vaccine against SARS-CoV-2 using both established and new vaccine production technologies (WHO, 2021), no one can yet claim what kind of therapy can be absolutely efficient for the treatment of or protection against COVID-19. Based on the involved pathophysiologic pathways, different kinds of therapeutic modalities have been conducted in numerous clinical trials with conflicting results.

SARS-CoV-2 is a spherical shaped virus with a diameter of about 60–140 nm with some pleomorphism belonging to Coronaviride family. The enveloped virus has distinctive spikes (Zhu et al., 2020). Its genome is around 29.8 kilobase with a single-stranded positive-sense RNA (Lu et al., 2020a; Chan et al., 2020; Zhou et al., 2020) encoding 12 putative structural and non-structural proteins; of which spike(S), envelope (E), membrane (M) and nuclecapsid (N) proteins are structural. The S protein complex has two subunits, the S1 subunit contains a single peptide, a receptor-binding domain (RBD) which mediates attachment of virion to host cell surface receptors, and an N-terminal domain (NTD). The S2 subunit mediates fusion between the viral and host cellular membranes which facilitates virus genome entry into the host cell (Gui et al., 2017; Kirchdoerfer et al., 2018; Song et al., 2018; Wan et al., 2020).

Although it has been discovered that the SARS-CoV-2 entry to cells strictly depends on Angiotensin-Converting Enzyme 2 (ACE2) (Li et al., 2003; Li and De Clercq, 2020), it has been shown in several studies that SARS-CoV-2 also needs transmembrane serine protease 2 (TMPRSS2) for S protein priming (Hoffmann et al., 2020a; Hoffmann et al., 2020b; Stopsack et al., 2020). Hoffmann et al. showed that the entrance of SARS-CoV-2 through ACE2 can be blocked using cellular TMPRSS2 blockers. The SARS-CoV-2 spike protein contains several arginine residues with a high cleavability action at the communication point in the S1/S2 cleavage site. It has been suggested that this zoonotic-origin cleavage site sequence is required for SARS-CoV-2 entrance into human cells. These findings are consistent with previous observations from several clinically relevant viruses such as MERS, other Coronaviruses, and Influenza A virus (Kim et al., 2006; Matsuyama et al., 2010; Glowacka et al., 2011; Shulla et al., 2011; Kawase et al., 2012; Gierer et al., 2013; Zhou et al., 2015; Shen et al., 2017; Iwata-Yoshikawa et al., 2019; Kleine-Weber et al., 2019). Moreover, Heurich et al. showed that TMPRSS2 and other potentially related proteases cleave the ACE2 and SARS-S protein leading to the SARS-CoV entry and fusion of the virus S protein with the host cell membrane, respectively (Heurich et al., 2014). Therefore blocking S protein priming by specific serine protease might have the potential to control SARS-CoV-2 infection. There are some therapeutic agents like bromehexine (Lucas et al., 2014), camostat mesylate (Shirato et al., 2013), and nafamostat mesylate (Yamamoto et al., 2016) that have been elucidated to be a good inhibitor of TMPRSS2.

Regarding TMPRSS2 inhibitors as potent candidates for anti SARS-CoV-2 infection, it may be promising to investigate natural resources to discover novel components with anti-TMPRSS2 activities. The extreme and unusual environment of the ocean has extraordinary organisms with astonishing properties, which can reveal new horizons for treatment in modern medicine including marine-derived secondary metabolites with evident anti-inflammatory, antitumor, antimicrobial, antiviral, antimalarial, and antioxidant activities (Molinski et al., 2009; Gogineni et al., 2015; Blunt et al., 2018; Riccio et al., 2020; Yi et al., 2020; Carroll et al., 2021). Based on the extraordinary self-defense capacities of marine organisms and the occurrence in them of some deadly viral infections, these organisms might be regarded as a source of novel antiviral agents, which may be able to combat a SARS-CoV-2 infection. Hence, marine-derived natural compounds should be considered in our efforts to overcome the challenges of COVID-19 treatments.

A promising method to investigate viral entry and proliferation is to apply a computer-aided active site directed inhibition study. To the best of our knowledge, there are no studies that have screened the marine natural products (MNPs) libraries specifically to find a blocker of TMPRSS2. However, several studies have screened the libraries of other natural products, particularly plant-derived compounds, to discover the potential inhibitors of TMPRSS2 (Chikhale et al., 2020; Da Silva Antonio et al., 2020; Idris et al., 2020; Rahman et al., 2020; Singh et al., 2020; Vivek-Ananth et al., 2020; Hu et al., 2021).

This study aimed to find potential inhibitor(s) of TMPRSS2 through virtual screening against a homology model of TMPRSS2 using the library of marine natural products (MNPs).

Materials and Methods

Homology Modeling

The essential step in the study of the structural and functional aspects of any protein is to have its suitable crystal structure. Unfortunately, the three-dimensional (3D) structure of TMPRSS2 had not been found at the time of the current study. Hence, in this case, the only option was to generate a 3D coordinate of TMPRSS2 by comparative prediction approach. Herein, the online server SWISS-MODEL (Guex et al., 2009; Bertoni et al., 2017; Waterhouse et al., 2018) (https://swissmodel.expasy.org/) was used to build the 3D structure of TMPRSS2. The amino acid sequence of human transmembrane protease serine 2, from the Universal Protein Resource “UniProtKB” (accession no: O15393) isoform-2,492 amino acids long (https://www.uniprot.org/uniprot/O15393), was selected for homology modeling that was performed by a template-based method. Then, the RAMPAGE server was used to validate the 3D modeled structure (http://mordred.bioc.cam.ac.uk/∼%7B%7Drapper/rampage.php).

Pharmacophore-Based Virtual Screening

The pharmacophore model was created using the Pharmit server (http://pharmit.csb.pitt.edu/). First, pharmacophore features were automatically extracted from the co-crystalized inhibitor of Serine protease Hepsin (PDB: 5CE1.A), and camostat mesylate was also used as another ligand for the pharmacophore modeling to obtain a more realistic model for homology structure of TMPRSS2 as input receptor. Hence, the pharmacophore model used for TMPRSS2 virtual screening was generated based on the predicted binding interactions of TMPRSS2 with these two inhibitors. The Pharmit parameters for 3D-pharmacophore research were changed according to these pharmacophore parameters. Then the MNP library, which contains 164,952 conformers from 14,064 molecules, was searched on this model. The hit compounds with an RMSD ≥4 Å and minimized affinity ≥ -6 were discarded. The remaining poses were minimized using functions of Pharmit. Finally, the pharmacophore-based minimized entries were further docked with the target protein to identify the lead compounds with the best docking scores.

Active Site Identification and Preparation of TMPRSS2 for Docking

The conserved domain (CD) search was done on the Fasta sequence of the homology model of TMPRSS2 using the NCBI’s conserved domain database (CDD/SPARCLE: https://www.ncbi.nlm.nih.gov/cdd/) (Lu et al., 2020b). Then, it was analyzed and its cleavage, active, and substrate binding site residues were predicted. The COACH-D server (https://yanglab.nankai.edu.cn/COACH-D/) (Yang et al., 2013; Wu et al., 2018) was also applied to predict TMPRSS2 putative ligand-binding sites. Then, in a model-template active site comparative study, the active and substrate binding residues of the (5CE1.A) template were similarly obtained.

Finally, the homology structure of TMPRSS2 was imported into Maestro Protein Wizard and sitemap analysis was done by SiteMap (Halgren, 2007; Halgren, 2009; Schrödinger, 2015b).

This software was used with default settings, in which the top five possible binding sites by a minimum of 15 points were identified while cropping site maps set at 4 A from the nearest site point. SiteMap used a more restrictive definition of hydrophobicity by standard grid. A SiteScore value above 0.80 is indicative of high druggability and promising drug-binding sites and is used in conjunction with Dscores, which serve as a measure of hydrophobicity. Dscore or druggability score penalizes increasing hydrophilicity and is thus used as a druggability measure for a pocket. In general, Dscore <0.83 is considered as “undruggable,” 0.83–0.98 as “difficult to drug” and >0.98 as “druggable” (Halgren, 2007; Halgren, 2009; Vidler et al., 2012).

Molecular Docking

Ligand Preparation

One conformation was generated per compound which was followed by geometry optimization with PM3 (Stewart, 1991), a semi-empirical method using Hyperchem release 7 for windows (HyperCube Inc., 2002). The geometry-optimized structures were retrieved in MOL type for further analysis with the LigPrep application which has been implemented in the Schrödinger 2015-2 suite of software (Schrödinger, 2015a). The ionization state was specified at pH = 7.00±2.0 using Epik (Epik, 2015) based on Hammett and Taft methodologies (Sastry et al., 2013). The desalt option was the same as the program default. All 32 possible conformations were produced for each compound at pH 7.00 in the OPLS3 (Harder et al., 2016). The obtained ligands were then used in the docking calculations.

Generation of the Grid

As adequate generation of the grid is a key step in the prediction of a ligand binding to a receptor, the 3D boundary for ligand binding was produced by Glide, version 10.2 of Mastero, Schrödinger (Glide, 2015). First, the protein preparation wizard was used with the following settings: 1) The original hydrogens removal and subsequent addition of hydrogens. 2) The atomic charges and bond orders were assigned. 3) The N and C termini were capped. 4) The disulfide bonds were generated between sulfur atoms (within 3.2 Å). 5) Epik was applied to generate possible protonation states at neutral pH. 6) The H-bonds were assigned, optimized by PROPKA (Olsson et al., 2011; Søndergaard et al., 2011) at pH 7.0, and then the structure was minimized with the OPLS3 force field. Glide was used to generate the grid on the catalytic domain of the receptor. The grid box size was set to 32*32*32 Å.

Molecular Interaction and Docking Studies

After grid generation, ligand docking was done according to the protocols in Glide version 10.2. The homology model of TMPRSS2 was used as the receptor, and the different internally produced conformations by the software were passed across some filters such as Euler angles, grid-based force field evaluation, and energy minimization by Monte Carlo. Finally, docking score is an important parameter for evaluating the conformations, and in this study, the output of standard precision (SP) docking was put forward in extra precision (XP) docking. The docked compounds were ranked based on their docking scores.

Pose Rescoring With Molecular Mechanics Generalized Born Surface Area

The binding energies of all docking poses were calculated using the molecular mechanics generalized Born surface area (MM-GBSA) approach implemented in the Prime program in the Schrödinger software suite (Prime, 2015). This approach employs a single minimized protein-ligand structure, and so is used as an efficient approach to rapidly refine and rescore docking results. A variable dielectric solvent model VSGB 2.0 (Li et al., 2011) was used, this solvent model contains several empirical corrections for modeling the directionality of hydrogen bond and π-stacking interactions. MM-GBSA has been shown to give good binding free energies for a wide range of protein-ligand complexes (Mulakala and Viswanadhan, 2013). It is also widely used to evaluate docking poses, to determine the stability of ligand-target complex for predicting binding affinity in drug design (Genheden and Ryde, 2011; Wang et al., 2019).

Absorption, Distribution, Metabolism, and Excretion and Drug-likeness Analysis

It is known that nearly 40% of drug candidates fail in clinical trials because of poor absorption, distribution, metabolism, and excretion (ADME). Hence, it is very crucial to recognize these problematic candidates at an early stage to avoid wasted time and resources. Accurate ADME prediction is based on full 3D molecular structures. Qikprop offers a set of several predictors including central nervous system (CNS) penetration, predicted apparent Caco-2 cell permeability across the gut-blood barrier in nm/sec (QPPCaco), apparent MDCK cell permeability (QPPMDCK), human oral absorption, Lipinski’s rule of five, and predicted maximum transdermal transport rate (JM). Another option of Qikprop (Qikprop, 2015) is to rank compounds based on how drug-like they are. In the current study, the ADME, drug-likeness, and medicinal chemistry parameters of these 11 compounds were predicted by QikProp (Qikprop, 2015).

Results and Discussion

Homology Modeling and Evaluation of Model Quality

Since the crystal structure of TMPRSS2 was unavailable at the time of this study, the 3D structure of TMPRSS2 was predicted using the online server SWISS-MODEL, as shown in Figure 1A. Its global model quality estimate (GMQE) score was 0.48, this score estimates the quality of the expected output model by a particular template. Its QMEAN Z-score was -1.47, with sequence coverage of 71%, sequence identity of 33.82%, and sequence similarity of 50% in comparison with the template (PDB ID: 5CE1.A). It has been shown that when the sequence similarity with the template is more than 30%, the obtained model can be considered reliable and suitable for further study (Xiang, 2006). Benkert et al. (2011) showed that QMEAN Z-score is an estimation of the degree of native likeness of the model, and a value close to 0 (and not lower than -4) could be an acceptable agreement criteria for the experimental structure of similar size Table 1.

FIGURE 1
www.frontiersin.org

FIGURE 1. (A) Three-dimensional structure of the modeled serine protease transmembrane protease serine 2 (TMPRSS2), (B) Ramachandran plot validation of the modeled 3D structure, (C) alignment of the target serine protease TMPRSS2 and the template serine protease hepsin (PDB ID: 5CE1.A.).

TABLE 1
www.frontiersin.org

TABLE 1. Homology modeling results and validation for predicting 3D structure of TMPRSS2.

By validating the obtained results of SWISS-MODEL and cross-checking in RAMPAGE, it was observed that there were 319 (92.7%) residues were in the favored region, 23 (6.7%) residues in the allowed region, and 2 (0.6%) residues in the outlier region. Moreover, there were no steric clashes or deviations in bond length or bond angle compared to the protein structure report (Supplementary Figure S1 and Supplementary Figure S2). These results indicated that the obtained model may have the correct geometry. The 3D arrangement of the model is shown in Figure 1B. The alignment of the template (PDB ID: 5CE1.A) and the target protein is shown in Figure 1C. The summary of obtained results is presented in Table 1.

Catalytic Site of the Homology Model of TMPRSS2

According to a conserved domain database (CDD) search, the protein classification of TMPRSS2 was Trypsin-like serine proteases with E-value 7.39e-100, and the residues of ILE293-GLN524 were involved in the characteristic domain of this protein. Based on the CDD algorithm, six amino acid residues have been identified as especially important in the active site of TMPRSS2, and one residue for its cleavage site. The active site included residue HIS333 to SER478 where (HIS333, ASP382, and SER478) were the three important amino acids at the catalytic site whereas (ASP472, SER497, and GLY499) residues were found to be the substrate binding site (Table 2). Moreover, the COACH-D results supported those findings. According to the COACH-D best prediction results, the protein template was a serine protease of the coagulation system (PDB ID: 5JB8) with a confidence score of 0.99.

TABLE 2
www.frontiersin.org

TABLE 2. Predicted binding site of homology model of TMPRSS2.

The corresponding docking energy of the template for its representative ligand (with pubchem CID: 137347860) was -6.0 kcal/mol. HIS333, LYS379, ASP472-GLY476, SER478, THR496-CYS502, and GLY509 were identified as the predicted binding residues, of which the CDD search showed HIS333, ASP472, SER478, SER497, and GLY499, as the same catalytic domain residues.

From the five binding sites which were calculated by sitemap, two sites had a Dscore above 1, and their locations were completely in the agreement with the CDD search and COACH-D results (Site 1: Sitescore 0.972, Dscore 1.001, size 87, and volume 150. Site 2: Sitescore 0.968, Dscore 1.026, size 91, and volume 228). The Dscore of site 3 was 1.044 with Sitescore 0.994 but its position is not in agreement with the catalytic site of Trypsin-like serine proteases active domain which is located in the base of its S1 pocket where it contains Asp472, and therefore is predicted to cleave after lysine or arginine residues (Hedstrom, 2002; Wilson et al., 2005; Blay and Pei, 2019). The predicted binding sites by all applied approaches are in agreement with the theoretical study by Idris et al. (2020). Two other predicted binding sites had a Dscore lower than 0.8, which indicates that they are not druggable. Hence, domains 1 and 2 were chosen for the docking study. Figures 2A,B illustrate the predicted binding site of TMPRSS2 from COACH-D and SiteMap analysis.

FIGURE 2
www.frontiersin.org

FIGURE 2. (A) The predicted binding site by CPACH-D server. The protein template is a serine protease of the coagulation system (PDB: 5jb8) and ligand is [(4S, 5S)-4-[[2-[[(2S)-2-amino-4-carboxybutanoyl]amino]acetyl]amino]-6-chloro-5-hydroxyhexyl]-(diaminomethylidene)azanium with pubchem CID: 137347860. (B) The binding sites predicted by Sitemap. Site1: Sitescore 0.972, Dscore 1.001, size 87, and volume 150. Site 2: Sitescore 0.968, Dscore 1.026, size 91, and volume 228. Site 3: Sitescore 0.994, Dscore 1.044, size 108, and volume 331.

Pharmacophore Model

Since a pharmacophore states the crucial features of interactions, such as the spatial arrangement of each interaction in the close contact of ligand and the target, its accurate setup is very important in binding site pharmacophore modeling. In this study, the 3D structure of TMPRSS2 homology was used to set pharmacophore using the Pharmit server; this server provides both pharmacophore and molecular shape search options and the results are ranked by the energy. The generated pharmacophore features were selected according to the co-crystal inhibitor of the template and camostat mesylate. In this modeling, the binding-site derived pharmacophore models include three subgroups of ligand binding sites: i) two amide nitrogen atoms were added to represent hydrogen bond donors (DON), ii) four negatively charged oxygen atoms (as in a carboxyl group) were added to represent a hydrogen bond acceptor (ACC), and iii) the two isopropyl group were added to represent a hydrophobic center (HYD) (Figure 3A). According to the generated pharmacophore model, a vast library of MNP (14,064 molecules, 164,952 conformers) was filtered. A total of 25,000 hits that met the criteria were minimized, resulting in 114 conformers. A total of 11 structures were retained by using one conformer for each molecule, with an RMSD lower than 4 Å and a binding score lower than -6 (Figure 3B and Table 3).

FIGURE 3
www.frontiersin.org

FIGURE 3. (A) Pharmacophore model generated by the Pharmit server, Two amide nitrogen atoms to represent hydrogen bond donors (DON) (green sphere), four negatively charged oxygen atoms (as in a carboxyl group) to represent a hydrogen bond acceptor (ACC) (orange sphere), and the two isopropyl group to represent a hydrophobic center (HYD) (yellow sphere) (B) superposition of all 11 aligned lead MNPs according to pharmacophore model.

TABLE 3
www.frontiersin.org

TABLE 3. CAS number, 2D and 3D presentations of investigated marine natural products.

Molecular Docking and Molecular Mechanics Generalized Born Surface Area Studies

The top 11 selected MNPs of the 114 structures which are shown in Table 3 were separately docked into the catalytic site of TMPRSS2. The results of docking of these structures, as well as camostat and nafamostat mesylate, as two standard inhibitors of TMPRSS2 (Scheme 1), are presented in Table 4.

SCHEME 1
www.frontiersin.org

SCHEME 1. Standard inhibitors (A) camostat mesylate and (B) nafamostat mesylate.

TABLE 4
www.frontiersin.org

TABLE 4. Glide Docking score, Glide energy, Glide emodel, and estimated free energy of binding for the best poses of investigated compounds in kcal/mol.

The molecular docking analysis revealed that all the studied compounds had comparable or lower docking scores than those of the standard inhibitors. The highest docking Glide score is −8.16 in compound 10, whereas these scores were −4.52 and −3.73 in camostat and nafamostat mesylate respectively. Also, the docking scores in compounds 3, 4, 6, 7, and 11 are lower than −7. Thus, these MNPs could be considered as the most potent inhibitors for TMPRSS2. Moreover, the highest Glide energy value was −59.21 for compound 10, and its Glide emodel value is −71.34, which are almost the highest values. The Glide energies were −42.47 and −38.27 for camostat mesylate and nafamostat mesylate respectively, and Glide emodel values were -53.19 and −45.75 respectively. Therefore, from the observed theoretical superiority of these 11 MNPs compared to the standard inhibitors (camostat and nafamostat mesylate), these compounds may be encouraging for further studies. Interestingly, among these selected compounds, compound 10 had the most promising results. Compound 10 (CAS number 107503-09-3; Watasenia Preluciferyl β-d-glucopyrasoiuronic acid) is a bioluminescent substance that was derived from the liver of myctophina fish, Diaphus elucens (Inoue et al., 1987; Blunt and Munro, 2007).

As it is presented in Table 5, the 2D template of best poses demonstrates the types of contacts formed between the ligands and target with cutoff 4.00 A. Remarkably, close contact/interactions within the catalytic domain were detected for all MNPs, however, all the important residues of catalytic domains significantly contributed in interactions with compound 10, such as HIS333 which was involved in π-π stacking interaction, and SER497 and GLY499 which were the important residues of substrate binding, in the close vicinity of phenol and iduronic acid moieties (Figure 4A). In addition, SER473, GLU426, LYS379, and GLU336 participated in intermolecular H-bond interaction with the active domain (Figure 4B). Taken together, compounds 3, 4, 10, and 11 had stronger interactions with TMPRSS2 than the standard ligands based on both their Glide and XP-pose emodel energies and were significantly involved in active domain contacts/interactions. On the other hand, the energy values of compounds 2 and 9 were also considerable and close to the references. Molecule 3, Downeyoside I, firstly isolated from starfish Henricia Downyae by Plagiano; molecule 2, Forbeside E, a sulfated sterol glycoside from starfish Asterias forbesi; molecule 4, Ulososide E, which has been derived from the sponge Ulosa sp; and molecule 9, Downeyoside E, sulfated steroid glycoside isolated from Henricia Downeyae (Blunt and Munro, 2007).

TABLE 5
www.frontiersin.org

TABLE 5. 2D presentation of best Glide docking pose at the catalytic domain of TMPRSS2 for investigated marine natural product.

FIGURE 4
www.frontiersin.org

FIGURE 4. (A) The best pose of docking (2D) of ligand 10 at the predicted catalytic domain of TMPRSS2. (B) The 3D presentation of the best docking pose of ligand 10.

In addition, the MM-GBSA calculations which estimated the values of ΔGbind are reported in Table 4. According to these results, the two standard inhibitors had values of −57.91 and −47.49 kcal/mol while eight molecules out of the 11 studied MNPs had more negative  ΔGbind than the standard inhibitors. The highest binding free energies were for compounds 8 and 10 (−76.00 kcal/mol) which is much higher than the standard inhibitors, and the values of compounds 2, 3, 4, 6, 9, and 11 were more negative than the standard inhibitors. After the rescoring of the 11 MNPs, compound 10 remained as a potent inhibitor of TMPRSS2, however other ligands might be also interesting for further studies. Although these compounds belong to various sources, starfish-derived products sound to be prominent. For example, 2, 3, and 9 have been isolated from starfish, whereas 7 has been isolated from the marine sponge Merriamum oxeato (Blunt and Munro, 2007). In this study, the marine natural products showed good measurable binding affinities for the TMPRSS2 residues. In other words, these binding affinities are indicative of the ligand’s contribution to ligand-target interactions and their sensible flexibility for this target. Based on XP Glide docking score, compounds 10 and 11 have strong interactions with the enzyme but rescoring by MM-GBSA suggests that other compounds such as 2, 3, 4, 8, and 9 are also important. Nevertheless, our proposed lead compound is compound 10, which has concurrently both high docking scores and a comparable  ΔGbind to standard inhibitors.

Absorption, Distribution, Metabolism, and Excretion and Drug-likeness Analysis

In this study, the ADME https://www.sciencedirect.com/topics/pharmacology-toxicology-and-pharmaceutical-science/admeproperties of 11 MNPs were analyzed using the QikProp tool. This analysis represents the physicochemical properties of chemical compounds along with their biological functions. The resulting physicochemical and biological properties are molecular formula, molecular weight, volume, SASA, acceptor H-bond, donor H-bond groups, the number of ring atoms, QPlogPw (−2–6.5), the percentage of human oral absorption, and CNS effects. However, it has been suggested that Lipinski’s rule of five (Zhang and Wilkinson, 2007) is not a strict criterion for natural compounds (Lipinski, 2003) and it has been revealed that natural compounds mostly do not follow Lipinski’s rule and they tend to keep their low hydrophobicity as well as their potential of donating the intermolecular H-bonds (Ganesan, 2008).

In general, Lipinski’s rule of five (Lipinski et al., 1997) is applied for predicting the drug-likeness with the following criteria: molecular mass less than 500 Da, up to 5 hydrogen bond donors, no more than 10 hydrogen bond acceptors, and an octanol-water partition coefficient (logPo/w) no higher than 5. The rule states that a molecule or an inhibitor can be orally absorbed/active if two or more of these thresholds are not violated. However, Jorgensen’s rule of three may also be used to evaluate the bioavailability of each marine natural product by estimating its solubility, permeability, and liver first-pass metabolism through the following rules: predicted aqueous solubility (logSwat) higher than -5.7 (with S in mol/dm3), predicted apparent Caco-2 cell rate permeability (BIPcaco-2) high than 22 nm/s, and number of primary metabolites up to 7 (Di and Kerns, 2015). In addition, the predicted qualitative human oral absorption (2 = medium and 3 = high) and the predicted skin permeability (logKp values between -8.0–1.0) are considered. Finally, the ADME-compliance score drug-likeness parameter (#star) was used to evaluate the pharmacokinetic of the studied compounds, including 25 different properties within the acceptable range of 95% of the known drugs. Herein, compounds 2, 5, 9, and 10 had the fewest violations when ADME-compliance score drug-likeness parameter (#stars) was considered less than 2 (Table 6), however, the recommended value is 0–5, and compounds 6 and 8 were also considered in the acceptable range of #stars. It is concluded that these MNPs may be proper candidate drugs for TMPRSS2 inhibition.

TABLE 6
www.frontiersin.org

TABLE 6. ADME parameters and drug similarity of investigated marine natural products.

Again, amongst the latter compounds, compound 10 showed a higher human oral absorption (23), as well as median aggregation to plasma proteins (Qlog k has serum protein binding: −0.13) and predicted aqueous solubility values (QPlog S: 3.22). Moreover, by considering the number of “stars” and the violations from the Lipsinki and Jorgensen rules, the obtained results for compound 10 indicated a high degree of reliability to be a drug candidate.

To find chemical similarity to the known drug molecules, the QikProp (Qikprop, 2015) software database identified five similar drug molecules for each entry according to its predicted descriptors. In this study, the results of two of the compounds were very promising. Accordingly, Azithromycin with a similarity of 65.70% was suggested for compound 6. It is interesting to know that in the early phase of COVID-19, azithromycin could reduce the need for hospitalization or duration of clinical recovery (Echeverría-Esnal et al., 2020; Million et al., 2020; Molina et al., 2020). There is also an opinion supporting the potential effectiveness of Azithromycin in SARS-CoV-2 infection, as well as its antiviral activity and immunomodulatory effects (Bleyzac et al., 2020). Interestingly, Amprenavir (Shen et al., 2010), as an antiretroviral protease inhibitor for HIV infection has been identified with 68.5% similarity to compound 10.

In conclusion, marine natural product 10, with high similarity to a known antiretroviral protease, might be considered as a potent inhibitor of TMPRSS2. However, molecular docking and MM-GBSA studies demonstrated that compounds 3, 4, 6, 7, 8, 9, and 11 were also able to inhibit TMPRSS2 as well. In a comparison study, these compounds showed better results than the standard TMPRSS2 inhibitors (camostat and nafamostat mesylate). Further computational, experimental, and clinical investigations are warranted to reveal their anti-SARS-CoV-2 activities.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.

Author Contributions

MF conceived and designed the format of the manuscript. All calculations have been done by MF. All the authors, MM, IN, and MK reviewed the manuscript and all of them contributed to the critical reading and discussion of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The authors wish to express their thanks for the partially financial support of Bushehr University of medical sciences.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2021.722633/full#supplementary-material

References

Astuti, I., and Ysrafil, Y. (2020). Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2): An Overview of Viral Structure and Host Response. Diabetes Metab. Syndr. Clin. Res. Rev. 14, 407–412. doi:10.1016/j.dsx.2020.04.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Benkert, P., Biasini, M., and Schwede, T. (2011). Toward the Estimation of the Absolute Quality of Individual Protein Structure Models. Bioinformatics 27, 343–350. doi:10.1093/bioinformatics/btq662

PubMed Abstract | CrossRef Full Text | Google Scholar

Bertoni, M., Kiefer, F., Biasini, M., Bordoli, L., and Schwede, T. (2017). Modeling Protein Quaternary Structure of Homo- and Hetero-Oligomers beyond Binary Interactions by Homology. Sci. Rep. 7, 10480–10495. doi:10.1038/s41598-017-09654-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Blay, V., and Pei, D. (2019). Serine Proteases: How Did Chemists Tease Out Their Catalytic Mechanism. ChemTexts 5, 1–7. doi:10.1007/s40828-019-0093-4

CrossRef Full Text | Google Scholar

Bleyzac, N., Goutelle, S., Bourguignon, L., and Tod, M. (2020). Azithromycin for COVID-19: More Than Just an Antimicrobial. Clin. Drug Investig. 40 (8), 683–686. doi:10.1007/s40261-020-00933-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Blunt, J. W., Carroll, A. R., Copp, B. R., Davis, R. A., Keyzers, R. A., and Prinsep, M. R. (2018). Marine Natural Products. Nat. Prod. Rep. 35, 8–53. doi:10.1039/c7np00052a

PubMed Abstract | CrossRef Full Text | Google Scholar

Blunt, J. W., and Munro, M. H. (2007). Dictionary of marine Natural Products with CD-ROM. London, United Kingdom: Chapman and Hall/CRC Press.

Google Scholar

Carroll, A. R., Copp, B. R., Davis, R. A., Keyzers, R. A., and Prinsep, M. R. (2021). Marine Natural Products. Nat. prod. Rep. 38 (2), 362–413. doi:10.1039/D0NP00089B

PubMed Abstract | CrossRef Full Text | Google Scholar

Chan, J. F.-W., Kok, K.-H., Zhu, Z., Chu, H., To, K. K.-W., Yuan, S., et al. (2020). Genomic Characterization of the 2019 Novel Human-Pathogenic Coronavirus Isolated from a Patient with Atypical Pneumonia after Visiting Wuhan. Emerging Microbes & Infections 9, 221–236. doi:10.1080/22221751.2020.1719902

CrossRef Full Text | Google Scholar

Chikhale, R. V., Gupta, V. K., Eldesoky, G. E., Wabaidur, S. M., Patil, S. A., and Islam, M. A. (2020). Identification of Potential Anti-TMPRSS2 Natural Products through Homology Modelling, Virtual Screening and Molecular Dynamics Simulation Studies. J. Biomol. Struct. Dyn. 39, 6660–6675. doi:10.1080/07391102.2020.1798813

CrossRef Full Text | Google Scholar

Da Silva Antonio, A., Wiedemann, L. S. M., and Veiga-Junior, V. F. (2020). Natural Products' Role against COVID-19. RSC Adv. 10, 23379–23393. doi:10.1039/d0ra03774e

CrossRef Full Text | Google Scholar

Di, L., and Kerns, E. H. (2015). Drug-Like Properties: Concepts, Structure Design and Methods from ADME to Toxicity Optimization. Cambridge: Academic Press.

Google Scholar

Echeverría-Esnal, D., Martin-Ontiyuelo, C., Navarrete-Rouco, M. E., De-Antonio Cuscó, M., Ferrández, O., Horcajada, J. P., et al. (2020). Azithromycin in the Treatment of COVID-19: A Review. Expert Rev. Anti Infect. Ther. 19 (2), 147–163. doi:10.1080/14787210.2020.1813024

PubMed Abstract | CrossRef Full Text | Google Scholar

Epik (2015). version 3. 2 ed. New York, NY: Schrödinger.

Ganesan, A. (2008). The Impact of Natural Products upon Modern Drug Discovery. Curr. Opin. Chem. Biol. 12, 306–317. doi:10.1016/j.cbpa.2008.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Genheden, S., and Ryde, U. (2011). Comparison of the Efficiency of the LIE and MM/GBSA Methods to Calculate Ligand-Binding Energies. J. Chem. Theor. Comput. 7, 3768–3778. doi:10.1021/ct200163c

PubMed Abstract | CrossRef Full Text | Google Scholar

Gierer, S., Bertram, S., Kaup, F., Wrensch, F., Heurich, A., Krämer-Kühl, A., et al. (2013). The Spike Protein of the Emerging Betacoronavirus EMC Uses a Novel Coronavirus Receptor for Entry, Can Be Activated by TMPRSS2, and Is Targeted by Neutralizing Antibodies. J. Virol. 87, 5502–5511. doi:10.1128/jvi.00128-13

CrossRef Full Text | Google Scholar

Glide (2015). version 6. 7 ed. New York, NY: Schrödinger.

Glowacka, I., Bertram, S., Müller, M. A., Allen, P., Soilleux, E., Pfefferle, S., et al. (2011). Evidence that TMPRSS2 Activates the Severe Acute Respiratory Syndrome Coronavirus Spike Protein for Membrane Fusion and Reduces Viral Control by the Humoral Immune Response. J. Virol. 85, 4122–4134. doi:10.1128/jvi.02232-10

CrossRef Full Text | Google Scholar

Gogineni, V., Schinazi, R. F., and Hamann, M. T. (2015). Role of Marine Natural Products in the Genesis of Antiviral Agents. Chem. Rev. 115, 9655–9706. doi:10.1021/cr4006318

PubMed Abstract | CrossRef Full Text | Google Scholar

Guan, W.-J., Ni, Z.-Y., Hu, Y., Liang, W.-H., Ou, C.-Q., He, J.-X., et al. (2020). Clinical Characteristics of Coronavirus Disease 2019 in China. N. Engl. J. Med. 382, 1708–1720. doi:10.1056/nejmoa2002032

PubMed Abstract | CrossRef Full Text | Google Scholar

Guex, N., Peitsch, M. C., and Schwede, T. (2009). Automated Comparative Protein Structure Modeling with SWISS-MODEL and Swiss-PdbViewer: A Historical Perspective. Electrophoresis 30, S162–S173. doi:10.1002/elps.200900140

PubMed Abstract | CrossRef Full Text | Google Scholar

Gui, M., Song, W., Zhou, H., Xu, J., Chen, S., Xiang, Y., et al. (2017). Cryo-electron Microscopy Structures of the SARS-CoV Spike Glycoprotein Reveal a Prerequisite Conformational State for Receptor Binding. Cell Res. 27, 119–129. doi:10.1038/cr.2016.152

PubMed Abstract | CrossRef Full Text | Google Scholar

Halgren, T. A. (2009). Identifying and Characterizing Binding Sites and Assessing Druggability. J. Chem. Inf. Model. 49, 377–389. doi:10.1021/ci800324m

CrossRef Full Text | Google Scholar

Halgren, T. (2007). New Method for Fast and Accurate Binding-Site Identification and Analysis. Chem. Biol. Drug Des. 69, 146–148. doi:10.1111/j.1747-0285.2007.00483.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Harder, E., Damm, W., Maple, J., Wu, C., Reboul, M., Xiang, J. Y., et al. (2016). OPLS3: A Force Field Providing Broad Coverage of Drug-Like Small Molecules and Proteins. J. Chem. Theor. Comput. 12, 281–296. doi:10.1021/acs.jctc.5b00864

PubMed Abstract | CrossRef Full Text | Google Scholar

Hedstrom, L. (2002). Serine Protease Mechanism and Specificity. Chem. Rev. 102, 4501–4524. doi:10.1021/cr000033x

PubMed Abstract | CrossRef Full Text | Google Scholar

Heurich, A., Hofmann-Winkler, H., Gierer, S., Liepold, T., Jahn, O., and Pöhlmann, S. (2014). TMPRSS2 and ADAM17 Cleave ACE2 Differentially and Only Proteolysis by TMPRSS2 Augments Entry Driven by the Severe Acute Respiratory Syndrome Coronavirus Spike Protein. J. Virol. 88, 1293–1307. doi:10.1128/jvi.02202-13

CrossRef Full Text | Google Scholar

Hoffmann, M., Kleine-Weber, H., Schroeder, S., Krüger, N., Herrler, T., Erichsen, S., et al. (2020a). SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. cell 181, 271–280. e278. doi:10.1016/j.cell.2020.02.052

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoffmann, M., Schroeder, S., Kleine-Weber, H., Müller, M. A., Drosten, C., and Pöhlmann, S. (2020b). Nafamostat Mesylate Blocks Activation of SARS-CoV-2: New Treatment Option for COVID-19. Antimicrob. Agents Chemother. 64 (6), e00754–20. doi:10.1128/AAC.00754-20

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, X., Shrimp, J. H., Guo, H., Xu, M., Chen, C. Z., Zhu, W., et al. (2021). Discovery of TMPRSS2 Inhibitors from Virtual Screening as a Potential Treatment of COVID-19. ACS Pharmacol. translational Sci. 4 (3), 1124–1135. doi:10.1021/acsptsci.0c00221

PubMed Abstract | CrossRef Full Text | Google Scholar

Hypercube Inc. (2002). Hyperchem 7.0. Toronto, C.: Hypercube Inc.

Idris, M. O., Yekeen, A. A., Alakanse, O. S., and Durojaye, O. A. (2020). Computer-Aided Screening for Potential TMPRSS2 Inhibitors: a Combination of Pharmacophore Modeling, Molecular Docking and Molecular Dynamics Simulation Approaches. J. Biomol. Struct. Dyn. 39 (15), 5638–5656. doi:10.1080/07391102.2020.1792346

CrossRef Full Text | Google Scholar

Inoue, S., Okada, K., Tanino, H., and Kakoi, H. (1987). Chemical Studies on Myctophina Fish Bioluminescence. Chem. Lett. 16, 417–418. doi:10.1246/cl.1987.417

CrossRef Full Text | Google Scholar

Iwata-Yoshikawa, N., Okamura, T., Shimizu, Y., Hasegawa, H., Takeda, M., and Nagata, N. (2019). TMPRSS2 Contributes to Virus Spread and Immunopathology in the Airways of Murine Models after Coronavirus Infection. J. Virol. 93, e01815–18. doi:10.1128/JVI.01815-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawase, M., Shirato, K., Van Der Hoek, L., Taguchi, F., and Matsuyama, S. (2012). Simultaneous Treatment of Human Bronchial Epithelial Cells with Serine and Cysteine Protease Inhibitors Prevents Severe Acute Respiratory Syndrome Coronavirus Entry. J. Virol. 86, 6537–6545. doi:10.1128/jvi.00094-12

CrossRef Full Text | Google Scholar

Kim, T. S., Heinlein, C., Hackman, R. C., and Nelson, P. S. (2006). Phenotypic Analysis of Mice Lacking the Tmprss2 -Encoded Protease. Mol. Cel. Biol. 26, 965–975. doi:10.1128/mcb.26.3.965-975.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Kirchdoerfer, R. N., Wang, N., Pallesen, J., Wrapp, D., Turner, H. L., Cottrell, C. A., et al. (2018). Stabilized Coronavirus Spikes Are Resistant to Conformational Changes Induced by Receptor Recognition or Proteolysis. Sci. Rep. 8, 15701–15711. doi:10.1038/s41598-018-34171-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Kleine-Weber, H., Elzayat, M. T., Wang, L., Graham, B. S., Müller, M. A., Drosten, C., et al. (2019). Mutations in the Spike Protein of Middle East Respiratory Syndrome Coronavirus Transmitted in Korea Increase Resistance to Antibody-Mediated Neutralization. J. Virol. 93, e01381–18. doi:10.1128/JVI.01381-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, G., and De Clercq, E. (2020). Therapeutic Options for the 2019 Novel Coronavirus (2019-nCoV). Nat. Rev. Drug Discov. 19, 149–150. doi:10.1038/d41573-020-00016-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Abel, R., Zhu, K., Cao, Y., Zhao, S., and Friesner, R. A. (2011). The VSGB 2.0 Model: A Next Generation Energy Model for High Resolution Protein Structure Modeling. Proteins 79, 2794–2812. doi:10.1002/prot.23106

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, W., Moore, M. J., Vasilieva, N., Sui, J., Wong, S. K., Berne, M. A., et al. (2003). Angiotensin-Converting Enzyme 2 Is a Functional Receptor for the SARS Coronavirus. Nature 426, 450–454. doi:10.1038/nature02145

PubMed Abstract | CrossRef Full Text | Google Scholar

Lipinski, C. A. (2003). Chris Lipinski Discusses Life and Chemistry after the Rule of Five. Drug Discov. Today 8, 12–16. doi:10.1016/s1359-6446(02)02556-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Lipinski, C. A., Lombardo, F., Dominy, B. W., and Feeney, P. J. (1997). Experimental and Computational Approaches to Estimate Solubility and Permeability in Drug Discovery and Development Settings. Adv. Drug Deliv. Rev. 23, 3–25. doi:10.1016/s0169-409x(96)00423-1

CrossRef Full Text | Google Scholar

Lu, R., Zhao, X., Li, J., Niu, P., Yang, B., Wu, H., et al. (2020a). Genomic Characterisation and Epidemiology of 2019 Novel Coronavirus: Implications for Virus Origins and Receptor Binding. The lancet 395, 565–574. doi:10.1016/s0140-6736(20)30251-8

CrossRef Full Text | Google Scholar

Lu, S., Wang, J., Chitsaz, F., Derbyshire, M. K., Geer, R. C., Gonzales, N. R., et al. (2020b). CDD/SPARCLE: The Conserved Domain Database in 2020. Nucleic Acids Res. 48, D265–D268. doi:10.1093/nar/gkz991

PubMed Abstract | CrossRef Full Text | Google Scholar

Lucas, J. M., Heinlein, C., Kim, T., Hernandez, S. A., Malik, M. S., True, L. D., et al. (2014). The Androgen-Regulated Protease TMPRSS2 Activates a Proteolytic cascade Involving Components of the Tumor Microenvironment and Promotes Prostate Cancer Metastasis. Cancer Discov. 4, 1310–1325. doi:10.1158/2159-8290.cd-13-1010

PubMed Abstract | CrossRef Full Text | Google Scholar

Madhavi Sastry, G., Adzhigirey, M., Day, T., Annabhimoju, R., and Sherman, W. (2013). Protein and Ligand Preparation: Parameters, Protocols, and Influence on Virtual Screening Enrichments. J. Comput. Aided Mol. Des. 27, 221–234. doi:10.1007/s10822-013-9644-8

CrossRef Full Text | Google Scholar

Matsuyama, S., Nagata, N., Shirato, K., Kawase, M., Takeda, M., and Taguchi, F. (2010). Efficient Activation of the Severe Acute Respiratory Syndrome Coronavirus Spike Protein by the Transmembrane Protease TMPRSS2. J. Virol. 84, 12658–12664. doi:10.1128/jvi.01542-10

CrossRef Full Text | Google Scholar

Million, M., Lagier, J.-C., Gautret, P., Colson, P., Fournier, P.-E., Amrane, S., et al. (2020). Early Treatment of COVID-19 Patients with Hydroxychloroquine and Azithromycin: A Retrospective Analysis of 1061 Cases in Marseille, France. Trav. Med. Infect. Dis. 35, 101738. doi:10.1016/j.tmaid.2020.101738

PubMed Abstract | CrossRef Full Text | Google Scholar

Molina, J. M., Delaugerre, C., Le Goff, J., Mela-Lima, B., Ponscarme, D., Goldwirt, L., et al. (2020). No Evidence of Rapid Antiviral Clearance or Clinical Benefit with the Combination of Hydroxychloroquine and Azithromycin in Patients with Severe COVID-19 Infection. Médecine et Maladies Infectieuses 50, 384. doi:10.1016/j.medmal.2020.03.006

CrossRef Full Text | Google Scholar

Molinski, T. F., Dalisay, D. S., Lievens, S. L., and Saludes, J. P. (2009). Drug Development from Marine Natural Products. Nat. Rev. Drug Discov. 8, 69–85. doi:10.1038/nrd2487

PubMed Abstract | CrossRef Full Text | Google Scholar

Mulakala, C., and Viswanadhan, V. N. (2013). Could MM-GBSA Be Accurate Enough for Calculation of Absolute Protein/Ligand Binding Free Energies. J. Mol. Graphics Model. 46, 41–51. doi:10.1016/j.jmgm.2013.09.005

CrossRef Full Text | Google Scholar

Olsson, M. H. M., Søndergaard, C. R., Rostkowski, M., and Jensen, J. H. (2011). PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. J. Chem. Theor. Comput. 7, 525–537. doi:10.1021/ct100578z

PubMed Abstract | CrossRef Full Text | Google Scholar

Prime (2015). version 4.0. New York, NY: Schrödinger, LLC.

Google Scholar

Qikprop (2015). version 4.4. New York, NY: Schrödinger, LLC.

Google Scholar

Rahman, N., Basharat, Z., Yousuf, M., Castaldo, G., Rastrelli, L., and Khan, H. (2020). Virtual Screening of Natural Products against Type II Transmembrane Serine Protease (TMPRSS2), the Priming Agent of Coronavirus 2 (SARS-CoV-2). Molecules 25, 2271. doi:10.3390/molecules25102271

PubMed Abstract | CrossRef Full Text | Google Scholar

Riccio, G., Ruocco, N., Mutalipassi, M., Costantini, M., Zupo, V., Coppola, D., et al. (2020). Ten-Year Research Update Review: Antiviral Activities from Marine Organisms. Biomolecules 10, 1007. doi:10.3390/biom10071007

PubMed Abstract | CrossRef Full Text | Google Scholar

Schrödinger, L. (2015a). LigPrep, Version 3.1. New York, NY: Schrödinger.

Google Scholar

Schrödinger, L. (2015b). LigPrep, Version 3. 5. New York, NY: Schrödinger, LLC.

Google Scholar

Shen, C.-H., Wang, Y.-F., Kovalevsky, A. Y., Harrison, R. W., and Weber, I. T. (2010). Amprenavir Complexes with HIV-1 Protease and its Drug-Resistant Mutants Altering Hydrophobic Clusters. FEBS J. 277, 3699–3714. doi:10.1111/j.1742-4658.2010.07771.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, L. W., Mao, H. J., Wu, Y. L., Tanaka, Y., and Zhang, W. (2017). TMPRSS2: A Potential Target for Treatment of Influenza Virus and Coronavirus Infections. Biochimie 142, 1–10. doi:10.1016/j.biochi.2017.07.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Shirato, K., Kawase, M., and Matsuyama, S. (2013). Middle East Respiratory Syndrome Coronavirus Infection Mediated by the Transmembrane Serine Protease TMPRSS2. J. Virol. 87, 12552–12561. doi:10.1128/jvi.01890-13

CrossRef Full Text | Google Scholar

Shulla, A., Heald-Sargent, T., Subramanya, G., Zhao, J., Perlman, S., and Gallagher, T. (2011). A Transmembrane Serine Protease Is Linked to the Severe Acute Respiratory Syndrome Coronavirus Receptor and Activates Virus Entry. J. Virol. 85, 873–882. doi:10.1128/jvi.02062-10

CrossRef Full Text | Google Scholar

Singh, R., Gautam, A., Chandel, S., Ghosh, A., Dey, D., Roy, S., et al. (2020). Protease Inhibitory Effect of Natural Polyphenolic Compounds on SARS-CoV-2: An In Silico Study. Molecules 25, 4604. doi:10.3390/molecules25204604

PubMed Abstract | CrossRef Full Text | Google Scholar

Søndergaard, C. R., Olsson, M. H., Rostkowski, M., and Jensen, J. H. (2011). Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of pKa Values. J. Chem. Theor. Comput. 7, 2284–2295. doi:10.1021/ct200133y

CrossRef Full Text | Google Scholar

Song, W., Gui, M., Wang, X., and Xiang, Y. (2018). Cryo-EM Structure of the SARS Coronavirus Spike Glycoprotein in Complex with its Host Cell Receptor ACE2. Plos Pathog. 14, e1007236. doi:10.1371/journal.ppat.1007236

PubMed Abstract | CrossRef Full Text | Google Scholar

Stewart, J. J. P. (1991). Optimization of Parameters for Semiempirical Methods. III Extension of PM3 to Be, Mg, Zn, Ga, Ge, As, Se, Cd, in, Sn, Sb, Te, Hg, Tl, Pb, and Bi. J. Comput. Chem. 12, 320–341. doi:10.1002/jcc.540120306

CrossRef Full Text | Google Scholar

Stopsack, K. H., Mucci, L. A., Antonarakis, E. S., Nelson, P. S., and Kantoff, P. W. (2020). TMPRSS2 and COVID-19: Serendipity or Opportunity for Intervention. Cancer Discov. 10, 779–782. doi:10.1158/2159-8290.cd-20-0451

PubMed Abstract | CrossRef Full Text | Google Scholar

Vidler, L. R., Brown, N., Knapp, S., and Hoelder, S. (2012). Druggability Analysis and Structural Classification of Bromodomain Acetyl-Lysine Binding Sites. J. Med. Chem. 55, 7346–7359. doi:10.1021/jm300346w

CrossRef Full Text | Google Scholar

Vivek-Ananth, R. P., Rana, A., Rajan, N., Biswal, H. S., and Samal, A. (2020). Silico Identification of Potential Natural Product Inhibitors of Human Proteases Key to SARS-CoV-2 Infection. Molecules 25, 3822. doi:10.3390/molecules25173822

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan, Y., Shang, J., Graham, R., Baric, R. S., and Li, F. (2020). Receptor Recognition by the Novel Coronavirus from Wuhan: An Analysis Based on Decade-Long Structural Studies of SARS Coronavirus. J. Virol. 94, e00127–20. doi:10.1128/JVI.00127-20

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, C., Horby, P. W., Hayden, F. G., and Gao, G. F. (2020). A Novel Coronavirus Outbreak of Global Health Concern. The lancet 395, 470–473. doi:10.1016/s0140-6736(20)30185-9

CrossRef Full Text | Google Scholar

Wang, E., Sun, H., Wang, J., Wang, Z., Liu, H., Zhang, J. Z. H., et al. (2019). End-Point Binding Free Energy Calculation with MM/PBSA and MM/GBSA: Strategies and Applications in Drug Design. Chem. Rev. 119, 9478–9508. doi:10.1021/acs.chemrev.9b00055

PubMed Abstract | CrossRef Full Text | Google Scholar

Waterhouse, A., Bertoni, M., Bienert, S., Studer, G., Tauriello, G., Gumienny, R., et al. (2018). SWISS-MODEL: Homology Modelling of Protein Structures and Complexes. Nucleic Acids Res. 46, W296–W303. doi:10.1093/nar/gky427

PubMed Abstract | CrossRef Full Text | Google Scholar

WHO (2021). Draft Landscape and Tracker of COVID-19. candidate vaccines-21 February [Online]. Available: https://www.who.int/publications/m/item/draft-landscape-of-covid-19-candidate-vaccines (Accessed October 1, 2021).

Google Scholar

Wilson, S., Greer, B., Hooper, J., Zijlstra, A., Walker, B., Quigley, J., et al. (2005). The Membrane-Anchored Serine Protease, TMPRSS2, Activates PAR-2 in Prostate Cancer Cells. Biochem. J. 388, 967–972. doi:10.1042/bj20041066

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Q., Peng, Z., Zhang, Y., and Yang, J. (2018). COACH-D: Improved Protein-Ligand Binding Sites Prediction with Refined Ligand-Binding Poses through Molecular Docking. Nucleic Acids Res. 46, W438–W442. doi:10.1093/nar/gky439

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang, Z. (2006). Advances in Homology Protein Structure Modeling. Cpps 7, 217–227. doi:10.2174/138920306777452312

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamamoto, M., Matsuyama, S., Li, X., Takeda, M., Kawaguchi, Y., Inoue, J.-I., et al. (2016). Identification of Nafamostat as a Potent Inhibitor of Middle East Respiratory Syndrome Coronavirus S Protein-Mediated Membrane Fusion Using the Split-Protein-Based Cell-Cell Fusion Assay. Antimicrob. Agents Chemother. 60, 6532–6539. doi:10.1128/aac.01043-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Roy, A., and Zhang, Y. (2013). Protein-ligand Binding Site Recognition Using Complementary Binding-Specific Substructure Comparison and Sequence Profile Alignment. Bioinformatics 29, 2588–2595. doi:10.1093/bioinformatics/btt447

PubMed Abstract | CrossRef Full Text | Google Scholar

Yi, M., Lin, S., Zhang, B., Jin, H., and Ding, L. (2020). Antiviral Potential of Natural Products from marine Microbes. Eur. J. Med. Chem. 207, 112790. doi:10.1016/j.ejmech.2020.112790

CrossRef Full Text | Google Scholar

Zhang, M.-Q., and Wilkinson, B. (2007). Drug Discovery Beyond the 'Rule-Of-Five'. Curr. Opin. Biotechnol. 18, 478–488. doi:10.1016/j.copbio.2007.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, P., Yang, X.-L., Wang, X.-G., Hu, B., Zhang, L., Zhang, W., et al. (2020). A Pneumonia Outbreak Associated with a New Coronavirus of Probable Bat Origin. nature 579, 270–273. doi:10.1038/s41586-020-2012-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y., Vedantham, P., Lu, K., Agudelo, J., Carrion, R., Nunneley, J. W., et al. (2015). Protease Inhibitors Targeting Coronavirus and Filovirus Entry. Antivir. Res. 116, 76–84. doi:10.1016/j.antiviral.2015.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, N., Zhang, D., Wang, W., Li, X., Yang, B., Song, J., et al. (2020). A Novel Coronavirus from Patients with Pneumonia in China. New Engl. J. Med. 382 (8), 727–733. doi:10.1056/NEJMoa2001017

CrossRef Full Text | Google Scholar

Keywords: SARS-CoV-2, COVID-19, tmprss2, molecular modeling, molecular docking

Citation: Mahmudpour M, Nabipour I, Keshavarz M and Farrokhnia M (2021) Virtual Screening on Marine Natural Products for Discovering TMPRSS2 Inhibitors. Front. Chem. 9:722633. doi: 10.3389/fchem.2021.722633

Received: 09 June 2021; Accepted: 17 September 2021;
Published: 12 October 2021.

Edited by:

Sylvain Broussy, University of Paris, France

Reviewed by:

Andrei I. Khlebnikov, Tomsk Polytechnic University, Russia
Luca Rastrelli, University of Salerno, Italy

Copyright © 2021 Mahmudpour, Nabipour, Keshavarz and Farrokhnia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Maryam Farrokhnia, farrokhnia.m@gmail.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.