Matlab code to run simulations of a batch reactor , MATLAB Programming

Background Protein therapeutics are a major component of the biotechnology industry, with sales estimated in the range of ~$99bn annually (2011) and steady market growth. Many pharmaceutical companies have been acquiring/developing intellectual property and the physical capacity to produce protein therapeutics, to offset revenues lost as traditional 'small molecule' drug revenues decline. As the development of protein therapeutics was fueled by advances in recombinant DNA technology that occurred between the late 1970's and late 1980's, the patents for some such drugs have now expired or are soon expiring, opening the market to generic versions, called "Bio-generics" or more often "Biosimilars."

In this project Your company is considering moving into the Biosimilars market and you need to design a system of reactors to produce a generic form of FancyMAb, commonly known as MAbZ1, which is secreted from cells. Populations of cells (CHO cells in this case) are grown from "clones," or individual cells (Nn=1), through division. However, individual clones of cells vary widely in their ability to produce (secrete) MAbZ1, and scaling up from the lab (batch reactors) to pilot plant (fed-batch reactors) often reveals differences in secretion potential between clones in each setting. Thus, you must test 8 different clones to find the one that performs best. The process of growing up single cells in the lab to large scale bioreactors is time consuming, so you need to plan out your experiments carefully and provide your boss with an approximate timeline for when you will need the large scale bioreactor for your studies and how long you will need it for. They will have to stop production of their high-value patent protected product while you conduct these experiments in the 5000 L reactor. You also have other work you would like to get back to, so you want to be done with this as soon as possible.

Under ideal conditions, the rate of cell growth at low density is given by:

{1 cell/ 0.1 ml < N < 5x105 cells/ml} dNcells / dt = k1 Ncells     k1 = 0.4 day-1

N = number of cells per ml

Cells under these low-density conditions are fed (i.e. the cells are collected and the growth substrate is replaced in the batch reactor) more often than needed to ensure that adequate nutrition is always present.

Once cells reach a critical density, their growth rate increases, and at such densities they consume substrate from the medium more rapidly, such that it typically begins to become depleted in a batch reactor. The substrate for cell growth is actually a mixture of many different nutrients, including glucose, proteins, amino acids, and dissolved oxygen, among others, but we will lump these together as "S" for simplicity. Too much substrate (> 3 M) is toxic to cells, so a culture is usually started off with 1 M substrate.

{5 x 105 cells/ml < N < 5 x 106 cells/ml}: dNcells / dt = {k6 Ncells [S]} / { k5 +[S]}

                                                                        d [S] / dt = - kNcells  

k5 =0.05 M

k6 = 1/day

k7 = 2 x 10-7 (M-ml)/(cells-day)

Finally, cells don't grow if they are too sparse or too dense:

{N < 1 cell / 0.1 ml; N > 5x106 cells/ml}   dNcells / dt =  0

The average amount of functional MAbZ1 secreted per cell depends on the cellular density, the time in the reactor (due to degradation), and the substrate concentration.

dE / dt = {k2 Ncells [S]} / { k5 +[S]} - k3 Ncells 2 - k4 Et       (Units of E: ng/ml)

k2=0.21 ng/cell-day for a 'typical' clone

k3=4 x 108 10-8 ng-ml/cell2/day


You must first estimate the optimal cell density for secretion of MAbZ1 using a small batch reactor, in which substrate can not easily be manipulated and will be consumed by the cells. You may wish to examine the effect of varying the starting cell density on the rate of MAbZ1 secretion to aid you in this estimate. However, you will only collect product from the final, largest reactor (5000 L). This final reactor is run as a fed-batch, starting at 20% capacity and optimal cell density, with water added continuously to keep the cell density optimal as the cells grow and substrate added as it is consumed to maintain the value of 3 M. Assume that MAbZ1 is only secreted in the presence of an inducing agent, which can be added to the medium but is not consumed. In the fed-batch reactor, the inducing agent will be added at t=0, and no MAbZ1 is present initially: E(t=0) = 0

You will need to:

Design a system of batch reactors to grow CHO cells from a single clone to a large population and produce maximal amounts of product in the 5000 L reactor (batch reactors available to for you to purchase range in volume from 25 μL - 5000 L, and you will need several).

Describe the number, type, and volume of each reactor needed to generate sufficient cells in a Table. Use a table to describe the timeline by which the 8 clones will be tested in each reactor. You will want to minimize the number of reactors in your system. For analysis of a single clone, create a multi-panel graphic of the system of reactors, and the cellular density as a function of time in each reactor (in one Figure).

Using MATLAB, run simulations of a batch reactor at different starting cellular densities to help determine/illustrate ideal cellular density needed for production of MAbZ1. Make a multi-panel figure of N, S, and dE/dt as a function of time and starting cell density and describe how you identified the optimal cell density you chose. 

In a separate figure, show how much S and total liquid volume (water and S) will need to be added to the fed batch reactor as a function of time to maintain the optimal cell density and the total amount of MAbZ1 secreted for a 'typical' clone.

Describe the results you obtained in text. The text should stand alone as a document without the figures, describing results concisely.

Describe your conclusions, the number and size of reactors used and the total amount of MAbZ1 produced per day in the 5000 L reactor.

Include your MATLAB code in an appendix after all display items.

Posted Date: 3/21/2013 3:09:20 AM | Location : United States

Related Discussions:- Matlab code to run simulations of a batch reactor , Assignment Help, Ask Question on Matlab code to run simulations of a batch reactor , Get Answer, Expert's Help, Matlab code to run simulations of a batch reactor Discussions

Write discussion on Matlab code to run simulations of a batch reactor
Your posts are moderated
Related Questions
1. Design a suitable phase-lead controller. You should explain in detail how your design has been developed. Your solution should include Bode plots constructed both manually and

Q. An expression for a sawtooth wave over the internal 0 ≤ t ≤ T0 is given by f(t) = At/T0. The student is encouraged to check the Fourier coefficients to be a 0 = A/2, an = 0 for

clear tic L=1; T=0.2; nust=2000; dt=T/nust; n=40; dx=L/n;   r=1;  omega=10:10:5000;%Store Range of Frequencies for Simulation u=zeros(n+1,nust+1);%

Customizing a Plot: Line Types, Color, Marker Types: Plots can be completed in the Command Window, if they are really simple. Though, at many times it is desirable to customiz

Design a controller which regulates flow and compensate pressure to my desire value in simscape.

Function Functions: The one reason for using function handles is to be able to pass functions to the other functions-these are known as function functions. For illustration

Format command: The format command can be used to control the spacing between the MATLAB command or expression and the outcome; it can be either compact or loose. >> format

Q.1: Consider the transmission of a sinusoid x(t) = cos(2f0t) through a channel a ected by multipath and Doppler. Let there be two paths, and assume the sinusoid is being sent fro

Relational Expressions: The Conditions in if statements use expressions which are theoretically, or logically, either true or false. These expressions are termed as relational

what equation equals 36