第12卷第6期北華大學學報(自然科學版)Vol. 12 No 62011年12月JOURNAL OF BEIHUA UNIVERSITY( Natural Science)Dec.2011文章編號:10094822(2011)0640621406HⅣ∨病理動力學模型分析蘇雙雙,王凱,夏米希努爾1.新疆大學數(shù)學與系統(tǒng)科學學院新疆烏魯木齊830046;2新疆醫(yī)科大學醫(yī)學工程與技術(shù)學院,新疆烏魯木齊830054)摘要:考慮了CDT細胞的全 logistic自我繁殖,構(gòu)造了具有治療的HV病理研究模型,研究了解的存在性、正性及穩(wěn)定性通過對每個感染T細胞釋放的HⅣ病毒平均數(shù)量N的討論,得到:如果NN灬,有無病平衡點和地方病平衡點兩個平衡點,但此時無病平衡點不穩(wěn)定,地方病平衡點是局部漸近穩(wěn)定的最后通過構(gòu)造 dulac- Benison函數(shù),得到了地方病平衡點全局穩(wěn)定的條件關鍵詞:HV;T細胞;漸近穩(wěn)定;平衡點中圖分類號:0175.12文獻標志碼:AAnalysis of the Dynamics Model of HIv Pathogenesis ModelSU Shuang-shuang, WANG Kai, Xamxinur(1. College of Mathematics and System Sciences, Xin jiang University, Urumqi 830046, China;2. Department of Medical Engineering and Technology, Xinjiang Medical University, Urumgi 830054, ChinaAbstract: HIV virus, which invade the human body, mostly parasitic in the T cells, where they transcribe RNAduplicate, then generate some new ones, at last kill the infected T cells, released from them. According to this wedivide the species in this paper into three kinds: uninfected T cells, infected T cells and HIV. Only uninfected Tcells can be created by body, others can't increase because of external input. Being predator, HIV just add frominfected T cells. We have proved that if NNit the infected steady state is locally asymptotically stable. Both of them are globally asymptotically stable ifone condition is satisfiedKey words HIV; T cells; asymptotically stable; equilibrium1引言自從20世紀80年代艾滋病被發(fā)現(xiàn)至今,已有許多關于HV病理學研究結(jié)果2,尤其是近些年將數(shù)學模型與HV病理學相結(jié)合3.我們更好地理解和治療HV提供了很大的幫助作為威脅人類最強的傳染病之一,對HIV病理研究是十分必要的而數(shù)學模型這一強大的工具,在史上幾次流行病研究中顯示出了其強大的威力眾所周知,HⅣV進入人體后主要攻擊T淋巴細胞(簡稱T細胞),寄宿在T細胞內(nèi),完成RNA反轉(zhuǎn)錄及復制,當產(chǎn)生新的HV病毒后再將T細胞徹底殺死『產(chǎn)中的*6,故對具有收稿日期:20110924CNMHG作者簡介:蘇雙雙(1986-),女,碩士研究生,主要從事常微分方程及其應用研究;王凱(1982-),男,講師博士,主要從事微分方程在生物傳染病、復雜網(wǎng)絡的應用研究622北華大學學報(自然科學版)第12卷藥物治療模型的T細胞及HIV的動力學性質(zhì)進行討論會對艾滋病的治療提供很大幫助.在這些模型中,人體免疫系統(tǒng)中其他對HV病毒具有殺傷力的,如CDT細胞也已考慮進去雖然完全殺死HⅣV是不可能的,但華裔科學家何大一發(fā)明的雞尾酒療法對控制HV病毒攜帶者體內(nèi)的病毒數(shù)量收到了意想不到的效果HIV病毒在人體內(nèi)可以達到一平衡狀態(tài),且可以局部漸近穩(wěn)定,但能否全局穩(wěn)定,我們將進行討論A S Perelson等7討論了如下基本模型:d-T-krv1.1)T(1)=kⅣV-δ7,V(t)=NδT’-cV但對系統(tǒng)(1.1)作者忽略了被感染的T細胞的自我繁殖能力HIV的RNA經(jīng)過反轉(zhuǎn)錄,形成DNA.因此,病毒DNA是通過整合到T細胞的DNA中才完成復制8],從這方面來看被感染T細胞應該也是具有自我繁殖能力的但因感染了HⅣV,治療中有抑制其繁殖的藥物,其活躍性要減弱通過以上分析,我們可做如下模型假設1)只有HⅣV病毒可以感染T細胞,被感染的T細胞不能傳染健康T細胞;2)T細胞在人體中的數(shù)量有最大值;3)基于以上理論及假設,我們可得如下模型7(t)=s-1T+r1T1i(t)=A1rv+r,n(1T+1)-a21(1.2)V(1)=Nu2/-u3V,其中:T,,V分別表示健康T細胞、被感染T細胞和HIV病毒數(shù)量;5表示人體淋巴系統(tǒng)產(chǎn)生健康T細胞的速率;Tm是T細胞在人體中含量的最大值;λ是HV病毒感染健康T細胞的速率;l1,u2,3分別表示健康T細胞感染T細胞和HV病毒的死亡率(u10,l(0)>0,V(0)>0(1.3)2解的正性及有界性分析如果沒有HIV病毒人侵,健康T細胞的動力系統(tǒng)為T=s-1+T1-可見健康T細胞在T時達到穩(wěn)定,其中TT4r且7滿足(u1-r1)T0下面我們研究模型(1.2)的正性定理2.1在初值條件(1.3)下,系統(tǒng)(1.2)的解都是正的證明反證法因為r(0)>0,假設正性不成立,則存在T(t)=0,其中to=inft|t>0,7(t)=0由系統(tǒng)(1.2)中的第1個等式可知T(4)=s>0.而根據(jù)導數(shù)的定義可知T()≤0.矛盾.所以T()>0恒成立中國煤化工現(xiàn)在我們證明(t)>0,V(t)>0恒成立若不成立,則存在CNMHG)I(t1)=0,W(t1)>0,且存在t∈(0,l1),有I(l)>0,V(t)>0,P(t1)≤0成立,而由系統(tǒng)(1.2)的第2個等式有6期蘇雙雙,等:HV病理動力學模型分析4)=(4)4)+6(4)(1-(4+()-(4)=7(44)>0,矛盾;i)I(41)>0,v(t1)=0,且存在t∈(0,),有I(t)>0,(t)>0,v(t1)≤0成立,而v(t1)=Nu2(t1)-u3V(1)>0,矛盾;ⅲ)I(t1)>0,(t1)=0,且存在t∈(0,t),有I(t)>0,v(t)>0,而l(t)=0,W(t)=0是系統(tǒng)(1.2)的一個解,且滿足初始條件l(t)=0,(t)=0,由解的唯一性,對任意t≥0,有l(wèi)(t)=0,V(t)=0恒成立,矛盾下面研究解的有界性定理2.2系統(tǒng)(1.2)滿足初始條件(1.3)的解有界證明由系統(tǒng)(1.2)的第1個等式得T≤-T+7(1-)所以,如果7(0)<70,那么 lim sup7(t)≤T,對所有t>0成立將系統(tǒng)(1.2)的前兩個等式相加得a(r+n=s-a1r+nr(1-T+1+n1(1-x+)-n1/≤r2(T+D)s+r1(T+1)-(T+D)r2(T+nT(n1-u1)(T+1)+s.(2.1)由式(21)可知,存在M1>0,使得T+/≤M1成立,又因為(t),W(t)常正,則存在M2>0,使得同樣,由系統(tǒng)(1.2)的第3個等式得那么N+(v(0)(2.2)通過不等式(2.2)可知v(t)恒有正上界(記為M3)3平衡點的存在性系統(tǒng)(1.2)的平衡點滿足的方程為A1v-u1T=0,A,TV +2112=03.1)-n3V=0.顯然系統(tǒng)(1.2)有無病平衡點E0=(7,0,0).由方程組(3.1)的第3個等式,有將式(3.2)代入式(3.1)的第2個等式得λ1L把式(3.1)的第1個等式改寫為A將式(3.2)和式(3.3)代入式(3.4)得中國煤化工(B1+B2V)(B3+B4V)CNMHG或B2B,V+(B,B4+B, B,)V+B,B624北華大學學報(自然科學版式中l(wèi)2r2lruat+ Tr1(u2-r2)u3TINAl1 TA,Nu2Tma(,-r2)u,+A, Nu2TB(NA,u,Tmas)系統(tǒng)(1.2)有正解當且僅當?shù)仁?3.5)有一個正解V使得B1+B2V>0.簡單起見,令-(B2B3+B1B)±√(B2B3+B1B)2-4B2B4(B1B3-s)2B.B如果u2>r2,B1B3-5<0等價于s(NAju2Tmax-r2u3)-(u2-r2)u3Tmax (u, -T(NAu2Tngx-r2u3)-(u2-r2)u3Tmar ,>0上述不等式的解為N>N(u2-r2)3r2A,u2To A,u2T(u2-r2)T其中:Nm=A1u2T0(r1T07T2s u,T2分3種情況討論:情形1.如果N>N_,那么B1>0,B2>0,B4>0,因而B2B4>0,B1B3-5<0,系統(tǒng)(3.5)有一個正解v,B1+B2V,>0,所以T=B1+B2V,=mV,此時V=V,所以系統(tǒng)(1.2)只有一個正平衡點E=(T”,I,V).情形2如果、20,B2>0,B,>0,進而BB4>0.當n10;當r1>u1時,由N≤Nn,有MA1u2mx-r2u3≤所以(ur-T,(NA u2T(u2 -r2)u3Tmax(u-r.)進而(u1-r1)(NA1a27十r;u2(2-1){(-x)La s所以B3>0.故有B1B3+B2B4>0,系統(tǒng)(1.2)無正平衡點情形3.如果0Nm時系統(tǒng)(1.2)只有一正平衡點;當NNa時系統(tǒng)(1.2)只有一正平衡點;當NN-,那么B1B3-s<0,B2>0,B4>0,從而B2B4>0,所以等式(3.5)只有一正根V.,B1+B,v。BB4-B2B+√(BB,-BB1)+4BBs、B1B、-B2B+|B1B1-B,B≥0,所以,系統(tǒng)2B4(1.2)只有一正平衡點E=(T,,V),此時V=Vⅱ)如科C0,B4>0.但B1<0,因為B1B3-<0,所以B3<0,得到兩個正根V和V,,但是B1+B2V=B,B,-B,B,+V(B,B4-B2 B,)+4B, B,BB1B4-B2B3+|B1B4-B2B3B,B4-B2 B,+V(B, BA-B2B3)+4B,BA2B≤0,B1+B2V,2B4B,B -B2 B,+B, B 4 -B, B3≥0,所以系統(tǒng)(1.2)只有一正平衡點E’=(T,,V),此時v=Vnuⅲ)如果Nan0,B1>0,B2<0,進而B3>0.當B4>0時,B2B4<0,V<0,v,>0,又B1+B2VBB-B2B+√(BB-B2B)+4BBs、BB-B+|BB1-BB1≥0,所以,系統(tǒng)(1.2)只有一正平衡點E=(T,',V),此時v=V;當B4<0時,B1+B2V=B,BA-B,B,+V(B, B4-B2B3)2+4B, B,s B, B,-B2B,+B,B,-B2B3≤0,B1BeV2B42BBB-BB+√(BA-B2)+4BB.、B-B1+1B1-B.≥0,所以,系統(tǒng)(12)只有2B42B4一正平衡點E*=(T”,,V),此時v=Viv)如果N0,B2<0,B4<0,那么B2B4>0,只有一正根V,但是B1+B, v-B, B4-B, B,+V(B, B4 -B2 B,)2+4B, BAS B,B4-B2B,+ B Bs-B2B3≤0,所以系統(tǒng)2B42B4(1.2)無正平衡點定理31如果NNm,系統(tǒng)(1.2)有無病平衡點E=(70,0,0)和地方病平衡點E=(T”,,V)4平衡點的穩(wěn)定性分析令E=(行,1,V)為任意平衡點,系統(tǒng)(1.2)在點E的線性化矩陣為2r1T+r27TT入λ1TA,V該矩陣的特征方程為A+Q1T 1, TT-n,vA+Q2-A=中國煤化工N4 +uCNMHG2r, Tr,T+2rI其中:Q1=u1-r1+即北華大學學報(自然科學版)第12卷(A+1-2+2x2)j2+cM+c1-c)=0,(4.1)其中:Co=a2{l1~n2)+(1-)+l3>0,C1=1-r2)T+2r27顯然A1=-u,2rT=-(u<0,方程(4.1)的其余特征根依賴于A2+C0A-C1-C2=0.如果N0,只有負根,E局部漸近穩(wěn)定;如果N>Nm,那么C1-C2<0,E。是鞍點,不穩(wěn)定因此,得到無病平衡點的局部穩(wěn)定性結(jié)果定理4.1如果NNom,E不穩(wěn)定.在點E處線性化后系統(tǒng)(1.2)的特征方程為A3+p1A2+P2A+p3=0,其中:rrT/AITV T27)+y(÷+)+(A,V令Q=P2-內(nèi)=(+r1T”A1Tv,r2IA,T"V+41,Ar4pA1V·」,得到E’的局部穩(wěn)定性結(jié)果定理42當N>Nm時,如果Q>0且7-A1V>0,E局部漸近穩(wěn)定最后我們用 Dulas準則來判斷系統(tǒng)有無周期解及極限環(huán).令B(T,1p)=7,可得a(BT+(BD)+0(Bn)=-s-+2所以,如果a2+u3>r2,則系統(tǒng)(1.2)無周期解定理4.3如果u2+u3>r2,系統(tǒng)(1.2)無周期解,則當NNn時,E·全局漸近穩(wěn)定參考文獻[1] A Perelson, A Neumann, M Markoitz, et al. HIV, Dynamics in Vivo: Virion Clearance Rate, Infected Cell Life-span, and ViralGeneration Time[J]. Science, 1996, 271: 1582-1586[2]R V Culshaw, SG Ruan. A Differential Equation Model for HIV Infection of CDA T-cells[ J]. Math Biosci, 2000(165): 27-39[3] P De Leeheer, H L Smith. Virus Dynamics: A Global Analysis[J]. SIAM I Appl Math, 2003 (631): 1313-1327[4 Patrick W Nelson, James D Murray, Alan S Perelson. A Model of HIV-1 Pathogenesis that Includes an Intracellular Delay[J]Mathematical Biosciences, 2000, 163(2): 201-215[5] LC Wang, M Y Li Mathematical Analysis of the Global Dynamics of a Model for HIV Infected of CD4 T-cells[J]. Appl Mathimput,2007(189):1331-13406] William A,O' Brien M d.艾滋病的治療與控制[J.中國性病艾滋病防中國煤化工[7 A Perelson, P W Nelson. HIV, Dynamics in Vivo[J]. SIAM REV, 1999CNMHG[8]池莉,水波,張永容艾滋病發(fā)病機制中細胞凋亡理論[J中國實驗動物學余志,0,0(4):449【責任編輯:伍林】