Mathematical models of dynamical systems facilitate the computation of characteristic properties that are not accessible experimentally. In cell biology, two main properties of interest are (1) the time-period a protein is accessible to other molecules in a certain state - its half-life - and (2) the time it spends when passing through a subsystem - its transit-time. We discuss two approaches to quantify the half-life, present the novel method of in silico labeling, and introduce the label half-life and label transit-time. The developed method has been motivated by laboratory tracer experiments. To investigate the kinetic properties and behavior of a substance of interest, we computationally label this species in order to track it throughout its life cycle. The corresponding mathematical model is extended by an additional set of reactions for the labeled species, avoiding any double-counting within closed circuits, correcting for the influences of upstream fluxes, and taking into account combinatorial multiplicity for complexes or reactions with several reactants or products. A profile likelihood approach is used to estimate confidence intervals on the label half-life and transit-time.
Application to the JAK-STAT signaling pathway in Epo-stimulated BaF3-EpoR cells enabled the calculation of the time-dependent label half-life and transit-time of STAT species. The results were robust against parameter uncertainties.
Our approach renders possible the estimation of species and label half-lives and transit-times. It is applicable to large non-linear systems and an implementation is provided within the PottersWheel modeling framework (http://www.potterswheel.de).