Compounds from the TCM Database@Taiwan were docked to H1 and N1 protein active sites reported in our previous study [91] . All procedures were completed under the forcefield of Chemistry at HARvard Molecular Mechanics (CHARMm) [105] . The virtual screening process was performed using LigandFit. The conformational search method was based on the Monte Carlo algorithm. Rigid body minimization following initial ligand placement was completed using Smart Minimizer. Scoring functions used by LigandFit were DockScore.