Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, August 13 until 12:00 AM ET on Friday, August 14 due to maintenance. We apologize for the inconvenience.


Title: Coupling compartmental models with Markov chains and measure evolution equations to capture virus mutability
The COVID-19 pandemic lit a fire under researchers who have subsequently raced to build models which capture various physical aspects of both the biology of the virus and its mobility throughout the human population. These models could include characteristics such as different genders, ages, frequency of interactions, mutation of virus, etc. Here, we propose two mathematical formulations to include virus mutation dynamics. The first uses a compartmental epidemiological model coupled with a discrete-time finite-state Markov chain. If one includes a nonlinear dependence of the transition matrix on current infected, the model is able to reproduce pandemic waves due to different variants. The second approach expands such an idea to a continuous state-space leveraging a combination of ordinary differential equations with an evolution equation for measure. This approach allows to include reinfections with partial immunity with respect to variants genetically similar to that of first infection.  more » « less
Award ID(s):
2033580
PAR ID:
10437523
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Mathematical Models and Methods in Applied Sciences
Volume:
32
Issue:
10
ISSN:
0218-2025
Page Range / eLocation ID:
2099 to 2119
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The emergence of new variants of the SARS-CoV-2 virus poses serious problems to the control of the current COVID-19 pandemic. Understanding how the variants originate is critical for effective control of the spread of the virus and the global pandemic. The study of the virus evolution so far has been dominated by phylogenetic tree analysis, which however is inappropriate for a few important reasons. Here we used phylogenetic network approach to study the origin of the VOC202012/01 (Alpha) or so-called UK variant (PANGO Lineage B.1.1.7). The multiple network analyses using different methods consistently revealed that the VOC202012/01 variant was a result of recombination, in contrast to the common assumption that the variant evolved from step-wise mutations in a linear order. The study provides an example for the power and application of phylogenetic network analysis in studying virus evolution, which can be applied to study the evolutionary processes leading to the emergence of other variants of the SARS-CoV-2 virus as well as many other viruses. 
    more » « less
  2. The objective is to understand the role of emerging variants, vaccination, and NPI policies on COVID-19 infections and deaths. We aim to identify scenarios in which COVID-19 can be managed such that the death rate from COVID-19 becomes comparable with the combined annual mortality rate from influenza and pneumonia. As a case study for a large urban area, we simulate COVID-19 transmission in King County, Washington, (greater Seattle) using an agent- based simulation model. Calibrated to local epidemiological data, our study uses detailed synthetic population data and includes interactions between vaccination and specific NPIs while considering waning immunity and emergence of variants. Virus mutation scenarios include 12 combinations of infectivity, disease severity, and immune evasiveness. A highly effective pancoronavirus vaccine that works against all strains is considered an optimistic scenario. Our findings highlight the potential benefits of pancoronavirus vaccines that offer enhanced and longer-lasting immunity. We emphasize the crucial role of nonpharmaceutical interventions in reducing COVID-19 deaths regardless of virus mutation scenarios. Owing to highly immune evasive and contagious SARS-CoV-2 variants, most scenarios in this study fail to reduce the mortality of COVID-19 to the level of influenza and pneumonia. However, our findings indicate that periodic vaccinations and a threshold nonpharmaceutical intervention policy may succeed in achieving this goal. This indicates the need for caution and vigilance in managing a continuing COVID-19 epidemic. 
    more » « less
  3. Compartmental models are essential for studying host-pathogen dynamics, evaluating intervention effectiveness, and predicting infection trends. However, the utility of these models for testing competing hypotheses is often overlooked. To address this, we propose a new model-based hypothesis testing (MBHT) approach, which uses compartmental models to evaluate the hypotheses in epidemiology. In our case, using the COVID-19 pandemic as a case study, we formulate hypotheses of SARS-CoV-2 mutation and construct a transmission model to test them. In addition to analyzing steady-state stability, deriving the basic reproduction number, and identifying a backward bifurcation, the model is fitted to seven peaks of U.S. COVID-19 data, each corresponding to periods of viral mutation and morbidity peaks. The estimated posterior probabilities reveal that Short-term within host selection primarily shaped mutations during the early pandemic stages, followed by immune selection driven by natural and vaccine-induced immunity. In later stages, mutations aligned with vaccination-induced virulence and transmission-virulence correlation, while the declining virulence and immune selection partially explained the final stages of SARS-CoV-2 mutation. In conclusion, model-based hypothesis testing offers a powerful yet underutilized approach to uncovering drivers of viral mutation and gaining deeper insights into pathogen evolution during outbreaks. 
    more » « less
  4. Summary We report a detailed analysis of the full-length SARS-CoV-2 spike dynamics within a native-like membrane environment and variants inaccessible to studies on soluble constructs by conducting hydrogen-deuterium exchange mass spectrometry (HDX-MS) on enveloped virus-like particles (eVLPs) displaying various spike constructs. We find that the previously identified open-interface trimer conformation is sampled in all eVLP-displayed spike variants studied including sequences from engineered vaccine constructs and native viral sequences. The D614G mutation, which arose early in the pandemic, favors the canonical ‘closed-interface’ prefusion conformation, potentially mitigating premature S1 shedding in the presence of a cleaved furin site and providing an evolutionary advantage to the virus. Remarkably, furin cleavage at the S1/S2 boundary allosterically increases the flexibility of the S2’ site, which may facilitate increased TMPRSS2 processing, enhancing viral infectivity. The use of eVLPs in HDX-MS studies provides a powerful platform for studying viral and membrane proteins in near-native environments. 
    more » « less
  5. The ongoing COVID-19 pandemic continues to infect people worldwide, and the virus continues to evolve in significant ways which can pose challenges to the efficiency of available vaccines and therapeutic drugs and cause future pandemic. Therefore, it is important to investigate the binding and interaction of ACE2 with different RBD variants. A comparative study using all-atom MD simulations was conducted on ACE2 binding with 8 different RBD variants, including N501Y, E484K, P479S, T478I, S477N, N439K, K417N and N501YE484K- K417N on RBD. Based on the RMSD, RMSF, and DSSP results, overall the binding of RBD variants with ACE2 is stable, and the secondary structure of RBD and ACE2 are consistent after the point mutation. Besides that, a similar buried surface area, a consistent binding interface and a similar amount of hydrogen bonds formed between RBD and ACE2 although the exact residue pairs on the binding interface were modified. The change of binding free energy from point mutation was predicted using the free energy perturbation (FEP) method. It is found that N501Y, N439K, and K417N can strengthen the binding of RBD with ACE2, while E484K and P479S weaken the binding, and S477N and T478I have negligible effect on the binding. Point mutations modified the dynamic correlation of residues in RBD based on the dihedral angle covariance matrix calculation. Doing dynamic network analysis, a common intrinsic network community extending from the tail of RBD to central, then to the binding interface region was found, which could communicate the dynamics in the binding interface region to the tail thus to the other sections of S protein. The result can supply unique methodology and molecular insight on studying the molecular structure and dynamics of possible future pandemics and design novel drugs. 
    more » « less