Strengths and Weaknesses of Docking Simulations in the SARS-CoV-2 Era: the Main Protease (Mpro) Case Study

Ensemble docking based virtual screening of SARS-CoV-2 main protease inhibitors

During the first years of COVID-19 pandemic, X-ray structures of the coronavirus drug targets were acquired at an unprecedented rate, giving hundreds of PDB depositions in less than a year. The main protease (Mpro) of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) is the primary validated target of direct-acting antivirals. The selection of the optimal ensemble of structures of Mpro for the docking-driven virtual screening campaign was thus non-trivial and required a systematic and automated approach. Here we report a semi-automated active site RMSD based procedure of ensemble selection from the SARS-CoV-2 Mpro crystallographic data and virtual screening of its inhibitors. The procedure was compared with other approaches to ensemble selection and validated with the help of hand-picked and peer-reviewed activity-annotated libraries. Prospective virtual screening of non-covalent Mpro inhibitors resulted in a new chemotype of thienopyrimidinone derivatives with experimentally confirmed enzyme inhibition.

Bioengineered amyloid peptide for rapid screening of inhibitors against main protease of SARS-CoV-2

The coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has evoked a worldwide pandemic. As the emergence of variants has hampered the neutralization capacity of currently available vaccines, developing effective antiviral therapeutics against SARS-CoV-2 and its variants becomes a significant challenge. The main protease (Mpro) of SARS-CoV-2 has received increased attention as an attractive pharmaceutical target because of its pivotal role in viral replication and proliferation. Here, we generated a de novo Mpro-inhibitor screening platform to evaluate the efficacies of Mpro inhibitors based on Mpro cleavage site-embedded amyloid peptide (MCAP)-coated gold nanoparticles (MCAP-AuNPs). We fabricated MCAPs comprising an amyloid-forming sequence and Mpro-cleavage sequence, mimicking in vivo viral replication process mediated by Mpro. By measuring the proteolytic activity of Mpro and the inhibitory efficacies of various drugs, we confirmed that the MCAP-AuNP-based platform was suitable for rapid screening potential of Mpro inhibitors. These results demonstrated that our MCAP-AuNP-based platform has great potential for discovering Mpro inhibitors and may accelerate the development of therapeutics against COVID-19.

An in silico approach to evaluate the bindings of natural flavonoids and RNA-DNA hybrids

Flavonoids, low molecular weight polyphenolic compounds, are important natural products that belong to plant secondary metabolites. They have diverse biomedical applications such as antioxidative, anti-inflammatory, enzyme inhibitory, antimutagenic, anticarcinogenic, aromatase inhibitory effects, etc. Some of the flavonoids have been exported for bindings with certain DNA and tRNA structures both experimentally and computationally. RNA-DNA hybrid (RDH) falls into an important category of noncanonical nucleic acid structures that have many important biological functions. We have investigated the interaction of RDH structures with some of the dietary flavonoids with the aid of computational methods such as docking and molecular dynamics simulation. The presence of the - OH group on the ligand and the availability of a proper binding pocket in the macromolecule are the two main factors driving the binding preference. Thus, this computationally guided report explains the binding of the flavonoids with RDH structures to assist the researchers in designing noncanonical nucleic acid-targeted drug molecules.Communicated by Ramaswamy H. Sarma.

ReGen-DTI: A novel generative drug target interaction model for predicting potential drug candidates against SARS-COV2

Covid-19 has caused massive numbers of infections and fatalities globally. In response, there has been a large-scale experimental and computational research effort to study and develop drugs. Towards this, Deep learning techniques are used for the generation of potential novel drug candidates that are proven to be effective against exploring large molecular search spaces. Recent advances in reinforcement learning in conjunction with generative techniques has proven to be a promising field in the area of drug discovery. In this regard, we propose a generative drug discovery approach using reinforcement techniques for sampling novel molecules that bind to the main protease of SARS-COV2. The generative method reported significant validity scores for the generated novel molecules and captured the underlying features of the training molecules. Further, the model is fine-tuned on existing re-purposed molecules which are active towards specific target proteins based on similarity metrics. Upon fine tuning the model generated 92.71% valid, 93.55% unique, and 100% novel molecules. Unlike previous methods which are dependent on docking procedures, we proposed a deep learning based novel drug target interaction (DTI) model to find the binding affinity between candidate molecules and target protease sequence. Finally, the binding affinity of the generated molecules is predicted against the 3CLPro main protease by using the proposed DTI model. Most of the generated molecules have shown binding affinity scores <100 nM (lower the better), which are significantly better compared to the existing commercial drugs including Remdesevir.

Garbage in, garbage out: how reliable training data improved a virtual screening approach against SARS-CoV-2 MPro

Introduction: The identification of chemical compounds that interfere with SARS-CoV-2 replication continues to be a priority in several academic and pharmaceutical laboratories. Computational tools and approaches have the power to integrate, process and analyze multiple data in a short time. However, these initiatives may yield unrealistic results if the applied models are not inferred from reliable data and the resulting predictions are not confirmed by experimental evidence. Methods: We undertook a drug discovery campaign against the essential major protease (MPro) from SARS-CoV-2, which relied on an in silico search strategy -performed in a large and diverse chemolibrary- complemented by experimental validation. The computational method comprises a recently reported ligand-based approach developed upon refinement/learning cycles, and structure-based approximations. Search models were applied to both retrospective (in silico) and prospective (experimentally confirmed) screening. Results: The first generation of ligand-based models were fed by data, which to a great extent, had not been published in peer-reviewed articles. The first screening campaign performed with 188 compounds (46 in silico hits and 100 analogues, and 40 unrelated compounds: flavonols and pyrazoles) yielded three hits against MPro (IC₅₀ ≤ 25 μM): two analogues of in silico hits (one glycoside and one benzo-thiazol) and one flavonol. A second generation of ligand-based models was developed based on this negative information and newly published peer-reviewed data for MPro inhibitors. This led to 43 new hit candidates belonging to different chemical families. From 45 compounds (28 in silico hits and 17 related analogues) tested in the second screening campaign, eight inhibited MPro with IC₅₀ = 0.12-20 μM and five of them also impaired the proliferation of SARS-CoV-2 in Vero cells (EC₅₀ 7-45 μM). Discussion: Our study provides an example of a virtuous loop between computational and experimental approaches applied to target-focused drug discovery against a major and global pathogen, reaffirming the well-known "garbage in, garbage out" machine learning principle.

Computer-assisted drug repurposing for thymidylate kinase drug target in monkeypox virus

Introduction

Monkeypox is a zoonotic disease caused by brick-shaped enveloped monkeypox (Mpox) virus that belongs to the family of ancient viruses known as Poxviridae. Subsequently, the viruses have been reported in various countries. The virus is transmitted by respiratory droplets, skin lesions, and infected body fluids. The infected patients experience fluid-filled blisters, maculopapular rash, myalgia, and fever. Due to the lack of effective drugs or vaccines, there is a need to identify the most potent and effective drugs to reduce the spread of monkeypox. The current study aimed to use computational methods to quickly identify potentially effective drugs against the Mpox virus.

Methods

In our study, the Mpox protein thymidylate kinase (A48R) was targeted because it is a unique drug target. We screened a library of 9000 FDA-approved compounds of the DrugBank database by using various in silico approaches, such as molecular docking and molecular dynamic (MD) simulation.

Results

Based on docking score and interaction analysis, compounds DB12380, DB13276, DB13276, DB11740, DB14675, DB11978, DB08526, DB06573, DB15796, DB08223, DB11736, DB16250, and DB16335 were predicted as the most potent. To examine the dynamic behavior and stability of the docked complexes, three compounds-DB16335, DB15796, and DB16250 -along with the Apo state were simulated for 300ns. The results revealed that compound DB16335 revealed the best docking score (-9.57 kcal/mol) against the Mpox protein thymidylate kinase.

Discussion

Additionally, during the 300 ns MD simulation period, thymidylate kinase DB16335 showed great stability. Further, in vitro and in vivo study is recommended for the final predicted compounds.

An overview of compound properties, multiparameter optimization, and computational drug design methods for PARP-1 inhibitor drugs

Breast cancer treatment with PARP-1 inhibitors remains challenging due to emerging toxicities, drug resistance, and unaffordable costs of treatment options. How do we invent strategies to design better anti-cancer drugs? A part of the answer is in optimized compound properties, desirability functions, and modern computational drug design methods that drive selectivity and toxicity and have not been reviewed for PARP-1 inhibitors. Nonetheless, comparisons of these compound properties for PARP-1 inhibitors are not available in the literature. In this review, we analyze the physchem, PKPD space to identify inherent desirability functions characteristic of approved drugs that can be valuable for the design of better candidates. Recent literature utilizing ligand, structure-based drug design strategies and matched molecular pair analysis (MMPA) for the discovery of novel PARP-1 inhibitors are also reviewed. Thus, this perspective provides valuable insights into the medchem and multiparameter optimization of PARP-1 inhibitors that might be useful to other medicinal chemists.

Lessons Learnt from COVID-19: Computational Strategies for Facing Present and Future Pandemics

Since its outbreak in December 2019, the COVID-19 pandemic has caused the death of more than 6.5 million people around the world. The high transmissibility of its causative agent, the SARS-CoV-2 virus, coupled with its potentially lethal outcome, provoked a profound global economic and social crisis. The urgency of finding suitable pharmacological tools to tame the pandemic shed light on the ever-increasing importance of computer simulations in rationalizing and speeding up the design of new drugs, further stressing the need for developing quick and reliable methods to identify novel active molecules and characterize their mechanism of action. In the present work, we aim at providing the reader with a general overview of the COVID-19 pandemic, discussing the hallmarks in its management, from the initial attempts at drug repurposing to the commercialization of Paxlovid, the first orally available COVID-19 drug. Furthermore, we analyze and discuss the role of computer-aided drug discovery (CADD) techniques, especially those that fall in the structure-based drug design (SBDD) category, in facing present and future pandemics, by showcasing several successful examples of drug discovery campaigns where commonly used methods such as docking and molecular dynamics have been employed in the rational design of effective therapeutic entities against COVID-19.

The diverse role of heparan sulfate and other GAGs in SARS-CoV-2 infections and therapeutics

In December 2019, the global coronavirus disease 2019 (COVID-19) pandemic began in Wuhan, China. COVID-19 is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which infects host cells primarily through the angiotensin-converting enzyme 2 (ACE2) receptor. In addition to ACE2, several studies have shown the importance of heparan sulfate (HS) on the host cell surface as a co-receptor for SARS-CoV-2-binding. This insight has driven research into antiviral therapies, aimed at inhibiting the HS co-receptor-binding, e.g., by glycosaminoglycans (GAGs), a family of sulfated polysaccharides that includes HS. Several GAGs, such as heparin (a highly sulfated analog of HS), are used to treat various health indications, including COVID-19. This review is focused on current research on the involvement of HS in SARS-CoV-2 infection, implications of viral mutations, as well as the use of GAGs and other sulfated polysaccharides as antiviral agents.

Molecular docking-based virtual screening: Challenges in hits identification for Anti-SARS-Cov-2 activity

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-Cov-2) requires finding new drugs or repurposing drugs for clinical use. Molecular docking belongs to structure-based drug design providing a fast method for identifying the hit compounds with antiviral activity against SARS-Cov-2. However, the weakness of the docking method is compounded by the limited crystallographic information and comparison drugs due to the novelty of this virus can present challenges in identifying hits of anti-SARS-Cov-2. In the current review, we highlighted several aspects, especially those related to the target structure, docking validation, and virtual hit selection, that need to be considered to obtain reliable docking results. Here, we discussed several cases pertaining to the issue highlighted and approaches that could be used to solve them.

Computational molecular interaction between SARS-CoV-2 main protease and theaflavin digallate using free energy perturbation and molecular dynamics

Our objective was to identify the molecule which can inhibit SARS-CoV-2 main protease and can be easily procured. Natural products may provide such molecules and can supplement the current custom chemical synthesis-based drug discovery for this objective. A combination of docking approaches, scoring functions, classical molecular dynamic simulation, binding pose metadynamics, and free energy perturbation calculations have been employed in this study. Theaflavin digallate has been observed in top-scoring compounds after the three independent virtual screening simulations of 598435 compounds (unique 27256 chemical entities). The main protease-theaflavin digallate complex interacts with critical active site residues of the main protease in molecular dynamics simulation independent of the explored computational framework, simulation time, initial structure, and force field used. Theaflavin digallate forms approximately three hydrogen bonds with Glutamate166 of main protease, primarily through hydroxyl groups in the benzene ring of benzo(7)annulen-6-one, along with other critical residues. Glu166 is the most critical amino acid for main protease dimerization, which is necessary for catalytic activity. The estimated binding free energy, calculated by Amber and Schrodinger MMGBSA module, reflects a high binding free energy between theaflavin digallate and main protease. Binding pose metadynamics simulation shows the highly persistent H-bond and a stable pose for the theaflavin digallate-main protease complex. Using method control, experimental controls, and test set, alchemical transformation studies confirm high relative binding free energy of theaflavin digallate with the main protease. Computational molecular interaction suggests that theaflavin digallate can inhibit the main protease of SARS-CoV-2.

Identification of new small molecule monoamine oxidase-B inhibitors through pharmacophore-based virtual screening, molecular docking and molecular dynamics simulation studies

The discovery of a safe and efficacious drug is a complex, time-consuming, and expensive process. Computational methodologies driven by cheminformatics tools play a central role in the high-throughput lead discovery and optimization process especially when the structure of the biological target is known. Monoamine oxidases are the membrane-bound FAD-containing enzymes and the isoform monoamine oxidase-B (MAO-B) is an attractive target for treating diseases like Alzheimer's disease, Parkinson's disease, glioma, etc. In the current study, we have used a pharmacophore-based virtual screening technique for the identification of new small molecule MAO-B inhibitors. Safinamide was used for building a pharmacophore model and the developed model was used to probe the ZINC database for potential hits. The obtained hits were filtered against drug-likeness and PAINS. Out of the hit's library, two compounds ZINC02181408, ZINC08853942 (most active), and ZINC53327382 (least active) were further subjected to molecular docking and dynamics simulation studies to assess their virtual binding affinities and stability of the resultant protein-ligand complex. The docking studies revealed that active ligands were well accommodated within the active site of MAO-B and interacted with both substrate and entrance cavity residues. MD simulation studies unveiled additional hydrogen bond interactions with the substrate cavity residues, Tyr398 and Tyr435 that are crucial for the catalytic role of MAO-B. Moreover, the predicted ADMET parameters suggest that the compounds ZINC08853942 and ZINC02181408 are suitable for CNS penetration. Thus, the attempted computational campaign yielded two potential MAO-B inhibitors that merit further experimental investigation.Communicated by Ramaswamy H. Sarma.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

Discovery of 2-Thiobenzimidazoles as Noncovalent Inhibitors of SARS-CoV-2 Main Protease

The discovery of antiviral agents against SARS-CoV-2 is an important step toward ending the COVID-19 pandemic and to tackle future outbreaks. In this context, the main protease (M^pro) represents an ideal target for developing coronavirus antivirals, being conserved among different strains and essential for survival. In this work, using in silico tools, we created and validated a docking protocol able to predict binders to the catalytic site of M^pro. The following structure-based virtual screening of a subset of the ZINC library (over 4.3 million unique structures), led to the identification of a hit compound having a 2-thiobenzimidazole scaffold. The inhibitory activity was confirmed using a FRET-based proteolytic assay against recombinant M^pro. Structure-activity relationships were obtained with the synthesis of a small library of analogs, guided by the analysis of the docking pose. Our efforts led to the identification of a micromolar M^pro inhibitor (IC₅₀ = 14.9 µM) with an original scaffold possessing ideal drug-like properties (predicted using the QikProp function) and representing a promising lead for the development of a novel class of coronavirus antivirals.

Design, Semisynthesis, and Estrogenic Activity of Lignan Derivatives from Natural Dibenzylbutyrolactones

Based on molecular docking studies on the ERα, a series of lignan derivatives (3–16) were designed and semisynthesized from the natural dibenzylbutyrolactones bursehernin (1) and matairesinol dimethyl ether (2). To examine their estrogenic and antiestrogenic potencies, the effects of these compounds on estrogen receptor element (ERE)-driven reporter gene expression and viability in human ER+ breast cancer cells were evaluated. Lignan compounds induced ERE-driven reporter gene expression with very low potency as compared with the pure agonist E2. However, coincubation of 5 μM of lignan derivatives 1, 3, 4, 7, 8, 9, 11, 13, and 14 with increasing concentrations of E2 (from 0.01 pM to 1 nM) reduced both the potency and efficacy of pure agonists. The binding to the rhERα-LBD was validated by TR-FRET competitive binding assay and lignans bound to the rhERα with IC₅₀ values from 0.16 μM (compound 14) to 6 μM (compound 4). Induced fit docking (IFD) and molecular dynamics (MD) simulations for compound 14 were carried out to further investigate the binding mode interactions. Finally, the in silico ADME predictions indicated that the most potent lignan derivatives exhibited good drug-likeness.

Finding a chink in the armor: Update, limitations, and challenges toward successful antivirals against flaviviruses

Flaviviruses have caused large epidemics and ongoing outbreaks for centuries. They are now distributed in every continent infecting up to millions of people annually and may emerge to cause future epidemics. Some of the viruses from this group cause severe illnesses ranging from hemorrhagic to neurological manifestations. Despite decades of research, there are currently no approved antiviral drugs against flaviviruses, urging for new strategies and antiviral targets. In recent years, integrated omics data-based drug repurposing paired with novel drug validation methodologies and appropriate animal models has substantially aided in the discovery of new antiviral medicines. Here, we aim to review the latest progress in the development of both new and repurposed (i) direct-acting antivirals; (ii) host-targeting antivirals; and (iii) multitarget antivirals against flaviviruses, which have been evaluated both in vitro and in vivo, with an emphasis on their targets and mechanisms. The search yielded 37 compounds that have been evaluated for their efficacy against flaviviruses in animal models; 20 of them are repurposed drugs, and the majority of them exhibit broad-spectrum antiviral activity. The review also highlighted the major limitations and challenges faced in the current in vitro and in vivo evaluations that hamper the development of successful antiviral drugs for flaviviruses. We provided an analysis of what can be learned from some of the approved antiviral drugs as well as drugs that failed clinical trials. Potent in vitro and in vivo antiviral efficacy alone does not warrant successful antiviral drugs; current gaps in studies need to be addressed to improve efficacy and safety in clinical trials.

Computational Repurposing of Drugs and Natural Products Against SARS-CoV-2 Main Protease (Mpro) as Potential COVID-19 Therapies

We urgently need to identify drugs to treat patients suffering from COVID-19 infection. Drugs rarely act at single molecular targets. Off-target effects are responsible for undesirable side effects and beneficial synergy between targets for specific illnesses. They have provided blockbuster drugs, e.g., Viagra for erectile dysfunction and Minoxidil for male pattern baldness. Existing drugs, those in clinical trials, and approved natural products constitute a rich resource of therapeutic agents that can be quickly repurposed, as they have already been assessed for safety in man. A key question is how to screen such compounds rapidly and efficiently for activity against new pandemic pathogens such as SARS-CoV-2. Here, we show how a fast and robust computational process can be used to screen large libraries of drugs and natural compounds to identify those that may inhibit the main protease of SARS-CoV-2. We show that the shortlist of 84 candidates with the strongest predicted binding affinities is highly enriched (≥25%) in compounds experimentally validated in vivo or in vitro to have activity in SARS-CoV-2. The top candidates also include drugs and natural products not previously identified as having COVID-19 activity, thereby providing leads for experimental validation. This predictive in silico screening pipeline will be valuable for repurposing existing drugs and discovering new drug candidates against other medically important pathogens relevant to future pandemics.

Predicting Antifouling Activity and Acetylcholinesterase Inhibition of Marine-Derived Compounds Using a Computer-Aided Drug Design Approach

Biofouling is the undesirable growth of micro- and macro-organisms on artificial water-immersed surfaces, which results in high costs for the prevention and maintenance of this process (billion €/year) for aquaculture, shipping and other industries that rely on coastal and off-shore infrastructure. To date, there are still no sustainable, economical and environmentally safe solutions to overcome this challenging phenomenon. A computer-aided drug design (CADD) approach comprising ligand- and structure-based methods was explored for predicting the antifouling activities of marine natural products (MNPs). In the CADD ligand-based method, 141 organic molecules extracted from the ChEMBL database and literature with antifouling screening data were used to build the quantitative structure–activity relationship (QSAR) classification model. An overall predictive accuracy score of up to 71% was achieved with the best QSAR model for external and internal validation using test and training sets. A virtual screening campaign of 14,492 MNPs from Encinar’s website and 14 MNPs that are currently in the clinical pipeline was also carried out using the best QSAR model developed. In the CADD structure-based approach, the 125 MNPs that were selected by the QSAR approach were used in molecular docking experiments against the acetylcholinesterase enzyme. Overall, 16 MNPs were proposed as the most promising marine drug-like leads as antifouling agents, e.g., macrocyclic lactam, macrocyclic alkaloids, indole and pyridine derivatives.

Computational Identification of Possible Allosteric Sites and Modulators of the SARS-CoV-2 Main Protease

In this study, we target the main protease (M^pro) of the SARS-CoV-2 virus as it is a crucial enzyme for viral replication. Herein, we report three plausible allosteric sites on M^pro that can expand structure-based drug discovery efforts for new M^pro inhibitors. To find these sites, we used mixed-solvent molecular dynamics (MixMD) simulations, an efficient computational protocol that finds binding hotspots through mapping the surface of unbound proteins with 5% cosolvents in water. We have used normal mode analysis to support our claim of allosteric control for these sites. Further, we have performed virtual screening against the sites with 361 hits from M^pro screenings available through the National Center for Advancing Translational Sciences (NCATS). We have identified the NCATS inhibitors that bind to the remote sites better than the active site of M^pro, and we propose these molecules may be allosteric regulators of the system. After identifying our sites, new X-ray crystal structures were released that show fragment molecules in the sites we found, supporting the notion that these sites are accurate and druggable.

Machine learning prediction of 3CLpro SARS-CoV-2 docking scores

Molecular docking results of two training sets containing 866 and 8,696 compounds were used to train three different machine learning (ML) approaches. Neural network approaches according to Keras and TensorFlow libraries and the gradient boosted decision trees approach of XGBoost were used with DScribe’s Smooth Overlap of Atomic Positions molecular descriptors. In addition, neural networks using the SchNetPack library and descriptors were used. The ML performance was tested on three different sets, including compounds for future organic synthesis. The final evaluation of the ML predicted docking scores was based on the ZINC in vivo set, from which 1,200 compounds were randomly selected with respect to their size. The results obtained showed a consistent ML prediction capability of docking scores, and even though compounds with more than 60 atoms were found slightly overestimated they remain valid for a subsequent evaluation of their drug repurposing suitability.

Discovery of C-12 dithiocarbamate andrographolide analogues as inhibitors of SARS-CoV-2 main protease: In vitro and in silico studies

•

Andrographolide analogues were found to inhibit SARS-CoV-2 main protease.

•

The compounds 3k, 3l, 3m and 3t showed promising in vitro inhibitory activity.

•

Most of the candidates could bind well to the SARS-CoV-2 main protease active site.

A global crisis of coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has impacted millions of people’s lives throughout the world. In parallel to vaccine development, identifying potential antiviral agents against SARS-CoV-2 has become an urgent need to combat COVID-19. One of the most attractive drug targets for discovering anti-SARS-CoV-2 agents is the main protease (M^pro), which plays a pivotal role in the viral life cycle. This study aimed to elucidate a series of twenty-one 12-dithiocarbamate-14-deoxyandrographolide analogues as SARS-CoV-2 M^pro inhibitors using in vitro and in silico studies. These compounds were initially screened for the inhibitory activity toward SARS-CoV-2 M^pro by in vitro enzyme-based assay. We found that compounds 3k, 3l, 3m and 3t showed promising inhibitory activity against SARS-CoV-2 M^pro with >50% inhibition at 10 μM. Afterward, the binding mode of each compound in the active site of SARS-CoV-2 M^pro was explored by molecular docking. The optimum docked complexes were then chosen and subjected to molecular dynamic (MD) simulations. The MD results suggested that all studied complexes were stable along the simulation time, and most of the compounds could fit well with the SARS-CoV-2 M^pro active site, particularly at S1, S2 and S4 subsites. The per-residue decomposition free energy calculations indicated that the hot-spot residues essential for ligand binding were T25, H41, C44, S46, M49, C145, H163, M165, E166, L167, D187, R188, Q189 and T190. Therefore, the obtained information from the combined experimental and computational techniques could lead to further optimization of more specific and potent andrographolide analogues toward SARS-CoV-2 M^pro.

Exploration of SARS-CoV-2 3CLpro Inhibitors by Virtual Screening Methods, FRET Detection, and CPE Assay

COVID-19 caused by a novel coronavirus (SARS-CoV-2) has been spreading all over the world since the end of 2019, and no specific drug has been developed yet. 3C-like protease (3CLpro) acts as an important part of the replication of novel coronavirus and is a promising target for the development of anticoronavirus drugs. In this paper, eight machine learning models were constructed using naïve Bayesian (NB) and recursive partitioning (RP) algorithms for 3CLpro on the basis of optimized two-dimensional (2D) molecular descriptors (MDs) combined with ECFP_4, ECFP_6, and MACCS molecular fingerprints. The optimal models were selected according to the results of 5-fold cross verification, test set verification, and external test set verification. A total of 5766 natural compounds from the internal natural product database were predicted, among which 369 chemical components were predicted to be active compounds by the optimal models and the EstPGood values were more than 0.6, as predicted by the NB (MD + ECFP_6) model. Through ADMET analysis, 31 compounds were selected for further biological activity determination by the fluorescence resonance energy transfer (FRET) method and cytopathic effect (CPE) detection. The results indicated that (+)-shikonin, shikonin, scutellarein, and 5,3',4'-trihydroxyflavone showed certain activity in inhibiting SARS-CoV-2 3CLpro with the half-maximal inhibitory concentration (IC50) values ranging from 4.38 to 87.76 μM. In the CPE assay, 5,3',4'-trihydroxyflavone showed a certain antiviral effect with an IC50 value of 8.22 μM. The binding mechanism of 5,3',4'-trihydroxyflavone with SARS-CoV-2 3CLpro was further revealed through CDOCKER analysis. In this study, 3CLpro prediction models were constructed based on machine learning algorithms for the prediction of active compounds, and the activity of potential inhibitors was determined by the FRET method and CPE assay, which provide important information for further discovery and development of antinovel coronavirus drugs.

A Review of the Current Landscape of SARS-CoV-2 Main Protease Inhibitors: Have We Hit the Bullseye Yet?

In this review, we collected 1765 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) M-pro inhibitors from the bibliography and other sources, such as the COVID Moonshot project and the ChEMBL database. This set of inhibitors includes only those compounds whose inhibitory capacity, mainly expressed as the half-maximal inhibitory concentration (IC50) value, against M-pro from SARS-CoV-2 has been determined. Several covalent warheads are used to treat covalent and non-covalent inhibitors separately. Chemical space, the variation of the IC50 inhibitory activity when measured by different methods or laboratories, and the influence of 1,4-dithiothreitol (DTT) are discussed. When available, we have collected the values of inhibition of viral replication measured with a cellular antiviral assay and expressed as half maximal effective concentration (EC50) values, and their possible relationship to inhibitory potency against M-pro is analyzed. Finally, the most potent covalent and non-covalent inhibitors that simultaneously inhibit the SARS-CoV-2 M-pro and the virus replication in vitro are discussed.

Epicatechin is a promising novel inhibitor of SARS-CoV-2 entry by disrupting interactions between angiotensin-converting enzyme type 2 and the viral receptor binding domain: A computational/simulation study

Angiotensin-converting enzyme 2 (ACE2) is the first target of SARS-CoV-2 and a key functional host receptor through which this virus hooks into and infects human cells. The necessity to block this receptor is one of the essential means to prevent the outbreak of COVID-19. This study was conducted to determine the most eligible natural compound to suppress ACE2 to counterfeit its interaction with the viral infection. To do this, the most known compounds of sixty-six Iraqi medicinal plants were generated and retrieved from PubChem database. After preparing a library for Iraqi medicinal plants, 3663 unique ligands’ conformers were docked to ACE2 using the GLIDE tool. Results found that twenty-three compounds exhibited the highest binding affinity with ACE2. The druglikeness and toxicity potentials of these compounds were evaluated using SwissADME and Protox servers respectively. Out of these virtually screened twenty-three compounds, epicatechin and kempferol were predicted to exert the highest druglikeness and lowest toxicity potentials. Extended Molecular dynamics (MD) simulations showed that ACE2-epicatechin complex exhibited a slightly higher binding stability than ACE2-kempferol complex. In addition to the well-known ACE2 inhibitors that were identified in previous studies, this study revealed for the first time that epicatechin from Hypericum perforatum provided a better static and dynamic inhibition for ACE2 with highly favourable pharmacokinetic properties than the other known ACE2 inhibiting compounds. This study entailed the ability of epicatechin to be used as a potent natural inhibitor that can be used to block or at least weaken the SARS-CoV-2 entry and its subsequent invasion. In vitro experiments are required to validate epicatechin effectiveness against the activity of the human ACE2 receptor.