Next Article in Journal
Market Segmentation by Motivations of Urban Forest Users and Differences in Perceived Effects
Next Article in Special Issue
Structural and Dynamic Insights into the W68L, L85P, and T87A Mutations of Mycobacterium tuberculosis Inducing Resistance to Pyrazinamide
Previous Article in Journal
Spatial-Temporal Pattern and Evolution Trend of the Cultivated Land Use Eco-Efficiency in the National Pilot Zone for Ecological Conservation in China
Previous Article in Special Issue
Determination of Novel Anti-Cancer Agents by Targeting OGG1 Enzyme Using Integrated Bioinformatics Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rescuing the Host Immune System by Targeting the Immune Evasion Complex ORF8-IRF3 in SARS-CoV-2 Infection with Natural Products Using Molecular Modeling Approaches

Department of Medical Biotechnology, College of Applied Medical Sciences, Qassim University, Buraydah 52571, Saudi Arabia
Int. J. Environ. Res. Public Health 2022, 19(1), 112; https://doi.org/10.3390/ijerph19010112
Submission received: 5 November 2021 / Revised: 20 December 2021 / Accepted: 20 December 2021 / Published: 23 December 2021

Abstract

:
The perennial emergence of SARS-CoV-2 and its new variants causing upper respiratory complexities since December 2019 has aggravated the pandemic situation around the world. SARS-CoV-2 encodes several proteins among which ORF8 is a novel factor that is unique to SARS-CoV-2 only and is reported to help the virus in disease severity and immune evasion. ORF8-IRF3 complex induces endoplasmic reticulum stress, thus helps in the evasion of immune response. Consequently, targeting the ORF8-IRF3 complex is considered as a prime target for the discovery of novel drugs against SARS-CoV-2. In this regard, computational methods are of great interest to fast track the identification and development of novel drugs. Virtual screening of South African Natural Compounds Database (SANCDB), followed by docking and molecular dynamics (MD) simulation analysis, were performed to determine novel natural compounds. Computational molecular search and rescoring of the SANCDB database followed by induced-fit docking (IFD) protocol identified Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside (SANC00850), Tribuloside (SANC01050), and Rutin (SANC00867) are the best scoring compounds. Structural-dynamic properties assessment revealed that these three compounds have stable dynamics, compactness, and a higher number of hydrogen bonds. For validation, we used MM/GBSA, in silico bioactivity estimation and dissociation constant (KD) approaches, which revealed that these compounds are the more potent inhibitors of the ORF8-IRF3 complex and would rescue the host immune system potentially. These compounds need further in vitro and in vivo validations to be used as therapeutics against SARS-CoV-2 to rescue the host immune system during COVID-19 infection.

1. Introduction

The perennial emergence of SARS-CoV-2 and its new variants causing upper respiratory complexities since December 2019 has aggravated the pandemic situation. This pandemic has caused 238,378,962 mortalities and 4,863,187 deaths worldwide as of 10 October 2021. This virus causing coronavirus Disease-19 (COVID-19) with a single strand RNA is continuously subjected to evolution by bringing new mutations [1,2]. The genomic organization of this virus shares commonalities with other beta coronaviruses. The six functional open reading frames (ORFs) are well-arranged in order from 5′ to 3′ replicas (ORF1a/ORF1b), membrane (M), envelope (E), spike (S), and nucleocapsid (N). In contrast, the nonstructural proteins include 3-chymotrypsin such as protease, papain-like protease, and RNA-dependent RNA polymerase, which are encoded by the ORF regions [3,4,5]. Among these, spike protein is essential which mediate cellular attachment and entry to the cell and starts the infections. ORF8 is a novel factor that is unique to SARS-CoV-2 only and is reported to help the virus in disease severity and immune evasion [6,7].
ORF8 with 121 amino acids display captivating features and is more susceptible to mutations. Among the key events regulated by ORF8 involve downregulation of type I interferon signalling and histocompatibility complex class I (MHC-I) pathways [8,9]. Many studies reported the interaction of ORF8 with IRF3 to hijack the host immune system and potentially evade [10]. Interaction of ORF8-IRF3 is reported to target different aspects of the IFN signaling cascade, thus suppressing the host immune system and proceed a successful infection. Many studies also reported mutations in ORF8 such as W45L, V62L, and L84S structurally characterized are reported to more robustly evade the host immune system [6]. ORF8 is also reported to form intracellular aggregates in the Epithelial Cells of the lungs. It is also reported to induce endoplasmic reticulum stress, thus helps in the evasion of immune response [7]. Thus, targeting the ORF8-IRF3 is considered as a prime target for the discovery of novel drugs against SARS-CoV-2 [11].
Due to the multi-faceted role of ORF8-IRF3 complex, it is important to abolish this complex either with small natural molecules, peptides inhibitors, or nanobodies. In this regard, computational methods are of great interest to fast track the identification and development of novel drugs [12,13,14,15,16,17,18]. Using computational methods various drugs for other targets, vaccines, peptides, and nanobodies are discovered against SARS-CoV-2 [19,20,21,22,23]. Thus, employing the same approach here would greatly assist the drug development cycle against this immune evasion complex.
Herein we used computational modeling and simulation approaches to target the interface of the ORF8-IRF3 complex with high-affinity natural products [24]. South African Natural Compounds Database (SANCDB) (http://sancdb.rubi.ru.ac.za/ (accessed on 10 October 2021) [25] containing compounds from plant sources and marine life was screened against the interface residues of IRF3. Computational molecular search and rescoring of the whole database followed by induced-fit docking (IFD) protocol were employed to identify the best scoring compounds. All-atoms simulation and binding free energy approaches were then applied to validate the potential of these compounds. Finally in silico bioactivity prediction and dissociation constant (KD) approaches validated the efficacy of these final hits. In sum the shortlisted compounds here provide basis for the development of next-generation drugs focusing on the host factors to rescue the host immune system.

2. Materials and Methods

2.1. IRF3, Database Retrieval and Preparation

The crystal structure of IRF3 available online was retrieved from RCSB using accession ID 1QWT. The retrieved structure was prepared, refined, and minimized using Chimera and GalaxyRefine (USCF, CA, USA) webserver [26,27,28].

2.2. Drugability Assessment of the Binding Interface

To find the druggable pockets in the binding interface of IRF3, the DoGSiteScorer web server is freely available at http://dogsite.zbh.uni-hamburg.de (accessed on 10 October 2021) was used [29]. The servers predicts the druggability potential of the interacting interface of the two proteins.

2.3. Molecular Search against ORF8-IRF3 Interface

Identification of potential natural compounds through a molecular search for IRF3, using SANCDB database. All the drugs, 1200 total including flavonoids, alkaloids, phenylpropanoids, polyketides, and terpenoids in the database were prepared and minimized using PyRx tool [30]. The interface residues 189–240 of IRF3 required for interaction with ORF8 were selected as reported by a previous literature [6]. For instance the role of medicinal natural compounds such as flavonoids, alkaloids, phenylpropanoids, polyketides, and terpenoids is very obvious against SARS-CoV-2 and has been previously reported by various studies [12,13,20,31]. Screening against the binding site residues was performed, and the best scoring hits were selected for re-docking and rescoring.

2.4. Rescoring and Ranking Using Induced-Fit Docking (IFD) of the Top Hits

The best scoring compounds were re-evaluated using IFD protocol using 64 exhaustiveness to remove and false-positive results and confirm the final hits. AutoDockFR–AutoDock for Flexible Receptors (ADFR) [32] based on AutoDock4 scoring function was used to re-score the best hits using IFD algorithm. This approach minimizes the internal energy of the receptor and achieves the receptor sidechain best conformation by employing 14 different sidechains to boost the frequency of accurate docking. In comparison with AutoDock Vina, AutoDockFR attained higher accuracy in a cross-validation check, and also, the speed of the docking is much higher.

2.5. Validation of the Top Hits

2.5.1. Molecular Dynamics Simulation of the Top Hits

Structural-dynamic features of the top hits were explored through all-atoms simulation using the AMBER18 package [33]. The Amber general force field (GAFF) and ff14SB forcefields were employed for drugs and IRF3 topologies. Solvation with TIP3P water box and CL- ions were added for neutralization of each system. Energy minimization of each complex was achieved in two steps, followed by heating and equilibration. For the annotation of long-range electrostatic interactions, a Particle Mesh Ewald (PME) algorithm [34] and 1.4 nm cut-off value for vdW and Columbic interactions of short-range was set. For temperature regulation, Langevin thermostat (300 K) and Berendsen barostate were used for pressure control. A time step of 2fs and a total simulation time of 200 ns for each complex were performed. The dynamic stability, flexibility, compactness to show the binding and unbinding events of the drugs in the cavity, and other features of the ligands-protein complexes were calculated by utilizing CPPTRAJ and PTRAJ modules [35].

2.5.2. The Binding Free Energy Calculations

The top hits in complex with IRF3 were subjected to binding free energy calculation using the script MMPBSA.PY by the whole trajectory [36,37,38,39]. Various studies have employed this method to compute the TBE of different compounds [40,41,42,43,44].
ΔGbind = ΔGcomplex − [ΔGreceptor + ΔGligand]
Here, ΔGbind stands for the total free binding energy, while others denote the free energy of the protein, the ligand and complex. The equation given below was then used to estimate each specific energy term that contributes to the total free energy:
G = Gbond + Gele + GvdW + Gpol + Gnpol
Bonded, electrostatic, polar, non-polar, and van der Waal energy terms are represented by the above equation.

2.5.3. In Silico Bioactivity Prediction

To validate the binding of the top hits, we predicted the bioactivity of each compound against different classes of proteins using the Molinspiration web tool. The scores between −0.5 to 0.5 demonstrate excellent inhibitory activity.

2.5.4. Dissociation Constant (KD) Evaluation

Binding kinetics is an important parameter to estimate the strength of the association of drug and a protein and is widely used by a larger number of studies. Herein, to estimate the dissociation constant of the top hits, we used PROtein binDIng enerGY prediction (PRODIGY) server (https://wenmr.science.uu.nl/prodigy/ (accessed on 18 October 2021)) [45].

3. Results and Discussion

3.1. Structural Retrieval and Analysis of Druggable Site

The pandemic caused by SARS-CoV-2 has troubled the world by spreading exponentially to every nook and cranny of the world. The pathophysiology of this virus demonstrates that it principally affects the upper respiratory canal by entering the host cell through spike protein. The adherence of the spike protein, receptor-binding domain to the host ACE2 protein, allows the entry to the host cell, and it is considered the primary drug target for novel drug discovery [46,47]. However, the role of other proteins, particularly the ORF8-IRF3 complex, is very important in the disease’s severity and immune evasion. Therefore, abrogation of this complex would rescue the host immune system and thus will ensure the efficient neutralization of the virus. Hence, drugs that can target this complex are required to rescue the host immune system. The current study uses state-of-the-art methods to search potential drugs from medicinal plants sources using genetic algorithm, IFD-algorithm, all-atoms simulation, binding free energy calculation, bioactivity prediction, and dissociation constant determination. Herein, we targeted the IRF3 interface residues (189–240) using SANCDB database, which contains drugs from natural plants and marine resources. The complex of ORF8-IRF3 was obtained through a protein-protein docking approach to confirm the binding interface with a previous study [6]. Our analysis revealed the same binding interface for both the structures shown in Figure 1A. The interface residues of IRF3 were then subjected to druggability assessment which revealed that the binding interface is a druggable site and can be targeted using drugs. The druggable residues are shown in Figure 1B,C.

3.2. Ranking the Best Compounds and Rescoring

Keeping in view the importance of ORF8-IRF3 complex in immune evasion, we screened phytomedicines from natural plants and marine sources. Using a multi-steps screening approach, we ranked the best compounds which could abrogate the IRF3 interaction with ORF8. In the first round, the whole database was screened using PyRx. The results revealed the docking scores range from −8.23 to −4.11 kcal/mol. Then we screened the best compounds with docking score >−6.0 kcal/mol. With this criterion, only 202 compounds were shortlisted for rescoring and analysis. In the second step, AutoDock Vina-based docking of these 202 compounds was performed, and the docking score range from −8.52 to −4.83 kcal/mol. Finally, using ADFR, IFD docking of the top 10% compounds were redocked and scored. The docking score for these 41 compounds using IFD approach range from −10.52 to −4.97 kcal/mol. Among these compounds with good docking scores were selected, which shortlisted only six compounds to have the best scores and interactions. These compounds include SANC00524 (−10.52 kcal/mol), SANC00526 (−9.53 kcal/mol), SANC00850 (−6.93 kcal/mol), SANC01085 (−6.37 kcal/mol), SANC00867 (−6.27 kcal/mol) and SANC00395 (−5.98 kcal/mol). The top six hits were then checked for Lipinski’s rule, which revealed that only three compounds SANC00850, SANC01085, and SANC00867 passed the criteria and were subjected to further analysis while the rest were discarded. The docking scores, database ID, 2D structure, and names of these compounds with the highest affinity towards the protein are given in Table 1. The binding mode of each compound is given in details below. These three compounds can be seen that they all possess almost similar scaffold with benzene rings surrounded by hydroxyl groups which make it possible for many interactions with the residues.

3.2.1. Binding Mode of Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside

Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside (SANC00850) isolated from Myrothamnus flabellifolia Welw is a flavonoid used against respiratory diseases and urinary tract infections in African countries [48]. It has also been shown that this compound induces antiviral effects against herpes simplex virus type 1. Here, with the docking score −6.93 kcal/mol, it has demonstrated the best activity against the SARS-CoV-2 ORF8-IRF3 complex. As shown in Figure 2, the compound interacts with the key interface residues and establishes several hydrogen bonds. Among the hydrogen bonding residues Glu201, Thr219, Ile220, Ser221, Glu293, ile395, Asp392, and Ser398 are involved. In this complex the [O2] of Glu201 interacts with the [O-] acting as acceptor here of the ligand. The [O2] of Thr219 amino acids contacted with [O2] at position 508 of the ligand atom. Furthermore the interactions formed by these atoms of amino acids and SANC00580 were reported to be responsible for the strongest binding 3722 [O3]...521 [O2], 551 [O2]...3724 [O2], 3724 [O2]...554 [O2], 3712 [O2]...3157 [O.co2], 3200 [O2]...3714 [O2], and 3714 [O2]...3247 [O2], respectively. For instance, these residues are previously reported to stabilize the IRF3-ORF8 complex, thus showing the inhibitory potential of this compound [6]. The binding mode of SANC00850 with IRF3 is shown in Figure 2.

3.2.2. Binding Mode of Rutin

Rutin (SANC01085) or Rutoside or Quercetin 3-O-rutinoside is a flavonoid isolated from various plants that exhibit hepatoprotective, anti-inflammatory, and anti-microbial properties. It was originally isolated from Aspalathus linearis, a south African tea widely used against cytochrome P450 [49]. This compound possesses a similar scaffold as Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside. The docking score for Rutin was −6.37 kcal/mol with seven hydrogen bonds established by 196Leu 3726 [O3]...137 [O2], 196Leu 3724 [O3]...137 [O2], 199Gly 183 [Nam]...3698 [O3], 201Glu 3712 [O2]...214 [O-], 223Pro 3720 [O3]...576 [O2], 224Glu 3718 [O3]...594 [O.co2], 224Glu 3716 [O3]...594 [O.co2], 276Arg 3722 [O3]...1372 [Ng+], 276Arg 1371 [Ng+]...3724 [O3], 276Arg 1372 [Ng+]...3724 [O3], 278Gly 1396 [Nam]...3722 [O], respectively. Among these, most of the interacting residues are conserved between Rutin and Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside, thus justify the inhibition potential of these compounds. The binding mode of SANC01085 with IRF3 is shown in Figure 3.

3.2.3. Binding Mode of Tribuloside

Tribuloside is also a flavonoid originally Tribulus terrestris L and characterized as an anti-mycobacterial compound. It is also known as kaempferol-3-β-d-(6″-p-coumaroyl) glucoside which is a derivate of Kaempferol [50]. Kaempferol has been reported to potentially inhibit the SARS-CoV-2 main protease (3CLpro), thus hinder the proteolysis of non-structural proteins [13,40]. Tribuloside with the docking score −6.27 kcal/mol formed nine hydrogen bonds with the key residues including Glu201, Glu203, Gln218, Thr219, Ile220, Ser221, Asp392, and Asn397. Deep interaction analysis revealed 201GLU 3708 [O2]...214 [O-], 203GLU 3716 [O3]...252 [O.co2], 203GLU 3718 [O3]...252 [O.co2], 218GLN 508 [Nam]...3712 [O2], 219THR 3712 [O2]...521 [O2], 221SER 551 [Nam]...3714 [O3], 221SER 3714 [O3]...556 [O3], 392ASP 3720 [O2]...3158 [O-] and 397ASN 3230 [Nam]...3696 [O2] formed interactions, respectively. This shows the inhibitory potential of these compounds against the IRF3-ORF8 complex. The binding mode of SANC01085 with IRF3 is shown in Figure 4.

3.3. Structural Dynamic Features of the Best Hits Complexes

3.3.1. Structural Stability Assessment

Dynamic stability of a ligand-protein complex is very important to demonstrate the inhibitory potential of a drug. Thus herein, we used root mean square deviation (RMSD) as a function of time to estimate the stability of each protein-ligand complex over 100 ns simulation time period. As given in Figure 5, it can be seen that all the complexes remained stable during the simulation with little deviations at different time intervals were experienced by the complexes. In the case of Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside-IRF3 complex, the system reached stability at 1.5 Å and at 5 ns. The structure did not face any convergence during the first 35 ns, however, a little deviation between 35–38 ns was showed then the structure remained stable. At different time intervals, particularly between 50–70 ns, faced a little perturbation, but overall, the complex remained stable. This shows the stable binding of Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside to the IRF3 interface. On the other side, the Tribuloside-IRF3 complex remained comparatively more stable than the Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside-IRF3 complex. Although the RMSD increased gradually but the average RMSD remained lower than 2.0 Å. A small deviation between 30–35 ns was experienced by Tribuloside-IRF3 complex but the complex dynamically remained stable. Moreover, the Rutin-IRF3 complex reported structural perturbation at different time intervals, particularly between 80–88 ns. Despite the structural deviation during the simulation, the average RMSD for this complex also remained lower (~2.2 Å). Consequently, this shows that all the compounds are stably bound to the target interface residues and exhibit the potential inhibitory properties. The RMSDs of all the complexes are shown in Figure 5.

3.3.2. Structural Compactness Calculation

Understanding the structural compactness in a dynamic environment using simulation helps in understanding the packing of protein-ligand complexes and demonstrates the binding and unbinding events that happened during the simulation. Herein, the radius of gyration (Rg) was estimated for each complex as a function of time. As given in Figure 6, the structural compactness results are in strong uniformity with the RMSD results. In the case of IRF3-Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside complex, the average Rg was reported to be 19.2 Å. Initially, the Rg increased until 10 ns but then decreased back and remained consistent until 30 ns. A significant drift was experienced between 40–50 ns, but then the structure remained compact. On the other hand, IRF3-Tribuloside demonstrated uniform results such as RMSD by remaining more compact throughout the simulation. The average Rg value was reported to be 19.0 Å for IRF3-Tribuloside complex. Furthermore, IRF3-Rutin complex displayed a higher Rg during the first 30 ns but then gradually decreased, and the average Rg was reported to be 19.1 Å. This shows that despite the binding and unbinding events, the top hits molecules stably bound the target protein steered by different bonds and biochemical events. The Rg(s) of all the complexes are shown in Figure 6.

3.3.3. Residual Flexibility Estimation

We computed the residual flexibility of all the complexes to provide an understanding of the strength contributed by each interacting residue in intermolecular association, the molecular recognition, and the probable effect on the global function of the biological molecules. Herein we calculated the residual flexibility as a root mean square fluctuation (RMSF). The results revealed that all the complexes exhibit minimal fluctuation stabilized by the binding of these molecules except at different regions where uniformly higher fluctuation was observed. These regions were witnessed as loops, and particularly, the terminal loop demonstrated higher fluctuation. Among the higher fluctuating regions, 40–50, 162–173, 190–200, and 235–240 demonstrated higher fluctuation. The RMSFs of all the complexes are shown in Figure 7.

3.3.4. Hydrogen Bonding Analysis

Protein-drugs association is primarily directed by a number of factors, among which hydrogen bonds and hydrophobic interactions are vital components. Hydrogen bonds in the protein interfaces are always engaged by water molecules which strive with the hydrogen bonding between the molecules [51]. The processes behind protein-protein/drugs coupling, as well as the implications to which hydrogen bonds play a role in this association, are very important [52]. Therefore, to comprehend the bonding pattern between the IRF3, Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside, Rutin and Tribuloside, post-simulation hydrogen bonding analysis of each complex was performed to inform the binding specificity and biochemical events steered by hydrogen bonding. Herein, we calculated the total number of hydrogen bonds in each complex as a function of time. In the case of Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside-IRF3 complex the average number of hydrogen bonds was reported to be 115, the average number of hydrogen bonds in Tribuloside-IRF3 complex were 113, while in Rutin-IRF3 complex, the average hydrogen bonds were also reported to be 115. This shows that hydrogen bonds formation and breakage events took place during the simulation. The total number of hydrogen bonds in all the complexes is shown in Figure 8.

3.4. Binding Free Energy Calculation

Estimation of binding free energy at the residue level of protein-ligand(s) complexes is the most frequently used approach in discovering small molecule inhibitors using in silico methods. This approach cross-validates the docking scores and re-evaluates the small molecule’s inhibitory potential interacting with a targeted protein. They are more comparatively more reliable than conventional docking scoring functions and computationally inexpensive compared to other methods such as Alchemical free energy approaches [53]. To estimate the binding free energies of all the complexes including vdW, electrostatic, SA, Gb and the total binding energy each trajectory was subjected to free energy calculation using MM/GBSA approach. The binding free energy calculation results revealed that for each complex the vdW was reported to be −34.86 kcal/mol (SANC00850-IRF3), −25.20 kcal/mol (SANC01085-IRF3) and −36.84 kcal/mol (SANC00867-IRF3). The electrostatic energy for each complex was reported to be −27.89 kcal/mol, −29.92 kcal/mol and −35.38 kcal/mol, respectively. The total binding energy given conclusive insights revealed that the total binding energy for SANC00850-IRF3 complex was −41.41 kcal/mol. For SANC01085-IRF3 the total binding energy was −40.33 kcal/mol while −40.62 kcal/mol was reported for SANC00867-IRF3. The fixation of the tails and backbone of the SANC00850 in the binding pocket gives it the advantage to bind robustly and inhibits the IRF3. It can be seen that 11 hydrogen bonds were more in number in Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside-IRF3 complex because of the binding of whole scaffold in the active site. Moreover, the whole ligand is bound by the hydrogen bonds such as the two tails and the backbone. On the hand, the two ligands in the active sites are free either at one end or at the central backbone. Thus demonstrating the differences in the binding. The hydrogen bonding analysis also reports that hydrogen bonding reprogramming took place during the simulation and the number of average bonds validate the docking predictions. This data shows that these compounds interact with the interface residues of IRF3 more robustly, thus hinder the binding of partner protein ORF8 and rescue the host immune system. The MM/GBSA results for all the complexes are given in Table 2.

3.5. Bioactivity Prediction

Molinspiration predicted the bioactivity of the top hits. For Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside the predicted bioactivity score was −0.05, for Tribuloside the predicted bioactivity score was 0.05 and for Rutin, the predicted bioactivity was 0.12. This shows the strongest inhibitory potential of these compounds because the server predicts scores for different compounds where the score ranges from −0.5 to 0.5 demonstrate excellent inhibitory activity. This approach has been previously used to estimate the bioactivity of potential compounds against SARS-CoV-2 [20].

3.6. Dissociation Constant Estimation (KD)

KD is an important parameter to compute the binding strength of a drug. Thus, employing the same approach, the dissociation constant of the top three hits was estimated. The results revealed that for the Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside-IRF3 complex the ∆G was −8.3 kcal/mol, for Rutin-IRF3 complex, the ∆G was predicted to be −7.7 kcal/mol, while for Tribuloside-IRF3 complex, the ∆G was predicted to be −7.8 kcal/mol. It can be seen that 11 hydrogen bonds were more in number in Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside-IRF3 complex because of the binding of whole scaffold in the active site. Moreover, the whole ligand is bound by the hydrogen bonds such as the two tails and the backbone. On the hand, the two ligands in the active sites are free either at one end or at the central backbone. Thus demonstrating the differences in the binding. This shows that these compounds potentially interact with the target protein residues and curtail their binding to the interface of the interacting partner.

4. Conclusions

In conclusion, the current study uses computational screening and simulation approaches to identify novel drugs against the ORF8-ORF3 complex. This complex has been reported to suppress the host immune system and increase the disease severity. Through molecular search, we reported three drugs Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside, Tribuloside and Rutin, potentially interact with the interface residues and abrogate the complex. Structural-dynamics analysis validated the inhibitory potential of these compounds. The total binding energy given conclusive insights revealed that the total binding energy for SANC00850-IRF3 complex was −41.41 kcal/mol. For SANC01085-IRF3 the total binding energy was −40.33 kcal/mol while −40.62 kcal/mol was reported for SANC00867-IRF3. This data shows that these compounds interact with the interface residues of IRF3 more robustly, thus hinder the binding of partner protein ORF8 and rescue the host immune system. Thus, these compounds need further in vitro and in vivo validations to be used as therapeutics against SARS-CoV-2.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available within the article.

Acknowledgments

The researcher would like to thank the Deanship of Scientific Research, Qassim University for funding the publication of this project.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Benvenuto, D.; Giovanetti, M.; Ciccozzi, A.; Spoto, S.; Angeletti, S.; Ciccozzi, M. The 2019-new coronavirus epidemic: Evidence for virus evolution. J. Med. Virol. 2020, 92, 455–459. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ul Qamar, M.T.; 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] [PubMed]
  3. Maison, D.P.; Ching, L.L.; Shikuma, C.M.; Nerurkar, V.R. Genetic Characteristics and Phylogeny of 969-bp S Gene Sequence of SARS-CoV-2 from Hawai ‘i Reveals the Worldwide Emerging P681H Mutation. Hawai’i J. Health Soc. Welf. 2021, 80, 52. [Google Scholar]
  4. Voloch, C.M.; Silva, R.d.F.; de Almeida, L.G.P.; Cardoso, C.C.; Brustolini, O.J.; Gerber, A.L.; Guimarães, A.P.d.C.; Mariani, D.; Costa, R.M.d.; Ferreira, O.C.; et al. Genomic characterization of a novel SARS-CoV-2 lineage from Rio de Janeiro, Brazil. J. Virol. 2021, 95, e00119-21. [Google Scholar] [CrossRef]
  5. Tang, D.; Comish, P.; Kang, R. The hallmarks of COVID-19 disease. PLoS Pathog. 2020, 16, e1008536. [Google Scholar] [CrossRef] [PubMed]
  6. Rashid, F.; Suleman, M.; Shah, A.; Dzakah, E.E.; Wang, H.; Chen, S.; Tang, S. Mutations in SARS-CoV-2 ORF8 Altered the Bonding Network With Interferon Regulatory Factor 3 to Evade Host Immune System. Front. Microbiol. 2021, 12. [Google Scholar] [CrossRef]
  7. Rashid, F.; Dzakah, E.E.; Wang, H.; Tang, S. The ORF8 protein of SARS-CoV-2 induced endoplasmic reticulum stress and mediated immune evasion by antagonizing production of interferon beta. Virus Res. 2021, 296, 198350. [Google Scholar] [CrossRef]
  8. Park, M.D. Immune evasion via SARS-CoV-2 ORF8 protein? Nat. Rev. Immunol. 2020, 20, 408. [Google Scholar] [CrossRef]
  9. Zhang, Y.; Zhang, J.; Chen, Y.; Luo, B.; Yuan, Y.; Huang, F.; Yang, T.; Yu, F.; Liu, J.; Liu, B. The ORF8 protein of SARS-CoV-2 mediates immune evasion through potently downregulating MHC-I. Proc. Natl. Acad. Sci. USA 2021, 118, e2024202118. [Google Scholar] [CrossRef]
  10. Hassan, S.S.; Aljabali, A.A.; Panda, P.K.; Ghosh, S.; Attrish, D.; Choudhury, P.P.; Seyran, M.; Pizzol, D.; Adadi, P.; Abd El-Aziz, T.M. A unique view of SARS-CoV-2 through the lens of ORF8 protein. Comput. Biol. Med. 2021, 133, 104380. [Google Scholar] [CrossRef]
  11. Zinzula, L. Lost in deletion: The enigmatic ORF8 protein of SARS-CoV-2. Biochem. Biophys. Res. Commun. 2021, 538, 116–124. [Google Scholar] [CrossRef]
  12. Humayun, F.; Khan, A.; Ahmad, S.; Yuchen, W.; Wei, G.; Nizam-Uddin, N.; Hussain, Z.; Khan, W.; Zaman, N.; Rizwan, M.; et al. Abrogation of SARS-CoV-2 interaction with host (NRP1) Neuropilin-1 receptor through high-affinity marine natural compounds to curtail the infectivity: A structural-dynamics data. Comput. Biol. Med. 2021, 141, 104714. [Google Scholar] [CrossRef] [PubMed]
  13. Khan, A.; Ali, S.S.; Khan, M.T.; Saleem, S.; Ali, A.; Suleman, M.; Babar, Z.; Shafiq, A.; Khan, M.; Wei, D.-Q. Combined drug repurposing and virtual screening strategies with molecular dynamics simulation identified potent inhibitors for SARS-CoV-2 main protease (3CLpro). J. Biomol. Struct. Dyn. 2020, 39, 4659–4670. [Google Scholar] [CrossRef] [PubMed]
  14. Khan, A.; Zia, T.; Suleman, M.; Khan, T.; Ali, S.S.; Abbasi, A.A.; Mohammad, A.; Wei, D.-Q. Higher infectivity of the SARS-CoV-2 new variants is associated with K417N/T, E484K, and N501Y mutants: An insight from structural data. J. Cell. Physiol. 2021, 236, 7045–7057. [Google Scholar] [CrossRef] [PubMed]
  15. Khan, A.; Khan, S.; Saleem, S.; Nizam-Uddin, N.; Mohammad, A.; Khan, T.; Ahmad, S.; Arshad, M.; Ali, S.S.; Suleman, M. Immunogenomics Guided Design of Immunomodulatory Multi-Epitope Subunit Vaccine against the SARS-CoV-2 new Variants, and its Validation through in Silico Cloning and Immune Simulation. Comput. Biol. Med. 2021, 133, 104420. [Google Scholar] [CrossRef] [PubMed]
  16. Altharawi, A.; Ahmad, S.; Alamri, M.A.; ul Qamar, M.T. Structural insight into the binding pattern and interaction mechanism of chemotherapeutic agents with Sorcin by docking and molecular dynamic simulation. Colloids Surf. B Biointerfaces 2021, 208, 112098. [Google Scholar] [CrossRef]
  17. Muneer, I.; Ahmad, S.; Naz, A.; Abbasi, S.W.; Alblihy, A.; Aloliqi, A.A.; Alkhayl, F.F.; Alrumaihi, F.; Ahmad, S.; El Bakri, Y. Discovery of Novel Inhibitors from Medicinal Plants for V-Domain Ig Suppressor of T-Cell Activation (VISTA). Front. Mol. Biosci. 2021, 8, 951. [Google Scholar] [CrossRef]
  18. Arif, R.; Ahmad, S.; Mustafa, G.; Mahrosh, H.S.; Ali, M.; Tahir ul Qamar, M.; Dar, H.R. Molecular Docking and Simulation Studies of Antidiabetic Agents Devised from Hypoglycemic Polypeptide-P of Momordica charantia. BioMed Res. Int. 2021, 2021, 5561129. [Google Scholar] [CrossRef]
  19. Khan, A.; Umbreen, S.; Hameed, A.; Fatima, R.; Zahoor, U.; Babar, Z.; Waseem, M.; Hussain, Z.; Rizwan, M.; Zaman, N. In Silico Mutagenesis-Based Remodelling of SARS-CoV-1 Peptide (ATLQAIAS) to Inhibit SARS-CoV-2: Structural-Dynamics and Free Energy Calculations. Interdiscip. Sci. Comput. Life Sci. 2021, 13, 521–534. [Google Scholar] [CrossRef]
  20. Khan, A.; Khan, M.; Saleem, S.; Babar, Z.; Ali, A.; Khan, A.A.; Sardar, Z.; Hamayun, F.; Ali, S.S.; Wei, D.-Q. Phylogenetic analysis and structural perspectives of RNA-dependent RNA-polymerase inhibition from SARS-CoV-2 with natural products. Interdiscip. Sci. Comput. Life Sci. 2020, 12, 335–348. [Google Scholar] [CrossRef]
  21. Alamri, M.A.; Tahir Ul Qamar, M.; Mirza, M.U.; Bhadane, R.; Alqahtani, S.M.; Muneer, I.; Froeyen, M.; Salo-Ahen, O.M. Pharmacoinformatics and molecular dynamics simulation studies reveal potential covalent and FDA-approved inhibitors of SARS-CoV-2 main protease 3CLpro. J. Biomol. Struct. Dyn. 2021, 39, 4936–4948. [Google Scholar] [CrossRef]
  22. Alamri, M.A.; ul Qamar, M.T.; Mirza, M.U.; Alqahtani, S.M.; Froeyen, M.; Chen, L.-L. Discovery of human coronaviruses pan-papain-like protease inhibitors using computational approaches. J. Pharm. Anal. 2020, 10, 546–559. [Google Scholar] [CrossRef]
  23. Mehmood, I.; Ijaz, M.; Ahmad, S.; Ahmed, T.; Bari, A.; Abro, A.; Allemailem, K.S.; Almatroudi, A.; Tahir ul Qamar, M. SARS-CoV-2: An update on genomics, risk assessment, potential therapeutics and vaccine development. Int. J. Environ. Res. Public Health 2021, 18, 1626. [Google Scholar] [CrossRef] [PubMed]
  24. Mumtaz, A.; Ashfaq, U.A.; ul Qamar, M.T.; Anwar, F.; Gulzar, F.; Ali, M.A.; Saari, N.; Pervez, M.T. MPD3: A useful medicinal plants database for drug designing. Nat. Prod. Res. 2017, 31, 1228–1236. [Google Scholar] [CrossRef] [PubMed]
  25. Diallo, B.N.t.; Glenister, M.; Musyoka, T.M.; Lobb, K.; Tastan Bishop, Ö. SANCDB: An update on South African natural compounds and their readily available analogs. J. Cheminformatics 2021, 13, 37. [Google Scholar] [CrossRef] [PubMed]
  26. Heo, L.; Park, H.; Seok, C. GalaxyRefine: Protein structure refinement driven by side-chain repacking. Nucleic Acids Res. 2013, 41, W384–W388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Goddard, T.D.; Huang, C.C.; Ferrin, T.E. Software extensions to UCSF chimera for interactive visualization of large molecular assemblies. Structure 2005, 13, 473–482. [Google Scholar] [CrossRef] [Green Version]
  28. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera—A visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef] [Green Version]
  29. 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] [Green Version]
  30. Dallakyan, S.; Olson, A.J. Small-molecule library screening by docking with PyRx. In Chemical Biology; Springer: Cham, Switzerland, 2015; pp. 243–250. [Google Scholar]
  31. Milanović, Ž.B.; Antonijević, M.R.; Amić, A.D.; Avdović, E.H.; Dimić, D.S.; Milenković, D.A.; Marković, Z.S. Inhibitory activity of quercetin, its metabolite, and standard antiviral drugs towards enzymes essential for SARS-CoV-2: The role of acid–base equilibria. Rsc Adv. 2021, 11, 2838–2847. [Google Scholar] [CrossRef]
  32. Ravindranath, P.A.; Forli, S.; Goodsell, D.S.; Olson, A.J.; Sanner, M.F. AutoDockFR: Advances in Protein-Ligand Docking with Explicitly Specified Binding Site Flexibility. PLoS Comput. Biol. 2015, 11, e1004586. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Case, D.A.; Cheatham, T.E., III; Darden, T.; Gohlke, H.; Luo, R.; Merz, K.M., Jr.; Onufriev, A.; Simmerling, C.; Wang, B.; Woods, R.J. The Amber biomolecular simulation programs. J. Comput. Chem. 2005, 26, 1668–1688. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Price, D.J.; Brooks, C.L., III. A modified TIP3P water potential for simulation with Ewald summation. J. Chem. Phys. 2004, 121, 10096–10103. [Google Scholar] [CrossRef]
  35. Roe, D.R.; Cheatham, T.E., III. PTRAJ and CPPTRAJ: Software for processing and analysis of molecular dynamics trajectory data. J. Chem. Theory Comput. 2013, 9, 3084–3095. [Google Scholar] [CrossRef] [PubMed]
  36. Sun, H.; Li, Y.; Tian, S.; Xu, L.; Hou, T. Assessing the performance of MM/PBSA and MM/GBSA methods. 4. Accuracies of MM/PBSA and MM/GBSA methodologies evaluated by various simulation protocols using PDBbind data set. Phys. Chem. Chem. Phys. 2014, 16, 16719–16729. [Google Scholar] [CrossRef]
  37. Hou, T.; Li, N.; Li, Y.; Wang, W. Characterization of domain–peptide interaction interface: Prediction of SH3 domain-mediated protein–protein interaction network in yeast by generic structure-based models. J. Proteome Res. 2012, 11, 2982–2995. [Google Scholar] [CrossRef]
  38. Chen, F.; Liu, H.; Sun, H.; Pan, P.; Li, Y.; Li, D.; Hou, T. Assessing the performance of the MM/PBSA and MM/GBSA methods. 6. Capability to predict protein–protein binding free energies and re-rank binding poses generated by protein–protein docking. Phys. Chem. Chem. Phys. 2016, 18, 22129–22139. [Google Scholar] [CrossRef]
  39. Miller, B.R., III; McGee, T.D., Jr.; Swails, J.M.; Homeyer, N.; Gohlke, H.; Roitberg, A.E. MMPBSA. py: An efficient program for end-state free energy calculations. J. Chem. Theory Comput. 2012, 8, 3314–3321. [Google Scholar] [CrossRef]
  40. Wang, Y.; Khan, A.; Chandra Kaushik, A.; Junaid, M.; Zhang, X.; Wei, D.-Q. The systematic modeling studies and free energy calculations of the phenazine compounds as anti-tuberculosis agents. J. Biomol. Struct. Dyn. 2019, 37, 4051–4069. [Google Scholar] [CrossRef]
  41. Khan, A.; Kaushik, A.C.; Ali, S.S.; Ahmad, N.; Wei, D.-Q. Deep-learning-based target screening and similarity search for the predicted inhibitors of the pathways in Parkinson’s disease. RSC Adv. 2019, 9, 10326–10339. [Google Scholar] [CrossRef] [Green Version]
  42. Alamri, M.A.; ul Qamar, M.T.; Afzal, O.; Alabbas, A.B.; Riadi, Y.; Alqahtani, S.M. Discovery of anti-MERS-CoV small covalent inhibitors through pharmacophore modeling, covalent docking and molecular dynamics simulation. J. Mol. Liq. 2021, 330, 115699. [Google Scholar] [CrossRef]
  43. Suleman, M.; ul Qamar, M.T.; Shoaib Saleem, S.A.; Ali, S.S.; Khan, H.; Akbar, F.; Khan, W.; Alblihy, A.; Alrumaihi, F.; Waseem, M. Mutational landscape of Pirin and elucidation of the impact of most detrimental missense variants that accelerate the breast cancer pathways: A computational modelling study. Front. Mol. Biosci. 2021, 8, 692835. [Google Scholar] [CrossRef] [PubMed]
  44. Durdagi, S.; ul Qamar, M.T.; Salmas, R.E.; Tariq, Q.; Anwar, F.; Ashfaq, U.A. Investigating the molecular mechanism of staphylococcal DNA gyrase inhibitors: A combined ligand-based and structure-based resources pipeline. J. Mol. Graph. Model. 2018, 85, 122–129. [Google Scholar] [CrossRef] [PubMed]
  45. Xue, L.C.; Rodrigues, J.P.; Kastritis, P.L.; Bonvin, A.M.; Vangone, A. PRODIGY: A web server for predicting the binding affinity of protein–protein complexes. Bioinformatics 2016, 32, 3676–3678. [Google Scholar] [CrossRef] [PubMed]
  46. Letko, M.; Marzi, A.; Munster, V. Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronaviruses. Nat. Microbiol. 2020, 5, 562–569. [Google Scholar] [CrossRef] [Green Version]
  47. Lan, J.; Ge, J.; Yu, J.; Shan, S.; Zhou, H.; Fan, S.; Zhang, Q.; Shi, X.; Wang, Q.; Zhang, L. Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor. Nature 2020, 581, 215–220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Engelhardt, C.; Petereit, F.; Lechtenberg, M.; Liefländer-Wulf, U.; Hensel, A. Qualitative and quantitative phytochemical characterization of Myrothamnus flabellifolia Welw. Fitoterapia 2016, 114, 69–80. [Google Scholar] [CrossRef]
  49. Fantoukh, O.I.; Dale, O.R.; Parveen, A.; Hawwal, M.F.; Ali, Z.; Manda, V.K.; Khan, S.I.; Chittiboyina, A.G.; Viljoen, A.; Khan, I.A. Safety Assessment of Phytochemicals Derived from the Globalized South African Rooibos Tea (Aspalathus linearis) through Interaction with CYP, PXR, and P-gp. J. Agric. Food Chem. 2019, 67, 4967–4975. [Google Scholar] [CrossRef]
  50. Pudumo, J.; Chaudhary, S.K.; Chen, W.; Viljoen, A.; Vermaak, I.; Veale, C.G.L. HPTLC fingerprinting of Croton gratissimus leaf extract with Preparative HPLC-MS-isolated marker compounds. S. Afr. J. Bot. 2018, 114, 32–36. [Google Scholar] [CrossRef]
  51. Chen, D.; Oezguen, N.; Urvil, P.; Ferguson, C.; Dann, S.; Savidge, T. Regulation of protein-ligand binding affinity by hydrogen bond pairing. Sci. Adv. 2016, 2, e1501240. [Google Scholar] [CrossRef] [Green Version]
  52. Chodera, J.D.; Mobley, D.L. Entropy-enthalpy compensation: Role and ramifications in biomolecular ligand recognition and design. Annu. Rev. Biophys. 2013, 42, 121–142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Tahir ul Qamar, M.; Ahmad, S.; Khan, A.; Mirza, M.U.; Ahmad, S.; Abro, A.; Chen, L.-L.; Almatroudi, A.; Wei, D.-Q. Structural probing of HapR to identify potent phytochemicals to control Vibrio cholera through integrated computational approaches. Comput. Biol. Med. 2021, 138, 104929. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Representation of the ORF8-IRF3 complex. (A) Surface representation of the ORF8-IRF3 complex, (B) interface residues of IRF3 while (C) Druggable site on IRF3 protein.
Figure 1. Representation of the ORF8-IRF3 complex. (A) Surface representation of the ORF8-IRF3 complex, (B) interface residues of IRF3 while (C) Druggable site on IRF3 protein.
Ijerph 19 00112 g001
Figure 2. Binding mode of Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside. The left panel show the surface representation of IRF3 and the binding of ligand. The right panel shows the interaction pattern and number of hydrogen bonds formed by each residue.
Figure 2. Binding mode of Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside. The left panel show the surface representation of IRF3 and the binding of ligand. The right panel shows the interaction pattern and number of hydrogen bonds formed by each residue.
Ijerph 19 00112 g002
Figure 3. Binding mode of Rutin (SANC01085). The left panel show the surface representation of IRF3 and the binding of ligand. The right panel shows the interaction pattern and number of hydrogen bonds formed by each residue.
Figure 3. Binding mode of Rutin (SANC01085). The left panel show the surface representation of IRF3 and the binding of ligand. The right panel shows the interaction pattern and number of hydrogen bonds formed by each residue.
Ijerph 19 00112 g003
Figure 4. Binding mode of Rutin (SANC00867). The left panel shows the surface representation of IRF3 and the binding of the ligand. The right panel shows the interaction pattern and number of hydrogen bonds formed by each residue.
Figure 4. Binding mode of Rutin (SANC00867). The left panel shows the surface representation of IRF3 and the binding of the ligand. The right panel shows the interaction pattern and number of hydrogen bonds formed by each residue.
Ijerph 19 00112 g004
Figure 5. Dynamic stability of IRF3-Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside, Tribuloside and Rutin complexes calculated as RMSD. The x-axis shows time in nanoseconds while y-axis show RMSD in Å.
Figure 5. Dynamic stability of IRF3-Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside, Tribuloside and Rutin complexes calculated as RMSD. The x-axis shows time in nanoseconds while y-axis show RMSD in Å.
Ijerph 19 00112 g005
Figure 6. Structural compactness of IRF3-Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside, Tribuloside and Rutin complexes calculated as Rg. The x-axis shows time in nanoseconds while y-axis show Rg in Å.
Figure 6. Structural compactness of IRF3-Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside, Tribuloside and Rutin complexes calculated as Rg. The x-axis shows time in nanoseconds while y-axis show Rg in Å.
Ijerph 19 00112 g006
Figure 7. Residual flexibility of IRF3-Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside, Tribuloside and Rutin complexes calculated as RMSF. The x-axis shows total residues while y-axis show RMSF in Å.
Figure 7. Residual flexibility of IRF3-Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside, Tribuloside and Rutin complexes calculated as RMSF. The x-axis shows total residues while y-axis show RMSF in Å.
Ijerph 19 00112 g007
Figure 8. Hydrogen bonds estimation in IRF3-Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside, Tribuloside and Rutin complexes. The x-axis shows time in nanoseconds while y-axis show hydrogen bonds count.
Figure 8. Hydrogen bonds estimation in IRF3-Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside, Tribuloside and Rutin complexes. The x-axis shows time in nanoseconds while y-axis show hydrogen bonds count.
Ijerph 19 00112 g008
Table 1. 2D structures, database ID, name and docking scores of the best hits.
Table 1. 2D structures, database ID, name and docking scores of the best hits.
2D StructureSANCDB IDNameIFD Scores
Ijerph 19 00112 i001SANC00850Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside−6.93
Ijerph 19 00112 i002SANC01085Rutin−6.37
Ijerph 19 00112 i003SANC00867Tribuloside−6.27
Table 2. MM/GBSA results for all the complexes including vdW, electrostatic and the total binding energy. All the energies are given in kcal/mol.
Table 2. MM/GBSA results for all the complexes including vdW, electrostatic and the total binding energy. All the energies are given in kcal/mol.
ComplexesvdWElectrostaticSAGBTotal
SANC00850−IRF3−34.86−27.89−5.3026.65−41.41
SANC01085−IRF3−25.20−29.92−2.9817.77−40.33
SANC00867−IRF3−36.84−35.38−3.1034.70−40.62
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Albutti, A. Rescuing the Host Immune System by Targeting the Immune Evasion Complex ORF8-IRF3 in SARS-CoV-2 Infection with Natural Products Using Molecular Modeling Approaches. Int. J. Environ. Res. Public Health 2022, 19, 112. https://doi.org/10.3390/ijerph19010112

AMA Style

Albutti A. Rescuing the Host Immune System by Targeting the Immune Evasion Complex ORF8-IRF3 in SARS-CoV-2 Infection with Natural Products Using Molecular Modeling Approaches. International Journal of Environmental Research and Public Health. 2022; 19(1):112. https://doi.org/10.3390/ijerph19010112

Chicago/Turabian Style

Albutti, Aqel. 2022. "Rescuing the Host Immune System by Targeting the Immune Evasion Complex ORF8-IRF3 in SARS-CoV-2 Infection with Natural Products Using Molecular Modeling Approaches" International Journal of Environmental Research and Public Health 19, no. 1: 112. https://doi.org/10.3390/ijerph19010112

APA Style

Albutti, A. (2022). Rescuing the Host Immune System by Targeting the Immune Evasion Complex ORF8-IRF3 in SARS-CoV-2 Infection with Natural Products Using Molecular Modeling Approaches. International Journal of Environmental Research and Public Health, 19(1), 112. https://doi.org/10.3390/ijerph19010112

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