The human immunodeficiency virus type 1 (HIV-1), one of the leading causes of infectious death globally, generates severe damages to people’s immune systems and makes them susceptible to serious diseases. To date, there are no drugs that completely remove HIV from the body. This paper focuses on screening 224,205 natural compounds of ZINC15 NPs subset to identify those with bioactivity similar to non-nucleoside reverse transcriptase inhibitors (NNRTIs) as promising candidates to treat HIV-1. To reach the goal, an in silico approach involving 3D-similarity search, ADMETox, HIV protein-inhibitor prediction, docking, and MM-GBSA free-binding energies was trained. The FDA-approved HIV drugs, efavirenz, etravirine, rilpivirine, and doravirine, were used as queries. The prioritized compounds were subjected to ADMETox, docking, and MM-GBSA studies against HIV-1 reverse transcriptase (RT). Lys101, Tyr181, Tyr188, Trp229, and Tyr318 residues and free-binding energies have proved that ligands can stably bind to HIV-1 RT. Three natural products (ZINC37538901, ZINC38321654, and ZINC67912677) containing oxan and oxolan rings with hydroxyl substituents and one (ZINC2103242) having 3,6,7,8-tetrahydro-2H-pyrido[1 ,2-a]pyrazine-1,4-dione core exhibited comparable profiles to etravirine and doravirine, with ZINC2103242 being the most promising anti-HIV candidate in terms of drug metabolism and safety profile. These findings may open new avenues to guide the rational design of novel HIV-1 NNRTIs.