Modelling the in-host dynamics of Neisseria gonorrhoeae infection

The bacterial species Neisseria gonorrhoeae (NG) has evolved to replicate effectively and exclusively in human epithelia, with its survival dependent on complex interactions between bacteria, host cells and antimicrobial agents. A better understanding of these interactions is needed to inform development of new approaches to gonorrhoea treatment and prevention but empirical studies have proven difficult, suggesting a role for mathematical modelling. Here, we describe an in-host model of progression of untreated male symptomatic urethral infection, including NG growth and interactions with epithelial cells and neutrophils, informed by in vivo and in vitro studies. The model reproduces key observations on bacterial load and clearance and we use multivariate sensitivity analysis to refine plausible ranges for model parameters. Model variants are also shown to describe mouse infection dynamics with altered parameter ranges that correspond to observed differences between human and mouse infection. Our results highlight the importance of NG internalisation, particularly within neutrophils, in sustaining infection in the human model, with ∼80% of the total NG population internalised from day 25 on. This new mechanistic model of in-host NG infection dynamics should also provide a platform for future studies relating to antimicrobial treatment and resistance and infection at other anatomical sites.