The influence of heterosubstituent and ring puckering angle (ϕ) on electronic properties (μ, α, β and γ ) of 16 four-membered heterocycles (1-16) were demonstrated theoretically using ab initio (HF and MP2), density functional theory (DFT) (PBE1PBE, SVWN5 and B3LYP) and semi-empirical (PM6 and PM3) treatments. Comparing the DFT and ab initio methods, these showed a similar description of electric dipole moment (μ) and polarisability (α) for each one of the studied heterocycles, with small differences in the prediction of α. The calculated hyperpolarisabilities (β and γ ) are very sensitive to the quality of level theory, especially when the heterocycles contain phosphorus and sulphur. On the other hand, comparing with the rest of methods used, the semi-empirical calculations show significant deviations in the prediction of the different properties, and only are in agreement with the calculations of μ for the heterocycles that do not contain phosphorus or sulphur. In general, β and γ are more sensitive to the molecular geometry than μ and α. The dipole moment is determined by the heteroatom (X) within the ring, where the oxetanes and thietanes exhibit a minimum value of μ at planar conformations, while the azetidines and phosphetanes at axial conformations. All heterocycles show a minimum value of α at planar conformation, and their largest values are associated to conformations highly puckered. β exhibits an implicit relationship with the C = Y bond according to the conformational states. In addition, the X and Y heterosubstituents affect significantly the values of α, β and γ , observing the largest values for the phosphorus-and sulphur-containing heterocycles. Finally, considering the values obtained for α, β and γ , the studied heterocycles can represent an interesting unit for the design of new promising push-pull systems as donor or acceptor group.