Exploring the possibilities of collectivity in small systems stands as one of the major goals for the upcoming p–O and O–O collisions at the LHC. Anisotropic flow coefficients, \(v_\textrm{n} \) which quantify the azimuthal momentum anisotropy of the produced hadrons are known to be sensitive to the initial-state nuclear configurations of the colliding nuclei. Since \(^{16}\mathrm O{ and}^{12}\mathrm C{ nucleiaretheorizedtopossessan}\alpha -\) clustered nuclear geometry, studying p–O and p–C collisions at TeV energies serve the purposes of exploring both small system collectivity as well as the effects of clustered nuclear geometry on the final-state flow coefficients. This work, for the first time, simulates p–O and p–C collisions using a multi-phase transport model, at \(\sqrt{s_\textrm{NN}}=9.9~{ TeVfor}\alpha \) -clustered and Sum-of-Gaussians nuclear density profiles for \(^{16}\mathrm O{ and}^{12}\mathrm C\) nuclei. We observe a very significant dependence of charged particle yield, eccentricities, and flow coefficients on nuclear density profile, collision species, and centrality.