skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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. 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
  4. The rapid spread of SARS-CoV-2 required immediate actions to control the transmission of the virus and minimize its impact on humanity. An extensive mutation rate of this viral genome contributes to the virus’ ability to quickly adapt to environmental changes, impacts transmissibility and antigenicity, and may facilitate immune escape. Therefore, it is of great interest for researchers working in vaccine development and drug design to consider the impact of mutations on virus-drug interactions. Here, we propose a multitarget drug discovery pipeline for identifying potential drug candidates which can efficiently inhibit the Receptor Binding Domain (RBD) of spike glycoproteins from different variants of SARS-CoV-2. Eight homology models of RBDs for selected variants were created and validated using reference crystal structures. We then investigated interactions between host receptor ACE2 and RBDs from nine variants of SARS-CoV-2. It led us to conclude that efficient multi-variant targeting drugs should be capable of blocking residues Q(R)493 and N487 in RBDs. Using methods of molecular docking, molecular mechanics, and molecular dynamics, we identified three lead compounds (hesperidin, narirutin, and neohesperidin) suitable for multitarget SARS-CoV-2 inhibition. These compounds are flavanone glycosides found in citrus fruits – an active ingredient of Traditional Chinese Medicines. The developed pipeline can be further used to (1) model mutants for which crystal structures are not yet available and (2) scan a more extensive library of compounds against other mutated viral proteins. 
    more » « less
  5. We present methods for building a Java Runtime-Alterable-Model Platform (RAMP) of complex dynamical systems. We illustrate our methods by building a multivariant SEIR (epidemic) RAMP. Underlying our RAMP is an individual-based model that includes adaptive contact rates, pathogen genetic drift, waning and cross-immunity. Besides allowing parameter values, process descriptions and scriptable runtime drivers to be easily modified during simulations, our RAMP can used within R-Studio and other computational platforms. Process descriptions that can be runtime altered within our SEIR RAMP include pathogen variant-dependent host shedding, environmental persistence, host transmission and within-host pathogen mutation and replication. They also include adaptive social distancing and adaptive application of vaccination rates and variant-valency of vaccines. We present simulation results using parameter values and process descriptions relevant to the current COVID-19 pandemic. Our results suggest that if waning immunity outpaces vaccination rates, then vaccination rollouts may fail to contain the most transmissible variants, particularly if vaccine valencies are not adapted to deal with escape mutations. Our SEIR RAMP is designed for easy use by others. More generally, our RAMP concept facilitates construction of highly flexible complex systems models of all types, which can then be easily shared as stand-alone application programs. 
    more » « less