Materials and Methods: A library of 120 phytochemical ligands was prepared, from which 5 were selected considering their absorption, distribution, metabolism, and excretion (ADMET) and quantitative structure-activity relationship (QSAR) profiles. The protein active sites and belonging quantum tunnels were defined to conduct supramolecular docking of the aforementioned ligands. The hydrogen bond formation and hydrophobic interactions between the ligand-receptor complexes were studied following the molecular docking steps. A comprehensive molecular dynamic simulation (MDS) was conducted for each of the ligand-receptor complexes to figure out the values - root mean square deviation (RMSD) (Å), root mean square fluctuation (RMSF) (Å), H-bonds, Cα, solvent accessible surface area (SASA) (Å2), molecular surface area (MolSA) (Å2), Rg (nm), and polar surface area (PSA) (Å). Finally, computational programming and algorithms were used to interpret the dynamic simulation outputs into their graphical quantitative forms.

Results: ADMET and QSAR profiles revealed that the most active candidates from the library to be used were apigenin, isovitexin, piperolactam A, and quercetin as test ligands, whereas serpentine as the control. Based on the binding affinities of supramolecular docking and the parameters of molecular dynamic simulation, the strength of the test ligands can be classified as isovitexin > quercetin > piperolactam A > apigenin when complexed with the hACE2 receptor. Surprisingly, serpentine showed lower affinity (-8.6 kcal/mol) than that of isovitexin (-9.9 kcal/mol) and quercetin (-8.9 kcal/mol). The MDS analysis revealed all ligands except isovitexin having a value lower than 2.5 Ǻ. All the test ligands exhibited acceptable fluctuation ranges of RMSD (Å), RMSF (Å), H-bonds, Cα, SASA (Å2), MolSA (Å2), Rg (nm), and PSA (Å) values.