Innovative Approximate Analytical Solution for Standard Model of Viral Dynamics: Hepatitis C with Direct-Acting Agents as an Implemented Case
In this article, a novel approximate analytical solution is presented for solving the standard viral dynamic model. Basically, the standard model is used to study viral dynamics in patients for a wide range of viruses like HIV, HPV, HBV, and HCV. In this research work, the standard model for hepatitis C virus (HCV) is considered in detail; however, the analysis and results can be applicable for all other viruses. This standard model is used to study viral dynamics in patients treated with direct-acting antiviral agents (DAAs). Power series solution combined with Laplace–Padé resummation method (PSLP) is used to obtain the approximate analytical solutions for the model. To test the capability as well as the validity of the proposed method, results are compared with available viral load data published in the literature and with published simulation results. Good fits are obtained in the comparison for all cases considered. Given that medical specialists and physicians are more interested in solutions that yield direct and simple predictions, it is expected that the proposed approximate analytical solution would be attractive to them and help them to obtain a straightforward and a proper estimation regarding the viral load due to variations in treatment and/or patient’s parameters.
Treatment for chronic hepatitis C infection began in the early 1990s with interferon-alpha . This injectable drug worked by improving the immune system, rather than by specifically attacking the virus. In 1998, the oral drug ribavirin was added to interferon . The development of the treatment occurred in 2002 with the approval of pegylated interferon-alpha, a process that makes interferon more durable and effective . In 2011, antiviral medications that stop the virus reproducing, which called direct-acting antivirals (DAAs), appeared . Sofosbuvir, commercial name Sovaldi, was approved in the United States in December 2013. Sovaldi-based oral therapy offers high cure rates for HCV infection with excellent tolerability . Currently, the recommended DAA regimens are combinations of multidrug treatments, containing in particular NS5A, NS5B, NS4A inhibitors, such as glecaprevir + pibrentasvir, sofosbuvir and ribavirin (SOF + RBV), telaprevir + pegylated interferon and ribavirin (TVR + PR), sofosbuvir + simeprevir (SOF + SIM), and sofosbuvir + velpatasvir. For more details, see [6, 7, 8]. These regimens offer high cure rates with a rate of curation near 100%.
The quick step of HCV drug development has led to the hopeful estimate that fully eradication of HCV is possible. Although there is still no vaccine for HCV, cure and eradication happened, but further study for the factors that increase this eradication rate is needed. The World Health Organization (WHO) has formulated the Global Health Sector Strategy on Viral Hepatitis, 2016–2021, and established service coverage targets to eliminate HCV as a public health threat by 2030 [9, 10]. Meanwhile, there remain many barriers that need to be overcome. Such barriers include the development of simplified and highly effective drug regimens, improving the rates of detection of infection, and the availability of funds including financial and medical expertise.
Mathematical modeling is a useful tool in the study of virus dynamics for many types of viruses such as HCV, HBV, HPV, and HIV. It can be used to predict behavior under certain conditions or decide which parameters enhance the spread of disease. It may also be used to calculate the medications required to eradicate the disease or, at least, get it under control. These models can be used also to understand the biological mechanisms and interpret the experimental results . Mathematical models in the form of a system of differential equations for the basic dynamics in vivo were developed and analyzed for HCV [8, 12–15], HBV [16, 17], HIV [18–22], and HPV . Models for HCV treatment with DAAs therapy are considered in [8, 14, 15]. Authors in [24–26] developed a multiscale model that considered both intracellular viral RNA replication and extracellular viral infection.
In general, numerical solutions were obtained for the models. However, analytical solutions could be useful for the estimation of parameters and for direct and simple predictions for the viral loads. Seldom, an analytical or an approximate analytical solution can be found in the literature. A simplified analytical solution for the standard model for HCV viral dynamics was constructed in . It was assumed that the number of target cells is constant; hence, the system of ordinary differential equations (ODE) is reduced, linearized, and solved. It can be noticed that this system of ODE contains 7 parameters while the obtained solutions contain only 3 parameters. An analytical solution was obtained for the multiscale model in [24, 26]. All new infections during therapy were ignored in that analysis.
In this research work, the approximate analytical solution is obtained for the standard model. Similar standard models have been used for HCV, HBV, and HIV in the literature [8, 12–23]. Nevertheless, the model considered is the HCV model. There is no reduction for the system of ODE or assumptions to simplify the equations of the model. Yet, the approximation is in the solution not in the system of differential equations. Therefore, it is more accurate and satisfactory than the previously mentioned analytical solutions. Power series solutions are first obtained for the system of ODE, and then, the Laplace–Padé resummation method (PSLP) is used to obtain the entailed approximate analytical solution. The technique used depends on the methodology developed in [27–34] for the general system of ODE. The solution is used for the analysis of the standard viral dynamic model of HCV after any type of DAA treatment initiation.
2. Standard Model
Consider the system of nonlinear ODEs for the standard viral dynamic mathematical model for HCV kinetics during treatment [11–14]:where is the target cells, which are produced at a constant rate and are assumed to die per capita rate , and are infected by virus, , at rate . Infected cells are assumed to die at per capita rate . Virions are generated at rate per infected cell and cleared from serum at rate per virion. Treatment is assumed to reduce the average viral production rate per infected cell from to , where is the in vivo antiviral effectiveness of therapy (). At treatment initiation , we used the standard assumption [10, 21, 22] that the system in the pretreatment is in steady state given by
3. Solution Method
3.1. Power Series Solution
Assuming the following solutions of order N,
Substituting equations (6)–(8) into equations (1)–(3) and equating terms having the same powers of . The , , and coefficients can be calculated. For instance, the power series solution for iswhere is the initial condition of and the coefficients , , and are as follows:
3.2. Laplace–Padé Resummation
The Padé approximant for a series function is a function with the same power series expansions as the original series function. It is suitable for approximating a divergent series function. The approximant is derived by expanding the function as a ratio of two power series and determining both the numerator and denominator coefficients [28–34]. If we have a function that can be represented in a power series form as
A Padé approximant is a rational function:which has a Maclaurin expansion that agrees with the series of equation (11) approximately. is chosen to be equal to . So there are independent numerator coefficients and independent denominator coefficients, making unknown coefficients. The should fit the power series of equation (11) through the orders 1, , , … , . Hence,
Up to the term of order , we can write
Denominator coefficients can be found from equation (14) by equating the coefficients of . Numerator coefficients can be found from equation (14) by equating the coefficients of . Hence, we have constructed the Padé approximant which agrees with through the order .
To extend the convergence zone of the power series solution, Laplace–Padé resummation can be introduced. In Laplace–Padé resummation, the Laplace transform is applied to the given series function, then Padé approximate is obtained for the resulted function, and consequently, the inverse Laplace transform is used to specify the final solution in the form of exponential functions. To apply Laplace–Padé resummation to the series solution of the viral load , the following procedure is employed:(1)From equation (9), the third-order power series solution of is(2)Taking Laplace transform to the series solution in equation (15) to obtain(3) is written instead of in equation (16) to obtain(4) Padé approximant for equation (17) can be written as
Then, equation (19) can be written as
Denominator coefficients and can be found by equating the coefficients of and . Numerator coefficients , and can be found by equating the coefficients of , and . Using Mathematica program, hence(5) is written instead of in ; hence, is obtained as(6)Using inverse Laplace transform for , the viral load is obtained aswhere
Equation (23) combined with equations (10) and (24) represents the general approximate analytical solution for the viral load for the standard dynamic model given in equations (1)–(3). It is worth to notice that this solution contains the 7 parameters given in the system of equations (1)–(3). Hence, it reflects properly the biological features inherited in the system. It can be noticed that the simplified approximate solution given in equation (23) contains only 3 parameters.
4. Study Cases
To illustrate the capabilities and competences of the proposed method, four cases of study are presented. These cases consider the analysis of viral kinetics using the standard viral dynamic mathematical model of HCV. The proposed PSLP method is applied to solve the nonlinear dynamic model of the viral kinetics for some patients after the initiation of treatment with DAAs.
For the first two study cases, the predictions of the PSLP method are compared, for each patient, with viral load data available in the literature [24, 25]. To provide best fits of data for each patient, parameters , and are estimated using the relation for initial condition and equation (23). In the first case, patients were infected with HCV genotype 1 and treated with danoprevir. The viral load for these patients is checked for 13 days after initiation of danoprevir, the data are available in , and the parameters are given in Table 1. In the second case, patients were treated with daclatasvir and the viral load of patients is checked for 2 days after initiation of daclatasvir, the data are available in , and the parameters are given in Table 2. Figures 1 and 2 demonstrate the comparison between the solution of the PSLP method and the corresponding viral load data for each patient.
Two other cases representing recently used combination of multidrug DAA treatments are considered. The predictions of the PSLP method are compared with the corresponding viral load obtained by simulation in . The used parameter values are given in Table 3, where have the same values assigned in  and has been chosen accordingly. The PSLP results and the published simulation results for 25 days after the initiation of treatment with TVR + PR are shown in Figure 3. The PSLP results and the published simulation results for 28 days after the initiation of treatment with SOF + SIM are shown in Figure 4. In both cases, the PSLP solution and the simulation results are almost identical. The comparison shows that the proposed PSLP method provides adequate approximate analytical solutions using the standard viral dynamic model for all the considered cases. It is worth to mention that, in these two cases, the rate of curation is near 100%.
Since it is highly desirable to predict the patient’s response to a specific treatment regimen prior to the starting of the treatment itself, the PSLP solution offers a simple and powerful tool for medical specialists and physicians to perform this task. They can use the patient’s parameters to calculate the constants in equations (10) and (24) and substitute them in equation (23) to get a closed-form solution for the viral load. Hence, the viral load can be plotted versus time or simply the viral load can be estimated at any instant by direct substitution in the closed-form solution.
Power series solution combined with the Laplace–Padé resummation method (PSLP) has been used to obtain a general approximate analytical solution for the nonlinear standard viral dynamic model of HCV for patients treated with DAAs. However, the solution is not limited to HCV model, and it can be applied to other viruses like HIV and HBV, for example, given that proper parameters are used. To test the applicability and accuracy of the proposed method, results have been compared with viral load data and with published simulated results. Satisfactory agreement between the PSLP solution and the corresponding viral load data has been found for all the considered cases. The PSLP results and the published simulated results are almost identical. The comparison proves that this innovative PSLP solution can be used with confidence for solving the nonlinear standard viral dynamic model. This solution can conveniently be used to fit patient data and estimate parameter values. So it would facilitate for physicians to monitor the changes in the viral load due to changes in treatment and to deal with changes in patient’s parameters.
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
T. H. T. Nguyen, J. Guedj, S. L. Uprichard, A. Kohli, S. Kottilil, and A. S. Perelson, “The paradox of highly effective sofosbuvir-based combination therapy despite slow viral decline: can we still rely on viral kinetics?” Scientific Reports, vol. 7, no. 1, Article ID 10233, 2017.View at: Publisher Site | Google Scholar
World Health Organization, Global Hepatitis Report 2017, Global Hepatitis Programme, Department of HIV/AIDS, WHO, Geneva, Switzerland, 2017.
World Health Organization, Combating Hepatitis B and C to Reach Elimination by 2030, World Health Organization, Geneva, Switzerland, 2016.
D. S. Jones and B. D. Sleeman, “Differential equations and mathematical biology,” in Mathematical Biology and Medicine Series, Chapman & Hall/CRC, Boca Raton, FL, USA, 2003.View at: Google Scholar
J. Guedj and A. U. Neumann, “Understanding hepatitis C viral dynamics with direct-acting antiviral agents due to the interplay between intracellular replication and cellular infection dynamics,” Journal of Theoretical Biology, vol. 267, no. 3, pp. 330–340, 2010.View at: Publisher Site | Google Scholar
J. Guedj, H. Dahari, L. Rong et al., “Modeling shows that the NS5A inhibitor daclatasvir has two modes of action and yields a shorter estimate of the hepatitis C virus half-life,” Proceedings of the National Academy of Sciences, vol. 110, no. 10, pp. 3991–3996, 2013.View at: Publisher Site | Google Scholar
H. A. Elkaranshawy, A. M. Abdelrazek, and H. M. Ezzat, “Power series solution to sliding velocity in three-dimensional multibody systems with impact and friction,” International Journal of Mathematical, Computational, Physical, Electrical and Computer Engineering, vol. 9, no. 10, 2015.View at: Google Scholar
I. M. Abdelrazik, “Improving the solutions of parker sochacki method for nonlinear ordinary differential equations,” Department of Engineering Mathematics and Physics, Faculty of Engineering, Alexandria University, Alexandria, Egypt, 2016, PhD thesis.View at: Google Scholar
G. A. Baker and P. R. Graves-Morris, Padé Approximants, vol. 59, Cambridge University Press, Cambridge, UK, 1996.