Although the π-stacking interaction between the substituted aromatic ring of small molecule inhibitors and the cytosine nucleobase present at the active site of HIV integrase (IN) are abundantly seen, no current study systematically covers differing positional substituents. Herewith, we report a comprehensive and systematic study aimed at unraveling the nature of substituent effects on model C||X-MB-Y (|| denotes π-stacking interaction, X is CHO, CN, NO 2 , COOH, Cl, SH, SCH 3 , F, NH 2 , OH, OCH 3 , CH 3 and H; MB denotes methylbenzene and Y typifies diverse substituent positioning which includes m 1 , m 2 , o 1 , o 2 and p) based on findings from the methodical incorporation of B3LYP-D3/6-311++G(d,p) binding energy (E bind) as well as symmetry-adopted perturbation theory (SAPT) results. Energy data indicated that stability of the C||X-MB-Y complexes increased in the presence of both the electron-withdrawing substituents (EWSs) and electron-donating substituents (EDSs) relative to the unsubstituted ones (X = H). Moreso, the lack of a good correlation between E bind and σ meta confirmed that the mechanism by which the various substituents facilitated the stability of each complex is not only by electrostatic interactions. In agreement with this finding, it was observed in this study that the combination of different substituent constant terms for all positions gave excellent correlations with E bind. In addition, the results indicated that the strength of π-stacking interactions depends on the corresponding orientations of interacting rings. Also, it was observed that the most stable complexes for EWSs and EDSs corresponded to positions o 1 and o 2 , respectively. Estimations of SAPT revealed the importance of the electrostatic term (E ele) in the prediction of E bind for each position, except at position m 2 , where all energy components played important roles in the stability of the complexes. Also, comparative studies of the different positions indicated that the stability of the complexes decreases by a decreasing E ele. Image 10895525 • The mechanisms by which the various substituents facilitated the stability of each complexes is not only by electrostatic interactions. • The strength of π-stacking interactions also depends on the corresponding orientations of interacting rings. • Estimations of SAPT revealed the importance of electrostatic term (E ele) in the prediction of E bind for each position, except at position m 2. [ABSTRACT FROM AUTHOR]