Cross sections of proton-induced nuclear reactions on enriched 86Sr target were measured by the activation technique up to proton energies of 44 MeV. The isomeric cross-section ratios for 86 m , gY and 85 m , gY as a function of projectile energy were deduced from their measured data. The present experimental data for the nuclear reaction products, namely 86 mY, 86 g + xmY, 85 mY, 85 gY, 84Rb and 83Rb were compared with the results of nuclear model calculations using the code TALYS, which combines the statistical, precompound, and direct interactions. In general, the experimental cross-section data as well as the isomeric cross-section ratios are reproduced well by the model calculations, provided the input model parameters are properly chosen and the level structure of the product nucleus is thoughtfully considered. The quality of the agreement between experimental data and model calculations was numerically quantified. For products formed via emission of a light complex particle as well as multi-nucleons (e.g., α and 2p2n), the contribution of the latter process starts increasing when its energy threshold is crossed.