Antiviral Efficacy of Selected Natural Phytochemicals against SARS-CoV-2 Spike Glycoprotein Using Structure-Based Drug Designing

SARS-CoV-2 is a highly virulent coronavirus that first surfaced in late 2019 and has since created a pandemic of the acute respiratory sickness known as “coronavirus disease 2019” (COVID-19), posing a threat to human health and public safety. S-RBD is a coronaviral protein that is essential for a coronavirus (CoV) to bind and penetrate into host cells. As a result, it has become a popular pharmacological target. The goal of this study was to find potential candidates for anti-coronavirus disease 2019 (COVID-19) drugs by targeting severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) S-RBD with novel bioactive compounds and molecular interaction studies of 15,000 phytochemicals belonging to different flavonoid subgroups. A spike protein crystal structure attached to the ACE2 structure was obtained from the PDB database. A library of 15,000 phytochemicals was made by collecting compounds from different databases, such as the Zinc-database, PubChem-database, and MPD3-database. This library was docked against a receptor binding domain of a spike glycoprotein through the Molecular Operating Environment (MOE). The top drug candidates Phylloflavan, Milk thistle, Ilexin B and Isosilybin B, after virtual screening, were selected on the basis of the least binding score. Phylloflavan ranked as the top compound because of its least binding affinity score of −14.09 kcal/mol. In silico studies showed that all those compounds showed good activity and could be used as an immunological response with no bioavailability issues. Absorption, distribution, metabolism, excretion and a toxicological analysis were conducted through SwissADME. Stability and effectiveness of the docked complexes were elucidated by performing the 100 ns molecular dynamic simulation through the Desmond package.     

Computational Screening of Phenylamino-Phenoxy-Quinoline Derivatives against the Main Protease of SARS-CoV-2 Using Molecular Docking and the ONIOM Method

In the search for new anti-HIV-1 agents, two forms of phenylamino-phenoxy-quinoline derivatives have been synthesized, namely, 2-phenylamino-4-phenoxy-quinoline and 6-phenylamino-4-phenoxy-quinoline. In this study, the binding interactions of phenylamino-phenoxy-quinoline derivatives and six commercially available drugs (hydroxychloroquine, ritonavir, remdesivir, S-217622, N3, and PF-07321332) with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (M^pro) were investigated using molecular docking and the ONIOM method. The molecular docking showed the hydrogen bonding and hydrophobic interactions of all the compounds in the pocket of SARS-CoV-2 main protease (M^pro), which plays an important role for the division and proliferation of the virus into the cell. The binding free energy values between the ligands and M^pro ranged from −7.06 to −10.61 kcal/mol. The molecular docking and ONIOM results suggested that 4-(2′,6′-dimethyl-4′-cyanophenoxy)-2-(4″-cyanophenyl)-aminoquinoline and 4-(4′-cyanophenoxy)-2-(4″-cyanophenyl)-aminoquinoline have low binding energy values and appropriate molecular properties; moreover, both compounds could bind to M^pro via hydrogen bonding and Pi-Pi stacking interactions with amino acid residues, namely, HIS41, GLU166, and GLN192. These amino acids are related to the proteolytic cleavage process of the catalytic triad mechanisms. Therefore, this study provides important information for further studies on synthetic quinoline derivatives as antiviral candidates in the treatment of SARS-CoV-2.     

Conducting the RBD of SARS-CoV-2 Omicron Variant with Phytoconstituents from Euphorbia dendroides to Repudiate the Binding of Spike Glycoprotein Using Computational Molecular Search and Simulation Approach

(1) Background: Natural constituents are still a preferred route for counteracting the outbreak of COVID-19. Essentially, flavonoids have been found to be among the most promising molecules identified as coronavirus inhibitors. Recently, a new SARS-CoV-2 B.1.1.529 variant has spread in many countries, which has raised awareness of the role of natural constituents in attempts to contribute to therapeutic protocols. (2) Methods: Using various chromatographic techniques, triterpenes (1–7), phenolics (8–11), and flavonoids (12–17) were isolated from Euphorbia dendroides and computationally screened against the receptor-binding domain (RBD) of the SARS-CoV-2 Omicron variant. As a first step, molecular docking calculations were performed for all investigated compounds. Promising compounds were subjected to molecular dynamics simulations (MD) for 200 ns, in addition to molecular mechanics Poisson–Boltzmann surface area calculations (MM/PBSA) to determine binding energy. (3) Results: MM/PBSA binding energy calculations showed that compound 14 (quercetin-3-O-β-D-glucuronopyranoside) and compound 15 (quercetin-3-O-glucuronide 6″-O-methyl ester) exhibited strong inhibition of Omicron, with ΔG_binding of −41.0 and −32.4 kcal/mol, respectively. Finally, drug likeness evaluations based on Lipinski’s rule of five also showed that the discovered compounds exhibited good oral bioavailability. (4) Conclusions: It is foreseeable that these results provide a novel intellectual contribution in light of the decreasing prevalence of SARS-CoV-2 B.1.1.529 and could be a good addition to the therapeutic protocol.     

Docking Analysis of Some Bioactive Compounds from Traditional Plants against SARS-CoV-2 Target Proteins

COVID-19 is still a global pandemic that has not been stopped. Many traditional medicines have been demonstrated to be incredibly helpful for treating COVID-19 patients while fighting the disease worldwide. We introduced 10 bioactive compounds derived from traditional medicinal plants and assessed their potential for inhibiting viral spike protein (S-protein), Papain-like protease (PLpro), and RNA dependent RNA polymerase (RdRp) using molecular docking protocols where we simulate the inhibitors bound to target proteins in various poses and at different known binding sites using Autodock version 4.0 and Chimera 1.8.1 software. Results found that the chicoric acid, quinine, and withaferin A ligand strongly inhibited CoV-2 S -protein with a binding energy of −8.63, −7.85, and −7.85 kcal/mol, respectively. Our modeling work also suggested that curcumin, quinine, and demothoxycurcumin exhibited high binding affinity toward RdRp with a binding energy of −7.80, −7.80, and −7.64 kcal/mol, respectively. The other ligands, namely chicoric acid, demothoxycurcumin, and curcumin express high binding energy than the other tested ligands docked to PLpro with −7.62, −6.81, and −6.70 kcal/mol, respectively. Prediction of drug-likeness properties revealed that all tested ligands have no violations to Lipinski’s Rule of Five except cepharanthine, chicoric acid, and theaflavin. Regarding the pharmacokinetic behavior, all ligand predicted to have high GI-absorption except chicoric acid and theaflavin. At the same way chicoric acid, withaferin A, and withanolide D predicted to be substrate for multidrug resistance protein (P-gp substrate). Caffeic acid, cepharanthine, chicoric acid, withaferin A, and withanolide D also have no inhibitory effect on any cytochrome P450 enzymes. Promisingly, chicoric acid, quinine, curcumin, and demothoxycurcumin exhibited high binding affinity on SARS-CoV-2 target proteins and expressed good drug-likeness and pharmacokinetic properties. Further research is required to investigate the potential uses of these compounds in the treatment of SARS-CoV-2.     

Virtual Screening and Molecular Dynamics Simulation Study of Influenza Polymerase PB2 Inhibitors

Influenza A virus is the main cause of worldwide epidemics and annual influenza outbreaks in humans. In this study, a virtual screen was performed to identify compounds that interact with the PB2 cap-binding domain (CBD) of influenza A polymerase. A virtual screening workflow based on Glide docking was used to screen an internal database containing 8417 molecules, and then the output compounds were selected based on solubility, absorbance, and structural fingerprints. Of the 16 compounds selected for biological evaluation, six compounds were identified that rescued cells from H1N1 virus-mediated death at non-cytotoxic concentrations, with EC50 values ranging from 2.5–55.43 μM, and that could bind to the PB2 CBD of H1N1, with Kd values ranging from 0.081–1.53 μM. Molecular dynamics (MD) simulations of the docking complexes of our active compounds revealed that each compound had its own binding characteristics that differed from those of VX-787. Our active compounds have novel structures and unique binding modes with PB2 proteins, and are suitable to serve as lead compounds for the development of PB2 inhibitors. An analysis of the MD simulation also helped us to identify the dominant amino acid residues that play a key role in binding the ligand to PB2, suggesting that we should focus on increasing and enhancing the interaction between inhibitors and these major amino acids during lead compound optimization to obtain more active PB2 inhibitors.     

Electron Density Analysis of SARS-CoV-2 RNA-Dependent RNA Polymerase Complexes

The work is devoted to the study of the complementarity of the electronic structures of the ligands and SARS-CoV-2 RNA-dependent RNA polymerase. The research methodology was based on determining of 3D maps of electron densities of complexes using an original quantum free-orbital AlteQ approach. We observed a positive relationship between the parameters of the electronic structure of the enzyme and ligands. A complementarity factor of the enzyme-ligand complexes has been proposed. The console applications of the AlteQ complementarity assessment for Windows and Linux (alteq_map_enzyme_ligand_4_win.exe and alteq_map_enzyme_ligand_4_linux) are available for free at the ChemoSophia webpage.     

VEGFR-2 与抑制剂Sunitinib 的分子对接及分子动力学研究

用分子对接方法研究了VEGFR-2 和抑制剂Sunitinib 的相互作用模式, 并对其复合物进行了10 ns 的分子动力学(Molecular Dynamics, MD)模拟. 结果表明, 抑制剂Sunitinib 能与VEGFR-2 中位于活性空腔的Glu885, Ile888, His1026,Asp1028, Asp1046 五个氨基酸残基形成疏水作用; 另外, VEGFR-2 中His1026, Cys1024, Asp1046 三个氨基酸残基能与Sunitinib 形成三个作用强度不同的氢键. 这些基团之间的相互作用是Sunitinib 抑制VEGFR-2 活性的关键因素. 研究结果可为VEGFR-2 抑制剂的结构改良、分子设计、合成提供理论参考, 并有助于寻找活性更高、效果更好的抗肿瘤药物.     

Computational identification of disulfiram and neratinib as putative SARS-CoV-2 main protease inhibitors

1. Download : Download high-res image (212KB)
2. Download : Download full-size image

     

In silico exploration of binding potentials of anti SARS-CoV-1 phytochemicals against main protease of SARS-CoV-2

The phytochemicals can play complementary medicine compared to synthetic drugs considering their natural origin, safety, and low cost. Phytochemicals hold a key position for the expansion of drug development against corona viruses and need better consideration to the agents that have already been shown to display effective activity against various strains of corona viruses. In this study, we performed molecular docking studies on potential forty seven phytochemicals which are SARS-CoV-1 M^pro inhibitors to identify potential candidate against the main proteins of SARS-CoV-2. In Silico Molecular docking studies revealed that phytochemicals 16 (Broussoflavan A), 22 (Dieckol), 31 (Hygromycin B), 45 (Sinigrin) and 46 (Theaflavin-3,3′-digallate) exhibited excellent SARS-CoV-2 M^pro inhibitors. Furthermore, supported by Molecular dynamics (MD) simulation analysis such as Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), Radius of gyration (Rg) and H-bond interaction analysis. We expect that our findings will provide designing principles for new corona virus strains and establish important frameworks for the future development of antiviral drugs.     

Insights into the Dynamics and Binding of Two Polyprotein Substrate Cleavage Points in the Context of the SARS-CoV-2 Main and Papain-like Proteases

It is well known that vital enzymes in the replication process of the coronavirus are the SARS-CoV-2 PLpro and SARS-CoV-2 3CLpro, both of which are important targets in the search for anti-coronavirus agents. These two enzymes are responsible for cleavage at various polyprotein sites in the SARS-CoV-2 lifecycle. Herein, the dynamics of the polyprotein cleavage sequences for the boundary between non-structural proteins Nsp1 and Nsp2 (CS1) and between Nsp2 and Nsp3 (CS2) in complex with both the papain-like protein PLpro and the main protease 3CLpro were explored using computational methods. The post dynamics analysis reveals that CS1 and CS2 both have greater stability when complexed with PLpro. Of these two, greater stability is observed for the CS1–PLpro complex, while destabilization resulting in loss of CS2 from the PLpro active site is observed for CS2-PLpro, suggesting the rate of exchange by the papain-like protease is faster for CS2 compared to CS1. On the other hand, the 3CLpro main protease also reveals stability for CS1 suggesting that the main protease could also play a potential role in the cleavage at point CS1. However, destabilization occurs early in the simulation for the complex CLpro–CS2 suggesting a poor interaction and non-plausible protease cleavage of the polyprotein at CS2 by the main protease. These findings could be used as a guide in the development and design of potent COVID-19 antiviral inhibitors that mimic the CS1 cleavage site.     

In Silico Identification and Validation of Organic Triazole Based Ligands as Potential Inhibitory Drug Compounds of SARS-CoV-2 Main Protease

The SARS-CoV-2 virus is highly contagious to humans and has caused a pandemic of global proportions. Despite worldwide research efforts, efficient targeted therapies against the virus are still lacking. With the ready availability of the macromolecular structures of coronavirus and its known variants, the search for anti-SARS-CoV-2 therapeutics through in silico analysis has become a highly promising field of research. In this study, we investigate the inhibiting potentialities of triazole-based compounds against the SARS-CoV-2 main protease (M^pro). The SARS-CoV-2 main protease (M^pro) is known to play a prominent role in the processing of polyproteins that are translated from the viral RNA. Compounds were pre-screened from 171 candidates (collected from the DrugBank database). The results showed that four candidates (Bemcentinib, Bisoctrizole, PYIITM, and NIPFC) had high binding affinity values and had the potential to interrupt the main protease (M^pro) activities of the SARS-CoV-2 virus. The pharmacokinetic parameters of these candidates were assessed and through molecular dynamic (MD) simulation their stability, interaction, and conformation were analyzed. In summary, this study identified the most suitable compounds for targeting Mpro, and we recommend using these compounds as potential drug molecules against SARS-CoV-2 after follow up studies.     

Aphrodisiac Performance of Bioactive Compounds from Mimosa pudica Linn.: In Silico Molecular Docking and Dynamics Simulation Approach

Plants and their derived molecules have been traditionally used to manage numerous pathological complications, including male erectile dysfunction (ED). Mimosa pudica Linn. commonly referred to as the touch-me-not plant, and its extract are important sources of new lead molecules in drug discovery research. The main goal of this study was to predict highly effective molecules from M. pudica Linn. for reaching and maintaining penile erection before and during sexual intercourse through in silico molecular docking and dynamics simulation tools. A total of 28 bioactive molecules were identified from this target plant through public repositories, and their chemical structures were drawn using Chemsketch software. Graph theoretical network principles were applied to identify the ideal target (phosphodiesterase type 5) and rebuild the network to visualize the responsible signaling genes, proteins, and enzymes. The 28 identified bioactive molecules were docked against the phosphodiesterase type 5 (PDE5) enzyme and compared with the standard PDE5 inhibitor (sildenafil). Pharmacokinetics (ADME), toxicity, and several physicochemical properties of bioactive molecules were assessed to confirm their drug-likeness property. Molecular dynamics (MD) simulation modeling was performed to investigate the stability of PDE5–ligand complexes. Four bioactive molecules (Bufadienolide (−12.30 kcal mol⁻¹), Stigmasterol (−11.40 kcal mol⁻¹), Isovitexin (−11.20 kcal mol⁻¹), and Apigetrin (−11.20 kcal mol⁻¹)) showed the top binding affinities with the PDE5 enzyme, much more powerful than the standard PDE5 inhibitor (−9.80 kcal mol⁻¹). The four top binding bioactive molecules were further validated for a stable binding affinity with the PDE5 enzyme and conformation during the MD simulation period as compared to the apoprotein and standard PDE5 inhibitor complexes. Further, the four top binding bioactive molecules demonstrated significant drug-likeness characteristics with lower toxicity profiles. According to the findings, the four top binding molecules may be used as potent and safe PDE5 inhibitors and could potentially be used in the treatment of ED.     

Marine Brominated Tyrosine Alkaloids as Promising Inhibitors of SARS-CoV-2

There have been more than 150 million confirmed cases of SARS-CoV-2 since the beginning of the pandemic in 2019. By June 2021, the mortality from such infections approached 3.9 million people. Despite the availability of a number of vaccines which provide protection against this virus, the evolution of new viral variants, inconsistent availability of the vaccine around the world, and vaccine hesitancy, in some countries, makes it unreasonable to rely on mass vaccination alone to combat this pandemic. Consequently, much effort is directed to identifying potential antiviral treatments. Marine brominated tyrosine alkaloids are recognized to have antiviral potential. We test here the antiviral capacity of fourteen marine brominated tyrosine alkaloids against five different target proteins from SARS-CoV-2, including main protease (M^pro) (PDB ID: 6lu7), spike glycoprotein (PDB ID: 6VYB), nucleocapsid phosphoprotein (PDB ID: 6VYO), membrane glycoprotein (PDB ID: 6M17), and non-structural protein 10 (nsp10) (PDB ID: 6W4H). These marine alkaloids, particularly the hexabrominated compound, fistularin-3, shows promising docking interactions with predicted binding affinities (S-score = −7.78, −7.65, −6.39, −6.28, −8.84 Kcal/mol) for the main protease (M^pro) (PDB ID: 6lu7), spike glycoprotein (PDB ID: 6VYB), nucleocapsid phosphoprotein (PDB ID: 6VYO), membrane glycoprotein (PDB ID: 6M17), and non-structural protein 10 (nsp10) (PDB ID: 6W4H), respectively, where it forms better interactions with the protein pockets than the native interaction. It also shows promising molecular dynamics, pharmacokinetics, and toxicity profiles. As such, further exploration of the antiviral properties of fistularin-3 against SARS-CoV-2 is merited.     

Naturally Available Flavonoid Aglycones as Potential Antiviral Drug Candidates against SARS-CoV-2

Flavonoids are important secondary plant metabolites that have been studied for a long time for their therapeutic potential in inflammatory diseases because of their cytokine-modulatory effects. Five flavonoid aglycones were isolated and identified from the hydrolyzed aqueous methanol extracts of Anastatica hierochuntica L., Citrus reticulata Blanco, and Kickxia aegyptiaca (L.) Nabelek. They were identified as taxifolin (1), pectolinarigenin (2), tangeretin (3), gardenin B (4), and hispidulin (5). These structures were elucidated based on chromatographic and spectral analysis. In this study, molecular docking studies were carried out for the isolated and identified compounds against SARS-CoV-2 main protease (Mpro) compared to the co-crystallized inhibitor of SARS-CoV-2 Mpro (α-ketoamide inhibitor (KI), IC₅₀ = 66.72 µg/mL) as a reference standard. Moreover, in vitro screening against SARS-CoV-2 was evaluated. Compounds 2 and 3 showed the highest virus inhibition with IC₅₀ 12.4 and 2.5 µg/mL, respectively. Our findings recommend further advanced in vitro and in vivo studies of the examined isolated flavonoids, especially pectolinarigenin (2), tangeretin (3), and gardenin B (4), either alone or in combination with each other to identify a promising lead to target SARS-CoV-2 effectively. This is the first report of the activity of these compounds against SARS-CoV-2.     

Identification of natural compounds targeting SARS-CoV-2 Mpro by virtual screening and molecular dynamics simulations

The SARS-CoV-2 main protease (Mpro) has been chosen as a conserved molecular target to develop broad-spectrum antiviral drugs. Using molecular docking and molecular dynamics (MD) simulations, a total of 5600 natural compounds available for virtual screening were tested to identify potential inhibitors of SARS-CoV-2 Mpro. As a result, three natural compounds (pentagalloylglucose, malonylawobanin and gnetin E dihydride) were found to be potential inhibitors of SARS-CoV-2, which confirms the theoretical and practical significance of this approach for the design of SARS-CoV-2 inhibitors.     

ZnCl2熔盐的分子动力学模拟

近年来，若干作者根据几种粒子间函数，对ZnCl2熔盐结构做过分子动力学模拟［1－3］其出的Zn－Cl和Cl－Cl离子间的偏径向分布函数与中子衍射实测值符合较好，但Zn－Zn离子间距以及Zn－Zn间配位数计算值多偏高．鉴于Busing势函数在多价卤化物馆盐的分子动力学计算中应用效果较好[4],我们试用Busing势函数为基础对ZnCl2熔盐结构和能量做分子动力学计算．1研究方法计算所用粒子势为Busing势函数此处，Zi为离子的电行数（ZZn。＋＝2，Zcl－＝1），几；是离子有效半径，人为＊离子的“硬度”参数·据文献问，f二0．005071，尸zn。十二0．…     

Computational Approaches to Discover Novel Natural Compounds for SARS-CoV-2 Therapeutics

Scientists all over the world are facing a challenging task of finding effective therapeutics for the coronavirus disease (COVID-19). One of the fastest ways of finding putative drug candidates is the use of computational drug discovery approaches. The purpose of the current study is to retrieve natural compounds that have obeyed to drug-like properties as potential inhibitors. Computational molecular modelling techniques were employed to discover compounds with potential SARS-CoV-2 inhibition properties. Accordingly, the InterBioScreen (IBS) database was obtained and was prepared by minimizing the compounds. To the resultant compounds, the absorption, distribution, metabolism, excretion and toxicity (ADMET) and Lipinski's Rule of Five was applied to yield drug-like compounds. The obtained compounds were subjected to molecular dynamics simulation studies to evaluate their stabilities. In the current article, we have employed the docking based virtual screening method using InterBioScreen (IBS) natural compound database yielding two compounds has potential hits. These compounds have demonstrated higher binding affinity scores than the reference compound together with good pharmacokinetic properties. Additionally, the identified hits have displayed stable interaction results inferred by molecular dynamics simulation results. Taken together, we advocate the use of two natural compounds, STOCK1N-71493 and STOCK1N-45683 as SARS-CoV-2 treatment regime.     

DCN1-UBE2M相互作用抑制剂的3D-QSAR、分子对接和分子动力学研究

类泛素化是一种蛋白质翻译后修饰,其异常会导致神经退行性疾病和多种肿瘤的发生,因此它被视为有希望的抗肿瘤靶标.研究表明,抑制DCN1-UBE2M相互作用可选择性阻遏类泛素化.本文基于哌啶基脲类DCN1-UBE2M相互作用抑制剂进行3D-QSAR、分子对接和分子动力学模拟研究.利用3D-QSAR中的CoMFA和CoMSIA...     

诺氟沙星-DNA复合物的分子动力学模拟

采用分子模建的方法构建了诺氟沙星-DNA复合物的初始结构, 通过2 ns的分子动力学(MD)模拟研究表明: 诺氟沙星能够和双螺旋d[ATATCGATAT]₂形成稳定的复合物, 药物分子可紧密结合在DNA的小沟区域, 并且能够与DNA的鸟嘌呤碱基形成两个稳定的氢键. 在分子水平上提供了诺氟沙星直接与双螺旋DNA相互作用的结构及复合物的动态变化情况.