Optimal control analysis of malaria-schistosomiasis co-infection dynamics

  • This paper presents a mathematical model for malaria-schistosom-iasis co-infection in order to investigate their synergistic relationship in the presence of treatment. We first analyse the single infection steady states, then investigate the existence and stability of equilibria and then calculate the basic reproduction numbers. Both the single-infection models and the co-infection model exhibit backward bifurcations. We carrying out a sensitivity analysis of the co-infection model and show that schistosomiasis infection may not be associated with an increased risk of malaria. Conversely, malaria infection may be associated with an increased risk of schistosomiasis. Furthermore, we found that effective treatment and prevention of schistosomiasis infection would also assist in the effective control and eradication of malaria. Finally, we apply Pontryagin's Maximum Principle to the model in order to determine optimal strategies for control of both diseases.

    Citation: Kazeem Oare Okosun, Robert Smith?. Optimal control analysis of malaria-schistosomiasis co-infection dynamics[J]. Mathematical Biosciences and Engineering, 2017, 14(2): 377-405. doi: 10.3934/mbe.2017024

    Related Papers:

    [1] Ivana Sirangelo, Margherita Borriello, Gaetano Irace, Clara Iannuzzi . Intrinsic blue-green fluorescence in amyloyd fibrils. AIMS Biophysics, 2018, 5(2): 155-165. doi: 10.3934/biophy.2018.2.155
    [2] Jonathan Trauth, Johannes Scheffer, Sophia Hasenjäger, Christof Taxis . Strategies to investigate protein turnover with fluorescent protein reporters in eukaryotic organisms. AIMS Biophysics, 2020, 7(2): 90-118. doi: 10.3934/biophy.2020008
    [3] Marek K. Korzeniowski, Barbara Baird, David Holowka . STIM1 activation is regulated by a 14 amino acid sequence adjacent to the CRAC activation domain. AIMS Biophysics, 2016, 3(1): 99-118. doi: 10.3934/biophy.2016.1.99
    [4] David E. Shoup . Diffusion-controlled reaction rates for clusters of binding sites on a cell. AIMS Biophysics, 2016, 3(4): 522-528. doi: 10.3934/biophy.2016.4.522
    [5] Vladimir B. Teif, Andrey G. Cherstvy . Chromatin and epigenetics: current biophysical views. AIMS Biophysics, 2016, 3(1): 88-98. doi: 10.3934/biophy.2016.1.88
    [6] Megan A. Ahern, Claudine P. Black, Gregory J. Seedorf, Christopher D. Baker, Douglas P. Shepherd . Hyperoxia impairs pro-angiogenic RNA production in preterm endothelial colony-forming cells. AIMS Biophysics, 2017, 4(2): 284-297. doi: 10.3934/biophy.2017.2.284
    [7] Ateeq Al-Zahrani, Natasha Cant, Vassilis Kargas, Tracy Rimington, Luba Aleksandrov, John R. Riordan, Robert C. Ford . Structure of the cystic fibrosis transmembrane conductance regulator in the inward-facing conformation revealed by single particle electron microscopy. AIMS Biophysics, 2015, 2(2): 131-152. doi: 10.3934/biophy.2015.2.131
    [8] Larisa A. Krasnobaeva, Ludmila V. Yakushevich . DNA kinks behavior in the potential pit-trap. AIMS Biophysics, 2022, 9(2): 130-146. doi: 10.3934/biophy.2022012
    [9] Yu-Ting Huang, Hui-Fen Liao, Shun-Li Wang, Shan-Yang Lin . Glycation and secondary conformational changes of human serum albumin: study of the FTIR spectroscopic curve-fitting technique. AIMS Biophysics, 2016, 3(2): 247-260. doi: 10.3934/biophy.2016.2.247
    [10] Larisa I. Fedoreyeva, Tatiana A. Dilovarova, Boris F. Vanyushin, Inna A. Chaban, Neonila V. Kononenko . Regulation of gene expression in Nicotiana tabacum seedlings by the MKASAA peptide through DNA methylation via the RdDM pathway. AIMS Biophysics, 2022, 9(2): 113-129. doi: 10.3934/biophy.2022011
  • This paper presents a mathematical model for malaria-schistosom-iasis co-infection in order to investigate their synergistic relationship in the presence of treatment. We first analyse the single infection steady states, then investigate the existence and stability of equilibria and then calculate the basic reproduction numbers. Both the single-infection models and the co-infection model exhibit backward bifurcations. We carrying out a sensitivity analysis of the co-infection model and show that schistosomiasis infection may not be associated with an increased risk of malaria. Conversely, malaria infection may be associated with an increased risk of schistosomiasis. Furthermore, we found that effective treatment and prevention of schistosomiasis infection would also assist in the effective control and eradication of malaria. Finally, we apply Pontryagin's Maximum Principle to the model in order to determine optimal strategies for control of both diseases.


    Some proteins are considered cytotoxic because their overexpression causes cellular growth defects. However, toxicity is a relative quantity that can only be defined by comparison. The expression limit of a protein, namely, the protein expression level that results in growth defects, may be considered as an index for protein cytotoxicity [1]. In the model eukaryote Saccharomyces cerevisiae, fluorescent proteins and glycolytic enzymes have the highest expression limits, which means that they are the least cytotoxic in this organism [2], while reports suggest that fluorescent proteins could be cytotoxic in other cells [3],[4]. The cytotoxicity of fluorescent proteins can be increased (i.e., their expression limits are decreased) by adding localization signals or mutations that cause misfolding [5],[6]. The mechanism determining the expression limit of a non-toxic protein is recognized as the “protein burden/cost effect,” in which cellular protein production resources are overloaded by extreme overexpression of the protein [1],[7]. In S. cerevisiae, the limit that triggers protein burden is about 15% of the total protein content [2].

    The nucleotide sequence that corresponds to the N-terminal of the protein strongly affects the expression level of the protein by modulating the efficiency of translation [8]. The N-terminal amino acid sequences of proteins also affect the stability and thus the expression levels of the protein [9]. N-terminal amino acid sequences also affect the cytotoxicity of proteins because they often transmit localization signals [5]. Although there are several systematic analyses revealing the relationship between amino acid sequences and expression levels [9],[10], predicting the effect of the N-terminal sequence on the expression level and cytotoxicity of the protein remains challenging.

    In this study, the effects of the N-terminal sequence on the expression level and cytotoxicity of green fluorescent protein (GFP) were systematically analyzed. Random 30 nucleotides (i.e., ten amino acids) were inserted into the coding sequence for the N-terminal region of GFP. This length was chosen because peptide sequences shorter than ten amino acid length could be linkers connecting different proteins [11], protein localization signals [12][14], and protein tags [15],[16] used for protein engineering. GFP consists of a beta-barrel structure as its main tertiary structure [17], and additional N-terminal sequences should be located outside of the main structure.

    The effects of the sequence on the expression level and GFP cytotoxicity were assessed using the expression per gene (EPG) and gene copy number limit (CNL) relationship. Several sequences that would decrease the expression level or increase the cytotoxicity were isolated. The lower expression level was associated with codon optimality, and cytotoxicity was increased by the existence of hydrophobic or cysteine residues.

    The S. cerevisiae strain BY4741 (MATa his3Δ1 leu2Δ0 met15Δ0 ura3Δ0) [18] was used in all experiments. pTOW40836-TDH3pro-yEGFP (Figure 1A) [19] was used as the base plasmid. The cultivation and transformation of yeast were performed as previously described [20]. Saccharomyces cerevisiae cells were grown in synthetic complete (SC) medium without uracil (−Ura) or without leucine and uracil (−Leu/Ura) conditions. The SC medium was made using the yeast nitrogen base (4027–512, MP Biomedicals), glucose (16806–25, Nacalai tesque), the –His/–Leu/–Ura drop-out mix (630423, Clontech), histidine (18116–92, Nacalai tesque), and leucine (20327–75, Nacalai tesque).

    Figure 1.  Experimental framework isolating N10 sequences affecting the expression level and cytotoxicity of N10-GFPs. (A) Plasmid construct used in this study. Thirty random nucleotides (ten amino acids) were inserted into the N-terminal of yEGFP. Constructed N10-GFPs were expressed under the control of the TDH3 promoter (TDH3pro). The base plasmid used was a multicopy plasmid pTOW40836, whose copy number may be controlled by culture conditions. (B) The experimental procedure of gTOW, with the principle of how copy number limit (CNL) of the target gene is determined by the gTOW. The detail of the figure is explained in the main text.

    Two DNA fragments were amplified with the primers OHM1002 (5′-CATTTTGTTTGTTTATGTGT-3′)–OSBI851 (5′-ATGAGTATTCAACATTTCCG-3′) and OHM600 (5′-ATGTCTAAAGGTGAAGAATT-3′)–OSBI850 (5′-TTACCAATGCTTAATCAGTG-3′) using pTOW40836-TDH3pro-yEGFP as a template. The N10 DNA fragment was prepared via single-cycle amplification of the template DNA (OHM1000: 5′-TAAACACACATAAACAAACAAAATGNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNTCTAAAGGTGAAGAATTATTCACTG-3′) using the primer OHM1001 (5′-CAGTGAATAATTCTTCACCTTTAGA-3′). The DNA fragments were amplified using KOD-Plus-NEO (KOD-401, TOYOBO). All three DNA fragments were simultaneously introduced into BY4741 cells to join them, and the transformants were cultured on an SC–Ura plate. The N10-GFP plasmids were constructed by homologous recombination-based cloning in vivo [21]. Randomly-chosen 94 colonies of the transformants were picked-up and used for further analysis.

    Each yeast strain was cultivated in 200 µL SC–Ura media in each well of a 96-well plate (161093, Thermo) at 30 °C for overnight to prepare the preculture. The precultured cells (5 µL each) were transferred into 200 µL of SC–Ura or SC–Leu/Ura medium in 96-well plates, and grown at 30 °C without agitation. The cell density (absorbance 595 nm) and GFP fluorescence (485 nm excitation – 535 emission) were measured every 30 min for 50 hours using the Infinite F200 microplate reader (TECAN). The maximum growth rate (min −1) during the 270 min interval was calculated from the cell density measurements, and the maximum GFP fluorescence (AU) was obtained from the GFP fluorescence measurements. After cell cultivation, cells were collected, and the plasmid copy number in the cell (copies per cell) was measured via real-time PCR as described previously [22]. The expression per gene (EPG, AU) of each N10-GFP was calculated by dividing the maximum GFP fluorescence in SC–Ura by the plasmid copy number in SC–Ura. The plasmid copy number in SC–Leu/Ura was considered as the copy number limit (CNL) of each N10-GFP as the reason described in Results.

    The folding energy of mRNA ΔG (kcal/mol) containing the N10 sequence (from −4 nt to +37 nt) was calculated using the mfold web server (http://unafold.rna.albany.edu/?q=mfold). Codon optimality was calculated as the tRNA adaptation index (tAI) using the stALcalc (http://tau-tai.azurewebsites.net/). The existence of a localization signal was predicted using DeepLoc-1.0 (http://www.cbs.dtu.dk/services/DeepLoc/). The GRAVY value was calculated using the GRAVY calculator (http://www.gravy-calculator.de/). The distribution of expected GRAVY values in the 10 amino acid sequences translated from random 30 nucleotide sequences (mean = −0.254 and standard deviation = 0.950) were estimated from the 100,000 randomly generated sequences. The existence of the transmembrane domain (TM) and secondary structure in the N-terminal 30 a.a. were predicted using TMHMM server v. 2.0 (http://www.cbs.dtu.dk/services/TMHMM/) and Jpred4 (http://www.compbio.dundee.ac.uk/jpred/), respectively. The statistical analyses were performed using Microsoft Excel (version 16.36).

    Figure 1A shows the plasmid used to assess the effects of the N-terminal sequences on the expression level and GFP cytotoxicity. A yeast codon-optimized enhanced GFP (yEGFP3) [23] was the GFP used. After the N-terminal methionine, random 30 nucleotides encoding 10 amino acids (N10 sequences) were inserted by homologous recombination-based cloning in S. cerevisiae. The resulting N10-GFPs were expressed under the control of the strong TDH3 promoter (TDH3pro). pTOW40836 was used as a backbone plasmid. The copy number of pTOW40836 in S. cerevisiae cells cultured without leucine and uracil (−Leu/Ura) reflects the CNL of the gene. Under such conditions, the overexpressed gene product inhibits the growth of the cell owing to the genetic tug-of-war (gTOW) effect [22],[24].

    The principle of how CNL is determined by the gTOW is explained as follows (Figure 1B). A target gene is cloned into a plasmid for gTOW (Step 1 in Figure 1B), and is used to transform an ura3Δ leu2Δ strain (Step 2 in Figure 1B). In this study, we used TDH3pro_N10-GFP, pTOW40836, and BY4741 as the target, the plasmid for gTOW, and the host strain, respectively. Because of 2µ ORI on the plasmid, the copy number of the plasmid becomes multicopy and diverse among the population. When the cells are transferred into the –Leu/Ura media, the plasmid copy number becomes >100 copies per cell (Step 3 in Figure 1B). This is performed due to leu2-89 on the plasmid, which is an allele of LEU2 harboring a large deletion within its promoter and thus weakly expressed. A higher copy number of leu2-89 is required to fulfill the demand of Leu2 protein to synthesize an essential nutrient leucine (red line in Figure 1B). However, if the target gene has the expression limit to cause cellular defects, the gene copy number and thus the plasmid copy number should be less than the limit, because overexpression of the target beyond the limit causes growth defects (blue line in Figure 1B). The resulting tug-of-war of selection biases between leu2-89 and the target gene (blue and red arrows in Figure 1B) determines the plasmid's copy number, which is determined by the real-time PCR (Step4 in Figure 1B). Because the bias of leu2-89 is always the same, the copy number should reflect the copy number limit (CNL) of the target gene.

    Our framework for identifying N10 sequences affecting the expression level and GFP cytotoxicity is explained as follows. The GFP EPG and CNL of the gene should be inversely proportional because the expression limit of GFP (EPG × CNL) is constant; for example, if there are three genes encoding the same GFP, but their EPGs are low, medium, and high, their CNLs should be high, medium, and low, respectively (Figure 2A). The inversely proportional curve is considered as the frontier line in the EGP–CNL space made by the expression limit of GFP (Figure 2B). If an N10 sequence simply changes the EPG of the N10-GFPs, but not its cytotoxicity and expression limit, the EPG–CNL relationship should be on the same frontier line. Conversely, if an N10 sequence increases the cytotoxicity, thus reducing the expression limit of N10-GFPs, the EPG–CNL relationship should shift to the left (Figure 2B, white arrow). Therefore, N10 sequences, which may alter the expression level and GFP cytotoxicity, can be systematically isolated by analyzing the EPG–CNL relationship.

    Figure 2.  The GFP expression level per gene (EPG) and gene copy number limit (CNL) are inversely proportional. See the main text for the details.

    Overall, 94 yeast colonies containing the N10-GFP plasmid were randomly isolated. Then, their growth rates, level of fluorescence, and plasmid copy numbers without uracil (−Ura) and −Leu/Ura conditions (Table S1) were measured. If analyzed N10-GFPs are under the gTOW regime described above, the growth rates and copy numbers in –Leu/Ura conditions should be on the linear red line created by the relationship between the growth rate and copy number of leu2-89 (Figure 1B). The relationship between the growth rates and copy numbers among the 94 isolates was r = 0.82 (Table S1), suggesting that they are under the gTOW regime, and thus the plasmid copy number under −Leu/Ura conditions was considered as the gene CNL. The GFP EPG was calculated by dividing the level of GFP fluorescence by the plasmid copy number. To calculate the EPG, the level of GFP fluorescence and plasmid copy number obtained under −Ura conditions were used, wherein GFP overexpression produces less growth defects. As shown in Figure 3A, the 94 isolates exhibited various EPG–CNL relationships, which indicates that the N10 sequences modulated the expression levels and toxicities of GFP. Importantly, the EPG–CNL relationships created a frontier line (dotted line in Figure 3A) as it is theoretically predicted in Figure 2B. Most of the N10-GFPs were located around the frontier line, suggesting that most N10 sequences just modulate the expression level of GFP but not the cytotoxicity. Sixteen representative isolates were selected (Figure 3A), and the N10 sequences of each were determined. Seven isolates (#6, #44, #2, #20, #90, #1, and #85) were on the frontier line, but their locations were different, which suggests that their N10 sequences modulated the expression levels of GFP. Four isolates (#15, #75, #80, and #82) were located inside the frontier line, which suggests that their N10 sequences were cytotoxic. Three toxic non-fluorescent isolates (#51, #52, and #64) and two non-toxic non-fluorescent isolates (#55 and #74) were also analyzed.

    Figure 3.  EPG–CNL relationships of isolated N10-GFPs. (A) Scatterplot between EPGs and CNLs of isolated N10-GFPs. Colored and numbered isolates were studied further. The dotted line indicates the predicted protein burden frontier line. (B) EPGs and CNLs of representative N10-GFP isolates. For each sample, the data from a single experiment are shown.

    Table 1 presents the sequencing results. First, toxic non-fluorescent isolates (#51, #52, and #64), but not N10-GFPs, contained the TDH3 gene. Therefore, the plasmids were probably constructed through homologous recombination between the base plasmid and TDH3 within the S. cerevisiae genome. Because Tdh3 has a lower expression limit than GFP [19], isolates expressing Tdh3 should have exhibited EPG–CNL relationships that differed from those of GFP. Second, non-toxic non-fluorescent isolates (#55 and #74) had nonsense codons within their N10 sequences, indicating that the genes do not encode any proteins. Third, isolates showing the highest EPGs (#6 and #44) did not contain the N10 sequence but rather GFP alone; these isolates were probably constructed by self-ligation of the base plasmid. Because #6 and #44 are located on the top of the frontier line, inserting an N10 sequence seemed to reduce the expression level or increase the GFP cytotoxicity. The other nine isolates contained various N10 sequences; among these, five (#2, #20, #90, #1, and #85) decreased the expression level, whereas four (#15, #75, #70, and #82) increased cytotoxicity as shown in Figure 3A. Figure 3B shows the CNLs and EPGs of the representative isolates.

    Table 1.  Isolated N10 sequences.
    Clone # Category Nucleotide sequence Amino acid sequence
    #6 High nontoxic without N10 sequence
    #44 High nontoxic without N10 sequence
    #2 middle nontoxic ATGAGCATGATACGGGAGGTTGATGAATGTTTA MSMIREVDECL
    #20 middle nontoxic ATG TGTACTGGACTCGATGTTGGCCCCGTCCAA MCTGLDVGPVQ
    #90 middle nontoxic ATGTCCGCAGCTGCGATCGGGTGCGGGTATGTA MSAAAIGCGYV
    #1 Low nontoxic ATGGGATGGCTGCCGAAAGGGTTCTTAGAAGAG MGWLPKGFLEE
    #85 Low nontoxic ATGGTGGTGGTGGGTTGGCAGCCCCCGGGTGGG MVVVGWQPPGG
    #15 toxic ATGTGCTTCGGGGTAATGAGGGGTATCTGGCTG MCFGVMRGIWL
    #75 toxic ATGTGCACTGAGCTGTGTTTTCTTCTCGGTGTC MCTELCFLLGV
    #80 toxic ATGCCGGTGATTGGGTATTTTGATGGCTTGTTA MPVIGYFDGLL
    #82 toxic ATGGATCCATTGTGCCAGTTGATGGGTGTGCTA MDPLCQLMGVL
    #51 No expression, toxic without N10 sequence, TDH3
    #52 no expression, toxic without N10 sequence, TDH3
    #64 no expression, toxic without N10 sequence, TDH3
    #55 no expression, nontoxic ATGTCCTCTGCTCATTCGTTGCCGGTATAGGGT MSSAHSLPV*G
    #74 no expression, nontoxic ATGTAGAAATTTACGGTTGCGATCGCCAGAGTG M*KFTVAIARV

     | Show Table
    DownLoad: CSV

    To reveal the sequence property that determines the expression level of N10-GFPs, the mRNA folding energy [25] and codon optimality (tRNA adaptation index; tAI) [26] were calculated (Table S2). The existence of localization signals [27] was also predicted. Although no significantly associated values with the expression level and cytotoxicity were noted, the tAIs of isolated N10 sequences were found to be lower than those of GFP (0.320). Their codon usage was thus optimized to that of the highly expressed TDH3. Codon optimization altered the EPG–CNL relationship profiles of N10-GFPs (Figure 4). Most importantly, the optimization increased the EPGs of isolates with low expression levels of N10-GFPs (#1 and #85), which suggests that codon optimality is responsible for their low EPGs. Thus, codon optimality of the N10 sequences seemed to be the most important contributing factor to the expression levels of N10-GFPs.

    Figure 4.  Effect of codon optimization on the EPGs and CNLs of N10-GFPs. Scatterplot between EPGs and CNLs of original and codon-optimized N10-GFPs. Arrows indicate the changes by codon optimization. The dotted line indicates the predicted protein burden frontier line the same as shown in Figure 3A. Averages of four biological replicates are shown.

    Next, the relationship between the amino acid composition in the N10 sequence and N10-GFP cytotoxicity was analyzed. First, the grand average of hydropathy (GRAVY) value of the amino acid sequence [28] (Table S2) was calculated, which has a positive relationship with hydrophobicity. The toxic N10 sequences had significantly higher GRAVY values than those of the non-toxic N10 sequences (p = 0.028, Student t-test, Figure 5A). The probability of obtaining four N10 sequences with GRAVY values of >0.696 is <0.0006, which suggests that hydrophobic residues contribute to the cytotoxicity of the N10 sequences. Next, the frequency of each amino acid in the four toxic N10 sequences (#15, #75, #80, and #82) was calculated and compared with the expected frequencies in the codon table (Figure 5B). The quantities of cysteine, methionine, and tryptophan residues were over two-fold of the expected frequencies. The p-value of the appearance of cysteine was 0.006 according to the results of the chi-square test, which was slightly insignificant under the false discovery rate <0.05 (p < 0.003). Therefore, we concluded that cysteine residues were over-represented in toxic N10-GFPs. Structural analyses of the N-terminal 30 a.a. containing N10 sequences and the N-terminal of GFP were also performed (Table S2). A transmembrane domain prediction by TMHMM (v2.0) [29] did not give any transmembrane domain in the sequences, and a secondary structure prediction by Jpred4 [30] gave some structures, but there was no structural consistency within the toxic and non-toxic N10 sequences.

    Figure 5.  Characteristics of toxic N10 sequences. (A) GRAVY values of the non-toxic and toxic N10 sequences. The p-value of the Student t-test between the GRAVY scores of the non-toxic and toxic N10 sequences is presented. (B) Frequency of appearance of each amino acid in the toxic N10 sequences. The frequency was calculated as the appeared over the expected frequency. The p-value of the chi-square test is shown.

    The nucleotide sequence encoding the N-terminal of a protein plays a substantial role in protein production because it is where ribosomes are assembled and translation is initiated. The secondary structure of mRNA and codon optimality affect the initiation of translation [8]. In this study, it was confirmed that codon optimality was a major factor in reducing the expression level of GFP, although mRNA stability did not affect the expression level. A question that was raised during this study focused on the mechanism causing cytotoxicity of N10 sequences. A possible explanation for this phenomenon is the artificial creation of localization signals; this occurs because the addition of localization signals at the N-terminal of GFP increases its cytotoxicity [5]. However, three out of the four toxic N10 sequences did not have any signature localization signal. Moreover, no localization of toxic N10-GFPs (data not shown) was observed. Agammaegative sequences also increase the cytotoxicity of proteins [31]; however, no agammaegation of toxic N10-GFPs (data not shown) was noted. Although a small subset of toxic N10 sequences was analyzed here, they showed statistically significant concentration of hydrophobic residues and a preferential appearance of cysteine residues. Overexpression of cysteine-containing glycolytic enzymes constitutes the agammaegates connected by the S–S bond [2],[32], and removal of cysteine from Tpi1 increases its expression limit [2]. Therefore, the existence of unimportant cysteines might increase the cytotoxicity of highly expressed proteins regardless of their positions. Hydrophobic residues seemed to increase the GFP cytotoxicity. This may be owing to non-specific hydrophobic interactions, and therefore, they might not be preferred in the outside of the main tertiary structure of proteins. Another possible mechanism to cause cytotoxicity is the perturbation of translation. Sequences rich in proline and positively-charged residues cause slowdown in translational elongation [33]. However, this might not be the case in this analysis because any enrichment of proline or positively-charged residues (lysine or arginine) was observed in the toxic N10 sequences (Figure 5B).

    [1] [ B. M. Adams,H. T. Banks,H. Kwon,H. T. Tran, Dynamic multidrug therapies for HIV: Optimal and STI control approaches, Mathematical Biosciences and Engineering, 1 (2004): 223-241.
    [2] [ F. B. Agusto, Optimal chemoprophylaxis and treatment control strategies of a tuberculosis transmission model, World Journal of Modelling and Simulation, 5 (2009): 163-173.
    [3] [ F. B. Agusto,K. O. Okosun, Optimal seasonal biocontrol for Eichhornia crassipes, International Journal of Biomathematics, 3 (2010): 383-397.
    [4] [ R. M. Anderson and R. M. May, Infectious Diseases of Humans: Dynamics and Control, Oxford University Press, 1991, Oxford.
    [5] [ K. W. Blayneh,Y. Cao,H. D. Kwon, Optimal control of vector-borne diseases: Treatment and Prevention, Discrete and Continuous Dynamical Systems Series B, 11 (2009): 587-611.
    [6] [ J. G. Breman,M. S. Alilio,A. Mills, Conquering the intolerable burden of malaria: What's new, what's needed: A summary, Am. J. Trop. Med. Hyg., 71 (2004): 1-15.
    [7] [ C. Castillo-Chavez,B. Song, Dynamical model of tuberculosis and their applications, Math. Biosci. Eng., 1 (2004): 361-404.
    [8] [ Z. Chen,L. Zou,D. Shen,W. Zhang,S. Ruan, Mathematical modelling and control of Schistosomiasis in Hubei Province, China, Acta Tropica, 115 (2010): 119-125.
    [9] [ E. T. Chiyaka,G. Magombedze,L. Mutimbu, Modelling within host parasite dynamics of schistosomiasis, Comp. Math. Meth. Med., 11 (2010): 255-280.
    [10] [ J. A. Clennon,C. G. King,E. M. Muchiri,U. Kitron, Hydrological modelling of snail dispersal patterns in Msambweni, Kenya and potential resurgence of Schistosoma haematobium transmission, Parasitology, 134 (2007): 683-693.
    [11] [ S. Doumbo,T. M. Tran,J. Sangala,S. Li,D. Doumtabe, Co-infection of long-term carriers of Plasmodium falciparum with Schistosoma haematobium enhances protection from febrile malaria: A prospective cohort study in Mali, PLoS Negl. Trop. Dis., 8 (2014): e3154.
    [12] [ M. Finkel, Malaria: Stopping a Global Killer, National Geographic, July 2007.
    [13] [ Z. Feng,A. Eppert,F. A. Milner,D. J. Minchella, Estimation of parameters governing the transmission dynamics of schistosomes, Appl. Math. Lett., 17 (2004): 1105-1112.
    [14] [ W. H. Fleming,R. W. Rishel, null, Deterministic and Stochastic Optimal Control, Springer Verlag, New York, 1975.
    [15] [ J. H. Ge,S. Q. Zhang,T. P. Wang,G. Zhang,C. Tao,D. Lu,Q. Wang,W. Wu, Effects of flood on the prevalence of schistosomiasis in Anhui province in 1998, Journal of Tropical Diseases and Parasitology, 2 (2004): 131-134.
    [16] [ P. J. Hotez,D. H. Molyneux,A. Fenwick,E. Ottesen, Ehrlich and S. Sachs et al., Incorporating a rapid-impact package for neglected tropical diseases with programs for HIV/AIDS, tuberculosis, and malaria, PLoS Med., 3 (2006): e102.
    [17] [ M. Y. Hyun, Comparison between schistosomiasis transmission modelings considering acquired immunity and age-structured contact pattern with infested water, Mathematical Biosciences, 184 (2003): 1-26.
    [18] [ H. R. Joshi, Optimal control of an HIV immunology model, Optimal Control Applications in Mathematics, 23 (2002): 199-213.
    [19] [ A. Kealey,R. J. Smith?, Neglected Tropical Diseases: Infection, modelling and control, J. Health Care for the Poor and Underserved, 21 (2010): 53-69.
    [20] [ J. Keiser,J. Utzinger,M. Caldas de Castro,T. A. Smith,M. Tanner,B. Singer, Urbanization in sub-Saharan Africa and implication for malaria control, Am. J. Trop. Med. Hyg., 71 (2004): 118-127.
    [21] [ D. Kirschner,S. Lenhart,S. Serbin, Optimal Control of the Chemotherapy of HIV, J. Math. Biol., 35 (1997): 775-792.
    [22] [ J. C. Koella,R. Anita, Epidemiological models for the spread of anti-malaria resistance, Malaria Journal, 2 (2003): p3.
    [23] [ C. M. Kribs-Zaleta,J. X. Velasco-Hernandez, A simple vaccination model with multiple endemic states, Math. Biosci., 164 (2000): 183-201.
    [24] [ V. Lakshmikantham, S. Leela and A. A. Martynyuk, Stability Analysis of Nonlinear Systems, Marcel Dekker, New York and Basel, 1989.
    [25] [ S. Lenhart and J. T. Workman, Control Applied to Biological Models, Chapman and Hall, London, 2007.
    [26] [ J. Li, D. Blakeley and R. J. Smith?, The failure of R0, Comp. Math. Meth. Med. , 2011 (2011), Article ID 527610, 17pp.
    [27] [ G. Li,Z. Jin, Global stability of a SEIR epidemic model with infectious force in latent, infected and immune period, Chaos, Solutions and Fractals, 25 (2005): 1177-1184.
    [28] [ Q. Longxing, J. Cui, T. Huang, F. Ye and L. Jiang, Mathematical model of schistosomiasis under flood in Anhui province Abstract and Applied Analysis, 2014(2014), Article ID 972189, 7pp.
    [29] [ A. D. Lopez,C. D. Mathers,M. Ezzati,D. T. Jamison,C. J. Murray, Global and regional burden of disease and risk factors, 2001: Systematic analysis of population health data, Lancet, 367 (2006): 1747-1757.
    [30] [ E. Mtisi,H. Rwezaura,J. M. Tchuenche, A mathematical analysis of malaria and Tuberculosis co-dynamics, Discrete and Continuous Dynamical Systems Series B, 12 (2009): 827-864.
    [31] [ Z. Mukandavire,A. B. Gumel,W. Garira,J. M. Tchuenche, Mathematical analysis of a model for HIV-Malaria co-infection, Mathematical Biosciences and Engineering, 6 (2009): 333-362.
    [32] [ S. Mushayabasa and C. P. Bhunu, Modeling Schistosomiasis and HIV/AIDS co-dynamics, Computational and Mathematical Methods in Medicine, 2011(2011), Article ID 846174, 15pp.
    [33] [ S. Mushayabasa,C. P. Bhunu, Is HIV infection associated with an increased risk for cholera? Insights from mathematical model, Biosystems, 109 (2012): 203-213.
    [34] [ I. S. Nikolaos,K. Dietz,D. Schenzle, Analysis of a model for the Pathogenesis of AIDS, Mathematical Biosciences, 145 (1997): 27-46.
    [35] [ K. O. Okosun,R. Ouifki,N. Marcus, Optimal control analysis of a malaria disease transmission model that includes treatment and vaccination with waning immunity, BioSystems, 106 (2011): 136-145.
    [36] [ K. O. Okosun,O. D. Makinde, Optimal control analysis of malaria in the presence of non-linear incidence rate, Appl. Comput. Math., 12 (2013): 20-32.
    [37] [ K. O. Okosun,O. D. Makinde, A co-infection model of malaria and cholera diseases with optimal control, Mathematical Biosciences, 258 (2014): 19-32.
    [38] [ L. S. Pontryagin, V. G. Boltyanskii, R. V. Gamkrelidze and E. F. Mishchenko, The Mathematical Theory of Optimal Processes, Wiley, New York, 1962.
    [39] [ R. Ross, The Prevention of Malaria, Murray, London, 1911.
    [40] [ P. Salgame,G. S. Yap,W. C. Gause, Effect of helminth-induced immunity on infections with microbial pathogens, Nature Immunology, 14 (2013): 1118-1126.
    [41] [ A. A. Semenya,J. S. Sullivan,J. W. Barnwell,W. E. Secor, Schistosoma mansoni Infection Impairs Antimalaria Treatment and Immune Responses of Rhesus Macaques Infected with Mosquito-Borne Plasmodium coatneyi, Infection and Immunity, 80 (2012): 3821-3827.
    [42] [ K. D. Silué,G. Raso,A. Yapi,P. Vounatsou,M. Tanner,E. Ńgoran,J. Utzinger, Spatially-explicit risk profiling of Plasmodium falciparum infections at a small scale: A geostatistical modelling approach, Malaria J., 7 (2008): p111.
    [43] [ R. J. Smith? and S. D. Hove-Musekwa, Determining effective spraying periods to control malaria via indoor residual spraying in sub-saharan Africa Journal of Applied Mathematics and Decision Sciences, 2008(2008), Article ID 745463, 19pp.
    [44] [ R. W. Snow,C. A. Guerra,A. M. Noor,H. Y. Myint,S. I. Hay, The global distribution of clinical episodes of Plasmodium falciparum malaria, Nature, 434 (2005): 214-217.
    [45] [ R. C. Spear,A. Hubbard,S. Liang,E. Seto, Disease transmission models for public health decision making: Toward an approach for designing intervention strategies for Schistosomiasis japonica, Environ. Health Perspect., 10 (2002): 907-915.
    [46] [ P. van den Driessche,J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002): 29-48.
    [47] [ R. B. Yapi,E. Hürlimann,C. A. Houngbedji,P. B. Ndri,K. D. Silué, Infection and Co-infection with Helminths and Plasmodium among School Children in Côte d'Ivoire: Results from a National Cross-Sectional Survey, PLoS Negl. Trop. Dis., 8 (2014): e2913.
    [48] [ X. N. Zhou,J. G. Guo,X. H. Wu, Epidemiology of schistosomiasis in the people's republic of China, 2004, Emerging Infectious Diseases, 13 (2007): 1470-1476.
  • This article has been cited by:

    1. Young Jin Lee, Joong Seob Lee, Olatunji Ajiteru, Ok Joo Lee, Ji Seung Lee, Hanna Lee, Seong Wan Kim, Jong Woo Park, Kee Young Kim, Kyu Young Choi, Heesun Hong, Tipu Sultan, Soon Hee Kim, Chan Hum Park, Biocompatible fluorescent silk fibroin bioink for digital light processing 3D printing, 2022, 213, 01418130, 317, 10.1016/j.ijbiomac.2022.05.123
    2. Shotaro Namba, Hisaaki Kato, Shuji Shigenobu, Takashi Makino, Hisao Moriya, C Hoffman, Massive expression of cysteine-containing proteins causes abnormal elongation of yeast cells by perturbing the proteasome, 2022, 12, 2160-1836, 10.1093/g3journal/jkac106
    3. Andreea Veronica Dediu Botezatu, Gabriela Elena Bahrim, Claudia Veronica Ungureanu, Anna Cazanevscaia Busuioc, Bianca Furdui, Rodica Mihaela Dinica, Green “one-pot” fluorescent bis-indolizine synthesis with whole-cell plant biocatalysis, 2023, 12, 2191-9550, 10.1515/gps-2023-0046
    4. Belen Torrado, Bruno Pannunzio, Leonel Malacrida, Michelle A. Digman, Fluorescence lifetime imaging microscopy, 2024, 4, 2662-8449, 10.1038/s43586-024-00358-8
    5. Zunlin Yang, Chunfu Zheng, 2025, Chapter 5, 978-1-0716-4614-4, 33, 10.1007/978-1-0716-4615-1_5
  • Reader Comments
  • © 2017 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4289) PDF downloads(740) Cited by(5)

Figures and Tables

Figures(10)  /  Tables(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog