Global Asymptotic Stability and Hopf Bifurcation for a Blood Cell Production Model

Fabien Crauste

Description

Niveau: Supérieur, Doctorat, Bac+8
Global Asymptotic Stability and Hopf Bifurcation for a Blood Cell Production Model? Fabien Crauste Year 2005 Laboratoire de Mathematiques Appliquees, UMR 5142, Universite de Pau et des Pays de l'Adour, Avenue de l'universite, 64000 Pau, France. ANUBIS project, INRIA–Futurs E-mail: Abstract We analyze the asymptotic stability of a nonlinear system of two differential equa- tions with delay describing the dynamics of blood cell production. This process takes place in the bone marrow where stem cells differentiate throughout divisions in blood cells. Taking into account an explicit role of the total population of hematopoietic stem cells on the introduction of cells in cycle, we are lead to study a characteristic equation with delay-dependent coefficients. We determine a necessary and sufficient condition for the global stability of the first steady state of our model, which describes the population's dying out, and we obtain the existence of a Hopf bifurcation for the only nontrivial pos- itive steady state, leading to the existence of periodic solutions. These latter are related to dynamical diseases affecting blood cells known for their cyclic nature. Keywords: asymptotic stability, delay differential equations, characteristic equation, delay- dependent coefficients, Hopf bifurcation, blood cell model, stem cells. 1 Introduction Blood cell production process is based upon the differentiation of so-called hematopoietic stem cells, located in the bone marrow.

• blood cell


• production model


• has been


• steady state


• unique solution


• stem cells


• been introduced


• cycle duration


• pos- itive steady


• crauste stability

Sujets

• Blood cell
• Production vehicle
• Has Been
• Steady state
• Célula madre

Informations

• Publié par: mijec
• Nombre de lectures: 24
• Langue: English

Extrait

GlobalAsymptoticStabilityandHopf BifurcationforaBloodCellProductionModel ∗ Fabien Crauste Year2005 LaboratoiredeMathe´matiquesApplique´es,UMR5142, Universite´dePauetdesPaysdel’Adour, Avenuedel’universite´,64000Pau,France. ANUBISproject,INRIA–Futurs E-mail:[e-mail] Abstract Weanalyzetheasymptoticstabilityofanonlinearsystemoftwodifferentialequa- tionswithdelaydescribingthedynamicsofbloodcellproduction.Thisprocesstakes placeinthebonemarrowwherestemcellsdifferentiatethroughoutdivisionsinblood cells.Takingintoaccountanexplicitroleofthetotalpopulationofhematopoieticstem cellsontheintroductionofcellsincycle,weareleadtostudyacharacteristicequation withdelay-dependentcoefficients.Wedetermineanecessaryandsufficientconditionfor theglobalstabilityofthefirststeadystateofourmodel,whichdescribesthepopulation’s dyingout,andweobtaintheexistenceofaHopfbifurcationfortheonlynontrivialpos- itivesteadystate,leadingtotheexistenceofperiodicsolutions.Theselatterarerelated todynamicaldiseasesaffectingbloodcellsknownfortheircyclicnature. Keywords:asymptoticstability,delaydifferentialequations,characteristicequation,delay- dependentcoefficients,Hopfbifurcation,bloodcellmodel,stemcells. 1Introduction Bloodcellproductionprocessisbaseduponthedifferentiationofso-calledhematopoietic stemcells,locatedinthebonemarrow.Theseundifferentiatedandunobservablecellshave uniquecapacitiesofdifferentiation(theabilitytoproducecellscommittedtooneofthethree bloodcelltypes:redbloodcells,whitecellsorplatelets)andself-renewal(theabilityto producecellswiththesameproperties). Mathematicalmodellingofhematopoieticstemcellsdynamicshasbeenintroducedat theendoftheseventiesbyMackey[21].Heproposedasystemoftwodifferentialequations withdelaywherethetimedelaydescribesthecellcycleduration.Inthismodel,hematopoi- eticstemcellsareseparatedinproliferatingandnonproliferatingcells,theselatterbeing ∗ToappearinMathematicalBiosciencesandEngineering 1

F.CrausteStabilityofabloodcellproductionmodel introducedintheproliferatingphasewithanonlinearratedependingonlyuponthenon- proliferatingcellpopulation.Theresultingsystemofdelaydifferentialequationsisthen uncoupled,withthenonproliferatingcellsequationcontainingthewholeinformationabout thedynamicsofthehematopoieticstemcellpopulation.Thestabilityanalysisofthemodel in[21]highlightedtheexistenceofperiodicsolutions,throughaHopfbifurcation,describing insomecasesdiseasesaffectingbloodcells,characterizedbyperiodicoscillations[19]. ThemodelofMackey[21]hasbeenstudiedbymanyauthors,mainlysincethebeginning ofthenineties.MackeyandRey[23,24,25]numericallystudiedthebehaviorofastruc- turedmodelbasedonthemodelin[21],stressingtheexistenceofstrangebehaviorsofthe cellpopulations(likeoscillations,orchaos).MackeyandRudnicky[26,27]developedthe descriptionofbloodcelldynamicsthroughanage-maturitystructuredmodel,stressingthe influenceofhematopoieticstemcellsonbloodproduction.Theirmodelhasbeenfurther developedbyDysonetal.[13,14,15],AdimyandPujo-Menjouet[7],AdimyandCrauste [2,3]andAdimyetal.[4].Recently,Adimyetal.[5,6]studiedthemodelproposedin [21]takingintoaccountthatcellsincycledivideaccordingtoadensityfunction(usually gammadistributionsplayanimportantroleincellcyclesdurations),contrarytowhathas beenassumedintheabove-citedworks,wherethedivisionhasalwaysbeenassumedtooccur atthesametime. Morerecently,Pujo-MenjouetandMackey[30]andPujo-Menjouetetal.[29]gavea betterinsightintothemodelofMackey[21],highlightingtheroleofeachparameterofthe modelontheappearanceofoscillationsand,moreparticularly,ofperiodicsolutions,when themodelisappliedtothestudyofchronicmyelogenousleukemia[16]. Contrarytotheassumptionusedinalloftheabove-citedworks,westudy,inthispaper, themodelintroducedbyMackey[21]consideringthattherateofintroductionintheprolifer- atingphase,whichcontainsthenonlinearityofthismodel,dependsuponthetotalpopulation ofhematopoieticstemcells,andnotonlyuponthenonproliferatingcellpopulation.Thein- troductionincellcycleispartlyknowntobeaconsequenceofanactivationofhematopoietic stemcellsduetomoleculesfixingonthem.Hence,theentirepopulationisincontactwith thesemoleculesanditisreasonabletothinkthatthetotalnumberofhematopoieticstem cellsplaysaroleintheintroductionofnonproliferatingcellsintheproliferatingphase. Thefirstconsequenceisthatthemodelisnotuncoupled,andthenonproliferatingcell populationequationdoesnotcontainthewholeinformationaboutthedynamicsofblood cellproduction,contrarytothemodelin[21,29,30].Therefore,weareleadtothestudyof amodifiedsystemofdelaydifferentialequations(system(3)–(4)),wherethedelaydescribes thecellcycleduration,withanonlinearpartdependingononeofthetwopopulations. Secondly,whilestudyingthelocalasymptoticstabilityofthesteadystatesofourmodel, wehavetodeterminerootsofacharacteristicequationtakingtheformofafirstdegree exponentialpolynomialwithdelay-dependentcoefficients.Forsuchequations,Berettaand Kuang[9]developedaveryusefulandpowerfultechnic,thatwewillapplytoourmodel. Ouraimistoshow,throughthestudyofthesteadystates’stability,thatourmodel, describedin(3)–(4),exhibitssimilarpropertiesthanthemodelin[21]andthatitcanbe usedtomodelbloodcellsproductiondynamicswithgoodresults,inparticularlywhenoneis interestedintheappearanceofperiodicsolutionsinbloodcelldynamicsmodels.Wewantto pointoutthattheusuallyacceptedassumptionabouttheintroductionratemaybelimitative andthatourmodelcandisplayinterestingdynamics,suchasstabilityswitches,thathave neverbeennotedbefore. Thepresentworkisorganizedasfollows.Inthenextsectionwepresentourmodel, statedinequations(3)and(4).Wethendeterminethesteadystatesofthismodel.In section3,welinearizethesystem(3)–(4)aboutasteadystateandwededucetheassociated characteristicequation.Insection4,weestablishnecessaryandsufficientconditionsforthe 2