Modelling Virus and Antibody Dynamics during Dengue Virus Infection Suggests a Role for Antibody in Virus Clearance

Within-host models unravelling the dynamics of dengue reinfections

Caused by four serotypes, dengue fever is a major public health concern worldwide. Current modeling efforts have mostly focused on primary and heterologous secondary infections, assuming that lifelong immunity prevents reinfections by the same serotype. However, recent findings challenge this assumption, prompting a reevaluation of dengue immunity dynamics. In this study, we develop a within-host modeling framework to explore different scenarios of dengue infections. Unlike previous studies, we go beyond a deterministic framework, considering individual immunological variability. Both deterministic and stochastic models are calibrated using empirical data on viral load and antibody (IgM and IgG) concentrations for all dengue serotypes, incorporating confidence intervals derived from stochastic realizations. With good agreement between the mean of the stochastic realizations and the mean field solution for each model, our approach not only successfully captures primary and heterologous secondary infection dynamics facilitated by antibody-dependent enhancement (ADE) but also provides, for the first time, insights into homotypic reinfection dynamics. Our study discusses the relevance of homotypic reinfections in dengue transmission at the population level, highlighting potential implications for disease prevention and control strategies.

How robust are estimates of key parameters in standard viral dynamic models?

Mathematical models of viral infection have been developed, fitted to data, and provide insight into disease pathogenesis for multiple agents that cause chronic infection, including HIV, hepatitis C, and B virus. However, for agents that cause acute infections or during the acute stage of agents that cause chronic infections, viral load data are often collected after symptoms develop, usually around or after the peak viral load. Consequently, we frequently lack data in the initial phase of viral growth, i.e., when pre-symptomatic transmission events occur. Missing data may make estimating the time of infection, the infectious period, and parameters in viral dynamic models, such as the cell infection rate, difficult. However, having extra information, such as the average time to peak viral load, may improve the robustness of the estimation. Here, we evaluated the robustness of estimates of key model parameters when viral load data prior to the viral load peak is missing, when we know the values of some parameters and/or the time from infection to peak viral load. Although estimates of the time of infection are sensitive to the quality and amount of available data, particularly pre-peak, other parameters important in understanding disease pathogenesis, such as the loss rate of infected cells, are less sensitive. Viral infectivity and the viral production rate are key parameters affecting the robustness of data fits. Fixing their values to literature values can help estimate the remaining model parameters when pre-peak data is missing or limited. We find a lack of data in the pre-peak growth phase underestimates the time to peak viral load by several days, leading to a shorter predicted growth phase. On the other hand, knowing the time of infection (e.g., from epidemiological data) and fixing it results in good estimates of dynamical parameters even in the absence of early data. While we provide ways to approximate model parameters in the absence of early viral load data, our results also suggest that these data, when available, are needed to estimate model parameters more precisely.

The Mathematical Modeling of the Host-Virus Interaction in Dengue Virus Infection: A Quantitative Study

Infectious diseases, such as Dengue fever, pose a significant public health threat. Developing a reliable mathematical model plays a crucial role in quantitatively elucidating the kinetic characteristics of antibody-virus interactions. By integrating previous models and incorporating the antibody dynamic theory, we have constructed a novel and robust model that can accurately simulate the dynamics of antibodies and viruses based on a comprehensive understanding of immunology principles. It explicitly formulates the viral clearance effect of antibodies, along with the positive feedback stimulation of virus-antibody complexes on antibody regeneration. In addition to providing quantitative insights into the dynamics of antibodies and viruses, the model exhibits a high degree of accuracy in capturing the kinetics of viruses and antibodies in Dengue fever patients. This model offers a valuable solution to modeling the differences between primary and secondary Dengue infections concerning IgM/IgG antibodies. Furthermore, it demonstrates that a faster removal rate of antibody-virus complexes might lead to a higher peak viral loading and worse clinical symptom. Moreover, it provides a reasonable explanation for the antibody-dependent enhancement of heterogeneous Dengue infections. Ultimately, this model serves as a foundation for constructing an optimal mathematical model to combat various infectious diseases in the future.

A Novel Mathematical Model That Predicts the Protection Time of SARS-CoV-2 Antibodies

Infectious diseases such as SARS-CoV-2 pose a considerable threat to public health. Constructing a reliable mathematical model helps us quantitatively explain the kinetic characteristics of antibody-virus interactions. A novel and robust model is developed to integrate antibody dynamics with virus dynamics based on a comprehensive understanding of immunology principles. This model explicitly formulizes the pernicious effect of the antibody, together with a positive feedback stimulation of the virus-antibody complex on the antibody regeneration. Besides providing quantitative insights into antibody and virus dynamics, it demonstrates good adaptivity in recapturing the virus-antibody interaction. It is proposed that the environmental antigenic substances help maintain the memory cell level and the corresponding neutralizing antibodies secreted by those memory cells. A broader application is also visualized in predicting the antibody protection time caused by a natural infection. Suitable binding antibodies and the presence of massive environmental antigenic substances would prolong the protection time against breakthrough infection. The model also displays excellent fitness and provides good explanations for antibody selection, antibody interference, and self-reinfection. It helps elucidate how our immune system efficiently develops neutralizing antibodies with good binding kinetics. It provides a reasonable explanation for the lower SARS-CoV-2 mortality in the population that was vaccinated with other vaccines. It is inferred that the best strategy for prolonging the vaccine protection time is not repeated inoculation but a directed induction of fast-binding antibodies. Eventually, this model will inform the future construction of an optimal mathematical model and help us fight against those infectious diseases.

Heterogeneity of Rift Valley fever virus transmission potential across livestock hosts, quantified through a model-based analysis of host viral load and vector infection

Quantifying the variation of pathogens’ life history traits in multiple host systems is crucial to understand their transmission dynamics. It is particularly important for arthropod-borne viruses (arboviruses), which are prone to infecting several species of vertebrate hosts. Here, we focus on how host-pathogen interactions determine the ability of host species to transmit a virus to susceptible vectors upon a potentially infectious contact. Rift Valley fever (RVF) is a viral, vector-borne, zoonotic disease, chosen as a case study. The relative contributions of livestock species to RVFV transmission has not been previously quantified. To estimate their potential to transmit the virus over the course of their infection, we 1) fitted a within-host model to viral RNA and infectious virus measures, obtained daily from infected lambs, calves, and young goats, 2) estimated the relationship between vertebrate host infectious titers and probability to infect mosquitoes, and 3) estimated the net infectiousness of each host species over the duration of their infectious periods, taking into account different survival outcomes for lambs. Our results indicate that the efficiency of viral replication, along with the lifespan of infectious particles, could be sources of heterogeneity between hosts. Given available data on RVFV competent vectors, we found that, for similar infectious titers, infection rates in the Aedes genus were on average higher than in the Culex genus. Consequently, for Aedes-mediated infections, we estimated the net infectiousness of lambs to be 2.93 (median) and 3.65 times higher than that of calves and goats, respectively. In lambs, we estimated the overall infectiousness to be 1.93 times higher in individuals which eventually died from the infection than in those recovering. Beyond infectiousness, the relative contributions of host species to transmission depend on local ecological factors, including relative abundances and vector host-feeding preferences. Quantifying these contributions will ultimately help design efficient, targeted, surveillance and vaccination strategies.

Viruses spread by mosquitoes present a major threat to animal and public health worldwide. When these pathogenic viruses can infect multiple species, controlling their spread becomes difficult. Rift Valley fever virus (RVFV) is such a virus. It spreads predominantly among ruminant livestock but can also spill over and cause severe disease in humans. Understanding which of these ruminant species are most important for the transmission of RVFV can help for effective control. One piece of this puzzle is to assess how effective infected animals are at transmitting RVFV to mosquitoes. To answer this question, we combine mathematical models with observations from experimental infections in cattle, sheep, and goats, and model changes in viremia over time within individuals. We then quantify the relationship between hosts’ viremia and the probability to infect mosquitoes. In combining these two analyses, we estimate the overall transmission potential of sheep, when in contact with mosquitoes, to be 3 to 5 times higher than that of goats and cattle. Further, sheep that experience a lethal infection have an even larger overall transmission potential. Once applied at the level of populations, with setting-specific herd composition and exposure to mosquitoes, these results will help unravel species’ role in RVF outbreaks.

Mathematical models for dengue fever epidemiology: A 10-year systematic review

Mathematical models have a long history in epidemiological research, and as the COVID-19 pandemic progressed, research on mathematical modeling became imperative and very influential to understand the epidemiological dynamics of disease spreading. Mathematical models describing dengue fever epidemiological dynamics are found back from 1970. Dengue fever is a viral mosquito-borne infection caused by four antigenically related but distinct serotypes (DENV-1 to DENV-4). With 2.5 billion people at risk of acquiring the infection, it is a major international public health concern. Although most of the cases are asymptomatic or mild, the disease immunological response is complex, with severe disease linked to the antibody-dependent enhancement (ADE) - a disease augmentation phenomenon where pre-existing antibodies to previous dengue infection do not neutralize but rather enhance the new infection. Here, we present a 10-year systematic review on mathematical models for dengue fever epidemiology. Specifically, we review multi-strain frameworks describing host-to-host and vector-host transmission models and within-host models describing viral replication and the respective immune response. Following a detailed literature search in standard scientific databases, different mathematical models in terms of their scope, analytical approach and structural form, including model validation and parameter estimation using empirical data, are described and analyzed. Aiming to identify a consensus on infectious diseases modeling aspects that can contribute to public health authorities for disease control, we revise the current understanding of epidemiological and immunological factors influencing the transmission dynamics of dengue. This review provide insights on general features to be considered to model aspects of real-world public health problems, such as the current epidemiological scenario we are living in.

Vaccine-Associated Enhanced Disease and Pathogenic Human Coronaviruses

Vaccine-associated enhanced disease (VAED) is a difficult phenomenon to define and can be confused with vaccine failure. Using studies on respiratory syncytial virus (RSV) vaccination and dengue virus infection, we highlight known and theoretical mechanisms of VAED, including antibody-dependent enhancement (ADE), antibody-enhanced disease (AED) and Th2-mediated pathology. We also critically review the literature surrounding this phenomenon in pathogenic human coronaviruses, including MERS-CoV, SARS-CoV-1 and SARS-CoV-2. Poor quality histopathological data and a lack of consistency in defining severe pathology and VAED in preclinical studies of MERS-CoV and SARS-CoV-1 vaccines in particular make it difficult to interrogate potential cases of VAED. Fortuitously, there have been only few reports of mild VAED in SARS-CoV-2 vaccination in preclinical models and no observations in their clinical use. We describe the problem areas and discuss methods to improve the characterisation of VAED in the future.

Modeling Dengue Immune Responses Mediated by Antibodies: A Qualitative Study

Simple Summary

With more than one-third of the world population at risk of acquiring the disease, dengue fever is a major public health problem. Caused by four antigenically distinct but related serotypes, disease severity is associated with the immunological status of the individual, seronegative or seropositive, prior to a natural dengue infection. While a primary natural dengue infection is often asymptomatic or mild, individuals experiencing a secondary dengue infection with a heterologous serotype have higher risk of developing the severe form of the disease, linked to the antibody-dependent enhancement (ADE) process. We develop a modeling framework to describe the dengue immune responses mediated by antibodies. Our model framework can describe qualitatively the dynamic of the viral load and antibodies production for scenarios of primary and secondary infections, as found in the empirical immunology literature. Studies such as the one described here serve as a baseline to further model extensions. Future refinements of our framework will be of use to evaluate the impact of imperfect dengue vaccines.

Abstract

Dengue fever is a viral mosquito-borne infection and a major international public health concern. With 2.5 billion people at risk of acquiring the infection around the world, disease severity is influenced by the immunological status of the individual, seronegative or seropositive, prior to natural infection. Caused by four antigenically related but distinct serotypes, DENV-1 to DENV-4, infection by one serotype confers life-long immunity to that serotype and a period of temporary cross-immunity (TCI) to other serotypes. The clinical response on exposure to a second serotype is complex with the so-called antibody-dependent enhancement (ADE) process, a disease augmentation phenomenon when pre-existing antibodies to previous dengue infection do not neutralize but rather enhance the new infection, used to explain the etiology of severe disease. In this paper, we present a minimalistic mathematical model framework developed to describe qualitatively the dengue immunological response mediated by antibodies. Three models are analyzed and compared: (i) primary dengue infection, (ii) secondary dengue infection with the same (homologous) dengue virus and (iii) secondary dengue infection with a different (heterologous) dengue virus. We explore the features of viral replication, antibody production and infection clearance over time. The model is developed based on body cells and free virus interactions resulting in infected cells activating antibody production. Our mathematical results are qualitatively similar to the ones described in the empiric immunology literature, providing insights into the immunopathogenesis of severe disease. Results presented here are of use for future research directions to evaluate the impact of dengue vaccines.

Optimal bang-bang control for variable-order dengue virus; numerical studies

Introduction

Dengue and Malaria are the most important mosquito-borne viral diseases affecting humans. Fever is transmitted between human hosts by infected female aedes mosquitoes. The modeling study of viral infections is very useful to show how the virus replicates in an infected individual and how the human antibody response acts to control that replication, which antibody playing a key role in controlling infection.

Objectives

Optimal control of a novel variable-order nonlinear model of dengue virus is studied in the present work. Bang-bang control is suggested to minimize the viral infection as well as quick clearance of the virus from the host. Necessary conditions for the control problem are given. The variable-order derivatives are given in the sense of Caputo. Moreover, the parameters of the proposed model are dependent on the same variable-order fractional power. Two numerical schemes are constructed for solving the optimality systems. Comparative studies and numerical simulations are implemented. The variable-order fractional derivative can be describe the effects of long variable memory of time dependent systems than the integer order and fractional order derivatives.

Methods

Both the nonstandard generalized fourth order Runge-Kutta and the nonstandard generalized Euler methods are presented.

Results

We have successfully applied a kind of Pontryagin's maximum principle with bang-bang control and were able to reduce the viraemia level by adding the dose of DI particles. The nonstandard generalized fourth order Runge-Kutta method has the best results than nonstandard generalized Euler method.

Conclusion

The combination of the variable-order fractional derivative and bang-bang control in the Dengue mathematical model improves the dynamics of the model. The nonstandard generalized Euler method and the nonstandard generalized fourth order Runge-Kutta method can be used to study the variable order fractional optimal control problem simply.

Efficacy profile of the CYD-TDV dengue vaccine revealed by Bayesian survival analysis of individual-level phase III data

Background:

Sanofi-Pasteur’s CYD-TDV is the only licensed dengue vaccine. Two phase three trials showed higher efficacy in seropositive than seronegative recipients. Hospital follow-up revealed increased hospitalisation in 2–5- year-old vaccinees, where serostatus and age effects were unresolved.

Methods:

We fit a survival model to individual-level data from both trials, including year 1 of hospital follow-up. We determine efficacy by age, serostatus, serotype and severity, and examine efficacy duration and vaccine action mechanism.

Results:

Our modelling indicates that vaccine-induced immunity is long-lived in seropositive recipients, and therefore that vaccinating seropositives gives higher protection than two natural infections. Long-term increased hospitalisation risk outweighs short-lived immunity in seronegatives. Independently of serostatus, transient immunity increases with age, and is highest against serotype 4. Benefit is higher in seropositives, and risk enhancement is greater in seronegatives, against hospitalised disease than against febrile disease.

Conclusions:

Our results support vaccinating seropositives only. Rapid diagnostic tests would enable viable ‘screen-then-vaccinate’ programs. Since CYD-TDV acts as a silent infection, long-term safety of other vaccine candidates must be closely monitored.

Funding:

Bill & Melinda Gates Foundation, National Institute for Health Research, UK Medical Research Council, Wellcome Trust, Royal Society.

Clinical trial number:

NCT01373281 and NCT01374516.

Modeling the Relationship Between Antibody-Dependent Enhancement and Disease Severity in Secondary Dengue Infection

Sequential infections with different dengue serotypes (DENV-1, 4) significantly increase the risk of a severe disease outcome (fever, shock, and hemorrhagic disorders). Two hypotheses have been proposed to explain the severity of the disease: (1) antibody-dependent enhancement (ADE) and (2) original T cell antigenic sin. In this work, we explored the first hypothesis through mathematical modeling. The proposed model reproduces the dynamic of susceptible and infected target cells and dengue virus in scenarios of infection-neutralizing and infection-enhancing antibody competition induced by two distinct serotypes of the dengue virus during secondary infection. The enhancement and neutralization functions are derived from basic concepts of chemical reactions and used to mimic binding to the virus by two distinct populations of antibodies. The analytic study of the model showed the existence of two equilibriums: a disease-free equilibrium and an endemic one. Using the concept of the basic reproduction number [Formula: see text], we performed the asymptotic stability analysis for the two equilibriums. To measure the severity of the disease, we considered the maximum value of infected cells as well as the time when this maximum is reached. We observed that it corresponds to the time when the maximum enhancing activity for the infection occurs. This critical time was calculated from the model to be a few days after the occurrence of the infection, which corresponds to what is observed in the literature. Finally, using as output [Formula: see text], we were able to rank the contribution of each parameter of the model. In particular, we highlighted that the cross-reactive antibody responses may be responsible for the disease enhancement during secondary heterologous dengue infection.

Dynamic causal modelling of immune heterogeneity

An interesting inference drawn by some COVID-19 epidemiological models is that there exists a proportion of the population who are not susceptible to infection-even at the start of the current pandemic. This paper introduces a model of the immune response to a virus. This is based upon the same sort of mean-field dynamics as used in epidemiology. However, in place of the location, clinical status, and other attributes of people in an epidemiological model, we consider the state of a virus, B and T-lymphocytes, and the antibodies they generate. Our aim is to formalise some key hypotheses as to the mechanism of resistance. We present a series of simple simulations illustrating changes to the dynamics of the immune response under these hypotheses. These include attenuated viral cell entry, pre-existing cross-reactive humoral (antibody-mediated) immunity, and enhanced T-cell dependent immunity. Finally, we illustrate the potential application of this sort of model by illustrating variational inversion (using simulated data) of this model to illustrate its use in testing hypotheses. In principle, this furnishes a fast and efficient immunological assay-based on sequential serology-that provides a (1) quantitative measure of latent immunological responses and (2) a Bayes optimal classification of the different kinds of immunological response (c.f., glucose tolerance tests used to test for insulin resistance). This may be especially useful in assessing SARS-CoV-2 vaccines.

Modeling the intracellular dynamics of the dengue viral infection and the innate immune response

The interplay between the dengue virus and the innate immune response is not fully understood. Here, we use deterministic and stochastic approaches to investigate the dynamics of the interaction between the interferon-mediated innate immune response and the dengue virus. We aim to develop a quantitative representation of these complex interactions and predict their system-level dynamics. Our simulation results predict bimodal and bistable dynamics that represent viral clearance and virus-producing states. Under normal conditions, we determined that the viral infection outcome is modulated by the innate immune response and the positive-strand viral RNA concentration. Additionally, we tested system perturbations by external stimulation, such as the direct induction of the innate immune response by interferon, and a therapeutic intervention consisting of the direct application of mRNA encoding for several interferon-stimulated genes. Our simulation results suggest optimal regimes for the studied intervention approaches.

Virological and Immunological Outcomes in Rhesus Monkeys after Exposure to Dengue Virus-Infected Aedes aegypti Mosquitoes

This study describes the natural history of dengue virus (DENV) infection in rhesus monkeys exposed to the bites of DENV-infected Aedes aegypti mosquitoes. Dengue virus–infected mosquitoes were generated by either intrathoracic inoculation or by oral feeding on viremic blood meals. Each of the six rhesus monkeys that were fed upon by intrathoracically infected mosquitoes developed non-structural protein 1 (NS1) antigenemia and an IgM response; viremia was detected in 4/6 individuals. No virological or immunological evidence of DENV infection was detected in the three monkeys exposed to mosquitoes that had been orally infected with DENV. These results demonstrate the utility of mosquito-borne challenge of rhesus monkeys with DENV.

A Coupled Mathematical Model of the Intracellular Replication of Dengue Virus and the Host Cell Immune Response to Infection

Dengue virus (DV) is a positive-strand RNA virus of the Flavivirus genus. It is one of the most prevalent mosquito-borne viruses, infecting globally 390 million individuals per year. The clinical spectrum of DV infection ranges from an asymptomatic course to severe complications such as dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS), the latter because of severe plasma leakage. Given that the outcome of infection is likely determined by the kinetics of viral replication and the antiviral host cell immune response (HIR) it is of importance to understand the interaction between these two parameters. In this study, we use mathematical modeling to characterize and understand the complex interplay between intracellular DV replication and the host cells' defense mechanisms. We first measured viral RNA, viral protein, and virus particle production in Huh7 cells, which exhibit a notoriously weak intrinsic antiviral response. Based on these measurements, we developed a detailed intracellular DV replication model. We then measured replication in IFN competent A549 cells and used this data to couple the replication model with a model describing IFN activation and production of IFN stimulated genes (ISGs), as well as their interplay with DV replication. By comparing the cell line specific DV replication, we found that host factors involved in replication complex formation and virus particle production are crucial for replication efficiency. Regarding possible modes of action of the HIR, our model fits suggest that the HIR mainly affects DV RNA translation initiation, cytosolic DV RNA degradation, and naïve cell infection. We further analyzed the potential of direct acting antiviral drugs targeting different processes of the DV lifecycle in silico and found that targeting RNA synthesis and virus assembly and release are the most promising anti-DV drug targets.

Modelling the impact of antibody-dependent enhancement on disease severity of Zika virus and dengue virus sequential and co-infection

Human infections with viruses of the genus Flavivirus, including dengue virus (DENV) and Zika virus (ZIKV), are of increasing global importance. Owing to antibody-dependent enhancement (ADE), secondary infection with one Flavivirus following primary infection with another Flavivirus can result in a significantly larger peak viral load with a much higher risk of severe disease. Although several mathematical models have been developed to quantify the virus dynamics in the primary and secondary infections of DENV, little progress has been made regarding secondary infection of DENV after a primary infection of ZIKV, or DENV-ZIKV co-infection. Here, we address this critical gap by developing compartmental models of virus dynamics. We first fitted the models to published data on dengue viral loads of the primary and secondary infections with the observation that the primary infection reaches its peak much more gradually than the secondary infection. We then quantitatively show that ADE is the key factor determining a sharp increase/decrease of viral load near the peak time in the secondary infection. In comparison, our simulations of DENV and ZIKV co-infection (simultaneous rather than sequential) show that ADE has very limited influence on the peak DENV viral load. This indicates pre-existing immunity to ZIKV is the determinant of a high level of ADE effect. Our numerical simulations show that (i) in the absence of ADE effect, a subsequent co-infection is beneficial to the second virus; and (ii) if ADE is feasible, then a subsequent co-infection can induce greater damage to the host with a higher peak viral load and a much earlier peak time for the second virus, and for the second peak for the first virus.

Infection severity across scales in multi-strain immuno-epidemiological Dengue model structured by host antibody level

Infection by distinct Dengue virus serotypes and host immunity are intricately linked. In particular, certain levels of cross-reactive antibodies in the host may actually enhance infection severity leading to Dengue hemorrhagic fever (DHF). The coupled immunological and epidemiological dynamics of Dengue calls for a multi-scale modeling approach. In this work, we formulate a within-host model which mechanistically recapitulates characteristics of antibody dependent enhancement in Dengue infection. The within-host scale is then linked to epidemiological spread by a vector-host partial differential equation model structured by host antibody level. The coupling allows for dynamic population-wide antibody levels to be tracked through primary and secondary infections by distinct Dengue strains, along with waning of cross-protective immunity after primary infection. Analysis of both the within-host and between-host systems are conducted. Stability results in the epidemic model are formulated via basic and invasion reproduction numbers as a function of immunological variables. Additionally, we develop numerical methods in order to simulate the multi-scale model and assess the influence of parameters on disease spread and DHF prevalence in the population.

Choices and trade-offs in inference with infectious disease models

Inference using mathematical models of infectious disease dynamics can be an invaluable tool for the interpretation and analysis of epidemiological data. However, researchers wishing to use this tool are faced with a choice of models and model types, simulation methods, inference methods and software packages. Given the multitude of options, it can be challenging to decide on the best approach. Here, we delineate the choices and trade-offs involved in deciding on an approach for inference, and discuss aspects that might inform this decision. We provide examples of inference with a dataset of influenza cases using the R packages pomp and rbi.

Antibody-Dependent Enhancement of Ebola Virus Infection by Human Antibodies Isolated from Survivors

Some monoclonal antibodies (mAbs) recovered from survivors of filovirus infections can protect against infection. It is currently unknown whether natural infection also induces some antibodies with the capacity for antibody-dependent enhancement (ADE). A panel of mAbs obtained from human survivors of filovirus infection caused by Ebola, Bundibugyo, or Marburg viruses was evaluated for their ability to facilitate ADE. ADE was observed readily with all mAbs examined at sub-neutralizing concentrations, and this effect was not restricted to mAbs with a particular epitope specificity, neutralizing capacity, or subclass. Blocking of specific Fcγ receptors reduced but did not abolish ADE that was associated with high-affinity binding antibodies, suggesting that lower-affinity interactions still cause ADE. Mutations of Fc fragments of an mAb that altered its interaction with Fc receptors rendered the antibody partially protective in vivo at a low dose, suggesting that ADE counteracts antibody-mediated protection and facilitates dissemination of filovirus infections.

A population of bang-bang switches of defective interfering particles makes within-host dynamics of dengue virus controllable

The titre of virus in a dengue patient and the duration of this viraemia has a profound effect on whether or not a mosquito will become infected when it feeds on the patient and this, in turn, is a key driver of the magnitude of a dengue outbreak. The assessment of the heterogeneity of viral dynamics in dengue-infected patients and its precise treatment are still uncertain. Infection onset, patient physiology and immune response are thought to play major roles in the development of the viral load. Research has explored the interference and spontaneous generation of defective virus particles, but have not examined both the antibody and defective particles during natural infection. We explore the intrinsic variability in the within-host dynamics of viraemias for a population of patients using the method of population of models (POMs). A dataset from 208 patients is used to initially calibrate 20,000 models for the infection kinetics for each of the four dengue virus serotypes. The calibrated POMs suggests that naturally generated defective particles may interfere with the viraemia, but the generated defective virus particles are not adequate to reduce high fever and viraemia duration. The effect of adding excess defective dengue virus interfering particles to patients as a therapeutic is evaluated using the calibrated POMs in a bang-bang (on-off or two-step) optimal control setting. Bang-bang control is a class of binary feedback control that turns either 'ON' or 'OFF' at different time points, determined by the system feedback. Here, the bang-bang control estimates the mathematically optimal dose and duration of the intervention for each model in the POM set.

Mathematical Analysis of Viral Replication Dynamics and Antiviral Treatment Strategies: From Basic Models to Age-Based Multi-Scale Modeling

Viral infectious diseases are a global health concern, as is evident by recent outbreaks of the middle east respiratory syndrome, Ebola virus disease, and re-emerging zika, dengue, and chikungunya fevers. Viral epidemics are a socio-economic burden that causes short- and long-term costs for disease diagnosis and treatment as well as a loss in productivity by absenteeism. These outbreaks and their socio-economic costs underline the necessity for a precise analysis of virus-host interactions, which would help to understand disease mechanisms and to develop therapeutic interventions. The combination of quantitative measurements and dynamic mathematical modeling has increased our understanding of the within-host infection dynamics and has led to important insights into viral pathogenesis, transmission, and disease progression. Furthermore, virus-host models helped to identify drug targets, to predict the treatment duration to achieve cure, and to reduce treatment costs. In this article, we review important achievements made by mathematical modeling of viral kinetics on the extracellular, intracellular, and multi-scale level for Human Immunodeficiency Virus, Hepatitis C Virus, Influenza A Virus, Ebola Virus, Dengue Virus, and Zika Virus. Herein, we focus on basic mathematical models on the population scale (so-called target cell-limited models), detailed models regarding the most important steps in the viral life cycle, and the combination of both. For this purpose, we review how mathematical modeling of viral dynamics helped to understand the virus-host interactions and disease progression or clearance. Additionally, we review different types and effects of therapeutic strategies and how mathematical modeling has been used to predict new treatment regimens.

Transmission-clearance trade-offs indicate that dengue virulence evolution depends on epidemiological context

An extensive body of theory addresses the topic of pathogen virulence evolution, yet few studies have empirically demonstrated the presence of fitness trade-offs that would select for intermediate virulence. Here we show the presence of transmission-clearance trade-offs in dengue virus using viremia measurements. By fitting a within-host model to these data, we further find that the interaction between dengue and the host immune response can account for the observed trade-offs. Finally, we consider dengue virulence evolution when selection acts on the virus’s production rate. By combining within-host model simulations with empirical findings on how host viral load affects human-to-mosquito transmission success, we show that the virus’s transmission potential is maximized at production rates associated with intermediate virulence and that the optimal production rate critically depends on dengue’s epidemiological context. These results indicate that long-term changes in dengue’s global distribution impact the invasion and spread of virulent dengue virus genotypes.

Theory predicts that pathogens will evolve towards intermediate virulence, yet the necessary trade-offs invoked by this theory have rarely been demonstrated empirically. Here, the authors show that dengue virus dynamics exhibit a trade-off between transmission and clearance rates.

Contributions from the silent majority dominate dengue virus transmission

Despite estimates that, each year, as many as 300 million dengue virus (DENV) infections result in either no perceptible symptoms (asymptomatic) or symptoms that are sufficiently mild to go undetected by surveillance systems (inapparent), it has been assumed that these infections contribute little to onward transmission. However, recent blood-feeding experiments with Aedes aegypti mosquitoes showed that people with asymptomatic and pre-symptomatic DENV infections are capable of infecting mosquitoes. To place those findings into context, we used models of within-host viral dynamics and human demographic projections to (1) quantify the net infectiousness of individuals across the spectrum of DENV infection severity and (2) estimate the fraction of transmission attributable to people with different severities of disease. Our results indicate that net infectiousness of people with asymptomatic infections is 80% (median) that of people with apparent or inapparent symptomatic infections (95% credible interval (CI): 0–146%). Due to their numerical prominence in the infectious reservoir, clinically inapparent infections in total could account for 84% (CI: 82–86%) of DENV transmission. Of infections that ultimately result in any level of symptoms, we estimate that 24% (95% CI: 0–79%) of onward transmission results from mosquitoes biting individuals during the pre-symptomatic phase of their infection. Only 1% (95% CI: 0.8–1.1%) of DENV transmission is attributable to people with clinically detected infections after they have developed symptoms. These findings emphasize the need to (1) reorient current practices for outbreak response to adoption of pre-emptive strategies that account for contributions of undetected infections and (2) apply methodologies that account for undetected infections in surveillance programs, when assessing intervention impact, and when modeling mosquito-borne virus transmission.

Most dengue virus infections result in either no perceptible symptoms or symptoms that are so mild that they go undetected by surveillance systems. It is unclear how much these infections contribute to the overall transmission and burden of dengue. At an individual level, we show that people with asymptomatic infections are approximately 80% as infectious to mosquitoes as their symptomatic counterparts. At a population level, we show that approximately 88% of infections result from people who display no apparent symptoms at the time of transmission. These results suggest that individuals undetected by surveillance systems may be the primary reservoir of dengue virus transmission and that policy for dengue control and prevention must be revised accordingly.

Challenges in dengue research: A computational perspective

The dengue virus is now the most widespread arbovirus affecting human populations, causing significant economic and social impact in South America and South-East Asia. Increasing urbanization and globalization, coupled with insufficient resources for control, misguided policies or lack of political will, and expansion of its mosquito vectors are some of the reasons why interventions have so far failed to curb this major public health problem. Computational approaches have elucidated on dengue's population dynamics with the aim to provide not only a better understanding of the evolution and epidemiology of the virus but also robust intervention strategies. It is clear, however, that these have been insufficient to address key aspects of dengue's biology, many of which will play a crucial role for the success of future control programmes, including vaccination. Within a multiscale perspective on this biological system, with the aim of linking evolutionary, ecological and epidemiological thinking, as well as to expand on classic modelling assumptions, we here propose, discuss and exemplify a few major computational avenues-real-time computational analysis of genetic data, phylodynamic modelling frameworks, within-host model frameworks and GPU-accelerated computing. We argue that these emerging approaches should offer valuable research opportunities over the coming years, as previously applied and demonstrated in the context of other pathogens.

The use of longitudinal cohorts for studies of dengue viral pathogenesis and protection

In this review, we describe how longitudinal prospective community-based, school-based, and household-based cohort studies contribute to improving our knowledge of viral disease, focusing specifically on contributions to understanding and preventing dengue. We describe how longitudinal cohorts enable measurement of essential disease parameters and risk factors; provide insights into biological correlates of protection and disease risk; enable rapid application of novel biological and statistical technologies; lead to development of new interventions and inform vaccine trial design; serve as sentinels in outbreak conditions and facilitate development of critical diagnostic assays; enable holistic studies on disease in the context of other infections, comorbidities, and environmental risk factors; and build research capacity that strengthens national and global public health response and disease surveillance.

Antibody-dependent enhancement of severe dengue disease in humans

For dengue viruses 1 to 4 (DENV1-4), a specific range of antibody titer has been shown to enhance viral replication in vitro and severe disease in animal models. Although suspected, such antibody-dependent enhancement of severe disease has not been shown to occur in humans. Using multiple statistical approaches to study a long-term pediatric cohort in Nicaragua, we show that risk of severe dengue disease is highest within a narrow range of preexisting anti-DENV antibody titers. By contrast, we observe protection from all symptomatic dengue disease at high antibody titers. Thus, immune correlates of severe dengue must be evaluated separately from correlates of protection against symptomatic disease. These results have implications for studies of dengue pathogenesis and for vaccine development, because enhancement, not just lack of protection, is of concern.

Modelling the effects of phylogeny and body size on within-host pathogen replication and immune response

Understanding how quickly pathogens replicate and how quickly the immune system responds is important for predicting the epidemic spread of emerging pathogens. Host body size, through its correlation with metabolic rates, is theoretically predicted to impact pathogen replication rates and immune system response rates. Here, we use mathematical models of viral time courses from multiple species of birds infected by a generalist pathogen (West Nile Virus; WNV) to test more thoroughly how disease progression and immune response depend on mass and host phylogeny. We use hierarchical Bayesian models coupled with nonlinear dynamical models of disease dynamics to incorporate the hierarchical nature of host phylogeny. Our analysis suggests an important role for both host phylogeny and species mass in determining factors important for viral spread such as the basic reproductive number, WNV production rate, peak viraemia in blood and competency of a host to infect mosquitoes. Our model is based on a principled analysis and gives a quantitative prediction for key epidemiological determinants and how they vary with species mass and phylogeny. This leads to new hypotheses about the mechanisms that cause certain taxonomic groups to have higher viraemia. For example, our models suggest that higher viral burst sizes cause corvids to have higher levels of viraemia and that the cellular rate of virus production is lower in larger species. We derive a metric of competency of a host to infect disease vectors and thereby sustain the disease between hosts. This suggests that smaller passerine species are highly competent at spreading the disease compared with larger non-passerine species. Our models lend mechanistic insight into why some species (smaller passerine species) are pathogen reservoirs and some (larger non-passerine species) are potentially dead-end hosts for WNV. Our techniques give insights into the role of body mass and host phylogeny in the spread of WNV and potentially other zoonotic diseases. The major contribution of this work is a computational framework for infectious disease modelling at the within-host level that leverages data from multiple species. This is likely to be of interest to modellers of infectious diseases that jump species barriers and infect multiple species. Our method can be used to computationally determine the competency of a host to infect mosquitoes that will sustain WNV and other zoonotic diseases. We find that smaller passerine species are more competent in spreading the disease than larger non-passerine species. This suggests the role of host phylogeny as an important determinant of within-host pathogen replication. Ultimately, we view our work as an important step in linking within-host viral dynamics models to between-host models that determine spread of infectious disease between different hosts.

Zika plasma viral dynamics in nonhuman primates provides insights into early infection and antiviral strategies

The recent outbreak of Zika virus (ZIKV) has been associated with fetal abnormalities and neurological complications, prompting global concern. Here we present a mathematical analysis of the within-host dynamics of plasma ZIKV burden in a nonhuman primate model, allowing for characterization of the growth and clearance of ZIKV within individual macaques. We estimate that the eclipse phase for ZIKV, the time between cell infection and viral production, is most likely short (∼4 h), the median within-host basic reproductive number R₀ is 10.7, the rate of viral production is rapid (>25,000 virions d^-1), and the lifetime of an infected cell while producing virus is ∼5 h. We also estimate that the minimum number of virions produced by an infected cell over its lifetime is ∼5,500. We assess the potential effect of an antiviral treatment that blocks viral replication, showing that the median time to undetectable plasma viral load (VL) can be reduced from ∼5 d to ∼3 d with a drug concentration ∼15 times the drug's EC₅₀ when treatment is given prophylactically starting at the time of infection. In the case of favipiravir, a polymerase inhibitor with activity against ZIKV, we predict a dose of 150 mg/kg given twice a day initiated at the time of infection can reduce the peak median VL by ∼3 logs and shorten the time to undetectable median VL by ∼2 d, whereas treatment given 2 d postinfection is mostly ineffective in accelerating plasma VL loss in macaques.

Drivers of Inter-individual Variation in Dengue Viral Load Dynamics

Dengue is a vector-borne viral disease of humans that endemically circulates in many tropical and subtropical regions worldwide. Infection with dengue can result in a range of disease outcomes. A considerable amount of research has sought to improve our understanding of this variation in disease outcomes and to identify predictors of severe disease. Contributing to this research, patterns of viral load in dengue infected patients have been quantified, with analyses indicating that peak viral load levels, rates of viral load decline, and time to peak viremia are useful predictors of severe disease. Here, we take a complementary approach to understanding patterns of clinical manifestation and inter-individual variation in viral load dynamics. Specifically, we statistically fit mathematical within-host models of dengue to individual-level viral load data to test virological and immunological hypotheses explaining inter-individual variation in dengue viral load. We choose between alternative models using model selection criteria to determine which hypotheses are best supported by the data. We first show that the cellular immune response plays an important role in regulating viral load in secondary dengue infections. We then provide statistical support for the process of antibody-dependent enhancement (but not original antigenic sin) in the development of severe disease in secondary dengue infections. Finally, we show statistical support for serotype-specific differences in viral infectivity rates, with infectivity rates of dengue serotypes 2 and 3 exceeding those of serotype 1. These results contribute to our understanding of dengue viral load patterns and their relationship to the development of severe dengue disease. They further have implications for understanding how dengue transmissibility may depend on the immune status of infected individuals and the identity of the infecting serotype.

Dengue is an important vector-borne disease that infects four-hundred million individuals annually. Infection results in a wide range of clinical symptoms. Though many risk factors of dengue are known, the mechanisms explaining why an individual will suffer severe symptoms are poorly understood. Clinical studies have shown characteristics of viral load kinetics of dengue-infected individuals may be indicators of disease severity. However, viral load measurements vary considerably by individual. Here we use statistical methods to empirically test hypotheses that may explain variation in dengue viral load patterns by clinical manifestation and by serotype. We show that there is statistical support for antibodies being responsible for higher disease severity during secondary dengue infections and for high viral infectivity rates of dengue serotypes 2 and 3 relative to dengue 1. These results further understanding of the relationship between viral load patterns and severe dengue disease and have important implications for dengue transmissibility.