Next Article in Journal
Longitudinal Consumption of Ergothioneine Reduces Oxidative Stress and Amyloid Plaques and Restores Glucose Metabolism in the 5XFAD Mouse Model of Alzheimer’s Disease
Next Article in Special Issue
Opportunities and Challenges for In Silico Drug Discovery at Delta Opioid Receptors
Previous Article in Journal
Rational Design by Structural Biology of Industrializable, Long-Acting Antihyperglycemic GLP-1 Receptor Agonists
Previous Article in Special Issue
In Silico Antiprotozoal Evaluation of 1,4-Naphthoquinone Derivatives against Chagas and Leishmaniasis Diseases Using QSAR, Molecular Docking, and ADME Approaches
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pharmaceuticals
Volume 15
Issue 6
10.3390/ph15060741
pharmaceuticals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Computational Investigations of Traditional Chinese Medicinal Compounds against the Omicron Variant of SARS-CoV-2 to Rescue the Host Immune System

by
Ziad Tareq Naman
1,
Salim Kadhim
2,
Zahraa J. K. Al-Isawi
3,
Christopher J. Butch
4,5,* and
Ziyad Tariq Muhseen
4,6,*
1
Department of Medical Laboratory Techniques, Al-Ma’Moon University College, Aladhamia, Baghdad 72029, Iraq
2
College of Pharmacy, University of Alkafeel, Najaf 61001, Iraq
3
Department of Pharmacology and Toxicology, College of Pharmacy, University of Kufa, Najaf 61001, Iraq
4
Department of Biomedical Engineering, College of Engineering and Applied Sciences, Nanjing University, Nanjing 210093, China
5
State Key Laboratory of Analytical Chemistry for Life Science, Jiangsu Key Laboratory of Artificial Functional Materials, Nanjing University, Nanjing 210093, China
6
Department of Pharmacy, Al-Mustaqbal University College, Hillah, Babylon 51001, Iraq
*
Authors to whom correspondence should be addressed.
Pharmaceuticals 2022, 15(6), 741; https://doi.org/10.3390/ph15060741
Submission received: 4 May 2022 / Revised: 27 May 2022 / Accepted: 31 May 2022 / Published: 13 June 2022
(This article belongs to the Special Issue In Silico Approaches in Drug Design)
Browse Figures
Versions Notes

Abstract

:
Macrodomain-I of the NSP3 (non-structural protein 3) is responsible for immune response hijacking in the SARS-CoV-2 infection known as COVID-19. In the omicron variant (B.1.1.529), this domain harbors a new mutation, V1069I, which may increase the binding of ADPr and consequently the infection severity. This macrodomain-I, due to its significant role in infection, is deemed to be an important drug target. Hence, using structural bioinformatics and molecular simulation approaches, we performed a virtual screening of the traditional Chinese medicines (TCM) database for potential anti-viral drugs. The screening of 57,000 compounds yielded the 10 best compounds with docking scores better than the control ADPr. Among the top ten, the best three hits—TCM42798, with a docking score of −13.70 kcal/mol, TCM47007 of −13.25 kcal/mol, and TCM30675 of −12.49 kcal/mol—were chosen as the best hits. Structural dynamic features were explored including stability, compactness, flexibility, and hydrogen bonding, further demonstrating the anti-viral potential of these hits. Using the MM/GBSA approach, the total binding free energy for each complex was reported to be −69.78 kcal/mol, −50.11 kcal/mol, and −47.64 kcal/mol, respectively, which consequently reflect the stronger binding and inhibitory potential of these compounds. These agents might suppress NSP3 directly, allowing the host immune system to recuperate. The current study lays the groundwork for the development of new drugs to combat SARS-CoV-2 and its variants.

Graphical Abstract

1. Introduction

Coronavirus disease-2019 (COVID-19) is defined as a disease caused by a novel coronavirus called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) [1]. SARS CoV-2 rapidly spread from person to person, and serious human disease was in recent times described in the city of Wuhan, Hubei Province, China. Coronavirus causes headache, fever, and respiratory disease, e.g., cough and shortness of breathing. SARS-CoV-2 infects humans through direct binding to host cell entry proteins (spike). 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), whereas nonstructural protein include3-chymotrypsin-like protease, papain-like protease, and RNA-dependent RNA polymerase, is encoded by the ORF region [2,3].
A SARS-CoV-2-associated global pandemic, with unceasing chaos, has now reported many variants of the virus [4]. Among the reported variants, many are particularly associated with immune response evasion, higher transmissibility, increased morbidity, and re-infection. Until now, these variants have been deemed either as variants of interest (VOIs) or variants of concern (VOCs) based on the threat they pose to public health. For instance, the alpha variant, also known as B.1.1.7, with 40–80% increased transmission, includes 69–70 del, N501Y, and P681H mutations, while beta, also called B.1.1.351, gamma, delta, and B.1.1.529 variants are classified as VOCs [5,6]. Moreover, a new version of B.1.1.7 supplemented with the E484K mutation was reported to be associated with 39 confirmed patients [7]. Moreover, B.1.351, which was reported in South Africa and harbors K417N, E484K, and N501Y variations, which were reported to increase the transmission and decrease the T cell-triggered immune response against COVID-19 infection, was recorded. In early 2021, P.1 variant was discovered in Brazil with K417T, E484K, and N501Y mutations in the RBD with transmission increased by 38% and mortality by 50% [8]. A molecular modelling study based on protein coupling deciphered the mechanism of higher infection associated with these variants of SARS-CoV-2 [9]. The binding variations were deeply discovered by demonstrating the interaction pattern and dynamic features of the complexes. In October 2020, India detected a more lethal strain of SARS-CoV-2, officially known as B.1.617.2. This variation has L452R, T478K, and P681R in the RBD, resulting in an 87 percent increase in transmissibility, an 85 percent increase in hospitalizations, and a 137 percent rise in death. Later, in India and the United Kingdom, a novel variant known as “Delta plus” was discovered with an additional mutation, K417N [10]. Similarly, a novel variant termed as Mu or B.1.621 was discovered to have spike protein alterations, some of which are shared with other VOCs. R346K, Y144T, Y145S, and insertion at position 146N are among the new mutations in this variant [11]. L452Q and F490S mutations in RBD are found in a unique VOI known as C.37 or its variant, which lowers antibody-mediated neutralization [12]. The B.1.617.1 strain with the L452R mutation in the (RBD) receptor-binding domain, on the other hand, has been attributed to decreased antibody neutralization by altering the corresponding conformational epitopes. Furthermore, the VOI Iota (lineage B.1.526) carries the mutation E484K, which confers resistance to therapeutic monoclonal antibodies and makes it more resistant to neutralization [13].
The most recent omicron strain, formally known as B.1.1.529, was discovered in November 2019 in South Africa and has 30 mutations in the spike protein including A67V, Δ69–70, T95I, G142D, Δ143–145, Δ211, L212I, ins214EPE, G339D, S371L, S373P, S375F, K417N, N440K, G446S, S477N, T478K, E484A, Q493R, G496S, Q498R, N501Y, Y505H, T547K, D614G, H655Y, N679K, P681H, N764K, D796Y, N856K, Q954H, N969K, and L981F, among which 15 mutations lie in the RBD [14]. Other proteins also reported mutations, such as NSP3 (K38R, V1069I, Δ1265, L1266I, A1892T), NSP4 (T492I), NSP5 (P132H), NSP6 (Δ105–107, A189V), and NSP12 (P323L), while I42V was reported in NSP14. Furthermore, sub-lineages, i.e., BA.1/B.1.1.529.1, BA.2/B.1.1.529.2, and BA.3/B.1.1.529.3, were reported by WHO (World Health Organization) as a threat to public health [15]. The variant is a major health concern all around the globe, and the therapeutic effectiveness of existing vaccines against it is yet unknown [16,17]. More research is needed to determine the molecular basis of pathogenicity in the omicron variant. Furthermore, new therapies against the recently emerging SARS-CoV-2 strains would require the use of advanced techniques.
The novel coronavirus-2019 genome encodes four structural proteins (S, M, E, and N), in which the S-protein gives the virus its corona-like shape, which is mainly responsible for the attachment to host cell receptors (ACE2s) or surface proteins, and 16 non-structural proteins (NSP1 to NSP16). The binding of spike protein to ACE2 (Angiotensin-Converting Enzyme 2) of the host initiates the infection in cells. The ACE2 is mainly expressed in the lungs, kidney, and small intestine, leading to serious illness [18]. During infection, the host cell protease cleaves the S-protein at the S1/S2 cleavage site. This priming (cleavage of S-protein) results in the division of protein into the S1 ectodomain at the N-terminal and the S2 membrane-anchored domain at the C-terminal. The S1 subunit recognizes the associated cell surface receptor, while the latter assists viral entry [19]. The SARS-CoV S1 subunit of the spike protein has conserved 14 aa in the RBD, which functions to recognize ACE2 and can infect both humans and bats. Among this conserved 14 aa in SARS-CoV, eight residues are highly conserved in 2019-nCoV, supporting the assumption that ACE2 is also the receptor of this new virus [20]. Among these proteins, the main proteins, PLpro, spike, and RdRp, are deemed to be direct drug and vaccine targets, but other proteins are also essential for drug design [21]. For instance, the Mac-I of the largest protein, NSP3, has also been deemed as an attractive drug target because of its role in disrupting the innate immune reaction and increasing of virulent properties of the virus [22]. NSP3 is comprised of three macrodomains and two adjacent SUD-M-like domains that have been reported to be associated with NSP3 functional modulation. Mac-I has a significant role in viral pathogenesis and increases interferons’ responses for viral neutralization [23]. Further reports have also disclosed the hijacking of the host immune response by Mac, which is carried out through interference with the IFN (Interferons) pathway and malfunctioning of STAT1 (signal transducer and activator of transcription 1) [24,25]. This hijacking of the IFN and STAT1 pathways has been reported to be possibly linked with the cytokines storm phenomena [26,27].
More research is needed to cope with this alarming pandemic situation in order to create safe and effective treatments swiftly. In this context, the macrodomain is thought to be the most druggable target for the development of COVID-19 therapy [28,29]. As a result, we used computational molecular screening to investigate the binding affinity of drugs against the macrodomain-I of the omicron (B.1.1.529) variant. The current study employs molecular docking and simulation-based methods [30,31,32] to identify potential anti-viral compounds from the TCM database against SARS-CoV-2. The findings provide crucial information on the antiviral effectiveness of the evaluated drugs against SARS-CoV-2. The findings will aid in the development and identification of potential medicinal solutions for the treatment of COVID-19.

2. Results and Discussion

2.1. Macrodomain of Omicron (B.1.1.529) and Structural Modelling

The continuously emerging variants of SARS-CoV-2 are a greater threat to public health [33]. These variants are associated with high infection, transmission, and re-infection. Moreover, some of the previous variants such as B.1.1.7, B.1.617, P.1, and A.30 have also been reported to escape the immune response. Recently, B.1.1.529, or the omicron variant, first identified in South Africa with the highest number of mutations ever, has been reported to be associated with higher transmissibility than other variants. All the proteins encoded by the genome of SARS-CoV-2 are responsible for essential roles in virulence and replication in different stages of infection [34,35,36]. Among these proteins, the largest and multi-domain protein called non-structural protein 3 or NSP3 is an essential factor for efficient translational and replication processes. NSP3 has eight functional domains that also include macrodomains, which play significant roles in viral escaping when attacked by the innate immune response. Among the macrodomains, macrodomain-I (X-domain) in particular also shares high sequence and structural conservation. Mac-I binds ADPr, and the strength of binding is directly associated with the pathogenicity level [37]. For instance, binding strength reduction or the complete loss of ADPr binding has been reported to be related to less or no infection in the IBV model [38]. Any mutation in the structure is associated with functional variations induced by conformation changes in the protein structures. For instance, K38R, V1069I, L1266I, and A1892T, as well as a deletion (Δ1266I) have recently been reported in the NSP3 protein of the omicron variant, as shown in Figure 1A. The Mac-I in the NSP3 starts from 1025–1194, and a single V1096I has been reported in this region, which may increase the binding of ADPr and consequently increase the pathogenesis by hijacking the innate immune response given in Figure 1A. Mac-I is also a potential therapeutic target; therefore, finding small molecule inhibitors for this domain might aid in the restoration of the host immunological innate response. For instance, a carbazole-based drug known as GeA-69 has been previously tested to target the human macrodomain to treat various types of cancers [39]. Similarly, structure-based approaches are being used to find drugs that target viral macrodomains [40]. New powerful therapies are needed to block Mac-I involvement in SARS-CoV-2 pathogenesis. Hence, our research combines structure-based drug discovery with molecular modelling to find new compounds that effectively target the Mac-I and salvage the innate response. The wild-type structure was used to construct a V1069I mutant of the Mac-I from the omicron variant (Figure 1B). The modelled V1069I-Mac-I, as shown in Figure 1C, demonstrated an RMSD (root mean square deviation) difference of 0.132 Å. The minimal RMSD difference demonstrated the accurate folding and distribution of secondary structural elements.

2.2. Virtual Screening and Re-Docking of TCM

Using EasyVS, an online web tool, the entire TCM database was screened against the Mac-I of the omicron variant. Prior to screening, the RO5 filter was applied, which reported 20,124 drugs that violated Lipinski’s rule of five. Among the total of 57,000 compounds, 36,876 compounds were screened against the binding site of ADPr. Among these, 14,255 compounds reported a docking score less than −4.0 kcal/mol, and 16,453 compounds reported a docking score less than −8.0 kcal/mol. Among the remaining compounds, 5221 compounds reported a docking score of less than −9.0 kcal/mol. Among the remaining 947 compounds, only 356 were reported to have docking scores greater than −9.5 kcal/mol. Compounds with docking scores greater than −9.50 kcal/mol were considered as the threshold, because for ADPr the docking score was previously reported to be −9.46 kcal/mol. The top scoring 50 compounds were selected for re-docking using AutoDock Vina. Among these 50 compounds, 10 compounds reported the best docking score and are given in Table 1. Among the top 10 compounds, TCM42798 had a docking score of −13.70 kcal/mol, TCM47007 of −13.25 kcal/mol, TCM30675 of −12.49 kcal/mol, TCM27763 of −11.93 kcal/mol, TCM33425 of −11.72 kcal/mol, TCM28788 of −11.46 kcal/mol, TCM42159 of −11.45 kcal/mol, TCM47184 of −11.36 kcal/mol, TCM31603 of −11.04 kcal/mol, and TCM31784 of −11.02 kcal/mol. Among the top 10 compounds reported to have higher docking scores than ADPr, only the top three were selected for further analysis, such as interactions, simulation, and post-simulation analysis. In addition to the top three compounds, TCM27763 and apigenin-bioside have been previously reported to have anti-adenovirus activity, which confirms the anti-viral activity of these compounds [41]. It was reported that apigenin reduces adenovirus replication and associated cellular toxicity. For instance, TCM28788 and TCM42159 have been reported to have anti-dengue viral properties, which demonstrate the anti-viral potential of these molecules [42]. The anti-herpes simplex virus type 2 infection activity of tibeticanol and the other compounds are well documented, and it has been reported that these compounds halt viral replication [43].
TCM42798, or mucic acid 1-methyl ester 2-O-gallate, was reported to have the best docking score of −13.70 kcal/mol. Mucic acid 1-methyl ester 2-O-gallate established a total of nine hydrogen bonds with Ala1060, Asn1062, Gly1068, Val1071, Ala1072, Leu1148, Ala1176, and Phe1178. For instance, these residues were previously reported to have an important role in the binding of ADPr. Comparatively, the docking score for ADPr has been reported to be −9.46 kcal/mol [22]. This shows the strongest binding of TCM42798 and consequently produces an inhibitory effect on Mac-I. The binding pattern of TCM42798 is given in Figure 2.
On the other hand, TCM47007, or (2S)-5,7,2’,5’-Tetrahydroxyflavanone 7-O--D-glucuronopyranoside, also demonstrated a good interaction profile by predicting the docking score of −13.25 kcal/mol. Unlike TCM42798, this complex reported only seven hydrogen bonds with the key residues required for interaction with ADPr. As shown in Figure 3, Ala1043, Ala1060, Gly1068, Val1071, Ala1072, and Ala1176 are involved in hydrogen bonding interactions. This shows the strongest binding of TCM47007 and consequently produces an inhibitory effect on Mac-I. Moreover, the interaction pattern for TCM30675, or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1, was also analyzed to demonstrate the binding mode of TCM30675. Residues such as Ala1060, Val1176, Phe1178, and Asp1179 are involved in direct interactions with the Mac-I. This shows the strongest binding of TCM30675 and consequently produces an inhibitory effect on Mac-I. The binding pattern of TCM42798 is given in Figure 4.

2.3. Dynamic Stability Analysis of the Top Compounds

Analysis of the drug-bound protein complexes to decipher the structural stability is a key process in determining the inhibitory potential of the interacting compound. To foresee the dynamic stability of each top-scoring compound, we calculated the root mean square deviation (RMSD) for all of the simulation trajectories. As given in Figure 5A–C, all the complexes attained equilibrium at the earlier simulation and reached stability at 1.0 Å. In the case of mucic acid 1-methyl ester 2-O-gallate, the complex demonstrated overall stable dynamics with no significant deviation, except at the start of the simulation (1–10 ns) when the structure reported minor acceptable deviation and then reached the equilibrium point with a uniform RMSD graph. The average RMSD reported for the mucic acid 1-methyl ester 2-O-gallate complex was reported to be 1.0 Å. The RMSD for mucic acid 1-methyl ester 2-O-gallate is given in Figure 5A. The RMSD for TCM47007 or (2S)-5,7,2’,5’-Tetrahydroxyflavanone 7-O--D-glucuronopyranoside reported comparatively small unstable dynamics until 15 ns, but then RMSD stabilized and attained equilibrium. As given in Figure 5B, after 15 ns the structure attained stability and demonstrated an average RMSD of 1.1 Å. Furthermore, the RMSD of TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 demonstrated more similar behavior to the mucic acid 1-methyl ester 2-O-gallate complex. The RMSD initially increased until 10 ns and then attained equilibrium. The complex reported no significant deviation except a small deviation between 10 and 20 ns, and then it stabilized again. An average RMSD of 1.2 Å was reported for TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1, as given in Figure 5C. Consequently, this shows the stronger binding and inhibitory potential of these compounds to rescue the host immune response against the COVID-19 infection.

2.4. Structural Compactness Analysis

Analysis of structural compactness in a dynamic environment is essential to understand the binding and unbinding events that happened during the simulation. These events demonstrated significant information regarding the binding stability and could be used to select the best compounds for the inhibition potential determination in the experimental setup. Thus, to determine the structural compactness, we calculated radius of gyration (Rg) over the simulation time as a function of time for each complex. As given in Figure 6A–C, all the complexes followed a similar pattern as RMSD. The TCM42798 complex demonstrated a uniform pattern of Rg with no increase or decrease over the simulation time. The structure reported no significant deviation, and the average Rg was reported to be 14.90 Å. The Rg for TCM42798 is shown in Figure 6A. On the other hand, TCM47007 demonstrated a small deviation between 23 and 28 ns and then reported a stable straight graph, which showed the stable binding steered by different kinds of bonds between TCM47007 and Mac-I. The TCM47007-Mac-I complex reported an average Rg of 14.96 Å. The Rg for TCM47007 is shown in Figure 6B. Furthermore, TCM30675 also reported a similar pattern of Rg as RMSD with no major significant deviation. A straight uniform Rg graph can be seen in Figure 6C, where a small fluctuation between 20 and 22 ns can be observed, while the Rg then stabilized and no significant deviation was reported. An average Rg of 15.90 Å was also reported for TCM30675. The Rg for TCM30675 is shown in Figure 6C.

2.5. Residues Flexibility Profiling

Residue flexibility is strongly correlated with the functional relevance of a protein, as it confers strength to the binding between the interacting molecules. Assessment of residue flexibility for the key residues demonstrates the impact of the small molecule binding. Thus, to determine how the residue flexibility is affected by the binding of these drugs, we calculated the root mean square fluctuation (RMSF) for each complex. As given in Figure 7, all the complexes demonstrated a similar pattern of residue flexibility. All the residues in each complex displayed minimal fluctuation, which showed the stronger binding of these compounds. The regions between 40 and 50, 90 and 100, and 120 and 140 demonstrated comparatively higher fluctuations due to the loop distribution.

2.6. Hydrogen Bonding Analysis

Hydrogen bonds confer strength to the binding of small molecules. The binding of small molecules is steered by hydrogen and many other bonds, which consequently show the inhibitory features. To determine the binding strength, we calculated the total number of hydrogen bonds and the population, and the simulation trajectories were analyzed. In TCM42798 or mucic acid 1-methyl ester 2-O-gallate, the average number of hydrogen bonds was reported to be 82, while in the case of TCM47007 or (2S)-5,7,2’,5’-Tetrahydroxyflavanone 7-O--D-glucuronopyranoside, the average number of hydrogen bonds was also reported to be 79, while TCM30675 or the (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 complex reported 78 average hydrogen bonds. The total number of hydrogen bonds in each complex is shown in Figure 8A–C. Moreover, the hydrogen bonding population for the key interacting residues was estimated and revealed in each complex, namely Ala1060 (63%, 74%, 58%), Asn1062 (66%, 53%, 61%), Gly1068 (51%, 43%, 34%) Val1071 (32%, 24%, 23%), Ala1072 (4%, 17%, 13%), Leu1148 (1%, 6%, 12%), Ala1176 (0%, 0%, 7%), and Phe1178 (0.36%, 2%, 11%) in TCM42798, TCM47007, and TCM30675 complexes, respectively. Together these results show that these compounds bind more strongly to the Mac-I and thus produce inhibitory properties that could help to rescue the host immune response against COVID-19 infection.

2.7. Binding Free Energy Estimation

The binding strength of small molecules using the binding free energy method, MM-GBSA, is a widely used method for re-demonstrating docking stability and correct binding. The above-mentioned MM-GBSA technique is more computationally economical than the more expensive alchemical free energy method. When compared to rational scoring functions, it is one of the most accurate techniques. Keeping in mind the implementation of this approach, we also employed the same method to compute the binding free energy for TCM42798, TCM47007, and TCM30675 complexes. The vdWs for these complexes were reported to be −84.26 kcal/mol, −59.79 kcal/mol, and −53.24 kcal/mol, respectively, while the electrostatic energies were reported to be −12.22 kcal/mol, −13.22 kcal/mol, and −15.66 kcal/mol, respectively. Moreover, the ESURF values were reported to be 17.45 kcal/mol, 14.68 kcal/mol, and 12.25 kcal/mol, respectively. The total binding free energy for these complexes was reported to be −69.78 kcal/mol, −50.11 kcal/mol, and −47.64 kcal/mol, respectively, which consequently reflected the stronger binding and inhibitory potential of these compounds. The MM/GBSA results are given in Table 2.

3. Material and Methods

3.1. Modelling of the Macrodomain-I (Mac-I) of B.1.1.529 Variant

The experimentally reported structure in the protein database was collected through accession number 6W02 to model the V1069I mutations reported in the Mac-I of NSP3 of the omicron variant [44]. For this purpose, the amino acid sequence using P0DTD1 accession number was obtained from UniProt, European union [45]. For modelling of the variant structure, AlphaFold 2.0 Seoul National University, Seoul, South Korea, was used, which is currently the best and most accurate approach for 3D structural modelling [46].

3.2. Virtual Screening of Traditional Chinese Medicine Database

The Traditional Chinese Medicine Database (TCM) currently holds 57,000 entries of different medicinal compounds isolated from Chinese herbs. It is considered as an important medicinal repository for discovering novel treatments that are safe and effective. Thus, we screened the complete TCM database against Mac-I of the B.1.1.529 variant. For screening, we used the EasyVS (http://biosig.unimelb.edu.au/easyvs, accessed on 10 April 2022) webserver [47]. Custom parameters were used and defined, such as the active site, and the RO5 filter was enabled for the filtration and removal of molecules that violated Lipinski’s rules. Furthermore, for the top 100 hits, a second round of screening was performed using the Auto Dock Vina algorithm [48]. Finally, the top-selected compounds were used for the Induced-Fit docking (IFD) approach to remove the false-positive results. The top hits were identified based on the docking score for ADPr, taken as a control, as previously reported [22].

3.3. Molecular Dynamics Simulation (MDS)

Characterization of dynamic features to accurately portray the inhibitory potential of the top hits was studied by performing molecular dynamics simulation [49]. The top three hits were investigated by using AMBER20 simulation software by adding an optimal point charge water model and sodium ions for neutralizing the effect of any charge. Protein and drug parameterization was achieved by recruiting FF19SB and GAFF forcefields. Each system after gentle minimization was subjected to heating followed by equilibration. The production runs for 50 ns each were completed. For any long-range electrostatic interactions (10.0 Å cutoff), the PME (particle mesh Ewald algorithm) was used, while covalent bonds were treated with the SHAKE algorithm.

3.4. Trajectories Analysis Using CPPTRAJ and PTRAJ

Evaluation of the simulation trajectories to forecast the dynamic stability, flexibility index, hydrogen bonding, structural compactness, and other parameters was completed by using CPPTRAJ and PTRAJ modules [50]. Root mean square deviation (RMSD) was estimated for stability, root mean square fluctuation (RMSF) to index the flexibility of each residue, radius of gyration (Rg) to foresee the compactness, and hydrogen bonding analysis to estimate the bonding population.

3.5. Estimation of Post-Simulation Binding Energy

The most generally utilized strategy in many related studies is to evaluate the strength of small molecule binding by employing the binding free energy (BFE) approach [51,52,53,54]. We also adopted the MMPBSA.py script to compute the binding free energy of the protein–ligand complexes by evaluating 2500 snapshots, keeping in mind the relevance of this strategy in re-ranking the binding conformations. For this purpose, the following equation was used to estimate the BFE:
“ Δ G b i n d = Δ G c o m p l e x − [ Δ G r e c e p t o r + Δ G l i g a n d ] ”  
The ΔGbind represents the total binding energy, while ΔGreceptor, ΔGligand, and ΔGcomplex represent the binding energy of protein, drug, and complex, respectively. The following equation was used to estimate individual binding energies such as bonded (Gbond), electrostatic (Gele), polar (Gpol), and non-polar (Gnpol), which contribute to the total binding free energy.
“ G = G b o n d G e l e c t r o s t a t i c − G v a n   d e r   W a a l − G p o l a r G n o n − p o l a r ”

4. Conclusions

The present study employs molecular modelling and MD simulation approaches to target the macrodomain-I of the B.1.1.529 variant of SARS-CoV-2. Three novel compounds, namely TCM42798, TCM47007, and TCM30675, were identified as potential inhibitors of Mac-I. These agents may suppress NSP3 directly, allowing the host immune system to recuperate. The current study lays the groundwork for the development of new drugs to combat SARS-CoV-2 and its variants.

Author Contributions

Conceptualization, C.J.B. and Z.T.M.; Data curation, Z.T.N.; S.K. and Z.J.K.A.-I.; Formal analysis, Z.T.N. and Z.T.M.; Funding acquisition, C.J.B.; Investigation, Z.T.N.; Methodology, S.K. and Z.J.K.A.-I.; Project administration, C.J.B. and Z.T.M.; Resources, C.J.B. and Z.T.M.; Software, Z.T.M.; Supervision, C.J.B.; Validation, S.K.; Z.J.K.A.-I. and C.J.B.; Visualization, Z.T.N. and Z.T.M.; Writing—original draft, Z.T.N. and Z.T.M.; Writing—review and editing, S.K.; Z.J.K.A.-I. and C.J.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research work was supported by the Nanjing University startup research funds.

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 authors would like to acknowledge Nanjing University, China, for providing facilities for this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 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]
  2. V’kovski, P.; Kratzel, A.; Steiner, S.; Stalder, H.; Thiel, V.J.N.R.M. Coronavirus biology and replication: Implications for SARS-CoV-2. Nat. Rev. Microbiol. 2021, 19, 155–170. [Google Scholar] [CrossRef]
  3. Fontanet, A.; Autran, B.; Lina, B.; Kieny, M.P.; Karim, S.S.A.; Sridhar, D.J.T.L. SARS-CoV-2 variants and ending the COVID-19 pandemic. Lancet 2021, 397, 952–954. [Google Scholar] [CrossRef]
  4. Tao, K.; Tzou, P.L.; Nouhin, J.; Gupta, R.K.; de Oliveira, T.; Kosakovsky Pond, S.L.; Fera, D.; Shafer, R.W. The biological and clinical significance of emerging SARS-CoV-2 variants. Nat. Rev. Genet. 2021, 22, 757–773. [Google Scholar] [CrossRef] [PubMed]
  5. Duong, D. Alpha, Beta, Delta, Gamma: What’s important to know about SARS-CoV-2 variants of concern? Can. Med. Assoc. 2021, 193, E1059–E1060. [Google Scholar] [CrossRef]
  6. Meng, B.; Kemp, S.A.; Papa, G.; Datir, R.; Ferreira, I.A.; Marelli, S.; Harvey, W.T.; Lytras, S.; Mohamed, A.; Gallo, G. Recurrent emergence of SARS-CoV-2 spike deletion H69/V70 and its role in the Alpha variant B. 1.1. 7. Cell Rep. 2021, 35, 109292. [Google Scholar] [CrossRef]
  7. Han, P.; Su, C.; Zhang, Y.; Bai, C.; Zheng, A.; Qiao, C.; Wang, Q.; Niu, S.; Chen, Q.; Zhang, Y. Molecular insights into receptor binding of recent emerging SARS-CoV-2 variants. Nat. Commun. 2021, 12, 6103. [Google Scholar] [CrossRef]
  8. Lauring, A.S.; Malani, P.N. Variants of SARS-CoV-2. JAMA 2021, 326, 880. [Google Scholar] [CrossRef]
  9. Khan, A.; Waris, H.; Rafique, M.; Suleman, M.; Mohammad, A.; Ali, S.S.; Khan, T.; Waheed, Y.; Liao, C.; Wei, D.Q. The Omicron (B.1.1.529) variant of SARS-CoV-2 binds to the hACE2 receptor more strongly and escapes the antibody response: Insights from structural and simulation data. Int. J. Biol. Macromol. 2022, 200, 438–448. [Google Scholar] [CrossRef]
  10. Mlcochova, P.; Kemp, S.A.; Dhar, M.S.; Papa, G.; Meng, B.; Ferreira, I.A.; Datir, R.; Collier, D.A.; Albecka, A.; Singh, S. SARS-CoV-2 B. 1.617. 2 Delta variant replication and immune evasion. Nature 2021, 599, 114–119. [Google Scholar] [CrossRef]
  11. Messali, S.; Bertelli, A.; Campisi, G.; Zani, A.; Ciccozzi, M.; Caruso, A.; Caccuri, F. A cluster of the new SARS-CoV-2 B. 1.621 lineage in Italy and sensitivity of the viral isolate to the BNT162b2 vaccine. J. Med. Virol. 2021, 93, 6468–6470. [Google Scholar] [CrossRef]
  12. Romero, P.E.; Dávila-Barclay, A.; Salvatierra, G.; González, L.; Cuicapuza, D.; Solís, L.; Marcos-Carbajal, P.; Huancachoque, J.; Maturrano, L.; Tsukayama, P. The emergence of SARS-CoV-2 variant lambda (C. 37) in South America. Microbiol. Spectr. 2021, 9, e00789-21. [Google Scholar] [CrossRef] [PubMed]
  13. Khan, A.; Khan, M.T.; Saleem, S.; Junaid, M.; Ali, A.; Ali, S.S.; Khan, M.; Wei, D.-Q. Structural Insights into the mechanism of RNA recognition by the N-terminal RNA-binding domain of the SARS-CoV-2 nucleocapsid phosphoprotein. Comput. Struct. Biotechnol. J. 2020, 18, 2174–2184. [Google Scholar] [CrossRef]
  14. He, X.; Hong, W.; Pan, X.; Lu, G.; Wei, X. SARS-CoV-2 Omicron variant: Characteristics and prevention. MedComm 2021, 2, 838–845. [Google Scholar] [CrossRef]
  15. Karim, S.S.A.; Karim, Q.A. Omicron SARS-CoV-2 variant: A new chapter in the COVID-19 pandemic. Lancet 2021, 398, 2126–2128. [Google Scholar] [CrossRef]
  16. Wilhelm, A.; Widera, M.; Grikscheit, K.; Toptan, T.; Schenk, B.; Pallas, C.; Metzler, M.; Kohmer, N.; Hoehl, S.; Helfritz, F.A. Reduced neutralization of SARS-CoV-2 Omicron variant by vaccine sera and monoclonal antibodies. MedRxiv 2021. [Google Scholar] [CrossRef]
  17. Roessler, A.; Riepler, L.; Bante, D.; von Laer, D.; Kimpel, J. SARS-CoV-2 B. 1.1. 529 variant (Omicron) evades neutralization by sera from vaccinated and convalescent individuals. MedRxiv 2021. [Google Scholar] [CrossRef]
  18. Duchene, S.; Featherstone, L.; Haritopoulou-Sinanidou, M.; Rambaut, A.; Lemey, P.; Baele, G. Temporal signal and the phylodynamic threshold of SARS-CoV-2. Virus Evol. 2020, 6, veaa061. [Google Scholar] [CrossRef] [PubMed]
  19. Wrapp, D.; Wang, N.; Corbett, K.S.; Goldsmith, J.A.; Hsieh, C.L.; Abiona, O.; Graham, B.S.; McLellan, J.S. Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science 2020, 367, 1260–1263. [Google Scholar] [CrossRef] [Green Version]
  20. Li, W.; Moore, M.J.; Vasilieva, N.; Sui, J.; Wong, S.K.; Berne, M.A.; Somasundaran, M.; Sullivan, J.L.; Luzuriaga, K.; Greenough, T.C. Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus. Nature 2003, 426, 450–454. [Google Scholar] [CrossRef] [Green Version]
  21. 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]
  22. Hussain, I.; Pervaiz, N.; Khan, A.; Saleem, S.; Shireen, H.; Wei, D.-Q.; Labrie, V.; Bao, Y.; Abbasi, A.A. Evolutionary and structural analysis of SARS-CoV-2 specific evasion of host immunity. Genes Immun. 2020, 21, 409–419. [Google Scholar] [CrossRef] [PubMed]
  23. Brosey, C.A.; Houl, J.H.; Katsonis, P.; Balapiti-Modarage, L.P.F.; Bommagani, S.; Arvai, A.; Moiani, D.; Bacolla, A.; Link, T.; Warden, L.S.; et al. Targeting SARS-CoV-2 Nsp3 macrodomain structure with insights from human poly(ADP-ribose) glycohydrolase (PARG) structures with inhibitors. Prog. Biophys. Mol. Biol. 2021, 163, 171–186. [Google Scholar] [CrossRef] [PubMed]
  24. Lin, M.-H.; Chang, S.-C.; Chiu, Y.-C.; Jiang, B.-C.; Wu, T.-H.; Hsu, C.-H. Structural, biophysical, and biochemical elucidation of the SARS-CoV-2 nonstructural protein 3 macro domain. ACS Infect. Dis. 2020, 6, 2970–2978. [Google Scholar] [CrossRef]
  25. Srinivasan, S.; Cui, H.; Gao, Z.; Liu, M.; Lu, S.; Mkandawire, W.; Narykov, O.; Sun, M.; Korkin, D. Structural genomics of SARS-CoV-2 indicates evolutionary conserved functional regions of viral proteins. Viruses 2020, 12, 360. [Google Scholar] [CrossRef] [Green Version]
  26. Hoch, N.C. Host ADP-ribosylation and the SARS-CoV-2 macrodomain. Biochem. Soc. Trans. 2021, 49, 1711–1721. [Google Scholar] [CrossRef]
  27. Molaei, S.; Dadkhah, M.; Asghariazar, V.; Karami, C.; Safarzadeh, E. The immune response and immune evasion characteristics in SARS-CoV, MERS-CoV, and SARS-CoV-2: Vaccine design strategies. Int. Immunopharmacol. 2021, 92, 107051. [Google Scholar] [CrossRef] [PubMed]
  28. Li, F. Structure, function, and evolution of coronavirus spike proteins. Annu. Rev. Virol. 2016, 3, 237–261. [Google Scholar] [CrossRef] [Green Version]
  29. Li, F.; Li, W.; Farzan, M.; Harrison, S.C. Structure of SARS coronavirus spike receptor-binding domain complexed with receptor. Science 2005, 309, 1864–1868. [Google Scholar] [CrossRef] [PubMed]
  30. Islam, S.; Hosen, M.A.; Ahmad, S.; ul Qamar, M.T.; Dey, S.; Hasan, I.; Fujii, Y.; Ozeki, Y.; Kawsar, S.M. Synthesis, antimicrobial, anticancer activities, PASS prediction, molecular docking, molecular dynamics and pharmacokinetic studies of designed methyl α-D-glucopyranoside esters. J. Mol. Struct. 2022, 1260, 132761. [Google Scholar] [CrossRef]
  31. Tahir ul Qamar, M.; Mirza, M.U.; Song, J.-M.; Rao, M.J.; Zhu, X.; Chen, L.-L. Probing the structural basis of Citrus phytochrome B using computational modelling and molecular dynamics simulation approaches. J. Mol. Liq. 2021, 340, 116895. [Google Scholar] [CrossRef]
  32. Ahmad, F.; Albutti, A.; Tariq, M.H.; Din, G.; Tahir ul Qamar, M.; Ahmad, S. Discovery of Potential Antiviral Compounds against Hendra Virus by Targeting Its Receptor-Binding Protein (G) Using Computational Approaches. Molecules 2022, 27, 554. [Google Scholar] [CrossRef]
  33. Abbas, K.; Khan, T.; Ali, S.; Aftab, S.; Wang, Y.; Qiankun, W.; Khan, M. SARS-CoV-2 new variants: Characteristic features and impact on the efficacy of different vaccines. Biomed. Pharmacother. 2021, 143, 112176. [Google Scholar]
  34. Magrane, M. UniProt Knowledgebase: A hub of integrated protein data. Database 2011, 2011, bar009. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Žídek, A.; Potapenko, A. Highly accurate protein structure prediction with AlphaFold. Nature 2021, 596, 583–589. [Google Scholar] [CrossRef]
  36. Pires, D.E.V.; Veloso, W.N.P.; Myung, Y.; Rodrigues, C.H.M.; Silk, M.; Rezende, P.M.; Silva, F.; Xavier, J.S.; Velloso, J.P.L.; da Silveira, C.H.; et al. EasyVS: A user-friendly web-based tool for molecule library selection and structure-based virtual screening. Bioinformatics 2020, 36, 4200–4202. [Google Scholar] [CrossRef]
  37. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef] [Green Version]
  38. 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, 716735. [Google Scholar] [CrossRef]
  39. 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]
  40. ul Qamar, M.T.; 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]
  41. 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] [PubMed]
  42. 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]
  43. 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]
  44. 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]
  45. 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]
  46. Muhseen, Z.T.; Hameed, A.R.; Al-Hasani, H.M.; ul Qamar, M.T.; Li, G. Promising terpenes as SARS-CoV-2 spike receptor-binding domain (RBD) attachment inhibitors to the human ACE2 receptor: Integrated computational approach. J. Mol. Liq. 2020, 320, 114493. [Google Scholar] [CrossRef] [PubMed]
  47. 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] [PubMed]
  48. Han, W.; Li, X.; Fu, X. The macro domain protein family: Structure, functions, and their potential therapeutic implications. Mutat. Res./Rev. Mutat. Res. 2011, 727, 86–103. [Google Scholar] [CrossRef] [PubMed]
  49. Leung, A.K.; McPherson, R.L.; Griffin, D.E. Macrodomain ADP-ribosylhydrolase and the pathogenesis of infectious diseases. PLoS Pathog. 2018, 14, e1006864. [Google Scholar] [CrossRef]
  50. Schuller, A.P.; Wu, C.C.-C.; Dever, T.E.; Buskirk, A.R.; Green, R. eIF5A functions globally in translation elongation and termination. Mol. Cell 2017, 66, 194–205.e5. [Google Scholar] [CrossRef] [Green Version]
  51. Babar, Z.; Khan, M.; Zahra, M.; Anwar, M.; Noor, K.; Hashmi, H.F.; Suleman, M.; Waseem, M.; Shah, A.; Ali, S. Drug similarity and structure-based screening of medicinal compounds to target macrodomain-I from SARS-CoV-2 to rescue the host immune system: A molecular dynamics study. J. Biomol. Struct. Dyn. 2020, 40, 523–537. [Google Scholar] [CrossRef] [PubMed]
  52. Kanerva, A.; Raki, M.; Ranki, T.; Särkioja, M.; Koponen, J.; Desmond, R.A.; Helin, A.; Stenman, U.H.; Isoniemi, H.; Höckerstedt, K.; et al. Chlorpromazine and apigenin reduce adenovirus replication and decrease replication associated toxicity. J. Gene Med. 2007, 9, 3–9. [Google Scholar] [CrossRef] [PubMed]
  53. Shimu, M.S.S.; Mahmud, S.; Tallei, T.E.; Sami, S.A.; Adam, A.A.; Acharjee, U.K.; Paul, G.K.; Emran, T.B.; Zaman, S.; Uddin, M.S.; et al. Phytochemical Compound Screening to Identify Novel Small Molecules against Dengue Virus: A Docking and Dynamics Study. Molecules 2022, 27, 653. [Google Scholar] [CrossRef] [PubMed]
  54. Zhang, C.-J.; Li, W.; Li, H.-Y.; Wang, Y.-L.; Yun, T.; Song, Z.-P.; Song, Y.; Zhao, X.-W. In vivo and in vitro antiviral activity of five Tibetan medicinal plant extracts against herpes simplex virus type 2 infection. Pharm. Biol. 2009, 47, 598–607. [Google Scholar] [CrossRef]
Pharmaceuticals 15 00741 g001 550
Figure 1. (A) NSP3 and reported mutations in the omicron variant. (B) Surface representation of Mac-I in complex with ADPr. (C) Superimposed wild type a V1069I mutation in Mac-I.
Figure 1. (A) NSP3 and reported mutations in the omicron variant. (B) Surface representation of Mac-I in complex with ADPr. (C) Superimposed wild type a V1069I mutation in Mac-I.
Pharmaceuticals 15 00741 g001
Pharmaceuticals 15 00741 g002 550
Figure 2. Binding of TCM42798 in the cavity of Mac-I. The right surface representation shows the binding mode of mucic acid 1-methyl ester 2-O-gallate (yellow sticks) inside the binding cavity. The right panel represent the 3D interaction pattern, where the TCM42798 is shown as yellow sticks, while the interacting residues are given in green sticks. The interactions (hydrogen bonds) are shown in pink color.
Figure 2. Binding of TCM42798 in the cavity of Mac-I. The right surface representation shows the binding mode of mucic acid 1-methyl ester 2-O-gallate (yellow sticks) inside the binding cavity. The right panel represent the 3D interaction pattern, where the TCM42798 is shown as yellow sticks, while the interacting residues are given in green sticks. The interactions (hydrogen bonds) are shown in pink color.
Pharmaceuticals 15 00741 g002
Pharmaceuticals 15 00741 g003 550
Figure 3. Binding of TCM47007 in the cavity of Mac-I. The right surface representation shows the binding mode of TCM47007 or (2S)-5,7,2’,5’-Tetrahydroxyflavanone 7-O--D-glucuronopyranoside (orange sticks) inside the binding cavity. The right panel represents the 3D interaction pattern, where the TCM47007 or (2S)-5,7,2’,5’-Tetrahydroxyflavanone 7-O--D-glucuronopyranoside is shown as orange sticks, while the interacting residues are given in green sticks. The interactions (hydrogen bonds) are shown in blue color.
Figure 3. Binding of TCM47007 in the cavity of Mac-I. The right surface representation shows the binding mode of TCM47007 or (2S)-5,7,2’,5’-Tetrahydroxyflavanone 7-O--D-glucuronopyranoside (orange sticks) inside the binding cavity. The right panel represents the 3D interaction pattern, where the TCM47007 or (2S)-5,7,2’,5’-Tetrahydroxyflavanone 7-O--D-glucuronopyranoside is shown as orange sticks, while the interacting residues are given in green sticks. The interactions (hydrogen bonds) are shown in blue color.
Pharmaceuticals 15 00741 g003
Pharmaceuticals 15 00741 g004 550
Figure 4. Binding of TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 in the cavity of Mac-I. The right surface representation shows the binding mode of TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 (blue sticks) inside the binding cavity. The right panel represent the 3D interaction pattern, where the TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 is shown as orange sticks, while the interacting residues are given in green sticks. The interactions (hydrogen bonds) are shown in blue color.
Figure 4. Binding of TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 in the cavity of Mac-I. The right surface representation shows the binding mode of TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 (blue sticks) inside the binding cavity. The right panel represent the 3D interaction pattern, where the TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 is shown as orange sticks, while the interacting residues are given in green sticks. The interactions (hydrogen bonds) are shown in blue color.
Pharmaceuticals 15 00741 g004
Pharmaceuticals 15 00741 g005 550
Figure 5. Dynamic stability analysis of the TCM42798, TCM47007, and TCM30675 complexes. (A) The RMSD graph for mucic acid 1-methyl ester 2-O-gallate. (B) The RMSD graph for TCM47007 or (2S)-5,7,2’,5’-Tetrahydroxyflavanone 7-O--D-glucuronopyranoside. (C) The RMSD graph for TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 complex.
Figure 5. Dynamic stability analysis of the TCM42798, TCM47007, and TCM30675 complexes. (A) The RMSD graph for mucic acid 1-methyl ester 2-O-gallate. (B) The RMSD graph for TCM47007 or (2S)-5,7,2’,5’-Tetrahydroxyflavanone 7-O--D-glucuronopyranoside. (C) The RMSD graph for TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 complex.
Pharmaceuticals 15 00741 g005
Pharmaceuticals 15 00741 g006 550
Figure 6. Structural compactness analysis in a dynamic environment of the TCM42798, TCM47007, and TCM30675 complexes calculated as Rg. (A) The Rg graph for mucic acid 1-methyl ester 2-O-gallate. (B) The Rg graph for TCM47007 or (2S)-5,7,2’,5’-Tetrahydroxyflavanone 7-O--D-glucuronopyranoside. (C) The Rg graph for TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 complex.
Figure 6. Structural compactness analysis in a dynamic environment of the TCM42798, TCM47007, and TCM30675 complexes calculated as Rg. (A) The Rg graph for mucic acid 1-methyl ester 2-O-gallate. (B) The Rg graph for TCM47007 or (2S)-5,7,2’,5’-Tetrahydroxyflavanone 7-O--D-glucuronopyranoside. (C) The Rg graph for TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 complex.
Pharmaceuticals 15 00741 g006
Pharmaceuticals 15 00741 g007 550
Figure 7. Residual flexibility analysis of the TCM42798, TCM47007, and TCM30675 complexes.
Figure 7. Residual flexibility analysis of the TCM42798, TCM47007, and TCM30675 complexes.
Pharmaceuticals 15 00741 g007
Pharmaceuticals 15 00741 g008 550
Figure 8. Hydrogen bonding analysis of the TCM42798, TCM47007, and TCM30675 complexes. (A) The H-bonds graph for mucic acid 1-methyl ester 2-O-gallate. (B) The H-bonds graph for TCM47007 or (2S)-5,7,2’,5’-Tetrahydroxyflavanone 7-O--D-glucuronopyranoside. (C) The H-bonds graph for TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 complex.
Figure 8. Hydrogen bonding analysis of the TCM42798, TCM47007, and TCM30675 complexes. (A) The H-bonds graph for mucic acid 1-methyl ester 2-O-gallate. (B) The H-bonds graph for TCM47007 or (2S)-5,7,2’,5’-Tetrahydroxyflavanone 7-O--D-glucuronopyranoside. (C) The H-bonds graph for TCM30675 or (5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-1 complex.
Pharmaceuticals 15 00741 g008
Table
Table 1. List of top 11 scoring compounds with their TCM IDs, compounds names, 3D structures, and docking scores.
Table 1. List of top 11 scoring compounds with their TCM IDs, compounds names, 3D structures, and docking scores.
TCM Database IDCompound Names2D StructuresDocking Scores
TCM42798Mucic_acid_1-methyl_ester_2-O-gallate Pharmaceuticals 15 00741 i001−13.70
TCM47007(2S)-5,7,2’,5’-Tetrahydroxyflavanone_7-O- -D-glucuronopyranoside Pharmaceuticals 15 00741 i002−13.25
TCM30675(5S,6S,7S,8R)-5,6,7,8-Tetrahydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethyl]-5,6,7,8-tetrahydro-4H-chromen-4-one Pharmaceuticals 15 00741 i003−12.49
TCM2776330389 Pharmaceuticals 15 00741 i004−11.93
TCM33425Apigenin-bioside Pharmaceuticals 15 00741 i005−11.72
TCM2878831943 Pharmaceuticals 15 00741 i006−11.46
TCM42159(-)-5’-Methoxyisolariciresinol-2-O-D-xylopyranoside_(D2) Pharmaceuticals 15 00741 i007−11.45
TCM47184Tibeticanol Pharmaceuticals 15 00741 i008−11.36
TCM3160336132 Pharmaceuticals 15 00741 i009−11.04
TCM3178436381 Pharmaceuticals 15 00741 i010−11.02
Table
Table 2. Binding free energy calculated as MM/GBSA. All the values are given in kcal/mol.
Table 2. Binding free energy calculated as MM/GBSA. All the values are given in kcal/mol.
MM/GBSATCM42798TCM47007TCM30675
vdW−84.26−59.79−53.24
electrostatic−12.22−13.22−15.66
ESURF17.4514.6812.25
EGB9.258.229.01
∆G Bind−69.78−50.11−47.64
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Naman, Z.T.; Kadhim, S.; Al-Isawi, Z.J.K.; Butch, C.J.; Muhseen, Z.T. Computational Investigations of Traditional Chinese Medicinal Compounds against the Omicron Variant of SARS-CoV-2 to Rescue the Host Immune System. Pharmaceuticals 2022, 15, 741. https://doi.org/10.3390/ph15060741

AMA Style

Naman ZT, Kadhim S, Al-Isawi ZJK, Butch CJ, Muhseen ZT. Computational Investigations of Traditional Chinese Medicinal Compounds against the Omicron Variant of SARS-CoV-2 to Rescue the Host Immune System. Pharmaceuticals. 2022; 15(6):741. https://doi.org/10.3390/ph15060741

Chicago/Turabian Style

Naman, Ziad Tareq, Salim Kadhim, Zahraa J. K. Al-Isawi, Christopher J. Butch, and Ziyad Tariq Muhseen. 2022. "Computational Investigations of Traditional Chinese Medicinal Compounds against the Omicron Variant of SARS-CoV-2 to Rescue the Host Immune System" Pharmaceuticals 15, no. 6: 741. https://doi.org/10.3390/ph15060741

APA Style

Naman, Z. T., Kadhim, S., Al-Isawi, Z. J. K., Butch, C. J., & Muhseen, Z. T. (2022). Computational Investigations of Traditional Chinese Medicinal Compounds against the Omicron Variant of SARS-CoV-2 to Rescue the Host Immune System. Pharmaceuticals, 15(6), 741. https://doi.org/10.3390/ph15060741

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