Next Article in Journal
Phage Therapy in Germany—Update 2023
Next Article in Special Issue
Use of Oral Antivirals Ritonavir-Nirmatrelvir and Molnupiravir in Patients with Multiple Myeloma Is Associated with Low Rates of Severe COVID-19: A Single-Center, Prospective Study
Previous Article in Journal
Bioinformatic Tools for NGS-Based Metagenomics to Improve the Clinical Diagnosis of Emerging, Re-Emerging and New Viruses
Previous Article in Special Issue
Understanding the SARS-CoV-2 Virus Neutralizing Antibody Response: Lessons to Be Learned from HIV and Respiratory Syncytial Virus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Department of Life Science, Dezhou University, Dezhou 253023, China
2
School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200030, China
3
Laboratory AGEIS EA 7407, Team Tools for e-Gnosis Medical, Faculty of Medicine, University Grenoble Alpes (UGA), 38700 La Tronche, France
*
Authors to whom correspondence should be addressed.
Viruses 2023, 15(2), 586; https://doi.org/10.3390/v15020586
Submission received: 30 January 2023 / Revised: 10 February 2023 / Accepted: 17 February 2023 / Published: 20 February 2023
(This article belongs to the Special Issue RNA Viruses and Antibody Response)

Abstract

:
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.

1. Introduction

The COVID-19 epidemic has caused more than 630 million infections and over 6.5 million deaths worldwide by the end of 2022, and it can hardly disappear in a short time based on observation [1,2,3]. Until now, vaccination has been the only approach to fighting SARS-CoV-2 infections. As the global promotion of the SARS-CoV-2 vaccine continues, we are gaining a better understanding of the antibodies triggered by vaccines and natural infections. It is gradually becoming apparent that antibodies against SARS-CoV-2, whether acquired from natural infection or triggered by vaccination, decay over time, with a concomitant decrease in protective efficiency. Unlike antibodies from vaccinations such as smallpox [4], antibodies against SARS-CoV-2 do not provide durable protection [5].
A significant portion of the population is reluctant to be vaccinated, according to a public poll survey [6,7]. The first concern comes from the side effects of vaccination, such as blood clots [8,9]. The second concern is this deficiency of the long-lasting protection effect brought by vaccination since many clinical reports indicate that the vaccine’s protection effect is declining over time [10,11,12,13]. Based on the clinical data, a significant decrease in the immune effect of the vaccine against reinfection has been noticed. The SARS-CoV-2 vaccine still shows an overall positive effect on protection currently. However, more people are concerned about how long the protection can be provided by the antibodies triggered by their immune response [5,14]. There is an urgent need to evaluate the protection time of neutralizing antibodies quantitatively.
Constructing a reliable mathematical model helps us quantitatively explain the virus dynamics in the host body, which could provide a reasonable prediction toward those sensitive concerns faced by the public. Researchers have conducted many studies on virus infections using computational approaches [15,16].
Specifically, depending on whether the immune response is integrated into the models, traditional cellular models can be divided into two categories. The first set contains models that do not include the immune response [17,18,19,20,21], that is, traditional viral kinetic models, whose following equation can be classically expressed as:
d T d t = s T d T β T V
d I d t = β T V δ I
d V d t = p I c V
where T denotes the target cells or the number of susceptible cells; I denotes the number of cells where the infection has occurred; and V denotes the number of viruses. This simplified model is wildly used to study the virus dynamics upon infection. While admitting its advantages in recapturing the viral dynamics in virus infection, concerns about this model are proposed. The limitation of this model originates from the absence of a connection to the body’s immune response. The main driving force behind eliminating the virus is activating the adaptive immune response by lymphocyte B and T cells, not the depletion of the susceptible cells. Therefore, the estimation based on those depletion mechanisms is skeptical, as it would return a much smaller initial susceptible cell number. The number of susceptible cells should be higher than the number estimated by several orders. Mandating the estimation would inevitably lead to underestimating T0, where T0 represents the number of initial susceptible cells. The T0 term is similar to the initial number of susceptible people when the SIR (susceptible, infection, recovery) model is used to predict the epidemic trend [22,23]. The estimated T0 values vary greatly due to the differences in the fitting data, which significantly reduces their reliability.
The second set consists of models that consider the immune response, where the classical expression of the model considering antibody binding is displayed as follows [24,25,26,27,28,29,30]:
d T d t = χ D β T V
d I 1 d t = β T V k I 1
d F d t = ω V α F
d D d t = δ I 2 χ D
d I 2 d t = k I 1 δ I 2
d A d t = f V h A  
d V d t = p 1 + ε 1 F   I 2 c V γ   T V κ A V
where T denotes the overall number of susceptible cells; A denotes the number of antibodies; V denotes the amount of virus; D denotes the number of dead cells; F denotes the amount of interferon in nonspecific immunity; I 1 denotes the number of cells infected in class 1 (cells infected without viral multiplication); and I 2 denotes the number of cells infected in class 2 (cells infected with viral multiplication). This model effectively combines antibody dynamics with viral changes. It performs well in recapturing and explaining the dynamic interaction between the virus and the host immune system. However, the kinetic properties of antibodies have not been explicitly integrated into this model. Different antigen-binding antibodies will exhibit different kinetics, and this diversity in binding behavior is an essential factor influencing the host-pathogen interaction. The difference in the antibody binding capacity would lead to a significant variation in the virus dynamics. However, this classical model does not explicitly formulize the antibody term.
Although mathematical modeling of antibody kinetics has been significantly advanced and developed in recent years, there is still great potential to be further improved. Therefore, a new mathematical model of antibody dynamics is proposed by us. This model displays good fitness in explaining some puzzling phenomena, especially in the case of SARS-CoV-2 infection, with improvements to previous models in several aspects.
First: New immunology responses are formulized and integrated into this model based on biochemistry principles. We explicitly formulize the virtual effect of antibodies on eliminating pathogenic microorganisms while noticing the stimulation effect of pathogenic antigens on antibodies [31]. This positive feedback stimulation is explicitly represented in a mathematical function as well.
Second: The dynamic model developed hypothesizes that environmental factors maintain memory cell levels. The natural attenuation of antibody levels can be recaptured in this model. Furthermore, their discontinuous decay trend can also be captured with the introduction of environmental antigens. Antigen-like substances in the environment can maintain a specific concentration of B cells or T cells at a certain level. This can explain why some vaccines can provide lifelong protection.
Third: The antigen–antibody binding process is represented as reversible in our model. Instead of using the equilibrium constant, the binding dynamic is described in both the binding and reverse dissociation reactions.
Based on the basic principles of immunology, we established the theoretical hypothesis of antibody kinetics. The results and methods sections illustrate specific principles and rationales in detail. In the end, this antibody kinetic model proposes some possible mechanisms underlying some real-world scenarios:
  • How are memory cells maintained?
  • How does our immune system screen for antibodies with a strong binding affinity?
  • Why do people who get influenza and other vaccines have a lower mortality rate from SARS-CoV-2?
  • How can we effectively calculate the duration of protection of a specific antibody?
  • Why are some recovered patients retested as positive cases without infections from other people?
  • Why do vaccinations show considerable differences in protection efficiency?
  • How can we improve the protective efficiency and duration of vaccines?

2. Materials and Methods

2.1. Mathematical Representation of the Antibody Production Process

A simple mathematical representation of the immune response is described in the diagram in Figure 1 below.
Here, x denotes the number of antibody–antigen (virus) complexes, y denotes the total number of antibodies, and z denotes the number of viruses. Six processes are displayed in our model. The first reaction represents the virus’s proliferation or replication with a reaction constant k 1 . The second reaction represents the binding reaction between virus and antibody, with a forward reaction constant k 2 and reverse constant k 2 . The third reaction represents removing the antibody–virus complex with a reaction constant k 3 with the help of natural killer (NK) cells [32]. The fourth reaction represents the induction of a new antibody by the antibody–virus complex with a kinetic constant k 4 . In immunology, those virus–antibody complexes are on the surface of B-cells since the antibodies are initially produced by B-cells and will attach to the plasma membrane of B-cells. Those complexes would further bind to the helper cells because the antibody has another structure binding region toward those receptors. Those helper cells will present the antigen part, a virus in this case, to the T-cells. The physical placement should be such that B-cells bind to those helper cells and further present themselves close to T-cells. The T-cells will handle those antigen substances; if those substances are not self-originated, they will secrete signal molecules to promote the proliferation or division of B-cells that attach to them. Therefore, B-cells finally proliferate, along with the antibodies their B-cells generate. The fifth reaction represents the degradation of the virus with a constant k 5 . The sixth reaction represents the degradation of antibodies with a rate constant k 6 .
The differential equations based on those reactions are derived as below:
We established the following equations to represent the proliferation process of antibodies:
d x d t = k 2 y z k 2 x k 3 x ,
d y d t = k 2 x k 2 y z + k 4 x k 6 y ,
d z d t = k 2 x k 2 y z k 5 z + k 1 z .

2.2. Mathematical Modeling including Environmental Antigens

The equations introduced in Section 2.1 describe antibody dynamics in the presence of the associated virus. In reality, other environmental antigens can induce antibody coupling and production. To further account for this possibility, we add a new set of equations as shown in Figure 2 below:
Where p denotes antigen-like substances in the environment that are supposed to be constant; q denotes antibodies bound to antigen-like substances in the environment. Reaction 7 represents the binding reaction between antibodies and environmental antigenic substances with a forward constant k 7   and a reverse constant k 7 . Reaction 8 represents the removal of the antibody–antigen complex q with a reaction rate k 3 . Reaction 9 represents the induction of a new antibody by q. Therefore, a new set of equations is derived as follows:
d x d t = k 2 y z k 2 x k 3 x ,
d y d t = k 2 x k 2 y z + k 4 x k 6 y k 7 p y + k 7 q + k 4 q ,
d z d t = k 2 x k 2 y z k 5 z + k 1 z ,
d p d t = 0 ,
d q d t = k 7 p y k 7 q k 3 q .
The antibody level will eventually drop to zero based on the first model described in Equations (11)–(13). Since the decaying coefficient k 6 is more significant than zero, the antibody will finally fade away quickly. However, as we all know, some antibodies can persist in the human body for a long time and provide lifelong protection. It forms the basis for a vaccine. Therefore, “environmental antigen-like substances” are introduced into the second model. They could be food-resourced, air-resourced, or even self-resourced. The presentation of those substances to T-cells would give weak signals to proliferate the memory B-cells or T-cells. The detailed physical background is introduced in Section 3.1.

2.3. A Simplified Model Simulating the Proliferation of Antibodies with Different Binding Kinetic Characteristics by the Immune System

d x i d t = K i on y i z K i off x i k 3 x i ,
d y i d t = K i o f f x i K i on y i z + k 4 x i k 6 ( y i w i ) ,
d z d t = i = 1 n ( K i off x i K i on y i z ) k 5 z + k 1 z .
Here, x i denotes the amount of the i-th antibody that binds to the viral antigen, y i denotes the total number of the i-th antibody, w i is a constant that indicates the threshold for maintaining the i-th antibody by the antigen-like substances in the environment. K i on denotes the binding constant of the free antibody i and the viral antigen. K i off   denotes the reaction coefficient for the dissociation reaction of the antibody i viral complex. After we finish the second model, we want to show how our immune system selects good-binding antibodies. Therefore, we introduce the third model. The third model is a combination of the first model and second model. The difference between the third and second models is that they do not explicitly integrate environmental factors. Since each kind of antibody will have a corresponding environmental counterpart, it would be too complicated if we explicitly integrated all of them. Therefore, we use a simplified term w i to represent the minimal antibody level threshold maintained by its corresponding environmental-antigenic substances. Each antibody has different binding kinetics toward the virus. The modeling result demonstrates the dynamic process by which our immune system selects those fast-binding antibodies over slower ones.
The above three models are numerically solved by MATLAB using the ode15s function [33].

3. Results

3.1. Physical Mechanism behind This Approach and the Underlying Relationships among the Three Models

Our antibody model is based on the following four rationales.
The first is the stimulus of the antibody–antigen complex on the regeneration of its corresponding antibody. The T-cells will induce the proliferation of the B-cells through the antigen-presenting process. Eventually, the corresponding antibody can be reproduced following the magnification of antibody-generating B-cells [34]. Although this regulation is critical in initializing the adaptive immune response, it was not explicitly represented in previous studies [24,25,26,27,28,29,30]. The k 4 x term in Equation (12) expresses this positive feedback regulation. Integrating feedback regulation into the classical representation of biological reactions allows us to simulate reactions more comprehensively. Specifically, when we refer to the concentration alternation of antibody and virus, the change of antibody in Equation (12) should be   k 2 x k 2 y z k 6 y , if we use chemical reactions to describe it. However, in this way, we ignore the virus–antibody complex’s stimulation of the regeneration of antibodies. In immunology, those virus–antibody complexes will proliferate the corresponding antibody with the help of T-cells after an antigen presentation [35]. Therefore, we use an external k 4 x term to model this positive feedback effect.
The second is the kinetic relationship between antibodies and antigens. This antigen–antibody binding is represented in the second process in the first diagram. The model is endowed with more dynamic characteristics with the introduction of a reversible reaction.
The third is the pernicious effect of antibodies on antigens, which is expressed as the k 3 x term in Equation (11). It has been recognized that the antibody–antigen complex can be efficiently removed with the help of functional immune cells such as NK (natural killer) cells [32].
The fourth is the introduction of environmental antigenic substances. We hypothesize that the environmental antigen contributes to the maintenance of memory cells. The declination of induced neutralizing antibodies is ubiquitous in almost all virus infections, such as Zika [36], Dengue [29], and SARS-CoV-2 [37]. However, this declination is always discontinuous after their concentration drops to a stable level. This is contributed by the existence of memory cells [38]. Immunologists regarded the long-lasting B-cells or T-cells as “memory cells.” Experiments gradually demonstrated that although the so-called “memory cells” are some specific forms of immune cells [39,40], they have similar half-lives as normal CD8+ cells [41]. Therefore, maintaining antibody levels in “memory cells” should be explained as a state of equilibrium between decay and regeneration. We suppose this final equilibrium state derives from environmental antigen-like substances. Due to the presence of cross-interaction [42], the neutralizing antibody, no matter how specific it is, could have weak binding with other macromolecules in the solution. Those weak binding partners are defined as “environmental antigen-like substances” in our model; they could be food-resourced, air-resourced, or self-resourced. The presentation of those substances to T cells would give weak proliferation signals to neighboring cells. Therefore, a lower bound for specific neutralizing antibodies is present due to the positive feedback regulation of environmental antigenic substances. This hypothesis could explain why some antibodies persist in the body at detectable levels for a long time. We propose that the presence of environmental antigenic substances allows us to produce lifelong immunity to some pathogens.
The antibody level will eventually drop to zero based on the first model described in Equations (11)–(13). This implies that the antibody would finally fade in a short time. However, some antibodies can persist in the human body for a long time and provide lifelong protection, which is why we use vaccines. Following the last rationale, a second model is constructed to integrate environmental antigenic substances. The protein sequence of environmental antigen-like substances is ordinarily close to our own body and less antigenic, leading to a weak T-cell proliferation signal. An equilibrium state can be reached when the decay of memory cells is equal to their new synthesis. A good fit between the two scenarios further improves the reliability of this model. The first is that those environmental antigen-like molecules would help maintain the antibody at a certain level. This scenario can be easily recaptured in our second model. Our model can also explain the second scenario: environmental antigen-like substances cannot naturally elevate the antibody level to the final equilibrium concentration. Environmental antigenic substances can only help maintain the antibody level instead of triggering the increment. An exception is the allergy reaction, which happens under extreme conditions with the massive presence of environmental antigenic substances, especially when they have a strong binding affinity. Our second model can also recapture this phenomenon by increasing p and k 7 values.
Although the absolute equilibrium of specific antibodies is challenging to simulate mathematically, the presence of environmental antigen-like substances dramatically attenuates the antibody decay rate. The simulation results are shown in Figure 3. The decay rate of antibodies is closely related to the concentration of the antigen-like substance in the environment. The high concentration of the environmental antigenic substance would lead to a slow decay speed. The antibody decay rate is significantly slower (shown in the solid yellow curve) when there are massive environmental antigen-like substances, as shown in Figure 3A. The antibody level attenuation is also significantly influenced by the binding kinetics between environmental antigen and its corresponding antibody. The antibody decay rate is significantly slower (shown in the dashed yellow line) when k 3 sets a considerable value in Figure 3B. It represents a better binding kinetics between environmental antigen-like substance and its corresponding antibody.
It is common sense that environmental antigenic substances alone can hardly trigger antibody proliferation without pathogenic antigens. Numerical simulations also indicate that antigen-like substances in the environment hardly stimulate antibody proliferation. One scenario in which environmental antigen-like substances do not trigger antibody growth is recaptured in Figure 4A. The insufficient capacity of environmental antigens to induce antibody magnification derives from their poor binding affinity toward antibodies. In most cases, the presence of antigen-like substances in the environment will not directly stimulate antibody proliferation but will significantly attenuate the decay rate after antibody proliferation, as shown in Figure 4A. The situation can be altered when the antigen-like substances are sufficient or their binding affinities toward some antibodies are powerful. Antigen-like substances in the environment can also stimulate antibody levels, as shown in Figure 4B. This might help elucidate the underlying mechanism of allergy. A significant difference between environmental antigen-like substances and pathogenic antigens is that the former are deficient in self-replication. However, their concentrations are maintained at a relatively stable level due to constant replenishment from the environment.

3.2. Characteristics of Immune Response after Infected with Different Virus Strains

For the same type of virus, highly virulent strains have high replication activity [43]. Infection with highly virulent variants will lead to a prominent peak viral load and can exhibit strong transmissibility per time unit. However, its infection cycle is typically short. These infection dynamics are defined as acute infections in clinical performance [44]. The cytokine storm induced by a quickly elevated antibody level is responsible for the acute response [45]. Weak virulent strains have the opposite behaviors, which could lead to chronic infection in clinical settings. It does not typically induce a quick immune response, resulting in a more extended infection cycle. All those situations can be recaptured using our model, which is displayed in Figure 5.
It can be noticed from Figure 5 that, under the exact initial viral invasion dosage, the infection caused by a less virulent strain has a more extended latency phase. Its peak viral load is significantly lower than a more virulent strain. The maximal level of the yellow line is significantly lower than that of the dashed blue line, indicating that the antibody production induced by the weak virus infection is significantly lower than that induced by a virulent strain. It is also reflected in Figure 5 that the infection cycle caused by the weak virus is remarkably longer than its strong counterpart.

3.3. How the Immune System Screens for Highly Binding Antibodies

Our immune system can automatically detect suitable binding antibodies and selectively proliferate them. Simulating the antibody screening process by mathematical modeling remains challenging for computational biology researchers. This study provides a preliminary attempt to elucidate this process. Our model suggests a strong correlation between the proliferation potential of a specific antibody and its binding kinetics. Antibodies with fast binding rates would obtain rapid proliferation, while those with low binding rates proliferate slower. Thus, the antibodies with fast binding rates can grow faster and dominate in the final antibodies’ composition, although various antibodies are triggered to proliferate during infection. Generally, antibodies with fast binding capacity have a more robust binding affinity. However, this relationship is not absolute. Binding affinity is the ratio of forward binding to reverse dissociation [46]. Infected individuals would have antibodies with good binding kinetics after infection [47,48]. Those antibodies display individual differences in binding affinity, with some antibodies having relatively stronger binding affinity [49]. We propose that the dominant factor affecting the efficiency of antibodies is not their absolute binding affinity K d but the binding reaction constant k 2 . Both   k 2 and k 2 have physical meaning in Equation (12): k 2 stands for the forward reaction constant, while k 2 represents the reaction constant of the reverse one. The dissociation constant K d = k 2 / k 2 can be used to evaluate the concentration of each component in an equilibrium state. K d can also be transformed into the binding energy between those two macromolecules and thus is an important indicator that describes the binding affinity between two molecules [50].
As shown in Figure 6, the antibody with the fastest binding rate (solid blue line in this figure) demonstrates a maximal proliferation magnitude. This selection will eventually lead to a high proportion of fast-binding antibodies in the final antibody composition. Although these five antibodies have the same binding affinity of K d = 1 × 10−9, their proliferations vary because of the differences in their binding speeds, and the antibody with the faster binding rate will gradually dominate. This example illustrates that the driving force in antibody reproduction is its binding speed, not its binding affinity. A different k 2 value and k 2 value is applied to keep the dissociation constant K d as a fixed number. It is subversive to claim that the immune system would like to select those fast-binding antibodies rather than tightly binding antibodies. Experimental researchers always seek an antibody with an excellent binding affinity [47,48,49]. However, our model indicates that a solid-binding antibody might not perform well in preventing reinfection. Its kinetic behavior is more critical than its static binding affinity.

3.4. High Concentrations of Weakly Binding Antibodies Can Provide Effective Protection

It has been statistically discovered that vaccination with other vaccines, such as the influenza vaccine, can also provide some degree of protection against SARS-CoV-2 [51,52,53,54]. The third model is utilized to explain this phenomenon. Those parameters used in Figure 7 are defined in Section 2.3.
Figure 7 illustrates that the presence of weakly binding antibodies can produce an inhibitory capacity against the virus. The concentration of the weakly binding antibody increases from one in case 2 to a thousand in case 1. The maximal value of the dashed yellow line is smaller than that of the solid blue line, indicating a more substantial inhibition effect on virus reproduction. It can be noticed that the maximal level of the solid-binding antibody also decreases with the enhanced initial level of the weakly binding antibody. A particular concentration of weakly binding antibodies can simultaneously inhibit the proliferation of strongly binding, efficient antibodies in the host.
This example suggests that vaccinated people will always have better protection against reinfection than non-vaccinated people, regardless of the mutation of the virus. However, when we aim to stimulate the production of high levels of firmly bound antibodies through vaccination, the presence of weakly bound antibodies might interfere with this process. The vaccination effect on each individual might vary due to the weak antibody inference. Some individuals can induce sufficient fast-binding antibodies, while others cannot, due to differences in their initial antibody library. The simulation results indicate antibody interferences might be greatly enhanced in people vaccinated against other viruses. This explains why some vaccines may occasionally lead to a higher overall mortality rate [14,55]. Keeping a relatively high antibody level for a specific pathogen, such as SARS-CoV-2, will shelter more people from infection. However, it might also have the side effect of decreasing the efficiency of other vaccines. The long-term effect of repeated booster vaccinations on SARS-CoV-2 needs a comprehensive study.

3.5. Calculation of the Protection Time Brought by Natural Infection

The prediction of protection time remains to be elucidated using a computational approach. The calculation of protection time is presented in this section with the application of a numerical approach. Figure 8 illustrates how we calculate the protection time upon natural infection. Details about how we calculate the protection time are described below:
First, we model the dynamic behaviors of all components (virus, virus–antibody complex, antibody). An initial number is assigned to each component before the first infection. For instance, z(0) = 1 stands for one invading virus. The virus number would vanish to zero soon after the antibody booms. The virus number dropped to zero after the 166th simulation time point, as illustrated in Figure 8.
After the antibody completely removes the virus, the second invasion is performed. The initial virus number would be manipulated to a non-zero integer at each following time unit, starting from 166 to 800. The early invasion would not lead to virus proliferation because the antibody level at that time was still high. However, the protection effect will gradually fade as the antibody level decreases. The earliest invasion that could lead to a virus proliferation is at the 326th time point in Figure 8. It is a time point when a breakthrough infection can happen. The time between the first infection and this time point is the protection duration brought by the first infection. Therefore, the protection time can be roughly estimated from 166th to 326th. We can also notice that the maximal viral load in the second infection is significantly lower than in the first infection. This indicates that the reinfection or breakthrough infection typically has a milder symptom than the first, although it will occur with a higher chance as time goes on. The protective effect will decrease as time increases, as shown in Figure 9 in the next section.

3.6. Factors Affecting the Duration of Antibody Protection: Concentration of the Environmental Antigen-like Substance, Viral Replication Capacity, and Antibody Binding Kinetics

Based on Equation (16), the virus will magnify itself when the k 2 x k 2 y z k 5 z + k 1 z , term is larger than zero. The analytic solution on the reinfection possibility is difficult to derive, although the antibody level and the virus replication capacity directly influence this term. A high concentration of environmental antigenic substances p and an excellent binding kinetic k 7 can help maintain the antibody at a relatively high level, thus providing an extended protection time. Meanwhile, the antibody kinetic features could also greatly influence the protection time, which is shown in Figure 9. A lifelong protective effect can happen when a virus has a weak replication capacity (small k 1 value), a high concentration of environmental antigen-like substances (large p-value), and a fast-binding antibody with an enormous k 2 value. Typical examples of this category are the smallpox virus [4], tetanus bacillus [56], and so on. In contrast, other viruses, such as influenza and SARS-CoV-2, may not be able to trigger a super-long protection time [5]. Meanwhile, for influenza and SARS-CoV-2 infections, there is a considerable variation in the protective time after natural infection or vaccination [57,58]. This variation is mainly attributed to the difference in their corresponding antibody kinetics.
As claimed in the previous section, the binding kinetic constant k 2 is much more important than the binding affinity K d . Figure 9 illustrates how k 2 influences the protection time in different infection cases.
The second infection time is marked in Figure 9A, which starts around the 275th time unit. The protection time can be roughly estimated to be between 125th and 275th time unit. The protection time is significantly extended (100th to 1080th time unit) in Figure 9C when a better binding antibody is introduced. A virus loading threshold for the severe case is arbitrarily set as 6.4 × 105 marked by the red arrow in Figure 9B. The protection time against severe infection can be derived based on Figure 9B,D, which is the time interval between those two green lines. Figure 9D demonstrates that the severe infection would happen after the 975th time unit with the same severe threshold. The protection time can be roughly estimated to range from 125th to 975th time unit. This protection phase is much longer compared to Figure 9B. This example demonstrates that an excellent binding antibody will protect people longer than the weak one. It indicates that the protection time brought by infection or vaccination could significantly vary for different individuals. The primary reason lies in the differences in the binding kinetics of antibodies–virus interaction. Another important discovery in this modeling is that people who generate suitable antibodies with fast binding kinetics rarely display significant symptoms in the first infection, which can be reflected in the smaller virus peak loading amount in the first infection in Figure 9C compared to Figure 9A. It is suggested that those fast-binding antibodies typically exist in patients without intense infection symptoms. Repeated vaccination might be limited in reshaping the composition of the final antibody reservoir, i.e., it does not significantly alter the binding activity of the antibodies. Therefore, targeted improvement of antibody binding activity against certain viruses is the key to extending vaccine protection in the long term. The heterogeneity of the individual antibody library might cause significant differences in the binding activity of antibodies to the same viral infection or vaccination, leading to a significant variation in protection duration. Therefore, we may need to change the vaccination strategy. Vaccination should induce the production of antibodies with high binding activity instead of a random incitation. The targeted induction of specific antibodies using gene editing provides a potential solution [59,60,61].

3.7. Parameter Estimation in Real Scenario

According to the data from clinical experiments, we can further fit the specific parameters to obtain the dynamic characteristics of antibodies within certain populations. The core content of this model is to accurately simulate the trend of mutual change between antibodies and viruses, so it is not required to add units to the mathematical model. By adjusting different parameters, especially the mean, and variance of the parameters, we can get the characteristics of the antibody dynamics behavior of a population. By comparing the available statistical data, we can estimate the distribution of antibody dynamics parameters of the whole population. Although these parameters have no direct physical meaning, they can more accurately reflect the decline of immunity in different individuals and the possibility of reinfection at different times.
We treat the concentration of the environmental antigenic substance as a constant value. The primary alternation is the difference in the binding capacity of the antibody in different individuals, which is defined as k 2 in our model. The binding kinetic between the environmental antigenic substance and the antibody, defined as k 7 in our model, varies among different people. The antibody kinetic performance characteristics of the population obtained using the above parameter combination are shown in Figure 10A. It was demonstrated that IgG antibody levels to the SARS-CoV-2 nucleocapsid waned within months, according to six months of data from a longitudinal seroprevalence study of 3217 UK healthcare workers [62]. The simulation result has a good match with this clinical report. The main factor determining the peak antibody concentration is k 2 , while the driving force determining the antibody waning speed comes from k 7 . The change in the overall antibody level of the population over time is shown in Figure 10A, and the change in the overall protective efficacy is shown in Figure 10B. According to this model, we can judge that in the absence of virus mutation, the overall protective efficacy of an initial 100% efficacy of the COVID-19 vaccine in a human population (10,000 human simulations) would drop to 97.21% after 100 days, 65.44% after 150 days, 39.28% after 200 days, and 28% after 240 days. The simulation results in Figure 8B are consistent with the clinical data on vaccine effectiveness. Based on a cohort study among US veterans, for the period 1 February 2021 to 1 October 2021, vaccine effectiveness against infection (VE-I) declined over time (p < 0.01 for time dependence), even after adjusting for age, sex, and comorbidity. VE-I declined for all vaccine types, with the most significant declines for Janssen, followed by Pfizer-BioNTech and Moderna. Specifically, in March, VE-I was 86.4% for Janssen, 89.2% for Moderna, and 86.9% for Pfizer-BioNTech. By September, VE-I had declined to 13.1% for Janssen, 58.0% for Moderna, and 43.3% for Pfizer-BioNTech [63]. According to a retrospective cohort study on the effectiveness of the mRNA BNT162b2 vaccine, effectiveness against infections declined from 88% (95% CI 86–89) during the first month after complete vaccination to 47% after five months. Among sequenced infections, vaccine effectiveness against infections of the delta variant was high during the first month after complete vaccination (93%) but declined to 53% after four months [64]. Similar trends were observed in the cohort study conducted in Qatar [65].

3.8. Recovered Patients with Retest Positive for SARS-CoV-2

In SARS-CoV-2 infections, we noticed some rare occurrences of reinfections by patients themselves [66]. This recurrent positivity is not caused by the invasion of environmental viruses but by the re-proliferation of internal viruses [67]. The self-reinfection may lead to a second epidemic outbreak, making it more challenging to prevent and control the epidemic. The scenario of self-reinfection is a frequent phenomenon in virus infection. Many people suffer from recurrent respiratory infections [68]. It is also typical when HBV or HCV infection occurs [69]. Our model suggests a possible mechanism underlying this phenomenon. The patient will experience reinfection under certain parameter sets, as shown in Figure 11.
This example above shows that it is difficult to completely eliminate a specific pathogen in the presence of fast-binding antibodies. However, the peak viral load during the infection is small, as shown in the solid red line in Figure 11A. This will lead to a relatively mild symptom. The peak viral load in Figure 11A is significantly lower than the typical infection cases illustrated in Figure 8 and Figure 9.
The low virus load does not have an intense stimulation on antibody production. Therefore, it is impossible to eliminate the virus from the body due to insufficient antibody quantity. The virus will regain the opportunity to proliferate when the antibody level decreases later. The patient may display a positive nucleic acid test result again. Multiple virus resurgences are marked in the dashed green cycle after the first infection, as displayed in Figure 11A. However, the infections in those cases are generally less symptomatic and even asymptomatic. A more extreme case is a long-term chronic infection, such as the appearance of a long-positive patient [70]. Like self-reinfection, an equilibrium state can be reached if the antibody–antigen interaction is moderate in those long-positive patients. In this case, pathogens would not be eliminated but maintained at a relatively stable level, forming a chronic infection. The low concentration of pathogenic antigens only provides a limited driving force for promoting antibody reproduction. The antibody and the virus will remain relatively low for a long time, as reflected in Figure 11B.
Pathogenic antigens and environmental antigenic substances all contribute to chronic inflammation. For chronic infections caused by pathogenic microorganisms, a short-term boost in antibody levels can be used to accomplish a complete clearance of pathogenic microorganisms [71]. Chronic infections could permanently disappear or significantly improve after they become acute infections in the clinic. The chronic symptom will disappear or become invisible after the healing from the acute infection [72]. Conversely, chronic inflammation caused by environmental factors can be removed by shielding environmental antigens for a certain period [73]. The blockade of antibody–antigen interaction will decrease the corresponding antibody level, thereby attenuating or eliminating the immune response. An allergic reaction can be significantly alleviated in the upcoming contact after this treatment.

4. Discussion

The application of mathematical modeling displays good epidemic forecast capacity at the population level [22,23,74,75]. Mathematical models are also helpful in quantitatively elucidating the immune process’s dynamics on an individual level. We are motivated to develop a new mathematical model that can explicitly model the antibody dynamics based on those pioneering modeling attempts [17,18,19,20,21,24,25,26,27,28,29,30]. A computational antibody dynamic model is finally proposed, which can help us explain some phenomena mentioned in the introduction.
I.
How are memory cells maintained?
The environmental antigen-like substances maintain memory cells in our model. It was demonstrated that, in the absence of a pathogen, the environmental antigen could not directly stimulate antibody proliferation in most cases except in allergic reactions. The concentration of memory cells is closely correlated to the property and concentration of its corresponding environmental antigenic substances. Some memory cells against a specific pathogen can be maintained at a high level for decades. The persistent stimulation by the environmental antigenic substances, rather than the eternal lifespan, leads to the long-term existence of those immune memory cells. The existence of super-binding antibodies and high concentrations of environmental antigens can also trigger allergic reactions.
II.
How does our immune system screen for antibodies with solid binding affinity?
Based on the model, antibodies with faster binding kinetics would increase proliferation. The binding constant is not linearly correlated to the binding affinity. The immune system tends to select antibodies with fast binding kinetics rather than solid binding affinity. The underlying mechanism comes from the positive feedback regulation of the virus–antibody complex. A good binding antibody would lead to a faster generation of the virus–antibody complex. This complex would further promote the proliferation of its corresponding antibody.
III.
Why do people vaccinated by the influenza vaccine or other vaccines have a lower mortality rate from SARS-CoV-2 infection?
An interesting phenomenon during the COVID-19 pandemic is that the mortality rate is significantly lower in people who have received the influenza vaccine and other vaccines than in their unvaccinated counterparts [51,52,53,54]. The non-specific binding antibody level can be elevated after non-specific vaccination. Although weakly binding antibodies inhibit the proliferation of firmly binding antibodies in vivo, they can significantly inhibit virus proliferation and reduce peak viral load. This would lead to a lower mortality rate and milder symptoms. The potential concerns are also discussed in Section 3.4. A particular concentration of weakly binding antibodies can inhibit the proliferation of strongly binding antibodies, as revealed in Figure 7. The rationality of repeated booster vaccination needs a comprehensive evaluation because it might weaken the immune response to other pathogens.
IV.
How could we effectively calculate the protection duration of a specific antibody?
Concerns are remarkably booming when people realize the protection of the SARS-CoV-2 vaccine has declined over time. Given the experimental data on viruses and antibodies, the kinetic parameters of this model can be inferred. The protection time of the vaccine or natural infection can be further deduced. Personalized prediction is also feasible with the availability of individual antibody behavior. The protection times of different individuals can be immensely varied, ranging from extraordinarily durable to very transient. Three factors that influence the protection time are discussed in Section 3.6. The virus’s attributes, such as the replication speed, influence the antibody production time. A faster-replicating virus, for instance, an RNA virus, inclines to have a shorter antibody protection time after natural infection. Virus immunogenicity against T-cells influences the duration of antibody protection. Viruses with strong T-cell immunogenicity provide a strong stimulus for antibody proliferation, exhibiting extended and even lifelong protection. The antibody binding kinetics also contribute to the difference in its protection performance, as illustrated in Section 3.6.
V.
Why are there cases of self-reinfection?
Self-reinfection prevails in less virulent strains such as the Omicron variant. Our model provides a plausible mechanism underlying this unusual phenomenon. The individual might exhibit self-reinfection with the presence of a good binding antibody; repeated infection can also be displayed when a less virulent strain invades the patient. This can help explain why Omicron infections are inclined to be asymptomatic and self-reinfect.
An extreme case of self-reinfection is a chronic status in which the individual displays positive nucleic acid results for a long time. The chronic infection of SARS-CoV-2, together with other chronic infections, can also be explained by this model. It was demonstrated that a sub-equilibrium state could be maintained under specific parameter combinations with a relatively low antibody level, thus exhibiting long-term chronic inflammation. The long-term chronic infection is maintained by the constantly active (even weakly) antigens, which can come (i) from the non-coding human genome, where there are, for example, proteins of ancient viruses such as HER-V that produce antibodies a long time after their passage in the host genome, providing cross-protection against other infectious agents [76,77,78,79]; (ii) from the V(D)J mechanism of the innate immunity which can also give birth to antibodies with a large spectrum of actively participating to non-specific defenses against pathogens [80,81,82,83]; and (iii) from other vaccines (it is for example known that BCG reprograms the innate immunity [84]).
VI.
Why do vaccinations show considerable differences in protection?
According to our antibody kinetic model, the antibody interference effect is the most important contributing factor to the differences in vaccine protection time, apart from viral mutations. Large amounts of interfering antibodies can inhibit the production of high concentrations of fast-binding antibodies after vaccination, causing a decrease in protection efficiency and protection period. The antibody reservoir’s heterogeneity in the human population also influences the generation of those fast-binding antibodies. The presence of good templates could promote the production of fast-binding mature antibodies after a few rounds of somatic hypermutation. Those fast-binding antibodies could provide extended protection. Environmental antigenic substances also influence the protection duration by influencing the antibodies’ decay rate. Clinical reports [85] have demonstrated that age and gender are statistically associated with differences in antibody response after vaccination. The IgG antibodies triggered by the SARS-CoV-2 BNT162b2 vaccine significantly varied among ages. Young people tend to generate more neutralizing antibodies compared to their elder counterparts. This mainly reflects a variation in k 4 values in our model. A faster antibody generation capacity, equal to a larger k 4 value, would confer a stronger immunity to younger people. However, it does not guarantee that the protection duration in young people would be longer than that in older people for each individual. As discussed in Section 3.6, the antibody attributes, together with environmental antigenic substances, strongly impact the protection time. Suitable antibodies with fast binding constant k 2 and high concentrations of environmental antigenic substances would prolong the protection time. The clinical report also indicates that females have stronger immunity than males. However, this trend is as insignificant as the influences of age. In this study, the forecast is performed to predict the population’s behavior, which does not consider the influences of age and gender. A more accurate and specific forecast can be performed in a future study when we integrate age and gender differences into the model. In this case, a different k 4 value would be assigned to each subgroup.
VII.
How can we improve the protective efficiency and duration of vaccines?
Persistent exposure to massive antigenic substances can boost the neutralizing antibody level in a short period. However, the repeated booster can hardly magnify the proportion of solid-binding antibodies in the final antibody composition. The targeted induction of antibodies with excellent binding kinetics could provide prolonged protection against reinfection. Researchers are paving the way in this direction through gene-editing techniques [59,60,61]. Another concern with the constant booster shots is that the antibody interference might diminish the effectiveness of other vaccines, including the expected vaccines against massively mutated strains of SARS-CoV-2. Therefore, instead of pursuing a short-term antibody surge, we suggest that researchers aim to induce antibodies with fast-binding activity.
Besides providing a quantitative explanation of virus–host dynamics during the infection cycle, this model has several practical applications. One application is to predict the evolutionary direction of SARS-CoV-2 mathematically. The relationship between virulence and transmissibility can be simulated. A delicate equilibrium point that optimizes the transmissibility can be numerically obtained. Based on this model, we predict that the virulence of SARS-CoV-2 might further decrease, accompanied by an enhancement of transmissibility [86]. The second application is to optimize the vaccine inoculum dose mathematically. The antigen level, which can be represented as the initial inoculum dose, significantly influences the vaccine’s efficiency. A low dose might not be able to induce sufficient IgG to provide durable protection. In contrast, a high dose is more inclined to trigger a high level of IgG. However, it also has a more substantial adverse effect, representing a high antibody–virus complex level in this model. The high dose also increases its production cost. Therefore, a mathematical optimization can be performed to evaluate the optimal inoculum dose before clinical trials. Thirdly, this study could theoretically pave the way for future vaccine development. The protection duration of antibodies generated by natural infection is explicitly modeled in this study. The modeling of vaccines is different because the vaccine, independent of its type, does not replicate but is injected with an extremely high initial viral concentration. Different vaccines would have a different protection duration against breakthrough infection. vaccination efficiency is very complicated and generally weaker than the protection brought by a natural infection. The inactivation process for the inactivated vaccine would truncate its original structure, especially the epitope spots of spikes and nucleocapsid proteins. Its efficiency might be less than the mRNA vaccine. The vaccine based on the viral vector also generated deficient antigens compared to the real virus. After all, two major factors influence the protection duration of different vaccines: the peak level of neutralizing antibodies and viral mutations. A large dose of inoculum would trigger more neutralizing antibodies, leading to more extended protection and stronger adverse effects caused by the immune response. All vaccines would display weakened or even zero protection against variants with tremendous mutations due to a reduced k 2 value in this model. mRNA vaccine is better than the traditional vaccine because it can gradually generate realistic antigens. The traditional concept for vaccine development, with the application of an inactivated virus, viral vector, or mRNA, strictly forbids the use of a live virus. Our model demonstrated that the attenuated virus with proliferation constraints (smaller k 1 value in our model) could induce the immune response and antibody level in a milder way (shown in Figure 5). Those attenuated live viruses form a solid basis for asymptomatic infections [87]. There are at least two remarkable advantages of attenuated live viruses in future vaccine development. Firstly, since our model demonstrates that the peak antibody level is independent of the initial virus concentration, a minimal dose of live viruses could elevate specific IgG to a high level. It is also more infectious than traditional vaccines. This means a large number of populations can be covered with a few vaccination attempts. Secondly, despite the unlikely inactivation process, the attenuated virus is structurally identical to the original virus. Therefore, the specificity of neutralizing antibodies triggered by those attenuated viruses is the same as that of those induced by natural infections. As demonstrated by our model, an excellent binding capacity could induce a more durable protection time. The safety of those attenuated live viruses can also be guaranteed based on the principle of RNA replication: truncated mRNA cannot produce fixed offspring (full genome size).
Nevertheless, we also have to admit that our model has many hypotheses and uncertainties. Furthermore, it still lacks an in vivo experiment data fitting process, while many of the predictions and theories in the article remain to be confirmed experimentally and statistically in the future.

Author Contributions

Conceptualization, Z.X.; methodology, Z.X. and J.D.; writing—original draft preparation, Z.X.; writing—review and editing, D.W. and H.Z.; funding acquisition, Z.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by DeZhou University, grant number 30101418.

Data Availability Statement

Matlab source codes are available at https://github.com/zhaobinxu23/antibody_dynamics, accessed on 1 January 2022.

Acknowledgments

We thank Zuyi Huang from Villanova University. Yushan Zhu from Tsinghua University, Jian Song, Peiyan Guan, XiangYong Li, and Zhenlin Wei from Dezhou University for helpful conversations, comments, and clarifications.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Skegg, D.; Gluckman, P.; Boulton, G.; Hackmann, H.; Karim, S.S.A.; Piot, P.; Woopen, C. Future scenarios for the COVID-19 pandemic. Lancet 2021, 397, 777–778. [Google Scholar] [CrossRef] [PubMed]
  2. Xu, Z.; Zhang, H. If we cannot eliminate them, should we tame them? Mathematics underpinning the dose effect of virus infection and its application on COVID-19 virulence evolution. medRxiv 2021. [Google Scholar] [CrossRef]
  3. Xu, Z.; Zhang, H.; Huang, Z. A Continuous Markov-Chain Model for the Simulation of COVID-19 Epidemic Dynamics. Biology 2022, 11, 190. [Google Scholar] [CrossRef] [PubMed]
  4. Parrino, J.; Graham, B.S. Smallpox vaccines: Past, present, and future. J. Allergy Clin. Immunol. 2006, 118, 1320–1326. [Google Scholar] [CrossRef] [PubMed]
  5. Sadoff, J.; Gray, G.; Vandebosch, A.; Cardenas, V.; Shukarev, G.; Grinsztejn, B.; Goepfert, P.A.; Truyers, C.; Fennema, H.; Spiessens, B.; et al. Safety and efficacy of single-dose Ad26. COV2. S vaccine against COVID-19. N. Engl. J. Med. 2021, 384, 2187–2201. [Google Scholar] [CrossRef] [PubMed]
  6. Allen, J.D.; Feng, W.; Corlin, L.; Porteny, T.; Acevedo, A.; Schildkraut, D.; King, E.; Ladin, K.; Fu, Q.; Stopka, T.J. Why are some people reluctant to be vaccinated for COVID-19? A cross-sectional survey among US Adults in May–June 2020. Prev. Med. Rep. 2021, 24, 101494. [Google Scholar] [CrossRef]
  7. Ruiz, J.B.; Bell, R.A. Predictors of intention to vaccinate against COVID-19: Results of a nationwide survey. Vaccine 2021, 39, 1080–1086. [Google Scholar] [CrossRef]
  8. Klugar, M.; Riad, A.; Mekhemar, M.; Conrad, J.; Buchbender, M.; Howaldt, H.P.; Attia, S. Side effects of mRNA-based and viral vector-based COVID-19 vaccines among German healthcare workers. Biology 2021, 10, 752. [Google Scholar] [CrossRef]
  9. Riad, A.; Pokorna, A.; Attia, S.; Klugarova, J.; Koscik, M.; Klugar, M. Prevalence of COVID-19 Vaccine Side Effects among Healthcare Workers in the Czech Republic. J. Clin. Med. 2021, 10, 1428. [Google Scholar] [CrossRef]
  10. Rzymski, P.; Camargo, C.A., Jr.; Fal, A.; Flisiak, R.; Gwenzi, W.; Kelishadi, R.; Leemans, A.; Nieto, J.J.; Ozen, A.; Perc, M.; et al. COVID-19 Vaccine Boosters: The good, the bad, and the ugly. Vaccines 2021, 9, 1299. [Google Scholar] [CrossRef]
  11. Gupta, R.K.; Topol, E.J. COVID-19 vaccine breakthrough infections. Science 2021, 374, 1561–1562. [Google Scholar] [CrossRef] [PubMed]
  12. Cevik, M.; Grubaugh, N.D.; Iwasaki, A.; Openshaw, P. COVID-19 vaccines: Keeping pace with SARS-CoV-2 variants. Cell 2021, 184, 5077–5081. [Google Scholar] [CrossRef] [PubMed]
  13. Krause, P.R.; Fleming, T.R.; Peto, R.; Longini, I.M.; Figueroa, J.P.; Sterne, J.A.C.; Cravioto, A.; Rees, H.; Higgins, J.P.T.; Boutron, I.; et al. Considerations in boosting COVID-19 vaccine immune responses. Lancet 2021, 398, 1377–1380. [Google Scholar] [CrossRef] [PubMed]
  14. Klein, S.L.; Shann, F.; Moss, W.J.; Benn, C.S.; Aabye, P. RTS, S malaria vaccine and increased mortality in girls. Mbio 2016, 7, e00514-16. [Google Scholar] [CrossRef] [Green Version]
  15. Bocharov, G.; Volpert, V.; Ludewig, B.; Meyerhans, A. Mathematical Immunology of Virus Infections; Springer International Publishing: Cham, Switzerland, 2018. [Google Scholar]
  16. Eftimie, R.; Gillard, J.J.; Cantrell, D.A. Mathematical models for immunology: Current state of the art and future research directions. Bull. Math. Biol. 2016, 78, 2091–2134. [Google Scholar] [CrossRef] [Green Version]
  17. Smith, A.M.; Adler, F.R.; McAuley, J.L.; Gutenkunst, R.N.; Ribeiro, R.M.; McCullers, J.A.; Perelson, A.S. Effect of 1918 PB1-F2 expression on influenza A virus infection kinetics. PLoS Comput. Biol. 2011, 7, e1001081. [Google Scholar] [CrossRef] [Green Version]
  18. Baccam, P.; Beauchemin, C.; Macken, C.A.; Hayden, F.G.; Perelson, A.S. Kinetics of influenza A virus infection in humans. J. Virol. 2006, 80, 7590–7599. [Google Scholar] [CrossRef] [Green Version]
  19. Beauchemin, C.A.; McSharry, J.J.; Drusano, G.L.; Nguyen, J.T.; Went, G.T.; Ribeiro, R.M.; Perelson, A.S. Modeling amantadine treatment of influenza A virus in vitro. J. Theor. Biol. 2008, 254, 439–451. [Google Scholar] [CrossRef] [Green Version]
  20. Handel, A.; Longini, I.M., Jr.; Antia, R. Neuraminidase inhibitor resistance in influenza: Assessing the danger of its generation and spread. PLoS Comput. Biol. 2007, 3, e240. [Google Scholar] [CrossRef]
  21. Canini, L.; Perelson, A.S. Viral kinetic modeling: State of the art. J. Pharmacokinet. Pharmacodyn. 2014, 41, 431–443. [Google Scholar] [CrossRef] [Green Version]
  22. Khajanchi, S.; Sarkar, K. Forecasting the daily and cumulative number of cases for the COVID-19 pandemic in India. Chaos Interdiscip. J. Nonlinear Sci. 2020, 30, 071101. [Google Scholar] [CrossRef]
  23. Mondal, J.; Khajanchi, S. Mathematical modeling and optimal intervention strategies of the COVID-19 outbreak. Nonlinear Dyn. 2022, 109, 177–202. [Google Scholar] [CrossRef] [PubMed]
  24. Hancioglu, B.; Swigon, D.; Clermont, G. A dynamical model of human immune response to influenza A virus infection. J. Theor. Biol. 2007, 246, 70–86. [Google Scholar] [CrossRef] [PubMed]
  25. Handel, A.; Longini, I.M., Jr.; Antia, R. Towards a quantitative understanding of the within-host dynamics of influenza a infections. J. R. Soc. Interface 2010, 7, 35–47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Lee, H.Y.; Topham, D.J.; Park, S.Y.; Hollenbaugh, J.; Treanor, J.; Mosmann, T.R.; Jin, X.; Ward, B.M.; Miao, H.; Holden-Wiltse, J.; et al. Simulation and prediction of the adaptive immune response to influenza A virus infection. J. Virol. 2009, 83, 7151–7165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Miao, H.; Hollenbaugh, J.A.; Zand, M.S.; Holden-Wiltse, J.; Mosmann, T.R.; Perelson, A.S.; Wu, H.; Topham, D.J. Quantifying the early immune response and adaptive immune response kinetics in mice infected with influenza A virus. J. Virol. 2010, 84, 6687–6698. [Google Scholar] [CrossRef] [Green Version]
  28. Tridane, A.; Kuang, Y. Modeling the interaction of cytotoxic T lymphocytes and influenza virus infected epithelial cells. Math. Biosci. Eng. 2010, 7, 171. [Google Scholar]
  29. Clapham, H.E.; Quyen, T.H.; Kien, D.T.H.; Dorigatti, I.; Simmons, C.; Ferguson, N. Modelling virus and antibody dynamics during dengue virus infection suggests a role for antibody in virus clearance. PLoS Comput. Biol. 2016, 12, e1004951. [Google Scholar] [CrossRef] [Green Version]
  30. Best, K.; Guedj, J.; Madelain, V.; de Lamballerie, X.; Lim, S.Y.; Osuna, C.E.; Whitney, J.B.; Perelson, A.S. Zika plasma viral dynamics in nonhuman primates provides insights into early infection and antiviral strategies. Proc. Natl. Acad. Sci. USA 2017, 114, 8847–8852. [Google Scholar] [CrossRef] [Green Version]
  31. Clark, E.A.; Ledbetter, J.A. How B and T cells talk to each other. Nature 1994, 367, 425–428. [Google Scholar] [CrossRef]
  32. Mujal, A.M.; Delconte, R.B.; Sun, J.C. Natural killer cells: From innate to adaptive features. Annu. Rev. Immunol. 2021, 39, 417–447. [Google Scholar] [CrossRef] [PubMed]
  33. Shampine, L.F.; Reichelt, M.W. The matlab ode suite. SIAM J. Sci. Comput. 1997, 18, 1–22. [Google Scholar] [CrossRef] [Green Version]
  34. Abbas, A.K.; Lichtman, A.H.; Pober, J.S. B cell activation and antibody production. Cell. Mol. Immunol. 2005, 1, 243–267. [Google Scholar]
  35. Guermonprez, P.; Valladeau, J.; Zitvogel, L.; Théry, C.; Amigorena, S. Antigen presentation and T cell stimulation by dendritic cells. Annu. Rev. Immunol. 2002, 20, 621–667. [Google Scholar] [CrossRef]
  36. Bonyah, E.; Okosun, K.O. Mathematical modeling of Zika virus. Asian Pac. J. Trop. Dis. 2016, 6, 673–679. [Google Scholar] [CrossRef]
  37. Yamayoshi, S.; Yasuhara, A.; Ito, M.; Akasaka, O.; Nakamura, M.; Nakachi, I.; Koga, M.; Mitamura, K.; Yagi, K.; Maeda, K.; et al. Antibody titers against SARS-CoV-2 decline, but do not disappear for several months. EClinicalMedicine 2021, 32, 100734. [Google Scholar] [CrossRef]
  38. Pennock, N.D.; White, J.T.; Cross, E.W.; Cheney, E.E.; Tamburini, B.A.; Kedl, R.M. T cell responses: Naive to memory and everything in between. Adv. Physiol. Educ. 2013, 37, 273–283. [Google Scholar] [CrossRef] [Green Version]
  39. Kurosaki, T.; Kometani, K.; Ise, W. Memory B cells. Nat. Rev. Immunol. 2015, 15, 149–159. [Google Scholar] [CrossRef]
  40. Inoue, T.; Moran, I.; Shinnakasu, R.; Phan, T.G.; Kurosaki, T. Generation of memory B cells and their reactivation. Immunol. Rev. 2018, 283, 138–149. [Google Scholar] [CrossRef]
  41. Van den Berg, S.P.H.; Derksen, L.Y.; Drylewicz, J.; Nanlohy, N.M.; Beckers, L.; Lanfermeijer, J.; Gessel, S.N.; Vos, M.; Otto, S.A.; de Boer, R.J.; et al. Quantification of T-cell dynamics during latent cytomegalovirus infection in humans. PLoS Pathog. 2021, 17, e1010152. [Google Scholar] [CrossRef]
  42. Aalberse, R.C.; Kleine, B.I.; Stapel, S.O.; van Ree, R. Structural aspects of cross-reactivity and its relation to antibody affinity. Allergy 2001, 56, 27–29. [Google Scholar] [CrossRef] [PubMed]
  43. Bull, J.J.; Lauring, A.S. Theory and empiricism in virulence evolution. PLoS Pathog. 2014, 10, e1004387. [Google Scholar] [CrossRef] [PubMed]
  44. Iqbal, F.M.; Lam, K.; Sounderajah, V.; Clarke, J.M.; Ashrafian, H.; Darzi, A. Characteristics and predictors of acute and chronic post-COVID syndrome: A systematic review and meta-analysis. EClinicalMedicine 2021, 36, 100899. [Google Scholar] [CrossRef] [PubMed]
  45. Fajgenbaum, D.C.; June, C.H. Cytokine storm. N. Engl. J. Med. 2020, 383, 2255–2273. [Google Scholar] [CrossRef] [PubMed]
  46. Kastritis, P.L.; Bonvin, A.M.J.J. On the binding affinity of macromolecular interactions: Daring to ask why proteins interact. J. R. Soc. Interface 2013, 10, 20120835. [Google Scholar] [CrossRef]
  47. Seydoux, E.; Homad, L.J.; MacCamy, A.J.; Parks, K.R.; Hurlburt, N.K.; Jennewein, M.F.; Akins, N.R.; Stuart, A.B.; Wan, Y.H.; Feng, J.; et al. Analysis of a SARS-CoV-2-infected individual reveals development of potent neutralizing antibodies with limited somatic mutation. Immunity 2020, 53, 98–105.e5. [Google Scholar] [CrossRef] [PubMed]
  48. Ju, B.; Zhang, Q.; Ge, J.; Wang, R.; Sun, J.; Ge, X.; Yu, J.; Shan, S.; Zhou, B.; Song, S.; et al. Human neutralizing antibodies elicited by SARS-CoV-2 infection. Nature 2020, 584, 115–119. [Google Scholar] [CrossRef]
  49. Muecksch, F.; Weisblum, Y.; Barnes, C.O.; Schmidt, F.; Schaefer-Babajew, D.; Wang, Z.; Lorenzi, J.C.; Flyak, A.I.; DeLaitsch, A.T.; Huey-Tubman, K.E.; et al. Affinity maturation of SARS-CoV-2 neutralizing antibodies confers potency, breadth, and resilience to viral escape mutations. Immunity 2021, 54, 1853–1868.e7. [Google Scholar] [CrossRef]
  50. Schneider, E.D.; Kay, J.J. Order from disorder: The thermodynamics of complexity in biology. In What Is Life? The Next Fifty Years: Speculations on the Future of Biology; Murphy, M.P., O’neill, L.A.J., Eds.; Cambridge University Press: Cambridge, UK, 1997; p. 161. [Google Scholar]
  51. Marín-Hernández, D.; Schwartz, R.E.; Nixon, D.F. Epidemiological evidence for association between higher influenza vaccine uptake in the elderly and lower COVID-19 deaths in Italy. J. Med. Virol. 2021, 93, 64. [Google Scholar] [CrossRef]
  52. Salem, M.L.; El-Hennawy, D. The possible beneficial adjuvant effect of influenza vaccine to minimize the severity of COVID-19. Med. Hypotheses 2020, 140, 109752. [Google Scholar] [CrossRef]
  53. Fink, G.; Orlova-Fink, N.; Schindler, T.; Grisi, S.; Ferrer, A.P.S.; Daubenberger, C.; Brentani, A. Inactivated trivalent influenza vaccination is associated with lower mortality among patients with COVID-19 in Brazil. BMJ Evid.-Based Med. 2021, 26, 192–193. [Google Scholar] [CrossRef] [PubMed]
  54. Escobar, L.E.; Molina-Cruz, A.; Barillas-Mury, C. BCG vaccine protection from severe coronavirus disease 2019 (COVID-19). Proc. Natl. Acad. Sci. USA 2020, 117, 17720–17726. [Google Scholar] [CrossRef] [PubMed]
  55. Aaby, P.; Mogensen, S.W.; Rodrigues, A.; Benn, C.S. Evidence of increase in mortality after the introduction of diphtheria-tetanus- pertussis vaccine to children aged 6–35 months in Guinea-Bissau: A time for reflection? Front. Public Health 2018, 6, 79. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Rhee, P.; Nunley, M.K.; Demetriades, D.; Velmahos, G.; Doucet, J.J. Tetanus and trauma: A review and recommendations. J. Trauma Acute Care Surg. 2005, 58, 1082–1088. [Google Scholar] [CrossRef] [PubMed]
  57. Milne, G.; Hames, T.; Scotton, C.; Gent, N.; Johnsen, A.; Anderson, R.M.; Ward, T. Does infection with or vaccination against SARS-CoV-2 lead to lasting immunity? Lancet Respir. Med. 2021, 9, 1450–1466. [Google Scholar] [CrossRef]
  58. Ferdinands, J.M.; Alyanak, E.; Reed, C.; Fry, A.M. Waning of influenza vaccine protection: Exploring the trade-offs of changes in vaccination timing among older adults. Clin. Infect. Dis. 2020, 70, 1550–1559. [Google Scholar] [CrossRef] [Green Version]
  59. Huang, D.; Tran, J.T.; Olson, A.; Vollbrecht, T.; Tenuta, M.; Guryleva, M.V.; Fuller, R.P.; Schiffner, T.; Abadejos, J.R.; Couvrette, L.; et al. Vaccine elicitation of HIV broadly neutralizing antibodies from engineered B cells. Nat. Commun. 2020, 11, 5850. [Google Scholar] [CrossRef]
  60. Voss, J.E.; Gonzalez-Martin, A.; Andrabi, R.; Fuller, R.P.; Murrell, B.; McCoy, L.E.; Porter, K.; Huang, D.; Li, W.; Sok, D.; et al. Reprogramming the antigen specificity of B cells using genome-editing technologies. Elife 2019, 8, e42995. [Google Scholar] [CrossRef]
  61. Moffett, H.F.; Harms, C.K.; Fitzpatrick, K.S.; Tooley, M.R.; Boonyaratanakornkit, J.; Taylor, J.J. B cells engineered to express pathogen-specific antibodies protect against infection. Sci. Immunol. 2019, 4, eaax0644. [Google Scholar] [CrossRef]
  62. Lumley, S.F.; Wei, J.; O’Donnell, D.; Stoesser, N.E.; Matthews, P.C.; Howarth, A.; Hatch, S.B.; Marsden, B.D.; Cox, S.; James, T.; et al. The duration, dynamics and determinants of SARS-CoV-2 antibody responses in individual healthcare workers. Clin. Infect. Dis. 2021. [Google Scholar] [CrossRef]
  63. Cohn, B.A.; Cirillo, P.M.; Murphy, C.C.; Krigbaum, N.Y.; Wallace, A.W. SARS-CoV-2 vaccine protection and deaths among US veterans during 2021. Science 2022, 375, 331–336. [Google Scholar] [CrossRef] [PubMed]
  64. Tartof, S.Y.; Slezak, J.M.; Fischer, H.; Hong, V.; Ackerson, B.K.; Ranasinghe, O.N.; Frankland, T.B.; Ogun, O.A.; Zamparo, J.M.; Gray, S.; et al. Effectiveness of mRNA BNT162b2 COVID-19 vaccine up to 6 months in a large integrated health system in the USA: A retrospective cohort study. Lancet 2021, 398, 1407–1416. [Google Scholar] [CrossRef] [PubMed]
  65. Chemaitelly, H.; Tang, P.; Hasan, M.R.; AlMukdad, S.; Yassine, H.M.; Benslimane, F.M.; Khatib, H.A.; Coyle, P.; Ayoub, H.H.; Kanaani, Z.; et al. Waning of BNT162b2 vaccine protection against SARS-CoV-2 infection in Qatar. N. Engl. J. Med. 2021, 385, e83. [Google Scholar] [CrossRef]
  66. Dao, T.L.; Hoang, V.T.; Gautret, P. Recurrence of SARS-CoV-2 viral RNA in recovered COVID-19 patients: A narrative review. Eur. J. Clin. Microbiol. Infect. Dis. 2021, 40, 13–25. [Google Scholar] [CrossRef]
  67. Azam, M.; Sulistiana, R.; Ratnawati, M.; Fibriana, A.I.; Bahrudin, U.; Widyaningrum, D.; Aljunid, S.M. Recurrent SARS-CoV-2 RNA positivity after COVID-19: A systematic review and meta-analysis. Sci. Rep. 2020, 10, 20692. [Google Scholar] [CrossRef]
  68. Bush, A. Recurrent respiratory infections. Pediatr. Clin. N. Am. 2009, 56, 67–100. [Google Scholar] [CrossRef] [PubMed]
  69. Midgard, H.; Weir, A.; Palmateer, N.; Lo Re, V., 3rd; Pineda, J.A.; Macias, J.; Dalgard, O. HCV epidemiology in high-risk groups and the risk of reinfection. J. Hepatol. 2016, 65, S33–S45. [Google Scholar] [CrossRef] [Green Version]
  70. Yang, J.R.; Deng, D.T.; Wu, N.; Yang, B.; Li, H.J.; Pan, X.B. Persistent viral RNA positivity during the recovery period of a patient with SARS-CoV-2 infection. J. Med. Virol. 2020, 92, 1681–1683. [Google Scholar] [CrossRef]
  71. Chan, C.E.; Chan, A.H.; Hanson, B.J.; Ooi, E.E. The use of antibodies in the treatment of infectious diseases. Singap. Med. J. 2009, 50, 663–672. [Google Scholar]
  72. Rehermann, B.; Nascimbeni, M. Immunology of hepatitis B virus and hepatitis C virus infection. Nat. Rev. Immunol. 2005, 5, 215–229. [Google Scholar] [CrossRef]
  73. Durham, S.R. Allergen avoidance measures. Respir. Med. 1996, 90, 441–445. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Rai, R.K.; Khajanchi, S.; Tiwari, P.K.; Venturino, E.; Misra, A.K. Impact of social media advertisements on the transmission dynamics of COVID-19 pandemic in India. J. Appl. Math. Comput. 2022, 68, 19–44. [Google Scholar] [CrossRef] [PubMed]
  75. Khajanchi, S.; Sarkar, K.; Mondal, J.; Nisar, K.S.; Abdelwahab, S.F. Mathematical modeling of the COVID-19 pandemic with intervention strategies. Results Phys. 2021, 25, 104285. [Google Scholar] [CrossRef] [PubMed]
  76. Laderoute, M.P.; Larocque, L.J.; Giulivi, A.; Diaz-Mitoma, F. Further evidence that human endogenous retrovirus K102 is a replication competent foamy virus that may antagonize HIV-1 replication. Open AIDS J. 2015, 9, 112. [Google Scholar] [CrossRef] [PubMed]
  77. Demongeot, J.; Seligmann, H. SARS-CoV-2 and miRNA-like inhibition power. Med. Hypotheses 2020, 144, 110245. [Google Scholar] [CrossRef] [PubMed]
  78. Srinivasachar Badarinarayan, S.; Sauter, D. Switching sides: How endogenous retroviruses protect us from viral infections. J. Virol. 2021, 95, e02299-20. [Google Scholar] [CrossRef] [PubMed]
  79. Liu, H.; Bergant, V.; Frishman, G.; Ruepp, A.; Pichlmair, A.; Vincendeau, M.; Frishman, D. Influenza A Virus Infection Reactivates Human Endogenous Retroviruses Associated with Modulation of Antiviral Immunity. Viruses 2022, 14, 1591. [Google Scholar] [CrossRef]
  80. Pasqual, N.; Gallagher, M.; Aude-Garcia, C.; Loiodice, M.; Thuderoz, F.; Demongeot, J.; Ceredig, R.; Marche, P.N.; Jouvin-Marche, E. Quantitative and qualitative changes in VJ α rearrangements during mouse thymocytes differentiation: Implication for a limited T cell receptor α chain repertoire. J. Exp. Med. 2002, 196, 1163–1174. [Google Scholar] [CrossRef] [Green Version]
  81. Baum, T.P.; Pasqual, N.; Thuderoz, F.; Hierle, V.; Chaume, D.; Lefranc, M.P.; Jouvin-Marche, E.; Marche, P.N.; Demongeot, J. IMGT/GeneInfo: Enhancing V (D) J recombination database accessibility. Nucleic Acids Res. 2004, 32, D51–D54. [Google Scholar] [CrossRef] [Green Version]
  82. Baum, T.P.; Hierle, V.; Pasqual, N.; Bellahcene, F.; Chaume, D.; Lefranc, M.P.; Jouvin-Marche, E.; Marche, P.N.; Demongeot, J. IMGT/GeneInfo: T cell receptor gamma TRG and delta TRD genes in database give access to all TR potential V (D) J recombinations. BMC Bioinform. 2006, 7, 1–7. [Google Scholar] [CrossRef]
  83. Thuderoz, F.; Simonet, M.A.; Hansen, O.; Pasqual, N.; Dariz, A.; Baum, T.P.; Hierle, V.; Demongeot, J.; Marche, P.N.; Jouvin-Marche, E. Numerical modelling of the VJ combinations of the T cell receptor TRA/TRD locus. PLoS Comput. Biol. 2010, 6, e1000682. [Google Scholar] [CrossRef] [PubMed]
  84. O’Neill, L.A.J.; Netea, M.G. BCG-induced trained immunity: Can it offer protection against COVID-19? Nat. Rev. Immunol. 2020, 20, 335–337. [Google Scholar] [CrossRef] [PubMed]
  85. Pellini, R.; Venuti, A.; Pimpinelli, F.; Abril, E.; Blandino, G.; Campo, F.; Conti, L.; De Virgilio, A.; De Marco, F.; Di Domenico, E.G.; et al. Initial observations on age, gender, BMI and hypertension in antibody responses to SARS-CoV-2 BNT162b2 vaccine. EClinicalMedicine 2021, 36, 100928. [Google Scholar] [CrossRef] [PubMed]
  86. Xu, Z.; Wei, D.; Zeng, Q.; Zhang, H.; Sun, Y.; Demongeot, J. More or less deadly? A mathematical model that predicts SARS-CoV-2 evolutionary direction. Comput. Biol. Med. 2023, 153, 106510. [Google Scholar] [CrossRef]
  87. Xu, Z.; Yang, D.; Wang, L.; Demongeot, J. Statistical analysis supports UTR (untranslated region) deletion theory in SARS-CoV-2. Virulence 2022, 13, 1772–1789. [Google Scholar] [CrossRef]
Figure 1. A simple diagram of host-virus interaction.
Figure 1. A simple diagram of host-virus interaction.
Viruses 15 00586 g001
Figure 2. The role of environmental antigens in immune response.
Figure 2. The role of environmental antigens in immune response.
Viruses 15 00586 g002
Figure 3. Antibody and virus dynamics modeling using different p(0) value. p(0) represents the concentration of environmental antigens. The virus-antibody dynamics are modeled under different environmental antigen concentrations (A) and different environmental antigen attributes (B). As shown in (A), the antibody decay rate is significantly slower (shown in solid yellow curve) when there is a large amount of environmental antigen-like substances. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 10, p(0) = 1000 or p(0) = 1,000,000, k 1 = 0.1 ,   k 2 = 1 × 10−5, k 2 = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−8, k 7 = 1 × 10−14. As shown in (B), the antibody decay rate is significantly slower (shown in solid yellow curve) when k 8 sets a large value which corresponds to a better binding kinetics between environmental antigen-like stuff and its corresponding antibody. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 10, p(0) = 1,000,000, k 1 = 0.1 ,   k 2 = 1 × 10−5, k 2   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−8 or k 7 = 1.8 × 10−8, k 7 = 1 × 10−14.
Figure 3. Antibody and virus dynamics modeling using different p(0) value. p(0) represents the concentration of environmental antigens. The virus-antibody dynamics are modeled under different environmental antigen concentrations (A) and different environmental antigen attributes (B). As shown in (A), the antibody decay rate is significantly slower (shown in solid yellow curve) when there is a large amount of environmental antigen-like substances. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 10, p(0) = 1000 or p(0) = 1,000,000, k 1 = 0.1 ,   k 2 = 1 × 10−5, k 2 = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−8, k 7 = 1 × 10−14. As shown in (B), the antibody decay rate is significantly slower (shown in solid yellow curve) when k 8 sets a large value which corresponds to a better binding kinetics between environmental antigen-like stuff and its corresponding antibody. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 10, p(0) = 1,000,000, k 1 = 0.1 ,   k 2 = 1 × 10−5, k 2   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−8 or k 7 = 1.8 × 10−8, k 7 = 1 × 10−14.
Viruses 15 00586 g003
Figure 4. Antibody dynamics modeling with different k 7   values. (A) One scenario where environmental antigen-like substances do not trigger antibody growth. As shown in (A), the antibody does not engage proliferation due to the presence of environment antigen-like molecules. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 0, p(0) = 1,000,000, k 1 = 0.1 ,   k 2 = 1 × 10−5, k 2 = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−8, k 7 = 1 × 10−14. (B) One scenario where environmental antigen-like substances do trigger antibody proliferation. As shown in (B), the antibody does engage proliferation due to the presence of environment antigen-like molecules. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 0, p(0) = 1,000,000, k 1 = 0.1 ,   k 2   = 1 × 10−5, k 2   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−7, k 7 = 1 × 10−14. The antibodies might significantly increase when the environmental antigenic substances bind strongly with them. This always induces allergic reactions.
Figure 4. Antibody dynamics modeling with different k 7   values. (A) One scenario where environmental antigen-like substances do not trigger antibody growth. As shown in (A), the antibody does not engage proliferation due to the presence of environment antigen-like molecules. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 0, p(0) = 1,000,000, k 1 = 0.1 ,   k 2 = 1 × 10−5, k 2 = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−8, k 7 = 1 × 10−14. (B) One scenario where environmental antigen-like substances do trigger antibody proliferation. As shown in (B), the antibody does engage proliferation due to the presence of environment antigen-like molecules. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 0, p(0) = 1,000,000, k 1 = 0.1 ,   k 2   = 1 × 10−5, k 2   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−7, k 7 = 1 × 10−14. The antibodies might significantly increase when the environmental antigenic substances bind strongly with them. This always induces allergic reactions.
Viruses 15 00586 g004
Figure 5. Different immune behaviors toward variants with different replication activities. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 10, p(0) = 1000,   k 1 = 0.06   o r   0.1 ,   k 2 = 1 × 10−5, k 2   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−8, k 7 = 1 × 10−14. As shown in this figure, the antibody response is milder when the host is infected by a less toxic strain (smaller   k 1 value 0.06). The peak viral load is also less compared to its high toxic counterpart (bigger   k 1 value 0.1). This indicates that even for the same virus infection, the immune response would vary greatly due to the differences in viral replication capacity of different variants.
Figure 5. Different immune behaviors toward variants with different replication activities. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 10, p(0) = 1000,   k 1 = 0.06   o r   0.1 ,   k 2 = 1 × 10−5, k 2   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−8, k 7 = 1 × 10−14. As shown in this figure, the antibody response is milder when the host is infected by a less toxic strain (smaller   k 1 value 0.06). The peak viral load is also less compared to its high toxic counterpart (bigger   k 1 value 0.1). This indicates that even for the same virus infection, the immune response would vary greatly due to the differences in viral replication capacity of different variants.
Viruses 15 00586 g005
Figure 6. Dynamics of different antibodies with different kinetic attributes. The parameter sets we used are: x(0) = 0, y(0) = 1, z(0) = 1, w = 1, k 1 = 0.1 ,   k 2   = 1 × 10−5, k 2   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, for antibody 1;   k 2 = 9 × 10−6, k 2 = 9 × 10−15 for antibody 2; k 2 = 8 × 10−6, k 2 = 8 × 10−15 for antibody 3; k 2 = 7 × 10−6, k 2 = 7 × 10−15 for antibody 4; k 2 = 6 × 10−6, k 2 = 6 × 10−15 for antibody 5. It is demonstrated in this figure that the faster-binding antibodies engage amplification at a greater magnitude. In this way, the immune system selects those specific neutralizing antibodies.
Figure 6. Dynamics of different antibodies with different kinetic attributes. The parameter sets we used are: x(0) = 0, y(0) = 1, z(0) = 1, w = 1, k 1 = 0.1 ,   k 2   = 1 × 10−5, k 2   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, for antibody 1;   k 2 = 9 × 10−6, k 2 = 9 × 10−15 for antibody 2; k 2 = 8 × 10−6, k 2 = 8 × 10−15 for antibody 3; k 2 = 7 × 10−6, k 2 = 7 × 10−15 for antibody 4; k 2 = 6 × 10−6, k 2 = 6 × 10−15 for antibody 5. It is demonstrated in this figure that the faster-binding antibodies engage amplification at a greater magnitude. In this way, the immune system selects those specific neutralizing antibodies.
Viruses 15 00586 g006
Figure 7. High concentrations of weakly binding antibodies can provide effective protection. A plot of the inhibitory capacity of a specific concentration of weakly binding antibodies against infection is presented. Two types of antibodies are presented in this figure: antibody 1 has a strong binding capacity ( K 1 on   = 1 × 10−5) while antibody 2 has a relatively weak binding capacity ( K 2 on   = 5 × 10−6). Two scenarios are simulated: both antibodies have low initial concentrations in case 1; weakly binding antibody has a high initial level in case 2. The parameter sets we used are: x 1 (0) = x 2 (0) = 0, y 1 (0) = y 2 (0) = 1, z(0) = 1, w 1 = w 2 = 1, K 1 on   = 1 × 10−5, K 1 off   = 1 × 10−14, K 2 on   = 5 × 10−6, K 2 off   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = k 6 = 0.02, k 1 = 0.1 for case 1; y 2 (0) = 100 for case 2. It can be seen that the peak viral load and antibody level are significantly lower in case 2, which corresponds to a milder immune response. It indicates that the elevated weakly binding antibodies could also provide protection against severe infection.
Figure 7. High concentrations of weakly binding antibodies can provide effective protection. A plot of the inhibitory capacity of a specific concentration of weakly binding antibodies against infection is presented. Two types of antibodies are presented in this figure: antibody 1 has a strong binding capacity ( K 1 on   = 1 × 10−5) while antibody 2 has a relatively weak binding capacity ( K 2 on   = 5 × 10−6). Two scenarios are simulated: both antibodies have low initial concentrations in case 1; weakly binding antibody has a high initial level in case 2. The parameter sets we used are: x 1 (0) = x 2 (0) = 0, y 1 (0) = y 2 (0) = 1, z(0) = 1, w 1 = w 2 = 1, K 1 on   = 1 × 10−5, K 1 off   = 1 × 10−14, K 2 on   = 5 × 10−6, K 2 off   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = k 6 = 0.02, k 1 = 0.1 for case 1; y 2 (0) = 100 for case 2. It can be seen that the peak viral load and antibody level are significantly lower in case 2, which corresponds to a milder immune response. It indicates that the elevated weakly binding antibodies could also provide protection against severe infection.
Viruses 15 00586 g007
Figure 8. An illustration of protection time calculation. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 10, p(0) = 10,000, k 1 = 0.1 ,   k 2 = 1 × 10−7, k 2 = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−8, k 7 = 1 × 10−14. Incubation time is calculated as the time interval between virus entrance and the production of antibodies. The viruses would engage a proliferation earlier than the antibodies. The patient is still asymptomatic even though the virus has reached a high level. Symptoms such as fever would appear when the antibody–virus complexes reach beyond a specific threshold. The second infection is marked with a green arrow. It can be seen in this figure that a second infection could occur when the IgG level drops below a certain level. It does not necessarily require a zero IgG level when a breakthrough infection happens.
Figure 8. An illustration of protection time calculation. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 10, p(0) = 10,000, k 1 = 0.1 ,   k 2 = 1 × 10−7, k 2 = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−8, k 7 = 1 × 10−14. Incubation time is calculated as the time interval between virus entrance and the production of antibodies. The viruses would engage a proliferation earlier than the antibodies. The patient is still asymptomatic even though the virus has reached a high level. Symptoms such as fever would appear when the antibody–virus complexes reach beyond a specific threshold. The second infection is marked with a green arrow. It can be seen in this figure that a second infection could occur when the IgG level drops below a certain level. It does not necessarily require a zero IgG level when a breakthrough infection happens.
Viruses 15 00586 g008
Figure 9. Protection time calculation when the antibody has a specific binding kinetic constant k 2 . It can be inferred from this figure that a fast-binding neutralizing antibody would provide a longer protection time when we compared (A) with (C). The protection time against severe infection can also be prolonged, given the faster binding kinetics when comparing (B) to (D). (A) Protection time calculation when the antibody has a weak binding kinetic constant k 2 = 1 × 10−6. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 10, q(0) = 1 × 106, k 1 = 0.1 ,   k 2 = 1 × 10−6, k 2 = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−9, k 7 = 1 × 10−14. (B) Maximal virus load at different infection points when the antibody has a binding kinetic constant k 1 = 1 × 10−6. The parameter set used is the same as (A). (C) Protection time calculation when the antibody has a strong binding kinetic constant k 2 = 1 × 10−5. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 10, q(0) = 1 × 106, k 1 = 0.1 ,   k 2 = 1 × 10−5, k 2   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−9, k 7 = 1 × 10−14. (D) Maximal virus load at different infection points when the antibody has a binding kinetic constant k 2 = 1 × 10−5. The parameter set used is the same as (C).
Figure 9. Protection time calculation when the antibody has a specific binding kinetic constant k 2 . It can be inferred from this figure that a fast-binding neutralizing antibody would provide a longer protection time when we compared (A) with (C). The protection time against severe infection can also be prolonged, given the faster binding kinetics when comparing (B) to (D). (A) Protection time calculation when the antibody has a weak binding kinetic constant k 2 = 1 × 10−6. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 10, q(0) = 1 × 106, k 1 = 0.1 ,   k 2 = 1 × 10−6, k 2 = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−9, k 7 = 1 × 10−14. (B) Maximal virus load at different infection points when the antibody has a binding kinetic constant k 1 = 1 × 10−6. The parameter set used is the same as (A). (C) Protection time calculation when the antibody has a strong binding kinetic constant k 2 = 1 × 10−5. The parameter set we used is: x(0) = 0, y(0) = 100, z(0) = 10, q(0) = 1 × 106, k 1 = 0.1 ,   k 2 = 1 × 10−5, k 2   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02, k 7 = 1 × 10−9, k 7 = 1 × 10−14. (D) Maximal virus load at different infection points when the antibody has a binding kinetic constant k 2 = 1 × 10−5. The parameter set used is the same as (C).
Viruses 15 00586 g009
Figure 10. (A) the dynamic behavior of antibodies in the overall population through time. (The blue zone around mean curve stands for 95% confidence interval). It can be seen that the IgG level would significantly decline after reaching a peak level. However, its degradation does not follow a simple mathematical formula. Its descent rate would gradually decline and be maintained at a relatively stable level after 200 days. (B) The protection performance of antibodies in the overall population through time. It can be seen in this figure that the protection efficiency of induced neutralizing antibodies would be maintained at a relatively high level in the first 100 days. Its protection efficiency would engage a rapid decline after the first three months.
Figure 10. (A) the dynamic behavior of antibodies in the overall population through time. (The blue zone around mean curve stands for 95% confidence interval). It can be seen that the IgG level would significantly decline after reaching a peak level. However, its degradation does not follow a simple mathematical formula. Its descent rate would gradually decline and be maintained at a relatively stable level after 200 days. (B) The protection performance of antibodies in the overall population through time. It can be seen in this figure that the protection efficiency of induced neutralizing antibodies would be maintained at a relatively high level in the first 100 days. Its protection efficiency would engage a rapid decline after the first three months.
Viruses 15 00586 g010
Figure 11. (A) Self-reinfection scenario. The parameter sets we used are: x(0) = 0, y(0) = 1000, z(0) = 1, w = 1000, k 1 = 0.1 ,   k 2   = 1 × 10−5, k 2   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02. Reinfections are represented as repeated waves in this figure. It indicates that the viruses could re-proliferate when the antibodies cannot completely eliminate them. The viruses start to proliferate when the antibodies drop to a certain level, leading to self-reinfection in this case. (B) Scenario of chronic infection. The parameter sets we used are: x(0) = 0, y(0) = 1000, z(0) = 1, w = 1000, k 1 = 0.1, k 2 = 3 × 10−5, k 2   = 1 × 10−5, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02. In this case, pathogens would not be eliminated but maintained at a relatively stable level, forming a chronic infection. The low concentration of pathogenic antigens only provides a limited driving force for promoting antibody reproduction.
Figure 11. (A) Self-reinfection scenario. The parameter sets we used are: x(0) = 0, y(0) = 1000, z(0) = 1, w = 1000, k 1 = 0.1 ,   k 2   = 1 × 10−5, k 2   = 1 × 10−14, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02. Reinfections are represented as repeated waves in this figure. It indicates that the viruses could re-proliferate when the antibodies cannot completely eliminate them. The viruses start to proliferate when the antibodies drop to a certain level, leading to self-reinfection in this case. (B) Scenario of chronic infection. The parameter sets we used are: x(0) = 0, y(0) = 1000, z(0) = 1, w = 1000, k 1 = 0.1, k 2 = 3 × 10−5, k 2   = 1 × 10−5, k 3 = 1, k 4 = 2, k 5 = 0.02, k 6 = 0.02. In this case, pathogens would not be eliminated but maintained at a relatively stable level, forming a chronic infection. The low concentration of pathogenic antigens only provides a limited driving force for promoting antibody reproduction.
Viruses 15 00586 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xu, Z.; Wei, D.; Zhang, H.; Demongeot, J. A Novel Mathematical Model That Predicts the Protection Time of SARS-CoV-2 Antibodies. Viruses 2023, 15, 586. https://doi.org/10.3390/v15020586

AMA Style

Xu Z, Wei D, Zhang H, Demongeot J. A Novel Mathematical Model That Predicts the Protection Time of SARS-CoV-2 Antibodies. Viruses. 2023; 15(2):586. https://doi.org/10.3390/v15020586

Chicago/Turabian Style

Xu, Zhaobin, Dongqing Wei, Hongmei Zhang, and Jacques Demongeot. 2023. "A Novel Mathematical Model That Predicts the Protection Time of SARS-CoV-2 Antibodies" Viruses 15, no. 2: 586. https://doi.org/10.3390/v15020586

APA Style

Xu, Z., Wei, D., Zhang, H., & Demongeot, J. (2023). A Novel Mathematical Model That Predicts the Protection Time of SARS-CoV-2 Antibodies. Viruses, 15(2), 586. https://doi.org/10.3390/v15020586

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop