The analysis of bioconvection flow nanofluids is the topic of concern in recent decades as it involves a variety of physical significance in biotechnology. Bioconvection has many applications in the interdisciplinary field of sciences such as in biomedical science, biofuel biotechnology, and enzyme-based biosensors, among others. The aim of the current work is to analyze the bioconvection phenomenon in the two-dimensional steady flow of viscoelastic nanofluid over a vertical surface. Here, the effects of activation energy, second-order slip, and nanoparticles zero mass flux conditions are considered to investigate the flow problem. Based on dimensionless variables, the governing partial differential equations (PDEs) are transformed into ordinary differential equations (ODEs) which are further solved numerically by using a built-in BVP4C approach in MATLAB software. Various controlling parameters like Hartman number, viscoelastic parameter, first and second-order slip factors, buoyancy ratio parameter, thermophoresis parameter, Brownian motion constant, bioconvection Lewis number and Peclet number are graphically illustrated for the distributions of velocity, temperature, concentration, and motile microorganism. Moreover, the variation of local Nusselt number, local Sherwood number, and motile density number are numerically investigated for the involved parameters.