This paper is about to formulate a design of predator–prey model with constant and time fractional variable order. The predator and prey act as agents in an ecosystem in this simulation. We focus on a time fractional order Atangana–Baleanu operator in the sense of Liouville–Caputo. Due to the nonlocality of the method, the predator–prey model is generated by using another FO derivative developed as a kernel based on the generalized Mittag-Leffler function. Two fractional-order systems are assumed, with and without delay. For the numerical solution of the models, we not only employ the Adams–Bashforth–Moulton method but also explore the existence and uniqueness of these schemes. We use the fixed point theorem which is useful in describing the existence of a new approach with a particular set of solutions. For the illustration, several numerical examples are added to the paper to show the effectiveness of the numerical method.