A Quantitative Structure -Activity Relationship (QSAR) analysis was carried out on a dataset of 135 six-and seven-membered cyclic urea-based Human Immunodeficiency Virus Type 1 (HIV-1) protease inhibitors. Using a larger and more diverse dataset than previous studies reported in literature allowed a more comprehensive analysis. A large set of molecular descriptors, calculated with CODESSA PRO, was used in Multiple Linear Regression (MLR) analysis and the resulting four-parameter model enabled accurate prediction of inhibitory activity for the structures both in training and external validation sets. The squared correlation coefficient for both sets combined is 0.824. Descriptors in the final model characterize the size and shape of the molecule and its charge distribution. Non-linear relationships between the HIV-1 protease inhibition activity and some descriptors are significant and described with simple non-linear transformations. After the model development, another dataset of 45 compounds was collected and additional external validation was performed.