The cytokine MCSF is a key activator of osteoclastogenesis, yet analysis of its effects via mathematical models of bone remodeling is more sparse than IL-6. is the first model that includes the role of MCSF. This model investigates the depth and duration of osteoclast erosion during resorption and models MCSF as a scalar variable but assumes it is always present. A sensitivity analysis of resorption depth indicates that changes in MCSF levels are equivalent to changes in maximum osteoclast activity, and both effects are minor compared to TGF-β. and are biomechanical models that include MCSF activation of uncommitted osteoclasts. However, they assume the MCSF concentration is constant, resulting in a continuous activator function term. investigates the kinetics of MCSF via network-based ODEs. As with and , this work investigates the effect of mechanical loading and the effects of PTH and circadian rhythm. In contrast, the network analysis includes multiple parameters that capture the role of MCSF in remodeling (outlined in Table ). Sensitivity analysis of the model shows that the secretion rate of MCSF by osteoblasts and preosteoclasts and the MCSF degradation rate results in a change of more than 5% in bone mass. When the rate of degradation doubles, bone mass increases by more than 60%. One limitation of this study is the assumption that MCSF secretion rates are considered equal across cell types. Still, these results warrant further mathematical investigation of MCSF in homeostatic and pathological bone remodeling. Ultimately, future mathematical models of bone biology should explore the complex and coordinated role of cytokines, growth factors, and hormones in bone remodeling.