 \tau_h$. As the number of hops becomes very large for $t \gg \tau_h$, the distribution becomes Gaussian again, due to central limit theorem. The kurtosis decays back to zero. Analytical results for the kurtosis behavior in this region will be reported elsewhere.

To conclude, we have proposed and analyzed a novel particle-tracking approach for identification of nonlinear interactions. Unlike common MSD techniques, our approach is based on the analysis of non-Gaussian characteristics of the particle dynamics, specifically the kurtosis of the displacement probability distribution. The functional dependence of the kurtosis on the measurement time carries a lot of information about the nonlinear interactions that contribute to the particle motion. E.g., if a particle is placed in a cage-type potential induced by cytoskeleton or transmembrane proteins, the resulting kurtosis of the displacement is a non-monotonous function with three distinct regions characterized by the change of the sign of the kurtosis slope. Specific structure of the kurtosis function depends on the characteristics of the potential: shape and size of the individual cells, heights and widths of the barriers but the change of sign feature is quite robust to the specifics of the potential. A number of other techniques have been proposed for identification of the potentials on biological membranes based on the analysis of individual trajectories. Most comprehensive approach is to solve the inverse problem of reconstruction of $U({\bf r})$ from the trajectories. However, such problem is generally ill-posed [@FokChouProcRoySoc2010] and needs very large statistics of trajectories. E.g., Ref.[@MassonVergassolaEtAlPRL2009] suggested to infer forces acting on the biomolecule [@AuthGovBiophysJ2009] requiring the multiple particle visits of each spatial location which is difficult to achieve experimentally. In addition, the potential $U$ in living cells can slowly change with time. This fact may limit the application of inverse problem approaches, that attempt to reconstruct the specific potential. Another approach [@KusumiSakoYamamotoBiophysJ1993; @KusumiAnnRevBio2005] focuses on identification of the potential barriers from MSD-based analysis. That approach is successful but requires very high temporal resolution of trajectories. One more method is based on the measurement of autocorrelation of SPT trajectories and recovering of the probability distribution of particle jumps [@YingHuertaSteinbergZunigaBullMathBiol2009] which indirectly displays the information about the inhomogeneity of the plasma membrane
\tau_h$. As the number of hops becomes very large for $ t \gg \tau_h$, the distribution become Gaussian again, due to cardinal limit theorem. The kurtosis decays back to zero. analytic solution for the kurtosis behavior in this region will be report elsewhere. 

 To reason, we have proposed and analyzed a novel particle - tracking overture for identification of nonlinear interactions. Unlike coarse MSD techniques, our approach is based on the analysis of non - Gaussian characteristics of the atom dynamics, specifically the kurtosis of the displacement probability distribution. The functional dependence of the kurtosis on the measurement fourth dimension carries a lot of data about the nonlinear interactions that put up to the particle movement. E.g., if a particle is placed in a cage - character potential induced by cytoskeleton or transmembrane proteins, the resulting kurtosis of the displacement is a non - monotonous function with three distinct regions characterized by the change of the sign of the kurtosis slope. Specific structure of the kurtosis function depends on the characteristics of the electric potential: human body and size of the individual cells, acme and width of the barriers but the change of sign feature is quite robust to the specifics of the potential. A number of early techniques have been proposed for identification of the potentials on biological membranes based on the analysis of individual trajectories. Most comprehensive approach is to solve the inverse problem of reconstruction of $ U({\bf r})$ from the trajectories. However, such trouble is generally ill - pose [ @FokChouProcRoySoc2010 ] and need very large statistics of trajectories. E.g., Ref.[@MassonVergassolaEtAlPRL2009 ] suggested to understand forces acting on the biomolecule [ @AuthGovBiophysJ2009 ] requiring the multiple particle visits of each spatial location which is difficult to achieve experimentally. In addition, the likely $ U$ in living cells can slowly change with time. This fact may restrict the application of inverse trouble approaches, that attempt to reconstruct the specific potential. Another overture [ @KusumiSakoYamamotoBiophysJ1993; @KusumiAnnRevBio2005 ] focuses on identification of the likely barriers from MSD - based analysis. That approach is successful but requires very high temporal resolution of trajectories. One more method is based on the measurement of autocorrelation of SPT trajectories and recovering of the probability distribution of atom jumps [ @YingHuertaSteinbergZunigaBullMathBiol2009 ] which indirectly displays the data about the inhomogeneity of the plasma membrane
 \tak_h$. As the number of hops becomes very lcege foc $t \gg \fau_h$, the distribution becomes Gaussien atain, eue to central limit tfeorem. Thv kurtosiw derays back to zero. Analytleal rsdultv for the kurtoxis behavimr in this reghov cill be reported elsewhere.

To concludq, we hafe proposed and wnaljzqd a novel particle-tracking approach fkr idenuification of nonlonear interactions. Unlike fommln MSD techniques, lur approacy is vased on the analysis of non-Gaussizn characteristics of the partizle dvnamics, spexidicwnly the kurvosis jf the displacement [robabikity distributlon. Tie fynctional dependence mf the kurtosis on the measgrzment time carries a oor of hnfosmatkin xboht tge nonpinxar interacfions that xontribute to the psrepvle motion. E.f., if a pwrticle is placed in a cage-type potentidl jnduced by cytoskeleton or transmembrane projeins, the wesulting kurtosis of the displacement is a non-mototonkjs ymkctiub aith three distinct regions characterized by ege cmange of the sigk of the kurtosis xllpr. Specific strocture of fhe kurtosis functlon depgnds ob the chawactrristics of the potential: syape and sizv of the individual cepls, heights and widtns of the barriers but che chznge of sigj feature ks quite robust go nhe vpecifics of the potential. A number of pther tdchnoques rave been oroposed for identificatlon oy the potentiald on biological membranes based on the analysos of individbal trsjectories. Moft comprehensire approcch is to solve nhe inverve problem jf reconstrucdlon of $U({\bf r})$ from thq trqjecrories. Fuwever, such prpblem is ytnerally iol-posed [@FokChouProgRoySud2010] and needs verv oarge statisticx ow twauertoriqv. E.g., Ref.[@MassmnVefgarxolaEgAlPRL2009] suggcstdd tp infer forces actinc on the biomolecule [@AitmGovBiophisJ2009] requiwing the multople particle visius of xach s'atial losation which is difficult to adhieve exoerlmentally. In wddibion, the potencial $U$ in living cells can slowly change with time. This faxt may limit the apkligation of interse [roblem a[proaches, that attempt to reconstruct the specific potentizl. Anodher wpproach [@KusumiSakoYamamotoBiophysJ1993; @KusumiAnnRevBio2005] focuses on identificarion oh ehe potentizl bsrriegs yrok MSD-bwsed enalysis. That appgoach is successful but requires very hich temporal resolution of trajrcgories. One mote method is based on the jeasuremrnt of autocorrelation of SPT ttajectories and recoverikg of the probzbilitj dmstribution of partivle jukps [@YungHuerraStcinbergZunigaBjllKabhBiol2009] whifh iidixectly displays the information avout tye inhomogeneity oh the plasma oembrzne
\tau_h$. As the number of hops becomes for \gg \tau_h$, distribution becomes Gaussian theorem. kurtosis decays back zero. Analytical results the kurtosis behavior in this region be reported elsewhere. To conclude, we have proposed and analyzed a novel particle-tracking for identification of nonlinear interactions. Unlike common MSD techniques, our approach is based the of characteristics the particle dynamics, specifically the kurtosis of the displacement probability distribution. The functional dependence of the on the measurement time carries a lot of about the nonlinear interactions contribute to the particle motion. if particle is in cage-type induced by cytoskeleton transmembrane proteins, the resulting kurtosis of the displacement is a non-monotonous function with three distinct regions characterized the change sign of kurtosis Specific of the kurtosis on the characteristics of the potential: of the individual cells, heights and widths of barriers but change of sign feature is quite to the specifics of the potential. A number other techniques have been proposed for identification of the potentials on biological membranes based on of individual trajectories. Most approach is to the problem reconstruction $U({\bf r})$ the trajectories. However, such problem is generally ill-posed [@FokChouProcRoySoc2010] and needs large statistics of trajectories. E.g., Ref.[@MassonVergassolaEtAlPRL2009] suggested to infer forces the [@AuthGovBiophysJ2009] requiring the particle visits of each location is difficult to achieve addition, potential cells slowly with time. This fact limit the application of inverse approaches, that attempt to approach [@KusumiSakoYamamotoBiophysJ1993; @KusumiAnnRevBio2005] focuses on identification of the barriers from MSD-based analysis. That approach is but requires very high temporal resolution of trajectories. One more method is on the autocorrelation of SPT trajectories and recovering of the distribution of particle jumps which indirectly displays the information about the inhomogeneity of plasma
 \tau_h$. As the number of hops becoMes very larGe for $T \gg \Tau_H$, tHe diStriBution becomes GAUssiAn again, due to central limIt theOrEM. The KUrTosis Decays bACk TO ZerO. ANaLytIcAL rEsultS foR the kurTosis behavIor In This region wiLL bE reported eLseWhere.

To conclUde, We have PrOpoSEd and AnaLyzed A novel PArticlE-tracking ApPRoach fOR identiFICaTion Of nonlinear interaCTiONs. Unlike common mSD tecHnIQuES, Our AppRoach is basEd On the ANalysis OF nON-gAusSIan characteriStics of the pARtiCle dynAmIcs, SPecifiCally ThE KurTosis of the dIsplAcement prObabilITy distrIBution. THe funcTioNal DepeNDeNcE of ThE KurTOsIs oN The MeasuremEnT tIme caRrieS A LOT of iNfoRmatIon abOut the nonlineAr iNterACtiOns thAt conTribUtE to thE partiCle moTiOn. E.g., if a particle Is plAced in a caGe-tYpE poTeNtial INduced By cYtoSkeletoN or tranSMemBrANE PrOteins, the resulting KuRTOsIs of the dIsplacEMeNt IS a non-monOtOnoUs fuNCTion wIth tHReE distincT regioNS cHaRacteriZeD by the ChAngE of The siGN of tHe kurtOsis slopE. SpecIFic structure of THe kurtosis funCTiON DePEnds On tHe characterIstiCS of tHe poTEnTiaL: Shape And siZe OF tHE individual cells, heiGhTs and wIdths Of the barriers But the chanGE OF sign feaTure IS qUIte robust to the SpeciFics of the pOTential. A NumbeR of other TechniqueS HAve been pRopOseD foR idENTiFication of the POTentIaLs on bioLogIcal memBraNes BasEd oN tHe analysiS of indivIdUaL tRaJecTorieS. most compReHenSiVe aPproaCH is to sOlve tHe inVeRsE ProBlem of rECoNSTrucTiOn Of $U({\bF r})$ fRoM the tRajeCTorIes. HoweVer, such prOblEM is gEnErAlly ill-Posed [@FokChouPRoCRoySoc2010] and NeEds Very laRGE statistIcs of trajectories. E.g., Ref.[@MASsonVerGasSolaETAlPrL2009] suggestEd tO infer ForCEs actiNg on thE biomOlEcuLE [@authGOVbiOphYsj2009] requiring THE muLtiplE pArtiCle visiTs of each spatial locATioN which is diffiCulT to aCHIeVe eXPeRImeNtALly. iN Addition, the poteNtial $U$ in liViNG cElls can sloWLy cHaNge with Time. ThiS fact MAy limit The applicAtion of inVeRse pROBleM approacheS, that attEmpt to recONstruCT tHe speCifIc poteNtIal. anothEr apprOAch [@kusumISakoYAmAmotoBIophySJ1993; @kusumiAnNRevBio2005] focuses on identifIcatioN of thE poTential baRriERs fRom MSD-basEd anAlysis. That AppRoaCh is sUccESsful But rEQuIreS Very hIgh tEMporal resOLuTioN OF tRajectories. oNE MorE methOd iS Based oN the Measurement of autoCOrrelation of SPt traJECtoRieS And rEcOvering of the prObaBiLITy distriBuTion of partiCle jumps [@yiNGHuerTaSteiNbergZUnigaBuLLmaTHBiol2009] wHich IndIrectly diSplAyS The infoRmAtIOn abouT the InHomogeNeity oF The pLASma membrane
 \tau_h$. As the number of  hops beco mes v ery  la rg e fo r $t  \gg \tau_h$,t he d istribution becomes Ga ussia na gain ,  d ue to  centra l  l i m itth eo rem .T he  kurt osi s decay s back tozer o.  Analyticalr es ults for t hekurtosis beh avi or inth isr egion  wi ll be  repor t ed els ewhere.
To conclu d e, we h a v eprop osed and analyzed an ovel particle- tracki ng ap p r oac h f or identif ic ation of nonl i ne a r int e ractions. Unl ike commonM SDtechni qu es, our ap proac hi s b ased on the  ana lysis ofnon-Ga u ssian c h aracter istics  of  th e pa r ti cl e d yn a mic s ,spe c ifi cally th eku rtosi s of t h e  dis pla ceme nt pr obability dis tri buti o n.The f uncti onal  d epend ence o f the  k urtosis on themeas urement t ime  c arr ie s a l o t of i nfo rma tion ab out the non li n e a rinteractions thatco n t ri bute tothe pa r ti cl e  motion.  E .g. , if a  part icle is  placedin a c a ge -t ype pot en tial i nd uce d b y cyt o skel eton o r transm embra n e proteins, th e  resulting ku r to s i so f th e d isplacement  isa  non -mon o to nou s  func tionwi t ht hree distinct regio ns  chara cteri zed by the ch ange of th e s ign of t he k u rt o sis slope. Spe cific  structure of the k urtos is funct ion depen d s  on thecha rac ter ist i c sof the potent i a l: s ha pe andsiz e of th e i ndi vid ual  c ells, hei ghts and  w id th softhe b a rriers b ut  th echa nge o f  signfeatu re i squ i terobustt ot h e sp ec if icsofth e pot enti a l.A numbe r of othe r t e chni qu es  have b een proposedfo r identifi ca tio n of t h e  potenti als on biological membr a nes bas edon th e an alysis of  in dividu alt raject ories.  Most  c omp r e hensi v e  a ppr oa ch is to s o l vethe i nv erse  proble m of reconstructio n  of  $U({\bf r})$  fr om t h e  t raj e ct o rie s. How e v er, such proble m is gener al l yill-posed[ @Fo kC houProc RoySoc2 010]a nd need s very la rge stati st icso f  tr ajectories . E.g.,Ref.[@Mas s onVer g as solaE tAl PRL200 9]  su ggest ed toi nfe r for ces ac ti ng onthe b io molecule  [@AuthGovBiophysJ2009]  requi ringthe  multiple  pa r tic le visits  ofeach spati alloc ation  wh i ch is  dif f ic ult to ac hiev e  experime n ta lly . In  addition,t h e  po tenti al$ U$ inlivi ng cells can slow l y change withtime . Thi s f a ct m ay  limit the app lic at i o n of inv er se problemapproach es ,  that  attem pt toreconst r u ct the sp ecif icpotential . A no t her app ro ac h  [@Kus umiS ak oYamam otoBio p hysJ 1 9 93; @KusumiAnnRe vBio2 0 0 5] fo c use s onid entific a tion  of the po tential bar riersfrom  MSD- based a na lysis.  Th at  approachi s success ful b ut requ ir es v ery  hightemp o r al re solu ti onof trajec t o ri e s.  O n e m oremetho dis b ased on t h e measur eme n t of au to cor r e lation of S PT traject ori es an d recovering of t h ep robab ilitydistri butiono f p ar ticle j ump s [@YingHue rtaSteinb e rgZ un igaB ullMathB io l200 9]  wh ic h indirectly di splays t he infor m a tion a bou t  the inhomo ge neity of the  plas ma membran e 
 \tau_h$._As the_number of hops becomes_very large_for_$t \gg_\tau_h$,_the distribution becomes_Gaussian again, due_to central limit theorem._The kurtosis decays_back_to zero. Analytical results for the kurtosis behavior in this region will be reported_elsewhere.

To_conclude, we_have_proposed_and analyzed a novel particle-tracking_approach for identification of nonlinear_interactions. Unlike_common MSD techniques, our approach is based on_the_analysis of non-Gaussian_characteristics of the particle dynamics, specifically the kurtosis of_the displacement probability distribution. The functional_dependence of the_kurtosis_on_the measurement time carries_a lot of information about the_nonlinear interactions that contribute to the_particle motion. E.g., if a particle is_placed in a cage-type potential induced_by cytoskeleton or transmembrane proteins,_the resulting_kurtosis of the displacement is_a non-monotonous function_with three_distinct regions characterized_by the change of the sign_of the kurtosis_slope. Specific structure of the kurtosis_function_depends on the_characteristics_of_the potential:_shape and size_of_the individual_cells,_heights and widths of the barriers_but_the change of sign feature is quite_robust to the specifics_of_the potential. A number_of other techniques have been_proposed for identification of the potentials_on biological_membranes based_on the analysis of individual trajectories. Most comprehensive approach is to_solve the inverse problem of reconstruction_of $U({\bf r})$ from_the trajectories._However,_such problem is_generally_ill-posed [@FokChouProcRoySoc2010]_and needs very large statistics of trajectories._E.g., Ref.[@MassonVergassolaEtAlPRL2009]_suggested to infer forces acting on_the biomolecule [@AuthGovBiophysJ2009] requiring_the_multiple particle visits of each spatial_location which is difficult to achieve_experimentally. In addition, the potential_$U$_in_living cells can slowly change_with time. This fact may limit_the application of_inverse problem approaches, that attempt to reconstruct_the_specific potential. Another approach [@KusumiSakoYamamotoBiophysJ1993; @KusumiAnnRevBio2005]_focuses_on identification of the potential barriers_from_MSD-based_analysis. That approach is successful_but requires very high temporal resolution_of trajectories. One more method is based on the_measurement of autocorrelation_of SPT trajectories and recovering_of_the_probability distribution of particle jumps [@YingHuertaSteinbergZunigaBullMathBiol2009] which indirectly displays the_information about_the inhomogeneity of_the plasma membrane