A model is developed and used to study within-human malaria parasite dynamics. The model integrates actors involved in the development-progression of parasitemia, gametocytogenesis and mechanisms for immune response activation. Model analyses under immune suppression reveal different dynamical behaviours for different healthy red blood cell (HRBC) generation functions. Existence of a threshold parameter determines conditions for HRBCs depletion. Oscillatory dynamics reminiscent of malaria parasitemia are obtained. A dependence exists on the type of recruitment function used to generate HRBCs, with complexities observed for a more nonlinear function. An upper bound that delimits the size of feasible parasitized steady-state solution exists for a logistic function but not a constant function. The upper bound is completely characterized and is affected by parameters associated with HRBCs recruitment, parasitized red blood cells generation and the release and time-to-release of free merozoites. A stable density size for mature gametocytes, the bridge to invertebrate hosts, is derived.
Keywords: Gametocytogenesis; Global stability; Innate and adaptive immune response; Malaria; Parasitemia; Recruitment; Red blood cells; Within-human–host dynamics.