Metabolic fluxes estimated from stable-isotope studies provide a key to understanding cell physiology and regulation of metabolism. A limitation of the classical method for metabolic flux analysis (MFA) is the requirement for isotopic steady state. To extend the scope of flux determination from stationary to nonstationary systems, we present a novel modeling strategy that combines key ideas from isotopomer spectral analysis (ISA) and stationary MFA. Isotopic transients of the precursor pool and the sampled products are described by two parameters, D and G parameters, respectively, which are incorporated into the flux model. The G value is the fraction of labeled product in the sample, and the D value is the fractional contribution of the feed for the production of labeled products. We illustrate the novel modeling strategy with a nonstationary system that closely resembles industrial production conditions, i.e. fed-batch fermentation of Escherichia coli that produces 1,3-propanediol (PDO). Metabolic fluxes and the D and G parameters were estimated by fitting labeling distributions of biomass amino acids measured by GC/MS to a model of E. coli metabolism. We obtained highly consistent fits from the data with 82 redundant measurements. Metabolic fluxes were estimated for 20 time points during course of the fermentation. As such we established, for the first time, detailed time profiles of in vivo fluxes. We found that intracellular fluxes changed significantly during the fed-batch. The intracellular flux associated with PDO pathway increased by 10%. Concurrently, we observed a decrease in the split ratio between glycolysis and pentose phosphate pathway from 70/30 to 50/50 as a function of time. The TCA cycle flux, on the other hand, remained constant throughout the fermentation. Furthermore, our flux results provided additional insight in support of the assumed genotype of the organism.