Modeling the kinetics of the neutralizing antibody response against SARS-CoV-2 variants after several administrations of Bnt162b2

Because SARS-CoV-2 constantly mutates to escape from the immune response, there is a reduction of neutralizing capacity of antibodies initially targeting the historical strain against emerging Variants of Concern (VoC)s. That is why the measure of the protection conferred by vaccination cannot solely rely on the antibody levels, but also requires to measure their neutralization capacity. Here we used a mathematical model to follow the humoral response in 26 individuals that received up to three vaccination doses of Bnt162b2 vaccine, and for whom both anti-S IgG and neutralization capacity was measured longitudinally against all main VoCs. Our model could identify two independent mechanisms that led to a marked increase in measured humoral response over the successive vaccination doses. In addition to the already known increase in IgG levels after each dose, we identified that the neutralization capacity was significantly increased after the third vaccine administration against all VoCs, despite large inter-individual variability. Consequently, the model projects that the mean duration of detectable neutralizing capacity against non-Omicron VoC is between 348 days (Beta variant, 95% Prediction Intervals PI [307; 389]) and 587 days (Alpha variant, 95% PI [537; 636]). Despite the low neutralization levels after three doses, the mean duration of detectable neutralizing capacity against Omicron variants varies between 173 days (BA.5 variant, 95% PI [142; 200]) and 256 days (BA.1 variant, 95% PI [227; 286]). Our model shows the benefit of incorporating the neutralization capacity in the follow-up of patients to better inform on their level of protection against the different SARS-CoV-2 variants. Trial registration: This clinical trial is registered with ClinicalTrials.gov, Trial IDs NCT04750720 and NCT05315583.