The emergence of SARS-CoV-2 in December
2019 has become a global
issue due to the continuous upsurge in patients and the lack of drug
efficacy for treatment. SARS-CoV-2 3CLPro is one of the
most intriguing biomolecular targets among scientists worldwide for
developing antiviral drugs due to its relevance in viral replication
and transcription. Herein, we utilized computer-assisted drug screening
to investigate 326 natural products from Thai traditional plants using
structure-based virtual screening against SARS-CoV-2 3CLPro. Following the virtual screening, the top 15 compounds based on
binding energy and their interactions with key amino acid Cys145 were
obtained. Subsequently, they were further evaluated for protein–ligand
complex stability via molecular dynamics simulation
and binding free energy calculation using molecular mechanics Poisson–Boltzmann
surface area (MM-PBSA) approaches. Following drug-likeness and ADME/Tox
assessments, seven bisbenzylisoquinolines were obtained, including
neferine (3), liensinine (4), isoliensinine
(5), dinklacorine (8), tiliacorinine (13), 2′-nortiliacorinine (14), and yanangcorinine
(15). These compounds computationally showed a higher
binding affinity than native N3 and GC-373 inhibitors and attained
stable interactions on the active site of 3CLpro during
100 ns in molecular dynamics (MD) simulation. Moreover, the in vitro enzymatic assay showed that most bisbenzylisoquinolines
could experimentally inhibit SARS-CoV-2 3CLPro. To our
delight, isoliensinine (5) isolated from Nelumbo
nucifera demonstrated the highest inhibition of protease
activity with the IC50 value of 29.93 μM with low
toxicity on Vero cells. Our findings suggested that bisbenzylisoquinoline
scaffolds could be potentially used as an in vivo model for the development of effective anti-SARS-CoV-2 drugs.