As the world becomes ever more connected, the chance of pandemics increases as well. The recent COVID-19 pandemic and the concurrent global mass vaccine roll-out provides an ideal setting to learn from and refine our understanding of infectious disease models for better future preparedness. In this review, we systematically analyze and categorize mathematical models that have been developed to design optimal vaccine prioritization strategies of an initially limited vaccine. As older individuals are disproportionately affected by COVID-19, the focus is on models that take age explicitly into account. The lower mobility and activity level of older individuals gives rise to non-trivial trade-offs. Secondary research questions concern the optimal time interval between vaccine doses and spatial vaccine distribution. This review showcases the effect of various modeling assumptions on model outcomes. A solid understanding of these relationships yields better infectious disease models and thus public health decisions during the next pandemic.