Identification of Potential SARS-CoV-2 Main Protease and Spike Protein Inhibitors from the Genus Aloe: An In Silico Study for Drug Development

Severe acute respiratory syndrome coronavirus (SARS-CoV-2) disease is a global rapidly spreading virus showing very high rates of complications and mortality. Till now, there is no effective specific treatment for the disease. Aloe is a rich source of isolated phytoconstituents that have an enormous range of biological activities. Since there are no available experimental techniques to examine these compounds for antiviral activity against SARS-CoV-2, we employed an in silico approach involving molecular docking, dynamics simulation, and binding free energy calculation using SARS-CoV-2 essential proteins as main protease and spike protein to identify lead compounds from Aloe that may help in novel drug discovery. Results retrieved from docking and molecular dynamics simulation suggested a number of promising inhibitors from Aloe. Root mean square deviation (RMSD) and root mean square fluctuation (RMSF) calculations indicated that compounds 132, 134, and 159 were the best scoring compounds against main protease, while compounds 115, 120, and 131 were the best scoring ones against spike glycoprotein. Compounds 120 and 131 were able to achieve significant stability and binding free energies during molecular dynamics simulation. In addition, the highest scoring compounds were investigated for their pharmacokinetic properties and drug-likeness. The Aloe compounds are promising active phytoconstituents for drug development for SARS-CoV-2.     

GC-MS,LC-MS/MS,Docking and Molecular Dynamics Approaches to Identify Potential SARS-CoV-2 3-Chymotrypsin-Like Protease Inhibitors from Zingiber officinale Roscoe

This study aims to identify and isolate the secondary metabolites of Zingiber officinale using GC-MS, preparative TLC, and LC-MS/MS methods, to evaluate the inhibitory potency on SARS-CoV-2 3 chymotrypsin-like protease enzyme, as well as to study the molecular interaction and stability by using docking and molecular dynamics simulations. GC-MS analysis suggested for the isolation of terpenoids compounds as major compounds on methanol extract of pseudostems and rhizomes. Isolation and LC-MS/MS analysis identified 5-hydro-7, 8, 2′-trimethoxyflavanone (9), (E)-hexadecyl-ferulate (1), isocyperol (2), N-isobutyl-(2E,4E)-octadecadienamide (3), and nootkatone (4) from the rhizome extract, as well as from the leaves extract with the absence of 9. Three known steroid compounds, i.e., spinasterone (7), spinasterol (8), and 24-methylcholesta-7-en-3β-on (6), were further identified from the pseudostem extract. Molecular docking showed that steroids compounds 7, 8, and 6 have lower predictive binding energies (MMGBSA) than other metabolites with binding energy of −87.91, −78.11, and −68.80 kcal/mole, respectively. Further characterization on the single isolated compound by NMR showed that 6 was identified and possessed 75% inhibitory activity on SARS-CoV-2 3CL protease enzyme that was slightly different with the positive control GC376 (77%). MD simulations showed the complex stability with compound 6 during 100 ns simulation time.     

Computational Simulation of HIV Protease Inhibitors to the Main Protease (Mpro) of SARS-CoV-2: Implications for COVID-19 Drugs Design

SARS-CoV-2 is highly homologous to SARS-CoV. To date, the main protease (Mpro) of SARS-CoV-2 is regarded as an important drug target for the treatment of Coronavirus Disease 2019 (COVID-19). Some experiments confirmed that several HIV protease inhibitors present the inhibitory effects on the replication of SARS-CoV-2 by inhibiting Mpro. However, the mechanism of action has still not been studied very clearly. In this work, the interaction mechanism of four HIV protease inhibitors Darunavir (DRV), Lopinavir (LPV), Nelfinavir (NFV), and Ritonavire (RTV) targeting SARS-CoV-2 Mpro was explored by applying docking, molecular dynamics (MD) simulations, and MM–GBSA methods using the broad-spectrum antiviral drug Ribavirin (RBV) as the negative and nonspecific control. Our results revealed that LPV, RTV, and NFV have higher binding affinities with Mpro, and they all interact with catalytic residues His41 and the other two key amino acids Met49 and Met165. Pharmacophore model analysis further revealed that the aromatic ring, hydrogen bond donor, and hydrophobic group are the essential infrastructure of Mpro inhibitors. Overall, this study applied computational simulation methods to study the interaction mechanism of HIV-1 protease inhibitors with SARS-CoV-2 Mpro, and the findings provide useful insights for the development of novel anti-SARS-CoV-2 agents for the treatment of COVID-19.     

A molecular docking study of potential inhibitors and repurposed drugs against SARS-CoV-2 main protease enzyme

COVID-19 has affected millions of people. Although many drugs are in use to combat disease, there is not any sufficient treatment yet. Having critical role in propagation of the novel coronavirus (SARS-CoV-2) works Main Protease up into a significant drug target. We have performed a molecular docking study to define possible inhibitor candidates against SARS-CoV-2 Main Protease enzyme. Besides docking Remdesivir, Ribavirin, Chloroquine and 28 other antiviral inhibitors (totally 31 inhibitors) to Main Protease enzyme, we have also performed a molecular docking study of 2177 ligands, which are used against Main Protease for the first time by using molecular docking program Autodock4. All ligands were successfully docked into Main Protease enzyme binding site. Among all ligands, EY16 coded ligand which previously used as EBNA1-DNA binding blocker candidate showed the best score for Main Protease with a binding free energy of −10.83 ​kcal/mol which was also lower than re-docking score of N3 ligand (−10.72 ​kcal/mol) contained in crystal structure of Main Protease. After analyzing the docking modes and docking scores we have found that our ligands have better binding free energy values than the inhibitors in use of treatment. We believe that further studies such as molecular dynamics or Molecular Mechanic Poisson Boltzmann Surface Area studies can make contribution that is more exhaustive to the docking results.     

Comparison of clinically approved molecules on SARS-CoV-2 drug target proteins: a molecular docking study

The new type of coronavirus, SARS-CoV-2 has affected more than 22.6 million people worldwide. Since the first day the virus was spotted in Wuhan, China, numerous drug design studies have been conducted all over the globe. Most of these studies target the receptor-binding domain of spike protein of SARS-CoV-2, which is known to bind to the human ACE2 receptor and SARS-CoV-2 main protease, vital for the virus’ replication. However, there might be a third target, human furin protease, which cleaves the virus’ S1-S2 domains playing an active role in its entry into the host cell. In this study, we docked five clinically used drug molecules, favipiravir, hydroxychloroquine, remdesivir, lopinavir, and ritonavir onto three target proteins, the receptor-binding domain of SARS-CoV-2 spike protein, SARS-CoV-2 main protease, and human furin protease. Results of molecular docking simulations revealed that human furin protease might be targeted by COVID-19. Remdesivir, a nucleic acid derivative, strongly bound to the active site of this protease, suggesting that this molecule can be used as a template for designing novel furin protease inhibitors to fight against the disease. Protein-drug interactions revealed in this study at the molecular level, can pave the way for better drug design for each specific target.     

Antiviral Activity of Metabolites from Peruvian Plants against SARS-CoV-2: An In Silico Approach

(1) Background: The COVID-19 pandemic lacks treatments; for this reason, the search for potential compounds against therapeutic targets is still necessary. Bioinformatics tools have allowed the rapid in silico screening of possible new metabolite candidates from natural resources or repurposing known ones. Thus, in this work, we aimed to select phytochemical candidates from Peruvian plants with antiviral potential against three therapeutical targets of SARS-CoV-2. (2) Methods: We applied in silico technics, such as virtual screening, molecular docking, molecular dynamics simulation, and MM/GBSA estimation. (3) Results: Rutin, a compound present in Peruvian native plants, showed affinity against three targets of SARS-CoV-2. The molecular dynamics simulation demonstrated the high stability of receptor–ligand systems during the time of the simulation. Our results showed that the Mpro-Rutin system exhibited higher binding free energy than PLpro-Rutin and N-Rutin systems through MM/GBSA analysis. (4) Conclusions: Our study provides insight on natural metabolites from Peruvian plants with therapeutical potential. We found Rutin as a potential candidate with multiple pharmacological properties against SARS-CoV-2.     

In silico studies on phytochemicals to combat the emerging COVID-19 infection

The current COVID-19 pandemic, caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and its variants, remains a serious health hazard globally. The SARS-CoV-2 Mpro and spike proteins, as well as the human ACE2 receptor, have previously been reported as good targets for the development of new drug leads to combat COVID-19. Various ligands, including synthetic and plant-derived small molecules, can interact with the aforementioned proteins. In this study, we investigated the interaction of eight phytochemicals, from selected medicinal plants (Aegle marmelos, Azadirachta indica, and Ocimum sanctum) commonly used in Indian traditional medicine, with SARS-CoV-2 Mpro (PDBID: 6LU7), SARS-CoV-2S spike protein (PDB ID: 6M0J) and the human ACE2 receptor (PDB ID: 6M18). All compounds were subjected to density functional theory (DFT) and frontier molecular orbitals (FMO) analysis to determine their geometry, and key electronic and energetic properties. Upon examining the interactions of the phytochemicals with the human ACE2 receptor and the SARS-CoV-2 Mpro, spike protein targets, two compounds (C-5 and C-8) were identified as the best binding ligands. These were further examined in MD simulation studies to determine the stability of the ligand–protein interactions. QSAR, pharmacokinetic and drug-likeness properties studies revealed that C-5 may be the best candidate to serve as a template for the design and development of new drugs to combat COVID-19.     

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.     

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.     

Prediction of COVID-19 manipulation by selective ACE inhibitory compounds of Potentilla reptant root: In silico study and ADMET profile

In the novel SARS-CoV-2 (COVID-19) as a global emergency event, the main reason of the cardiac injury from COVID-19 is angiotensin-converting enzyme 2 (ACE2) targeting in SARS-CoV-2 infection. The inhibition of ACE2 induces an increase in the angiotensin II (Ang II) and the angiotensin II receptor type 1 (AT1R) leading to impaired cardiac function or cardiac inflammatory responses. The ethyl acetate fraction of Potentilla reptans L. root can rescue heart dysfunction, oxidative stress, cardiac arrhythmias and apoptosis. Therefore, isolated components of P. reptans evaluated to identify natural anti-SARS-CoV-2 agents via molecular docking.In silico molecular docking study were carried out using the Auto Dock software on the isolated compounds of Potentilla reptans root. The protein targets of selective ACE and others obtained from Protein Data Bank (PDB). The best binding pose between amino acid residues involved in active site of the targets and compounds was discovered via molecular docking. Furthermore, ADMET properties of the compounds were evaluated.The triterpenoids of P. reptans showed more ACE inhibitory potential than catechin in both domains. They were selective on the nACE domain, especially compound 5. Also, the compound 5 & 6 had the highest binding affinity toward active site of nACE, cACE, AT1R, ACE2, and TNF-α receptors. Meanwhile, compound 3 showed more activity to inhibit TXA2. Drug likeness and ADMET analysis showed that the compounds passed the criteria of drug likeness and Lipinski rules. The current study depicted that P. reptans root showed cardioprotective effect in COVID-19 infection and manipulation of angiotensin II-induced side effects.     

Interaction Between Tangeretin and Ovalbumin to Reduce the Allergic Effects of Ovalbumin

Food allergy has become a worldwide problem. To find an effective way to reduce the allergic effect of allergen is an interesting topic. In this work, the interaction between ovalbumin(OVA) and tangeretin was studied by means of spectroscopy, molecular docking calculation and animal experiment. As the results show, static quenching of fluorescence intensity happened due to the interaction between OVA and tangeretin. Binding constants, binding sites, and thermometric parameters were obtained. The conformation changing of OVA occurred due to the interaction between it and tangeretin. The molecular docking calculation results show that the epitopes Arginine 84(ARG84), Leucine 87(LEU87), Asparagine 88(ASN88), Serine 103(SER103), Arginine 104(ARG104), Leucine 124(LEU124), Arginine 126(ARG126), Glycine 127(GLY127), Glycine 128(GLY128) and Tryptohan 148(TRP148) of OVA were occupied by tangeretin after the docking, which means tangeretin may have some inhibition effects. As the results of animal experiments, the amount of IgE and lung tissue sections of female BALB/c mice in different groups show that tangeretin with different concentrations has better effects than dexamethasone. In conclusion, the tangeretin may have the effects to reduce the allergic effects of OVA via occupying the IgE-binding epitopes of OVA.     

Targeting 3CLpro and SARS-CoV-2 RdRp by Amphimedon sp. Metabolites: A Computational Study

Since December 2019, novel coronavirus disease 2019 (COVID-19) pandemic has caused tremendous economic loss and serious health problems worldwide. In this study, we investigated 14 natural compounds isolated from Amphimedon sp. via a molecular docking study, to examine their ability to act as anti-COVID-19 agents. Moreover, the pharmacokinetic properties of the most promising compounds were studied. The docking study showed that virtually screened compounds were effective against the new coronavirus via dual inhibition of SARS-CoV-2 RdRp and the 3CL main protease. In particular, nakinadine B (1), 20-hepacosenoic acid (11) and amphimedoside C (12) were the most promising compounds, as they demonstrated good interactions with the pockets of both enzymes. Based on the analysis of the molecular docking results, compounds 1 and 12 were selected for molecular dynamics simulation studies. Our results showed Amphimedon sp. to be a rich source for anti-COVID-19 metabolites.     

In Silico Mining of Terpenes from Red-Sea Invertebrates for SARS-CoV-2 Main Protease (Mpro) Inhibitors

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent for the COVID-19 pandemic, which generated more than 1.82 million deaths in 2020 alone, in addition to 83.8 million infections. Currently, there is no antiviral medication to treat COVID-19. In the search for drug leads, marine-derived metabolites are reported here as prospective SARS-CoV-2 inhibitors. Two hundred and twenty-seven terpene natural products isolated from the biodiverse Red-Sea ecosystem were screened for inhibitor activity against the SARS-CoV-2 main protease (M^pro) using molecular docking and molecular dynamics (MD) simulations combined with molecular mechanics/generalized Born surface area binding energy calculations. On the basis of in silico analyses, six terpenes demonstrated high potency as M^pro inhibitors with ΔG_binding ≤ −40.0 kcal/mol. The stability and binding affinity of the most potent metabolite, erylosides B, were compared to the human immunodeficiency virus protease inhibitor, lopinavir. Erylosides B showed greater binding affinity towards SARS-CoV-2 M^pro than lopinavir over 100 ns with ΔG_binding values of −51.9 vs. −33.6 kcal/mol, respectively. Protein–protein interactions indicate that erylosides B biochemical signaling shares gene components that mediate severe acute respiratory syndrome diseases, including the cytokine- and immune-signaling components BCL2L1, IL2, and PRKC. Pathway enrichment analysis and Boolean network modeling were performed towards a deep dissection and mining of the erylosides B target–function interactions. The current study identifies erylosides B as a promising anti-COVID-19 drug lead that warrants further in vitro and in vivo testing.     

In Silico Evaluation of Iranian Medicinal Plant Phytoconstituents as Inhibitors against Main Protease and the Receptor-Binding Domain of SARS-CoV-2

The novel coronavirus disease 2019 (COVID-19) is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which initially appeared in Wuhan, China, in December 2019. Elderly individuals and those with comorbid conditions may be more vulnerable to this disease. Consequently, several research laboratories continue to focus on developing drugs to treat this infection because this disease has developed into a global pandemic with an extremely limited number of specific treatments available. Natural herbal remedies have long been used to treat illnesses in a variety of cultures. Modern medicine has achieved success due to the effectiveness of traditional medicines, which are derived from medicinal plants. The objective of this study was to determine whether components of natural origin from Iranian medicinal plants have an antiviral effect that can prevent humans from this coronavirus infection using the most reliable molecular docking method; in our case, we focused on the main protease (M^pro) and a receptor-binding domain (RBD). The results of molecular docking showed that among 169 molecules of natural origin from common Iranian medicinal plants, 20 molecules (chelidimerine, rutin, fumariline, catechin gallate, adlumidine, astragalin, somniferine, etc.) can be proposed as inhibitors against this coronavirus based on the binding free energy and type of interactions between these molecules and the studied proteins. Moreover, a molecular dynamics simulation study revealed that the chelidimerine–M^pro and somniferine–RBD complexes were stable for up to 50 ns below 0.5 nm. Our results provide valuable insights into this mechanism, which sheds light on future structure-based designs of high-potency inhibitors for SARS-CoV-2.     

Ligand release mechanisms and channels in histone deacetylases

Exploring the molecular channels of class I histone deacetylases (HDACs) with buried active sites are important to understand their structures and functionalities. In this work, we perform hybrid classical molecular dynamics and random acceleration molecular dynamics simulations to explore the B3N [i.e., (4‐(dimethylamino)N‐[7(hydroxyamino)‐7‐oxoheptyle] benzamide)] exit channels in the x‐ray crystal structures of HDAC3 and HDAC8 enzymes. Our simulations identify B3N release through four different channels in HDAC3 (denoted as A1, A2, B1, and B2) and HDAC8 (referred as A1, B1, B2, and B3) enzymes, among which egression through channel A1 is more predominant in both the enzymes. This mechanism is similar to ligand release in HDAC1 and HDAC2 described in our previous study and can be the fingerprint ligand release mechanisms in class I HDACs. Ligand release events through B channels, on the other hand, are different among HDAC3 and HDAC8, highlighting the significances of substituted residues in controlling the access to these channels This study reveals a novel aromatic gating mechanism elicited by TYR154‐TRP141‐TYR111 that controls the B3N access to all the B channels in HDAC8. The TRP141 in HDAC8 is substituted by LEU133 in HDAC3, which do not hinder the access to B channels in HDAC3. However, two hydrogen bonded barricades formed as ARG28‐GLY297‐GLY295‐GLY131 and TRP129‐ARG28‐ALA130‐LEU29‐TRP129 obstruct the B3N from exploring the B channels in HDAC3. The structural and dynamical characterizations of molecular channels and ligand unbinding mechanisms reported in this study provide novel structural insights and atomic level perspectives on HDAC3 and HDAC8 enzymes, thereby potentially aiding in the design of more specific HDAC inhibitors.Copyright © 2013 Wiley Periodicals, Inc.     

分子对接和动力学模拟提高嗜热蛋白酶PhpI的活力

通过分子对接和动力学模拟对嗜热蛋白酶的分子进行改造, 确定蛋白酶PH1704(PhpI)定点突变残基, 并通过分子生物学实验进行验证. 突变体K43C的蛋白酶活力提高了5.8倍. 分子动力学模拟结果表明, 经过8 ns的动力学模拟后, K43C突变体二级结构由野生型的S2片层(F11-E12-D13)变成环状结构. E12和K43均是活性位点的重要残基, 这种变化将导致活性位点的柔性增强, 有利于催化反应的发生.     

Molecular Docking and Dynamics Investigations for Identifying Potential Inhibitors of the 3-Chymotrypsin-like Protease of SARS-CoV-2: Repurposing of Approved Pyrimidonic Pharmaceuticals for COVID-19 Treatment

This study demonstrates the inhibitory effect of 42 pyrimidonic pharmaceuticals (PPs) on the 3-chymotrypsin-like protease of SARS-CoV-2 (3CL^pro) through molecular docking, molecular dynamics simulations, and free binding energies by means of molecular mechanics–Poisson Boltzmann surface area (MM-PBSA) and molecular mechanics–generalized Born surface area (MM-GBSA). Of these tested PPs, 11 drugs approved by the US Food and Drug Administration showed an excellent binding affinity to the catalytic residues of 3CL^pro of His41 and Cys145: uracil mustard, cytarabine, floxuridine, trifluridine, stavudine, lamivudine, zalcitabine, telbivudine, tipiracil, citicoline, and uridine triacetate. Their percentage of residues involved in binding at the active sites ranged from 56 to 100, and their binding affinities were in the range from −4.6 ± 0.14 to −7.0 ± 0.19 kcal/mol. The molecular dynamics as determined by a 200 ns simulation run of solvated docked complexes confirmed the stability of PP conformations that bound to the catalytic dyad and the active sites of 3CL^pro. The free energy of binding also demonstrates the stability of the PP–3CL^pro complexes. Citicoline and uridine triacetate showed free binding energies of −25.53 and −7.07 kcal/mol, respectively. Therefore, I recommend that they be repurposed for the fight against COVID-19, following proper experimental and clinical validation.     

Site mapping and small molecule blind docking reveal a possible target site on the SARS-CoV-2 main protease dimer interface

The SARS-CoV-2 virus is causing COVID-19 resulting in an ongoing pandemic with serious health, social, and economic implications. Much research is focused in repurposing or identifying new small molecules which may interact with viral or host-cell molecular targets. An important SARS-CoV-2 target is the main protease (M^pro), and the peptidomimetic α-ketoamides represent prototypical experimental inhibitors. The protease is characterised by the dimerization of two monomers each which contains the catalytic dyad defined by Cys¹⁴⁵ and His⁴¹ residues (active site). Dimerization yields the functional homodimer. Here, our aim was to investigate small molecules, including lopinavir and ritonavir, α-ketoamide 13b, and ebselen, for their ability to interact with the M^pro. The sirtuin 1 agonist SRT1720 was also used in our analyses. Blind docking to each monomer individually indicated preferential binding of the ligands in the active site. Site-mapping of the dimeric protease indicated a highly reactive pocket in the dimerization region at the domain III apex. Blind docking consistently indicated a strong preference of ligand binding in domain III, away from the active site. Molecular dynamics simulations indicated that ligands docked both to the active site and in the dimerization region at the apex, formed relatively stable interactions. Overall, our findings do not obviate the superior potency with respect to inhibition of protease activity of covalently-linked inhibitors such as α-ketoamide 13b in the M^pro active site. Nevertheless, along with those from others, our findings highlight the importance of further characterisation of the M^pro active site and any potential allosteric sites.     

Pharmacophore model-aided virtual screening combined with comparative molecular docking and molecular dynamics for identification of marine natural products as SARS-CoV-2 papain-like protease inhibitors

Targeting SARS-CoV-2 papain-like protease using inhibitors is a suitable approach for inhibition of virus replication and dysregulation of host anti-viral immunity. Engaging all five binding sites far from the catalytic site of PLpro is essential for developing a potent inhibitor. We developed and validated a structure-based pharmacophore model with 9 features of a potent PLpro inhibitor. The pharmacophore model-aided virtual screening of the comprehensive marine natural product database predicted 66 initial hits. This hit library was downsized by filtration through a molecular weight filter of ≤ 500 g/mol. The 50 resultant hits were screened by comparative molecular docking using AutoDock and AutoDock Vina. Comparative molecular docking enables benchmarking docking and relieves the disparities in the search and scoring functions of docking engines. Both docking engines retrieved 3 same compounds at different positions in the top 1 % rank, hence consensus scoring was applied, through which CMNPD28766, aspergillipeptide F emerged as the best PLpro inhibitor. Aspergillipeptide F topped the 50-hit library with a pharmacophore-fit score of 75.916. Favorable binding interactions were predicted between aspergillipeptide F and PLpro similar to the native ligand XR8-24. Aspergillipeptide F was able to engage all the 5 binding sites including the newly discovered BL2 groove, site V. Molecular dynamics for quantification of Cα-atom movements of PLpro after ligand binding indicated that it exhibits highly correlated domain movements contributing to the low free energy of binding and a stable conformation. Thus, aspergillipeptide F is a promising candidate for pharmaceutical and clinical development as a potent SARS-CoV-2 PLpro inhibitor.     

Docking and molecular dynamics study on the inhibitory activity of coumarins on aldose reductase

In order to explore the inhibitory mechanism of coumarins toward aldose reductase (ALR2), AutoDock and Gromacs software were used for docking and molecular dynamics studies on 14 coumarins (CM) and ALR2 protease. The docking results indicate that residues TYR48, HIS110, and TRP111 construct the active pocket of ALR2 and, besides van der Waals and hydrophobic interaction, CM mainly interact with ALR2 by forming hydrogen bonds to cause inhibitory behavior. Except for CM1, all the other coumarins take the lactone part as acceptor to build up the hydrogen bond network with active-pocket residues. Unlike CM3, which has two comparable binding modes with ALR2, most coumarins only have one dominant orientation in their binding sites. The molecular dynamics calculation, based on the docking results, implies that the orientations of CM in the active pocket show different stabilities. Orientation of CM1 and CM3a take an unstable binding mode with ALR2; their conformations and RMSDs relative to ALR2 change a lot with the dynamic process. While the remaining CM are always hydrogen-bonded with residues TYR48 and HIS110 through the carbonyl O atom of the lactone group during the whole process, they retain the original binding mode and gradually reach dynamic equilibrium.