Efficient optimization of an inactive 2-anilinyl-5-benzyloxadiazole core has been guided by free energy perturbation (FEP) calculations to provide potent non-nucleoside inhibitors of human immunodeficiency virus (HIV) reverse transcriptase (NNRTIs). An FEP "chlorine scan" was performed to identify the most promising sites for substitution of aryl hydrogens. This yielded NNRTIs 8 and 10 with activities (EC50) of 820 and 310 nM for protection of human T-cells from infection by wild-type HIV-1. FEP calculations for additional substituent modifications and change of the core heterocycle readily led to oxazoles 28 and 29, which were confirmed as highly potent anti-HIV agents with activities in the 10-20 nM range. The designed compounds were also monitored for possession of desirable pharmacological properties by use of additional computational tools. Overall, the trends predicted by the FEP calculations were well borne out by the assay results. FEP-guided lead optimization is confirmed as a valuable tool for molecular design including drug discovery; chlorine scans are particularly attractive since they are both straightforward to perform and highly informative.