Backward bifurcation and optimal control in a co-infection model for SARS-CoV-2 and ZIKV

In co-infection models for two diseases, it is mostly claimed that, the dynamical behavior of the sub-models usually predict or drive the behavior of the complete models. However, under a certain assumption such as, allowing incident co-infection with both diseases, we have a different observation. In this paper, a new mathematical model for SARS-CoV-2 and Zika co-dynamics is presented which incorporates incident co-infection by susceptible individuals. It is worth mentioning that the assumption is missing in many existing co-infection models. We shall discuss the impact of this assumption on the dynamics of a co-infection model. The model also captures sexual transmission of Zika virus. The positivity and boundedness of solution of the proposed model are studied, in addition to the local asymptotic stability analysis. The model is shown to exhibit backward bifurcation caused by the disease-induced death rates and parameters associated with susceptibility to a second infection by those singly infected. Using Lyapunov functions, the disease free and endemic equilibria are shown to be globally asymptotically stable for [Formula: see text] , respectively. To manage the co-circulation of both infections effectively, under an endemic setting, time dependent controls in the form of SARS-CoV-2, Zika and co-infection prevention strategies are incorporated into the model. The simulations show that SARS-CoV-2 prevention could greatly reduce the burden of co-infections with Zika. Furthermore, it is also shown that prevention controls for Zika can significantly decrease the burden of co-infections with SARS-CoV-2.