Pyrrolyl aryl sulfones (PASs) have been recently reported as a new class of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) inhibitors acting at the non-nucleoside binding site of this enzyme (Artico, M.; et al. J. Med. Chem. 1996, 39, 522-530). Compound 3, the most potent inhibitor within the series (EC(50) = 0.14 microM, IC(50) = 0.4 microM, and SI > 1429), was then selected as a lead compound for a synthetic project based on molecular modeling studies. Using the three-dimensional structure of RT cocrystallized with the alpha-APA derivative R95845, we derived a model of the RT/3 complex by taking into account previously developed structure-activity relationships. Inspection of this model and docking calculations on virtual compounds prompted the design of novel PAS derivatives and related analogues. Our computational approach proved to be effective in making qualitative predictions, that is in discriminating active versus inactive compounds. Among the compounds synthesized and tested, 20 was the most active one, with EC(50) = 0.045 microM, IC(50) = 0.05 microM, and SI = 5333. Compared with the lead 3, these values represent a 3- and 8-fold improvement in the cell-based and enzyme assays, respectively, together with the highest selectivity achieved so far in the PAS series.