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

Abstract

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.

Keywords

• COVID-19
• density functional theory (DFT)
• molecular docking (MD)
• molecular dynamics simulation
• structure–activity relationship (SAR)
• pharmacokinetics