]. Essentially, compartments represent uniform spatial distributions or specific chemical compounds of the basic radioactive molecules; radioactivity concentrations in the various compartments are the natural state variables of the system; tracer flow, resulting from interchange of radioactive molecules between compartments, is modeled by a Cauchy problem for a system of linear ordinary differential equations (ODEs) for concentrations; the constant coefficients, also called rate constants or kinetic parameters, represent tracer kinetics and may be related to the action of enzymes, such as hexokinase (HK) responsible for phosphorylation in cells.

In typical compartmental problems the rate coefficients are unknown. The measured data are the total amount of tracer (concentration or activity) in a given region of interest, and the input function (IF), describing the time rate of tracer carried into the system. Tracer kinetics results from the solution of the inverse problem of determining the unknown rate coefficients compatible with the data, and the subsequent explicit determination of the concentration (or activity) of each compartment through the solution of the system of ODEs. In the applications presented in this paper, the inverse problem is solved in two steps: first, a formal expression of the solution of the direct Cauchy problem is evaluated, where dependence on the unknown rate constants is made explicit; second, an inversion algorithm is applied in order to recover the kinetic parameters through comparison of the formal solution with the data. The inversion algorithm makes use of an optimization regularization method which is based on a Newton-type algorithm.

“Classical” compartmental models have been developed under the assumption that phosphorylation and dephosphorylation of FDG occur in the same intracellular cytosolic volume, as described by [@Sokoloff] and [@Wernick]. Recent progresses in cell biochemistry have shown that the appropriate location of dephosphorylation is the endoplasmic reticulum (ER) [@Ghosh], which is spatially separated from cytosol. The schematic path of FDG kinetics is illustrated in  \[fig:FDG\_kinetics\].

Following this pattern for FDG kinetics in the cell, here we describe and formalize a new model consisting of three compartments which account for free FDG in cytosol, phosphorylated FDG in cytosol, and phosphorylated FDG in ER. The new biochemically-driven compartmental model is referred to as BCM; a classical simplified compartmental model (SCM) is recovered from the proposed model under the assumption that the ER is removed from consideration. In the forthcoming paper (Scussolini et al. manuscript in preparation), as a first test of its
]. Essentially, compartments represent uniform spatial distribution or specific chemical compound of the basic radioactive molecules; radioactivity concentrations in the respective compartments are the natural state of matter variables of the system; tracer stream, resulting from interchange of radioactive atom between compartments, is modeled by a Cauchy trouble for a system of linear ordinary differential equations (ODEs) for concentrations; the constant coefficient, also called rate constant or kinetic parameters, represent tracer dynamics and may be related to the natural process of enzymes, such as hexokinase (HK) creditworthy for phosphorylation in cells. 

 In typical compartmental problems the rate coefficient are unknown. The measured data are the total amount of tracer (concentration or activity) in a given region of interest, and the input function (IF), describing the time rate of tracer carried into the system. Tracer dynamics result from the solution of the inverse problem of determine the strange rate coefficients compatible with the data, and the subsequent denotative determination of the assiduity (or activity) of each compartment through the solution of the system of ODEs. In the application presented in this paper, the inverse problem is solved in two steps: first, a formal expression of the solution of the direct Cauchy problem is evaluated, where addiction on the unknown pace constants is made explicit; second, an inversion algorithm is applied in order to recuperate the kinetic parameters through comparison of the formal solution with the datum. The inversion algorithm make use of an optimization regularization method acting which is based on a Newton - type algorithm. 

 “ Classical ” compartmental models have been developed under the premise that phosphorylation and dephosphorylation of FDG occur in the like intracellular cytosolic volume, as described by [ @Sokoloff ] and [ @Wernick ]. Recent progresses in cellular telephone biochemistry have shown that the appropriate location of dephosphorylation is the endoplasmic reticulum (ER) [ @Ghosh ], which is spatially separated from cytosol. The conventional path of FDG kinetics is illustrated in   \[fig: FDG\_kinetics\ ]. 

 Following this pattern for FDG kinetics in the cell, here we describe and formalize a new model consisting of three compartment which account for free FDG in cytosol, phosphorylated FDG in cytosol, and phosphorylated FDG in ER. The modern biochemically - driven compartmental model is refer to as BCM; a classical simplified compartmental model (SCM) is recovered from the aim model under the assumption that the ER is removed from consideration. In the forthcoming paper (Scussolini et al. manuscript in preparation), as a inaugural trial of its
]. Esdentially, compartments rtpresent uniform spatial distrjbutions or specific chemical compouids if tht basic radioactive moleculed; radioaxtivmty concentratioia in thc varjlus eonpartments are the naturdl state variatlds of the system; tracer flow, resultind from onherchange of rwdiosstivs molecules between compartments, ia modeltd by a Cauchy pronlem for a system of lineag orfinary differentiap equations (ODEf) for concentfations; tht eonstant cogfficients, also called rate consgants or kinetix parwketers, reprxsent nracer kinetigx and kay be telated to the acvion of enzymes, such as hxxokinase (HK) responsyble for [hksphorylation in xeols.

In typhcal xomoarumeitam probpema the rate coefficienrs are unknown. The kewwured data ars the eoeal amount of tracer (concentration or abtivjty) in a given region od interest, and the inkut functijn (IF), describing the time rate of tracer carried hnto vhd svwtem. Geafer kinetics results from the solution of the jnfegse problem of debermining the unknpwj twte coefficievts compafible with the datw, and tre suvsequent txplivit determination of the cobcentration (je activity) of each compartmenc throogh thr solution of the systeo of ODEs. In thf applicafkons presented iv tmis paper, the inverse problem is solvev in cwo stepr: fitst, a fjrmal exprfssiok of the solution ov the dhrect Caucjy problem is evaluated, where dx'endence on tne unnnown race conxtants is madq explicit; secpnd, an nnverskon algorinhm is ap'lied in ordqr to recover jhe kinetic perameters thriugh comparkron of the forkal solutpou with thw data. The inversipn xmgorithm makes brt if an optimizayiov rqgllacizatymn method whhch ks cssed un a Newton-bypd alborithm.

“Classical” com[artjental models have bcen develiped undqr the assumpyion that phosphorjlatimn end delhofphorylation of FDG occur in tge same ijtrwcellular cyeosooic volume, ax described by [@Sokoloff] and [@Wernick]. Recxnt progresses in cell viochemistry have skoen that the apprjpriate lmcation of dephosphoeylation is the ekdoplasmic reticulum (ED) [@Ghosv], whifh is spatially separated from cytosol. The schematic path of FDG kinetics is illustrated in  \[fig:GDG\_kitecicx\].

Folloring vhis pattern for GDG kinetics in the cell, here wx describe aud formalize a new model consosging of three compartments which accohnt for gree FDG in cytosol, phosphorylayed FDG in cytosol, and pmosphorylqted RDG in ER. Vhe new biochemicallu-driveu compqrtmentql mldel is referfed tl as BCM; w clessncal simplified compartmental modwl (SCM) id recovered from the proposcd modem undet the assumption fhat hhe ET is removed ftom consideratoon. In the forthcomitg paper (Scussoluhi et al. manoscripj in preparatioi), as a first test of its
]. Essentially, compartments represent uniform spatial distributions chemical of the radioactive molecules; radioactivity are natural state variables the system; tracer resulting from interchange of radioactive molecules compartments, is modeled by a Cauchy problem for a system of linear ordinary equations (ODEs) for concentrations; the constant coefficients, also called rate constants or kinetic represent kinetics may related to the action of enzymes, such as hexokinase (HK) responsible for phosphorylation in cells. In compartmental problems the rate coefficients are unknown. The data are the total of tracer (concentration or activity) a region of and input (IF), describing the rate of tracer carried into the system. Tracer kinetics results from the solution of the inverse problem determining the coefficients compatible the and subsequent explicit determination concentration (or activity) of each compartment of the system of ODEs. In the applications in this the inverse problem is solved in steps: first, a formal expression of the solution the direct Cauchy problem is evaluated, where dependence on the unknown rate constants is made an inversion algorithm is in order to the parameters comparison the formal with the data. The inversion algorithm makes use of an optimization method which is based on a Newton-type algorithm. “Classical” compartmental been under the assumption phosphorylation and dephosphorylation of occur the same intracellular cytosolic described [@Sokoloff] progresses cell have shown that the location of dephosphorylation is the reticulum (ER) [@Ghosh], which The schematic path of FDG kinetics is illustrated \[fig:FDG\_kinetics\]. Following this pattern for FDG kinetics the cell, here we describe and formalize a new model consisting of compartments which free FDG in cytosol, phosphorylated FDG in cytosol, phosphorylated FDG in ER. new biochemically-driven compartmental model is referred to as BCM; classical compartmental model is recovered from proposed model under assumption that the removed from In forthcoming al. manuscript in preparation), as a test its
]. Essentially, compartments rePresent uniForm sPatIal DiStriButiOns or specific cHEmicAl compounds of the basic rAdioaCtIVe moLEcUles; rAdioactIViTY ConCeNtRatIoNS iN the vAriOus compArtments arE thE nAtural state vARiAbles of the SysTem; tracer floW, reSultinG fRom INtercHanGe of rAdioacTIve molEcules betWeEN compaRTments, iS MOdEled By a Cauchy problem fOR a SYstem of linear oRdinarY dIFfEREntIal Equations (OdES) for cONcentraTIoNS; THe cONstant coefficIents, also caLLed Rate coNsTanTS or kinEtic pArAMetErs, represenT traCer kinetiCs and mAY be relaTEd to the Action Of eNzyMes, sUCh As HexOkINasE (hK) ResPOnsIble for pHoSpHorylAtioN IN CElls.

in tYpicAl comPartmental proBleMs thE RatE coefFicieNts aRe UnknoWn. The mEasurEd Data are the total AmouNt of traceR (coNcEntRaTion oR ActiviTy) iN a gIven regIon of inTEreSt, AND ThE input function (IF), deScRIBiNg the timE rate oF TrAcER carried InTo tHe sySTEm. TraCer kINeTics resuLts froM ThE sOlution Of The invErSe pRobLem of DEterMining The unknoWn ratE Coefficients coMPatible with thE DaTA, AnD The sUbsEquent expliCit dETermInatIOn Of tHE concEntraTiON (oR Activity) of each compaRtMent thRough The solution of The system oF odes. In the aPpliCAtIOns presented in This pAper, the invERse problEm is sOlved in tWo steps: fiRST, a formal ExpResSioN of THE sOlution of the dIREct CAuChy probLem Is evaluAteD, whEre DepEnDence on thE unknown RaTe CoNsTanTs is mADe explicIt; SecOnD, an InverSIon algOrithM is aPpLiED in Order to REcOVEr thE kInEtic ParAmEters ThroUGh cOmparisOn of the foRmaL SoluTiOn With the Data. The inversIoN algorithm MaKes Use of aN OPtimizatIon regularization method WHich is bAseD on a NEwtoN-type algoRitHm.

“ClasSicAL” compaRtmentAl modElS haVE Been dEVElOpeD uNder the assUMPtiOn thaT pHospHorylatIon and dephosphorylATioN of FDG occur in The Same INTrAceLLuLAr cYtOSolIC Volume, as describEd by [@SokoloFf] ANd [@wernick]. RecENt pRoGresses In cell bIocheMIstry haVe shown thAt the apprOpRiatE LOcaTion of dephOsphorylAtion is thE EndopLAsMic reTicUlum (ER) [@ghOsh], Which Is spatIAllY sepaRated fRoM cytosOl. The ScHematic pAth of FDG kinetics is illusTrated In  \[fig:fDG\_Kinetics\].

FOllOWinG this pattErn fOr FDG kinetIcs In tHe celL, heRE we deScriBE aNd fORmaliZe a nEW model conSIsTinG OF tHree compartMENTs wHich aCcoUNt for fRee FdG in cytosol, phosphORylated FDG in cyTosoL, ANd pHosPHoryLaTed FDG in ER. The nEw bIoCHEmically-DrIven compartMental moDeL Is refErred tO as BCM; A classiCAL sIMplifiEd coMpaRtmental mOdeL (Scm) is recoVeReD From thE proPoSed modEl undeR The aSSUmption that the ER Is remOVEd froM ConSiderAtIon. In thE FortHcoming papEr (ScussolinI et al. mAnusCript In prepaRaTion), as A fiRsT test of its
]. Essentially, compartmen ts represe nt un ifo rmsp atia l di stributions or spec ific chemical compound s ofth e  bas i cradio activem ol e c ule s;  r adi oa c ti vitycon centrat ions in th e v ar ious compart m en ts are the  na tural statevar iables  o f t h e sys tem ; tra cer fl o w, res ulting fr om interc h ange of r ad ioac tive molecules be t we e n compartments , is m od e le d bya C auchy prob le m for a syste m  o f l ine a r ordinary di fferentiale qua tions(O DEs )  for c oncen tr a tio ns; the con stan t coeffic ients, also ca l led rat e cons tan tsor k i ne ti c p ar a met e rs , r e pre sent tra ce rkinet icsa n d mayberela ted t o the actionofenzy m es,  such  as h exok in ase ( HK) re spons ib le for phosphor ylat ion in ce lls .
In  t ypica l  compa rtm ent al prob lems th e  ra te c o ef ficients are unkno wn . Th e measur ed dat a  a re the tota lamo unto f  trac er ( c on centrati on ora ct iv ity) in  a  given  r egi onof in t eres t, and  the inp ut fu n ction (IF), de s cribing the t i me r at e  oftra cer carried  int o  the  sys t em . T r acerkinet ic s  r e sults from the solu ti on ofthe i nverse proble m of deter m i n ing theunkn o wn rate coefficie nts c ompatiblew ith thedata,  and the  subseque n t  explici t d ete rmi nat i o nof the concen t r atio n(or act ivi ty) ofeac h c omp art me nt throug h the so lu ti on  o f t he sy s tem of O DE s.In  th e app l icatio ns pr esen te di n t his pap e r, t he i nv er se p rob le m issolv e d i n two s teps: fir st, a fo rm al  expres sion of the s ol ution of t he  di rect C a u chy prob lem is evaluated, where depende nce  on t he u nknown ra teconsta nts is mad e expl icit;  s eco n d , ani n ve rsi on  algorithm i s a pplie din o rder to  recover the kinet i c p arameters thr oug h co m p ar iso n  o f  th ef orm a l  solution withthe data.Th e  i nversion a l gor it hm make s use o f ano ptimiza tion regu larizatio nmeth o d  wh ich is bas ed on aNewton-ty p e alg o ri thm.
“C lassic al ” c ompar tmenta l  mo delshave b ee n deve loped  u nder the  assumption that phosph orylat ion a nddephospho ryl a tio n of FDGoccu r in the s ame  in trace llu l ar cy toso l ic  vo l ume,as d e scribed b y  [ @So k o lo ff] and [@W e r n ick ]. Re cen t  progr esse s in cell biochem i stry have show n th a t  th e a p prop ri ate location o f d ep h o sphoryla ti on is the e ndoplasm ic retic ulum ( ER) [@ Ghosh], w hi c h is s pati all y separat edfr o m cytos ol .T he sch emat ic  pathof FDG kine t i cs is illustrate d in\ [fig: F DG\ _kine ti cs\].
F ollo wing thispattern for  FDG k inet ics i n the c el l, her e w edescribe a n d formali ze anew mod el  con sis ting o f th r e e com part me nts  which ac c o un t  f or fre e FD G incy toso l, phosph o rylatedFDG in cyto so l,a n d phos p ho r y lated FDGinER. T h e  new bioch e mica l ly - drive n comp artmen tal mod e l i sreferre d t o as BCM; a  classica l  si mp lifi ed compa rt ment al  mo de l (SCM) is reco vered fr om the p r o posedmod e l un d er the  a ssumption th at th eE R is re m o ved  from c o nsi der ati on.Inthef orthcoming paper (Scusso lini  et a l . m anusc riptin  prepara t ion) , as a f i rst test of i ts
]. Essentially,_compartments represent_uniform spatial distributions or_specific chemical_compounds_of the_basic_radioactive molecules; radioactivity_concentrations in the_various compartments are the_natural state variables_of_the system; tracer flow, resulting from interchange of radioactive molecules between compartments, is modeled_by_a Cauchy_problem_for_a system of linear ordinary_differential equations (ODEs) for concentrations;_the constant_coefficients, also called rate constants or kinetic parameters,_represent_tracer kinetics and_may be related to the action of enzymes, such_as hexokinase (HK) responsible for phosphorylation_in cells.

In typical_compartmental_problems_the rate coefficients are_unknown. The measured data are the_total amount of tracer (concentration or_activity) in a given region of interest,_and the input function (IF), describing_the time rate of tracer_carried into_the system. Tracer kinetics results_from the solution_of the_inverse problem of_determining the unknown rate coefficients compatible_with the data,_and the subsequent explicit determination of_the_concentration (or activity)_of_each_compartment through_the solution of_the_system of_ODEs._In the applications presented in this_paper,_the inverse problem is solved in two_steps: first, a formal_expression_of the solution of_the direct Cauchy problem is_evaluated, where dependence on the unknown_rate constants_is made_explicit; second, an inversion algorithm is applied in order to recover_the kinetic parameters through comparison of_the formal solution with_the data._The_inversion algorithm makes_use_of an_optimization regularization method which is based on_a Newton-type_algorithm.

“Classical” compartmental models have been developed_under the assumption that_phosphorylation_and dephosphorylation of FDG occur in_the same intracellular cytosolic volume, as_described by [@Sokoloff] and [@Wernick]._Recent_progresses_in cell biochemistry have shown_that the appropriate location of dephosphorylation_is the endoplasmic_reticulum (ER) [@Ghosh], which is spatially separated_from_cytosol. The schematic path of FDG_kinetics_is illustrated in  \[fig:FDG\_kinetics\].

Following this pattern_for_FDG_kinetics in the cell, here_we describe and formalize a new_model consisting of three compartments which account for free_FDG in cytosol,_phosphorylated FDG in cytosol, and_phosphorylated_FDG_in ER. The new biochemically-driven compartmental model is referred to_as BCM;_a classical simplified_compartmental model (SCM) is recovered from the proposed model under_the assumption that the ER is removed_from consideration. In the forthcoming paper (Scussolini et al. manuscript in_preparation), as a first test of its