Johne’s disease (JD) is a chronic disease in ruminants and is caused by infection with Mycobacterium avium subspecies paratuberculosis (MAP). At late stages of the disease, MAP bacilli are shed via feces excretion and in turn create the potential for oral-fecal transmission. The role of the host immune response in MAP bacteria shedding patterns at different stages of JD is still unclear. We employed mathematical modeling to predict if the variation in MAP shedding could be correlated to the immune response in infected animals. We used a novel inverse modeling approach that assumed biological interactions among the antigen-specific lymphocyte proliferation response (cell-mediated response), antibody/humoral immune responses, and MAP bacteria. The modeling framework was used to predict and test possible biological interactions between the measured variables and returns only the essential interactions that are relevant in explaining the observed cattle MAP experimental infection data. Through confronting the models with data, we predicted observed effects (enhancement or suppression) and extents of interactions among the three variables. This analysis enabled classification of the infected cattle into three different groups that correspond to the unique predicted immune responses that are essential to explain the data from cattle within these groups. Our analysis highlights the strong and weak points of the modeling approach, as well as the key immune mechanisms predicted to be expressed in all animals and those that were different between animals, hence giving insight into how animals exhibit different disease dynamics and bacteria shedding patterns.