The role of cells refractory to productive infection in acute hepatitis B viral dynamics

Application of MIDD to accelerate the development of anti-infectives: Current status and future perspectives

This review examines the role of model-informed drug development (MIDD) in advancing antibacterial and antiviral drug development, with an emphasis on the inclusion of host system dynamics into modeling efforts. Amidst the growing challenges of multidrug resistance and diminishing market returns, innovative methodologies are crucial for continuous drug discovery and development. The MIDD approach, with its robust capacity to integrate diverse data types, offers a promising solution. In particular, the utilization of appropriate modeling and simulation techniques for better characterization and early assessment of drug resistance are discussed. The evolution of MIDD practices across different infectious disease fields is also summarized, and compared to advancements achieved in oncology. Moving forward, the application of MIDD should expand into host system dynamics as these considerations are critical for the development of "live drugs" (e.g. chimeric antigen receptor T cells or bacteriophages) to address issues like antibiotic resistance or latent viral infections.

A multiscale model of the action of a capsid assembly modulator for the treatment of chronic hepatitis B

Chronic hepatitis B virus (HBV) infection is strongly associated with increased risk of liver cancer and cirrhosis. While existing treatments effectively inhibit the HBV life cycle, viral rebound occurs rapidly following treatment interruption. Consequently, functional cure rates of chronic HBV infection remain low and there is increased interest in a novel treatment modality, capsid assembly modulators (CAMs). Here, we develop a multiscale mathematical model of CAM treatment in chronic HBV infection. By fitting the model to participant data from a phase I trial of the first-generation CAM vebicorvir, we estimate the drug's dose-dependent effectiveness and identify the physiological mechanisms that drive the observed biphasic decline in HBV DNA and RNA, and mechanistic differences between HBeAg-positive and negative infection. Finally, we demonstrate analytically and numerically that HBV RNA is more sensitive than HBV DNA to increases in CAM effectiveness.

Incorporating Intracellular Processes in Virus Dynamics Models

In-host models have been essential for understanding the dynamics of virus infection inside an infected individual. When used together with biological data, they provide insight into viral life cycle, intracellular and cellular virus-host interactions, and the role, efficacy, and mode of action of therapeutics. In this review, we present the standard model of virus dynamics and highlight situations where added model complexity accounting for intracellular processes is needed. We present several examples from acute and chronic viral infections where such inclusion in explicit and implicit manner has led to improvement in parameter estimates, unification of conclusions, guidance for targeted therapeutics, and crossover among model systems. We also discuss trade-offs between model realism and predictive power and highlight the need of increased data collection at finer scale of resolution to better validate complex models.

Mathematical Models of Early Hepatitis B Virus Dynamics in Humanized Mice

Analyzing the impact of the adaptive immune response during acute hepatitis B virus (HBV) infection is essential for understanding disease progression and control. Here we developed mathematical models of HBV infection which either lack terms for adaptive immune responses, or assume adaptive immune responses in the form of cytolytic immune killing, non-cytolytic immune cure, or non-cytolytic-mediated block of viral production. We validated the model that does not include immune responses against temporal serum hepatitis B DNA (sHBV) and temporal serum hepatitis B surface-antigen (HBsAg) experimental data from mice engrafted with human hepatocytes (HEP). Moreover, we validated the immune models against sHBV and HBsAg experimental data from mice engrafted with HEP and human immune system (HEP/HIS). As expected, the model that does not include adaptive immune responses matches the observed high sHBV and HBsAg concentrations in all HEP mice. By contrast, while all immune response models predict reduction in sHBV and HBsAg concentrations in HEP/HIS mice, the Akaike Information Criterion cannot discriminate between non-cytolytic cure (resulting in a class of cells refractory to reinfection) and antiviral block functions (of up to 99 % viral production 1-3 weeks following peak viral load). We can, however, reject cytolytic killing, as it can only match the sHBV and HBsAg data when we predict unrealistic levels of hepatocyte loss.

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.

Multiscale modeling of HBV infection integrating intra- and intercellular viral propagation to analyze extracellular viral markers

Chronic infection with hepatitis B virus (HBV) is caused by the persistence of closed circular DNA (cccDNA) in the nucleus of infected hepatocytes. Despite available therapeutic anti-HBV agents, eliminating the cccDNA remains challenging. Thus, quantifying and understanding the dynamics of cccDNA are essential for developing effective treatment strategies and new drugs. However, such study requires repeated liver biopsy to measure the intrahepatic cccDNA, which is basically not accepted because liver biopsy is potentially morbid and not common during hepatitis B treatment. We here aimed to develop a noninvasive method for quantifying cccDNA in the liver using surrogate markers in peripheral blood. We constructed a multiscale mathematical model that explicitly incorporates both intracellular and intercellular HBV infection processes. The model, based on age-structured partial differential equations, integrates experimental data from in vitro and in vivo investigations. By applying this model, we roughly predicted the amount and dynamics of intrahepatic cccDNA within a certain range using specific viral markers in serum samples, including HBV DNA, HBsAg, HBeAg, and HBcrAg. Our study represents a significant step towards advancing the understanding of chronic HBV infection. The noninvasive quantification of cccDNA using our proposed method holds promise for improving clinical analyses and treatment strategies. By comprehensively describing the interactions of all components involved in HBV infection, our multiscale mathematical model provides a valuable framework for further research and the development of targeted interventions.

Global stability of secondary DENV infection models with non-specific and strain-specific CTLs

Dengue virus (DENV) is a highly perilous virus that is transmitted to humans through mosquito bites and causes dengue fever. Consequently, extensive efforts are being made to develop effective treatments and vaccines. Mathematical modeling plays a significant role in comprehending the dynamics of DENV within a host in the presence of cytotoxic T lymphocytes (CTL) immune response. This study examines two models for secondary DENV infections that elucidate the dynamics of DENV under the influence of two types of CTL responses, namely non-specific and strain-specific responses. The first model encompasses five compartments, which consist of uninfected monocytes, infected monocytes, free DENV particles, non-specific CTLs, and strain-specific CTLs. In the second model, latently infected cells are introduced into the model. We posit that the CTL responsiveness is determined by a combination of self-regulating CTL response and a predator-prey-like CTL response. The model's solutions are verified to be nonnegativity and bounded and the model possesses two equilibrium states: the uninfected equilibrium EQ 0 and the infected equilibrium EQ ⁎ . Furthermore, we calculate the basic reproduction number R 0 , which determines the existence and stability of the model's equilibria. We examine the global stability by constructing suitable Lyapunov functions. Our analysis reveals that if R 0 ≤ 1 , then EQ 0 is globally asymptotically stable (G.A.S), and if R 0 > 1 , then EQ 0 is unstable while EQ ⁎ is G.A.S. To illustrate our findings analytically, we conduct numerical simulations for each model. Additionally, we perform sensitivity analysis to demonstrate how the parameter values of the proposed model impact R 0 given a set of data. Finally, we discuss the implications of including the CTL immune response and latently infected cells in the secondary DENV infection model. Our study demonstrates that incorporating the CTL immune response and latently infected cells diminishes R 0 and enhances the system's stability around EQ 0 .

Multiscale modeling of HBV infection integrating intra- and intercellular viral propagation for analyzing extracellular viral markers

Chronic infection of hepatitis B virus (HBV) is caused by the persistence of closed circular DNA (cccDNA) in the nucleus of infected hepatocytes. Despite available therapeutic anti-HBV agents, eliminating the cccDNA remains challenging. The quantifying and understanding dynamics of cccDNA are essential for developing effective treatment strategies and new drugs. However, it requires a liver biopsy to measure the intrahepatic cccDNA, which is basically not accepted because of the ethical aspect. We here aimed to develop a non-invasive method for quantifying cccDNA in the liver using surrogate markers present in peripheral blood. We constructed a multiscale mathematical model that explicitly incorporates both intracellular and intercellular HBV infection processes. The model, based on age-structured partial differential equations (PDEs), integrates experimental data from in vitro and in vivo investigations. By applying this model, we successfully predicted the amount and dynamics of intrahepatic cccDNA using specific viral markers in serum samples, including HBV DNA, HBsAg, HBeAg, and HBcrAg. Our study represents a significant step towards advancing the understanding of chronic HBV infection. The non-invasive quantification of cccDNA using our proposed methodology holds promise for improving clinical analyses and treatment strategies. By comprehensively describing the interactions of all components involved in HBV infection, our multiscale mathematical model provides a valuable framework for further research and the development of targeted interventions.

Global dynamics of multiple delays within-host model for a hepatitis B virus infection of hepatocytes with immune response and drug therapy

In this paper, a mathematical model describing the hepatitis B virus (HBV) infection of hepatocytes with the intracellular HBV-DNA containing capsids, cytotoxic T-lymphocyte (CTL), antibodies including drug therapy (blocking new infection and inhibiting viral production) with two-time delays is studied. It incorporates the delay in the productively infected hepatocytes and the delay in an antigenic stimulation generating CTL. We verify the positivity and boundedness of solutions and determine the basic reproduction number. The local and global stability of three equilibrium points (infection-free, immune-free, and immune-activated) are investigated. Finally, the numerical simulations are established to show the role of these therapies in reducing viral replication and HBV infection. Our results show that the treatment by blocking new infection gives more significant results than the treatment by inhibiting viral production for infected hepatocytes. Further, both delays affect the number of infections and duration i.e. the longer the delay, the more severe the HBV infection.

Advances in Parameter Estimation and Learning from Data for Mathematical Models of Hepatitis C Viral Kinetics

Mathematical models, some of which incorporate both intracellular and extracellular hepatitis C viral kinetics, have been advanced in recent years for studying HCV–host dynamics, antivirals mode of action, and their efficacy. The standard ordinary differential equation (ODE) hepatitis C virus (HCV) kinetic model keeps track of uninfected cells, infected cells, and free virus. In multiscale models, a fourth partial differential equation (PDE) accounts for the intracellular viral RNA (vRNA) kinetics in an infected cell. The PDE multiscale model is substantially more difficult to solve compared to the standard ODE model, with governing differential equations that are stiff. In previous contributions, we developed and implemented stable and efficient numerical methods for the multiscale model for both the solution of the model equations and parameter estimation. In this contribution, we perform sensitivity analysis on model parameters to gain insight into important properties and to ensure our numerical methods can be safely used for HCV viral dynamic simulations. Furthermore, we generate in-silico patients using the multiscale models to perform machine learning from the data, which enables us to remove HCV measurements on certain days and still be able to estimate meaningful observations with a sufficiently small error.

An age-structured model of hepatitis B viral infection highlights the potential of different therapeutic strategies

Hepatitis B virus (HBV) is a global health threat, and its elimination by 2030 has been prioritised by the World Health Organisation. Here we present an age-structured model for the immune response to an HBV infection, which takes into account contributions from both cell-mediated and humoral immunity. The model has been validated using published patient data recorded during acute infection. It has been adapted to the scenarios of chronic infection, clearance of infection, and flare-ups via variation of the immune response parameters. The impacts of immune response exhaustion and non-infectious subviral particles on the immune response dynamics are analysed. A comparison of different treatment options in the context of this model reveals that drugs targeting aspects of the viral life cycle are more effective than exhaustion therapy, a form of therapy mitigating immune response exhaustion. Our results suggest that antiviral treatment is best started when viral load is declining rather than in a flare-up. The model suggests that a fast antibody production rate always leads to viral clearance, highlighting the promise of antibody therapies currently in clinical trials.

Modeling the dynamics of Usutu virus infection in birds

Usutu virus is an emerging zoonotic flavivirus causing high avian mortality rates and occasional severe neurological disorders in humans. Several virus strains are co-circulating and the differences in their characteristics and avian pathogenesis levels are still unknown. In this study, we use within-host mathematical models to characterize the mechanisms responsible for virus expansion and clearance in juvenile chickens challenged with four Usutu virus strains. We find heterogeneity between the virus strains, with the time between cell infection and viral production varying between 16 h and 23 h, the infected cell lifespan varying between 48 min and 9.5 h, and the basic reproductive number R₀ varying between 12.05 and 19.49. The strains with high basic reproductive number have short infected cell lifespan, indicative of immune responses. The virus strains with low basic reproductive number have lower viral peaks and longer lasting viremia, due to lower infection rates and high infected cell lifespan. We discuss how the host and virus heterogeneities may differently impact the public health threat presented by these virus strains.

Towards a Mathematical Model for the Viral Progression in the Pharynx

In this work, a comprehensive model for the viral progression in the pharynx has been developed. This one-dimension model considers both Fickian diffusion and convective flow coupled with chemical reactions, such as virus population growth, infected and uninfected cell accumulation as well as virus clearance. The effect of a sterilizing agent such as an alcoholic solution on the viral progression in the pharynx was taken into account and a parametric analysis for the effect of kinetic rate parameters on virus propagation was made. Moreover, different conditions caused by further medical treatment, such as a decrease in virus yield per infected cell, were examined. It is shown that the infection fails to establish by decreasing the virus yield per infected cell. It is believed that this work could be used to further investigate the medical treatment of viral progression in the pharynx.

A mathematical model to evaluate the role of memory B and T cells in heterologous secondary dengue infection

We propose a mathematical model to investigate the antibody-dependent enhancement (ADE) phenomenon during secondary dengue infection. The model consists of an ODE system that describes the interaction of the dengue virus with macrophages and memory B and T cell role during the infection. The qualitative model analysis is done in terms of memory B and T cell cloning parameters and the basic reproduction number R₀. In the absence of memory B and T cell cloning, if R₀<1 the dengue virus population extinguishes, while for R₀>1, it tends asymptotically to a positive equilibrium. However, when we consider the memory B cell cloning, it is possible to occur dengue infection even when R₀<1. Memory T cells have an essential role in eliminating the possibility of ADE occurrence when R₀<1.

A quantitative systems pharmacology model for acute viral hepatitis B

•

Mechanistic model characterizing acute immune response and HBV system interactions.

•

Key role of the cellular and regulatory response triggering hepatitis B chronicity.

•

Modelling framework to easily incorporate and explore additional biological mechanisms.

Hepatitis B liver infection is caused by hepatitis B virus (HBV) and represents a major global disease problem when it becomes chronic, as is the case for 80–90% of vertical or early life infections. However, in the vast majority (>95%) of adult exposures, the infected individuals are capable of mounting an effective immune response leading to infection resolution. A good understanding of HBV dynamics and the interaction between the virus and immune system during acute infection represents an essential step to characterize and understand the key biological processes involved in disease resolution, which may help to identify potential interventions to prevent chronic hepatitis B.

In this work, a quantitative systems pharmacology model for acute hepatitis B characterizing viral dynamics and the main components of the innate, adaptive, and tolerant immune response has been successfully developed. To do so, information from multiple sources and across different organization levels has been integrated in a common mechanistic framework. The final model adequately describes the chronology and plausibility of an HBV-triggered immune response, as well as clinical data from acute patients reported in the literature. Given the holistic nature of the framework, the model can be used to illustrate the relevance of the different immune pathways and biological processes to ultimate response, observing the negligible contribution of the innate response and the key contribution of the cellular response on viral clearance. More specifically, moderate reductions of the proliferation of activated cytotoxic CD8+ lymphocytes or increased immunoregulatory effects can drive the system towards chronicity.

Modeling Within-Host Dynamics of SARS-CoV-2 Infection: A Case Study in Ferrets

The pre-clinical development of antiviral agents involves experimental trials in animals and ferrets as an animal model for the study of SARS-CoV-2. Here, we used mathematical models and experimental data to characterize the within-host infection dynamics of SARS-CoV-2 in ferrets. We also performed a global sensitivity analysis of model parameters impacting the characteristics of the viral infection. We provide estimates of the viral dynamic parameters in ferrets, such as the infection rate, the virus production rate, the infectious virus proportion, the infected cell death rate, the virus clearance rate, as well as other related characteristics, including the basic reproduction number, pre-peak infectious viral growth rate, post-peak infectious viral decay rate, pre-peak infectious viral doubling time, post-peak infectious virus half-life, and the target cell loss in the respiratory tract. These parameters and indices are not significantly different between animals infected with viral strains isolated from the environment and isolated from human hosts, indicating a potential for transmission from fomites. While the infection period in ferrets is relatively short, the similarity observed between our results and previous results in humans supports that ferrets can be an appropriate animal model for SARS-CoV-2 dynamics-related studies, and our estimates provide helpful information for such studies.

Comparing antiviral strategies against COVID-19 via multiscale within-host modelling

Within-host models of COVID-19 infection dynamics enable the merits of different forms of antiviral therapy to be assessed in individual patients. A stochastic agent-based model of COVID-19 intracellular dynamics is introduced here, that incorporates essential steps of the viral life cycle targeted by treatment options. Integration of model predictions with an intercellular ODE model of within-host infection dynamics, fitted to patient data, generates a generic profile of disease progression in patients that have recovered in the absence of treatment. This is contrasted with the profiles obtained after variation of model parameters pertinent to the immune response, such as effector cell and antibody proliferation rates, mimicking disease progression in immunocompromised patients. These profiles are then compared with disease progression in the presence of antiviral and convalescent plasma therapy against COVID-19 infections. The model reveals that using both therapies in combination can be very effective in reducing the length of infection, but these synergistic effects decline with a delayed treatment start. Conversely, early treatment with either therapy alone can actually increase the duration of infection, with infectious virions still present after the decline of other markers of infection. This suggests that usage of these treatments should remain carefully controlled in a clinical environment.

Early events in hepatitis B infection: the role of inoculum dose

The relationship between the inoculum dose and the ability of the pathogen to invade the host is poorly understood. Experimental studies in non-human primates infected with different inoculum doses of hepatitis B virus have shown a non-monotonic relationship between dose magnitude and infection outcome, with high and low doses leading to 100% liver infection and intermediate doses leading to less than 0.1% liver infection, corresponding to CD4 T-cell priming. Since hepatitis B clearance is CD8 T-cell mediated, the question of whether the inoculum dose influences CD8 T-cell dynamics arises. To help answer this question, we developed a mathematical model of virus-host interaction following hepatitis B virus infection. Our model explains the experimental data well, and predicts that the inoculum dose affects both the timing of the CD8 T-cell expansion and the quality of its response, especially the non-cytotoxic function. We find that a low-dose challenge leads to slow CD8 T-cell expansion, weak non-cytotoxic functions, and virus persistence; high- and medium-dose challenges lead to fast CD8 T-cell expansion, strong cytotoxic and non-cytotoxic function, and virus clearance; while a super-low-dose challenge leads to delayed CD8 T-cell expansion, strong cytotoxic and non-cytotoxic function, and virus clearance. These results are useful for designing immune cell-based interventions.

Understanding the antiviral effects of RNAi-based therapy in HBeAg-positive chronic hepatitis B infection

The RNA interference (RNAi) drug ARC-520 was shown to be effective in reducing serum hepatitis B virus (HBV) DNA, hepatitis B e antigen (HBeAg) and hepatitis B surface antigen (HBsAg) in HBeAg-positive patients treated with a single dose of ARC-520 and daily nucleosidic analogue (entecavir). To provide insights into HBV dynamics under ARC-520 treatment and its efficacy in blocking HBV DNA, HBsAg, and HBeAg production we developed a multi-compartmental pharmacokinetic-pharamacodynamic model and calibrated it with frequent measured HBV kinetic data. We showed that the time-dependent single dose ARC-520 efficacies in blocking HBsAg and HBeAg are more than 96% effective around day 1, and slowly wane to 50% in 1-4 months. The combined single dose ARC-520 and entecavir effect on HBV DNA was constant over time, with efficacy of more than 99.8%. The observed continuous HBV DNA decline is entecavir mediated, the strong but transient HBsAg and HBeAg decays are ARC-520 mediated. The modeling framework may help assess ongoing RNAi drug development for hepatitis B virus infection.

Stochastic dynamics in a time-delayed model for autoimmunity

In this paper we study interactions between stochasticity and time delays in the dynamics of immune response to viral infections, with particular interest in the onset and development of autoimmune response. Starting with a deterministic time-delayed model of immune response to infection, which includes cytokines and T cells with different activation thresholds, we derive an exact delayed chemical master equation for the probability density. We use system size expansion and linear noise approximation to explore how variance and coherence of stochastic oscillations depend on parameters, and to show that stochastic oscillations become more regular when regulatory T cells become more effective at clearing autoreactive T cells. Reformulating the model as an Itô stochastic delay differential equation, we perform numerical simulations to illustrate the dynamics of the model and associated probability distributions in different parameter regimes. The results suggest that even in cases where the deterministic model has stable steady states, in individual stochastic realisations, the model can exhibit sustained stochastic oscillations, whose variance increases as one gets closer to the deterministic stability boundary. Furthermore, in the regime of bi-stability, whereas deterministically the system would approach one of the steady states (or periodic solutions) depending on the initial conditions, due to the presence of stochasticity, it is now possible for the system to reach both of those dynamical states with certain probability. Biological significance of this result lies in highlighting the fact that since normally in a laboratory or clinical setting one would observe a single individual realisation of the course of the disease, even for all parameters characterising the immune system and the strength of infection being the same, there is a proportion of cases where a spontaneous recovery can be observed, and similarly, where a disease can develop in a situation that otherwise would result in a normal disease clearance.

Within-host mathematical models of hepatitis B virus infection: Past, present, and future

Mathematical modeling has been instrumental in enhancing our understanding of the viral dynamics of hepatitis B virus (HBV) infection. We give a primer on HBV infection in humans and a brief overview of the development of within-host mathematical models of HBV infection. In the last decade, models have advanced from considering chronic HBV infections under therapy to the pathogenesis of infection. We also summarize estimates of key viral dynamic parameters that have varied greatly among studies, and show that they impact model predictions. Future directions for mathematical modeling of HBV infection are proposed to better understand emerging therapies, the HBV life cycle, predicting cure, and the mechanisms involved in the immune response to HBV infection.

Acute hepatitis B virus infection model within the host incorporating immune cells and cytokine responses

We formulate and analyze a within-host hepatitis B viral mathematical model for hepatitis B in the acute phase of infection. The model incorporates hepatocytes, hepatitis B virus, immune system cells and cytokine dynamics using a system of ordinary differential equations. We use the model to demonstrate the trends of the hepatitis B infection qualitatively without the effects of immune cells and cytokines. Using these trends, we tested the effects of incorporating the immune cells only and immune cells with cytokine responses at low and high inhibitions on the hepatitis B virus infection. Our results showed that it is impossible to have the immune cells work independently from cytokines when there is an acute hepatitis B virus infection. Therefore, our results suggest that incorporating immune cells and cytokine dynamics in the acute hepatitis B virus infection stage delays infection in the hepatocytes and excluding such dynamics speeds up infection during this phase. Results from this study are useful in developing strategies for control of hepatocellular carcinoma which is caused by hepatitis B virus infection.

A Parameter Estimation Method for Multiscale Models of Hepatitis C Virus Dynamics

Mathematical models that are based on differential equations require detailed knowledge about the parameters that are included in the equations. Some of the parameters can be measured experimentally while others need to be estimated. When the models become more sophisticated, such as in the case of multiscale models of hepatitis C virus dynamics that deal with partial differential equations (PDEs), several strategies can be tried. It is possible to use parameter estimation on an analytical approximation of the solution to the multiscale model equations, namely the long-term approximation, but this limits the scope of the parameter estimation method used and a long-term approximation needs to be derived for each model. It is possible to transform the PDE multiscale model to a system of ODEs, but this has an effect on the model parameters themselves and the transformation can become problematic for some models. Finally, it is possible to use numerical solutions for the multiscale model and then use canned methods for the parameter estimation, but the latter is making the user dependent on a black box without having full control over the method. The strategy developed here is to start by working directly on the multiscale model equations for preparing them toward the parameter estimation method that is fully coded and controlled by the user. It can also be adapted to multiscale models of other viruses. The new method is described, and illustrations are provided using a user-friendly simulator that incorporates the method.

Modeling the dynamics of hepatitis B infection, immunity, and drug therapy

Hepatitis B virus infection is the cause of liver diseases such as cirrhosis and liver cancer. Understanding the host-virus mechanisms that mediate virus pathogenesis can help design better preventive measures for disease control. Mathematical models have been used alongside experimental data to provide insight into the role of immune responses during the acute and chronic hepatitis B infections as well as virus dynamics following administration of combined drug therapy. In this paper, we review several modeling studies on virus-host interactions during acute infection, the virus-host characteristics responsible for transition to chronic disease, and the efficacy and optimal control measures of drug therapy. We conclude by presenting our opinion on the future directions of the field.

Superinfection and cure of infected cells as mechanisms for hepatitis C virus adaptation and persistence

RNA viruses exist as a genetically diverse quasispecies with extraordinary ability to adapt to abrupt changes in the host environment. However, the molecular mechanisms that contribute to their rapid adaptation and persistence in vivo are not well studied. Here, we probe hepatitis C virus (HCV) persistence by analyzing clinical samples taken from subjects who were treated with a second-generation HCV protease inhibitor. Frequent longitudinal viral load determinations and large-scale single-genome sequence analyses revealed rapid antiviral resistance development, and surprisingly, dynamic turnover of dominant drug-resistant mutant populations long after treatment cessation. We fitted mathematical models to both the viral load and the viral sequencing data, and the results provided strong support for the critical roles that superinfection and cure of infected cells play in facilitating the rapid turnover and persistence of viral populations. More broadly, our results highlight the importance of considering viral dynamics and competition at the intracellular level in understanding rapid viral adaptation. Thus, we propose a theoretical framework integrating viral and molecular mechanisms to explain rapid viral evolution, resistance, and persistence despite antiviral treatment and host immune responses.

Numerical schemes for solving and optimizing multiscale models with age of hepatitis C virus dynamics

Age-structured PDE models have been developed to study viral infection and treatment. However, they are notoriously difficult to solve. Here, we investigate the numerical solutions of an age-based multiscale model of hepatitis C virus (HCV) dynamics during antiviral therapy and compare them with an analytical approximation, namely its long-term approximation. First, starting from a simple yet flexible numerical solution that also considers an integral approximated over previous iterations, we show that the long-term approximation is an underestimate of the PDE model solution as expected since some infection events are being ignored. We then argue for the importance of having a numerical solution that takes into account previous iterations for the associated integral, making problematic the use of canned solvers. Second, we demonstrate that the governing differential equations are stiff and the stability of the numerical scheme should be considered. Third, we show that considerable gain in efficiency can be achieved by using adaptive stepsize methods over fixed stepsize methods for simulating realistic scenarios when solving multiscale models numerically. Finally, we compare between several numerical schemes for the solution of the equations and demonstrate the use of a numerical optimization scheme for the parameter estimation performed directly from the equations.

Mathematical model of immune response to hepatitis B

A new detailed mathematical model for dynamics of immune response to hepatitis B is proposed, which takes into account contributions from innate and adaptive immune responses, as well as cytokines. Stability analysis of different steady states is performed to identify parameter regions where the model exhibits clearance of infection, maintenance of a chronic infection, or periodic oscillations. Effects of nucleoside analogues and interferon treatments are analysed, and the critical drug efficiency is determined.

A multi-scale spatial model of hepatitis-B viral dynamics

Chronic hepatitis B viral infection (HBV) afflicts around 250 million individuals globally and few options for treatment exist. Once infected, the virus entrenches itself in the liver with a notoriously resilient colonisation of viral DNA (covalently-closed circular DNA, cccDNA). The majority of infections are cleared, yet we do not understand why 5% of adult immune responses fail leading to the chronic state with its collateral morbid effects such as cirrhosis and eventual hepatic carcinoma. The liver environment exhibits particularly complex spatial structures for metabolic processing and corresponding distributions of nutrients and transporters that may influence successful HBV entrenchment. We assembled a multi-scaled mathematical model of the fundamental hepatic processing unit, the sinusoid, into a whole-liver representation to investigate the impact of this intrinsic spatial heterogeneity on the HBV dynamic. Our results suggest HBV may be exploiting spatial aspects of the liver environment. We distributed increased HBV replication rates coincident with elevated levels of nutrients in the sinusoid entry point (the periportal region) in tandem with similar distributions of hepatocyte transporters key to HBV invasion (e.g., the sodium-taurocholate cotransporting polypeptide or NTCP), or immune system activity. According to our results, such co-alignment of spatial distributions may contribute to persistence of HBV infections, depending on spatial distributions and intensity of immune response as well. Moreover, inspired by previous HBV models and experimentalist suggestions of extra-hepatic HBV replication, we tested in our model influence of HBV blood replication and observe an overall nominal effect on persistent liver infection. Regardless, we confirm prior results showing a solo cccDNA is sufficient to re-infect an entire liver, with corresponding concerns for transplantation and treatment.

The Role of Infected Cell Proliferation in the Clearance of Acute HBV Infection in Humans

Around 90-95% of hepatitis B virus (HBV) infected adults do not progress to the chronic phase and, instead, recover naturally. The strengths of the cytolytic and non-cytolytic immune responses are key players that decide the fate of acute HBV infection. In addition, it has been hypothesized that proliferation of infected cells resulting in uninfected progeny and/or cytokine-mediated degradation of covalently closed circular DNA (cccDNA) leading to the cure of infected cells are two major mechanisms assisting the adaptive immune response in the clearance of acute HBV infection in humans. We employed fitting of mathematical models to human acute infection data together with physiological constraints to investigate the role of these hypothesized mechanisms in the clearance of infection. Results suggest that cellular proliferation of infected cells resulting in two uninfected cells is required to minimize the destruction of the liver during the clearance of acute HBV infection. In contrast, we find that a cytokine-mediated cure of infected cells alone is insufficient to clear acute HBV infection. In conclusion, our modeling indicates that HBV clearance without lethal loss of liver mass is associated with the production of two uninfected cells upon proliferation of an infected cell.

A Robust and Efficient Numerical Method for RNA-Mediated Viral Dynamics

The multiscale model of hepatitis C virus (HCV) dynamics, which includes intracellular viral RNA (vRNA) replication, has been formulated in recent years in order to provide a new conceptual framework for understanding the mechanism of action of a variety of agents for the treatment of HCV. We present a robust and efficient numerical method that belongs to the family of adaptive stepsize methods and is implicit, a Rosenbrock type method that is highly suited to solve this problem. We provide a Graphical User Interface that applies this method and is useful for simulating viral dynamics during treatment with anti-HCV agents that act against HCV on the molecular level.

Understanding the Complex Patterns Observed during Hepatitis B Virus Therapy

Data from human clinical trials have shown that the hepatitis B virus (HBV) follows complex profiles, such as bi-phasic, tri-phasic, stepwise decay and rebound. We utilized a deterministic model of HBV kinetics following antiviral therapy to uncover the mechanistic interactions behind HBV dynamics. Analytical investigation of the model was used to separate the parameter space describing virus decay and rebound. Monte Carlo sampling of the parameter space was used to determine the virological, pharmacological and immunological factors that separate the bi-phasic and tri-phasic virus profiles. We found that the level of liver infection at the start of therapy best separates the decay patterns. Moreover, drug efficacy, ratio between division of uninfected and infected cells, and the strength of cytotoxic immune response are important in assessing the amount of liver damage experienced over time and in quantifying the duration of therapy leading to virus resolution in each of the observed profiles.

In-host modeling

Understanding the mechanisms governing host-pathogen kinetics is important and can guide human interventions. In-host mathematical models, together with biological data, have been used in this endeavor. In this review, we present basic models used to describe acute and chronic pathogenic infections. We highlight the power of model predictions, the role of drug therapy, and advantage of considering the dynamics of immune responses. We also present the limitations of these models due in part to the trade-off between the complexity of the model and their predictive power, and the challenges a modeler faces in determining the appropriate formulation for a given problem.

Modelling the Impact of Cell-To-Cell Transmission in Hepatitis B Virus

Cell-free virus is a well-recognized and efficient mechanism for the spread of hepatitis B virus (HBV) infection in the liver. Cell-to-cell transmission (CCT) can be a more efficient means of virus propagation. Despite experimental evidence implying CCT occurs in HBV, its relative impact is uncertain. We develop a 3-D agent-based model where each hepatocyte changes its viral state according to a dynamical process driven by cell-free virus infection, CCT and intracellular replication. We determine the relative importance of CCT in the development and resolution of acute HBV infection in the presence of cytolytic (CTL) and non-CTL mechanisms. T cell clearance number is defined as the minimum number of infected cells needed to be killed by each T cell at peak infection that results in infection clearance within 12 weeks with hepatocyte turnover (HT, number of equivalent livers) ≤3. We find that CCT has very little impact on the establishment of infection as the mean cccDNA copies/cell remains between 15 to 20 at the peak of the infection regardless of CCT strength. In contrast, CCT inhibit immune-mediated clearance of acute HBV infection as higher CCT strength requires higher T cell clearance number and increases the probability of T cell exhaustion. An effective non-CTL inhibition can counter these negative effects of higher strengths of CCT by supporting rapid, efficient viral clearance and with little liver destruction. This is evident as the T cell clearance number drops by approximately 50% when non-CTL inhibition is increased from 10% to 80%. Higher CCT strength also increases the probability of the incidence of fulminant hepatitis with this phenomenon being unlikely to arise for no CCT. In conclusion, we report the possibility of CCT impacting HBV clearance and its contribution to fulminant hepatitis.

A Computational Model of Inhibition of HIV-1 by Interferon-Alpha

Type 1 interferons such as interferon-alpha (IFNα) inhibit replication of Human immunodeficiency virus (HIV-1) by upregulating the expression of genes that interfere with specific steps in the viral life cycle. This pathway thus represents a potential target for immune-based therapies that can alter the dynamics of host-virus interactions to benefit the host. To obtain a deeper mechanistic understanding of how IFNα impacts spreading HIV-1 infection, we modeled the interaction of HIV-1 with CD4 T cells and IFNα as a dynamical system. This model was then tested using experimental data from a cell culture model of spreading HIV-1 infection. We found that a model in which IFNα induces reversible cellular states that block both early and late stages of HIV-1 infection, combined with a saturating rate of conversion to these states, was able to successfully fit the experimental dataset. Sensitivity analysis showed that the potency of inhibition by IFNα was particularly dependent on specific network parameters and rate constants. This model will be useful for designing new therapies targeting the IFNα network in HIV-1-infected individuals, as well as potentially serving as a template for understanding the interaction of IFNα with other viruses.

Dynamics of in vivo hepatitis D virus infection

Hepatitis-D virus (HDV) is a satellite virus of hepatitis-B virus (HBV) whose intracellular products are required for the completion of the HDV life cycle. HDV can replicate in a cell without the presence of HBV but needs hepatitis B surface antigen (HBsAg) to complete virus assembly and packaging. In order to better understand HDV dynamics, we developed a mathematical model and successfully simulated HBV and HDV data under a range of scenarios. Compared to HBV mono-infection, dual HDV infection resulted in lower chronic HBV DNA levels, with more marked suppression for coinfection (1 logs HBV DNA copies/ml lower) compared to superinfection (0.6 logs HBV DNA copies/ml). Although they have no effect on HBV, prenylation inhibitors may provide the best therapy for reducing HDV viremia irrespective of the stage in which they are commenced. We found that highly effective long term pegylated interferon (IFN) therapy (99.99%) eliminates HBV and HDV viremia while less effective long term IFN therapy (99%) will only produce approximately 2.03 logs and no decrease in HBV and HDV viremia respectively in both coinfection and superinfection settings. Our study also suggests that there is a substantial difference in the outcome of therapies depending upon the time of commencement.

CONCLUSION

Mathematical modeling of HDV infection can describe the complex interplay between this virus and HBV. Simulations suggest that HDV impacts on the feedback mechanisms that maintain cccDNA levels and that targeting these mechanisms may result in new therapeutic agents for both viruses.

Multiscale modeling of virus replication and spread

Replication and spread of human viruses is based on the simultaneous exploitation of many different host functions, bridging multiple scales in space and time. Mathematical modeling is essential to obtain a systems-level understanding of how human viruses manage to proceed through their life cycles. Here, we review corresponding advances for viral systems of large medical relevance, such as human immunodeficiency virus-1 (HIV-1) and hepatitis C virus (HCV). We will outline how the combination of mathematical models and experimental data has advanced our quantitative knowledge about various processes of these pathogens, and how novel quantitative approaches promise to fill remaining gaps.

A THEORETICAL ASSESSMENT OF THE EFFECTS OF DISTRIBUTED DELAY ON THE TRANSMISSION DYNAMICS OF HEPATITIS B

In this paper, we investigate the global dynamics of a system of delay differential equations which describes the interaction of hepatitis B virus (HBV) with both liver and blood cells. The model has two distributed time delays describing the time needed for infection of cell and virus replication. We also include the efficiency of drug therapy in inhibiting viral production and the efficiency of drug therapy in blocking new infection. We compute the basic reproduction number and find that increasing delays will decrease the value of the basic reproduction number. We study the sensitivity analysis on the key parameters that drive the disease dynamics in order to determine their relative importance to disease transmission and prevalence. Our analysis reveals that the model exhibits the phenomenon of backward bifurcation (where a stable disease-free equilibrium (DFE) co-exists with a stable endemic equilibrium when the basic reproduction number is less than unity). Numerical simulations are presented to evaluate the impact of time-delays on the prevalence of the disease.

Mathematical models of the interrelated dynamics of hepatitis D and B

The hepatitis delta virus (HDV) is a rarest form of viral hepatitis, but has the worst outcomes for patients.It is a subviral satellite dependent on coinfection with hepatitis B (HBV) to replicate within the host liver.To date, there has been little to no modeling effort for HDV. Deriving and analyzing such a mathematical model poses difficulty as it requires the inclusion of (HBV). Here we begin with a well-studied HBV model from the literature and expand it to incorporate HDV. We investigate two models, one with and one without infected hepatocyte replication. Additionally, we consider treatment by the drug lamivudine. Comparison of model simulations with experimental results of lamivudine treatment indicate that infected cell proliferation may play a significant role in chronic HDV infection. Our results also shed light on several questions surrounding HDV and illustrate the need for more data.

In silico single cell dynamics of hepatitis B virus infection and clearance

The progression of acute hepatitis B virus (HBV) to chronic infection or clearance is highly dependent on the host immune response composed of cytolytic (CTL) and non-cytolytic (non-CTL) effects. Cytolytic processes induce hepatocyte killing while non-CTL processes inhibit intracellular replication. Both effects are widely recognized and accepted. However, there are uncertainties about the assistance provided by either the loss of covalently circular closed DNA (cccDNA) during cell proliferation or the emergence of refractory cells to immune mediated clearance. We developed an agent-based mathematical model and tested the relative roles of different mechanisms of the immune system in the clearance of acute HBV infection. HBV viremia clearance time and hepatocyte turnover (HT) were used as the two major criteria in determining reasonable outcomes. Modelling results in 90% of cells containing between 1 and 17 cccDNA copies and normally distributed at the peak of infection. Variations in p36 levels, responsible for determining export of virions or recirculation to amplify cccDNA numbers, have a much greater impact on mean cccDNA level/cell at peak viremia than virus infectivity and cccDNA half-life. A strong CTL effect alone failed to clear infection with HT ≈ 10. Acute infection clearance was possible with combined CTL and non-CTL effects along with complete loss of intracellular viral components during cell proliferation resulting in the desired range of HT (0.7-1). The emergence of cells refractory to infection can reduce HT by up to 90%. However their impact was less effective than complete loss of intracellular viral components during cell proliferation.

CONCLUSION

the existence of refractory cells is not necessary when there is complete loss of intracellular quantities during cell proliferation but is essential with only partial clearance.

Understanding early serum hepatitis D virus and hepatitis B surface antigen kinetics during pegylated interferon-alpha therapy via mathematical modeling

There is little information on the early kinetics of hepatitis delta virus (HDV) and hepatitis B surface antigen (HBsAg) during interferon-α therapy. Here a mathematical model was developed and fitted to frequent HDV and HBsAg kinetic data from 10 patients during the first 28 weeks of pegylated-interferon-α2a (peg-IFN) therapy. Three patients achieved a complete virological response (CVR), defined as undetectable HDV 6 months after treatment stopped with loss of HBsAg and anti-HBsAg seroconversion. After initiation of therapy, a median delay of 9 days (interquartile range [IQR]: 5-15) was observed with no significant changes in HDV level. Thereafter, HDV declined in a biphasic manner, where a rapid first phase lasting for 25 days (IQR: 23-58) was followed by a slower or plateau second phase. The model predicts that the main effect of peg-IFN is to reduce HDV production/release with a median effectiveness of 96% (IQR: 93-99.8). Median serum HDV half-life (t1/2 ) was estimated as 2.9 days (IQR: 1.5-5.3) corresponding to a pretreatment production and clearance of about 10(10) (IQR: 10(9.7) -10(10.7) ) virions/day. None of the patients with flat second phase in HDV achieved CVR. HBsAg kinetics of decline paralleled the second phase of HDV decline consistent with HBsAg-productive-infected cells being the main source of production of HDV, with a median t1/2 of 135 days (IQR: 20-460). The interferon lambda-3 polymorphism (rs12979860) was not associated with kinetic parameters.

CONCLUSION

Modeling results provide insights into HDV-host dynamics, the relationship between serum HBsAg levels and HBsAg-infected cells, IFN's mode of action, and its effectiveness. The observation that a flat second phase in HDV and HBsAg kinetics was associated with failure to achieve CVR provides the basis to develop early stopping rules during peg-IFN treatment in HDV-infected patients.

Antibody responses during hepatitis B viral infection

Hepatitis B is a DNA virus that infects liver cells and can cause both acute and chronic disease. It is believed that both viral and host factors are responsible for determining whether the infection is cleared or becomes chronic. Here we investigate the mechanism of protection by developing a mathematical model of the antibody response following hepatitis B virus (HBV) infection. We fitted the model to data from seven infected adults identified during acute infection and determined the ability of the virus to escape neutralization through overproduction of non-infectious subviral particles, which have HBs proteins on their surface, but do not contain nucleocapsid protein and viral nucleic acids. We showed that viral clearance can be achieved for high anti-HBV antibody levels, as in vaccinated individuals, when: (1) the rate of synthesis of hepatitis B subviral particles is slow; (2) the rate of synthesis of hepatitis B subviral particles is high but either anti-HBV antibody production is fast, the antibody affinity is high, or the levels of pre-existent HBV-specific antibody at the time of infection are high, as could be attained by vaccination. We further showed that viral clearance can be achieved for low equilibrium anti-HBV antibody levels, as in unvaccinated individuals, when a strong cellular immune response controls early infection.

Viral kinetic modeling: state of the art

Viral kinetic (VK) modeling has led to increased understanding of the within host dynamics of viral infections and the effects of therapy. Here we review recent developments in the modeling of viral infection kinetics with emphasis on two infectious diseases: hepatitis C and influenza. We review how VK modeling has evolved from simple models of viral infections treated with a drug or drug cocktail with an assumed constant effectiveness to models that incorporate drug pharmacokinetics and pharmacodynamics, as well as phenomenological models that simply assume drugs have time varying-effectiveness. We also discuss multiscale models that include intracellular events in viral replication, models of drug-resistance, models that include innate and adaptive immune responses and models that incorporate cell-to-cell spread of infection. Overall, VK modeling has provided new insights into the understanding of the disease progression and the modes of action of several drugs. We expect that VK modeling will be increasingly used in the coming years to optimize drug regimens in order to improve therapeutic outcomes and treatment tolerability for infectious diseases.

Human systems immunology: hypothesis-based modeling and unbiased data-driven approaches

Systems immunology is an emerging paradigm that aims at a more systematic and quantitative understanding of the immune system. Two major approaches have been utilized to date in this field: unbiased data-driven modeling to comprehensively identify molecular and cellular components of a system and their interactions; and hypothesis-based quantitative modeling to understand the operating principles of a system by extracting a minimal set of variables and rules underlying them. In this review, we describe applications of the two approaches to the study of viral infections and autoimmune diseases in humans, and discuss possible ways by which these two approaches can synergize when applied to human immunology.

Quantification of viral infection dynamics in animal experiments

Analyzing the time-course of several viral infections using mathematical models based on experimental data can provide important quantitative insights regarding infection dynamics. Over the past decade, the importance and significance of mathematical modeling has been gaining recognition among virologists. In the near future, many animal models of human-specific infections and experimental data from high-throughput techniques will become available. This will provide us with the opportunity to develop new quantitative approaches, combining experimental and mathematical analyses. In this paper, we review the various quantitative analyses of viral infections and discuss their possible applications.

Optimal Antiviral Treatment Strategies of HBV Infection Model with Logistic Hepatocyte Growth

This study considers an optimal therapy strategy for HBV infection by incorporating two controls laws into a previous hepatitis B viral infection model with logistic hepatocyte growth. Our goal is to maximize the number of healthy cells and to minimize the cost of the therapy. In this context, the existence of an optimal control is proved. The optimal control is obtained by solving the optimality system which was composed of three nonlinear ODEs with initial conditions and three nonlinear adjoint ODEs with transversality conditions. The results were analysed and interpreted numerically using MATLAB.

Chronic hepatitis B: what should be the goal for new therapies?

Chronic hepatitis B can currently be medically managed with either pegylated interferon-alpha (pegIFN-α) or one of the five nucleos(t)ide analog Direct Acting Antivirals (DAAs) that inhibit the hepatitis B virus (HBV) DNA polymerase. While pegIFN-α is effective in approximately one-third of the treated patients, the polymerase inhibitors significantly reduce viral load in the vast majority of those treated. However, neither pegIFN-α nor nucleosi(t)de analogs are capable of reliably eliminating the virus and achieving a cure. Moreover, the interferons and polymerase inhibitors are recommended by US, European and Asian professional society practice guidelines for use in only a subset of those infected with HBV. This subset is the population with the greatest levels of circulating viral DNA and abnormal liver function. Although this is the population at the highest risk for cirrhosis and liver cancer, those who fall outside the treatment guidelines, with low levels of viral replication and normal serum ALTs, may also benefit from antiviral therapy. The questions are thus: are new classes of drugs needed to manage chronic hepatitis B? Is a cure possible? Is a cure even necessary? It is therefore important to define the meaning of a cure and determine what the goals of new therapies should be. In this article, we address those questions and propose two operational definitions of medically attainable cures. The first is a "functional cure" based on the clinical outcome, in which the patient's life expectancy becomes the same as that of an individual who has resolved his HBV infection without therapy. Because such an outcome cannot be measured over the short term, we also define an "apparent virological cure," based on the stable off-drug suppression of HBV viremia and antigenemia and the normalization of ALTs and other laboratory tests. We suggest that such a virological cure should be the goal of future therapeutics in all patients with chronic hepatitis B. The extent to which a virological cure predicts a functional cure will only be determined by long-term follow-up.

Modeling the acute and chronic phases of Theiler murine encephalomyelitis virus infection

Theiler murine encephalomyelitis virus (TMEV) infection of a mouse's central nervous system is biphasic: first the virus infects motor neurons (acute phase), and this is followed by a chronic phase in which the virus infects glial cells (primarily microglia and macrophages [M]) of the spinal cord white matter, leading to inflammation and demyelination. As such, TMEV-induced demyelinating disease in mice provides a highly relevant experimental animal model for multiple sclerosis. Mathematical models have proven valuable in understanding the in vivo dynamics of persistent virus infections, such as HIV-1, hepatitis B virus, and hepatitis C virus infections. However, viral dynamic modeling has not been used for understanding TMEV infection. We constructed the first mathematical model of TMEV-host kinetics during acute and early chronic infections in mice and fit measured viral kinetic data with the model. The data fitting allowed us to estimate several unknown parameters, including the following: the rate of infection of neurons, 0.5 × 10(-8) to 5.6 × 10(-8) day(-1); the percent reduction of the infection rate due to the presence of virus-specific antibodies, which reaches 98.5 to 99.9% after day 15 postinfection (p.i.); the half-life of infected neurons, 0.1 to 1.2 days; and a cytokine-enhanced macrophage source rate of 25 to 350 M/day into the spinal cord starting at 10.9 to 12.9 days p.i. The model presented here is a first step toward building a comprehensive model for TMEV-induced demyelinating disease. Moreover, the model can serve as an important tool in understanding TMEV infectious mechanisms and may prove useful in evaluating antivirals and/or therapeutic modalities to prevent or inhibit demyelination.