Andrew's Substrate Inhibition Model of Cell Growth

Biochemical Engineering


Problem Statement: The growth of an organism is limited by the availability of a particular nutrient, which at high concentrations also inhibits cell growth (Andrew's inhibition kinetics). We wish to cultivate this culture in a continuous bioreactor.

  1. Write the dynamic equations that describe the biomass and substrate concentrations in the continuous bioreactor.

    Solution:

    dx/dt = (μ-D)*x
    ds/dt = D*(sf-s) - μ*x/Y
    where μ = μm*s/(K+s+Ki*s2)
    

  2. The model from Part a) is: (circle applicable description)
      Structured/Unstructured
      Segregated/Nonsegregated
      Stochastic/Deterministic
      

    Solution: Unstructured, nonsegregated, & deterministic.

  3. Identify the operating parameters in your model. In other words, what major handles does an operator have control over?

    Solution: D and sf

  4. Find all the steady-states. How many steady-states exist for this system? What factor(s) determine which steady-state the bioreactor eventually approach?

    Solution:

    0 = dx/dt = (μ-D)*x
    0 = ds/dt = D*(sf-s) - μ*x/Y
    There are three distinct regions:
    1. For Dmax<D, the only steady-state is washout: x=0 & s=sf
         where Dmax=um/(2*sqrt(K+Ki)+1) is the maximum in the specific growth rate
         This occurs at smax=sqrt(K/Ki)
    2. For D<μsf<Dmax and sf<smax, there are two steady-states.
         Washout steady-state     x=0 & s=sf
         Growth steady-state      x=Y*(sf-s) & s=μ-1(D)  du(s)/ds>0
    3. Otherwise, there are three steady states:
         Washout steady-state     x=0 & s=sf
         Growth steady-state #1   x=Y*(sf-s) & s=μ-1(D)  du(s)/ds>0
         Growth steady-state #2   x=Y*(sf-s) & s=μ-1(D)  du(s)/ds<0
    
    See the following file for detail. The eventual steady-state is determined by the following three factors: D, sf, and the initial condition.

  5. Outline the major steps in localized stability analysis.

    Solution:

    Step 1. Derive dynamic equations for the system: dx/dt=f(t,x)
    Step 2. Find steady-state by setting dx/dt=0
    Step 3. Linearize the dynamic equation around the steady-state
    Step 4. Find the eigenvalues and eigenvectors
    Step 5. Determine localized stability based on the signs of the eigenvalues;
            determine trends based on the directions given by the eigenvectors
    

  6. Without actually performing the detailed mathematics, which steady-state(s) are stable? Give a brief reasoning for your answer. Which steady-state(s) are attainable in practice?

    Solution:

    Stable ones:
      For D > Dmax,  x=0 & s=sf
      For D<μsf,     x=Y*(sf-s) & s=μ-1(D)  dμ(s)/ds>0
    Others are unstable.
    

  7. How would you find the dilution rate that maximizes the rate of biomass production from your model? (If you wish, write the equation(s) that must be solved, but do not actually solve them.) From a purely operational consideration (not an economic one, which is discussed in the next part), would you operate at the dilution rate you have just calculated?

    Solution:

    Find the point at which D*x is maximum. Substrate inhibition is characterized by a point where the specific growth rate exhibits a maximum as the substate concentration increases. Let's call this point smax and μmax. When there is substrate inhibition, the point at which D*x is maximum coincides with the maximum in the specific growth rate. If sf<smax, this point is the "hill-top" of the specific growth rate versus s curve, i.e., μmax. If sf>smax, this point is μ(sf). We do not necessariy want to operate at this point, because it sits on the edge of a "cliff", too close to the washout condition. The biomass productivity declines very abruptly once we operate past that point.

  8. When you consider the cost of the substrate, do you expect the optimum dilution rate to shift to a higher value or lower value? Why?

    Solution: Shift to a lower value because the term -D*sf will pull the original biomass productivity curve D*x downward, more the higher the D. The same curve is pulled downward further the higher the cost of the substrate.

  9. Find the substrate concentration in the feed (sf) that maximizes the rate of biomass production. (Do not solve any equations. In view of the problem statement that high substrate level inhibits cell growth, should you limit sf to a low level in maximizing the rate of biomass production?)

    Solution:

    Operate at a high sf. In practice, the upper limit will be limited, among other factors, by the substrate's solubility. Several files in the class' Handout section deal with optimization with respect to sf. For example, see:

  10. Which of the following mode(s) of bioreactor operation would you definitely not consider for biomass production? Why?
      ---------------------------------------------------------
      Mode                       Consider?   Productivity?
      ---------------------------------------------------------
      Batch                         Y/N          H/L
      Fed-batch                     Y/N          H/L
      Continuous                    Y/N          --- (baseline)
      Cell recycle                  Y/N          H/L
      Reactors-in-series            Y/N          H/L
      Immobilized cell bioreactor   Y/N          H/L
      ---------------------------------------------------------
    

    Solution: Other than an immobilized cell bioreactor, which is totally unsuited for biomass production, all other modes of operation, each with its own advantages and disadvantages, can be considered. The solution is bold-faced below.

      ---------------------------------------------------------
      Mode                       Consider?   Productivity?
      ---------------------------------------------------------
      Batch                         Y/N          H/L
      Fed-batch                     Y/N          H/L
      Continuous                    Y/N          --- (baseline)
      Cell recycle                  Y/N          H/L
      Reactors-in-series            Y/N          H/L
      Immobilized cell bioreactor   Y/N          H/L (not even a candidate!)
      ---------------------------------------------------------
    

  11. Which of the above mode(s) of bioreactor operation do you expect to yield a higher biomass productivity compared to the continuous mode?

    Solution: Cell recycle.


Return to Prof. Nam Sun Wang's Home Page
Return to Biochemical Engineering (ENCH482)

Biochemical Engineering -- Andrew's Substrate Inhibition Model of Cell Growth
Forward comments to:
Nam Sun Wang
Department of Chemical & Biomolecular Engineering
University of Maryland
College Park, MD 20742-2111
301-405-1910 (voice)
301-314-9126 (FAX)
e-mail: nsw@umd.edu ©1996-2007 by Nam Sun Wang
UMCP logo