Computational Evaluation and Multi-Criteria Optimization of Natural Compound Analogs Targeting SARS-CoV-2 Proteases

TL;DR

This study uses computational methods to identify and optimize natural compound analogs that could target SARS-CoV-2 proteases, offering potential new therapeutic options for treating COVID-19.

Contribution

The study introduces a multi-criteria in silico framework for evaluating and optimizing natural compound analogs against SARS-CoV-2 proteases.

Findings

CHEMBL1720210 showed strong interaction with PLpro with a binding score of −9.34 kcal/mol.

CHEMBL1495225 exhibited high affinity for 3CLpro with a binding score of −8.04 kcal/mol.

CHEMBL4069090 emerged as a lead compound with favorable drug-likeness and predicted binding to PLpro.

Abstract

The global impact of the COVID-19 crisis has underscored the need for novel therapeutic candidates capable of efficiently targeting essential viral proteins. Existing therapeutic strategies continue to encounter limitations such as reduced efficacy against emerging variants, safety concerns, and suboptimal pharmacodynamics, which emphasize the potential of natural-origin compounds as supportive agents with immunomodulatory, anti-inflammatory, and antioxidant benefits. The present study significantly advances prior molecular docking research through comprehensive virtual screening of structurally related analogs derived from antiviral phytochemicals. These compounds were evaluated specifically against the SARS-CoV-2 main protease (3CLpro) and papain-like protease (PLpro). Utilizing chemical similarity algorithms via the ChEMBL database, over 600 candidate molecules were retrieved and subjected to automated docking, interaction pattern analysis, and comprehensive ADMET profiling. Several analogs showed enhanced binding scores relative to their parent scaffolds, with CHEMBL1720210 (a shogaol-derived analog) demonstrating strong interaction with PLpro (−9.34 kcal/mol), and CHEMBL1495225 (a 6-gingerol derivative) showing high affinity for 3CLpro (−8.04 kcal/mol). Molecular interaction analysis revealed that CHEMBL1720210 forms hydrogen bonds with key PLpro residues including GLY163, LEU162, GLN269, TYR265, and TYR273, complemented by hydrophobic interactions with TYR268 and PRO248. CHEMBL1495225 establishes multiple hydrogen bonds with the 3CLpro residues ASP197, ARG131, TYR239, LEU272, and GLY195, along with hydrophobic contacts with LEU287. Gene expression predictions via DIGEP-Pred indicated that the top-ranked compounds could influence biological pathways linked to inflammation and oxidative stress, processes implicated in COVID-19’s pathology. Notably, CHEMBL4069090 emerged as a lead compound with favorable drug-likeness and predicted binding to PLpro. Overall, the applied in silico framework facilitated the rational prioritization of bioactive analogs with promising pharmacological profiles, supporting their advancement toward experimental validation and therapeutic exploration against SARS-CoV-2.

Figures