Next Article in Journal
Changes in Porcine Corpus Luteum Proteome Associated with Development, Maintenance, Regression, and Rescue during Estrous Cycle and Early Pregnancy
Next Article in Special Issue
Combined Pharmacophore and Grid-Independent Molecular Descriptors (GRIND) Analysis to Probe 3D Features of Inositol 1,4,5-Trisphosphate Receptor (IP3R) Inhibitors in Cancer
Previous Article in Journal
Chenodeoxycholic Acid Has Non-Thermogenic, Mitodynamic Anti-Obesity Effects in an In Vitro CRISPR/Cas9 Model of Bile Acid Receptor TGR5 Knockdown
Previous Article in Special Issue
Anti-Oxidation and Anti-Inflammatory Potency Evaluation of Ferulic Acid Derivatives Obtained through Virtual Screening
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Natural Products-Based Drug Design against SARS-CoV-2 Mpro 3CLpro

1
Graduate Program on Medicinal Chemistry and Molecular Modeling, Institute of Health Science, Federal University of Pará, Augusto Corrêa 01-Guamá, Belém 66075-110, PA, Brazil
2
Departamento de Química, Faculdade de Filosofia, Ciências e Letras de Ribeirão Preto, Universidade de São Paulo, Ribeirão Preto 14040-901, SP, Brazil
3
Graduate Program on Pharmacy (PPGFAR-UFBA), Pharmacy College, Federal University of Bahia, Salvador 40170-115, BA, Brazil
4
Laboratory of Bioinformatics and Molecular Modeling (LaBiMM), Pharmacy College, Federal University of Bahia, Av. Barão de Jeremoabo, 147 Ondina, Salvador 40170-115, BA, Brazil
5
Department of Pharmaceutical and Organic Chemistry, Faculty of Pharmacy, Campus of Cartuja, University of Granada, 18071 Granada, Spain
6
Biosanitary Institute of Granada (ibs.GRANADA), University of Granada, 18071 Granada, Spain
7
Department of Physical Sciences, University of Embu, Embu 6-60100, Kenya
8
School of Pharmaceutical Sciences of Ribeirao Preto, University of Sao Paulo, Ribeirao Preto 14040-903, SP, Brazil
9
Laboratory of Modeling and Computational Chemistry, Department of Biological and Health Sciences, Federal University of Amapá, Macapá 68902-280, AP, Brazil
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(21), 11739; https://doi.org/10.3390/ijms222111739
Submission received: 14 September 2021 / Revised: 15 October 2021 / Accepted: 18 October 2021 / Published: 29 October 2021
(This article belongs to the Special Issue Drug Design and Virtual Screening)

Abstract

:
Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has received global attention due to the serious threat it poses to public health. Since the outbreak in December 2019, millions of people have been affected and its rapid global spread has led to an upsurge in the search for treatment. To discover hit compounds that can be used alone or in combination with repositioned drugs, we first analyzed the pharmacokinetic and toxicological properties of natural products from Brazil’s semiarid region. After, we analyzed the site prediction and druggability of the SARS-CoV-2 main protease (Mpro), followed by docking and molecular dynamics simulation. The best SARS-CoV-2 Mpro complexes revealed that other sites were accessed, confirming that our approach could be employed as a suitable starting protocol for ligand prioritization, reinforcing the importance of catalytic cysteine-histidine residues and providing new structural data that could increase the antiviral development mainly against SARS-CoV-2. Here, we selected 10 molecules that could be in vitro assayed in response to COVID-19. Two compounds (b01 and b02) suggest a better potential for interaction with SARS-CoV-2 Mpro and could be further studied.

1. Introduction

Although coronaviruses (CoVs) have been known since 1940 [1], reports on human infections causing mild respiratory infections were first reported in the 1960s. In December 2019 [2,3,4,5,6], a new viral respiratory disease emerged from Wuhan, China, and rapidly spread globally. This virus was named SARS-CoV-2 because its genomic RNA is about 82% identical to the coronavirus that causes Severe Acute Respiratory Syndrome (SARS-CoV). The disease caused by SARS-CoV-2 is called COVID-19 [4]. On 11 March 2020, the World Health Organization (WHO) declared the COVID-19 outbreak a pandemic [3,4,5] and data indicate that the cumulative number of cases reported globally now exceeds 183 million, with fatalities of almost 4 million (https://www.who.int/emergencies/diseases/novel-coronavirus-2019, accessed on 9 June 2021).
Therefore, the rapid discovery of safe, effective, and broad-spectrum anti-COVID-19 drugs is urgent. Currently, there is no specific therapy for COVID-19 patients and drugs are limited [7], despite more than 80 clinical trials employing distinct drugs such as Chloroquine, Hydroxychloroquine, Arbidol, Remdesivir, Favipiravir, Lopinavir, Ritonavir, Oseltamivir, Methylprednisolone, and Bevacizumab [1,8]. Due to its relevance to viral replication, the Mpro enzyme, which plays a central role in mediating viral replication and transcription [4,7,9,10,11], was the main target investigated in this study and the compounds identified here will be evaluated in future assays.
Individuals infected with COVID-19 have flu-like symptoms with fever, headache, muscular soreness, cough, and dyspnea, and in a few cases, were reported vomiting and diarrhea. People with asymptomatic COVID-19 behaviors may have a longer incubation period, thus it is indicated that quick tracking of the infected individual is essential to limit person-to-person transmissions. The quarantine period of 14 days for travelers was considered recommendable. The reduction of social interaction by limiting activities at workplaces, education centers, and public functions were indicated to be more effective in preventing and reducing the rapid spread of infection [12,13,14,15].
Recent studies assume that supplementing one’s diet has been widely reported to improve the immunity, health, and microbiota of humans, which play a key role in the treatment of COVID-19, even in infected cases; healthy eating habits and symbiosis decrease the chances of acute inflammation and severe cases [16].
Natural products research and development (R&D) potentially plays a pivotal role in innovative drug discovery [17]. Medicinal plants have been employed by humans for thousands of years to treat illnesses, health disorders [18], and to identify some plants useful as antiviral drugs. Plants produce a high diversity of secondary metabolites with interesting biological activities and both their natural products and natural product-inspired agents have attracted significant attention because they have played an integral role in the treatment of many different conditions, such as viral infection, especially SARS-CoV-2 [18,19,20,21,22,23,24]. Carolacton, Homoharringtonine, Emetine, and Cepharanthine, as well as natural product-inspired small molecules (Ivermectin, GS-5734, EIDD-2801, and Ebselen), are potential anti-SARS-CoV-2 agents that have attracted significant attention due to their broad-spectrum antiviral activities [24].
A good starting point to find new antiviral natural products would be traditional medicinal plants. Brazil hosts about 20% of all the world’s biodiversity, being considered the country with the most endemic species worldwide [25]. The Natural Products Database of the Bahia Semi-Arid region (NatProDB, http://natprodb.uefs.br/, accessed on 7 May 2021) contains naturally occurring compounds with a wide chemical diversity. These compounds are yet to be fully explored for the discovery of bioactive molecules [26]. Thus, we searched for new molecular entities that could be potent inhibitors of SARS-CoV-2 Mpro in the NatProDB database.
Medicinal chemistry approaches have seen a great technological advancement over the years and have contributed immensely to the discovery of promising molecules. However, many compounds discovered by these techniques have shown unsatisfactory absorption, distribution, metabolism, excretion, and toxicity (ADMET) properties in in vivo tests [27].
In this study, we propose to employ a structure-based drug design (SBDD) and ligand-based drug design (LBDD) approach against the main sites of the SARS-CoV-2 main protein (Mpro) through different computer-aided methods. We explore ADMET properties using two distinct tools, namely QikPropTM and DerekTM, intending to reduce the selection of false-positive ligands and obtain hits that can be subjected to in vitro tests. We further evaluated the most likely corresponding bioactive conformations through tactics such as cavity detection, docking, and molecular dynamics (MD) simulations. All steps performed in this study were as described in Figure 1.

2. Results

In our previous work, we analyzed the NatProDB compound [28] interactions with a Trypanosoma cruzi enzyme. Here, we decided to analyze these data in further detail, employing additional distinct ADMET in silico methods: QikPropTM and DerekTM [29].

2.1. ADMET Analysis by QikPropTM

ADMET prediction was performed through the QikPropTM module. This tool predicts some relevant properties (e.g., water/gas log Ps, log S, log BB, overall CNS activity, Caco-2, and MDCK cell permeabilities; log Khsa for human serum albumin-binding and log IC50 for HERG K+-channel blockage) based on the full 3D molecular structure over twenty physical descriptors [30].
First, we collected 555 molecules from the NatProDB using the filter #star to calculate the number of properties that were outside the range by comparing a particular molecule’s properties with those of the 95% of known drugs via QikPropTM [30,31,32]. Although it is acceptable to have at least five properties per hit outside the range, in this study, all hits with #star values of ≥1 were eliminated following the E64 (positive control), yielding 376 compounds. Subsequently, we evaluated the predicted activity on the central nervous system (CNS), which ranged from −2 (inactive) to +2 (active) [32,33,34]. Only compounds with values of <0 were kept (203 hits) considering that the antiviral action mechanism proposed in this study does not necessarily require an effect on the CNS.
Next, the human oral absorption properties [35] were predicted for the remaining compounds using the values of a percentage ranging from >80% (high and desirable) to ≤25% (low absorption) as standard criterion [32,33]. Aiming to search for good absorption potential, molecules with values of ≥60% were maintained since E-64 presented 14.41%, obtaining 201 candidates.
To assess the apparent permeability on the Caco-2 cell membrane, the Boehringer-Ingelheim scale (nm/s) was employed. Compounds with low permeability presented <5 nm/s and those with high permeability showed values of >500 nm/s [30,31,32]. In this process, only molecules with values above 100 nm/s were considered and yielded in the 153 molecules with majorly pharmacokinetic desirable properties for the top 10 hits (Supplementary Table S1).

2.2. Toxicological Analysis by Derek NexusTM

For the toxicity analysis predictions, we sought to identify and eliminate molecules that exhibited any toxicity alert (carcinogenicity, mutagenicity, genotoxicity, skin sensitization, teratogenicity, nephrotoxicity, hepatotoxicity, irritation, respiratory sensitization, and reproductive effects) [29] for specific chemical groups in the early stages of the process for developing prototype drugs, making it possible to optimize costs and time during research. Thus, the evaluation of chemical groups of the remaining 153 NatProDB compounds resulted in 47 candidates with low toxicological potential and adequate pharmacokinetic profiles for an oral route, i.e., good intestinal absorption, high cellular permeability, and low excretion, thereby promising antiviral candidate drugs.
The polar molecular surface (PSA) area is a property based on polar atom surface that was shown to correlate with passive molecular transport through membranes and allow for the prediction of the transport properties of drugs [36,37]. Its desirable value ranges from 7.0 to 200.0 Å2 [30,31,32]. Molecule 6 had the highest PSA index of 107.35 Å2, while compound 10 had the lowest value at 46.80 Å2. This indicates a good permeability of the selected molecules through the lipid membrane. It is noteworthy that a PSA value of <75 Å2 has been associated with an increased risk of adverse effects when combined with high lipophilicity (logP greater than 4) and only one matched that situation, which was excluded [36].
A drug bound to a protein is inactive and it is only the free drug that can bind with a receptor, providing pharmacological action [38]. The free drug concentration is affected by binding on human serum albumin (HSA), therefore this pharmacokinetic property was also predicted for the compounds [39,40]. This was evaluated by its logarithmic constant and has a recommended range of −1.5 to 1.5 [30,32,41]. Values outside the standard range may indicate non-enzymatic glycosylation of HSA interference on drug-binding potency, and may influence the distribution and excretion of molecules affecting the duration and intensity of the biologic effects [40,42,43]. All molecules selected showed satisfactory values that were better than the E64, which was the smallest in the series, specifically logKHSA = −1.31 (Table S1).
Finally, the hERG encodes the inward rectifying voltage-gated potassium channel in the heart (IKr), which is associated with cardiac repolarization. Inhibition of the hERG current induces QT interval prolongation and could represent fatal ventricular arrhythmia called Torsade de Pointes [44]. In fact, recent studies pointed out that some medications such as Hydroxychloroquine and azithromycin have facilitated the risk of fatal arrhythmia in patients with COVID-19 infection [45,46]. In our study, molecules that presented QPlogHERG values below −5 (2 = −5.14; 4 = −5.53 and 5 = −6.48) required greater attention and further in vitro analysis should be done to validate their pharmacological safety in the treatment against COVID-19.

2.3. Site Prediction and Druggability Analysis of SARS-CoV-2 Mpro

We noted that all programs correctly identified the SARS-CoV-2 Mpro active site (Figure 2A–C) and additionally indicated that this site was druggable, e.g., with a druggability probability of >0.5 in PockDrug [47] as well as a strong primary hot spot with ≥16 probes clusters and hot spots connected with a center-to-center distance of <8 Å in FTMap [48,49]. These agree with all experimental SARS-CoV-2 Mpro structures [4,7,9,50,51], where distinct ligands were tested and crystallized.
After the protein preparation step, pocket predictions and the druggability analysis of the SARS-CoV-2 Mpro were done using three different methods: Difference of Gaussian (DoG) [52,53], PockDrug [47], and FTMap [49]. The main pockets determined in this protein are shown below (Figure 1)
Other sites were also signed as druggable in our analysis: SeeSAR (II and III—Figure 2A); PockDrug (green, yellow, and gray sites—Figure 2B); and FTMap (yellow—Figure 2C). Since two SARS-CoV-2 Mpro sites were located in the same locations and were specified by distinct methods (II—yellow surface and yellow mesh; III—green surface and yellow mesh, Figure 2), we hypothesized that these two regions could be explored in further studies.
The probes posed by FTMap revealed that the druggable site (magenta sticks, Figure 2C) accommodates polar acyclic compounds (ethanol, dimethylformamide, and acetamide) and some rings (benzene, benzaldehyde, and cyclohexane), and the borderline site (yellow sticks, Figure 2C) lodges ethanol, methanamine, benzene, and cyclohexane. This information is helpful for fragment-based drug design. In another further step, gray sites predicted by PockDrug (Figure 2B), which had higher druggability scores, could also be explored.

2.4. Molecular Docking Simulations

The docking calculations were done by SeeSAR version 9.2 (BioSolveIT GmbH) and focused on the active site of SARS-CoV-2 Mpro to better sample the interaction energy (ΔG) between NatProDB ligands and this enzyme. The hydrogen bond and dehydration (HYDE) energy scoring function in SeeSAR provides a range of affinities, spanning an upper and lower limit.
We, therefore, used this descending range as an affinity parameter to compare distinct molecules and their interactions on the Mpro active site. For the sake of clarity, we only show the 10 best-docked compounds with higher estimated affinities (EA and µM) in Table 1 instead of 46 NatProDB and E-64 molecules as previously filtered from the ADMET analysis (see preceding discussions).
SeeSAR version 9.2 (BioSolveIT GmbH) analyzes the estimated affinities’ (EA) energies and the best poses selection was based on both the HYDE scoring function [61,62,63] and torsional analysis [64]. The docked molecules at the SARS-CoV-2 Mpro active site were ordered by their EA on SARS-CoV-2 Mpro-NatProDB complexes and the “top 10” were selected for further analysis.
In comparing the docking results presented in Table 1, we noted that NatProDB compounds’ energy ranked in the “top 10” were better ranked than the Mpro crystallographic ligand (RZS, 255th better energy). These results indicate a good interaction between NatProDB and SARS-CoV-2 Mpro residues, mainly on the active site. Additionally, the “top 10” NatProDB compounds presented better energy when compared to a positive control (E-64), which was a well-known cysteine protease inhibitor (562th better energy, Table 1) [65,66].
The best-scored poses were selected by visual inspection and, since they interacted better with Mpro active site residues, we decided to compare their interactions by analyzing each complex (Mpro- RZS, E-64, and NatProDB), aiming to find some structural information that could help regarding the drug discovery against COVID-19 (Figures S1–S4).
These hits occupied the active site in a varied manner and interacted with Mpro through hydrogen bonds with catalytic residues (His41 and Cys145), as well as made additional hydrophobic contacts as shown from the docking results (Figures S1–S4).
The main interactions were described by the Protein–Ligand Interaction Profiler server (https://plip-tool.biotec.tu-dresden.de/plip-web/plip/index, accessed on 27 May 2021) [67] at each pose (Table S2), where we noted that the two better-docked NatProDB compounds (b01 and b02) presented better contacts than the crystallographic ligand (RZS).

2.5. Molecular Dynamics (MD) Simulations

To understand the behavior of SARS-CoV-2 Mpro complexes, we analyzed them through molecular dynamics (MD) simulations. These data allowed us to analyze the time influence on these interactions’ patterns. First, we analyzed the stability of our eight systems: apo and complexes: holo (PDB ID: 5R82), as well as E-64 and the five better-ranked NatProDB compounds (Table S3). From our root-mean square deviation (RMSD) data, we noted that all the systems equilibrated after 40 ns (Figure S5A). Then, we defined our productive phase into time intervals from 40 to 100 ns for all the simulations (Figure S5B,C).
In comparing the Mpro structure fluctuations (3D RMSF) within all the complexes, we distinguished their overall protein-folding stability, which was mainly at the active site represented by catalytic residues (His41 and Cys145, Figure 3). These data correlated with the Mpro crystallographic structure, where active site residues were well solved in the electron density maps from 2.1 Å resolution [7].
Additionally, we analyzed their secondary structure stability by the DSSP 3.1.4 module [68,69,70], which was installed on GROMACS 5.1.4 [71,72,73,74,75,76]. All Mpro complexes maintained their secondary structure stability during our simulation time (Figure S6). Some regions presented higher fluctuations in Mpro structure, mostly at some loops and at both N and C-terminals (Figure 3). The SARS-CoV-2 Mpro crystal had a long loop region of domain II (residues 185–200) connecting to domain III and this region is highly variable as emphasized by the superposition of 12 Mpro crystal structures [7].
Our MD data revealed that these loops on domain II fluctuated higher than other secondary protein structures, corroborating with this experimental study. Previous results also cited the interface region among Mpro protomers composed by domain II residues as unstable [4,7,9] and these data were also visualized in our simulations (Figure 3).
Many studies described Mpro subsites expanding from S6 to S1 with catalytic sites located on S1 [4,7,9,51]. Based on these results, we analyzed the interaction pattern on SARS-CoV-2 Mpro of NatProDB derivatives compared with the crystallographic ligand (RZS) and E-64. Additionally, we evaluated the hydrogen bond pattern, through the GROMACS H-bond module [77] and HbMap2Grace program [78], and the molecular surface area, by the SurfinMD program [79]. The molecular dynamics (MD) data for the hydrogen bond (H-bond) analysis showed that NatProDB compounds exhibited more interactions than RZS, and b04 and b05 were similar to E-64, which has the highest H-bond number (Table S3). Since these interaction patterns could favor their inhibitory behavior, we analyzed them in detail. Initially, we defined and calculated the hydrogen bond capability for our ligands (Table S3), (Hbondcapac.), during the MD simulation. This measure could aid us in quantifying how “tightly” the ligands interacted with Mpro since Hbondcapac. ≥ 1 indicates that each ligand’s hydrogen bond atoms (donors and acceptors) interacted with protein residues.
Our results from the MD productive phase correlates well with the estimated binding energy obtained before from the docking results (Table 1), i.e., the ligands with the best docking energies presented higher Hbondcapac., which revealed to us that the hydrogen bond is a favorable interaction in developing Mpro inhibitors. Additionally, we noted that our ligand accessed distinct sub-sites on SARS-CoV-2 Mpro, interacting beyond the catalytic ones—S1 (Figure 4).
The H-bond pattern presented interactions with His41 and Cys145 (catalytic residues); Thr25 and 26 from S1′; Gly143 and Ser144 from the canonical oxyanion role on S1 [7]; Glu166 from S1 and Gln189; and Thr190 from S4 [4,7,9,51]. Since the H-bond is considered as the “driving force” for Mpro inhibition, we could map and analyze them through our MD simulations.
We also calculated the atomic contacts involving SARS-CoV-2 Mpro and NatProDB compounds (Figure 5). The contact surface area revealed additional interactions with apolar residues on the same sites described before [4,7,9]: S1′ (Thr25, 26, Leu 27, and Cys145), S2 (His41, Thr45, Ser46, Asp48, and Met49), S4 (Met165, Glu166, Leu 167, and Gln189), and S5 (Pro168).
Another measure to evaluate Mpro interaction with NatProDB ligands concerns their active site occupancy. Ideally, it is expected that potent inhibitors should present reversible action mechanisms and tight binding [80]. Since Mpro active sites have a well-defined area (335.9 Å2) and volume (364.1 Å3) [81], we compared these results with our MD data. The Fpocket program [82,83] was employed to calculate the volume variation throughout the simulation. We noted that the active site volume oscillates from 461.8 Å3 of E-64 complexes to 262.9 Å3 for b04 (Table S3), which agreed well with previous data [4,50,81] and with both the molecular area and volume of NatProDB ligands (Table 1). These results could emphasize that SARS-CoV-2 Mpro and NatProDB ligands show induced-fit mechanisms, correlating with their dynamics and molecular interactions at the active site.
Based on the interaction results of complexes after MD simulations, we assumed that new Mpro sites were pointed to: (i) V42, C44, T45 and S46, and both near and on (ii) Y118 and N119 (Table S3). The former residues made hydrophobic interactions with alpha-ketoamide inhibitors as described before [4], including FDA-approved drugs, medicinal plant compounds [81], and ZINC-15-docked compounds [84]. Our results (Figure 6) reinforce the importance of these non-polar interactions and opens up a new region (together with latter residues) for drug discovery development, aiming to search for new Mpro inhibitors.
Additionally, we calculated the binding free energy of all Mpro complexes through MM-PBSA methods [85,86]. The binding energy (ΔEbinding) calculated by the solvent accessible surface area showed that all compounds interacted favorably with Mpro (Table S3). Since these values were directly correlated to interacting protein residues, we decided to discriminate how amino acids presented better contacts with ligands. These residue decomposition energy analyses selected residues near the ligand (<5 Å) during the MD simulation and also those that participated actively in complex stabilization (ΔEbinding > ±5 kJ/mol), as shown in Figure 7.
We noted that ligand interactions with catalytic residues (His41 and Cys145) are highly favorable (negative values), as expected for reversible inhibitors. This behavior remained for two residues of the new binding site (Arg188 and Gln189), as first described in this work. Other residues (Val42, Cys44, Ser46, and Glu48) did not interact favorably with the ligands (positive values—Figure 6). This could be related to their position on the Mpro site, i.e., their side chains were pointed outwards instead of towards the active site cavity [4,9].

3. Discussion

Natural products represent an important source in discovering new and effective therapies to treat patients with COVID-19. Plant extracts are used around the planet as natural immune stimulants and anti-infective agents [87]. A good starting point to find new antiviral natural products would be traditional medicinal plants, such as those from Asia, Africa, or America that have been employed to treat infections. Progressive understanding of the development of efficient antiviral drugs and their bioactivity mechanisms has encouraged the exploration of natural products as an important way to indicate effective treatments for patients with COVID-19. Computational studies obtain and provide comparative approaches in a fast and economical way [88,89].
Carolacton is an antibacterial macrolide keto-carboxylic acid synthetized by the myxobacterium Sorangium cellulosum [90] and is reported as an inhibitor of Methylenetetrahydrofolate dehydrogenase 1 (MTHFD1), which is a potential target for developing anti-SARS-CoV-2 agents. The Carolacton strongly inhibited SARS-CoV-2 replication in Vero cells at a half-maximal inhibitory concentration (IC50) and exhibited a moderate cytotoxicity (CC50) profile [24,91].
Tetrahydroisoquinoline alkaloid emetine, a known natural product that was isolated from the plant Psychotria ipecacuanha [92], showed potential cardiotoxicity and acted as a protein synthesis inhibitor [93]; notably, this natural product has very recently been recognized as a promising antiviral drug with in vitro activity against MHV-A59 [94], severe acute respiratory syndrome coronavirus (SARS-CoV), and Middle East respiratory syndrome coronavirus (MERS-CoV) [95]. Recently, researchers performed a virtual screening against the SARS-CoV-2 Mpro binding site using the library of Marine Natural Products. Through many comparative prediction approaches, they selected seventeen most promising compounds in terms of being potential SARS-CoV-2 inhibitors [96].
The Natural Products Database of the Bahia Semi-Arid region (NatProDB, http://natprodb.uefs.br/, accessed on 7 May 2021) contains naturally occurring compounds with a wide chemical diversity. NatProDB phytocompounds were explored for their antiviral potential using a ligand and structure-based drug design approach against the main sites on the SARS-CoV-2 main protein (Mpro) through different computer-aided methods. Many bioactive compounds do not progress to clinical phase trials due to undesirable pharmacokinetic properties identified in advanced stages of the drug design, representing a significant loss of money, equipment, and time. Thus, it is essential to understand the pharmacokinetic profile of potential drug candidates at an early stage of research and development [41,97,98].
Compound b05 had the highest PSA index (107.35 Å2), listed in Table S1, compared to both the b08 and E-64 control (46.80 Å2), which were the lowest values of the series. Such evidence suggests good permeability of the molecule in a cellular lipid membrane [35]. A low PSA (PSA < 75 Å2) has been associated with an increased risk of adverse events due to non-specific toxicity, particularly when combined with high lipophilicity (log P > 4) [36]. A high PSA value has been associated with poor membrane permeability and, in particular, with low blood–brain barrier penetration, which denotes a specific structure that protects the central nervous system (CNS) [99]. For good oral bioavailability, an upper limit of 140 Å2 has been suggested, particularly when associated with a large number of rotatable bond surfaces belonging to polar atoms [100].
Cerebrovascular diseases are among the comorbidities of patients with confirmed COVID-19 who develop severe respiratory complications [101]. The intense systemic inflammatory response linked to viral infection can lead to blood–brain barrier (BBB) breakdown. This in turn can allow peripheral cytokines to gain access to the CNS, where they may trigger or exacerbate neuroinflammation, leading to encephalitis [102,103]. Recently, studies estimated the incidence of neurological syndromes to be about 37% [104,105,106]. In Table S1, it is possible to visualize that the b01, b04, b05, and b08 compounds presented a CNS activity value of −1, suggesting that they can be more active than the E-64 control (−2). For a drug with biological activity in the CNS region, a value higher than −2 is required [107,108].
The angiotensin-converting enzyme (ACE2) is expressed in Caco-2, Calu-3, and Vero-6 cells on the apical membrane domains [109,110]. A study investigated a library of thousands of marketed drugs and prototypes to evaluate the inhibition of virus-induced cytotoxicity using Caco-2 cell lines from human epithelial colorectal adenocarcinoma and SARS-CoV-2 obtained from a patient originally infected with the virus in Wuhan, China. Thus, 19 compounds with IC50 < 1 μM were identified, including the Remdesivir [111], exhibiting the same inhibition of SARS-CoV-2 in Caco-2 and Vero-6 cells (IC50 0.77 μM), and surprisingly, promising potency (IC50 0.07 μM) in SARS-CoV and MERS [112]. In our analysis, only hits which presented high Caco-2 values considered. These included compounds b02 (759.23 nm/s) and b04 (1101.54 nm/s). The others exhibited values between 5 nm/s and 500 nm/s. All hits had better values than the E-64 control (1.52 nm/s), considering the low apparent permeability on the Caco-2 cell membrane. This suggests that the compounds may be transported across the host cellular membrane or may block the virus–host membrane interaction that is critical for the progression of COVID-19 infection.
To achieve better permeability, the drug should have a moderate logP value (−2 to 6.5). In this range, there is likely to be a good balance between permeability and solubility [113,114]. LogP values in Table S1 show that three bioactive compounds had higher affinity for the organic phase, characterized as more hydrophobic, in the descending order of b01, b04, and b02. All bioactive compounds presented optimum logP values within the range, favoring good permeability and solubility for drugs, with molecule b03 (0.15) being the most hydrophilic. In our analysis, logP values were directly related to the values of human oral absorption (HOA%). The E-64 control presented a value (−3.93) much lower than necessary for what is considered adequate for logarithmic partition coefficients (n-octanol and water phases).
The site prediction and druggability analysis of SARS-CoV-2 Mpro was identified by the probes posed by FTMap and revealed the druggable site (magenta sticks, Figure 2C) which accommodates polar acyclic compounds (ethanol, dimethylformamide, and acetamide) and some rings (benzene, benzaldehyde, and cyclohexane), as well as the borderline site (yellow sticks, Figure 2C) which accommodates ethanol, methanamine, benzene, and cyclohexane. This information is helpful for fragment-based drug design. In a further step, gray sites predicted by PockDrug (Figure 2B), which had higher druggability scores, could also be explored. These agreed with all experimental SARS-CoV-2 Mpro structures, where distinct ligands were tested and crystallized [4,7,9,50,51].
The initial molecular docking analysis revealed that all the selected NatProDB phytocompounds against the SARS CoV-2 3CLpro occupied the active site in a varied manner, interacted with Mpro through hydrogen bonds with catalytic residues (His41 and Cys145), and made additional hydrophobic contacts as shown from the docking results in Figures S1–S4. The main interactions were described by the Protein–Ligand Interaction Profiler server (https://plip-tool.biotec.tu-dresden.de/plip-web/plip/index) [67] at each pose (Table S2), where we noted the two best-docked compounds, thus only the interactions of both with 3CLpro were analyzed in detail and showed satisfactory intervals between the lower and upper limits of the estimated affinities b01 (3.29–326.4 µM) and b02 (4.42–439.6 µM). In addition, both had better contacts than the crystallographic ligand (RZS), as depicted in Figures S1–S4.
The atomic clashes between the protein and docked ligand atoms (inter-clash type), as predicted by SeeSAR9.2 (BioSolveIT, Sankt Augustin, Germany, 2018, www.biosolveit.de/SeeSAR, accessed on 25 July 2021) were classified as + (no clashes), 0 (few clashes), and (many clashes) Table 1. All hits exhibited few or no clashes, including the RZS and E-64 controls. This suggests the appropriateness of the docked conformations generated from the docking tool and that these docked conformations can have equivalent intermolecular electron density overlaps similar to the experimental SARS-CoV-2 3CLpro data. Further validation of these conformations in MD simulations showed that the b01, b02, and RZS conformations show a stable interaction with SARS-CoV-2 3CLpro. Our results from the MD productive phase correlated well with the estimated binding energy obtained before from the docking results (Table 1). These ligands with the best docking energies presented higher Hbondcapac., which revealed to us that the hydrogen bond is a favorable interaction in developing Mpro inhibitors.
Ideally, it is expected that potent inhibitors should present reversible action mechanisms and tight binding [80]. Since the Mpro active site has a well-defined area (335.9 Å2) and volume (364.1 Å3) [81], we compared these results with our MD data. The Fpocket program [82,83] was employed to calculate the volume variation throughout the simulation. We noted that the active site volume oscillates from 461.8 Å3 of E-64 complexes to 262.9 Å3 for b04 (Table S3), which agreed well with previous data [4,50,81] and with both the molecular area and volume of NatProDB ligands (Table 1). These results could emphasize that SARS-CoV-2 Mpro and NatProDB ligands show induced-fit mechanisms, correlating with their dynamics and molecular interactions at the active site.
The H-bond pattern presented interactions with His41 and Cys145 (catalytic residues); Thr25 and 26 from S1′; Gly143 and Ser144 from the canonical oxyanion role on S1 [7]; Glu166 from S1 and Gln189; and Thr190 from S4 [4,7,9,51]. Hence, it can be concluded that the docking procedure generated correct poses for the tested compounds. Finally, binding free energy calculations and the decomposition analysis showed a good binding affinity of b01 and b02. These ligands interact with the key 3CLpro residues, including a strong interaction with the Cys145–His41 catalytic dyad (Figures S3 and S4), which is crucial for the 3CLpro function [115,116,117,118].
Altogether, the present investigation provides a comprehensive overview of Bahia Semi-Arid region (NatProDB) compounds and shows their potential against SARS-CoV-2 3CLpro. The in silico analysis elucidated two potential candidates, namely b01 and b02 as potential 3CLpro inhibitors. However, both are subject to some uncertainty and additional in vitro and in vivo studies are essential to validate their efficacy as good drug candidates for the inhibition of 3CLpro activity [119].

4. Materials and Methods

Molecular modelling calculations, estimations, and visualizations of SARS-CoV-2 Mpro binding affinities were carried out using the SeeSAR version 9.2 software package (BioSolveIT GmbH) with the HYDE visual affinity assessment [61,62,63] carried on a personal computer with 4 logical processors of an Intel® Core TM i5-6200U and 2.3 GHz processor using the Windows 10 Home operational system.

4.1. SARS-CoV-2 Mpro Protein Preparation

The high-resolution X-ray crystallographic structure of SARS-CoV-2 Mpro (1.31 Å) was obtained from the Protein Data Bank (PDB ID: 5R82, [120]) and was imported into SeeSAR version 9.2 (BioSolveIT, Sankt Augustin, Germany, 2018, www.biosolveit.de/SeeSAR, accessed on 25 July 2021). The structures were then prepared using the Protein Editor tool as implemented in the BioSolveIT Suite. Briefly, the raw PDB structure was processed by automatically assigning bond orders, adding hydrogens, adding missing side-chains, creating possible disulfide bridges, deleting waters, and generating hetero protonation states at pH 7.4. Residues with alternate positions were locked in the conformations with the highest average occupancy. Small ligand (6-(ethylamino) pyridine-3-carbonitrile-RZS) and dimethyl sulfoxide (DMS) originating from the crystallization buffer were removed. The hydrogen-bonding networks were optimized automatically by sampling water orientations and optimizations of hydroxyls, Asn, Gln, and His residue states using the Protein Editor.

4.2. Site Prediction and Druggability Analysis of SARS-CoV-2 Mpro

4.2.1. SeeSAR

The protein structures prepared as described above were exported as PDB files and inserted on the “Binding Site” tool using BioSolveIT Suite. The site predictions and druggability analysis were determined by applying the Difference of Gaussian (DoG) site pocket finder [52,53,121]. DoGSite provides the functionality to detect potential binding pockets and sub-pockets of a protein of interest. Subsequently, it analyzes the geometric and physico-chemical properties of these pockets and estimates druggability with the aid of a support vector machine (SVM) [121]. Thus, binding pockets and their druggability were extracted and analyzed, employing the default parameters.

4.2.2. PockDrug

PockDrug predicts pocket druggability on both pockets guided and not guided by the ligand proximity [35] using different estimation methods based on 36 physicochemical and 16 geometrical descriptors to characterize each estimated pocket [122]. The prepared PDB file was uploaded to the PockDrug server (http://pockdrug.rpbs.univ-paris-diderot.fr/cgi-bin/index.py?page=Home, accessed on 25 July 2021) and assessed for binding sites regarding their 6 pocket descriptors (volume of convex hull (Å3); Kyte scale of residues’ hydrophobicity [123]; frequency of polar residues; frequency of aromatic residues; oxygen atoms’ frequency from Tyr residues in pockets [48]; and number of pocket residues) and corresponding average druggability probability, as well as its associated standard deviation. For a probability greater than 0.5, pockets are considered as druggable according to the published data [47]. Local pocket properties, sizes, shapes, and hydrophobicity were extracted for all the identified sites and annotated to pocket numbers.

4.2.3. FTMap

The FTMap identifies binding hot spots of macromolecules (regions of the surface with major contributions to the ligand-binding free energy) using a computational mapping approach [54]. It uses small organic molecules as probes to sample the protein surface and scores them using the interaction energy [124]. Protein sites that bind different probes indicate this region as binding hot spots [54]. The prepared PDB file was uploaded to the FTMap server (https://ftmap.bu.edu/, accessed on 14 July 2021) and examined for the number of probes per cluster found according to the published protocols [64,125,126].
Similar to Borrel et al. [122], in this study, the predicted pocket of interest from the entire set of candidate cavities detected was selected as the one that best overlaps the crystallographic ligand-binding site and presented the best druggability index on all three programs.

4.2.4. Molecular Docking on SARS-CoV-2 Mpro

Docking experiments were performed on previously prepared SARS-CoV-2 Mpro structures with the SeeSAR version 9.2 through sampling their spatial coordinates on the best interaction sites, which were defined by both the DoG site pocket finder [52,121] and druggability analysis from SeeSAR, PockDrug, and FTMap. Then, the docking library was loaded, including a set of NatProDB ligands previously filtered from QikPropTM and DerekTM, in addition to E64, which is a well know cysteine protease inhibitor, as a positive control.
Docking calculations were done for each compound, generating 20 poses where each one has its binding affinity estimated from µM into the SARS-CoV-2 Mpro site. The best poses were selected based on their estimated affinity from the HYDE score function, while also considering the Torsion angle values to the binding conformation of the protein–ligand [61,121,127]. HYDE’s empirical score function (Equation (1)) relies on intrinsically balanced terms of atom-specific desolvation, hydration, and hydrogen bonds (H-bond; i.e., free of weighting parameters) based on the logP atomic increment system [128]. The binding affinity of a ligand, which is described as the quotient of ΔG and the non-hydrogen numbers of atoms in the compound (Equation (2)) [129], is expressed as follows.
ΔGHyde = Σatom i [ΔGi Dehydn + ΔGi h -bonds]
LE = ΔG⁄N
where ΔG = -RT ln Ki and N = number of non-hydrogen atoms.
The selection of the best poses was based on their visual HYDE scores while also considering a statistics-based torsional analysis [127]. SeeSAR version 9.2 enables an interactive assessment of torsions and energies (in KJ·mol−1), including the desolvation (dehydration, −TΔS) and enthalpic (interaction, ΔH) contributions to the binding for both the protein and ligand. Furthermore, it also quantitatively reports the energy contributions for all heavy atoms (with a united atom approach for bounded H-atoms) and allows for a semi-quantitative estimation of the thermodynamic profile for all the tested compounds. Analyses, comparisons, and visualizations of the predicted binding pose predictions and all images were done by educational Pymol version 2.4.1 [55,130].

4.2.5. SARS-CoV-2 Ligand Preparation

The Natural Products Database of the Bahia Semi-Arid region (NatProDB, http://natprodb.uefs.br/, accessed on 7 May 2021) contains 555 molecular structures publicly available [26]. All these molecules were automatically employed to perform an ADMET property analysis by QikPropTM and DerekTM.

4.2.6. ADMET Analysis by QikPropTM and DerekTM

The predictions described below were made on a computer with 6 logical processors of an Intel® Core TM i5-9400 and 4.10 GHz processor using the Windows 10 Professional operational system.
The QikProp program provides ranges to compare the physically significant descriptors and important ADMET properties of a specific compound with those of the 95% of approved drugs. Predicted human oral absorption on a 0 to 100% scale (HOA%) predicted apparent Caco-2 cell permeability in nm/sec (QPPCaco); predicted apparent Maden Darby Canine Kidney cell permeability (QPPMDCK); predicted the brain–blood partition coefficient (QPlogBB); and predicted the octanol–water partition coefficient (QPlogPo/w), among others. The program analyzes the chemical structure as fragments or as a whole and supports predictions about the structure in three dimensions. In addition, it flags 30 types of reactive functional groups that can cause false positives in high-throughput screening (HTS) assays [131,132].
Derek Nexus is an “expert rule-based” system for the prediction of toxicity [48,49,50] that uses a virtual base of molecules with reported toxicity alerts, including, for example, mutagenicity, genotoxicity, teratogenicity, carcinogenicity, and hepatotoxicity, for comparison with the potentially toxic group present in each query molecule [39]. The program generates an alert based on literature evidence and describes the potential toxicity for the complete structure [133].

4.2.7. Molecular Dynamics (MD) Simulations on SARS-CoV-2 Mpro-NatProDB Complexes

The SARS-CoV-2 Mpro systems apo, holo, and five complexed with better-docked NatProDB compounds were submitted to molecular dynamic (MD) simulations using GROMACS 5.1.4 [71,72,73,74,75,76,134], available at the National Center for High-Performance Computing in São Paulo (CENAPAD-SP), employing the following parameters: 1 atm, 310 K, pH 7.4, GROMOS54A7 force field-updated [135], and PME [136] for electrostatic treatment with cut-off = 1.0 nm in a dodecahedral box solvated with water model SPC/E [137] with periodic boundary conditions (PBC). Na+ and Cl− ions were added to maintain the physiological salt concentration (0.15M) and to neutralize the residual system charge at pH = 7.4. At first, the system was energy-minimized (steepest descent/conjugate gradient) until forces reached ≤ 10 kJ·mol−1 and nm−1, followed by a pre-equilibrium simulation step (heavy atoms’ position restrained for 1 ns).
Finally, 100 ns of the unrestrained simulation was performed for all SARS-CoV-2 Mpro systems, where atomic coordinates were recorded every 10 ps for later analyses in an NPT ensemble with a V-rescale thermostat [138] and Berendsen barostat [139] with the SETTLE [140] algorithm for solvent bonds and the LINCS [141] algorithm for other bonds. The crystallographic ligand (RZS) and NatProDB compounds had their topology coordinates built in the Automated Topology Builder (ATB) version 3.0 server http://compbio.biosci.uq.edu.au/atb/, accessed on 9 May 2021 [142,143,144].
The 3D root-mean square fluctuation (RMSF) analysis converted RMSF data to B-factor from the RMSF module of GROMACS 5.1.4 [71,72,73,74,75,76,134]. The hydrogen bond analysis was conducted in the GROMACS H-bond module [77] and HbMap2Grace program [78], available on http://lmdm.biof.ufrj.br/software/hbmap2grace/index.html, accessed on 22 July 2021 and the molecular surface area (Å2; Table 1) analysis was conducted by the SurfinMD program [79], accessible on http://lmdm.biof.ufrj.br/software/surfinmd/index.html, accessed on 9 May 2021.

4.2.8. Hydrogen Bond Capacity Analysis

Hydrogen bonds have been calculated between the number of donors/acceptors and the distance/angle between them, and there are few discussions about this relationship in drug discovery. Hydrogen bond capacity was quantitatively described [145,146] through an empirical definition. Wang et al. [56] described the hydrogen bond occurrence probability from MD simulation results and, based on their definition, we calculated the hydrogen bond capacity (Hbondcapac.) as defined by:
H b o n d c a p a c . = H b o n d H B D , H B A
where <Hbond> is the average hydrogen bond number during the MD simulation, ∑HBD, while HBA is the sum of the hydrogen bond donor (HBD) and acceptor (HBA) of the molecule.
All steps performed in this study are described in Figure 1.

5. Conclusions

In this study, diverse and robust computer-aided drug discovery methods were sequentially employed against the main protease Mpro of SARS-CoV-2 (SARS-CoV-2 Mpro), aiming to select the most active molecules from the Natural Products Database of the Bahia Semi-Arid region (NatProDB, http://natprodb.uefs.br/, accessed on 7 May 2021). From the “top 10” docked compounds, some common scaffolds were selected (secotrachylobanoic acid, hexahydronaphthalen-1-yl, indolinedione, and benzopyran-4-one) indicating that Mpro could privilege the hydrophobic chains and (poly)aromatic rings. The best interactions with SARS-CoV-2 Mpro revealed that some enzyme sites were accessed, thereby confirming that this method can be employed as a suitable starting method for the identification of novel SARS-CoV-2 Mpro inhibitors.
The pharmacokinetic (ADMET) and toxicological analysis confirmed that nine out of the ten best molecules did not have any issues and could be employed in subsequent in vitro assays. Finally, this study should emphasize that the sequential application of in silico methods (ADMET analysis as well as docking and molecular dynamics) are valuable tools for searching the chemical space and selecting the best ligand structures that could be employed in experimental tests, minimizing the costs of testing many compounds. We also hope that the natural products from the Bahia Semi-arid region chosen here could be safe and effective for treating SARS-CoV-2 infection.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms222111739/s1.

Author Contributions

Conceptualization, R.C.S., C.H.T.P.S., S.S.R.P. and C.B.R.S.; formal analysis, R.C.S., H.F.F., S.S.R.P. and R.S.B.; funding acquisition, R.S.B.; investigation, J.M.C., S.S.R.P. and C.B.R.S.; methodology, H.F.F., N.M.K. and S.S.R.P.; project administration, S.S.R.P.; supervision, S.S.R.P. and C.B.R.S.; visualization, J.M.C., N.M.K., C.H.T.P.S. and R.S.B.; writing—original draft preparation, R.C.S., H.F.F., S.S.R.P. and C.B.R.S.; writing—review and editing, J.M.C., N.M.K., C.H.T.P.S., S.S.R.P. and R.S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Research Dean and Graduate Studies of the Federal University of Pará (PROPESP/UFPA); Brazilian National Council for Scientific and Technological Development (CNPq); Brazilian Coordination for Improvement of Personnel Higher Education (CAPES); and Bahia Research Foundation (FAPESB, grant numbers APP071/2011, JCB-0039/2013, and RED-008/2013).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research project (project IDs: proj651 and proj822) was developed with the help of CENAPAD-SP (National Center for High Performance Computing in São Paulo), UNICAMP/FINEP-MCT project in collaboration with Laboratory of Bioinformatics and Molecular Modeling (LaBiMM) on the Pharmacy College at Federal University of Bahia (UFBA) and the Laboratory of Modeling and Computational Chemistry of the Department of Biological and Health Sciences at the Federal University of Amapá (UNIFAP).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

COVID-19Coronavirus disease 2019
SARS-CoV-2severe acute respiratory syndrome coronavirus 2
ACE2angiotensin-converting enzyme
ADMETabsorption, distribution, metabolism, and excretion/toxicity
MDmolecular dynamics
NatProDBNatural Products Database of the Bahia Semi-Arid region
3CLproSARS-CoV-2 Main Protease
CoVscoronaviruses
PDBProtein Data Bank
PSApolar surface area
%HOAhuman oral absorption in percentage
QPlogPo/wlogarithm of the partition coefficient in 1-octanol/water predicted by QikProp
QPLogBBlogarithm of the blood–brain barrier predicted by QikProp
QPPCacopermeability across Caco-2 cells predicted by QikProp
CNScentral nervous system
RMSFroot-mean square fluctuation
LogKhsalogarithmic human serum albumin-binding predicted by QikProp
EAestimated affinities
Hbondcapac.hydrogen bond capability
ΔEbindingbinding energy
QPPMDCKpermeability across Madin-Darby Canine Kidney cells predicted by QikProp
MWmolecular weight
HTShigh-throughput screening

References

  1. Pillaiyar, T.; Meenakshisundaram, S.; Manickam, M. Recent Discovery and Development of Inhibitors Targeting Coronaviruses. Drug Discov. Today 2020, 25, 668–688. [Google Scholar] [CrossRef] [PubMed]
  2. Guan, W.; Ni, Z.; Hu, Y.; Liang, W.; Ou, C.; He, J.; Liu, L.; Shan, H.; Lei, C.; Hui, D.S.C.; et al. Clinical Characteristics of Coronavirus Disease 2019 in China. N. Engl. J. Med. 2020, 382, 1708–1720. [Google Scholar] [CrossRef] [PubMed]
  3. Wu, F.; Zhao, S.; Yu, B.; Chen, Y.-M.; Wang, W.; Song, Z.-G.; Hu, Y.; Tao, Z.-W.; Tian, J.-H.; Pei, Y.-Y.; et al. A New Coronavirus Associated with Human Respiratory Disease in China. Nature 2008, 579, 265–269. [Google Scholar] [CrossRef] [Green Version]
  4. Zhang, L.; Lin, D.; Sun, X.; Curth, U.; Drosten, C.; Sauerhering, L.; Becker, S.; Rox, K.; Hilgenfeld, R. Crystal Structure of SARS-CoV-2 Main Protease Provides a Basis for Design of Improved α-Ketoamide Inhibitors. Science 2020, 368, 409–412. [Google Scholar] [CrossRef] [Green Version]
  5. Zhou, P.; Yang, X.-L.; Wang, X.-G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.-R.; Zhu, Y.; Li, B.; Huang, C.-L.; et al. A Pneumonia Outbreak Associated with a New Coronavirus of Probable Bat Origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef] [Green Version]
  6. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. N. Engl. J. Med. 2020, 382, 727–733. [Google Scholar] [CrossRef]
  7. Jin, Z.; Du, X.; Xu, Y.; Deng, Y.; Liu, M.; Zhao, Y.; Zhang, B.; Li, X.; Zhang, L.; Peng, C.; et al. Structure of Mpro from SARS-CoV-2 and Discovery of Its Inhibitors. Nature 2020, 582, 289–293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Rosa, S.G.V.; Santos, W.C. Clinical Trials on Drug Repositioning for COVID-19 Treatment. Rev. Panam. Salud Públ. 2020, 44, e40. [Google Scholar] [CrossRef]
  9. Anand, K. Coronavirus Main Proteinase (3CLpro) Structure: Basis for Design of Anti-SARS Drugs. Science 2003, 300, 1763–1767. [Google Scholar] [CrossRef] [Green Version]
  10. Ziebuhr, J. Molecular Biology of Severe Acute Respiratory Syndrome Coronavirus. Curr. Opin. Microbiol. 2004, 7, 412–419. [Google Scholar] [CrossRef]
  11. Gorbalenya, A.E.; Snijder, E.J.; Ziebuhr, J. Virus-Encoded Proteinases and Proteolytic Processing in the Nidovirales. J. Gen. Virol. 2000, 81, 853–879. [Google Scholar] [CrossRef]
  12. Bajaj, A.; Purohit, H.J. Understanding SARS-CoV-2: Genetic Diversity, Transmission and Cure in Human. Indian J. Microbiol. 2020, 60, 398–401. [Google Scholar] [CrossRef] [Green Version]
  13. Sun, J.; He, W.-T.; Wang, L.; Lai, A.; Ji, X.; Zhai, X.; Li, G.; Suchard, M.A.; Tian, J.; Zhou, J.; et al. COVID-19: Epidemiology, Evolution, and Cross-Disciplinary Perspectives. Trends Mol. Med. 2020, 26, 483–495. [Google Scholar] [CrossRef] [Green Version]
  14. Han, Q.; Lin, Q.; Jin, S.; You, L. Coronavirus 2019-NCoV: A Brief Perspective from the Front Line. J. Infect. 2020, 80, 373–377. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Patel, S.K.S.; Lee, J.-K.; Kalia, V.C. Deploying Biomolecules as Anti-COVID-19 Agents. Indian J. Microbiol. 2020, 60, 263–268. [Google Scholar] [CrossRef]
  16. Rishi, P.; Thakur, K.; Vij, S.; Rishi, L.; Singh, A.; Kaur, I.P.; Patel, S.K.S.; Lee, J.-K.; Kalia, V.C. Diet, Gut Microbiota and COVID-19. Indian J. Microbiol. 2020, 60, 420–429. [Google Scholar] [CrossRef] [PubMed]
  17. Thomford, N.; Senthebane, D.; Rowe, A.; Munro, D.; Seele, P.; Maroyi, A.; Dzobo, K. Natural Products for Drug Discovery in the 21st Century: Innovations for Novel Drug Discovery. Int. J. Mol. Sci. 2018, 19, 1578. [Google Scholar] [CrossRef] [Green Version]
  18. Wink, M. Medicinal Plants: A Source of Anti-Parasitic Secondary Metabolites. Molecules 2012, 17, 2771. [Google Scholar] [CrossRef] [Green Version]
  19. Fakhar, Z.; Faramarzi, B.; Pacifico, S.; Faramarzi, S. Anthocyanin Derivatives as Potent Inhibitors of SARS-CoV-2 Main Protease: An in-Silico Perspective of Therapeutic Targets against COVID-19 Pandemic. J. Biomol. Struct. Dyn. 2020, 39, 6171–6183. [Google Scholar] [CrossRef] [PubMed]
  20. Mattio, L.M.; Catinella, G.; Pinto, A.; Dallavalle, S. Natural and Nature-Inspired Stilbenoids as Antiviral Agents. Eur. J. Med. Chem. 2020, 202, 112541. [Google Scholar] [CrossRef] [PubMed]
  21. Newman, D.J.; Cragg, G.M. Natural Products As Sources of New Drugs over the 30 Years from 1981 to 2010. J. Nat. Prod. 2012, 75, 311–335. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Olubiyi, O.O.; Olagunju, M.; Keutmann, M.; Loschwitz, J.; Strodel, B. High Throughput Virtual Screening to Discover Inhibitors of the Main Protease of the Coronavirus SARS-CoV-2. Molecules 2020, 25, 3193. [Google Scholar] [CrossRef]
  23. Sillapachaiyaporn, C.; Chuchawankul, S. HIV-1 Protease and Reverse Transcriptase Inhibition by Tiger Milk Mushroom (Lignosus Rhinocerus) Sclerotium Extracts: In Vitro and in Silico Studies. J. Tradit. Complement. Med. 2020, 10, 396–404. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, Z.; Yang, L. Turning the Tide: Natural Products and Natural-Product-Inspired Chemicals as Potential Counters to SARS-CoV-2 Infection. Front. Pharmacol. 2019, 11, 1013. [Google Scholar] [CrossRef] [PubMed]
  25. Lewinsohn, T.M.; Prado, P.I. Biodiversity of Brazil: A synthesis of the current state of knowledge. In Evaluation of the State of Knowledge on Biological Diversity in Brazil; Baumgarten, L.C., de Andrade, L.A.Z., Cariello, M.O., Eds.; CID Ambiental: Brasilia, Brazil, 2003; pp. 139–144. [Google Scholar]
  26. Lucchese, A.M.; Vale, A.E. Plantas Da Caatinga: Perfil Botânico, Fitoquímica e Atividade Biológica; Associação Plantas do Nordeste: Recife, Brazil, 2006; Volume 4. [Google Scholar]
  27. Pires, D.E.V.; Blundell, T.L.; Ascher, D.B. PkCSM: Predicting Small-Molecule Pharmacokinetic and Toxicity Properties Using Graph-Based Signatures. J. Med. Chem. 2015, 58, 4066–4072. [Google Scholar] [CrossRef] [PubMed]
  28. Da Paixão, V.G.; Pita, S.S.D.R. In Silico Identification and Evaluation of New Trypanosoma Cruzi Trypanothione Reductase (TcTR) Inhibitors Obtained from Natural Products Database of the Bahia Semi-Arid Region (NatProDB). Comput. Biol. Chem. 2019, 79, 36–47. [Google Scholar] [CrossRef]
  29. Macmillan, D.S.; Chilton, M.L. A Defined Approach for Predicting Skin Sensitisation Hazard and Potency Based on the Guided Integration of in Silico, in Chemico and in Vitro Data Using Exclusion Criteria. Regul. Toxicol. Pharmacol. 2019, 101, 35–47. [Google Scholar] [CrossRef]
  30. Silva, R.; Poiani, J.; Ramos, R.; Costa, J.; Silva, C.; Brasil, D.; Santos, C. Ligand- and Structure-Based Virtual Screening from 16-((Diisobutylamino)Methyl)-6α-Hydroxyivouacapane-7β,17β-Lactone a Compound with Potential Anti-Prostate Cancer Activity. J. Serbian Chem. Soc. 2019, 84, 153–174. [Google Scholar] [CrossRef] [Green Version]
  31. Duffy, E.M.; Jorgensen, W.L. Prediction of Properties from Simulations: Free Energies of Solvation in Hexadecane, Octanol, and Water. J. Am. Chem. Soc. 2000, 122, 2878–2888. [Google Scholar] [CrossRef]
  32. Jorgensen, W.L.; Duffy, E.M. Prediction of Drug Solubility from Monte Carlo Simulations. Bioorg. Med. Chem. Lett. 2000, 10, 1155–1158. [Google Scholar] [CrossRef]
  33. Silva, R.C.; Ferreira, I.M.; Federico, L.B.; Hage-Melim, L.I.d.S.; Macêdo, W.J.C.; Porto, A.L.M.; Taft, C.A.; Silva, C.H.T.P.; dos Santos, C.B.R. ADME/Tox Study and Molecular Dynamics Simulations Applied in the Design of New Potential GABA-AT Inhibitors. In Functional Properties of Advanced Engineering Materials and Biomolecules; Springer: Cham, Switzerland, 2021. [Google Scholar]
  34. Cruz, J.V.; Giuliatti, S.; Alves, L.B.; Silva, R.C.; Ferreira, E.F.B.; Kimani, N.M.; Silva, C.H.T.P.; Souza, J.S.N.d.; Espejo-Román, J.M.; Santos, C.B.R. Identification of Novel Potential Cyclooxygenase-2 Inhibitors Using Ligand- and Structure-Based Virtual Screening Approaches. J. Biomol. Struct. Dyn. 2021, 1–23. [Google Scholar] [CrossRef] [PubMed]
  35. Sutter, A.; Amberg, A.; Boyer, S.; Brigo, A.; Contrera, J.F.; Custer, L.L.; Dobo, K.L.; Gervais, V.; Glowienke, S.; van Gompel, J.; et al. Use of in Silico Systems and Expert Knowledge for Structure-Based Assessment of Potentially Mutagenic Impurities. Regul. Toxicol. Pharmacol. 2013, 67, 39–52. [Google Scholar] [CrossRef]
  36. Ertl, P.; Rohde, B.; Selzer, P. Fast Calculation of Molecular Polar Surface Area as a Sum of Fragment-Based Contributions and Its Application to the Prediction of Drug Transport Properties. J. Med. Chem. 2000, 43, 3714–3717. [Google Scholar] [CrossRef] [PubMed]
  37. Palm, K.; Stenberg, P.; Luthman, K.; Artursson1, P. Polar Molecular Surface Properties Predict the Intestinal Absorption of Drugs in Humans. Pharm. Res. 1997, 14, 568–571. [Google Scholar] [CrossRef]
  38. Johnson-Davis, K.L.; Dasgupta, A. Special Issues in Therapeutic Drug Monitoring in Patients. In Clinical Challenges in Therapeutic Drug Monitoring; Elsevier: Amsterdam, The Netherlands, 2016. [Google Scholar]
  39. Braggio, S.; Corsi, M.; Feriani, A.; Fontana, S.; Marocchio, L.; Virginio, C. CHAPTER 15. Discovery Toxicology In Lead Optimisation. In The Handbook of Medicinal Chemistry; Royal Society of Chemistry: Cambridge, UK, 2015. [Google Scholar]
  40. Hanai, T.; Koseki, A.; Yoshikawa, R.; Ueno, M.; Kinoshita, T.; Homma, H. Prediction of Human Serum Albumin–Drug Binding Affinity without Albumin. Anal. Chim. Acta 2002, 454, 101–108. [Google Scholar] [CrossRef]
  41. Jorgensen, W.L.; Duffy, E.M. Prediction of Drug Solubility from Structure. Adv. drug Deliv. Rev. 2002, 54, 355–366. [Google Scholar] [CrossRef]
  42. Noctor, T.A.G.; Wainer, I.W.; Hage, D.S. Allosteric and Competitive Displacement of Drugs from Human Serum Albumin by Octanoic Acid, as Revealed by High-Performance Liquid Affinity Chromatography, on a Human Serum Albumin-Based Stationary Phase. J. Chromatogr. B Biomed. Sci. Appl. 1992, 577, 305–315. [Google Scholar] [CrossRef]
  43. Urien, S.; Bre, F.; Breillout, F.; Bastian, G.; Krikorian, A.; Tillement, J.P. Vinorelbine High-Affinity Binding to Human Platelets and Lymphocytes: Distribution in Human Blood. Cancer Chemother. Pharmacol. 1993, 32, 231–234. [Google Scholar] [CrossRef]
  44. Haverkamp, W. The Potential for QT Prolongation and Proarrhythmia by Non-Antiarrhythmic Drugs: Clinical and Regulatory Implications. Report on a Policy Conference of the European Society of Cardiology. Eur. Heart J. 2000, 21, 1216–1231. [Google Scholar] [CrossRef]
  45. Gautret, P.; Lagier, J.-C.; Parola, P.; Hoang, V.T.; Meddeb, L.; Mailhe, M.; Doudier, B.; Courjon, J.; Giordanengo, V.; Vieira, V.E.; et al. Hydroxychloroquine and Azithromycin as a Treatment of COVID-19: Results of an Open-Label Non-Randomized Clinical Trial. Int. J. Antimicrob. Agents 2020, 56, 105949. [Google Scholar] [CrossRef] [PubMed]
  46. Malviya, A. Ventricular Arrhythmia Risk Due to Chloroquine/Hydroxychloroquine Treatment for COVID-19: Should It Be Given. Indian Heart J. 2020, 72, 131–132. [Google Scholar] [CrossRef]
  47. Hussein, H.A.; Borrel, A.; Geneix, C.; Petitjean, M.; Regad, L.; Camproux, A.-C. PockDrug-Server: A New Web Server for Predicting Pocket Druggability on Holo and Apo Proteins. Nucleic Acids Res. 2015, 43, W436–W442. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Milletti, F.; Vulpetti, A. Predicting Polypharmacology by Binding Site Similarity: From Kinases to the Protein Universe. J. Chem. Inf. Model. 2010, 50, 1418–1431. [Google Scholar] [CrossRef] [PubMed]
  49. Kozakov, D.; Grove, L.E.; Hall, D.R.; Bohnuud, T.; Mottarella, S.E.; Luo, L.; Xia, B.; Beglov, D.; Vajda, S. The FTMap Family of Web Servers for Determining and Characterizing Ligand-Binding Hot Spots of Proteins. Nat. Protoc. 2015, 10, 733–755. [Google Scholar] [CrossRef] [Green Version]
  50. Tan, J.; Verschueren, K.H.G.; Anand, K.; Shen, J.; Yang, M.; Xu, Y.; Rao, Z.; Bigalke, J.; Heisen, B.; Mesters, J.R.; et al. PH-Dependent Conformational Flexibility of the SARS-CoV Main Proteinase (Mpro) Dimer: Molecular Dynamics Simulations and Multiple X-Ray Structure Analyses. J. Mol. Biol. 2005, 354, 25–40. [Google Scholar] [CrossRef]
  51. Yang, H.; Xie, W.; Xue, X.; Yang, K.; Ma, J.; Liang, W.; Zhao, Q.; Zhou, Z.; Pei, D.; Ziebuhr, J.; et al. Design of Wide-Spectrum Inhibitors Targeting Coronavirus Main Proteases. PLoS Biol. 2005, 3, e324. [Google Scholar] [CrossRef]
  52. Volkamer, A.; Griewel, A.; Grombacher, T.; Rarey, M. Analyzing the Topology of Active Sites: On the Prediction of Pockets and Subpockets. J. Chem. Inf. Model. 2010, 50, 2041–2052. [Google Scholar] [CrossRef]
  53. Volkamer, A.; Kuhn, D.; Rippmann, F.; Rarey, M. DoGSiteScorer: A Web Server for Automatic Binding Site Prediction, Analysis and Druggability Assessment. Bioinformatics 2012, 28, 2074–2075. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Kozakov, D.; Hall, D.R.; Napoleon, R.L.; Yueh, C.; Whitty, A.; Vajda, S. New Frontiers in Druggability. J. Med. Chem. 2015, 58, 9063–9088. [Google Scholar] [CrossRef] [PubMed]
  55. Schrödinger LLC. The Pymol Molecular Graphics System. 2013. Available online: http://www.pymol.org (accessed on 18 October 2021).
  56. Wang, T.; Xie, H.-B.; Song, Z.; Niu, J.; Chen, D.-L.; Xia, D.; Chen, J. Role of Hydrogen Bond Capacity of Solvents in Reactions of Amines with CO2: A Computational Study. J. Environ. Sci. 2020, 91, 271–278. [Google Scholar] [CrossRef]
  57. Lagorce, D.; Sperandio, O.; Galons, H.; Miteva, M.A.; Villoutreix, B.O. FAF-Drugs2: Free ADME/Tox Filtering Tool to Assist Drug Discovery and Chemical Biology Projects. BMC Bioinform. 2008, 9, 396. [Google Scholar] [CrossRef] [Green Version]
  58. Lagorce, D.; Oliveira, N.; Miteva, M.A.; Villoutreix, B.O. Pan-Assay Interference Compounds (PAINS) That May Not Be Too Painful for Chemical Biology Projects. Drug Discov. Today 2017, 22, 1131–1133. [Google Scholar] [CrossRef] [PubMed]
  59. Lagorce, D.; Sperandio, O.; Baell, J.B.; Miteva, M.A.; Villoutreix, B.O. FAF-Drugs3: A Web Server for Compound Property Calculation and Chemical Library Design. Nucleic Acids Res. 2015, 43, W200–W207. [Google Scholar] [CrossRef] [PubMed]
  60. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and Computational Approaches to Estimate Solubility and Permeability in Drug Discovery and Development Settings. Adv. Drug Deliv. Rev. 2001, 46, 3–26. [Google Scholar] [CrossRef]
  61. Reulecke, I.; Lange, G.; Albrecht, J.; Klein, R.; Rarey, M. Towards an Integrated Description of Hydrogen Bonding and Dehydration: Decreasing False Positives in Virtual Screening with the HYDE Scoring Function. ChemMedChem 2008, 3, 885–897. [Google Scholar] [CrossRef]
  62. Schneider, N.; Hindle, S.; Lange, G.; Klein, R.; Albrecht, J.; Briem, H.; Beyer, K.; Claußen, H.; Gastreich, M.; Lemmen, C.; et al. Substantial Improvements in Large-Scale Redocking and Screening Using the Novel HYDE Scoring Function. J. Comput. Aided. Mol. Des. 2012, 26, 701–723. [Google Scholar] [CrossRef]
  63. Schneider, N.; Lange, G.; Hindle, S.; Klein, R.; Rarey, M. A Consistent Description of HYdrogen Bond and DEhydration Energies in Protein–Ligand Complexes: Methods behind the HYDE Scoring Function. J. Comput. Aided. Mol. Des. 2013, 27, 15–29. [Google Scholar] [CrossRef] [PubMed]
  64. Vajda, S.; Beglov, D.; Wakefield, A.E.; Egbert, M.; Whitty, A. Cryptic Binding Sites on Proteins: Definition, Detection, and Druggability. Curr. Opin. Chem. Biol. 2018, 44, 1–8. [Google Scholar] [CrossRef]
  65. Barrett, A.J.; Kembhavi, A.A.; Brown, M.A.; Kirschke, H.; Knight, C.G.; Tamai, M.; Hanada, K. L-Trans-Epoxysuccinyl-Leucylamido(4-Guanidino)Butane (E-64) and Its Analogues as Inhibitors of Cysteine Proteinases Including Cathepsins B, H and L. Biochem. J. 1982, 201, 189–198. [Google Scholar] [CrossRef] [Green Version]
  66. Matsumoto, K.; Mizoue, K.; Kitamura, K.; Tse, W.-C.; Huber, C.P.; Ishida, T. Structural Basis of Inhibition of Cysteine Proteases by E-64 and Its Derivatives. Biopolymers 1999, 51, 99–107. [Google Scholar] [CrossRef]
  67. Salentin, S.; Schreiber, S.; Haupt, V.J.; Adasme, M.F.; Schroeder, M. PLIP: Fully Automated Protein–Ligand Interaction Profiler. Nucleic Acids Res. 2015, 43, W443–W447. [Google Scholar] [CrossRef]
  68. Kabsch, W.; Sander, C. Dictionary of Protein Secondary Structure: Pattern Recognition of Hydrogen-Bonded and Geometrical Features. Biopolymers 1983, 22, 2577–2637. [Google Scholar] [CrossRef] [PubMed]
  69. Sander, C.; Schneider, R. Database of Homology-Derived Protein Structures and the Structural Meaning of Sequence Alignment. Proteins Struct. Funct. Genet. 1991, 9, 56–68. [Google Scholar] [CrossRef]
  70. Touw, W.G.; Baakman, C.; Black, J.; te Beek, T.A.H.; Krieger, E.; Joosten, R.P.; Vriend, G. A Series of PDB-Related Databanks for Everyday Needs. Nucleic Acids Res. 2015, 43, D364–D368. [Google Scholar] [CrossRef]
  71. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers. SoftwareX 2015, 1–2, 19–25. [Google Scholar] [CrossRef] [Green Version]
  72. Berendsen, H.J.C.; van der Spoel, D.; van Drunen, R. GROMACS: A Message-Passing Parallel Molecular Dynamics Implementation. Comput. Phys. Commun. 1995, 91, 43–56. [Google Scholar] [CrossRef]
  73. Lindahl, E.; Hess, B.; van der Spoel, D. GROMACS 3.0: A Package for Molecular Simulation and Trajectory Analysis. J. Mol. Model. 2001, 7, 306–317. [Google Scholar] [CrossRef]
  74. Pronk, S.; Páll, S.; Schulz, R.; Larsson, P.; Bjelkmar, P.; Apostolov, R.; Shirts, M.R.; Smith, J.C.; Kasson, P.M.; van der Spoel, D.; et al. GROMACS 4.5: A High-Throughput and Highly Parallel Open Source Molecular Simulation Toolkit. Bioinformatics 2013, 29, 845–854. [Google Scholar] [CrossRef]
  75. Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J.C. GROMACS: Fast, Flexible, and Free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef]
  76. Bekker, H.; Berendsen, H.J.C.; Dijkstra, E.J.; Achterop, S.; Van Drunen, R.; Van der Spoel, D.; Berendsen, H.; Van der Westhuyzen, D.R. Gromacs-a Parallel Computer for Molecular-Dynamics Simulations. In Proceedings of the 4th International Conference on Computational Physics (PC 92); DeGroot, R.A., Nadrchal, J., Eds.; World Scientific Publishing: Singapore, 1993; pp. 252–256. [Google Scholar]
  77. Van der Spoel, D.; van Maaren, P.J.; Larsson, P.; Tîmneanu, N. Thermodynamics of Hydrogen Bonding in Hydrophilic and Hydrophobic Media. J. Phys. Chem. B 2006, 110, 4393–4398. [Google Scholar] [CrossRef]
  78. Gomes, D.E.B.; Silva, A.W.; Lins, R.D.; Pascutti, P.G.; Soares, T.A. HbMap2Grace 2002. Available online: http://lmdm.biof.ufrj.br/software/hbmap2grace/index.html (accessed on 18 October 2021).
  79. Gomes, D.E.B.; Sousa, G.L.S.C.; Silva, A.W.S.D.; Pascutti, P.G. SurfinMD 2012. Available online: http://lmdm.biof.ufrj.br/software/surfinmd/index.html (accessed on 18 October 2021).
  80. Patil, R.; Das, S.; Stanley, A.; Yadav, L.; Sudhakar, A.; Varma, A.K. Optimized Hydrophobic Interactions and Hydrogen Bonding at the Target-Ligand Interface Leads the Pathways of Drug-Designing. PLoS ONE 2010, 5, e12029. [Google Scholar] [CrossRef]
  81. Sharma, P.; Vijayan, V.; Pant, P.; Sharma, M.; Vikram, N.; Kaur, P.; Singh, T.P.; Sharma, S. Identification of Potential Drug Candidates to Combat COVID-19: A Structural Study Using the Main Protease (Mpro) of SARS-CoV-2. J. Biomol. Struct. Dyn. 2020, 6649–6659. [Google Scholar] [CrossRef] [PubMed]
  82. Le Guilloux, V.; Schmidtke, P.; Tuffery, P. Fpocket: An Open Source Platform for Ligand Pocket Detection. BMC Bioinformatics 2009, 10, 168. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Schmidtke, P.; Le Guilloux, V.; Maupetit, J.; Tuffery, P. Fpocket: Online Tools for Protein Ensemble Pocket Detection and Tracking. Nucleic Acids Res. 2010, 38, W582–W589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Ton, A.; Gentile, F.; Hsing, M.; Ban, F.; Cherkasov, A. Rapid Identification of Potential Inhibitors of SARS-CoV-2 Main Protease by Deep Docking of 1.3 Billion Compounds. Mol. Inform. 2020, 39, 2000028. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Baker, N.A.; Sept, D.; Joseph, S.; Holst, M.J.; McCammon, J.A. Electrostatics of Nanosystems: Application to Microtubules and the Ribosome. Proc. Natl. Acad. Sci. USA 2001, 98, 10037–10041. [Google Scholar] [CrossRef] [Green Version]
  86. Kumari, R.; Kumar, R.; Lynn, A. G_mmpbsa —A GROMACS Tool for High-Throughput MM-PBSA Calculations. J. Chem. Inf. Model. 2014, 54, 1951–1962. [Google Scholar] [CrossRef]
  87. Gbadamosi, I.T. Stay Safe: Helpful Herbal Remedies in Covid-19 Infection. African J. Biomed. Res. 2020, 23, 131–133. [Google Scholar]
  88. Chou, C.-Y.; Chang, H.-C.; Hsu, W.-C.; Lin, T.-Z.; Lin, C.-H.; Chang, G.-G. Quaternary Structure of the Severe Acute Respiratory Syndrome (SARS) Coronavirus Main Protease. Biochemistry 2004, 43, 14958–14970. [Google Scholar] [CrossRef]
  89. Zhong, N.; Zhang, S.; Zou, P.; Chen, J.; Kang, X.; Li, Z.; Liang, C.; Jin, C.; Xia, B. Without Its N-Finger, the Main Protease of Severe Acute Respiratory Syndrome Coronavirus Can Form a Novel Dimer through Its C-Terminal Domain. J. Virol. 2008, 82, 4227–4234. [Google Scholar] [CrossRef] [Green Version]
  90. Jansen, R.; Irschik, H.; Huch, V.; Schummer, D.; Steinmetz, H.; Bock, M.; Schmidt, T.; Kirschning, A.; Müller, R. Carolacton—A Macrolide Ketocarbonic Acid That Reduces Biofilm Formation by the Caries- and Endocarditis-Associated Bacterium Streptococcus Mutans. Eur. J. Org. Chem. 2010, 2010, 1284–1289. [Google Scholar] [CrossRef]
  91. Anderson, D.E.; Cui, J.; Ye, Q.; Huang, B.; Zu, W.; Gong, J.; Liu, W.; Young Kim, S.; Guo Yan, B.; Sigmundsson, K.; et al. Orthogonal Genome-Wide Screenings in Bat Cells Identify MTHFD1 as a Target of 1 Broad Antiviral Therapy 2 3. bioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  92. Wiegrebe, W.; Kramer, W.J.; Shamma, M. The Emetine Alkaloids. J. Nat. Prod. 1984, 47, 397–408. [Google Scholar] [CrossRef]
  93. Akinboye, E.S.; Rosen, M.D.; Denmeade, S.R.; Kwabi-Addo, B.; Bakare, O. Design, Synthesis, and Evaluation of PH-Dependent Hydrolyzable Emetine Analogues as Treatment for Prostate Cancer. J. Med. Chem. 2012, 55, 7450–7459. [Google Scholar] [CrossRef] [PubMed]
  94. Shen, L.; Niu, J.; Wang, C.; Huang, B.; Wang, W.; Zhu, N.; Deng, Y.; Wang, H.; Ye, F.; Cen, S.; et al. High-Throughput Screening and Identification of Potent Broad-Spectrum Inhibitors of Coronaviruses. J. Virol. 2019, 93, 93. [Google Scholar] [CrossRef] [Green Version]
  95. Dyall, J.; Coleman, C.M.; Hart, B.J.; Venkataraman, T.; Holbrook, M.R.; Kindrachuk, J.; Johnson, R.F.; Olinger, G.G.; Jahrling, P.B.; Laidlaw, M.; et al. Repurposing of Clinically Developed Drugs for Treatment of Middle East Respiratory Syndrome Coronavirus Infection. Antimicrob. Agents Chemother. 2014, 58, 4885–4893. [Google Scholar] [CrossRef] [Green Version]
  96. Gentile, D.; Patamia, V.; Scala, A.; Sciortino, M.T.; Piperno, A.; Rescifina, A. Putative Inhibitors of SARS-CoV-2 Main Protease from A Library of Marine Natural Products: A Virtual Screening and Molecular Modeling Study. Mar. Drugs 2020, 18, 225. [Google Scholar] [CrossRef] [Green Version]
  97. Borges, R.S.; Palheta, I.C.; Ota, S.S.B.; Morais, R.B.; Barros, V.A.; Ramos, R.S.; Silva, R.C.; da Silva Costa, J.; Silva, C.H.T.P.; Campos, J.M.; et al. Toward of Safer Phenylbutazone Derivatives by Exploration of Toxicity Mechanism. Molecules 2018, 23, 143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  98. Cunha, E.L.; Santos, C.F.; Braga, F.S.; Costa, J.S.; Silva, R.C.; Favacho, H.A.S.; Hage-Melim, L.I.S.; Carvalho, J.C.T.; da Silva, C.H.T.P.; Santos, C.B.R. Computational Investigation of Antifungal Compounds Using Molecular Modeling and Prediction of ADME/Tox Properties. J. Comput. Theor. Nanosci. 2015, 12, 3682–3691. [Google Scholar] [CrossRef]
  99. Stenberg, P.; Norinder, U.; Luthman, K.; Artursson, P. Experimental and Computational Screening Models for the Prediction of Intestinal Drug Absorption. J. Med. Chem. 2001, 44, 1927–1937. [Google Scholar] [CrossRef] [PubMed]
  100. Shityakov, S.; Förster, C. Advances and Applications in Bioinformatics and Chemistry Dovepress in Silico Predictive Model to Determine Vector-Mediated Transport Properties for the Blood-Brain Barrier Choline Transporter. Adv. Appl. Bioinforma. Chem. 2014, 7, 7–23. [Google Scholar] [CrossRef] [Green Version]
  101. Yang, X.; Yu, Y.; Xu, J.; Shu, H.; Xia, J.; Liu, H.; Wu, Y.; Zhang, L.; Yu, Z.; Fang, M.; et al. Clinical Course and Outcomes of Critically Ill Patients with SARS-CoV-2 Pneumonia in Wuhan, China: A Single-Centered, Retrospective, Observational Study. Lancet Respir. Med. 2020, 8, 475–481. [Google Scholar] [CrossRef] [Green Version]
  102. Platt, M.P.; Bolding, K.A.; Wayne, C.R.; Chaudhry, S.; Cutforth, T.; Franks, K.M.; Agalliu, D. Th17 Lymphocytes Drive Vascular and Neuronal Deficits in a Mouse Model of Postinfectious Autoimmune Encephalitis. Proc. Natl. Acad. Sci. USA 2020, 117, 6708–6716. [Google Scholar] [CrossRef] [PubMed]
  103. Han, H.; Yang, L.; Liu, R.; Liu, F.; Wu, K.; Li, J.; Liu, X.; Zhu, C. Prominent Changes in Blood Coagulation of Patients with SARS-CoV-2 Infection. Clin. Chem. Lab. Med. 2020, 58, 1116–1120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  104. Mao, L.; Jin, H.; Wang, M.; Hu, Y.; Chen, S.; He, Q.; Chang, J.; Hong, C.; Zhou, Y.; Wang, D.; et al. Neurologic Manifestations of Hospitalized Patients With Coronavirus Disease 2019 in Wuhan, China. JAMA Neurol. 2020, 77, 683. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Mannan Baig, A.; Khaleeq, A.; Ali, U.; Syeda, H. Evidence of the COVID-19 Virus Targeting the CNS: Tissue Distribution, Host−Virus Interaction, and Proposed Neurotropic Mechanisms. ACS Chem. Neurosci. 2020, 11, 995–998. [Google Scholar] [CrossRef] [Green Version]
  106. De Felice, F.G.; Tovar-Moll, F.; Moll, J.; Munoz, D.P.; Ferreira, S.T. Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) and the Central Nervous System. Trends Neurosci. 2020, 43, 355–357. [Google Scholar] [CrossRef] [PubMed]
  107. Rojas, H.; Ritter, C.; Pizzol, F.D. Mechanisms of Dysfunction of the Blood-Brain Barrier in Critically Ill Patients: Emphasis on the Role of Matrix Metalloproteinases. Rev. Bras. Ter. Intensiva 2011, 23, 222–227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  108. Wang, J.; Hou, T. Drug and Drug Candidate Building Block Analysis. J. Chem. Inf. Model. 2010, 50, 55–67. [Google Scholar] [CrossRef]
  109. Ren, X.; Glende, J.; Al-Falah, M.; de Vries, V.; Schwegmann-Wessels, C.; Qu, X.; Tan, L.; Tschernig, T.; Deng, H.; Naim, H.Y.; et al. Analysis of ACE2 in Polarized Epithelial Cells: Surface Expression and Function as Receptor for Severe Acute Respiratory Syndrome-Associated Coronavirus. J. Gen. Virol. 2006, 87, 1691–1695. [Google Scholar] [CrossRef]
  110. Uemura, K.; Sasaki, M.; Sanaki, T.; Toba, S.; Takahashi, Y.; Orba, Y.; Hall, W.W.; Maenaka, K.; Sawa, H.; Sato, A. MRC5 Cells Engineered to Express ACE2 Serve as a Model System for the Discovery of Antivirals Targeting SARS-CoV-2. Sci. Rep. 2021, 11, 5376. [Google Scholar] [CrossRef]
  111. Warren, T.K.; Jordan, R.; Lo, M.K.; Ray, A.S.; Mackman, R.L.; Soloveva, V.; Siegel, D.; Perron, M.; Bannister, R.; Hui, H.C.; et al. Therapeutic Efficacy of the Small Molecule GS-5734 against Ebola Virus in Rhesus Monkeys. Nature 2016, 531, 1–9. [Google Scholar] [CrossRef] [PubMed]
  112. Ellinger, B.; Bojkova, D.; Zaliani, A.; Cinatl, J.; Claussen, C.; Westhaus, S.; Keminer, O.; Reinshagen, J.; Kuzikov, M.; Wolf, M.; et al. A SARS-CoV-2 Cytopathicity Dataset Generated by High-Content Screening of a Large Drug Repurposing Collection. Sci. Data 2021, 8, 1–10. [Google Scholar] [CrossRef] [PubMed]
  113. Hansch, C.; Fujita, T. P -σ-π Analysis. A Method for the Correlation of Biological Activity and Chemical Structure. J. Am. Chem. Soc. 1964, 86, 1616–1626. [Google Scholar] [CrossRef]
  114. Box, K.; Comer, J. Using Measured PKa, LogP and Solubility to Investigate Supersaturation and Predict BCS Class. Curr. Drug Metab. 2008, 9, 869–878. [Google Scholar] [CrossRef]
  115. Khalifa, I.; Zhu, W.; Hamed, H.; Mohammed, H.; Dutta, K.; Chunmei, L.; Li, C. Tannins Inhibit SARS-CoV-2 through Binding with Catalytic Dyad Residues of 3CL pro: An in Silico Approach with 19 Structural Different Hydrolysable Tannins Practical Applications. J. Food Biochem. 2020, 44, 13432. [Google Scholar] [CrossRef]
  116. Tahir ul Qamar, M.; Alqahtani, S.M.; Alamri, M.A.; Chen, L.-L. Structural Basis of SARS-CoV-2 3CLpro and Anti-COVID-19 Drug Discovery from Medicinal Plants. J. Pharm. Anal. 2020, 10, 313–319. [Google Scholar] [CrossRef]
  117. Ferreira, J.C.; Fadl, S.; Villanueva, A.J.; Rabeh, W.M. Catalytic Dyad Residues His41 and Cys145 Impact the Catalytic Activity and Overall Conformational Fold of the Main SARS-CoV-2 Protease 3-Chymotrypsin-Like Protease. Front. Chem. 2021, 9, 491. [Google Scholar] [CrossRef]
  118. Alamri, M.A.; Tahir ul Qamar, M.; Mirza, M.U.; Bhadane, R.; Alqahtani, S.M.; Muneer, I.; Froeyen, M.; Salo-Ahen, O.M.H. Pharmacoinformatics and Molecular Dynamics Simulation Studies Reveal Potential Covalent and FDA-Approved Inhibitors of SARS-CoV-2 Main Protease 3CL Pro. J. Biomol. Struct. Dyn. 2021, 39, 4936–4948. [Google Scholar] [CrossRef] [PubMed]
  119. Xu, T.; Zheng, W.; Huang, R. High-throughput Screening Assays for SARS-CoV-2 Drug Development: Current Status and Future Directions. Drug Discov. Today 2021, 26, 2439–2444. [Google Scholar] [CrossRef]
  120. Douangamath, A.; Fearon, D.; Gehrtz, P.; Krojer, T.; Lukacik, P.; Owen, C.D.; Resnick, E.; Strain-Damerell, C.; Aimon, A.; Ábrányi-Balogh, P.; et al. Crystallographic and Electrophilic Fragment Screening of the SARS-CoV-2 Main Protease. Nat. Commun. 2020, 11, 1–11. [Google Scholar] [CrossRef]
  121. Volkamer, A.; Kuhn, D.; Grombacher, T.; Rippmann, F.; Rarey, M. Combining Global and Local Measures for Structure-Based Druggability Predictions. J. Chem. Inf. Model. 2012, 52, 360–372. [Google Scholar] [CrossRef] [PubMed]
  122. Borrel, A.; Regad, L.; Xhaard, H.; Petitjean, M.; Camproux, A.-C. PockDrug: A Model for Predicting Pocket Druggability That Overcomes Pocket Estimation Uncertainties. J. Chem. Inf. Model. 2015, 55, 882–895. [Google Scholar] [CrossRef]
  123. Kyte, J.; Doolittle, R.F. A Simple Method for Displaying the Hydropathic Character of a Protein. J. Mol. Biol. 1982, 157, 105–132. [Google Scholar] [CrossRef] [Green Version]
  124. Ngan, C.H.; Bohnuud, T.; Mottarella, S.E.; Beglov, D.; Villar, E.A.; Hall, D.R.; Kozakov, D.; Vajda, S. FTMAP: Extended Protein Mapping with User-Selected Probe Molecules. Nucleic Acids Res. 2012, 40, W271–W275. [Google Scholar] [CrossRef] [PubMed]
  125. Brenke, R.; Kozakov, D.; Chuang, G.-Y.; Beglov, D.; Hall, D.; Landon, M.R.; Mattos, C.; Vajda, S. Fragment-Based Identification of Druggable ‘Hot Spots’ of Proteins Using Fourier Domain Correlation Techniques. Bioinformatics 2009, 25, 621–627. [Google Scholar] [CrossRef] [Green Version]
  126. Kozakov, D.; Hall, D.R.; Chuang, G.-Y.; Cencic, R.; Brenke, R.; Grove, L.E.; Beglov, D.; Pelletier, J.; Whitty, A.; Vajda, S. Structural Conservation of Druggable Hot Spots in Protein-Protein Interfaces. Proc. Natl. Acad. Sci. USA 2011, 108, 13528–13533. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  127. Schärfer, C.; Schulz-Gasch, T.; Ehrlich, H.-C.; Guba, W.; Rarey, M.; Stahl, M. Torsion Angle Preferences in Druglike Chemical Space: A Comprehensive Guide. J. Med. Chem. 2013, 56, 2016–2028. [Google Scholar] [CrossRef]
  128. Brethon, A.; Chantalat, L.; Christin, O.; Clary, L.; Fournier, J.-F.; Gastreich, M.; Harris, C.S.; Isabet, T.; Pascau, J.; Thoreau, E.; et al. New Caspase-1 Inhibitor by Scaffold Hopping into Bio-Inspired 3D-Fragment Space. Bioorg. Med. Chem. Lett. 2017, 27, 5373–5377. [Google Scholar] [CrossRef]
  129. Hopkins, A.L.; Groom, C.R.; Alex, A. Ligand Efficiency: A Useful Metric for Lead Selection. Drug Discov. Today 2004, 9, 430–431. [Google Scholar] [CrossRef]
  130. Macêdo, W.J.C.; Silva, R.C.; Taft, C.A.; Silva, C.H.T.P.; Rodriguez, A.F.R.; Campos, J.M.; dos Santos, C.B.R. High-Throughput-Based Virtual Screening via Molecular Docking for Oxidative Stress Mediated by ROS Enzyme. In Functional Properties of Advanced Engineering Materials and Biomolecules; Engineering Materials; La Porta, F.A., Taft, C.A., Eds.; Springer: Cham, Switzerland; Hillerød, Denmark, 2021; pp. 489–513. [Google Scholar] [CrossRef]
  131. Schrödinger QikProp: Rapid ADME Predictions of Drug Candidates. Available online: https://www.schrodinger.com/products/qikprop (accessed on 18 October 2021).
  132. Hou, T.; Wang, J. Structure—ADME Relationship: Still a Long Way to Go? Expert Opin. Drug Metab. Toxicol. 2008, 4, 759–770. [Google Scholar] [CrossRef]
  133. Greene, N.; Judson, P.N.; Langowski, J.J.; Marchant, C.A. Knowledge-Based Expert Systems for Toxicity and Metabolism Prediction: DEREK, StAR and METEOR. SAR QSAR Environ. Res. 1999, 10, 299–314. [Google Scholar] [CrossRef] [PubMed]
  134. Hess, B.; Kutzner, C.; van der Spoel, D.; Lindahl, E. GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J. Chem. Theory Comput. 2008, 4, 435–447. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  135. Schmid, N.; Eichenberger, A.P.; Choutko, A.; Riniker, S.; Winger, M.; Mark, A.E.; van Gunsteren, W.F. Definition and Testing of the GROMOS Force-Field Versions 54A7 and 54B7. Eur. Biophys. J. 2011, 40, 843–856. [Google Scholar] [CrossRef]
  136. Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald: An N log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98, 10089–10092. [Google Scholar] [CrossRef] [Green Version]
  137. Berendsen, H.J.C.; Grigera, J.R.; Straatsma, T.P. The Missing Term in Effective Pair Potentials. J. Phys. Chem. 1987, 91, 6269–6271. [Google Scholar] [CrossRef]
  138. Bussi, G.; Donadio, D.; Parrinello, M. Canonical Sampling through Velocity Rescaling. J. Chem. Phys. 2007, 126, 014101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  139. Berendsen, H.J.C.; Postma, J.P.M.; van Gunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular Dynamics with Coupling to an External Bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef] [Green Version]
  140. Miyamoto, S.; Kollman, P.A. Settle: An Analytical Version of the SHAKE and RATTLE Algorithm for Rigid Water Models. J. Comput. Chem. 1992, 13, 952–962. [Google Scholar] [CrossRef]
  141. Hess, B.; Bekker, H.; Berendsen, H.J.C.; Fraaije, J.G.E.M. LINCS: A Linear Constraint Simulations. J. Comput. Chem. 1997, 18, 1463–1472. [Google Scholar] [CrossRef]
  142. Koziara, K.B.; Stroet, M.; Malde, A.K.; Mark, A.E. Testing and Validation of the Automated Topology Builder (ATB) Version 2.0: Prediction of Hydration Free Enthalpies. J. Comput. Aided. Mol. Des. 2014, 28, 221–233. [Google Scholar] [CrossRef]
  143. Malde, A.K.; Zuo, L.; Breeze, M.; Stroet, M.; Poger, D.; Nair, P.C.; Oostenbrink, C.; Mark, A.E. An Automated Force Field Topology Builder (ATB) and Repository: Version 1.0. J. Chem. Theory Comput. 2011, 7, 4026–4037. [Google Scholar] [CrossRef] [PubMed]
  144. Stroet, M.; Caron, B.; Visscher, K.M.; Geerke, D.P.; Malde, A.K.; Mark, A.E. Automated Topology Builder Version 3.0: Prediction of Solvation Free Enthalpies in Water and Hexane. J. Chem. Theory Comput. 2018, 14, 5834–5845. [Google Scholar] [CrossRef]
  145. Tayar, N.E.; Tsai, R.-S.; Testa, B.; Carrupt, P.-A.; Leo, A. Partitioning of Solutes in Different Solvent Systems: The Contribution of Hydrogen-Bonding Capacity and Polarity. J. Pharm. Sci. 1991, 80, 590–598. [Google Scholar] [CrossRef] [PubMed]
  146. Seiler, P. Interconversion of Lipophilicities from Hydrocarbon/Water Systems into the Octanol/Water System. Eur. J. Med. Chem. 1974, 9, 473–479. [Google Scholar]
Figure 1. Virtual screening steps applied in this study to identify new compounds from NatProdB (555 hits) with potential antiviral activity against the SARS-CoV-2 Main protease, including: ADME property filters; toxicity alert filter for chemical groups; and docking simulations and molecular dynamics filters.
Figure 1. Virtual screening steps applied in this study to identify new compounds from NatProdB (555 hits) with potential antiviral activity against the SARS-CoV-2 Main protease, including: ADME property filters; toxicity alert filter for chemical groups; and docking simulations and molecular dynamics filters.
Ijms 22 11739 g001
Figure 2. Three-dimensional structure of the SARS-CoV-2 main protease (SARS-CoV-2 Mpro, PDB ID: 5R82, blue cartoon). (A) Binding sites (surface and pink sticks) found by SeeSAR version 9.2 (BioSolveIT GmbH) in the active site (I) are colored in pink (25 residues), at site II (19 residues) in yellow, and at site III (15 residues) in green. This image was generated by the SeeSAR program. (B) Main pockets (color surfaces) predicted by PockDrug server [47]. The active site (I) is colored in pink (volume = 1892 Å3, druggability score = 0.63 ± 0.09), at site II (volume = 643 Å3, druggability score = 0.35 ± 0.08) in yellow, and at site III (volume = 760 Å3, druggability score = 0.33 ± 0.03) in green. The remaining gray sites presented a druggability score of 0.55 ± 0.10 (upper), 0.98 ± 0.01 (middle), and 0.61 ± 0.04 (lower). (C) Main hot spots (mesh and sticks) calculated by FTMap server [49]. The active site is colored in pink and was classified as druggable (≥16 probes) [54], the second site (yellow) is borderline (13 ≤ probes < 16), and pink and yellow sticks show the location of the probe. The B and C images were generated by educational Pymol 2.4.1 [55].
Figure 2. Three-dimensional structure of the SARS-CoV-2 main protease (SARS-CoV-2 Mpro, PDB ID: 5R82, blue cartoon). (A) Binding sites (surface and pink sticks) found by SeeSAR version 9.2 (BioSolveIT GmbH) in the active site (I) are colored in pink (25 residues), at site II (19 residues) in yellow, and at site III (15 residues) in green. This image was generated by the SeeSAR program. (B) Main pockets (color surfaces) predicted by PockDrug server [47]. The active site (I) is colored in pink (volume = 1892 Å3, druggability score = 0.63 ± 0.09), at site II (volume = 643 Å3, druggability score = 0.35 ± 0.08) in yellow, and at site III (volume = 760 Å3, druggability score = 0.33 ± 0.03) in green. The remaining gray sites presented a druggability score of 0.55 ± 0.10 (upper), 0.98 ± 0.01 (middle), and 0.61 ± 0.04 (lower). (C) Main hot spots (mesh and sticks) calculated by FTMap server [49]. The active site is colored in pink and was classified as druggable (≥16 probes) [54], the second site (yellow) is borderline (13 ≤ probes < 16), and pink and yellow sticks show the location of the probe. The B and C images were generated by educational Pymol 2.4.1 [55].
Ijms 22 11739 g002
Figure 3. SARS-CoV-2 Mpro structures are shown in B-factor putty mode in PyMOL. The highest B-factor in each structure is colored in red and the lowest in dark blue, as indicated by the B-factor scale bar. The thickness of the protein backbone also is proportional to the B-factors. The catalytic residues (His41 and Cys145) are displayed as sticks. This image was generated by educational Pymol 2.4.1 [64].
Figure 3. SARS-CoV-2 Mpro structures are shown in B-factor putty mode in PyMOL. The highest B-factor in each structure is colored in red and the lowest in dark blue, as indicated by the B-factor scale bar. The thickness of the protein backbone also is proportional to the B-factors. The catalytic residues (His41 and Cys145) are displayed as sticks. This image was generated by educational Pymol 2.4.1 [64].
Ijms 22 11739 g003
Figure 4. Hydrogen bond stability in SARS-CoV-2 Mpro complexes as calculated by HbMap2Grace after 100 ns of molecular dynamics simulation.
Figure 4. Hydrogen bond stability in SARS-CoV-2 Mpro complexes as calculated by HbMap2Grace after 100 ns of molecular dynamics simulation.
Ijms 22 11739 g004
Figure 5. Surface molecular area (Å2) of SARS-CoV-2 Mpro complexes as calculated by SurfinMD after 100 ns of molecular dynamics simulation.
Figure 5. Surface molecular area (Å2) of SARS-CoV-2 Mpro complexes as calculated by SurfinMD after 100 ns of molecular dynamics simulation.
Ijms 22 11739 g005
Figure 6. SARS-CoV-2 Mpro main residues that interacted with the best-docked NatProDB compounds throughout the productive phase of the molecular dynamics simulation. Mpro is shown as a cartoon (blue) with its main interacting residues labeled (blue sticks); b01 is shown in green sticks; b02 in cyan; RZS in magenta; and E-64 in orange. Polar interactions are depicted as yellow dashed lines. This image was generated by educational Pymol 2.4.1 [64].
Figure 6. SARS-CoV-2 Mpro main residues that interacted with the best-docked NatProDB compounds throughout the productive phase of the molecular dynamics simulation. Mpro is shown as a cartoon (blue) with its main interacting residues labeled (blue sticks); b01 is shown in green sticks; b02 in cyan; RZS in magenta; and E-64 in orange. Polar interactions are depicted as yellow dashed lines. This image was generated by educational Pymol 2.4.1 [64].
Ijms 22 11739 g006
Figure 7. Residue contributions to the binding energy of SARS-CoV-2 Mpro complexes. The main residues with energy interactions (ΔEbinding > ±5 kJ/mol) are highlighted. For Mpro sub-site definitions, see the text.
Figure 7. Residue contributions to the binding energy of SARS-CoV-2 Mpro complexes. The main residues with energy interactions (ΔEbinding > ±5 kJ/mol) are highlighted. For Mpro sub-site definitions, see the text.
Ijms 22 11739 g007
Table 1. Best docking results carried out by SeeSAR version 9.2 (BioSolveIT GmbH) on the active site of the SARS-CoV-2 main protease (SARS-CoV-2 Mpro) using the Natural Products Database of the Bahia Semi-Arid region (NatProDB), E-64 (positive control), and RZS (crystallographic ligand, PDB ID: 5R82).
Table 1. Best docking results carried out by SeeSAR version 9.2 (BioSolveIT GmbH) on the active site of the SARS-CoV-2 main protease (SARS-CoV-2 Mpro) using the Natural Products Database of the Bahia Semi-Arid region (NatProDB), E-64 (positive control), and RZS (crystallographic ligand, PDB ID: 5R82).
#CompoundsChemical StructureHydrogen Bond Donors (HBD) 1Hydrogen Bond Acceptors (HBA)Lipinki’s Rule 2 ViolationMolecular Surface Area (Å2) 3Solvent Accessible Surface (Å3) 4Estimated Affinity Range 5 (µM)Inter Clash Type 6
b01VE0DIA0AF Ijms 22 11739 i001121134.92478.23.29–326.4+
b02VE0PPA0AF Ijms 22 11739 i002550150.25481.34.42–439.60
b03VE0ISA0AF Ijms 22 11739 i003150152.87645.37.24–719.2+
b04VE0FKA0AF Ijms 22 11739 i004250150.34576.217.25–1714.3+
b05VE0FEA0SF Ijms 22 11739 i005460150.17519.518.18–1806.50
b06VE0JDA0SF Ijms 22 11739 i006350121.72355.920.82–2068.70
b07VE0ZDA0AF Ijms 22 11739 i007150151.74652.241.3–4103.0+
b08VE0NCA0SF Ijms 22 11739 i00812072.59309.251.97–5164.2+
b09VE0KJA0SI Ijms 22 11739 i009380148.67599.852.40–5206.50
b10VE0NHA0SF Ijms 22 11739 i010360122.04543.653.85–5350.2+
255Control(RZS)- Ijms 22 11739 i01113065.69350.4504.29–50,103.8+
562Positive control(E-64)- Ijms 22 11739 i0127101145.91638.396,164.969, 455,201.90
1 The number of hydrogen bond donors (HBD), acceptors (HBA), and Lipinski’s rule violation were determined from the FAFF-Drugs server (https://fafdrugs4.rpbs.univ-paris-diderot.fr/, accessed on 21 April 2021) [56,57,58,59]. 2 Lipinski’s rule was described in [60]. 3 The molecular surface area were obtained from the pkCSM server (http://biosig.unimelb.edu.au/pkcsm/, accessed on 15 April 2021) [27]. 4 The solvent accessible surface (Å3) was calculated from the Marvin Sketch 18.8.0 Geometry plugin (http://www.chemaxon.com) employing the solvent radius 1.4 Å and pH = 7.4. 5 The estimated affinity (EA) ranges mark the lower and upper bound of the estimated affinities (±2 log-units). 6 Inter clash type: atomic clashes between the protein and docked ligand atoms as predicted by SeeSAR9.2. Inter clash type are classified as + (no clashes), 0 (few clashes), and − (many clashes).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Silva, R.C.; Freitas, H.F.; Campos, J.M.; Kimani, N.M.; Silva, C.H.T.P.; Borges, R.S.; Pita, S.S.R.; Santos, C.B.R. Natural Products-Based Drug Design against SARS-CoV-2 Mpro 3CLpro. Int. J. Mol. Sci. 2021, 22, 11739. https://doi.org/10.3390/ijms222111739

AMA Style

Silva RC, Freitas HF, Campos JM, Kimani NM, Silva CHTP, Borges RS, Pita SSR, Santos CBR. Natural Products-Based Drug Design against SARS-CoV-2 Mpro 3CLpro. International Journal of Molecular Sciences. 2021; 22(21):11739. https://doi.org/10.3390/ijms222111739

Chicago/Turabian Style

Silva, Rai C., Humberto F. Freitas, Joaquín M. Campos, Njogu M. Kimani, Carlos H. T. P. Silva, Rosivaldo S. Borges, Samuel S. R. Pita, and Cleydson B. R. Santos. 2021. "Natural Products-Based Drug Design against SARS-CoV-2 Mpro 3CLpro" International Journal of Molecular Sciences 22, no. 21: 11739. https://doi.org/10.3390/ijms222111739

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop