Modeling Suggests a Mechanism of Synergy Between Hepatitis C Virus Entry Inhibitors and Drugs of Other Classes

Delay induced stability switch in a mathematical model of CD8 T-cell response to SARS-CoV-2 mediated by receptor ACE2

The pathogen SARS-CoV-2 binds to the receptor angiotensin-converting enzyme 2 (ACE2) of the target cells and then replicates itself through the host, eventually releasing free virus particles. After infection, the CD8 T-cell response is triggered and appears to play a critical role in the defense against virus infections. Infected cells and their activated CD8 T-cells can cause tissue damage. Here, we established a mathematical model of within-host SARS-CoV-2 infection that incorporates the receptor ACE2, the CD8 T-cell response, and the damaged tissues. According to this model, we can get the basic reproduction number R0 and the immune reproduction number R1. We provide the theoretical proof for the stability of the disease-free equilibrium, immune-inactivated equilibrium, and immune-activated equilibrium. Finally, our numerical simulations show that the time delay in CD8 T-cell production can induce complex dynamics such as stability switching. These results provide insights into the mechanisms of SARS-CoV-2 infection and may help in the development of effective drugs against COVID-19.

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

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

Modelling how increased Cathepsin B/L and decreased TMPRSS2 usage for cell entry by the SARS-CoV-2 Omicron variant may affect the efficacy and synergy of TMPRSS2 and Cathepsin B/L inhibitors

The SARS-CoV-2 Omicron variant harbours many mutations in its spike protein compared to the original SARS-CoV-2 strain, which may alter its ability to enter cells, cell tropism, and response to interventions blocking virus entry. To elucidate these effects, we developed a mathematical model of SARS-CoV-2 entry into target cells and applied it to analyse recent in vitro data. SARS-CoV-2 can enter cells via two pathways, one using the host proteases Cathepsin B/L and the other using the host protease TMPRSS2. We found enhanced entry efficiency of the Omicron variant in cells where the original strain preferentially used Cathepsin B/L and reduced efficiency where it used TMPRSS2. The Omicron variant thus appears to have evolved to use the Cathepsin B/L pathway better but at the expense of its ability to use the TMPRSS2 pathway compared to the original strain. We estimated >4-fold enhanced efficiency of the Omicron variant in entry via the Cathepsin B/L pathway and >3-fold reduced efficiency via the TMPRSS2 pathway compared to the original or other strains in a cell type-dependent manner. Our model predicted that Cathepsin B/L inhibitors would be more efficacious and TMPRSS2 inhibitors less efficacious in blocking Omicron variant entry into cells than the original strain. Furthermore, model predictions suggested that drugs simultaneously targeting the two pathways would exhibit synergy. The maximum synergy and drug concentrations yielding it would differ for the Omicron variant compared to the original strain. Our findings provide insights into the cell entry mechanisms of the Omicron variant and have implications for intervention targeting these mechanisms.

Global asymptotic stability of a delay differential equation model for SARS-CoV-2 virus infection mediated by ACE2 receptor protein

The COVID-19 pandemic has brought a serious threat to human life safety worldwide. SARS-CoV-2 virus mainly binds to the target cell surface receptor ACE2 (Angiotensin-converting enzyme 2 ) through the S protein expressed on the surface of the virus, resulting in infection of target cells. During this infection process, the target cell ACE2 receptor plays a very important mediating role. In this paper, a delay differential equation model containing the mediated effect of target cell receptor is established according to the mechanism of SARS-CoV-2 virus invasion of target cells, and the global stability of the infection-free equilibrium and the infected equilibrium of the model is obtained by using the basic reproduction number  ℛ 0  and constructing the appropriate Lyapunov functional. The expression of the basic reproduction number  ℛ 0  intuitively gives the dependence on the expression ratio of the target cell surface ACE2 receptor, which is helpful for the understanding of the mechanism of SARS-CoV-2 virus infection.

Mechanistic Models of COVID-19: Insights into Disease Progression, Vaccines, and Therapeutics

The COVID-19 pandemic has severely impacted health systems and economies worldwide. Significant global efforts are therefore ongoing to improve vaccine efficacies, optimize vaccine deployment, and develop new antiviral therapies to combat the pandemic. Mechanistic viral dynamics and quantitative systems pharmacology models of SARS-CoV-2 infection, vaccines, immunomodulatory agents, and antiviral therapeutics have played a key role in advancing our understanding of SARS-CoV-2 pathogenesis and transmission, the interplay between innate and adaptive immunity to influence the outcomes of infection, effectiveness of treatments, mechanisms and performance of COVID-19 vaccines, and the impact of emerging SARS-CoV-2 variants. Here, we review some of the critical insights provided by these models and discuss the challenges ahead.

Modeling how antibody responses may determine the efficacy of COVID-19 vaccines

Predicting the efficacy of COVID-19 vaccines would aid vaccine development and usage strategies, which is of importance given their limited supplies. Here we develop a multiscale mathematical model that proposes mechanistic links between COVID-19 vaccine efficacies and the neutralizing antibody (NAb) responses they elicit. We hypothesized that the collection of all NAbs would constitute a shape space and that responses of individuals are random samples from this space. We constructed the shape space by analyzing reported in vitro dose-response curves of ~80 NAbs. Sampling NAb subsets from the space, we recapitulated the responses of convalescent patients. We assumed that vaccination would elicit similar NAb responses. We developed a model of within-host SARS-CoV-2 dynamics, applied it to virtual patient populations and, invoking the NAb responses above, predicted vaccine efficacies. Our predictions quantitatively captured the efficacies from clinical trials. Our study thus suggests plausible mechanistic underpinnings of COVID-19 vaccines and generates testable hypotheses for establishing them.

Life cycle process dependencies of positive-sense RNA viruses suggest strategies for inhibiting productive cellular infection

Life cycle processes of positive-strand (+)RNA viruses are broadly conserved across families, yet they employ different strategies to grow in the cell. Using a generalized dynamical model for intracellular (+)RNA virus growth, we decipher these life cycle determinants and their dependencies for several viruses and parse the effects of viral mutations, drugs and host cell permissivity. We show that poliovirus employs rapid replication and virus assembly, whereas the Japanese encephalitis virus leverages its higher rate of translation and efficient cellular reorganization compared to the hepatitis C virus. Stochastic simulations demonstrate infection extinction if all seeding (inoculating) viral RNA degrade before establishing robust replication critical for infection. The probability of this productive cellular infection, 'cellular infectivity', is affected by virus-host processes and defined by early life cycle events and viral seeding. An increase in cytoplasmic RNA degradation and delay in vesicular compartment formation reduces infectivity, more so when combined. Synergy among these parameters in limiting (+)RNA virus infection as predicted by our model suggests new avenues for inhibiting infections by targeting the early life cycle bottlenecks.

Flavonoid-triazolyl hybrids as potential anti-hepatitis C virus agents: Synthesis and biological evaluation

A series of flavonoid-triazolyl hybrids were synthesized and evaluated as novel inhibitors of hepatitis C virus (HCV). The results of anti-HCV activity assays showed that most of the synthesized derivatives at a concentration of 100 μg/mL inhibited the generation of progeny virus. Among these derivatives, 10m and 10r exhibited the most potent anti-HCV activity and inhibited the production of HCV in a dose-dependent manner. Interestingly, 10m and 10r had no significant inhibitory effect on viral translation or replication. Additional action mechanism studies revealed that the most potent compounds, 10m and 10r, significantly inhibited viral entry to 34.0% and 52.0%, respectively, at 10 μM. These results suggest further effective application of 10m and 10r as potential HCV preventive agents.

Targeting TMPRSS2 and Cathepsin B/L together may be synergistic against SARS-CoV-2 infection

The entry of SARS-CoV-2 into target cells requires the activation of its surface spike protein, S, by host proteases. The host serine protease TMPRSS2 and cysteine proteases Cathepsin B/L can activate S, making two independent entry pathways accessible to SARS-CoV-2. Blocking the proteases prevents SARS-CoV-2 entry in vitro. This blockade may be achieved in vivo through 'repurposing' drugs, a potential treatment option for COVID-19 that is now in clinical trials. Here, we found, surprisingly, that drugs targeting the two pathways, although independent, could display strong synergy in blocking virus entry. We predicted this synergy first using a mathematical model of SARS-CoV-2 entry and dynamics in vitro. The model considered the two pathways explicitly, let the entry efficiency through a pathway depend on the corresponding protease expression level, which varied across cells, and let inhibitors compromise the efficiency in a dose-dependent manner. The synergy predicted was novel and arose from effects of the drugs at both the single cell and the cell population levels. Validating our predictions, available in vitro data on SARS-CoV-2 and SARS-CoV entry displayed this synergy. Further, analysing the data using our model, we estimated the relative usage of the two pathways and found it to vary widely across cell lines, suggesting that targeting both pathways in vivo may be important and synergistic given the broad tissue tropism of SARS-CoV-2. Our findings provide insights into SARS-CoV-2 entry into target cells and may help improve the deployability of drug combinations targeting host proteases required for the entry.

You Cannot Have Your Synergy and Efficacy Too

Synergistic drugs are preferred in combination therapies for many diseases, including viral infections and cancers. Maximizing synergy, however, may come at the cost of efficacy. This synergy-efficacy trade-off appears to be widely prevalent and independent of the specific drug interactions yielding synergy. We present examples of the trade-off in drug combinations used in hepatitis C, HIV, and cancer therapies and believe that screens for optimal drug combinations that presently seek to maximize synergy may be improved by considering the trade-off.

Interferon at the cellular, individual, and population level in hepatitis C virus infection: Its role in the interferon-free treatment era

The advent of powerful direct-acting antiviral agents (DAAs) has revolutionized the treatment of hepatitis C. DAAs cure nearly all patients with short duration, oral treatments. Significant efforts are now underway to optimize DAA-based treatments. We discuss the potential role of interferon in this optimization. Clinical studies present compelling evidence that DAAs perform better in treatment-naive individuals than in individuals who previously failed treatment with interferon, a surprising correlation because interferon and DAAs are thought to act independently. Recent mathematical models explore a mechanistic hypothesis underlying this correlation. The hypothesis invokes the action of interferon at the cellular, individual, and population levels. Strong interferon responses prevent the productive infection of cells, reduce viral replication, and impede the development of resistance to DAAs in infected individuals and improve cure rates elicited by DAAs in treated populations. The models develop descriptions of these processes, integrate them into a comprehensive framework, and capture clinical data quantitatively, providing a successful test of the hypothesis. Individuals with strong endogenous interferon responses thus present a promising subpopulation for reducing DAA treatment durations. This review discusses the conceptual advances made by the models, highlights the new insights they unravel, and examines their applicability to optimize DAA-based treatments.

Modelling how responsiveness to interferon improves interferon-free treatment of hepatitis C virus infection

Direct-acting antiviral agents (DAAs) for hepatitis C treatment tend to fare better in individuals who are also likely to respond well to interferon-alpha (IFN), a surprising correlation given that DAAs target specific viral proteins whereas IFN triggers a generic antiviral immune response. Here, we posit a causal relationship between IFN-responsiveness and DAA treatment outcome. IFN-responsiveness restricts viral replication, which would prevent the growth of viral variants resistant to DAAs and improve treatment outcome. To test this hypothesis, we developed a multiscale mathematical model integrating IFN-responsiveness at the cellular level, viral kinetics and evolution leading to drug resistance at the individual level, and treatment outcome at the population level. Model predictions quantitatively captured data from over 50 clinical trials demonstrating poorer response to DAAs in previous non-responders to IFN than treatment-naïve individuals, presenting strong evidence supporting the hypothesis. Model predictions additionally described several unexplained clinical observations, viz., the percentages of infected individuals who 1) spontaneously clear HCV, 2) get chronically infected but respond to IFN-based therapy, and 3) fail IFN-based therapy but respond to DAA-based therapy, resulting in a comprehensive understanding of HCV infection and treatment. An implication of the causal relationship is that failure of DAA-based treatments may be averted by adding IFN, a strategy of potential use in settings with limited access to DAAs. A second, wider implication is that individuals with greater IFN-responsiveness would require shorter DAA-based treatment durations, presenting a basis and a promising population for response-guided therapy.