Predicting drug-Protein interaction with deep learning framework for molecular graphs and sequences: Potential candidates against SAR-CoV-2.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused the COVID-19 disease, which represents a new life-threatening disaster. Regarding viral infection, many therapeutics have been investigated to alleviate the epidemiology such as vaccines and receptor decoys. However, the continuous mutating coronavirus, especially the variants of Delta and Omicron, are tended to invalidate the therapeutic biological product. Thus, it is necessary to develop molecular entities as broad-spectrum antiviral drugs. Coronavirus replication is controlled by the viral 3-chymotrypsin-like cysteine protease (3CLpro) enzyme, which is required for the virus's life cycle. In the cases of severe acute respiratory syndrome coronavirus (SARS-CoV) and middle east respiratory syndrome coronavirus (MERS-CoV), 3CLpro has been shown to be a promising therapeutic development target. Here we proposed an attention-based deep learning framework for molecular graphs and sequences, training from the BindingDB 3CLpro dataset (114,555 compounds). After construction of such model, we conducted large-scale screening the in vivo/vitro dataset (276,003 compounds) from Zinc Database and visualize the candidate compounds with attention score. geometric-based affinity prediction was employed for validation. Finally, we established a 3CLpro-specific deep learning framework, namely GraphDPI-3CL (AUROC: 0.958) achieved superior performance beyond the existing state of the art model and discovered 10 molecules with a high binding affinity of 3CLpro and superior binding mode.

Drug repurposing: identification of SARS-CoV-2 potential inhibitors by virtual screening and pharmacokinetics strategies.

INTRODUCTION: The coronavirus disease 2019 (COVID-19) pandemic caused global health, economic, and population loss. Variants of the coronavirus contributed to the severity of the disease and persistent rise in infections. This study aimed to identify potential drug candidates from fifteen approved antiviral drugs against SARS-CoV-2 (6LU7), SARS-CoV (5B6O), and SARS-CoV-2 spike protein (6M0J) using virtual screening and pharmacokinetics to gain insights into COVID-19 therapeutics. METHODOLOGY: We employed drug repurposing approach to analyze binding performance of fifteen clinically approved antiviral drugs against the main protease of SARS-CoV-2 (6LU7), SARS-CoV (5B6O), and SARS-CoV-2 spike proteins bound to ACE-2 receptor (6M0J), to provide an insight into the therapeutics of COVID-19. AutoDock Vina was used for docking studies. The binding affinities were calculated, and 2-3D structures of protein-ligand interactions were drawn. RESULTS: Rutin, hesperidin, and nelfinavir are clinically approved antiviral drugs with high binding affinity to proteins 6LU7, 5B6O, and 6M0J. These ligands have excellent pharmacokinetics, ensuring efficient absorption, metabolism, excretion, and digestibility. Hesperidin showed the most potent interaction with spike protein 6M0J, forming four H-bonds. Nelfinavir had a high human intestinal absorption (HIA) score of 0.93, indicating maximum absorption in the body and promising interactions with 6LU7. CONCLUSIONS: Our results indicated that rutin, hesperidin, and nelfinavir had the highest binding results against the proposed drug targets. The computational approach effectively identified SARS-CoV-2 inhibitors. COVID-19 is still a recurrent threat globally and predictive analysis using natural compounds might serve as a starting point for new drug development against SARS-CoV-2 and related viruses.

SARS-CoV-2 Mpro responds to oxidation by forming disulfide and NOS/SONOS bonds.

The main protease (Mpro) of SARS-CoV-2 is critical for viral function and a key drug target. Mpro is only active when reduced; turnover ceases upon oxidation but is restored by re-reduction. This suggests the system has evolved to survive periods in an oxidative environment, but the mechanism of this protection has not been confirmed. Here, we report a crystal structure of oxidized Mpro showing a disulfide bond between the active site cysteine, C145, and a distal cysteine, C117. Previous work proposed this disulfide provides the mechanism of protection from irreversible oxidation. Mpro forms an obligate homodimer, and the C117-C145 structure shows disruption of interactions bridging the dimer interface, implying a correlation between oxidation and dimerization. We confirm dimer stability is weakened in solution upon oxidation. Finally, we observe the protein's crystallization behavior is linked to its redox state. Oxidized Mpro spontaneously forms a distinct, more loosely packed lattice. Seeding with crystals of this lattice yields a structure with an oxidation pattern incorporating one cysteine-lysine-cysteine (SONOS) and two lysine-cysteine (NOS) bridges. These structures further our understanding of the oxidative regulation of Mpro and the crystallization conditions necessary to study this structurally.

Synthesis and in silico inhibitory action studies of azo-anchored imidazo[4,5-b]indole scaffolds against the COVID-19 main protease (Mpro).

The present work elicits a novel approach to combating COVID-19 by synthesizing a series of azo-anchored 3,4-dihydroimidazo[4,5-b]indole derivatives. The envisaged methodology involves the L-proline-catalyzed condensation of para-amino-functionalized azo benzene, indoline-2,3-dione, and ammonium acetate precursors with pertinent aryl aldehyde derivatives under ultrasonic conditions. The structures of synthesized compounds were corroborated through FT-IR, 1H NMR, 13C NMR, and mass analysis data. Molecular docking studies assessed the inhibitory potential of these compounds against the main protease (Mpro) of SARS-CoV-2. Remarkably, in silico investigations revealed significant inhibitory action surpassing standard drugs such as Remdesivir, Paxlovid, Molnupiravir, Chloroquine, Hydroxychloroquine (HCQ), and (N3), an irreversible Michael acceptor inhibitor. Furthermore, the highly active compound was also screened for cytotoxicity activity against HEK-293 cells and exhibited minimal toxicity across a range of concentrations, affirming its favorable safety profile and potential suitability. The pharmacokinetic properties (ADME) of the synthesized compounds have also been deliberated. This study paves the way for in vitro and in vivo testing of these scaffolds in the ongoing battle against SARS-CoV-2.

SARS-CoV-2 Mpro oligomerization as a potential target for therapy.

The main protease (Mpro) of SARS-CoV-2 is critical in the virus's replication cycle, facilitating the maturation of polyproteins into functional units. Due to its conservation across taxa, Mpro is a promising target for broad-spectrum antiviral drugs. Targeting Mpro with small molecule inhibitors, such as nirmatrelvir combined with ritonavir (Paxlovid™), which the FDA has approved for post-exposure treatment and prophylaxis, can effectively interrupt the replication process of the virus. A key aspect of Mpro's function is its ability to form a functional dimer. However, the mechanics of dimerization and its influence on proteolytic activity remain less understood. In this study, we utilized biochemical, structural, and molecular modelling approaches to explore Mpro dimerization. We evaluated critical residues, specifically Arg4 and Arg298, that are essential for dimerization. Our results show that changes in the oligomerization state of Mpro directly affect its enzymatic activity and dimerization propensity. We discovered a synergistic relationship influencing dimer formation, involving both intra- and intermolecular interactions. These findings highlight the potential for developing allosteric inhibitors targeting Mpro, offering promising new directions for therapeutic strategies.

Statine-based peptidomimetic compounds as inhibitors for SARS-CoV-2 main protease (SARS-CoV­2 Mpro).

COVID-19 is a multisystemic disease caused by the SARS-CoV-2 airborne virus, a member of the Coronaviridae family. It has a positive sense single-stranded RNA genome and encodes two non-structural proteins through viral cysteine-proteases processing. Blocking this step is crucial to control virus replication. In this work, we reported the synthesis of 23 statine-based peptidomimetics to determine their ability to inhibit the main protease (Mpro) activity of SARS-CoV-2. Among the 23 peptidomimetics, 15 compounds effectively inhibited Mpro activity by 50% or more, while three compounds (7d, 8e, and 9g) exhibited maximum inhibition above 70% and IC50 < 1 µM. Compounds 7d, 8e, and 9g inhibited roughly 80% of SARS-CoV-2 replication and proved no cytotoxicity. Molecular docking simulations show putative hydrogen bond and hydrophobic interactions between specific amino acids and these inhibitors. Molecular dynamics simulations further confirmed the stability and persisting interactions in Mpro's subsites, exhibiting favorable free energy binding (ΔGbind) values. These findings suggest the statine-based peptidomimetics as potential therapeutic agents against SARS-CoV-2 by targeting Mpro.

Deep profiling of potential substrate atlas of porcine epidemic diarrhea virus 3C-like protease.

Coronavirus (CoV) 3C-like protease (3CLpro) is essential for viral replication and is involved in immune escape by proteolyzing host proteins. Deep profiling the 3CLpro substrates in the host proteome extends our understanding of viral pathogenesis and facilitates antiviral drug discovery. Here, 3CLpro from porcine epidemic diarrhea virus (PEDV), an enteropathogenic CoV, was used as a model which to identify the potential 3CLpro cleavage motifs in all porcine proteins. We characterized the selectivity of PEDV 3CLpro at sites P5-P4'. We then compiled the 3CLpro substrate preferences into a position-specific scoring matrix and developed a 3CLpro profiling strategy to delineate the protein substrate landscape of CoV 3CLpro. We identified 1,398 potential targets in the porcine proteome containing at least one putative cleavage site and experimentally validated the reliability of the substrate degradome. The PEDV 3CLpro-targeted pathways are involved in mRNA processing, translation, and key effectors of autophagy and the immune system. We also demonstrated that PEDV 3CLpro suppresses the type 1 interferon (IFN-I) cascade via the proteolysis of multiple signaling adaptors in the retinoic acid-inducible gene I (RIG-I) signaling pathway. Our composite method is reproducible and accurate, with an unprecedented depth of coverage for substrate motifs. The 3CLpro substrate degradome establishes a comprehensive substrate atlas that will accelerate the investigation of CoV pathogenicity and the development of anti-CoV drugs.IMPORTANCECoronaviruses (CoVs) are major pathogens that infect humans and animals. The 3C-like protease (3CLpro) encoded by CoV not only cleaves the CoV polyproteins but also degrades host proteins and is considered an attractive target for the development of anti-CoV drugs. However, the comprehensive characterization of an atlas of CoV 3CLpro substrates is a long-standing challenge. Using porcine epidemic diarrhea virus (PEDV) 3CLpro as a model, we developed a method that accurately predicts the substrates of 3CLpro and comprehensively maps the substrate degradome of PEDV 3CLpro. Interestingly, we found that 3CLpro may simultaneously degrade multiple molecules responsible for a specific function. For instance, it cleaves at least four adaptors in the RIG-I signaling pathway to suppress type 1 interferon production. These findings highlight the complexity of the 3CLpro substrate degradome and provide new insights to facilitate the development of anti-CoV drugs.

In vitro biological evaluation and in silico insights into the antiviral activity of standardized olive leaves extract against SARS-CoV-2.

There is still a great global need for efficient treatments for the management of SARS-CoV-2 illness notwithstanding the availability and efficacy of COVID-19 vaccinations. Olive leaf is an herbal remedy with a potential antiviral activity that could improve the recovery of COVID-19 patients. In this work, the olive leaves major metabolites were screened in silico for their activity against SARS-CoV-2 by molecular docking on several viral targets such as methyl transferase, helicase, Plpro, Mpro, and RdRp. The results of in silico docking study showed that olive leaves phytoconstituents exhibited strong potential antiviral activity against SARS-CoV-2 selected targets. Verbacoside demonstrated a strong inhibition against methyl transferase, helicase, Plpro, Mpro, and RdRp (docking scores = -17.2, -20, -18.2, -19.8, and -21.7 kcal/mol.) respectively. Oleuropein inhibited 5rmm, Mpro, and RdRp (docking scores = -15, -16.6 and -18.6 kcal/mol., respectively) respectively. Apigenin-7-O-glucoside exhibited activity against methyl transferase and RdRp (docking score = -16.1 and -19.4 kcal/mol., respectively) while Luteolin-7-O-glucoside inhibited Plpro and RdRp (docking score = -15.2 and -20 kcal/mol., respectively). The in vitro antiviral assay was carried out on standardized olive leaf extract (SOLE) containing 20% oleuropein and IC50 was calculated. The results revealed that 20% SOLE demonstrated a moderate antiviral activity against SARS-CoV-2 with IC50 of 118.3 µg /mL. Accordingly, olive leaf could be a potential herbal therapy against SARS-CoV-2 but more in vivo and clinical investigations are recommended.

Anti-Coronavirus Potential of Polyether Ionophores: The New Application of Veterinary Antibiotics in Livestock.

Coronaviruses have consistently posed a major global concern in the field of livestock industry and public health. However, there is currently a lack of efficient drugs with broad-spectrum antiviral activity to address the challenges presented by emerging mutated strains or drug resistance. Additionally, the method for identifying multitarget drugs is also insufficient. Aminopeptidase N (APN) and 3C-like proteinase (3CLpro) represent promising targets for host-directed and virus-directed strategies, respectively, in the development of effective drugs against various coronaviruses. In this study, maduramycin ammonium demonstrated a broad-spectrum antiviral effect by targeting both of the proteins. The binding domains 4 Å from the ligand of both target proteins shared a structural similarity, suggesting that screening and designing drugs based on these domains might exhibit broad-spectrum and highly effective antiviral activity. Furthermore, it was identified that the polyether ionophores' ability to carry zinc ion might be one of the reasons why they were able to target APN and exhibit antiviral effect. The findings of this experiment provide novel perspectives for future drug screening and design, while also offering valuable references for the utilization of polyether ionophores in the management of livestock health.

De novo design of potential peptide analogs against the main protease of Omicron variant using in silico studies.

SARS-CoV-2 and its variants are crossing the immunity barrier induced through vaccination. Recent Omicron sub-variants are highly transmissible and have a low mortality rate. Despite the low severity of Omicron variants, these new variants are known to cause acute post-infectious syndromes. Nowadays, novel strategies to develop new potential inhibitors for SARS-CoV-2 and other Omicron variants have gained prominence. For viral replication and survival the main protease of SARS-CoV-2 plays a vital role. Peptide-like inhibitors that mimic the substrate peptide have already proved to be effective in inhibiting the Mpro of SARS-CoV-2 variants. Our systematic canonical amino acid point mutation analysis on the native peptide has revealed various ways to improve the native peptide of the main protease. Multi mutation analysis has led us to identify and design potent peptide-analog inhibitors that act against the Mpro of the Omicron sub-variants. Our in-depth analysis of all-atom molecular dynamics studies has paved the way to characterize the atomistic behavior of Mpro in Omicron variants. Our goal is to develop potent peptide-analogs that could be therapeutically effective against Omicron and its sub-variants.

Discovery of First-in-Class PROTAC Degraders of SARS-CoV-2 Main Protease.

We have witnessed three coronavirus (CoV) outbreaks in the past two decades, including the COVID-19 pandemic caused by SARS-CoV-2. Main protease (MPro), a highly conserved protease among various CoVs, is essential for viral replication and pathogenesis, making it a prime target for antiviral drug development. Here, we leverage proteolysis targeting chimera (PROTAC) technology to develop a new class of small-molecule antivirals that induce the degradation of SARS-CoV-2 MPro. Among them, MPD2 was demonstrated to effectively reduce MPro protein levels in 293T cells, relying on a time-dependent, CRBN-mediated, and proteasome-driven mechanism. Furthermore, MPD2 exhibited remarkable efficacy in diminishing MPro protein levels in SARS-CoV-2-infected A549-ACE2 cells. MPD2 also displayed potent antiviral activity against various SARS-CoV-2 strains and exhibited enhanced potency against nirmatrelvir-resistant viruses. Overall, this proof-of-concept study highlights the potential of targeted protein degradation of MPro as an innovative approach for developing antivirals that could fight against drug-resistant viral variants.

Exploring the Efficacy of Noncovalent SARS-CoV-2 Main Protease Inhibitors: A Computational Simulation Analysis Study.

The SARS-CoV-2 main protease, as a key target for antiviral therapeutics, is instrumental in maintaining virus stability, facilitating translation, and enabling the virus to evade innate immunity. Our research focused on designing non-covalent inhibitors to counteract the action of this protease. Utilizing a 3D-QSAR model and contour map, we successfully engineered eight novel non-covalent inhibitors. Further evaluation and comparison of these novel compounds through methodologies including molecular docking, ADMET analysis, frontier molecular orbital studies, molecular dynamics simulations, and binding free energy revealed that the inhibitors N02 and N03 demonstrated superior research performance (N02 ΔGbind=-206.648 kJ/mol, N03 ΔGbind=-185.602 kJ/mol). These findings offer insightful guidance for the further refinement of molecular structures and the development of more efficacious inhibitors. Consequently, future investigations can draw upon these findings to unearth more potent inhibitors, thereby amplifying their impact in the treatment and prevention of associated diseases.

In Silico and In Vitro Screening of Some Pregnane Glycosides Isolated from Certain Caralluma Species as SARS-COV-2 Main Protease Inhibitors.

SARS-CoV-2 caused pandemic represented a major risk for the worldwide human health, animal health and economy, forcing extraordinary efforts to discover drugs for its prevention and cure. Considering the extensive interest in the pregnane glycosides because of their diverse structures and excellent biological activities, we investigated them as antiviral agents against SARS-COV-2. We selected 21 pregnane glycosides previously isolated from the genus Caralluma from Asclepiadaceae family to be tested through virtual screening molecular docking simulations for their potential inhibition of SARS-CoV-2 Mpro. Almost all target compounds showed a more or equally negative docking energy score relative to the co-crystallized inhibitor X77 (S=-12.53 kcal/mol) with docking score range of (-12.55 to -19.76 kcal/mol) and so with a potent predicted binding affinity to the target enzyme. The activity of the most promising candidates was validated by in vitro testing. Arabincoside C showed the highest activity (IC50=35.42 µg/ml) and the highest selectivity index (SI=9.9) followed by Russelioside B (IC50=50.80 µg/ml), and Arabincoside B (IC50=53.31 µg/ml).

Synthesis, molecular docking analysis, molecular dynamic simulation, ADMET, DFT, and drug likeness studies: Novel Indeno[1,2-b]pyrrol-4(1H)-one as SARS-CoV-2 main protease inhibitors.

BACKGROUND: The COVID-19 pandemic began in 2019 as a result of the advent of a novel coronavirus, SARS-CoV-2. At present, there are a limited number of approved antiviral agents for the treatment of COVID-19. Remdesivir, Molnupiravir, and Paxlovid have been approved by the FDA to treat COVID-19 infections. Research has shown that the main protease enzyme (Mpro) of SARS-CoV-2 plays a crucial role in the enzymatic processing of viral polyproteins. This makes Mpro an interesting therapeutic target for combating infections caused by emerging coronaviruses. METHODS: The pharmacological effects of pyrroles and their derivatives have a wide range of applications. In our study, we focused on synthesizing nine novel derivatives of 2-arylamino-dihydro-indeno[1,2-b] pyrrol-4(1H)-one, with a particular emphasis on their antiviral properties. Using in silico studies involving molecular docking and DFT analyses in the gas phase using the B3LYP/6-31++G(d,p) basis set, we studied these compounds with respect to their interactions with the Mpro of SARS-CoV-2. The results of the docking analysis revealed that the synthesized compounds exhibited favorable inhibitory effects. Notably, compound 5f demonstrated the highest effectiveness against the target protein. Furthermore, the pharmacokinetic and drug-like properties of the synthesized derivatives of 2-arylamino-dihydroindeno[1,2-b] pyrrol-4(1H)-one indicated their potential as promising candidates for further development as inhibitors targeting SARS-CoV-2. However, it is imperative to determine the in vitro efficacy of these compounds through comprehensive biochemical and structural analyses.

Sensitive detection of SARS-CoV-2 main protease 3CLpro with an engineered ribonuclease zymogen.

Alongside vaccines and antiviral therapeutics, diagnostic tools are a crucial aid in combating the COVID-19 pandemic caused by the etiological agent SARS-CoV-2. All common assays for infection rely on the detection of viral sub-components, including structural proteins of the virion or fragments of the viral genome. Selective pressure imposed by human intervention of COVID-19 can, however, induce viral mutations that decrease the sensitivity of diagnostic assays based on biomolecular structure, leading to an increase in false-negative results. In comparison, mutations are unlikely to alter the function of viral proteins, and viral machinery is under less selective pressure from vaccines and therapeutics. Accordingly, diagnostic assays that rely on biomolecular function can be more robust than ones that rely on biopolymer structure. Toward this end, we used a split intein to create a circular ribonuclease zymogen that is activated by the SARS-CoV-2 main protease, 3CLpro . Zymogen activation by 3CLpro leads to a >300-fold increase in ribonucleolytic activity, which can be detected with a highly sensitive fluorogenic substrate. This coupled assay can detect low nanomolar concentrations of 3CLpro within a timeframe comparable to that of common antigen-detection protocols. More generally, the concept of detecting a protease by activating a ribonuclease could be the basis of diagnostic tools for other indications.

Prediction of Binding Pose and Affinity of Nelfinavir, a SARS-CoV-2 Main Protease Repositioned Drug, by Combining Docking, Molecular Dynamics, and Fragment Molecular Orbital Calculations.

A novel in silico drug design procedure is described targeting the Main protease (Mpro) of the SARS-CoV-2 virus. The procedure combines molecular docking, molecular dynamics (MD), and fragment molecular orbital (FMO) calculations. The binding structure and properties of Mpro were predicted for Nelfinavir (NFV), which had been identified as a candidate compound through drug repositioning, targeting Mpro. Several poses of the Mpro and NFV complexes were generated by docking, from which four docking poses were selected by scoring with FMO energy. Then, each pose was subjected to MD simulation, 100 snapshot structures were sampled from each of the generated MD trajectories, and the structures were evaluated by FMO calculations to rank the pose based on binding energy. Several residues were found to be important in ligand recognition, including Glu47, Asp48, Glu166, Asp187, and Gln189, all of which interacted strongly with NFV. Asn142 is presumably regarded to form hydrogen bonds or CH/π interaction with NFV; however, in the present calculation, their interactions were transient. Moreover, the tert-butyl group of NFV had no interaction with Mpro. Identifying such strong and weak interactions provides candidates for maintaining and substituting ligand functional groups and important suggestions for drug discovery using drug repositioning. Besides the interaction between NFV and the amino acid residues of Mpro, the desolvation effect of the binding pocket also affected the ranking order. A similar procedure of drug design was applied to Lopinavir, and the calculated interaction energy and experimental inhibitory activity value trends were consistent. Our approach provides a new guideline for structure-based drug design starting from a candidate compound whose complex crystal structure has not been obtained.

Exploiting high-energy hydration sites for the discovery of potent peptide aldehyde inhibitors of the SARS-CoV-2 main protease with cellular antiviral activity.

Small-molecule antivirals that prevent the replication of the SARS-CoV-2 virus by blocking the enzymatic activity of its main protease (Mpro) are and will be a tenet of pandemic preparedness. However, the peptidic nature of such compounds often precludes the design of compounds within favorable physical property ranges, limiting cellular activity. Here we describe the discovery of peptide aldehyde Mpro inhibitors with potent enzymatic and cellular antiviral activity. This structure-activity relationship (SAR) exploration was guided by the use of calculated hydration site thermodynamic maps (WaterMap) to drive potency via displacement of waters from high-energy sites. Thousands of diverse compounds were designed to target these high-energy hydration sites and then prioritized for synthesis by physics- and structure-based Free-Energy Perturbation (FEP+) simulations, which accurately predicted biochemical potencies. This approach ultimately led to the rapid discovery of lead compounds with unique SAR that exhibited potent enzymatic and cellular activity with excellent pan-coronavirus coverage.

Improved fluorescence-based assay for rapid screening and evaluation of SARS-CoV-2 main protease inhibitors.

The outbreak of coronavirus disease 2019 (COVID-19) pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a global threat to human health. In parallel with vaccines, efficacious antivirals are urgently needed. SARS-CoV-2 main protease (Mpro) is an attractive drug target for antiviral development owing to its key roles in virus replication and host immune evasion. Due to the limitations of currently available methods, the development of novel high-throughput screening assays is of the highest importance for the discovery of Mpro inhibitors. In this study, we first developed an improved fluorescence-based assay for rapid screening of Mpro inhibitors from an anti-infection compound library using a versatile dimerization-dependent red fluorescent protein (ddRFP) biosensor. Utilizing this assay, we identified MG-101 as a competitive Mpro inhibitor in vitro. Moreover, our results revealed that ensitrelvir is a potent Mpro inhibitor, but baicalein, chloroquine, ebselen, echinatin, and silibinin are not. Therefore, this robust ddRFP assay provides a faithful avenue for rapid screening and evaluation of Mpro inhibitors to fight against COVID-19.

SAR, Molecular Docking and Molecular Dynamic Simulation of Natural Inhibitors against SARS-CoV-2 Mpro Spike Protein.

The SARS-CoV-2 virus and its mutations have affected human health globally and created significant danger for the health of people all around the world. To cure this virus, the human Angiotensin Converting Enzyme-2 (ACE2) receptor, the SARS-CoV-2 main protease (Mpro), and spike proteins were found to be likely candidates for the synthesis of novel therapeutic drug. In the past, proteins were capable of engaging in interaction with a wide variety of ligands, including both manmade and plant-derived small molecules. Pyrus communis L., Ginko bibola, Carica papaya, Syrian rue, and Pimenta dioica were some of the plant species that were studied for their tendency to interact with SARS-CoV-2 main protease (Mpro) in this research project (6LU7). This scenario investigates the geometry, electronic, and thermodynamic properties computationally. Assessing the intermolecular forces of phytochemicals with the targets of the SARS-CoV-2 Mpro spike protein (SP) resulted in the recognition of a compound, kaempferol, as the most potent binding ligand, -7.7 kcal mol-1. Kaempferol interacted with ASP-187, CYS-145, SER-144, LEU 141, MET-165, and GLU-166 residues. Through additional molecular dynamic simulations, the stability of ligand-protein interactions was assessed for 100 ns. GLU-166 remained intact with 33% contact strength with phenolic OH group. We noted a change in torsional conformation, and the molecular dynamics simulation showed a potential variation in the range from 3.36 to 7.44 against a 45-50-degree angle rotation. SAR, pharmacokinetics, and drug-likeness characteristic investigations showed that kaempferol may be the suitable candidate to serve as a model for designing and developing new anti-COVID-19 medicines.