Team:USTC CHINA/Modeling/DesignsofImmuneExperiments

From 2013.igem.org

Introduction

Our mice experiment has primarily proven the validity of our project. However, just like most scientific immune experiments on animals, the aim of our mice experiment was verification instead of exploring the optimal conditions for the production of our T-vaccine. In fact, few optimization experiments have been done by pure scientific researches, as most scientists care about facts and theories only, whereas exploring the optimal conditions is often viewed as the task of pharmaceutical factories. Yet iGEM itself frequently involves industrial fields, which makes itself seems like more an engineered competition than a science competition sometimes.
We investigated the methodology of Design of Experiment (DOE) in our project, and realized although most papers propagate the wide application of DOE, the popularity of DOE is much lower in scientific fields compared with that in engineered fields. Perhaps the disparity of ideology between science and engineered cause this puzzling phenomenon. Consider the dual qualities of igem, we decided to explore on DOE further and design further experiments for pharmaceutical factories. We have also utilized DOE on our optimization of B.subtilis medium experiment. Various designs have been made, and the overall runs of them all are too large for our laboratory.
But general factories have adequate time and equipment to fulfill our designs, and different designs will meet the requirement of different situations.

Sweeping factors

The final effects of our T-vaccine hinge on various factors, in fact over ten factors. The more factors, the more runs. In our laboratory, any experiment involving over five factors is hard to design, whatever the methodology. Fortunately, we were designing experiments for pharmaceutical factories, which enabled us to take more factors into account without sacrificing accuracy too much. The first step of any method in DOE is to make a list of controllable factors, and then find out levels of each factors. In our design, we finally selected eight factor as follows:
The rate of four engineered bacteria, which produce antigen, LTB, KNFα and reporter respectively;(We selected the concentration of antigen as our standard, fixed at 1, and the rates of other three bacteria to engineered bacteria produced antigen provides three independent factors);

  • The area of the sticky patch;
  • The concentration of bacteria per unit area;
  • The body temperature of the vaccines;
  • The time consumed for culturing the bacteria;
  • The molecule weight of the antigen;

The ranges of these factor are given as follows:

Factor

Level Values

LTB

-1 0 1

KNFα

-1 0 1

Reporter

-4 -3

Temperature/℃

35.5 36 36.5 37 37.5

Time/h

4 5 6 7

Area/

1 3 7 10

Concentration(the number of engineered bacteria per square centimeter )

7 8 9

Molecule Weight(K D)

10 20 40 80


Note:

The ranges of rates and concentration of engineered bacteria were too large, and thus we used the common logarithms instead of the original values. For example, the low level of LTB was -1, meaning the lowest rate of LTB to antigen was 0.1.

Abstract of DOE methods

The classification standards of DOE methods are not unified, and according to one classification the DOE methods can be classified into three plots:
Factorial Design: Factorial Design is the most traditional method of DOE, and theoretically all other plots origin from it. Factorial Design is recommended when the ranges of factors is too large.
Response Surface Design: Response Surface Design utilizes response surface and excels in data analysis.
Taguchi Design: Taguchi Design utilizes orthogonal tables to decrease runs, and emphasizes the stability of qualities. Some mathematicians doubt the accuracy of this method, but its wide success has proven its power.
We have tried them all in our project.

Factorial Designs

To some extent, all DOE methods are branches of Factorial Design. The easiest subplot of Factorial Designs is Full Factorial Design, which means making a list of all combinations of all levels, which in fact tries nothing to minimize the runs. Surely the overall runs of Full Factorial Design is larger than any other method, but it does provide the most detailed information, so it is recommended when the factory does not care about money and time.
Generally Full Factorial Design has nothing mathematically sophisticated, all required is to list the specific values of all factors without any limitation on levels, which grants us more flexibility and freedom. Here is our table of levels of factors:

Factor

Level Values

LTB

-1 0 1

KNFα

-1 0 1

Reporter

-4 -3

Temperature/℃

35.5 36 36.5 37 37.5

Time/h

4 5 6 7

Area/

1 3 7 10

Concentration(the number of engineered bacteria per square centimeter )

7 8 9

Molecule Weight(K D)

10 20 40 80


And we got our first design, whose number of overall runs is 17820!
Full Factorial Designs 17280 runs
In reality we did not deem this level values table was detailed enough, but the number of runs was already enormous. Perhaps only the biggest pharmaceutical factory can afford this design.

Next we turned to traditional Factional Factorial Design. To minimize the runs, the levels of all factors were fixed at 2. A general 2-level-8-factor Full Factorial design contains 2^8=256 treatments, but we can further decrease the runs by defining alias. That is to say, define some specific factors as logical operation results of others.
Here we got a half Factional Factorial Design and a quater one, and the numbers of runs are 128 and 64, separately.
Factorial Designs 64runs
Factorial Designs 128 runs
Any effort trying to decrease runs will inevitably lower the cogency of the experiments, and this influence is irreversible whatever the methods. Factories are supposed to strike a balance between the accuracy of experiments and the costs they can afford when designing experiments.

Plackett-Burman Design

As an important subplot of Factorial Design, Plackett-Burman Design is excellent in dealing with mass factors. Generally it was applied in the primary experiments to select the key factors for further experiments. The number of runs can be controlled at very low values, although it is hard to get the best treatment from Plackett-Burman Design.
Naturally the levels of all factors were two. On most occasions it works as preparation for other DOE methods, like RSM. In our project, we made three Plackett-Burman Designs of 12 runs, 20 runs and 48 runs. The more runs, the more reliable results will be get, but even the last one still requires further designs.
Plackett-Burman 20 runs
Plackett-Burman 12 runs
Plackett-Burman 48 runs

Response Surface Design

Utilized response surface and gradient, Response Surface Design excels in data analysis, which makes it more mathematically gracefully than Taguchi Design, which partly accounts for why we selected it for our experiments on the optimization of medium. The most widespread subplots of Response Surface Design is Central Composite Design and Box-Behnken Design, both of which were considered when we designed our experiments on medium. The number of factors of Box-Behnken Designs is fixed on some given values, which does not include eight, therefore we had to turn to Central Composite Design (CCD). CCD itself contains three subplots, namely Central Composite Circumscribed Design (CCC), Central Composite Inscribed Design (CCI) and Central Composite Face-centered Design (CCF). CCC is mathematically preferred for it is rotatable, and we designed the experiments on CCC and CCF. The numbers of runs in half CCC and CCF designs were 154, whereas in quarter designs 90.
CCC 90runs
CCC 154runs
CCF 90runs
CCF 154runs
In spite of the mathematical advantages of CCC, the alpha value, which means the distance from axial point to the center point, is larger than one, and some absurd treatment might be yielded. In our half CCC design the alpha value was 3.364, while in quarter CCC design 2.828. In both designs, some treatments are irrational, because their area or concentration were negative, which obviously contradicts the common sense. However, factories can still adopt these CCC designs by giving up the irrational treatments.

Taguchi Design

Taguchi Design uses orthogonal tables to decrease the runs. Created by Doctor Taguchi, it has obtained wide success all over the world, especially in Asia. Different from Response Surface Design, it does not aim to calculate a fitting surface or function but just find out the best level value of each factor. Generally the number of runs is smaller compared with RSM, while the ranges of factors in Taguchi Design is relatively smaller.
We tried to make Taguchi Designs but our tool software is unable to make the design with eight factor. Additionally, our ranges of factors were too large for Taguchi Design, therefore we had to give up this method in our design.
If our T-vaccine is lucky enough to be produced at mass scale, we hope our designs could help these pharmaceutical factories.