The key objective of the present paper is to propose a numerical scheme based on the homotopy analysis transform technique to analyze the time-fractional nonlinear predator-prey population model. The population model is coupled fractional order nonlinear partial differential equations often employed to narrate the dynamics of biological systems in which two species interact, first is a predator and the second is a prey. The proposed scheme provides the series solution with great freedom and flexibility by choosing appropriate parameters. The convergence of results is free from small or large parameters. Three examples are discussed to demonstrate the correctness and efficiency of the used computational approach.