Diagnostic and therapeutic purpose on human healthcare is moving to a trend of combined engineering principles and design concepts. It is therefore beneficial to refine the study on thermodynamic properties of bionanofluid towards a better state of diagnoses and treatments. This paper aims to determine the characteristics of a magnetohydrodynamics (MHD) bionanofluid slip flow model in a fluid-permeable plate. A Buongiorno mathematical model describing steady two-dimensional MHD bionanofluid slip flow across a permeable plate is constructed. Numerical analysis by shooting technique and sensitivity analysis under the principle of RSM is configured to determine the correlation between input parameters and output responses. The relationship between the four parameters of interest namely skin friction coefficient, local Nusselt number, local Sherwood number and local density of motile microorganism with respective velocity, temperature, nanoparticle concentration and microorganism density distributions to the selected input parameters are investigated. Skin friction coefficient and local density of motile microorganism increase with suction and Schmidt number respectively, signifying positive sensitivity. Local Sherwood number is conditionally more sensitive to Brownian motion parameter or else Lewis number. While, local Nusselt number is always gaining negative sensitivity from Lewis number, thermophoresis and Brownian motion parameters. Extensive numerical studies on boundary layer problem have been done; yet this is the first study invoking sensitivity analysis to discuss the sensitivity of parameters of interest to the input thermophysical parameters.