Sildenafil, one of selective phosphodiesterase 5 (PDE5) inhibitors, is a widely used oral agent for the treatment of erectile dysfunction. To develop new PDE5 inhibitors with improved therapeutic efficacy, a series of sildenafil analogues have been prepared and their in vitro PDE5 inhibitory activities were evaluated. Their IC(50) values ranged from 423 to 0.05 nM. Herein, the results of 3D-QSAR (CoMFA and CoMSIA) analyses on these inhibitors are reported. Both CoMFA and CoMSIA gave reliable models with q(2) values >0.75 and r(2) values >0.99. The resulting CoMFA and CoMSIA models reveal a good correlation between the contour maps and the active site residues critical for the interaction with inhibitor, and nicely predict the key structural features of new analogues with improved activity and selectivity.