In silico investigation of potential inhibitors to main protease and spike protein of SARS-CoV-2 in propolis

Docking analysis of propolis's natural compound was successfully performed against SARS-CoV-2 main protease (Mpro) and spike protein subunit 2 (S2). Initially, the propolis's protein was screened using chromatography analysis and successfully identified 22 compounds in the propolis. Four compounds were further investigated, i.e., neoblavaisoflavone, methylophiopogonone A, 3′-Methoxydaidzin, and genistin. The binding affinity of 3′-Methoxydaidzin was ?7.7 kcal/mol, which is similar to nelfinavir (control), while the others were ?7.6 kcal/mol. However, we found the key residue of Glu A:166 in the methylophiopogonone A and genistin, even though the predicted binding energy slightly higher than nelfinavir. In contrast, the predicted binding affinity of neoblavaisoflavone, methylophiopogonone A, 3′-Methoxydaidzin, and genistin against S2 were ?8.1, ?8.2, ?8.3, and ?8.3 kcal/mol, respectively, which is far below of the control (pravastatin, ?7.3 kcal/mol). Instead of conventional hydrogen bonding, the π bonding influenced the binding affinity against S2. The results reveal that this is the first report about methylophiopogonone A, 3′-Methoxydaidzin, and genistin as candidates for anti-viral agents. Those compounds can then be further explored and used as a parent backbone molecule to develop a new supplementation for preventing SARS-CoV-2 infections during COVID-19 outbreaks.     

An in silico approach unveils the potential of antiviral compounds in preclinical and clinical trials as SARS-CoV-2 omicron inhibitors

The increased transmissibility and highly infectious nature of the new variant of concern (VOC) that is severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron and lack of effective therapy need the rapid discovery of therapeutic antivirals against it. The present investigation aimed to identify antiviral compounds that would be effective against SARS-CoV-2 Omicron. In this study, molecular docking experiments were carried out using the recently reported experimental structure of omicron spike protein in complex with human angiotensin-converting enzyme 2 (ACE2) and various antivirals in preclinical and clinical trial studies. Out of 36 tested compounds, Abemaciclib, Dasatinib and Spiperone are the three top-ranked molecules which scored binding energies of ?10.08 kcal/mol, ?10.06 kcal/mol and ?9.54 kcal/mol respectively. Phe338, Asp339, and Asp364 are crucial omicron receptor residues involved in hydrogen bond interactions, while other residues were mostly involved in hydrophobic interactions with the lead molecules. The identified lead compounds also scored well in terms of drug-likeness. Molecular dynamics (MD) simulation, essential dynamics (ED) and entropic analysis indicate the ability of these molecules to modulate the activity of omicron spike protein. Therefore, Abemaciclib, Dasatinib and Spiperone are likely to be viable drug-candidate molecules that can block the interaction between the omicron spike protein and the host cellular receptor ACE2. Though our findings are compelling, more research into these molecules is needed before they can be employed as drugs to treat SARS-CoV-2 omicron infections.     

Determination of the interaction between the receptor binding domain of 2019-nCoV spike protein,TMPRSS2, cathepsin B and cathepsin L,and glycosidic and aglycon forms of some flavonols

The novel coronavirus (COVID-19, SARS-CoV-2) is a rapidly spreading disease with a high mortality. In this research, the interactions between specific flavonols and the 2019-nCoV receptor binding domain (RBD), transmembrane protease, serine 2 (TMPRSS2), and cathepsins (CatB and CatL) were analyzed. According to the relative binding capacity index (RBCI) calculated based on the free energy of binding and calculated inhibition constants, it was determined that robinin (ROB) and gossypetin (GOS) were the most effective flavonols on all targets. While the binding free energy of ROB with the spike glycoprotein RBD, TMPRSS2, CatB, and CatL were –5.02, –7.57, –10.10, and –6.11 kcal/mol, the values for GOS were –4.67, –5.24, –8.31, and –6.76, respectively. Furthermore, both compounds maintained their stability for at least 170 ns on respective targets in molecular dynamics simulations. The molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) calculations also corroborated these data. Considering Lipinski’s rule of five, ROB and GOS exhibited 3 (MW>500, N or O>10, NH or OH>5), and 1 (NH or OH>5) violations, respectively. Neither ROB nor GOS showed AMES toxicity or hepatotoxicity. The LD50 of these compounds in rats were 2.482 and 2.527 mol/kg, respectively. Therefore, we conclude that these compounds could be considered as alternative therapeutic agents in the treatment of COVID-19. However, the possible inhibitory effects of these compounds on cytochromes (CYPs) should be verified by in vitro or in vivo tests and their adverse effects on cellular energy metabolism should be minimized by performing molecular modifications if necessary.     

In silico analyses of major active constituents of fingerroot (Boesenbergia rotunda) unveils inhibitory activities against SARS-CoV-2 main protease enzyme

Boesenbergia rotunda (L.) Mansf., commonly known as fingerroot is a perennial herb in the Zingiberaceae family with anticancer, anti-leptospiral, anti-inflammatory, antioxidant, antiulcer, and anti-herpes viral activities. While the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) inhibitory activity of B. rotunda extract has been recently found, the active compounds contributing to this activity are yet unknown. The main protease (M^pro) enzyme is one of the most well established therapeutic targets among coronaviruses which plays a vital role in the maturation and cleavage of polyproteins during viral replication. The current work aims to identify active phytochemical substances from B. rotunda extract that can inhibit the replication of SARS-CoV-2 by using a combined molecular docking and dynamic simulation approaches. The virtual screening experiment revealed that fifteen molecules out of twenty-three major active compounds in the plant extract have acceptable drug-like characteristics. Alpinetin, Pinocembrin, and Pinostrobin have binding energies of ?7.51 kcal/mol, ?7.21 kcal/mol, and ?7.18 kcal/mol, respectively, and can suppress M^pro activity. The stability of the simulated complexes of the lead compounds with the drug-receptor is demonstrated by 100-ns MD simulations. The binding free energies study utilizing molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and molecular mechanics generalized Born surface area (MM-GBSA) show that the compounds and M^pro enzyme have favourable thermodynamic interactions, which are majorly driven by van der Waals forces. Thus, the selected bioactive phytochemicals from B. rotunda might be used as anti-SARS-CoV-2 candidates that target the M^pro enzyme.     

Potential targets of severe acute respiratory syndrome coronavirus 2 of clinical drug fluvoxamine: Docking and molecular dynamics studies to elucidate viral action

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continued evolving for survival and adaptation by mutating itself into different variants of concern, including omicron. Several studies and clinical trials found fluvoxamine, an Food and Drug Administration-approved antidepressant drug, to be effective at preventing mild coronavirus disease 2019 (COVID-19) from progressing to severe diseases. However, the mechanism of fluvoxamine's direct antiviral action against COVID-19 is still unknown. Fluvoxamine was docked with 11 SARS-CoV-2 targets and subjected to stability, conformational changes, and binding free energy analyses to explore its mode of action. Of the targets, nonstructural protein 14 (NSP14), main protease (Mpro), and papain-like protease (PLpro) had the best docking scores with fluvoxamine. Consistent with the docking results, it was confirmed by molecular dynamics simulations that the NSP14 N7-MTase ((N7-guanine)-methyltransferase)–fluvoxamine, Mpro–fluvoxamine, and PLpro–fluvoxamine complexes are stable, with the lowest binding free energies of −105.1, −82.7, and − 38.5 kJ/mol, respectively. A number of hotspot residues involved in the interaction were also identified. These include Glu166, Asp187, His41, and Cys145 in Mpro, Gly163 and Arg166 in PLpro, and Glu302, Gly333, and Phe426 in NSP14, which could aid in the development of better antivirals against SARS-CoV-2.     

Penta-peptide ATN-161 based neutralization mechanism of SARS-CoV-2 spike protein

SARS-CoV-2 has become a big challenge for the scientific community worldwide. SARS-CoV-2 enters into the host cell by the spike protein binding with an ACE2 receptor present on the host cell. Developing safe and effective inhibitor appears an urgent need to interrupt the binding of SARS-CoV-2 spike protein with ACE2 receptor in order to reduce the SARS-CoV-2 infection. We have examined the penta-peptide ATN-161 as potential inhibitor of ACE2 and SARS-CoV-2 spike protein binding, where ATN-161 has been commercially approved for the safety and possess high affinity and specificity towards the receptor binding domain (RBD) of S1 subunit in SARS-CoV-2 spike protein. We carried out experiments and confirmed these phenomena that the virus bindings were indeed minimized. ATN-161 peptide can be used as an inhibitor of protein-protein interaction (PPI) stands as a crucial interaction in biological systems. The molecular docking finding suggests that the binding energy of the ACE2-spike protein complex is reduced in the presence of ATN-161. Protein-protein docking binding energy (-40.50 kcal/mol) of the spike glycoprotein toward the human ACE2 and binding of ATN-161 at their binding interface reduced the biding energy (-26.25 kcal/mol). The finding of this study suggests that ATN-161 peptide can mask the RBD of the spike protein and be considered as a neutralizing candidate by binding with the ACE2 receptor. Peptide-based masking of spike S1 protein (RBD) and its neutralization is a highly promising strategy to prevent virus penetration into the host cell. Thus masking of the RBD leads to the loss of receptor recognition property which can reduce the chance of infection host cells.     

D614G substitution at the hinge region enhances the stability of trimeric SARS-CoV-2 spike protein

Mutations in the spike protein of SARS-CoV-2 are the major causes for the modulation of ongoing COVID-19 infection. Currently, the D614G substitution in the spike protein has become dominant worldwide. It is associated with higher infectivity than the ancestral (D614)variant. We demonstrate using Gaussian network model-based normal mode analysis that the D614G substitution occurs at the hinge region that facilitates domain-domain motions between receptor binding domain and S2 region of the spike protein. Computer-aided mutagenesis and inter-residue energy calculations reveal that contacts involving D614 are energetically frustrated. However, contacts involving G614 are energetically favourable, implying the substitution strengthens residue contacts that are formed within as well as between protomers. We also find that the free energy difference (ΔΔG) between two variants is -2.6 kcal/mol for closed and -2.0 kcal/mol for 1-RBD up conformation. Thus, the thermodynamic stability has increased upon D614G substitution. Whereas the reverse mutation in spike protein structures having G614 substitution has resulted in the free energy differences of 6.6 kcal/mol and 6.3 kcal/mol for closed and 1-RBD up conformations, respectively, indicating that the overall thermodynamic stability has decreased. These results suggest that the D614G substitution modulates the flexibility of spike protein and confers enhanced thermodynamic stability irrespective of conformational states. This data concurs with the known information demonstrating increased availability of the functional form of spikeprotein trimer upon D614G substitution.     

Transcriptome network analyses in human coronavirus infections suggest a rational use of immunomodulatory drugs for COVID-19 therapy

The recent outbreak of coronavirus disease 2019 (COVID-19) by SARS-CoV-2 has led to uptodate 24.3 M cases and 0.8 M deaths. It is thus in urgent need to rationalize potential therapeutic targets against the progression of diseases. An effective, feasible way is to use the pre-existing ΔORF6 mutant of SARS-CoV as a surrogate for SARS-CoV-2, since both lack the moiety responsible for interferon antagonistic effects. By analyzing temporal profiles of upregulated genes in ΔORF6-infected Calu-3 cells, we prioritized 55 genes and 238 ligands to reposition currently available medications for COVID-19 therapy. Eight of them are already in clinical trials, including dexamethasone, ritonavir, baricitinib, tofacitinib, naproxen, budesonide, ciclesonide and formoterol. We also pinpointed 16 drug groups from the Anatomical Therapeutic Chemical classification system, with the potential to mitigate symptoms of SARS-CoV-2 infection and thus to be repositioned for COVID-19 therapy.     

Analysis of the docking property of host variants of hACE2 for SARS-CoV-2 in a large cohort

The recent novel coronavirus disease (COVID-19) outbreak, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is threatening global health. However, an understanding of the interaction of SARS-CoV-2 with human cells, including the physical docking property influenced by the host’s genetic diversity, is still lacking. Here, based on germline variants in the UK Biobank covering 502,543 individuals, we revealed the molecular interactions between human angiotensin-converting enzyme 2 (hACE2), which is the representative receptor for SARS-CoV-2 entry, and COVID-19 infection. We identified six nonsense and missense variants of hACE2 from 2585 subjects in the UK Biobank covering 500000 individuals. Using our molecular dynamics simulations, three hACE2 variants from 2585 individuals we selected showed higher levels of binding free energy for docking in the range of 1.44–3.69 kcal/mol. Although there are diverse contributors to SARS-CoV-2 infections, including the mobility of individuals, we analyzed the diagnosis records of individuals with these three variants of hACE2. Our molecular dynamics simulations combined with population-based genomic data provided an atomistic understanding of the interaction between hACE2 and the spike protein of SARS-CoV-2.     

Caffeic acid derivatives (CAFDs) as inhibitors of SARS-CoV-2: CAFDs-based functional foods as a potential alternative approach to combat COVID-19

BackgroundSARS-CoV-2, an emerging strain of coronavirus, has affected millions of people from all the continents of world and received worldwide attention. This emerging health crisis calls for the urgent development of specific therapeutics against COVID-19 to potentially reduce the burden of this emerging pandemic.PurposeThis study aims to evaluate the anti-viral efficacy of natural bioactive entities against COVID-19 via molecular docking and molecular dynamics simulation.MethodsA library of 27 caffeic-acid derivatives was screened against 5 proteins of SARS-CoV-2 by using Molegro Virtual Docker 7 to obtain the binding energies and interactions between compounds and SARS-CoV-2 proteins. ADME properties and toxicity profiles were investigated via www.swissadme.ch web tools and Toxtree respectively. Molecular dynamics simulation was performed to determine the stability of the lead-protein interactions.ResultsOur obtained results has uncovered khainaoside C, 6-O-Caffeoylarbutin, khainaoside B, khainaoside C and vitexfolin A as potent modulators of COVID-19 possessing more binding energies than nelfinavir against COVID-19 M^pro, Nsp15, SARS-CoV-2 spike S2 subunit, spike open state and closed state structure respectively. While Calceolarioside B was identified as pan inhibitor, showing strong molecular interactions with all proteins except SARS-CoV-2 spike glycoprotein closed state. The results are supported by 20 ns molecular dynamics simulations of the best complexes.ConclusionThis study will hopefully pave a way for development of phytonutrients-based antiviral therapeutic for treatment or prevention of COVID-19 and further studies are recommended to evaluate the antiviral effects of these phytochemicals against SARS-CoV-2 in in vitro and in vivo models.     

Prolonged humoral and cellular immunity in COVID-19-recovered patients

By the beginning of 2021, the battle against coronavirus disease 2019 (COVID-19) remains ongoing. Investigating the adaptive immune response against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which causes COVID-19, in patients who have recovered from this disease could contribute to our understanding of the natural host immune response. We enrolled 38 participants in this study. 7 healthy participants and 31 COVID-19 patients who had recovered from COVID-19 and categorized them into 3 groups according to their previous clinical presentations: 10 moderate, 9 mild, and 12 asymptomatic. Flow cytometry analysis of peripheral lymphocyte counts in recovered patients showed significantly increased levels of CD4⁺ T cells in patients with a history of mild and moderate COVID-19 symptoms compared with those healthy individuals (p < 0.05 and p < 0.0001 respectively). whereas no significant difference was observed in the CD8⁺ T cell percentage in COVID-19-recovered patients compared with healthy individuals. Our study demonstrated that antibodies against the SARS-CoV-2 spike protein (anti-S) IgG antibody production could be observed in all recovered COVID-19 patients, regardless of whether they were asymptomatic (p < 0.05)or presented with mild (p < 0.0001) or moderate symptoms (p < 0.01). Anti-S IgG antibodies could be detected in participants up to 90 days post-infection. In conclusion, the lymphocyte levels in recovered patients were associated with the clinical presentation of the disease, and further analysis is required to investigate relationships between different clinical presentations and lymphocyte activation and function.     

Identification of SARS-CoV-2 RNA-dependent RNA polymerase inhibitors from the major phytochemicals of Nigella sativa: An in silico approach

The coronavirus disease 2019 (COVID-19), which emerged in December 2019, continues to be a serious health concern worldwide. There is an urgent need to develop effective drugs and vaccines to control the spread of this disease. In the current study, the main phytochemical compounds of Nigella sativa were screened for their binding affinity for the active site of the RNA-dependent RNA polymerase (RdRp) enzyme of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The binding affinity was investigated using molecular docking methods, and the interaction of phytochemicals with the RdRp active site was analyzed and visualized using suitable software. Out of the nine phytochemicals of N. sativa screened in this study, a significant docking score was observed for four compounds, namely α-hederin, dithymoquinone, nigellicine, and nigellidine. Based on the findings of our study, we report that α-hederin, which was found to possess the lowest binding energy (–8.6 kcal/mol) and hence the best binding affinity, is the best inhibitor of RdRp of SARS-CoV-2, among all the compounds screened here. Our results prove that the top four potential phytochemical molecules of N. sativa, especially α-hederin, could be considered for ongoing drug development strategies against SARS-CoV-2. However, further in vitro and in vivo testing are required to confirm the findings of this study.     

In silico evaluation of flavonoids as effective antiviral agents on the spike glycoprotein of SARS-CoV-2

The novel coronavirus pandemic has spread over in 213 countries as of July 2020. Approximately 12 million people have been infected so far according to the reports from World Health Organization (WHO). Preventive measures are being taken globally to avoid the rapid spread of virus. In the current study, an in silico approach is carried out as a means of inhibiting the spike protein of the novel coronavirus by flavonoids from natural sources that possess both antiviral and anti-inflammatory properties. The methodology is focused on molecular docking of 10 flavonoid compounds that are docked with the spike protein of SARS-CoV-2, to determine the highest binding affinity at the binding site. Molecular dynamics simulation was carried out with the flavonoid-protein complex showing the highest binding affinity and highest interactions. The flavonoid naringin showed the least binding energy of −9.8 Kcal/mol with the spike protein which was compared with the standard drug, dexamethasone which is being repurposed to treat critically ill patients. MD simulation was carried out on naringin-spike protein complex for their conformational stability in the active site of the novel coronavirus spike protein. The RMSD of the complex appeared to be more stable when compared to that of the protein from 0.2 nm to 0.4 nm. With the aid of this in silico approach further in vitro studies can be carried out on these flavonoids against the novel coronavirus as a means of viral protein inhibitors.     

Glycyrrhizic acid exerts inhibitory activity against the spike protein of SARS-CoV-2

Coronavirus causes a disease with high infectivity and pathogenicity, especially SARS in 2003, MERS in 2012, and COVID-2019 currently. The spike proteins of these coronaviruses are critical for host cell entry by receptors. Thus, searching for broad-spectrum anti-coronavirus candidates, such as spike protein inhibitors, is vital and desirable due to the mutations in the spike protein. In this study, a combination of computer-aided drug design and biological verification was used to discover active monomers from traditional Chinese medicine. Surface plasmon resonance (SPR) assays and NanoBit assays were used to verify the predicated compounds with their binding activities to spike proteins and inhibitory activities on the SARS-CoV-2 RBD/ACE2 interaction, respectively. Furthermore, an MTT assay was used to evaluate the cell toxicities of active compounds. As a result, glycyrrhizic acid (ZZY-44) was found to be the most efficient and nontoxic broad-spectrum anti-coronavirus molecule in vitro, especially, the significant effect on SARS-CoV-2, which provided a theoretical basis for the study of the pharmacodynamic material basis of traditional Chinese medicine against SARS-CoV-2 and offered a lead compound for further structural modification in order to obtain more effective candidate drugs against SARS-CoV-2.     

新型冠状病毒抗体和小分子抑制剂的研究现状

世界卫生组织已宣布新型冠状病毒感染（coronavirus disease 2019,COVID-19）的爆发为全球大流行。中和抗体和小分子抑制剂在预防及治疗COVID-19中发挥重要作用。尽管已开发出了多种中和抗体以及疫苗,但是随着病原体严重急性呼吸综合征冠状病毒2(severe acute respiratory syndrome coronavirus 2,SARS-CoV-2)的不断变异,现有的抗体及疫苗面临巨大的挑战。小分子抑制剂主要通过干扰病毒与宿主的结合以及病毒自身的复制达到消灭病毒以及抑制病毒感染的作用,并且对SARS-CoV-2突变株具有广谱抑制作用,是当前研究的热点。近年来国内外学者对SARS-CoV-2的小分子抑制剂做了大量的研究工作,本文根据中和抗体识别的抗原表位以及小分子抑制剂的作用机制分别对用于预防及治疗COVID-19的中和抗体和小分子抑制剂进行综述,讨论其研究现状,并展望小分子抑制剂的应用前景,以期为该领域的进一步研究提供参考。     

Multi-targeted molecular docking,pharmacokinetics, and drug-likeness evaluation of coumarin based compounds targeting proteins involved in development of COVID-19

COVID-19 is a progressing pandemic of coronavirus disease-2019, which had drowned the whole world in a deep sorrow sea. Uncountable deaths were extending the list of deaths every single day. The present research was aimed to study the multi-target interaction of coumarins against COVID-19 using molecular docking analysis. The structure of coumarin compounds was checked for ADME and Lipinski rule of five by using SwissADME, an online tool. SARS-CoV-2 proteins such as RdRp, PLpro, Mpro and spike protein were collected from the Protein Data Bank. The molecular docking study was performed in the PyRx tool, and the molecular interactions were visualised by Discovery Studio Visualizer. All the coumarin compounds used in the study were obeyed Lipinski’s rule of 5 without any violations. All the three designed derivatives of phenprocoumon, hymecromone, and psoralen were showed high binding affinity and prominent interactions with the drug target. The presence of –OH groups in the compound, His41, a catalytic dyad in Mpro, number of and the distance of hydrogen bond interactions with SARS-CoV-2 targets was accountable for the high binding attractions. The modified drug structures possess better binding efficacy towards at least three targets compared to their parent compounds. Further, molecular dynamic studies can be suggested to find the ligand–protein complex stability. The present study outcome reveals that the designed coumarins can be synthesised and examined as a potent inhibitory drug of SARS-CoV-2.     

Design,synthesis and biological evaluation of novel hydroxamic acids bearing artemisinin skeleton

A series of novel hydroxamic acids bearing artemisinin skeleton was designed and synthesized. Some compounds in this series exhibited moderate inhibition against the whole cell HDAC enzymes. Especially, compound 6g displayed potent cytotoxicity against three human cancer cell lines, including HepG2 (liver cancer), MCF-7 (breast cancer) and HL-60 (leukemia cancer), with IC₅₀ values of 2.50, 2.62 and 1.28 μg/mL, respectively. Docking studies performed with two potent compounds 6a and 6g using Autodock Vina showed that both compounds bound to HDAC2 with relatively high binding affinities from −7.1 to 7.0 kcal/mol compared to SAHA (−7.4 kcal/mol). It was found in this research that most of the target compounds seemed to be more cytotoxic toward blood cancer cells (HL-60) than liver (HepG2), and breast (MCF-7) cancer cells.     

A multiplex antigen microarray for simultaneous IgG and IgM detection against SARS-CoV-2 reveals higher seroprevalence than reported

The surge of SARS-CoV-2 has challenged health systems worldwide and efficient tests to detect viral particles, as well as antibodies generated against them, are needed. Specificity, sensitivity, promptness or scalability are the main parameters to estimate the final performance, but rarely all of them match in a single test. We have developed SCOVAM, a protein microarray with several viral antigens (spike, nucleocapsid, main protease Nsp5) as capturing probes in a fluorescence immunoassay for COVID-19 serological testing. SCOVAM depicts IgG and IgM antibody responses against each of these proteins of 22 individuals in a single microscope slide. It detects specific IgM (0.094 μg ml^-1) and IgG (~0.017 μg ml^-1) and is scalable and cost-effective. We validated SCOVAM by comparing with a widely used chemiluminescent commercial serological test (n = 742). SCOVAM showed twice the sensitivity and allowed following seroconversion in a single assay. By analysing the prevalence 4 months later in a subset of 76 positive sera, we still detected 93.42% of positives, almost doubling the detection of the commercial assay. The higher sensitivity of SCOVAM is especially relevant to screen sera for convalescent plasma-based treatments, high-throughput antibody response monitoring after vaccination or evaluation of vaccine efficiency.