In this work, a theoretical model with a numerical solution is brought forward for a bio-nanofluid with varying fluid features over a slippery sheet. The partial differential equations (PDEs) involving temperature-dependent quantities have been translated into ordinary differential equations (ODEs) by using similarity variables. Numerical verifications have been done in three different methods: finite difference method, shooting method, and bvp4c. To figure out the influence of parameters on the flows, the graphs are plotted for the velocity, temperature, concentration, and microorganism curves. The boundary layer thickness of the microorganism profile reduces with the Schmidt number and Peclet number. In addition to adding radiative heat flux, we added heat generation, rate of chemical reaction, and first-order slip. Adding these parameters brought new aspects to the underlying profiles. Moreover, the obtained data of the skin friction coefficient, the local Nusselt number, the local Sherwood number, and the local density of motile microorganisms are tabulated against various parameters for the physical parameters. From the results, it is apparent that the local Nusselt number decreases with the Brownian and thermophoretic parameters. The data obtained for physical parameters have a close agreement with the published data. Finally, the graphs for slip conditions are significantly different when the comparison is drawn with no-slip condition.