A lattice Boltzmann (LB) interfacial gas-solid three-dimensional model is developed for isothermal multicomponent flows with strongly nonequimolar catalytic reactions, further accounting for the presence of velocity slips and concentration jumps. The model includes diffusion coefficients of all reactive species in the calculation of the catalytic reaction rates as well as an updated velocity at the reactive boundary node. Lattice Boltzmann simulations are performed in a catalytic channel-flow geometry under a wide range of Knudsen (Kn) and surface Damköhler (Da_s) numbers. Comparisons with simulations from a computational fluid dynamics (CFD) Navier-Stokes solver show good agreement in the continuum regime (Kn < 0.01) in terms of flow velocity and reactive species distributions, while comparisons with direct simulation Monte Carlo results from the literature attest to the model's applicability in capturing the correct slip velocity at Kn as high as 0.1, even with a significantly reduced number of grid points (N = 10) in the cross-flow direction. Theoretical and numerical results demonstrate that the term Da_s × Kn × A₂ (where A₂ is a function of the mass accommodation coefficient) determines the significance of the concentration jump on the catalytic reaction rate. The developed model is applicable for many catalytic microflow systems with complex geometries (such as reactors with porous networks) and large velocity/concentration slips (such as catalytic microthrusters for space applications).