Large-scale ab initio calculations of the electronic contribution to the electric quadrupole hyperfine constant B were performed for the 5s25p6s1,3P1o excited states of neutral tin. To probe the sensitivity of B to different electron correlation effects, three sets of variational multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction calculations employing different strategies were carried out. In addition, a fourth set of calculations was based on the configuration interaction Dirac-Fock-Sturm theory. For the 5s25p6sP1o1 state, the final value of B/Q=703(50) MHz/b differs by 0.4% from the one recently used by Yordanov et al. [Commun. Phys. 3, 107 (2020)10.1038/s42005-020-0348-9] to extract the nuclear quadrupole moments Q for tin isotopes in the range Sn117-131 from collinear laser spectroscopy measurements. Efforts were made to provide a realistic theoretical uncertainty for the final B/Q value of the 5s25p6sP1o1 state based on statistical principles and on correlation with the electronic contribution to the magnetic dipole hyperfine constant A.