When Modern Technology Meets Ancient Traditional Chinese Medicine
View this Special IssueResearch Article  Open Access
Modeling and Simulating Dynamics of Complete and PoorResponse Chronic Hepatitis B Chinese Patients for Adefovir and Traditional Chinese Medicine Plus Adefovir Therapy
Abstract
ChiCTRTRC11001263 study was the first largescale doubleblind randomized placebocontrolled traditional Chinese medicines (TCMs) and adefovir (ADV) antihepatitis B virus (HBV) infection trial in the world. A total of 560 hepatitis B e antigen (HBeAg) positive Chinese patients with chronical HBV were randomly classified, in 1 : 1 ratio, into two groups: experimental group (EXG) receiving TCMs + ADV and controlled group (CTG) receiving ADV + TCMplacebo treatment for 48 weeks. This paper introduces two models to model and simulate the evolutions of dynamics for the completeresponse patients and the poorresponse patients in EXG and CTG, respectively. The stimulated mean HBV DNA and alanine aminotransferase (ALT) levels were close to the patients’ experimental data. Analysis and simulations suggest that the activated patients’ immune functions by TCMs + ADV may not only clear infected hepatocytes, but also clear HBV, which made the completeresponse patients’ mean serum HBV DNA levels in EXG reduce rapidly 12 weeks’ earlier than the ones in CTG. One can assume that both the TCMs and ADV have the function of preventing completeresponse patients’ infected hepatocytes from being injured by cytotoxic T lymphocytes (CTLs); the patients’ activated immune cells may also block HBV replications.
1. Author Summary
Nucleoside analogues (NAs), such as lamivudine, adefovir, entecavir, and telbivudine, suppress HBV replication and result in the improvement of the liver architecture. Some TCMs are able to activate patients’ immune function because patients’ serum HBeAg levels may reduce rapidly much earlier before their serum HBV DNA levels decrease significantly. ChiCTRTRC11001263 was the first international registered ADV + TCMplacebo (control group CTG) and TCM + ADVplacebo switching to TCM + ADV (experimental group EXG) antiHBV infection therapy trial. Based on Nowak et al.’s uninfected cellinfected cellfree virus basic virus infection model, this paper introduces two models with additional immune variable and alanine aminotransferase loads to describe and understand the two group patients’ dynamics for antiHBV infection therapy. The results include the determinations of the model parameters, predicting the outcome of the longterm treatment, finding that both the TCMs and ADV may have the function of preventing completeresponse patients’ infected hepatocytes from being injured by CTLs; activated CTLs may also play the role of blocking HBV replications; HBeAg seroconversion may be defined as a predictor that patients can keep their activated immune function via oneyear additional treatment, then ending their therapy.
2. Introduction
Hepatitis B is a lifethreatening liver infection caused by hepatitis B virus (HBV), which can cause chronic liver disease and make people die of cirrhosis of the liver and liver cancer. Two billion people worldwide have been infected with HBV and more than 400 million have chronic (longterm) liver infections. An estimated 1 million people die every year due to the consequences of hepatitis B [1].
The goal of antiCHB infection treatment is to achieve sustained suppression of HBV DNA and remission of liver disease [2]. Nucleoside analogues, such as lamivudine, adefovir, entecavir, and telbivudine, are popular drugs to treat HBV infection. The main role of nucleoside analogues is to block the replication of HBV DNA in vivo.
Some TCMs antiHBV infection therapies have the advantages of rare viral mutation, rare side, and cheap price. Lines of evidence show that TCMs can regulate CHB patients’ immune functions [3].
Monotherapy may have low response rates. Most CHB patients need longterm medication, which can maintain a low response rate after withdrawing drugs and result in higher rate of drug resistance [4–8]. The disadvantages of some NA monotherapies limit the clinical application of CHB patients’ treatments.
NA + TCM therapy has better efficacy than monotone treatment, which is able to increase proportion of patients’ achieving HBeAg loss and clear HBV directly without damaging patients’ hepatocytes [3, 9].
Modelling the dynamics of HBV infection and other virus infections has attracted considerable attentions. Mathematical models play a significant role in improving understanding of the dynamics of the HBV infections in vivo. The models typically used to study HBV dynamics in vivo tend to focus on healthy cells, free virus, and infected cells [10, 11].
The basic viral infection dynamic mathematical model (BVIM) proposed by Nowak et al. [12, 13] has been widely used in the study of the dynamic of infectious agents such as hepatitis B, C and HIV. The BVIM has the following form [12]: where , , and are the numbers of uninfected cells, virusinfected cells, and free virus, respectively. Uninfected cells are produced at a constant rate , die at rate , and become infected at rate . Virusinfected cells are produced at rate and die at rate . Free virus is produced from virusinfected cells at rate and is removed at rate .
Equation (1) has an infection freesteady state : representing an infected person’s complete recovery. Equation (1) has also an endemic steady state : representing an infected person’s persistent infection. Here, is called the basic virus reproductive number of model (1).
It has been proved that if , then the infectionfree steady state of the model (1) is globally attractive; otherwise the endemic steady state of the model (1) is globally attractive [14].
Since in represents the total number of uninfected cells of the patient’s target organ, this implies that an individual with a larger liver will be more difficult to be cured than a person with a smaller one. The meaning of is questionable. Recently, some amended basic viral infection models (ABVIM) [10, 15, 16] are established. One of them takes the following form [15]: where the meanings of the variables , and and the parameters , and are the same as those given in model (1). is the viral infected rate of uninfected cells by free virus and produced rate of virus from virusinfected cells.
The ABVIM has a basic virus reproductive number , which is independent on the total number of cells of the patient’s target organ. It has been proved that if , then the infection free steady state is globally attractive; otherwise the endemic steady state is globally attractive [15, 17].
During the process of viral infections, the immune response has been shown to be universal and necessary to eliminate or control the disease [18, 19]. Actually, in most virus infections, cytotoxic T lymphocytes (CTLs) play a critical role in antiviral defense by attacking virusinfected cells [20].
Therefore, many viral infection dynamic mathematical models with immune response have been studied in recent years [3, 13, 20–22]. One of them has the following form [13]: where the meanings of the variables , and and the parameters , and are the same as those given in model (1). The variable represents the number of cytotoxic T lymphocytes (CTLs). CTLs are produced at rate and die at rate . The term is the death rate of virusinfected cells caused by immune response. Model (5) has a basic virus reproductive number , which is also dependent on the total number of cells of the patient’s target organ. The infection freesteady state of model (5) is independent on the parameters and which relate to the production of CTLs.
Based on the experimental data and previous researches on the dynamics of virus infection model [3, 10, 12, 13, 15, 16, 20–22], this paper introduces two mathematical models to model, simulate, and analyze the dynamics of the evolutions of patients’ mean serum HBV DNA and ALT levels and make longterm prediction for the completeresponse patients and the poorresponse patients for ADV monotherapy and TCM + ADV combination therapy.
3. Methods
3.1. Experiment
ChiCTRTRC11001263 study was a doubleblind randomized placebocontrolled trial. ADV and two kinds of TCMs, Tiaoganjianpihuoxue grain (TCM1) and Tiaoganjieduhuashi grain (TCM2), were used in the trial. TCM1 consists of 13 herbal ingredients, and TCM2 consists of 15 herbal ingredients. A total of 560 Chinese HBeAgpositive CHB patients were randomly classified into, in 1 : 1 ratio, two groups: control group (CTG) and experimental group (EXG).
The patients’ plasma HBV DNA level baselines were 3 by PCR assay. And the alanine aminotransferase (ALT) level baselines were 2 ULN~12 ULN (Disease: ULN, upper limit of normal), where the abbreviation ULN represents “upper limit of normal”. Total bilirubin (TBIL) load baselines were less or equal to 3ULN.
The patients in CTG received ADV (10 mg, once daily) + TCMplacebo (twice daily) for 48 weeks. The patients in EXG were divided into 3 subgroups as follows.(a)Group EXG1 has 207 patients whose ALT levels were larger than 2 ULN and less than 6 ULN.(b)Group EXG2 has 39 patients whose ALT levels were larger than 6 ULN and less than 12 ULN.(c)Group EXG3 has 34 patients whose HBeAg levels were less than 60 S/CO.
The patients’ numbers of the corresponding three subgroups in CTG are 206, 35, and 39, respectively. The control group and the experimental three subgroups have the same characteristics (ITT).
The experimental schemes of the three subEXGs were designed as follows.(a)The patients in EXG1 received TCM1 (twice daily) + ADVplacebo (10 mg once daily) for the first 24 weeks and then switched to TCM2 + AD for additional 24 weeks.(b)The patients in EXG2 received TCM2 (twice daily) + ADVplacebo (10 mg once daily) for the first 24 weeks and then switched to receive TCM2 + AD for additional 24 weeks.(c)The patients in EXG3 continuously received TCM2 (twice daily) + ADV (10 mg once daily) for 48 weeks.(d)During the first 24week therapy if a patient’s ALT level in EXG2 was larger than , then the patient switched to receive TCM2 + ADVplacebo until the 24th week and then switched to receive TCM2 + ADV for additional 24 weeks.(e)For any one in EXG1 or EXG2, if a patient’s HBeAg level was less than 60 S/CO, or ALT level was larger than , or TBIL level was larger than , then the patient switched to receive the scheme of EXG3 for therapy.
The main function of TCM1 is to regulate patients’ immune abilities, and the main role of TCM2 is to block the repletion of HBV. However, efficacy of TCM2 is limited, and it needs NA (e.g., ADV) for combination therapy to increase its efficacy.
The therapy scheme suggests that the patients with lower immune abilities whose ALT levels were less than 6 ULN or HBeAg levels were larger than 60 S/Co or TBIL levels were less than 3 ULN should receive only TCM1 therapy to regulate their immune functions for the first 24 weeks if their tested items did not change to the levels given in item (e).
The conditions item (e) may be a criterion which makes corresponding patients switch to use TCM2 + ADC scheme for further therapy.
Consequently, the purpose of the above therapy scheme was to expect that the 24week therapy would make some patients in groups EXG1 and EXG2 achieve the conditions in item (e). And then the patients in the three groups received the TCM2 + ADV combination treatment for additional 24 weeks.
Some virologic and biochemical responses of the six subgroups are listed in Table 1 (also see [9]). The results show that TCM + ADV antiHBV combination therapy resulted in increased proportion of patients achieving HBeAg loss in the EXG1 and EXG3 versus the CTG1 and CTG3 at week 48. The other virologic and biochemical responses of the controlled group and the experimental group had not significant differences.

At week 48, there were 28 and 31 patients in EXG and CTG who achieved complete response (denoted by CEXG and CCTG), respectively. Meanwhile, there were 42 and 55 patients in EXG and CTG responded poorly (denoted by PEXG and PCTG). Here complete response is defined as HBV DNA level being lower than undetectable level (1000 copies/mL) and HBeAg seroconversion (HBeAg 1 and antiHBe 1); poor response is defined as less than 1 log_{10} copies/mL decrease in HBV DNA level from the baseline at the 48th week.
The outcomes of the patients’ ALT loads in CEXG and CCTG are shown in Figures 1, 2, 3, and 4, respectively. The patients’ mean HBV DNA levels, ALT loads, and HBeAg levels are listed in Table 2 (also see [23]).

Figures 1–4 and the data in Table 2 suggest that the main function of the TCMs is to regulate the patients’ immune functions. The additional 24week TCM + ADV therapy speeded up the patients’ enhancement of immune functions. This observation motivates us to introduce two models to describe the dynamics of antiHBV infection with ADV and TCMs + ADV in the next section.
3.2. Models
Based on the previous work [3, 10, 12, 13, 15, 16, 20–22] and the above analysis, we use model (6) to describe the dynamics of the CTG with the ADV antiHBV infection therapy (similar to that proposed in [22] which does not include the ATL level variable ):
By the similar reasons, we introduce model (7) to describe the dynamics of the EXG with the TCM + ADV antiHBV infection therapy (similar to [3, 24]):
Here, the meanings of the variables , and and the parameters , and are the same as those given in model (1); represents the number of CTLs which are produced at rate and die at rate ; is the death rate of virusinfected cells generated by immune killing; is the clearing rate of virus generated by some specific immune abilities activated via antivirus infection therapy. The variable represents the serum ALT levels. A liver without immune attacking produces ALT at rate , and ALT dies at rate . A CHB patient’s liver produces ALT at rate . () are the treatment efficacy variables during the antiHBV treatment.
Model (6) and model (7) both have the same infectionfree equilibrium :
Model (6) has an endemic equilibrium : representing persistent virus infection, where
Model (7) has an endemic equilibrium :
where
For model (6), the basic virus reproductive number is
For model (7), the basic virus reproductive number is
Similar to [22], we can prove the following theorems.
Theorem 1. Let be defined by (13). If , then the infectionfree equilibrium of (6) is locally stable.
Theorem 2. If then the infectionfree equilibrium of (6) is globally asymptotically stable.
Similar to [3, 24], we can prove the following theorems.
Theorem 3. Let be defined by (14). If , then the infectionfree equilibrium of (7) is locally stable.
Theorem 4. If then the infection free equilibrium of (7) is globally asymptotically stable.
4. Numerical Simulation
In this subsection, in order to interpret more clearly the specific role of TCM in the antiHBV infection therapy, we use model (7) and model (6) to simulate the dynamics of the evolutions of mean serum HBV DNA and ALT levels and make longterm prediction for the CRP and PRP in CTG and EXG antiHBV infection therapies, respectively.
4.1. Simulations for CompleteResponse Patients’ Dynamics in EXG
Model (7) is used to simulate completeresponse patients’ dynamics for TCM + ADV antiHBV infection therapy. Use the methods in [3, 15, 17, 22] to determine approximately the parameters in model (7) as follows.
Because a human liver contains about hepatocytes [13], we obtain
Since the halflife of a hepatocyte is about half a year [25], we get
Assuming the natural death rate of infected hepatocytes is the same as that of uninfected hepatocytes, hence we obtain
A CHB patient typically has between and infected hepatocytes [13]. Different CHB patient’s serum HBV DNA load varies ranging from cps/mL to cps/mL. Hence, we assume that and have the following relation:
We can calculate and via the following equations:
Consequently, we obtain
In the completeresponse patients in EXG for TCM + ADV antiHBV infection therapy, cps/mL. Hence, we calculate
Chronic HBV infection makes some infected hepatocytes undergo apoptosis and be replaced by hepatic stellate cells [26]. Define a parameter , and the patient’s hepatocytes are reduced by percent. Hence, we get where we choose . Furthermore, we obtain
Assuming the halflife of a virus is about one day [13], we obtain
Assuming that the baseline , we get
Since the halflife of CTLs is about 77 days [27], we obtain
Because the halflife of ALT is about 2~3 days [28], we select
We assume that 22 U/L is the mean normal ALT level because the complete response patients’ mean ALT level was about 24 U/L at the week 48. When a human is healthy, in model (7). Hence, we can obtain that because an individual has about about 3liter serum.
Solving the equilibrium point equation gives
Because a healthy Chinese has about counts/L CD8 + T cells, we assume that
Select the mean serum HBV DNA level copies/mL at week 0 as the initial value, and an individual have 3liter serum; hence we determine
Select the mean serum ALT level value 180.06 U/L (see Table 2) at week 0 as the initial value, and an individual have 3liter serum; hence we obtain
Select since none of the available nucleoside analogues inhibitors have been shown to prevent infection of uninfected hepatocytes [29].
In order to agree with the experimental data, the parameters , and need to be changed during the treatment. Their values are listed in Table 3.

Selecting the following initial condition then one can simulate the dynamics of the completeresponse group in EXG for TCM + ADV antiHBV infection therapy.
4.2. Simulations for CompleteResponse Patients’ Dynamics in CTG
Model (6) is used to simulate the completeresponse patients’ dynamics in CTG for ADV antiHBV infection therapy. The parameters , and in model (6) have the same values as those given in ~ and ~ in the above section. Consider the following.(1)Substituting cps/mL into formula (22) gives . (2)Since the mean serum HBV DNA load value at week 0 is , it follows that via in the above section.(3)Since the mean serum ALT level value at week 0 is 198.02 U/L (see Table 2), it follows that via in the above section.
The parameters , and have changed during the treatment. Their values are listed in Table 4.

Selecting the following initial condition: then one can simulate the completeresponse patients’ dynamics in CTG.
4.3. Simulations for PoorResponse Patients’ Dynamics in EXG
Model (7) is used to simulate the poorresponse patients’ dynamics of in EXG for TCM + ADV antiHBV infection therapy. The parameters in model (7) have the same values as those given in and ~ in the above section. Consider the following.(1)Substituting cps/mL into formula (22) gives . (2)Since the mean serum HBV DNA load value at week 0 is , it follows that . (3)Since the mean serum ALT level value at week 0 is U/L, it follows that .
The parameters , and have changed during the treatment. Their values are listed in Table 5.

Selecting the following initial condition: then one can simulate the poorresponse patients’ dynamics in EXG for TCM + ADV antiHBV infection therapy.
4.4. Simulations for PoorResponse Patients’ Dynamics in CTG
Model (6) is used to simulate the poorresponse patients’ dynamics in CTG for ADV antiHBV infection therapy. The parameters , and in model (6) have the same values as those given in and in the above section. Consider the following.(1)Substituting cps/mL into formula (22) gives . (2)Since the mean serum HBV DNA level at week 0 is , it follows that via in the above section. (3)Since the serum ALT level value at week 0 is 128.36 U/L, it follows that .
The parameters , and have changed during the treatment. Their values are listed in Table 6.

Selecting the following initial condition: then one can simulate the dynamics of the poor response patients in CTG for ADV antiHBV infection therapy.
5. Results
The numerical simulations of the evolution dynamics of patients mean serum HBV DNA and ALT levels for the four subgroups are shown in Figures 5 and 6, respectively. Observe that the stimulated evolutions of the mean serum HBV DNA levels and ALT levels are close to the experimental data.
(a)
(b)
(c)
(a)
(b)
(c)
The numerical simulations of the patients’ dynamics of the antiHBV infection therapies give the following results.
At the week 0 (baseline), the basic virus reproductive numbers and of the 4 subgroups CEXG, CCTG, PEXG, and CTG are 1.0526, 1.0527, 1.0527, and 1.0527, respectively. This can interpret why they become viruspersistent CHB patients.
After the 48week therapy, the basic virus reproductive numbers and of the two subgroups CEXG and CCTG were reduced to and , respectively. Further simulations show that it needs about 6.5 and 6.8 years of treatment to make all infected hepatocytes be replaced by normal ones, respectively.
After the 48week therapy, the basic virus reproductive numbers and of the two subgroups PEXG and PCTG were only reduced to 0.84085 and 0.87121, respectively. Further simulations show that the poorresponse patients in EXG and CTG cannot be recovered completely until 20 years of treatments.
Figure 5 and Table 2 show that the completeresponse patients’ mean serum HBV DNA levels in CEXG have reduced rapidly during week 24 and week 36, which was 12 weeks earlier than the complete response patients in CCTG.
In order to model this phenomenon, the parameter in model (7) related to the production rate of CTLs has increased from to while the parameter in model (7) has been designed to increase form to during week 24 and week 36. represents the clearing rate of virus generated by some specific immune abilities activated via antiHBV infection therapy, which can clear HBV directly.
Combining the trial data listed in Tables 1 and 2 and the modeling data given in Tables 3 and 4 follows that the TCMs + ADV combination therapy may offer superior efficacy for suppressing HBV replications than monotone ADV therapy.
Comparing the model parameters given in Tables 3–6 gives the following.(a)At week 48, the poorresponse patients’ parameter value on is much smaller than the completeresponse patients’ one.(b)During weeks 25 to 48, the poorresponse patients’ parameter value on increased while the completeresponse patients’ one kept unchanged.
The above results imply that for the poorresponse patients, the drug resistance made the therapy efficacy (parameter ) reduce rapidly, and the virus replication rate (parameter ) was increased quickly.
6. Discussion
Based on the experimental data of CHB patients’ serum HBV DNA levels and ALT levels, this paper introduces two differential equation models (6) and (7) to describe the CHB patients’ dynamics for the ADV monotone treatment and the TCMs + ADV combination therapy. An amended term related to the ability for killing infected hepatocytes is included in the models to describe the evolution of the patients’ ALT levels.
Making some simplified assumptions, one can determine 11 of the 13 parameters in (6) and (7). The simulation results are close to the patients’ mean HBV DNA levels and mean ALT levels.
Based on the experimental data (see Table 2) and the simulation results, one can propose the following hypotheses.
Hypothesis (a). Both the TCMs and ADV have the function that prevents completeresponse patients’ infected hepatocytes from being injured by CTLs; that is, the killing parameter in (6) and (7) becomes smaller than its baseline values.
This hypothesis may interpret why the completeresponse patients’ ALT loads decreased quickly while their HBV DNA levels decreased slowly or increased during the first 12week therapy (see Tables 2, 3, and 4 and Figure 5).
Clinically, some patients with NA or TCM treatments may show serum HBV DNA levels to rebound higher than their baseline levels after cessation treatments. Hypothesis (a) may interpret that the patients kept the function of preventing infected hepatocytes from being injured by CTLs after stopping therapy.
The experimental data (see Table 2) and the simulation results (see Figure 6) suggest that the Chinese patients with high baseline HBV DNA levels, as well as HBeAg loads, and low baseline ALT levels may not obtain complete responses for the ADV or the CTMs + ADV treatments in 48 weeks.
The numerical simulations show that for the completeresponse patients in EXG and CTG, it needs about 6.5 and 6.8 years of treatment to make all infected hepatocytes be replaced by normal if no virus mutations will appear and the efficacy of the therapy will be kept.
Clinically, a complete response CHB patient with nucleoside analogues treatment usually needs much longer times to obtain hepatitis B surface antigen loss. This fact suggests that patients’ activated immune abilities may decrease as patients’ HBV DNA levels decrease to very low levels.
At week 48, the very high efficacy (see given in Tables 3 and 4 and Figure 5) of suppressing HBV replications makes us propose the second hypothesis.
Hypothesis (b). The efficacy of blocking HBV replications is not generated via TCMs and/or ADV alone. The CTLs (represented by variable in (6) and (7)) efficiently control HBV replication by noncytolytic mechanisms [30] contributing also to block HBV replications.
This hypothesis may interpret why some patients’ HBV DNA levels reduced rapidly at some specific time during their treatments because the activated noncytolytic mechanisms of CTL may play roles.
Based on a review article on the endpoints of hepatitis B treatment [31] and hypothesis (b), one can propose the following.
Hypothesis (c). For completeresponse CHB patients with nucleoside analogues treatments, additional oneyear consolidation therapy can make most patients keep their activated immune abilities (parameter in (6) and (7)), contributing to the treatment efficacy (parameter ) after finishing the consolidation treatment.
This hypothesis may interpret Recommendation 9 given in the Chronic Hepatitis B Guideline of the AsianPacific Association for the Study of the Liver.
For oral antiviral agents, in HBeAgpositive patients, treatment can be stopped when HBeAg seroconversion with undetectable HBVDNA has been documented on 2 separate occasions at least 6 months apart [32].
Based on hypothesis (c), we assume that after finishing twoyear treatment, the completeresponse patients in EXG and CTG keep their immune parameter unchanged, returns to baseline, , and efficacy parameter reduces to . The simulated evolutions dynamics of HBV DNA levels and ALT levels are shown in Figure 7. Observe that the treatment benefits are kept.
(a)
(b)
Modeling the dynamics for the ADV monotone treatment and the TCMs + ADV combination therapy may also provide some theoretical interpretation for the medical statistic results (see Table 1).
Since the TCMs + ADV therapy made the patients have an additional immune term (see the third equation in (7)), the TCMs + ADV therapy significantly resulted in increased proportion of the patients achieving HBeAg loss in the experimental group (see Table 1).
The modeling analysis with the experimental data analysis motivates to propose the previous three hypotheses, which may interpret some clinical experience judgements. The dynamics of antiHBV infection therapy are very complex. It is difficult to set up mathematical model to describe them accurately. However, modelling dynamics of antiHBV infection therapy would enable a better understanding, prediction, and design of antiHBV infection treatments.
Conflict of Interests
The authors have no conflict of interests to declare.
Acknowledgments
The authors would like to thank the referees for their valuable comments. The authors wish also to thank the following people who performed the experiments: Fengzhen Zhao (The First Hospital affiliated to Tianjin University of Chinese Medicine), Daqiao Zhou (Shenzhen College of Chinese Medicine), Mingxiang Zhang (Shenyan 6th People Hospital), Jingdong Xue (Shanxi Provincial Hospital of Chinese Medicine), Tiejun Liu (Hospital affiliated to Changchun University of Chinese Medicine), Xiaoling Chi (Guangdong Provincial Hospital of Chinese Medicine), Xianbo Wang (Ditan Hospital Affiliated to Capital University of Medicine), Bingjiu Lu (Liaoning Provincial Hospital of Chinese Medicine), Jun Li (302 Military Hospital of China), Qin Li (Fuzhou Infection Disease Hospital), Dewen Mao (The First Hospital affiliated to Guangxi College of Chinese Medicine), Huasheng Yang (Youan Hospital affiliated to Capital University of Medicine), Hongzhi Yang (The Third Hospital affiliated to Zhongshan University), Wenxia Zhao (The First Hospital affiliated to Henan College of Chinese Medicine), Yong Li (The Hospital Affiliated to Shandong University of Chinese Medicine), Guoliang Zhang (Anhui Provincial Hospital of Chinese Medicine), and Feng Jiang (Dongzhimen Hospital affiliated to Beijing University of Chinese Medicine). This work is jointly supported by the 11th 5Year Plan Key Research Project of China (no. 2008ZX10005006), the NNSF of China (Grant no. 61074192).
References
 World health Organization, “Hepatitis B: Fact Sheet,” 2013, http://www.who.int/mediacentre/factsheets/fs204/en/. View at: Google Scholar
 A. S. Lok and B. J. McMahon, “Chronic hepatitis B,” Hepatology, vol. 45, no. 2, pp. 507–539, 2007. View at: Publisher Site  Google Scholar
 X. Chen, L. Min, Y. Ye, and Y. Zheng, “Modeling and simulation for dynamics of antiHBV infection therapy,” Lecture Notes in Electrical Engineering, vol. 123, no. 2, pp. 557–566, 2011. View at: Publisher Site  Google Scholar
 E. B. Keeffe, S. Zeuzem, R. S. Koff et al., “Report of an International workshop: roadmap for management of patients receiving oral therapy for chronic hepatitis B,” Clinical Gastroenterology and Hepatology, vol. 5, no. 8, pp. 890–897, 2007. View at: Publisher Site  Google Scholar
 F. Zoulim and R. Perrillo, “Hepatitis B: reflections on the current approach to antiviral therapy,” Journal of Hepatology, vol. 48, supplement 1, pp. S2–S19, 2008. View at: Publisher Site  Google Scholar
 Y. F. Liaw, “On treatment outcome prediction and adjustment during chronic hepatitis B therapy: now and future,” Antiviral Therapy, vol. 14, no. 1, pp. 13–22, 2009. View at: Google Scholar
 A. S. Lok and B. J. McMahon, “Chronic hepatitis B: update 2009,” Hepatology, vol. 50, no. 3, pp. 661–662, 2009. View at: Publisher Site  Google Scholar
 European Association for The Study of The Liver, “EASL clinical practice guidelines: management of chronic hepatitis B,” Journal of Hepatology, vol. 50, no. 2, pp. 227–242, 2009. View at: Publisher Site  Google Scholar
 Y. Ye, L. Min, Q. Zhang et al., “Evaluation of 48 week adefovior dipvoxl (AD) and Chinese hebral medicine plus AD treatment in HBeAg(+) chronic hepatitis B Chinese patients: a doublebind randomized trial,” in Proceedings of the 62st Annual Meeting of the American Association for the Study of Liver Diseases: The Liver Meeting 2011, vol. 54, p. 1047A, San Francisco, NC, USA, 2011. View at: Google Scholar
 S. Eikenberry, S. Hews, J. D. Nagy, and Y. Kuang, “The dynamics of a delay model of hepatitis B virus infection with logistic hepatocyte growth,” Mathematical Biosciences and Engineering, vol. 6, no. 2, pp. 283–299, 2009. View at: Publisher Site  Google Scholar
 S. Hews, S. Eikenberry, J. D. Nagy, and Y. Kuang, “Rich dynamics of a hepatitis B viral infection model with logistic hepatocyte growth,” Journal of Mathematical Biology, vol. 60, no. 4, pp. 573–590, 2010. View at: Publisher Site  Google Scholar
 M. A. Nowak, S. Bonhoeffer, A. M. Hill et al., “Viral dynamics in hepatitis B virus infection,” Proceedings of the National Academy of Sciences of the United States of America, vol. 93, no. 9, pp. 4398–4402, 1996. View at: Publisher Site  Google Scholar
 M. A. Nowak and R. M. May, Virus Dynamics: Mathematical Principles of Immunology and Virology, Oxford University, Oxford, UK, 2000.
 P. D. Leenheer and H. L. Smith, “Virus dynamics: a global analysis,” SIAM Journal on Applied Mathematics, vol. 63, no. 4, pp. 1313–1327, 2003. View at: Publisher Site  Google Scholar
 L. Min, Y. Su, and Y. Kuang, “Mathematical analysis of a basic virus infection model with application to HBV infection,” Rocky Mountain Journal of Mathematics, vol. 38, no. 5, pp. 1573–1585, 2008. View at: Publisher Site  Google Scholar
 S. A. Gourley, Y. Kuang, and J. D. Nagy, “Dynamics of a delay differential model of hepatitis B virus infection,” Journal of Biological Dynamics, vol. 2, no. 2, pp. 140–153, 2008. View at: Publisher Site  Google Scholar
 Y. Zheng, L. Min, Y. Ji, Y. Su, and Y. Kuang, “Global stability of endemic equilibrium point of basic virus infection model with application to HBV infection,” Journal of Systems Science and Complexity, vol. 23, no. 6, pp. 1221–1230, 2010. View at: Publisher Site  Google Scholar
 D. Kägi and H. Hengartner, “Different roles for cytotoxic T cells in the control of infections with cytopathic versus noncytopathic viruses,” Current Opinion in Immunology, vol. 8, no. 4, pp. 472–477, 1996. View at: Publisher Site  Google Scholar
 J. E. Schmitz, M. J. Kuroda, S. Santra et al., “Control of viremia in simian immunodeficiency virus infection by CD8+ lymphocytes,” Science, vol. 283, no. 5403, pp. 857–860, 1999. View at: Publisher Site  Google Scholar
 Q. Z. Xie, D. W. Huang, S. D. Zhang, and J. Cao, “Analysis of a viral infection model with delayed immune response,” Applied Mathematical Modelling, vol. 34, no. 9, pp. 2388–2395, 2010. View at: Publisher Site  Google Scholar
 D. Wodarz, J. P. Christensen, and A. R. Thomsen, “The importance of lytic and nonlytic immune responses in viral infections,” Trends in Immunology, vol. 23, no. 4, pp. 194–200, 2002. View at: Publisher Site  Google Scholar
 Y. Zheng, L. Min, X. Chen, and Y. Ye, “Dynamic analysis of HBV infection model with simulations for antiHBV infection therapy,” in 2011 IEEE International conference on Intelligent Computation and BioMedical Instrumentation, ICBMI 2011, pp. 291–295, chn, December 2011. View at: Publisher Site  Google Scholar
 L. Min, X. Chen, Y. Ye et al., “Modeling and simulating dynamics of complete and poor response hbeag positive chronic hepatitis B chinese patients for adefovir (adv) and traditional chinese medicine plus adv therapy,” in Proceedings of the 22st Conference of the Asian Pacific Association for the Study of the Liver, Hepatology Intenational, vol. 6, p. 62, Taipei, Taiwan, February 2012. View at: Google Scholar
 X. Chen, L. Min, Y. Zheng et al., “Modeling and simulation for dynamics of antiHBV infection therapy,” Computer Engineering and Applications, vol. 48, pp. 20–27, 2012 (Chinese). View at: Google Scholar
 A. S. Perelson, “Modelling viral and immune system dynamics,” Nature Reviews Immunology, vol. 2, no. 1, pp. 28–36, 2002. View at: Google Scholar
 R. Bataller and D. A. Brenner, “Liver fibrosis,” Journal of Clinical Investigation, vol. 115, no. 2, pp. 209–218, 2005. View at: Publisher Site  Google Scholar
 M. Hellerstein, M. B. Hanley, D. Cesar et al., “Directly measured kinetics of circulating T lymphocytes in normal and HIV1infected humans,” Nature Medicine, vol. 5, no. 1, pp. 83–89, 1999. View at: Publisher Site  Google Scholar
 S. L. Stockham and M. A. Scott, Fundamentals of Veterinary Clinical Pathology, Iowa State University Press, Ames, Iowa, USA, 2002.
 F. Zoulim, M. Buti, and A. S. Lok, “Antiviralresistant hepatitis B virus: can we prevent this monster from growing?” Journal of Viral Hepatitis, vol. 14, no. 1, pp. 29–36, 2007. View at: Publisher Site  Google Scholar
 J. J. Chang, A. J. V. Thompson, K. Visvanathan et al., “The phenotype of hepatitis B virusspecific T cells differ in the liver and blood in chronic hepatitis B virus infection,” Hepatology, vol. 46, no. 5, pp. 1332–1340, 2007. View at: Publisher Site  Google Scholar
 W. Chotiyaputta and A. S. F. Lok, “Endpoints of hepatitis B treatment,” Journal of Viral Hepatitis, vol. 17, no. 10, pp. 675–684, 2010. View at: Publisher Site  Google Scholar
 Y. Liaw, N. Leung, J. Kao et al., “AsianPacific consensus statement on the management of chronic hepatitis B: a 2008 update,” Hepatology International, vol. 2, no. 3, pp. 263–283, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Lequan Min et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.