A search for charmless four-body decays of $Lambda_{b}^{0}$ and $Xi_{b}^{0}$ baryons with a proton and three charged mesons (either kaons or pions) in the final state is performed. The data sample used was recorded in 2011 and 2012 with the LHCb experiment and corresponds to an integrated luminosity of 3 fb$^{-1}$. Six decay modes are observed, among which $Lambda_{b}^{0} to pK^{-}pi^{+}pi^{-}$, $Lambda_{b}^{0} to pK^{-}K^{+}K^{-}$, $Xi_{b}^{0} to pK^{-}pi^{+}pi^{-}$ and $Xi_{b}^{0} to pK^{-} pi^{+}K^{-}$ are established for the first time. Their branching fractions (including the ratio of hadronisation fractions in the case of the $Xi_{b}^{0}$ baryon) are determined relative to the $Lambda_{b}^{0} to Lambda_{c}^{+}pi^{-}$ decay.