We present a mathematical model that explores the progression of Alzheimer’s disease, with a particular focus on the involvement of disease-related proteins and astrocytes. Our model consists of a coupled system of differential equations that delineates the dynamics of amyloid beta plaques, amyloid beta protein, tau protein, and astrocytes. Amyloid beta plaques can be considered fibrils that depend on both the plaque size and time. We change our mathematical model to a temporal system by applying an integration operation with respect to the plaque size. Theoretical analysis including existence, uniqueness, positivity, and boundedness is performed in our model. We extend our mathematical model by adding two populations, namely prion protein and amyloid beta-prion complex. We characterize the system dynamics by locating biologically feasible steady states and their local stability analysis for both models. The characterization of the proposed model can help inform in advancing our understanding of the development of Alzheimer’s disease as well as its complicated dynamics. We investigate the global stability analysis around the interior equilibrium point by constructing a suitable Lyapunov function. We validate our theoretical analysis with the aid of extensive numerical illustrations.