Respiratory syncytial virus (RSV) is a common virus that can have varying effects ranging from mild cold-like symptoms to mortality depending on the age and immune status of the individual. We combined mathematical modelling using ordinary differential equations (ODEs) with measurement of RSV infection kinetics in primary well-differentiated human bronchial epithelial cultures in vitro and in immunocompetent and immunosuppressed cotton rats to glean mechanistic details that underlie RSV infection kinetics in the lung. Quantitative analysis of viral titre kinetics in our mathematical model showed that the elimination of infected cells by the adaptive immune response generates unique RSV titre kinetic features including a faster timescale of viral titre clearance than viral production, and a monotonic decrease in the peak RSV titre with decreasing inoculum dose. Parameter estimation in the ODE model using a nonlinear mixed effects approach revealed a very low rate (average single-cell lifetime > 10 days) of cell lysis by RSV before the adaptive immune response is initiated. Our model predicted negligible changes in the RSV titre kinetics at early times post-infection (less than 5 dpi) but a slower decay in RSV titre in immunosuppressed cotton rats compared to that in non-suppressed cotton rats at later times (greater than 5 dpi) in silico. These predictions were in excellent agreement with the experimental results. Our combined approach quantified the importance of the adaptive immune response in suppressing RSV infection in cotton rats, which could be useful in testing RSV vaccine candidates.
Keywords: adaptive immunity; cotton rats; nonlinear mixed effects modelling; ordinary differential equations; respiratory syncytial virus; viral kinetics.