In heavy-ion fusion reactions, the enhancement in the sub-barrier fusion cross-section has been observed as compared to the 1-Dimensional barrier penetration model due to the coupling of many degrees of freedom to the relative motion. This enhancement can be explained theoretically by including nuclear structure effects like deformation and the coupling of relative motion among two colliding nuclei. The present work aims to investigate the effect of individual rotational energy levels on the fusion cross-sections for 16O-based reaction systems, namely, 16O + 182,184,186W, 16O + 176,180Hf, 16O + 174,176Yb, 16O + 166Er, 16O + 148,152,154Sm, 16O + 150Nd at energies below the fusion barrier. Using the CCFULL code, the effect of low-lying rotational energy levels on the fusion cross-section for 16O induced reactions has been investigated at energies below and around the Coulomb barrier. The calculations are performed by assuming the fixed value of diffuseness parameter a 0 = 0.65 fm in the Woods-Saxon nuclear potential and the other two parameters are optimised by fitting the experimental data at the above barrier. Here we have determined the V 0 and r 0 as a function of Z P Z T , where experimental cross-sections are available. From our calculations, it is observed that the hexadecapole deformation (β 4) with different magnitudes has a significant influence on the fusion cross sections. For the case of the +ve value of β 4, beyond 10+, the rotational levels cease to contribute significantly and also there is a significant difference between the contribution of sequential channels. On the other hand, in the case of -ve β 4, up to 6+ levels contribute significantly. Furthermore, we have established an algebraic systematic of fitting, which one can use to determine the parameters V 0, r 0 of Woods-Saxon nuclear potential within the range of Z P Z T lie in between 480 ≤ Z P Z T ≤ 592.