(Advances in Biochemical Engineering Biotechnology) R. Aarts, M. Aynsley, J.E. Bailey, P.M. Doran, I.J. Dunn, K. Friehs, E. Heinzle, A. Hofland
(Advances in Biochemical Engineering Biotechnology) R. Aarts, M. Aynsley, J.E. Bailey, P.M. Doran, I.J. Dunn, K. Friehs, E. Heinzle, A. Hofland
(Advances in Biochemical Engineering Biotechnology) R. Aarts, M. Aynsley, J.E. Bailey, P.M. Doran, I.J. Dunn, K. Friehs, E. Heinzle, A. Hofland
48 Biotechnology
With contributions by
R. Aarts, M. Aynsley, J. E. Bailey, R M. Doran,
I. J. Dunn, K. Friehs, E. Heinzle, A. Hofland,
K. B. Konstantinov, C. Di Massimo,
G. A. Montague, A. J. Morris, K. E Reardon,
G. B. Ryhiner, T. Yoshida
Springer-Verlag
Berlin Heidelberg NewYork
London Paris Tokyo
Hong Kong Barcelona Budapest
ISBN 3-540-56315-6 Springer-Verlag Berlin Heidelberg NewYork
ISBN 0-3 87-56315-6 Springer-Verlag NewYork Berlin Heidelberg
This work is subject to copyright. All rights are reserved, whether the whole or part of
the material is concerned, specifically the rights of translation, reprinting, re-use of
illustrations, recitation, broadcasting, reproduction on microfilms or in other ways, and
storage in data banks. Duplication of this publication or parts thereof is only permitted
under the provisions of the German Copyright Law of September 9, 1965, in its current
version, and a copyright fee must always be paid.
The use of registered names, trademarks, etc. in this publication does not imply, even
in the absence of a specific statement, that such names are exempt from the relevant
protective laws and regulations and therefore free for general use.
Editorial Board
The file has been compressed using the popular shareware program
“PKZIP” (Trademark of PKware INC., PKZIP is available from most
BBS and shareware distributors).
A u t h o r I n d e x Volumes 1 - 48 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 193
M . A y n s l e y , A. H o f l a n d , A.J. M o r r i s , G . A . M o n t a g u e a n d C. D i M a s s i m o
Department of Chemical and Process Engineering, University of Newcastle,
Newcastle upon Tyne, NE1 7RU, England
Abbreviations. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2 The Bioprocesses Studied . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.1 An Industrial Penicillin Process . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.2 An Industrial Mycelial Process . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
3 Real-Time Knowledge-Based Supervision of Bioprocesses . . . . . . . . . . . . . . . . . . . . . . . . . . 6
3.1 The Real-Time Knowledge-Based System Shell . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
3.2 Bioprocess Supervision Control and Analysis (Bio-SCAN) . . . . . . . . . . . . . . . . . . . . . 7
3.3 Knowledge-Based Systems for Bioprocess Scheduling Operations . . . . . . . . . . . . . . . . 8
3.4 Bioprocess Database System . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3.5 System Configuration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3.6 System Operation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
3.7 Integration with Process Analysis and Control Software . . . . . . . . . . . . . . . . . . . . . . . 14
4 Application of Real-Time Knowledge-Based Bioprocess Supervision . . . . . . . . . . . . . . . . . 14
5 Bioprocess State Estimation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
6 Artificial Neural Network Based Estimation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
7 Applications of Neural Networks to Bioprocess Supervision . . . . . . . . . . . . . . . . . . . . . . . 20
8 Discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
9 References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
The ability to supervise and control a highly non-linear and time variant bioprocess is of
considerable importance to the biotechnological industries which are continually striving to obtain
higher yields and improved uniformity of production. Two AI methodologies aimed at contributing
to the overall intelligent monitoring and control of bioprocess operations are discussed.
The development and application of a real-time knowledge-based system to provide supervisory
control of fed-batch bioprocesses is reviewed. The system performs sensor validation, fault detection
and diagnosis and incorporates relevant expertise arid experience drawn from both bioprocess
engineering and control engineering domains. A complementary approach, that of artificial neural
networks is also addressed. The development of neural network modelling tools for use in bioprocess
state estimation and inferential control are reviewed. An attractive characteristic of neural networks
is that with the appropriate topology any non-linear functional relationship can be modelled, hence
significantly reducing model-process mismatch. Results from industrial applications are presented.
Advancesin BiochemicalEngineering
Biotechnology,Vol.48
ManagingEditor:A. Fiechter
Springer-VerlagBerlinHeidelberg1993
2 M. Aynsley et al.
Abbreviations
AI artificial intelligence
CER carbon dioxide evolution rate
ICP intelligent communication protocol
KBS knowledge-based system
OUR oxygen uptake rate
RQ respiratory quotient
RTKBS real-time knowledge-based system
SQL Structural query language
TCS Turnbull control system
Artificial Intelligenceand the Supervision of Bioprocesses
1 Introduction
- better sensors
-better sampling and automatic analysis systems
-on-line estimation of difficult to measure, or unmeasurable, variables.
Although all of these are attracting research effort, the techniques of on-line
state and parameter estimation are receiving the greatest interest.
Whilst developments in biosensors are in some cases helping to contribute to
the on-line determination of bioprocess parameters, they are by no means at the
stage where general applicability has been achieved. Until this is the case the
development of on-line estimation techniques will provide the major means by
which on-line closed loop feedback control of critical bioprocess variables can
be achieved. (The need to obtain information from process variables which are
difficult to measure is not restricted to any one type of process. Problems can be
identified in such diverse industries as for example the processing of food,
chemicals, polymers, minerals, etc.)
Several different approaches can be adopted in the development of bio-
process estimation algorithms. The easiest, but potentially least accurate, way of
obtaining an estimate of a "difficult to measure" primary process variable is to
establish a correlation with a measurable secondary process variable, whilst
ignoring measurement errors, noise, etc. A more robust approach might be to
adopt a numerical estimation technique, either for estimating the parameters of
a pre-defined model structure, usually of a generalised linear time series form, or
directly to obtain an estimate of a difficult to measure primary process variable
(state estimation or soft-sensor). An alternative approach is to exploit the
nonlinear mechanistic structure of the bioprocess, if known, in the form of a
state observer.
A relatively new development in artificial intelligence is the area of neural
computing. Neural networks are dynamic systems composed of highly intercon-
nected layers of "simple" neurone-like processing elements. In the field of
process engineering it is their ability to capture process non-linearities which
offers the potential benefits. Other useful characteristics are their ability to
adjust dynamically to environmental and time-variant changes, to infer general
rules from specific examples, and to recognize invariances from complex high-
dimensional data. These properties provide neural networks with the potential
to out-perform other "learning" techniques. Given a series of examples the
Artificial Intelligence and the Supervision of Bioprocesses 5
network is able to establish the governing relationships in the training data. This
ability can be exploited to aid the nonlinear modelling, control and optimisation
of complex processes, as well as being used in existing predictive control
methods.
Two industrial systems are used to demonstrate the application of the AI
methodologies reviewed the fed-batch production of penicillin G and a large
scale mycelial process which operates in a continuous mode.
The diagnostic knowledge base has been structured into a number of discrete
rule-sets with meta rules again determining which rule-set should be assessed
next. The decision to invoke a particular set of rules is based upon a comparison
of the current raw, or processed, data with that permitted or expected in the
current phase of the bioreaction (characteristic trajectories of any number of
variables of interest can be incorporated into the system in addition to con-
straints). This ensures that knowledge can be focused in response to the
occurrence of specific events.
Once a fault has been detected the operator is informed of the problem on a
message board and additional rule-sets are invoked in an attempt to establish
the nature of the fault. A list of possible causes is then presented to the operator
together with advice on how to correct the fault. If Bio-SCAN cannot solve the
problem automatically, the operator is requested to confirm a diagnosis (e.g.
contamination). Such advice is readily achieved when the process is sufficiently
overdetermined, since the knowledge-based system capitalises on redundancies
to test for data consistency and error identification. The source of these
redundancies is usually balances or load cells, although software sensor estima-
tion of variables reduces the need for hardware redundancy [-14, 15].
Since structured, generic rule-sets have been used in the development of the
system, coupled with an object-orientated approach to the description of the
bioprocess, it is hoped that much of the knowledge base will be generally
applicable across a wide range of bioprocesses.
manipulated, to a certain extent, during the first few hours but then becomes
committed to transfer its culture at a fixed time. An appropriate production
vessel must therefore be available at this time. Production batches then typically
run for around 200 h before harvesting and recovery of product. In order to
achieve maximum plant efficiency, production vessels must be utilized in such a
way as to ensure high potency media are always available at regular intervals so
that downstream processing facilities are continuously supplied. In theory this
simply requires an orderly sequence of operations involving seed inoculations,
growth of seed, transfer of culture to production vessels and subsequent growth
and product formation. In practice, however, the ideal schedule may not be
realized because of "upsets" caused by factors such as contaminations of seed or
production vessels, maintenance requirements, plant problems, etc. The prob-
lem is compounded in that there are different sizes of both seed and production
vessels and there may be constraints imposed on where a seed vessel can transfer
its culture. This can result in lengthy downtimes for production units if an
existing schedule cannot be changed by redirection of previously committed
seed cultures.
Further disturbances to the "orderly" harvesting of batches arises as a
consequence of the inherent batch to batch variability in production rate. At any
one time several batches on the plant may have achieved a sufficient titre to
make them candidates for harvesting. However, it is often found that total plant
productivity can be improved by harvesting a vessel out of its inoculated
sequence. For this reason an on-line forecaster has been developed to identify
poor performing batches based on off-line titre information. One method of
forecasting is based upon the classification of historical titre profiles and
identification of the current batch type to predict future behaviour. An alter-
native strategy is to employ a curve-fitting approach to a number of possible
model structures. Whatever strategy is used, the forecaster seeks to identify
which bioreactors should be dispatched to downstream processing as fast as
possible and which should be left to run-on. The scheduler must then decide if
harvesting a bioreactor out of sequence is workable and economically feasible
with regard to the current schedule in operation.
The (simplified) problem description outlined above is not of the conven-
tional flow or job-shop type described in the literature, see e.g. Ku et al. [18] and
is essentially one of on-line rescheduling. A solution to the problem using a KBS
approach is being developed which combines a model of the plant with a module
which will use rule-based heuristics to examine operating strategies. This does
not, however preclude the use of algorithmic methods if a mixed approach is
found to be needed. The model is being de.signed to incorporate a detailed
knowledge of plant structure and operation and is based on the activity-event
representations used in Critical Path Analysis [19]. In this approach however,
both activities and events exist as connected objects. Knowledge of how each
object exists in the context of plant operation is contained in an associated frame
representation. The model will be used to test what-if scenarios by process
operators in a manner similar to the system proposed by Bernstein et al. [20].
10 M. Aynsleyet al.
off-gas
Air flow
measurement
] C02
analyser
~ 1~. ample
_J off-gas
I
I
i
. . . . . . . . . ]
i
i
TCS Serial
link Ethernet link
signal -.4F -- --I~-
processor
Serial link
I
T
Air Feed I
pump
I- i Balance
3 It is not strictly necessary to run an extension process on the same machine as the G2 real-time
expert system. G2 and GSI use the Intelligent Communications Protocol (ICP) which allows
extension processes to be run on different host computers and to communicate via ethernet.
4 RPC is a trademark of Sun Micro Systems.
Artificial Intelligence and the Supervision of Bioprocesses 13
State Estimation] ~ , , . I
,( S o f .t S e n s i n. g ,~ t~ - - [ ~S i m u. l a t i o n n~ IOptimisationl
P r e d i c t i o n of' | ~
Future Control [ ~
#
Bio - SCAN
Monitoring and Supervisory Control
- Phase Identification
- Operating Strategy
- I n t e r n a l Models
- Expected Trajectories
-
-
- Fault Detection & Diagnosis
Sensor Validation
R e c o v e r y Advice ~ l n t e llige n t /
User
Interface
~,9 ..................
[
Local C o n t r o l C o m p u t e r ] Off-Line
Contaminated Batch
,=,,
Feed Pump Failure
Z
8
(3
at
i i 410 i i i i f i i t i i i i I t l
TIME (hi
will fire, how these limits could change during the course of a reaction and
perhaps the probability of a process fault. The CO 2 must also be reconciled with
other process measurements such as the oxygen uptake rate and dissolved
oxygen concentration before diagnostic rule-sets are subsequently invoked. In
the contaminated bioreaction shown in Fig. 3, Bio-SCAN detected the in-
creased CO2 production rate early in the batch. Since this reading was found to
be reconcilable with other measurements and no other abnormalities in the
operation of the process could be detected, the system concluded that a
contamination had most likely occurred and informed the operator after 17 h. It
is interesting to note that this contamination was actually detected by the actual
process operators after 95 h when laboratory analysis identified the presence of
a yeast and terminated the reaction. Had the reaction been terminated earlier a
considerable saving in substrate feed costs and plant operation costs would have
been realized.
In the case of the second abnormal reaction shown in Fig. 3, Bio-SCAN
detected an atypical fall in the CO2 measurement after 150 h which was not
reconcilable with any mini-harvests. As all other operational parameters ap-
peared normal, the system invoked a rule-set which established from balance
measurements that substrate flow into the vessel was below normal. The
operator was advised to investigate the cause of this problem (e.g. a blocked line
or pump failure). Once the problem has been corrected Bio-SCAN is "informed"
via the intelligent user interface. Supervisory rules then modify the feeding
strategy in an attempt to bring microbial growth back onto a pre-defined
trajectory and recover the productivity of the reaction.
16 M. Aynsleyet al.
Analyser Recalibration
F-
z
8
o
[ I I I I I ~ I I I I I I I I I I
TIME (h)
Figure 4 compares the original two normal bioreactions with a batch during
which the CO2 analyser was detected to have drifted off calibration. Once the
suspected calibration error had been confirmed it was corrected at time 81 h and
normal bioprocess operation resumed.
bioprocess operator. The philosophy being that more frequent and more up-to-
date information enables improved process operation.
Numerous techniques for the estimation of important "hidden" bioprocess
variables papers have appeared in the literature since the early 1980s. The
methodologies developed have, however, yet to be proved completely reliable
and find acceptance with the industrial community. Industrial systems are
inherently non-linear and it is expected that the use of a process model which
reflects this essential non-linear structure would prove to be most beneficial. As
a consequence the estimation algorithms developed should either be inherently
based upon the non-linear structure or alternatively approximate the non-
linearities using a linear model with adaptation. Standard non-linear estimation
techniques are available, such as the extended Kalman filter [24]. They can,
however, suffer from numerical problems and convergence difficulties. Despite
these complications many successful applications have been reported, see e.g.
Svrcek et al. [25]; Stephanopoulos and San [26, 27].
Several examples were presented at the 1st IFAC Workshop on Modelling
and Control of Biotechnical Processes [28] in 1982, e.g. Dekkers [29];
Swiniarski et al. [30], and the following IFAC Symposium held in the Nether-
lands in 1985 [31], e.g. Ghoul et al. [32]; Montague et al. [33]. A clear message
from these early studies was that although model based estimation methods
could achieve quite reasonable performance, varying growth regimes and ill-
known process disturbance and noise characteristics could invalidate the pro-
cess description. Much of this early work concentrated on an extended Kalman
filter approach coupled with a mechanistic model, the accuracy of which played
a major determining factor in the quality of estimation. Deficiencies in the
process model, such as unmodelled process dynamics and model parameter
variations, have to a certain extent been overcome by applying adaptive
estimation techniques, e.g. Leigh and Ng [34]; Shioya et al. [35]; Chattaway and
Stephanopoulos [36]. In some cases estimation of the model parameters and
process state variables simultaneously is often necessary in order to obtain
reliable results.
A prerequisite for any model based estimation technique is the development
of a balanced mathematical model of the system under study. The model must
be balanced in the sense that a compromise must be made between the
conflicting requirements of dynamic complexity and the ability to obtain the
parameters of the model either through estimation or experimental measure-
ment. Several techniques have been suggested by Dochain [37] for the develop-
ment of bioprocess observers, as well as controllers [38]. The observers de-
veloped are either of fixed parameter form, partially adaptive form or fully
adaptive. The experimental validation of an observer for the on-line state
estimation of bioreactor performance (biomass and product) has recently been
presented by Dochain et al. [39].
Whilst a large number of process models are available, two major problems
are faced. Firstly, the parameters of the models are usually difficult to determine
experimentally. Secondly, and perhaps more importantly, the simplifications
18 M. Aynsleyet al.
made in any process model usually result in a need to adapt the model
parameters to compensate for dynamic model inaccuracies (model-process
mismatch). One approach that has been applied to the penicillin process in order
to provide an estimate of biomass concentration [40], is to estimate the
parameters of a CER model off-line and use the resulting model in conjunction
with an on-line estimator. This results in a partially adaptive observer. That is, a
priori knowledge of some of the model parameter values is required for observer
application, although it is possible to estimate some of the parameters on-line. In
such an observer, the error in the prediction of CER acts as the driving force for
state and parameter updating. Although convergence and stability proofs have
been developed for observers of the form described, providing certain key
criteria are met, partially adaptive observers tend to be insensitive to initial
conditions [37]. Adaptive observers that do not require any prior knowledge of
the model parameters are known as being fully adaptive. Here the parameters
are assumed to be unknown and are identified along with the state estimation.
An alternative starting point for the development of adaptive estimation
algorithms is to adopt the philosophy that the bioprocess dynamics can be
represented by a general linear input-output model. This approach is similar to
that adopted for the development of well known adaptive control laws, except
that an additional term representing the secondary process measurement is now
included. In the bioprocess plant the secondary measured variable, v, is CO2 in
bioreactor off-gas, the infrequently measured controlled variable, y, is biomass
concentration assay, and the manipulated variable variables, u, are the bio-
reactor feeds. The algorithm thus provides a means by which frequent estimates
Disturbances
Point
+
I ~ . J P,ooess
I y
i
j . . . . I
I Parameter i~1 _ 1
I Estimator I~1. . . . . .
19
I I Parameter
,~/I z estimates
Estimator ~
quality of the model fit produced by the neural network approach can be seen
from the plots in Fig. 6, where the step-like curve is the laboratory (off-line)
biomass assay. The bias that can be observed between the biomass assay results
and the neural network estimate, from around 25 h onwards, is due to a drift in
the on-line CO2 analyser calibration. This was corrected at around 170 h. At
around 210 h a major plant steriliser failure occurs which takes about 5 h to
correct. It is encouraging to observe that the neural network estimator remains
reliable over the period that this large process disturbance is effective. The
variable frequency of the laboratory assays is also quite evident.
The continuous bioprocess considered is required to operate as closely as
possible to pre-defined "steady state" operating conditions. In this respect a
linear model approximation was acceptable in the vicinity of normal process
operating conditions. A fed-batch bioprocess, however, presents more of a
problem in terms of model development; the system passes through a wide range
of dynamic situations never achieving a steady state.
The measurement problems experienced in controlling the penicillin process
are similar to those experienced in the continuous bioprocess. It is desirable to
control the biomass growth rate at a low level to optimise the production of
penicillin. Whilst the growth rate can be adjusted by varying the rate of
substrate (feed) addition, it is unfortunately not possible to measure the rate of
growth on-line. The aim therefore is again to develop a model which relates an
on-line bioprocess measurement to the process variable of interest. In this case
CER was used as the major variable from which to infer biomass concentration.
It is worthy of note, however, that in other studies the oxygen uptake rate
(OUR) has been used equally successfully. (In the penicillin process after the first
20 h or so, the ratio of CER to OUR, known as the respiratory quotient (RQ),
remains constant at around 0.9.)
z
_o
I,,.--
Z
I..IJ
O
Z
O
O
69
<
_o Estimate
-- Measured
4-
4-
,4- 94- 4--I-
Z
o 4-§ -4- 4- 4-
4-
I,Ll
0
Z
0
0
g
4--4-
O
N
Estimated
Measured
z
o_
t~ .4-
uJ ,~. "4"
o
z
0
o
4-
,<
0
Estimated
+ + + + Measured
200
T~ME(~,l
z
o
~
LI3
o
z
0
(,3
,<
+//~+ Estimated
+ + + + Measured
180
T~ME (h)
from the laboratory assays. The operators, having established confidence in the
network estimator, were led to check the calibration of the on-line CO z
analyser. The analyser calibration had drifted resulting in a significant error in
the predicted biomass. This was corrected at around the 100 h point. Towards
the end of the reaction (140 h) the estimated biomass concentration again
diverges. In this case off-line laboratory assays confirmed a contamination. This
type of behaviour, considered alongside other recent work in process engineer-
ing, suggests that with appropriate reconciliation of data there is potential for
24 M. A y n s l e y et al.
q"~ ,4.
<
1.1/
z
0
o
g
<
0
Estimated
+ + + + Measured
2012
TIME (h)
use in fault detection. Bioprocess faults can of course be associated with reactor
and sensor hardware as well as the bioreaction itself.
Acceptable biomass concentration estimates have been achieved without the
need for corrective action from off-line assays. It is possible, however, to make
such updates. In a similar vein, it is also possible to make use of a first-principles
model of the penicillin process in concert with a neural network. In such a
scenario the network is used to take account of the process/first-principles-
model mismatch. Studies are on-going in these areas.
8 Discussion
The initial objective of the RTKBS study outlined in this paper has been the
development of a generally applicable real-time knowledge-based system for the
supervisory control of bioprocesses. The supervisory strategy also performs
sensor validation, fault detection and diagnostic tasks and provides advice to aid
productivity recovery following fault correction. In order to achieve this aim,
generic-type rules have been used as much as possible, particularly in the case of
hardware fault detection and diagnosis. Process specific knowledge can be
incorporated within the existing framework relatively easily owing to the
modularity and object orientated approach adopted in the development of this
system.
Further work is underway aimed at incorporating on-line modules to
perform process simulation, soft-sensing and control. Intelligent use of this
information should significantly improve the overall performance of the system
Artificial Intelligence and the Supervision of Bioproccsses 25
9 References
8. Chantler MJ (1988) Real-time aspects of expert systems in process control. In: Expert systems in
process control, IEE Colloqium, March 1988
9. Stephanopoulos G, Stephanopoulos G (1986) Artificial intelligence in the development and
design of biochemical processes. Trends in Biotechnology, p 241
10. Moore RL, Hawkinson LB, Levin M, Hofman AG, Matthews BL, David MH (1987) Expert
systems methodology for real-time process control. Proceedings of the 10th IFAC World
Congress, Munich, FRG, July
11. Qi C, Wang S-Q, Wang J-C (1988) Application of expert system to the operation and control of
industrial antibiotic fermentation process. Preprints of 4th International Congress on computer
applications in fermentation technology, Cambridge, UK, September
12. Moore R (1988) The G2 real-time expert system for process control. International symposium
on advanced process supervision and real-time knowledge-based control, University of New-
castle, UK, November
13. Aynsley M, Peel D, Morris AJ (1989) A real-time knowledge-based system for fermentation
control. Proc ACC, Pittsburgh, USA, June
14. Guilandoust MT, Morris AJ, Tham MT (1987) Adaptive inferential control. Proc IEE, vol 134,
PtD
15. Tham MT, Montague GA, Morris AJ, Lant PA (1991) Soft-sensors for process estimation and
inferential control. J Proc Cont, vol 13-14
16. D.T.I Report (1989) Batch process scheduling computer software. The R&D Clearing House
17. Bruno G, Elia A, Laface P (1986) IEEE Computer 32:32
18. Ku H-M, Rajagopalan D, Karimi I (1987) Scheduling in batch processes. Chem Eng Prog,
August, p 35
19. Barnetson P (1969) Critical path plann!ng. Butterworth, London
20. Bernstein G, Carlson EC, Felder RM (1989) Development of a simulation-based decision
support system for a multipurpose pharmaceutical plant. AIChE Annual Meeting, San
Francisco, CA, November
21. Carleysmith SW (1988) Data handling for fermentation development - An industrial approach.
Preprints, 4th Int Cong on computer applications in fermentation technology, Cambridge, UK,
September
22. Flynn DS (1982) Instrumentation for fermentation processes. Proc 1st IFAC Conf, modelling
and control of biotechnical processes, Helsinki, Finland
23. Codd EF (1970) Comm ACM 13:377
24. Anderson BDO, Moore JL (1979) Optimal filtering. Prentice Hall, New Jersey
25. Svrcek WY, Elliot RF, Zajic JE (1974) Biotech Bioeng 16:827
26. Stephanopoulos G, San KY (1981) State estimation for computer control of biochemical
reactors. Advances in biotechnology, vol 1, p 399
27. Stephanopoulos G, San KY (1984) Biotech Bioeng 26:1176
28. Halme A (ed) (1982) Proceedings 1st IFAC workshop on modelling and control of biotechnical
processes, Helsinki, August, Pergamon, Oxford
29. Dekkers RM (1982) State estimation of a fed-batch bakers yeast fermentation. Proc 1st IFAC
workshop, modelling and control of biotechnical processes. Helsinki, Finland, August
30. Swiniarski R, Lesniewski A, Dewshi MAM, Ng MH, Leigh JR (1982) Progress towards
estimation of biomass in a batch fermentation process. Proc 1st IFAC workshop, modelling and
control of biotechnical processes, Helsinki, Finland
31. Johnson A (ed) (1985) Proceedings 1st IFAC symposium on modelling and control of biotechni-
cal processes, Noordwijkerhout, December, Pergamon, Oxford
32. Ghoul M, Pons MN, Engasser JM, Border J (1985) Extended Kalman filtering technique for the
control of Candida utilis production. 1st IFAC symposium modelling and control of biotechni-
cal processes, Noordwijkerhout, p 165
33. Montague GA, Morris AJ, Wright AR, Aynsley M, Ward A (1985) Parameter adaptive control
of the fed-batch penicillin fermentation. 1st IFAC symposium modelling and control of
biotechnological processes, Noordwijkerhout, p 39
34. Leigh JR, Ng MH (1984) Estimation of biomass and secondary product in batch fermentation.
6th International Conference on analysis and optimisation of systems, Nice, France, p t9
35. Shioya S, Shimizu H, Ogata M, Takamatsu T (1985) Measurement of state variables and
controlling biochemical reaction processes. 1st IFAC symp modelling and control of biotechni-
cal processes, Noordwijkerhout, p 49
36. Chattaway T, Stephanopoulos G (1987) A new technique for estimation in bioreactors: adaptive
Artificial Intelligence and the Supervision of Bioprocesses 27
James E. Bailey
Department of Chemical Engineering, California Institute of Technology,
Pasadena, CA 91125, USA
Introduction of a DNA vector into E. coli for the purposes of cloned gene expression can perturb
native cell functions at many levels. The presence of foreign D N A can alter regulation of
chromosomal D N A replication, regulation of transcription of chromosomal genes, ribosome
functions and RNA turnover, activities of regulatory proteins, chaperone and protease levels,
membrane energetics and protein post-translational processing, as well as energy and intermediary
metabolism of the cell. The combined effects of these interactions on the metabolic characteristics of
the host-vector system have major implications for yields of cloned biotechnological products and
interactions of genetically engineered organisms with their environment.
Advancesin BiochemicalEngineering
Biotechnology,Vol. 48
ManagingEditor: A. Fiechter
9 Springer-VerlagBerlinHeidelberg1993
30 J. E. Bailey
1 Introduction
The genome of an organism contains the genes and regulatory elements which
interact and function to establish the chemical, physical,.and functional charac-
teristics of that organism. "Wild-type" or "normal" attributes of the organism
are a consequence of these interactions and functions. In seeking to clone a
segment of DNA or to express a cloned gene, introduction of recombinant
vectors into the bacterium Escherichia coli is a common procedure. The
insertion of new DNA into the cell can perturb and in some instances severely
disrupt the normal balance of activities found in a wild-type cell. This review
considers experiments and mathematical modeling which have been conducted
to elucidate the types of perturbations which can arise.
Studies of the interlocking influences of vector and host cell upon each other
have major implications in biotechnology and basic science. For example, when
manufacturing a cloned therapeutic protein such as granulocyte colony stimu-
lating factor using genetically engineered E. coli, the reduction in specific growth
rate caused by high-level expression of the product must be taken into account
when designing the host-vector expression system and when operating the
process. From a basic science perspective, careful study of the metabolic and
regulatory characteristics of cells containing a recombinant DNA vector offers
insights into native cell function in the sense of a stimulus-response experiment.
Alterations in host cell metabolism which result from introduction of a cloning
vector are also of basic importance in the study of metabolic engineering. When
seeking to discern the influence on cell metabolism of expressing an active
cloned molecule in the cell, it is essential to recognize and to control for other
influences on cell functions which result from all other interactions with the
cloning vector.
This review will emphasize general points and will not endeavor to catalog
every experiment or study involving a different vector or host strain. Because all
data on interactions between hosts and vectors and their consequences are
obtained from in vivo studies, it is difficult if not impossible to dissect clearly the
most critical loci of the interactions. Many different types of perturbations and
changes in cell functions occur simultaneously. Nevertheless, in an effort to
emphasize the diverse ways by which introduction of new DNA can and does
alter cell function, the following comments will begin with the discussion of
perturbation of functions of host cell DNA, progress through perturbations in
RNA and protein functions and processing, and subsequently turn to other
special areas of metabolism before arriving at the most common subject of
vector effects on specific growth rates.
The key concept is the parasitic and perturbing influence of the added vector.
However, it is critically important that this parasitism be recognized as a
complex network of interconnections with many aspects of cell function and not
something as simple as "a demand for extra DNA synthesis" or "a demand for
more ATP". Considering alterations in the bacterial genome in a somewhat
32 J.E. Bailey
The most direct possible interaction of a cloning vector with the DNA of the
host and its functions occurs when introduction of the vector modifies the host
genome. This occurs for integrating vectors or if the autonomously replicating
vector contains a transposon or a region of homology with the host genome
such that some insertion is made into the host DNA. Such insertions can have
many effects, of course, ranging from silence to lethality by disrupting a function
essential for viability. Many other types of influences are possible by inactivating
important genes or regulatory elements. Such interactions are commonplace
with expression vectors for mammalian cells such that the resulting clones differ
widely in their growth characteristics, cell size, productivity, and other pheno-
typic characteristics, presumably due to different integration sites and their
subsequent consequences for host cell function. Such difficulties can be avoided
in E. coli, for example, by using vector and selection technology which localizes
the integration site such that a function unimportant for the growth conditions
of interest is targeted. Such a manipulation, for example, would involve integra-
ting a gene into the lacZ locus for a strain which would be grown on a glucose-
containing medium. Other types of interactions which can also exert profound
effects on the characteristics of the host-vector system are indirect in the sense
that the presence of the vector alters the native distribution of DNA-binding
proteins on the host DNA.
2.1 Replication
Bacterial plasmid vectors contain origins of replication consisting of regulatory
sequences which serve to control initiation of replication of the plasmid and also
the site of replication initiation. However, the vector typically does not encode
the necessary DNA polymerase and associated enzymes for DNA replication,
depending upon the host enzymes for these functions. Further, the plasmid may
depend upon other host proteins involved in regulation of chromosomal DNA
Host-Vector Interactions in Escherichia coil 33
ORI 0RI
-_ pr
Barn HI digestion
and ligation
Fig. 1. Geneology of the plasmids used in many studies of host-plasmid interactions. Indicated in
parentheses is the copy number (plasmids/cell) of each plasmid in E. coli HB101 grown in Luria-
Bertani (LB) medium. Alternative nomenclature for these plasmids is based on these reference copy
numbers (e.g., P24 denotes pDM246, P60 denotes RSF 1050, etc.). Reprinted from [24]
34 J.E. Bailey
Table 1. Estimated C and D periods for recombinant E. coil HB101 strains. Estimates were obtained
by two independent methods. Flow cytometry analysis of the single-cellDNA content distribution
following chloramphenicol (CM) addition combined with theory enabled estimation of the number
of active replication forks in an exponentiallygrowing population and thereby,based on population
theory, the C and D times. Direct comparison results were obtained by least squares optimization of
C and D in the theoretical steady-state single-cell DNA distribution compared to corresponding
flow cytometry data. Growth conditions: LB medium at 37 ~ at pH 7.0. Reprinted from [1]
CM-treatment Direct comparison
C D C+D C D C+D
Strain (min) (rain) (min) (min) (rain) (min)
HB101 only 36.8 19.9 56.7 37.3 20.5 57.8
HB101 :pDM247 30.3 15.4 45.7 31.0 15.5 46,5
HB101:pDM246 30.9 16.7 47.6 31.8 15.0 46.8
HB101 :RSFI050 33.3 16.2 49.5 33.8 13.7 47.6
HB101:pDM248 32.1 15.0 47.1 33.5 14.6 48,1
to expect increases in these cell cycle times rather than decreases. Furthermore,
it is extremely interesting that the plasmid copy number exerts almost no effect
on the magnitude of the change of C and D; merely the presence of plasmids, in
this case all plasmids with a ColEl-type origin of replication, is sufficient to give
rise to these responses in the host-vector construct.
2.2 Transcription
E. coli RNA polymerase apoenzyme and holoenzyme bind nonspecifically to
DNA. Accordingly, introduction of additional D N A into the cell as a multi-
copy plasmid can be expected to cause significant parasitic partitioning of RNA
polymerase onto the plasmid DNA. Of course by simple mass action relation-
ships, this parasitic partitioning would not be significant for a small, low copy
plasmid but could be substantial for a large or high copy vector. A model of
cloned gene expression and its regulation which considers the influences of
nonspecific binding on vector D N A on availability of polymerase for gene
transcription clearly indicates substantial effects when plasmid D N A content of
the cell approaches the D N A content of the host cell genome (corresponding to
around 1000 copies of a vector of 4 kb [2]).
Moreover, plasmids carry several promoters. For example, in pBR322 there
are two promoters associated with the origin of replication (for the small
transcripts RNAI and R N A I I ) as well as the promoters associated with the
ampicillin resistance and tetracycline resistance genes. In a construct designed
for high level expression of a cloned protein, promoters of such great activity are
employed that they may, even when present at relatively low levels, compete
strongly for available RNA polymerase holoenzyme.
These disturbances in the distribution of RNA polymerase relative to that in
the wild-type cell are expected to influence the level of expression of various host
Host-Vector Interactions in Escherichia coli 35
Table 2. Determinations of relative levels of total RNA and ribosome components in E. coli HB101
carrying different plasmids from Fig. 1 (P0 denotes the plasmid-freestrain). Reprinted from [3]
RNA"/protein 70Sb 50S~ 30Sc
Strain (lag lag- 1) (AUFS per 103 dpm) (70S- 1) (70S- 1)
P0 0.26 - -
P24 0.25 3.6 0.23 0.16
P60 0.28 3.7 0.19 0.14
P120 0.23 2.3 0.14 0.09
a Calculated with Ab26 o of 1.0 to be equivalent to 40 p.g RNA per ml.
b Normalized with dpm of 3H-labelled cells and Ab42 o of P24 AUFS denotes absorption units full
scale.
c Height ratio from absorbance (254 nm) scan of ribosome fractions separated in a sucrose gradient
cell genes. Experiments to examine ribosomal RNA contents and relative levels
of ribosome components in cells carrying different members of the family of
plasmids shown in Fig. 1 have clearly shown the consequences of such perturba-
tions [3]. As-plasmid copy number increases, RNA and active 70S ribosome
levels first decline slowly, then drop rapidly (Table 2). These trends parallel
those for specific growth rate discussed below.
W o o d and Peretti also observed a decline in ribosomal RNA (rRNA) levels
for increasing copy number in experiments employing the plasmids in Fig. 1 in
E. coli HB101 [4]. Interestingly, their detailed study showed that this decline in
steady-state rRNA levels with increasing plasmid content was a consequence of
an increase in rRNA degradation which more than compensated for an increase
in rRNA synthesis rate with increasing copy number. A different trend was
observed in studies of the consequences of induction of the strong tac promoter
in a pUC-based expression plasmid [5]. In this case rRNA content increased as
the cloned promoter was induced because the increase in rRNA synthesis rate
with increasing induction dominated a smaller increase in rRNA degradation
rate.
If the pattern of individual proteins found in the cell is regarded as indicative,
at least in part, of the corresponding pattern of transcription, the alterations in
cellular content of m a n y individual host cell proteins which arise upon introduc-
tion of multicopy plasmids indicate significant perturbations in host cell gene
expression patterns which presumably result from shifts in distribution of RNA
polymerase. Based on computer image analysis of two-dimensional gel electro-
phoresis patterns of recombinant E. coli HB 101 extracts, levels of many host cell
proteins are altered by presence of multicopy plasmids (Table 3) [3]. Levels of
stress response proteins detected are generally higher, while many metabolic
enzymes exhibit reduced levels at the highest copy number studied. Lower
activities of glucose-6-phosphate dehydrogenase, fructose 1,6-diphosphate al-
dolase, and fructose 1,6-diphosphotase were observed in E. coli D H 5 ~ carrying
the higher copy members of the plasmid family in Fig. 1 [6].
The appearance of foreign or misfolded native proteins in the cell has been
suggested to be the primary signal for the stress or "heat shock" response of the
36 J.E. Bailey
Table 3. Relative amounts of individual proteins in plasmid-free (P0) E. coli HB101 and in HB101
carrying P60 and PI20 as determined by analytical two-dimensional protein gel electrophoresis.
MW is molecular weight (kDa) and pl is isoelectric pH. Results for P0 indicate the fraction of that
protein relative to all proteins detected. P0 numbers denote 3sS spot intensity on the basis of a total
intensity of 106 for all spots detected on the P0 gel. The values for P60 and P120 strains were
normalized by the corresponding values for P0 to indicate changes in protein levels relative to the
plasmid-free host. A value of 1.0 thus indicates no change of protein level in comparison to the
plasmid-free strain P0. Reprinted from [3]
Strain
Mr
Name (kDa) pI P0 a P60 b P120 b
Carbon metabolism
Oxoglutarate dehydrogenase
complex
Oxoglutarate dehydrogenase 110.0 6.36 1314 2.68 0.92
Dihydrolipoate succinyl-
transferase 48.0 5.74 2239 2.55 1.14
Dihydrolipoate dehydrogenase 49.3 6.25 2716 1.72 0.81
PEP Carboxylase 92.9 5.88 215l 0.99 1.64
Pyruvate dehydrogenase complex
Pyruvate dehydrogenase 99.6 5.70 4652 1.13 0.60
Dihydrolipoate transacetylase 74.7 5.12 10,593 0.90 0.75
Dihydrolipoate dehydrogenase 49.3 6.25 2716 1.72 0.81
Pyruvate kinase I 53.8 6.23 8631 0.97 0.65
Succinate dehydrogenase 71.3 6.17 199 3.73 2.20
Succinate thiokinase (~) 41.8 5.66 1270 2.26 1.09
Amino acid metabolism
Aspartate aminotransferase 41.7 5.93 945 0.75 1.35
Glutamine synthetase,
Adenylated 50.7 5.52 2721 0.93 1.51
Unadenylated 50.8 5.57 1167 1.23 0.34
Nucleotide metabolism
Aspartate transcarbamoylase 33.9 7.00 1556 2.03 0.57
Carbamoyl-phosphate
synthetase (~) 136.2 5.57 1754 1.20 1.77
Protein synthesis
translation
Aspartyl-tRNA synthetase 66.1 5.75 1278 0.96 0.25
Elongation factor
G 87.5 5.60 11,475 1.35 0.67
Tu 43.2 5.71 103,503 0.71 0.63
Ribosomal subunit proteins
S1 69.1 5.00 8608 - 0.58
$5 17.2 - 21,326 1.15 0.45
S6a 15.0 5.70 2195 2.02 0.08
S6b 15.1 5.63 10,567 1.07 0.31
$6c 15.2 5.53 7674 0.60
L1 34.4 - 6358 1.14 0.45
L3 28.7 - 10,163 1.12 0.26
L7/12 12.2 4.91 28,168 0.92 0.84
Lll 14.7 - 26,411 0.85 0.28
L25 11.0 - 12,972 0.81 0.69
Heat shock proteins
Dna K 71.3 4.96 7506 1.14 1.12
Htp,
G 69.3 5.24 563 1.52 2.05
Host-Vector Interactions in Escherichia coli 37
Table 3. (continued)
Strain
Mr
Name (kDa) pI P0 a P60 b P120b
" This is ppm as dpm (specific protein) per dpm total valid protein spots.
b This is ppm per ppm P0.
c This is normalized with respect to P24 which has 2889 ppm
Fig. 2A,B. Two dimensional gel electrophoresis of E. coli BGF1 before (A) and after (B) the
induction of the CSH 11 mutant 13-galactosidase.Cells were grown at 37 ~ and pulse-labeled before
(A) and 30 rain after (B) IPTG addition. The following stress proteins were identified according to
known databases (Phillips et al., 1987): (1) Protease La (H94.0, H94.1); (2) HtpM (F84.1); (3) DnaK
(B66.0); (4) HtpG (C62.5); (5) GroEL (B56.5);(6) GrpE (B25.3);and (7) GroES (C15.4). The position of
the abnormal [3-galactosidase is indicated with an arrow. The right side of the gels corresponds to
lower pH values. Samples loaded each had 100000 epm total. Reprinted from [9]
copy of a lacZ gene containing the CSH 11 mutation [8]. Addition of IPTG to a
culture of BGF1 resulted in elevated expression of several stress response
proteins as determined by analytical two-dimensional gel electrophoresis
(Fig. 2) [9]. As a warning on the difficulties of designing experiments which
clearly delineate only host-plasmid interactions, the recent discovery that
addition of IPTG to wild type E. coli alters metabolism and protein pattern
should be noted [10].
One possible strategy for reducing alterations in host cell gene transcription
resulting from parasitic partitioning of host cell RNA polymerase is introduc-
tion of a separate polymerase-promoter pair for cloned gene expression. Indeed,
this strategy is embodied to some extent in the expression strategy which
employs phage T7 polymerase and the corresponding phage T7 promoter for
cloned gene expression [11]. However, this promoter is so strong that cell
specific growth rate is reduced following induction, presumably as a direct and
indirect consequence of the high level of message produced.
plasmid message in principle could reduce the deleterious effects of high levels of
plasmid message and of strong ribosome binding sites on those messages on
translation of host-cell message. An extensive investigation of this strategy and
its consequences has recently been presented by W o o d and Peretti [14]. They
observed a 30% reduction in specific growth rate of the culture following
expression of the specialized 16S RNA; however, no further decline in growth
rate was observed following induction of transcription of cloned [3-galactosidase
message bearing the modified SD sequence.
3.2 R N A Degradation
Table 4. Kinetic data for synthesis and degradation of !3-1actamasemessenger RNA in E. coil HB101
carrying the different plasmids shown in Fig. 1. The first column of data shows relative synthesis
rates. Division of these by the corresponding plasmid copy numbers gives the relative transcription
efficiencies rhx indicated. The relative [3-1actamasemessage specific degradation rate constants kdm!3L
were estimated from Eq. (1) making use of other measurements of relative steady-state ~3-1actamase
mRNA levels. Reprinted from [15]
Strain Counts per min a (10 -6) Yltx (10-8 cpm per gene kdm[3L (h 1)
per rain)
PI2 0.17 0.89 21
P60 0.72 0.75 2.8
P120 6.6 3.4 79
P408 690 110 330
aAmount of labeled [Mactamase mRNA on a per cell basis. Total RNA was pulse-labeled with [aH]
uracil, hybridized to nitrocellulose bound ~-lactamase DNA and digested with both RNase A and
RNase T 1 prior to scintillation counting
40 J.E. Bailey
multicopy plasmids or a small quantity of mutant protein into the cell activates
expression of stress response proteins including several chaperones. Besides
perturbing host cell protein metabolism, these interactions may also influence
the folding and targeting of the cloned protein.
Cloned proteins may be targeted for export through the cytoplasmic membrane
by gene fusion of a signal sequence with the product gene. Although this strategy
has not proven generally effective (only those proteins secreted in their native
hosts appear to be reasonable candidates for secretion in a heterologous host),
many such constructs have been employed effectively to obtain export of a
cloned protein into the periplasm of E. coll. In some cases, use of such fusion
constructs with strong promoters has resulted in lethality following induction.
This so called "overproduction lethality" has been attributed to interference of
the cloned preprotein with processing of host preproteins. That is, by mass
action, the large quantity of cloned preprotein occupies much of the limited
capacity of the host cell secretory apparatus, thereby impeding normal rates of
host cell preprotein processing and transport. Thus, there is again a parasitic
interaction and competition between the newly expressed cloned preprotein and
the normal processing of host cell secreted proteins. Another consequence of this
type of interaction may be reduction in integrity of the outer membrane ofE. coli
since its proteins are not processed sufficiently rapidly (relative to the cell
division rate) at the cytoplasmic membrane. This can result in leakage of
periplasmic proteins into the medium [see, for example, 19]. The latter phenom-
enon can be exploited to advantage in practice to produce cloned products
42 J.E. Bailey
Axe and Bailey have reported an in vivo nuclear magnetic resonance phospho-
rous-31 spectroscopy study of host and p!asmid-carrying E. coli [20]. Observa-
tions of intracellular and extracellular pH transients following addition of
glucose suggest a difference in membrane energetics between the transformed
and unmodified cells. Whether these changes result from differences in proton
pumping activities or differences in susceptibility to passive leakage of proteins
is not known. Also, it should be noted that the plasmid involved in these
experiments encoded pre-[3-1actamase, a protein new to the cell which must be
processed at the cytoplasmic membrane. Whether the observed effects on
transmembrane pH difference is connected with the expression of this cloned
preprotein is not established but worth considering.
Table5. Theoretical estimates of the maximum cell yield from ATP (Y.... ATP) and of ATP
requirements for synthesis of different classes of cellular components. Six cases are considered for
recombinant cells. The parameters n and 6 denote the number of plasmids per cell and the fraction of
total protein synthesized which is encoded on the plasmid, respectively. The "product retained"
figures include plasmid-encoded protein as part of total cellular protein; the "product secreted"
figures do not. Thus the latter figures may also be interpreted as requirements and yields for active
biomass. Reprinted from [23]
Recombinant cell
v e r y h i g h c o p y n u m b e r v e c t o r s d o e s n o t m a k e m u c h difference o n g l o b a l e n e r g y
m e t a b o l i s m . H o w e v e r , a c c o r d i n g to this s t o i c h i o m e t r i c analysis, p e r t u r b a t i o n s
in A T P r e q u i r e m e n t s o f the cell c a n b e c o m e i m p o r t a n t w h e n the e x p r e s s i o n level
o f the c l o n e d p r o t e i n s b e c o m e s significant r e l a t i v e to t h a t for all h o s t cell
p r o t e i n s . T h e fact t h a t s u b s t a n t i a l g r o w t h rate r e d u c t i o n s o c c u r w i t h h i g h Copy
p l a s m i d s a n d w i t h m o d e r a t e to l o w o v e r a l l e x p r e s s i o n o f c l o n e d p r o t e i n s [24] is
a clear i n d i c a t i o n t h a t the i n t e r a c t i o n s d i s c u s s e d in p r e v i o u s s e c t i o n s m a y
d o m i n a t e p l a s m i d effects o n c e l l u l a r A T P utilization,
T h i s t y p e o f a n a l y s i s of c o u r s e d o e s n o t i n c l u d e a n d c a n n o t e s t i m a t e the
effects of p l a s m i d p r e s e n c e o n m a i n t e n a n c e . As p o i n t e d o u t by S t o u t h a m e r a n d
others, m a i n t e n a n c e is a k i n e t i c r a t h e r t h a n a s t o i c h i o m e t r i c p r o c e s s and, in
E, coli, c o n s i s t s p r i m a r i l y of m a i n t a i n i n g the e n e r g e t i c c h a r g e of the c y t o p l a s m i c
m e m b r a n e [22]. U n f o r t u n a t e l y , there h a v e b e e n few definitive studies o n
44 J.E. Bailey
Table 6. Plasmid content effects on recombinant E. eoli HB101 maximum specific growth rates."
Reprinted from 1-24]
" Values shown are the ratio of measured specific growth rate to specific growth rate of the plasmid-
fiee host in the same medium. E. coli HB101 specific growth rates: LB medium: 1.38 h - l ; M9C
medium: 0.79 h 1; M9 medium: 0.41 h - i
Fig. 3A,B. Photomicrographs of E. coli JM101:pRED2 cells. A: phase contrast micrograph (bar is
1 pro). B: transmission electron micrograph (bar is 1 gm). In addition to refractile-type inclusion
bodies, cells also contain another inclusion body morphologydesignated floccule-type.Reprinted
from 1-31]
small amounts so greatly interfere with cell growth is also not known. Some
hypothetical causes include binding of the heterologous protein to crucial DNA
or RNA sequences, interaction of the heterologous DNA with chaperones or
proteases in the cell, or catalysis by the heterologous protein of reactions which
produce toxic components.
Host-Vector Interactions in Escherichia coli 47
In any description of cell kinetics, simple models which represent the cell by a
small number of pseudocomponents are convenient and often, by judicious
choice of a small number of parameters, can be adjusted to describe the growth
or product formation trajectories of interest. Although the components of such
small models are often assigned a loose physical interpretation such as total
nucleic acids or total protein or "total synthetic component", these labels are
more useful for conceptual guidance in proposing reasonable kinetic forms than
they are rigorously valid interpretations. A small structured model is basically
an empirical fit to the observed phenomena in which most if not all model
parameters have no clear physical relationship with the fundamental parameters
characterizing the reactions and equilibria actually occurring within the cell,
However, as a trade-off, the numbers of parameters involved and the com-
putational complexity of such models is minimized. Therefore, these types of
models are particularly useful for implementation in control schemes and for
representing experimental observations of the overall characteristics of the
cultures.
The four compartment model of Nielsen et al. includes variables for "active
biomass" interpreted to be primarily ribosomes, a lumped component consisting
of structural material and chromosomal DNA, plasmid DNA, and cloned
protein [33]. This model simulates many of the important qualitative features of
host-plasmid interactions including decreasing maximum specific growth rate
with increasing plasmid content and variations in RNA concentration and
plasmid content with culture specific growth rate. These trends are also well
simulated by an eight-compartment model formulated by Bentley and Kompala
[34].
Large-scale computer models containing many variables (say of order 30) and
parameters (in excess of 100) serve a different purpose. These models represent
an attempt to calculate, in a systematic and coordinated fashion, the conse-
quences of many simultaneous interactions within the cell, some of which are
represented in great detail, perhaps even at the level of the known molecular
mechanism. These models have an important place for hypothesis testing in the
following sense: the model in itself is an embodiment of a certain set of variables
and corresponding interactions precisely and rigorously defined by the model
equations. The model thus defines a self-consistent "computer cell" which
possesses its own intrinsic characteristics based entirely on the model's structure
and parameter values. Therefore, the model calculations with the computer cell
show the extent to which such a structure is capable of manifesting the behavior
observed experimentally for biological cells. As such, agreement between model
and experiment is an indication that the structure implemented in the model is a
plausible one for interpretation of the experimentally observed phenomenon.
Host-Vector Interactions in Escherichia coli 49
tRNA , ~
rRNA
mRNA c
TRANSCRIPTION
" RNA Polymerase
- Ribonucleotides
~ RNA p
mRNA p
Chromosomal ~ TRANSLATION
Plasmid
Protein " Ribosomes Protein
- AminoAcids
,o/ ,o/
PROMOTER PROMOTER
Fig. 5. RNA polymerase exists in many conformations. Binding of the sigma subunit to core
polymerase alters binding selectivityof enzyme.Resultant holoenzymehas much greater affinityfor
promoter sites and decreased affinity for nonpromoter regions relative to core polymerase. These
protein and DNA species and their relative affinitiesare explicitlycalculated in the Peretti-Bailey
single-cellmodel for E. coll. Reprinted from [35]
Table 7. Comparison of experimental and model simulation results for the specific growth rate and
the amount of cloned gene product per plasmid as a function of the relative number of plasmids per
cell. Reprinted from [12]
Cloned gene product, plasmid
Number of Relative specific growth rates specific synthesis
plasmids
per cell Experimental Calculated ExperimentaP Calculated
0 1 1
1 0.92 0.94 1 1
2 0.91 0.90 2.1 1.9
5 0.87 0.89 4.2 3.6
10 0182 0.87 5.4 4.1
20 0.66 - 5.2
34 0.77 - 7.0 -
Experimental values are taken from Ref. [24] for E. coli HB101 with the plasmid in Fig. 1 grown at
37 ~ in EB enriched with leucine, proline and thiamine and supplemented with 0.2% glucose
11 Concluding Remarks
References
Dedicated to Professor Dr. Karl Schiigerl on the occasion of his 65th birthday
1 P a r a m e t e r s Relating to D N A . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
1.1 P l a s m i d C o p y N u m b e r . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
1.1.1 Increased C o p y N u m b e r by P l a s m i d Design . . . . . . . . . . . . . . . . . . 55
1.1.2 Influence of C u l t i v a t i o n C o n d i t i o n s . . . . . . . . . . . . . . . . . . . . . . 56
1.2 P l a s m i d Stability . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
1.2.1 Segregational Genetic Stability . . . . . . . . . . . . . . . . . . . . . . . . 57
1.2.2 Structural Genetic Stability . . . . . . . . . . . . . . . . . . . . . . . . . . 59
2 P a r a m e t e r s R e l a t i n g to P r o t e i n Synthesis . . . . . . . . . . . . . . . . . . . . . . . . 60
2.1 P r o m o t e r s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
2.1.1 P r o m o t e r Strength . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
2.1.2 I n d u c t i o n . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
2.2 T e r m i n a t o r s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
2.3 m R N A . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
2.4 R i b o s o m a l Binding Sites . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
2.5 Stop C o d o n s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
2.6 C o d o n U s a g e . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
3 P a r a m e t e r s Relating to P r o t e i n s . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
3.1 Proteolysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
3.1.1 Protease Deficient Strains . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
3.1.2 I n h i b i t i o n of Proteases . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
3.1.3 F u s i o n Proteins . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
3.1.4 Protein E x p o r t . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
3.1.5 R a p i d P r o d u c t F o r m a t i o n . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
3.2 Inclusion Bodies . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
3,2.l A d v a n t a g e s and D i s a d v a n t a g e s of Inclusion Body F o r m a t i o n . . . . . . . . . . 65
3.2.2 Factors Influencing F o r m a t i o n . . . . . . . . . . . . . . . . . . . . . . . . 66
4 P a r a m e t e r s Relating to D o w n s t r e a m Processing . . . . . . . . . . . . . . . . . . . . . 67
4.1 Cell H a r v e s t and Cell D i s r u p t i o n . . . . . . . . . . . . . . . . . . . . . . . . . . 67
4.1.1 Cell Harvest . . : . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
4.1.2 Cell D i s r u p t i o n . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
4.2 Protein T r a n s p o r t . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
4.2.1 Signal Sequences . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
4.2.2 Secretion Systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
4.3 Protein F o l d i n g . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.4 Protein S e p a r a t i o n and Purification . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.5 I n a c t i v a t i o n of Biological W a s t e . . . . . . . . . . . . . . . . . . . . . . . . . . 73
5 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
6 References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
In the past 10 to 15 years, many of the promises of microbial genetic engineering have been realized:
the use of recombinant Escherichia coli has moved from the laboratory to the production facility,
and the manufacture of therapeutic recombinant proteins such as human growth hormone and
interleukins is a rapidly growing industry.
Along with this progress, however, have come new problems to solve: bioreactor operators have
discovered that large-scale cultivations of plasmid-containing bacteria do not behave in exactly the
same way as those of plasmid-free cells, plasmid stability has been recognized as a major hurdle, and
the protein product might not be present in a soluble form but rather as intracellular granules that
resist solubilization. These and other difficulties represent a new generation of challenges for genetic
engineering.
However, genetic engineering can do more than solve these problems. Molecular biological
techniques also have the ability to create new opportunities: to produce new compounds, to use
cheaper substrates, to facilitate downstream processing, and to optimize production in new ways.
The productivity of a cultivation can generally be expressed as the product of the cell density and
the specific biological activity. Both of these parameters are influenced by a variety of factors. For
recombinant cultivations, though, the level of biological activity, a reflection of the plasmid copy
number and expression efficiency, is the more interesting and important consideration and will
therefore be given more attention in our review. In this contribution, our general goal is to discuss
the factors that influence the productivity of recombinant E. coli cultivations, covering
--parameters relating to DNA;
- - parameters relating to protein synthesis;
- parameters
- relating to proteins; and
- parameters
- relating to downstream processing.
The object is not to tell the reader how to choose the perfect plasmid, host, and cultivation
conditions, but to make known the many variables involved in designing a recombinant process and
to point out recent and potential advances made possible by genetic engineering. The discussion
focuses on the production of a protein, but many of the same concepts apply to other cultivations of
recombinant E. coli, including cases in which the desired product is not a protein or the ceils have
been designed for a special metabolic capability such as pollutant biodegradation.
Parameters Influencingthe Productivityof Recombinant Cultivations 55
Since the copy number is closely related to plasmid replication, it has been
possible to design higher copy number plasmids by modifying replication
functions. Examples of this include:
- - c o p and tom mutations (in ColE1 derivatives)
- - Multiple tandem gene repeats
--Runaway replication plasmids
- - Dual low/high copy number origins of replication
The first two of these provide a high plasmid copy number throughout a
cultivation. The cop and rom mutations in ColEl-type plasmids influence the
regulation of plasmid replication [2, 3]. These types of mutations have been
used to produce ColE1 derivatives maintained at a level of 500 copies per cell
[4]. Similar opportunities for exploiting mutations in replication regulation are
available for plasmids such as R1, R6K, and pSC101 [5, 6, 7]. The use of tandem
gene copies in a plasmid has also been successful; for example, E. coli that
contained a plasmid with four copies of the chloramphenicol acetyltrans-
ferase gene produced four times as much product as cells harboring a plasmid
with a single gene copy [8]. However, tandem gene copies are prone to
homologous recombination.
Since continuous maintenance of plasmids at high copy numbers is a burden
to the host cell, resulting in significantly reduced growth rates as well as
increased plasmid instability, plasmid designs allowing one to increase the copy
number at an appropriate point in a cultivation are often desirable. One such
design is a temperature-sensitive "runaway" replication mutant, with which the
copy number of a plasmid is low at 25~ but rapidly increases when the
temperature is shifted to 37 ~ [9]. A different method involves the use of both
high and low copy number origins of replication on the same plasmid. Here, the
idea is to insert a controllable promoter in front of the replication primer of the
high copy number plasmid. An example of this is pMG411, which is maintained
at 4 copies per cell at 30~ and at 140 per cell when the culture temperature is
increased to 42 ~ [10].
56 K. Friehs and K. F. Reardon
Cultivation parameters also play a role in determining the copy number; thus, it
is also possible to change the copy number of a plasmid by manipulating the
cultivation conditions. The effects vary with the plasmid being used, but some
trends are evident.
Seo and Bailey [11] utilized different media to change the growth rate (and
perhaps other influential parameters) in batch cultures, and observed that the
copy number increased with decreasing growth rate. The influences of nutrient
limitation in continuous culture have also been studied. When a minimal
medium was used in a chemostat with glucose as the limiting substrate, Jones
et al. [12] noted significant decreases in copy number after one week of
operation.
When using runaway replication plasmids, the culture conditions have a
high impact on the achievable copy number. For example, the influence of
substrate feeding during runaway induction of pOU140 is shown in Fig. 1 [13].
Clearly, the addition of the carbon source, lactose, led to a significantly greater
number of plasmid copies than in the case in which no substrate was added.
100
Induction
O 80
~' 60
E
:~ 40
O
o 20
0
0 5 10 15 20 25
Time (h)
Fig. 1. Runaway replication with addition of lactose medium at point of induction. ( 4 1 - with
addition; --t-- without addition)
Although some wild type, low copy number plasmids exhibit high stability, most
plasmids of industrial interest (high copy number and expressing heterologous
proteins) are lost at frequencies of 10 -2 to 10-s per cell per generation [14]. A
recombinant host-vector system with an appreciable level of segregational
plasmid instability will be less productive than desired for two reasons:
- - t h e lowered average copy number will generally result in lower specific
productivity; and
- - plasmid-free cells which eventually emerge have a higher specific growth rate
(since they are no longer metabolically burdened by the plasmids) and will
thus reach a higher concentration than the plasmid-containing cells.
In a model developed by Imanaka and Aiba [15], cells were grouped as either
plasmid-containing or plasmid-free (i.e. copy number was not considered). Two
parameters were used to determine the fraction F of plasmid-containing cells
remaining after N generations: the frequency of plasmid loss (formation of
plasmid-free cells) and the ratio of the growth rates of plasmid-free to plasmid-
containing cells. The value of F was much more sensitive to the growth rate
ratio than to the loss frequency.
Strategies for overcoming segregational plasmid instability are based on the
alteration of one of these key parameters. Several types of approaches can be
identified. These have been classified as either selective (based on eliminating
plasmid-free cells) vs non-selective [16], or cellular/molecular vs bioprocess
methods [17]. In the following paragraphs, methods for preventing or overcom-
ing plasmid instability problems will be discussed according to the level at which
they are implemented, i.e.
- plasmid construction
-
- bioreactor configuration
-
maintained for more than 200 generations in the valS system [19]. However,
problems can arise when the mutated host cell is unable to grow rapidly or when
spontaneous reverse mutations occur.
Selection for plasmid-containing cells can also be accomplished by killing
plasmid-free cells. In one approach, the parB locus of plasmid R1 is incorporated
into the plasmid. This sequence includes two genes: hok, which produces a
potent bacteriocin, and sok, which forms an mRNA product that prevents the
translation of hok mRNA. Since the sok mRNA has a much shorter lifetime than
hok mRNA, ceils that lose the plasmid after division are killed by the bacteriocin
[21]. An alternative method of killing plasmid-free cells is to include the
bacteriophage ~ repressor on the plasmid and to infect the host cells with the
phage. Any plasmid-free cells that appear are then lysed [22].
The use of adjustable copy number (runaway replication) plasmids is a
different type of method that allows both good plasmid stability and high
productivity. Plasmid stability is maintained during an initial low copy number
growth phase due to a favorable ratio of growth rates between plasmid-
containing and plasmid-free cells. Following the growth phase, the copy number
is increased (e.g. by a temperature shift) and a period of high rates of product
formation occurs [9]. This approach can be applied to continuous cultivation
by utilizing a two-stage chemostat [23].
Most of the preceding methods are similar in that plasmid-containing cells
are given a growth advantage over plasmid-free cells without inherently increas-
ing the segregational plasmid stability. Another problem with selective tech-
niques is that they do nothing to maintain plasmid copy number but rather
require only that one plasmid per cell be present.
Non-selective, or genetic, approaches are more promising for ensuring that
high copy numbers are stably maintained. An example of these methods is the
inclusion of strong terminator sequences on the plasmid to prevent stability
problems posed by the strong promoters that are frequently used [24].
The par locus of pSC101 is responsible for stable maintenance of that
plasmid. When par has been cloned into other plasmids, their loss frequencies
have decreased [25, 26]. Meacock and Cohen [26] determined that the par locus
effectively increases the stability of low copy number plasmids but had little
impact on high copy number vectors. It should also be noted that insertion of
the par sequence can reduce the copy number [27].
Another example of a genetic method for enhanced plasmid stability is the
incorporation of the cer locus of the plasmid ColE1. When plasmids containing
cer are grown in a xer + (genes for site-specific recombination) host, multimeriz-
ation of plasmids is reduced, leading to increased plasmid stability [14].
Influences of plasmid copy number. The number of plasmids per cell has an
influence on the segregational plasmid stability of the culture. Although low
copy number plasmids with an active partition mechanism (e.g. par) are usually
very stable, high copy numbers generally result in greater stability when random
Parameters Influencing the Productivity of Recombinant Cultivations 59
2.1 Promoters
The synthesis of a protein starts with the promoter. The initiation of trans-
cription is a rate-limiting process for mRNA synthesis. Comparisons of more
than 100 promoters of E. coli have shown that there are two regions of
conserved DNA sequences [41, 42]. These regions, located 10 and 35 base pairs
upstream from the transcription initiation site, strongly influence the strength of
a promoter, which in turn determines the rate of transcription initiation. During
the development of the molecular biology of E. coli, many different promoters
were investigated; some of them are currently used in biotechnology and others
may be useful after further study.
2.1.2 Induction
2.2 Terminators
The first step of protein biosynthesis ( m R N A formation) ends with the termin-
ation of D N A transcription, when the R N A polymerase reaches the t e r m i n a t o r
sequence. W i t h o u t correct t e r m i n a t i o n , the m R N A would form the w r o n g
proteins, leading to an overall decrease in productivity. The use of strong
62 K. Friehs and K. F. Reardon
2.3 mRNA
The translation of mRNA, the second step in protein biosynthesis, starts with
the binding of the ribosomes at specific binding sites (RBS) on the mRNA
molecule. Weak binding sites lead to a low level of expression. In many systems,
the expression of foreign genes utilizes the native RBS. An increase in the
expression level can be achieved by replacing this natural RBS with an altered,
more efficient sequence [793. For example, Olsen et al. were able to enhance the
expression of bovine growth hormone in E. coli by enriching the sequence
flanking the RBS with A and T nucleotides [80]. The distance between the RBS
and the AUG start codon also has an impact on the rate of initiation of
translation [81]. A more complete understanding of the influence of the
secondary structure of the RBS on the rate of translation should provide further
opportunities to optimize foreign gene expression {-82].
Parameters Influencingthe Productivityof RecombinantCultivations 63
When using synthetic genes, the amino acid sequence of a protein is used to
develop a DNA sequence. Due to the nature of the genetic code, a choice of
codons is usually available. It has been shown that the use of particular codons
can play a role in gene expression. For example, highly expressed genes in
several species show a bias for certain synonymous codons [84]. Differences in
codon usage can also influence mRNA lifetimes.
3.1 Proteolysis
The levels of protease activity vary among E. coli strains, and this may represent
one criterion for host strain selection. More directed efforts have focused on
mutants deficient in the production of one or more of the known proteases. The
best known examples are the lon- mutants, which cannot form the cytoplasmic
64 K. Friehs and K. F. Reardon
Finally, producing the cloned gene product at high rates so as to saturate the
proteolytic enzymes is an effective method of minimizing degradation. This can
be accomplished with runaway replication vectors or with plasmids containing
inducible promoters. Since this strategy has other benefits (e.g. increased
plasmid stability), it is frequently implemented.
Parameters Influencing the Productivityof RecombinantCultivations 65
If the product is not secreted into the medium or at least exported into the
periplasmic space, the cells must be harvested and disrupted to release the
product. This is the case for both inclusion bodies and soluble intracellular
proteins.
Both mechanical and non-mechanical methods are used to disrupt cells in large-
scale processing [105, 107]. Mechanical methods such as high pressure homo-
genization and bead mills often result in product losses through inefficient
disruption or thermal denaturation of proteins. Non-mechanical methods in-
clude those based on physical processes (e.g. sonication) and chemical effects
(e.g. organic solvents and enzymatic lysis). Each of these have disadvantages; for
68 K. Friehs and K. F. Reardon
In many cases, it would be desirable to secrete the recombinant product into the
medium. Purification would be simpler than for an intracellular protein since
the product would not be contaminated with cytoplasmic components. (It
would be necessary to handle larger volumes, but this problem has been lessened
by newer chromatographic methods.) In addition, the formation of inclusion
bodies could be avoided and the toxic effects of some protein products on the
host cell could also be reduced. Protein secretion would also reduce proteolysis,
unless exoproteases are produced by the transport system.
Excretion into the periplasmic space of E. coli also provides many of these
benefits. Proteins can be freed from the periplasmic space with gentle treatments
that remove the cell wall and outer membrane.
A significant disadvantage facing protein secretion is that protein folding
may be incorrect or may not occur. As discussed in the following section, this is a
problem in the production of all recombinant proteins. The investigation of
special reactors for promoting folding is underway in several laboratories.
Two important factors influence the transport of proteins: the type of cellular
transport system and the nature of the signal (or leader) sequences that allow a
protein to use that transport system.
Parameters Influencingthe Productivityof RecombinantCultivations 69
Signal sequences are short peptides that allow the protein to be transported. In
bacteria, these signal sequences are typically 15 to 30 amino acids long [114,
1t5] and are located at the N-terminus of proteins, although some, like that of
hemolysin, are positioned at the C-terminus. Signal sequences form positively
charged heads that are able to pass through the membrane. Between the signal
peptide and the preprotein, a cleavage site for specific signal peptidases is
located to allow formation of the mature protein after transport.
In addition to the leader sequences, transported proteins often have other
sequences that are important for secretion. One type, located inside the protein
sequence, aids in the formation of a transport competent structure. The function
of other sequences, found in membrane proteins, is to stop transport. Thus, the
fusion of a gene to a signal sequence does not guarantee that the product will be
transported [116].
Some signal sequences that have been used to transport recombinant
proteins are shown in Table 2. Well known signal sequences like bla (~-
lactamase), malE (maltose binding protein), ompA (outer membrane protein),
and phoA (alkaline phosphatase) come from proteins that are exported into the
periplasm. Fusion proteins formed with one of these are therefore mainly
transported into the periplasm, but some secretion has also been observed. Of
special interest are the signal sequences spA (Staphylococcus protein A) and
malE, since both can be used as purification tags in affinity chromatography (see
Sect. 4.4).
Table 2. Signalsequences and target proteins for transport in E. coil; the final protein location is
shown as periplasmic space (P) or medium (M)
Other abbreviations: amy = amylase, bla = [Mactamase, BRP = bacteriocine release protein, hly
= hemolysin,M = medium, malE = maltosebinding protein, ompA = outer membrane protein A,
ompF = outer membrane protein F, P = periplasmic space, phoA = alk. phosphatase, phoS
= phosphate binding protein, spA = Staphylococcus aureus protein A.
The correct folding of a protein is essential for its biological activity. Although
issues of recombinant protein structure and overproduction receive a great
amount of attention, relatively little work has been done on protein folding,
despite its immense commercial impact and status as an important fundamental
question [ 152, 153]. The factors influencing protein folding are diverse and must
be considered when the production of active proteins is desired 1-154, 155, 156].
The production of recombinant proteins leads to at least two additional
problems with regard to protein folding: incorrect folding and the need to
renature proteins produced as inclusion bodies.
Some proteins, especially those from eukaryotic sources, are not folded
correctly in the oxidative milieu of the bacterial cytoplasm. One way to address
this problem is to export the protein into the periplasmic space of E. coli (via
fusion to a signal sequence). The periplasmic space is a reducing environment
and supports correct folding of eukaryotic recombinant proteins.
The need to refold proteins from inclusion bodies can be obviated by
preventing their formation. As discussed earlier (Sect. 3.2.2), this can be achieved
by using low-expression systems, by changing the cultivation conditions, or by
fusion to a signal peptide, among other methods.
A good technique for improving protein purification is the use of tags fused
to the protein product. The tag (e.g. Staphylococcusprotein A) is chosen to allow
the use of efficient separation techniques such as affinity chromatography. In
order to facilitate removal of the tag following purification, a specific cleavage
site can be placed between the tag and the desired protein. In Table 3, a number
of different tags and their chromatographic ligands are presented.
An example of the power of this method is presented in Fig. 2. In this case,
Staphylococcus protein A (SPA) was used as a fusion tag for the production of
101 15
8 4"--
E
t-
O o
00 6 3 ~
Cq o
"r"
c-
O v
4 2
L
0 . M
r
..Q
< 2
0
2 4 6 8 10 12 14 16 18 20 22 24 26 28 30
Fraction
Fig. 2. Affinity chromatography of EcoRI-SpA fusion on IgG 9 Total protein (Abs.) [] EcoRI
activity
Parameters Influencing the Productivity of Recombinant Cultivations 73
0.1
[////
r////
er i////
0.01
V/A
////
i/ill
I/// i/ill
J / / / / ~
0.001
Zephirol Formal- Formic acid Phenol Tego
dehyde
the restriction enzyme EcoR1 [1583. The fusion protein could be enriched and
purified in one step using affinity chromatography with IgG as the ligand.
For safety reasons and for public acceptance, it is necessary to inactivate the
biological waste from the production of recombinant proteins. While this is also
done for many other cultivations, inactivation of recombinant microorganisms
involves not only killing whole cells but also significantly reducing the level of
recombinant DNA (such as transformable plasmid DNA).
A series of experiments was performed to determine the efficiency of plasmid
inactivation by different methods. The intact biological activity of the plasmids
was measured by comparing the transformation rate of isolated plasmid DNA
before and after treating the cells. The results of these experiments (shown in
Fig. 3) demonstrate that several commonly used compounds, including formal-
dehyde, formic acid, phenol, Tego, and Zephirol cannot inactivate plasmids
when used in normal concentrations [1643.
The only truly efficient method to inactivate plasmid DNA with respect to
transformation is thermal incubation. This can be accomplished continuously,
just as sterilization is.
5 Conclusion
In this review, we have presented a wide variety of factors and techniques that
influence the productivity of recombinant E. coli cultivations. These parameters
74 K. Friehs and K. F. Reardon
6 References
38. Henry KL, Davis RH, Taylor AL (1990) Biotechnol Prog 6:7
39. Berry F, Sayadi S, Nasri M, Thomas D, Barbotin JN (1990) J Biotechnol 16:199
40. Godwin D, Slater H (1979) J Gen Microbiol 111:201
41. Hawley DK, McClure R (1983) Nuc Acd Res 1l:2237
42. Horwitz MSZ, Loeb LA (1990) Prog. Nuc. Acd, Res. Mol. Biol. 38:137
43. Kobayashi M, Nagata K, lshihama A (1990) Nuc Acd Res 18:7367
44. Reinikainen P, L/ihde M, Karp M, Suominen I, Markkanen P, M/ints/il/i P (1988) Biotechnol
Lett 10:149
45. Botterman J, H6fte H, Zabeau M (1987) J Biotechnol 6:71
46. Elvin CM, Thompson PR, Argall ME, Hendry P, Stamford NPJ, Lilley PE, Dixon NE (1990)
Gene 87:123
47. Shapira SK, Chou J, Richaud FV, Casadaban MJ (1983) Gene 25:71
48. Amann E, Brosius J, Ptashne M, (1983) Gene 25:167
49. Makoff AJ, Oxer MD (1991) Nuc Acd Res 19:2417
50. Perez L, Vega J, Chuay C, Menendez A, Ubieta R, Montero M, Padron G, Silva A, Santizo C,
Besada V, Herrera L (1990) Appl Microbiol Biotechnol 33:429
51. Vidal-Ingigliardi D, Raibaud O (1985) Nuc Acd Res 13:1163
52. Pistillo, JM, Vishwanatha JK (1990) Biochem Biophys Res Comm 169:1129
53. Lukacsovich T, Orosz A, Baliko G, Venetianer P (1990) J Biotechnol 16:49
54. Pharmacia (1991) Pharmacia LKB Biotechnology, Uppsala, Sweden
55. Deng T, Noel JP, Tsai MD (1990) Gene 93:229
56. Yasukawa K, Saito T (1990) Biotechnol Lett 12:419
57. Owolabi JB, Rosen BP (1990) J Bacteriol 172:2367
58. Flores N, de Anda R, Guerec L, Cruz N, Antonio S, Balbas P, Bolivar F, Valle F (1986) Appl
Microbiol Biotechnol 25:267
59. Dikshit KL, Dikshit RP, Webster DA (1990) Nuc Acd Res 18:4149
60. Scanlan DJ, Bloye SA, Mann NH, Hodgson DA, Carr NG (1990) Gene 90:43
61. Berg BL, Stewart V (1990) Genetics 125:691
62. Su TZ, Schweizer H, Oxender DL (1990) Gene 90:129
63. Metcalf WW, Steed PM, Wanner BL (1990) J Bacteriol 172:3191
64. Claverie-Martin F, Magasanik B (1991) Proc Natl Acad Sci USA 88:1631
65. Bingham RJ, Hall KS, Slonczewski JL (1990) J Bacteriol 172:2184
66. Poindexter K0 Gayle III RB (1991) Gene 97:125
67. Repoila F, Gutierrez C (1991) Mol Micorbiol 5:747
68. Lucht' JM, Bremer E (1991) J Bacteriol 173:801
69. Unden G, Trageser M, Duchene A (1990) Mol Microbiol 4:315
70. Terasawa N, Masayuki I, Uchida Y, Kobayashi M, Kurusu Y, Yukawa H (1991) Appl
Microbiol Biotechnol 34:623
71. Fischer M, Fytlovich S, Amit B, Wortzel A, Beck Y (1990) Appl Microbiol Biotechnol 33:424
72. Mott JE, Grant RA, Ho Y, Platt T (1985) Proc Natl Acad Sci USA 82:88
73. Iijima S, Kawai S, Mizutani S, Taniguchi M, Kobayashi T (1987) Appl Microbiol Biotechnol
26:542
74. Stanley KK0 Luzio JP (1984) EMBO J 3:1429
75. Chan WKY, Belfort G, Belfort M (1988) Gene 73:295
76. Chan WKY, Belfort G, Belfort M (1991) J Biotechnol 18:225
77. Claverie-Martin F, Diaz-Torres MR, Yancey SD, Kushner SR (1991) J Biol Chem 266:2843
78. Donovan WP, Kushner SR (1986) Proc Natl Acad Sci USA 83:120
79. Schmucker R, Gfilland U, Will M, Hitlen W (1989) Appl Microbiol Bioteehnol 30:509
80. Olsen MK, Rockenbach SK, Curry KA, Tomich C-SC (1989) J Biotechnol 9:179
81. Surek B, Wilhelm M, Hillen W (1991) Appl Microbiol Biotechnol 34:488
82. Spanjaard RA, van Dijk MCM, Turion AJ, van Duin J (1989) Gene 80:345
83. Br6ker M, Amann E (1986) Appl Microbiol Biotechnol 23:294
84. Ernst JF (1988) Trends Biotechnol 6:196
85. Maurizi MR, Trisler P, Gottesman S (1985) J Bacteriol 164:1124
86. Boss MA, Kenton JH, Wood CR, Emtage JS (1984) Nucl Acids Res 12:3791
87. Simon LD, Randolph B, Irwin N, Binowski G (1983) Proc Natl Acad Sci USA 80:2059
88. Kadokura H, Yoda K, Mitsunobu I, Yamasaki M (1990) Appl Environ Microbiol 56:2742
89. Wittliff JL, Wenz LL, Dong J0 Nawaz Z, Butt TR (1990) J Biol Chem 265:35
90. Quaas R, McKeown Y, Stanssens P, Frank R, B16cker H, Hahn U (1988) Eur J Biochem
173:617
76 K. Friehs and K. F. Reardon
137. Anba J, Baty D, Lloubes R, Pages J-M, Joseph-Liauzum E, Shire D, Roskam W, Lazdunski C
(1987) Gene 53:219
138. Hogset A, Blinsmo OR, Saether O, Gautvik VT, Holmgren E, Hartmanis M, Josephson S,
Gabrielsen OS, Gordelaze JO, Alestrom P, Gautvik KM (1990) J Biol Chem 265:7338
139. Wadenstein H, Ekebacke A, Hammarberg B, Holmgren E, Kalderen C, Tally M, Moks T,
Uhlen M, Josephson S, Hartmanis M (1991) Bioteehnol Appl Bioehem 13:412
140. Fraser TH, Bruce TJ (1987) Proc Natl Aead Sci USA 75:5936
141. Kornacker MG, Pugsley AP (1990) Molec Microbiol 4:1101
142. Talmadge K, Brosius J, Gilbert W (1981) Nature 294:176
143. Morioka-Fujimoto K, Marumoto R, Fukuda T (1991) J Biol Chem 266:1728
144. Legoux R, Leptlatois P, Joseph-Liauzun, E (1991) US Patent 4945047
145. Letoffe S, Delepelaire P, Wandersman C (1991) J Bacteriol 173:2160
146. van Putten AJ, de Graaf FK, Oudega B (1987) Proc 4th Europ Con Biotechnol 4:593
147. Mackman N, Nicaud JM, Gray L, Holland IB (1986) Cur Top Microbiol Immunol 125:159
148. Lazzaroni JC, Portalier RC (1982) Eur J Appl Microb Biotechnol 16:146
149. Arian D, Portalier RC (1984) Eur J Appl Microb Biotechnol 19:5
150. Mosbach K, Birnbaum S, Hardy K, Davies J, Bulow L (1983) Nature (London) 302:543
151. Murakami Y, Furusato Y, Kato C, Habuka N, Kudo T, Horikoshi K (1989) Appl Microbiol
Biotechnol 30:619
152. Bialy H (1987) Bio/Technology 5:883
153. Thatcher DR (1990) Biochem Soc Trans 18:234
154. Fischer G, Schmid FX (1990) Biochem 29:2205
155. Dill KA (1990) Biochem 29:7133
156. Jaenicke R (1987) Prog Biophys Molec Biol 49:117
157. Hanada K, Yamamoto I, Anraku Y (1987) J Bio Chem 262:14100
158. Maschke H-E, Hebenbrock K, Friehs K (1990) DECHEMA Biotechnol Conf 4:379
159. Ullman A (1984) Gene 29:27
160. Maina CV, Riggs PD, Grandea AG, Slatko BE, Moran LS, Tagliamonte JA, McReynolds LA,
di Guan C (1988) Gene 74:365
161. Hochuli E, Bannwarth W, D6beli H, Gentz R, Stfiber D (1988) Bio/Technology 11:1321
162. Smith JC, Derbyshire RB, Cook E, Dunthorne L, Viney J, Brewer SJ, Sassenfeld HM, Bell LD
(1984) Gene 32:321
163. Persson M, Bergstrand MG, Billow L, Mosbach K (1988) Anal Biochem 172:330/337
164. Stiegelmeier C, Friehs K, Schilgerl K (1992) Acta Biotechnol (in press)
Modeling and Control for Anaerobic
Wastewater Treatment
Dedicated to Professor Dr. Karl Schfigerl on the occasion o f his 65th birthday
List of Symbols . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
2 Recent Reviews of the Field . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
3 Reactions in Anaerobic Digestion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
4 Dynamic Mass Balance Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
4.1 Biomass Balance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
4.2 Substrates and Product Balances . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
5 Stoichiometry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
6 Kinetics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
7 Ion Charge Balance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
8 Modeling Biofilm Kinetics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
9 Gas Phase Modeling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
10 Comparison of Simulations with Experiment . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
10.1 Steady State Comparison of Experimental and Simulated Carbon Balances . . . . . . 96
10.2 Comparison of Dynamic Experimental and Simulation Results . . . . . . . . . . . . . . . . 96
11 On-line Measurement and Observation Methods . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
12 Control of Anaerobic Processes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
12.1 Simple Controller Design Using Simulations of a Complex Process Model . . . . . . . 100
12.2 Proportional-Integral-Differential Control of pH . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
12.3 PID Control of Organic Acid Concentrations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
12.4 Control of Dissolved Hydrogen Concentration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
12.5 Control Based on the Methane Content of the Gas . . . . . . . . . . . . . . . . . . . . . . . . . 105
13 Adaptive Control and Optimization of Anaerobic Processes . . . . . . . . . . . . . . . . . . . . . . 106
14 Concluding Remarks . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 ll
15 References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1t 2
The recent literature on the modeling and control of anaerobic wastewater treatment processes is
reviewed. An example from the author's personal work is used to describe how a dynamic simulation
model can be developed from the basis of multi-organism growth kinetics and mass balancing
techniques. This included consideration of the organic acid dissociation equilibria important for pH
calculation and the thermodynamic influence of hydrogen on the reactions involving propionic and
butyric acids. Liquid phase balances were linked to gas phase balances by gas-liquid transfer
considerations. It is shown in detail how the model was applied to one and two stage experimental
reactors for the design and tuning of controllers. Both conventional PID controllers and adaptive
optimizing controllers, employing simple input-output models with an objective function, were
tested.
List of Symbols
Greek symbols:
[*] Gain adaptive optimizer
[3 [,] Working point optimizer
Modeling and Control for Anaerobic Wastewater Treatment 81
Indices:
AH acid
A- dissociated AH
Ac acetic acid
An- anion
B base
Bu butyric acid
D differential
g gas phase
HAc undissociated acetic acid
HBu undissociated butyric acid
HPr undissociated propionic acid
i refers to component
I integral
K+ cation
L liquid phase
P product
Pr propionic acid
S substrate
sl solid
titr titrator base
tot total
W water
X biomass
Z surplus cations
0 reactor feed stream
Abbreviations:
COD Chemical oxygen demand
GC Gas chromatograph(y)
P Proportional
PI Proportional-integral
PID Proportional-integral-differential
VFA Volatile fatty acids
82 E. Heinzle,I. J. Dunn and G. B. Ryhiner
1 Introduction
Table 1. Rate parameters as reported from the literature by Pavlostathis and Giraldo-Gomez [41].
VSS - volatile suspended solids
Substrate
~maxCs~
k Ks gmax Y
KsY /
g COD g COD g VSS
d-1
g VSS d m3 g COD
Polymers (Proteins,
fats, carbohydrates ] '~
Monomers (sugars, [
amino acids,...) [
IHBu - =
Fig. 1. Reaction scheme of anaerobic degradation. Dashed arrows indicate gaseous compounds
transferred out of the liquid phase. HPr - propionic acid; Pr- - propionate; HBu - butyric acid;
Bu - butyrate; HAc - acetic acid; Ac- - acetate
Table 2. Most important reactions in the anaerobic degradation of whey. Reaction R1 is given in a
C-mol basis (one carbon atom per molecule); vi stoichiometric coefficients.
(AGo' f o r T = 2 5 ~ p H = 7 , p = l b a r , xn2o = 1, c = l m o l k g - 1 )
Acetification
CH3(CH2)2COO- 4- 2HzO ,~-2CH3COO- + 2 H z + H + AG ~ + 71.7kJ (R2)
CH3CH2COO + 2 H20 ~- CH3COO- + CO 2 -~ 3 H 2 AG O' = + 48.3 kJ (R3)
Methanogenesis
CH3COO- + H + CH 4 4- CO 2 AG O' = - 35.8 kJ (R4)
CO 2 + 4 H 2 ~- CH 4 + 2 H20 AG O' = - 130.7 kJ (R5)
R e f e r r i n g t o T a b l e 2, t w o t y p i c a l p r o c e s s c a s e s m a y b e d i s t i n g u i s h e d .
If solid matter has to be degraded, then the hydrolysis step may be rate
l i m i t i n g , a n d o v e r l o a d i n g will n o t b e i m p o r t a n t .
- O v e r l o a d i n g o f t h e r e a c t o r m a y o c c u r if e a s i l y d e g r a d a b l e m a t e r i a l is t o b e
t r e a t e d . I t is a g r e e d t h a t h y d r o g e n p l a y s a n i m p o r t a n t r o l e i n t h e p r o c e s s
s t a b i l i t y . D u e t o t h e p o s i t i v e free e n e r g y o f r e a c t i o n , r e a c t i o n s R 2 a n d R 3 in
86 E. Heinzle,I. J. Dunn and G. B. Ryhiner
Table 2 can only proceed to the right, if the products, especially H2, are
continuously removed. Thus if the final methanogenesis steps are inhibited due
to toxicity or pH then the hydrogen and acetic acid accumulation will stop the
reactions, R2 and R3. Because the concentration of H 2 influences the equilib-
rium relation to the second (R2) and even third power (R3), it is evident that
the hydrogen partial pressure must be kept sufficiently low. This is accom-
plished by the methanogenic scavengers, whose rates cannot keep up in the
event of substrate overloads. H2 consuming organisms can remove H 2 locally
to concentration levels which make the reaction proceed exergonically [11].
because only the Syntrophobacter can use it. One thing is clear from the new
modeling concepts presented in this paper: many mysteries still remain with
regard to the microbiology and its kinetic modeling.
Smith and McCarty [57] used a four bacteria-type model for mixed
substrates ethanol and propionate. Two types were used for oxidation of the two
substrates to acetate and the other two types for methanogenesis. Thus four
component balances and four biomass balances were used.
Costello et al. [12] used the Mosey model [-39] as a basis but also included
lactic acid as an important intermediate. Six groups of bacteria were considered.
One group formed acetic, butyric or lactic acid and a second degraded lactic
acid into propionic or acetic acid. Two groups were acetogenic and formed
acetic acid from butyric and from propionic acid. Methane was formed by two
groups, one by hydrogen-consuming bacteria reducing carbon dioxide and
another from acetate.
A four organism Mosey-type model was used by Dochain et al. [19] to
develop an adaptive control algorithm for the control of hydrogen concentra-
tion. A quasi-steady state assumption regarding the glucose concentration
allowed relating the accumulation rate for hydrogen to the glucose inflow rate
and the hydrogen outflow gas rate. A further assumption allowed neglecting the
propionate reaction path. For the simulated results both hydrogen and acetate
were assumed to have inhibiting effects on the methanogenic steps and the
originally neglected propionate reaction was included. Without control the
hydrogen-consuming methanogens would be washed out of the reactor.
Ryhiner [52] used the reaction scheme shown in Fig. 1 and Table 2 including
a six organism group model to describe anaerobic digestion of whey.
For the purpose of modeling, many reactors used in anaerobic degradation can
be considered as well mixed. The basic balance equation for a well mixed liquid
phase is
mass transfer from -
F accumulation ~ = input - output + reaction + the gas to the
[_of component i j
liquid phase
d C L.
V L d ~ = FL (CLio -- CLi) + ri VL + Ni VL (1)
dX i Xi
dt - rx, - Xi'kd, -- - - (2)
Zsl
where
rxl = ~iXi (3)
where i indicates populations feeding on the various substrates (whey, butyric,
propionic and acetic acid, hydrogen). The term Xikd~ also accounts for mainte-
nance and endogenous reactions.
Each substrate may be consumed and produced by several reactions which are
taken into account by the summation. The corresponding balances for volatile
components (H2, CO2 and CH4)include mass transfer term as shown in Eq. (1).
The relation of substrate consumption is
-- ~iXi (5)
rs~ -- Yxi/si
t'tiXiYPjsi (6)
rvj -- Yx~/si
5 Stoichiometry
hydrolysis reactions are not known a priori. Experiments under defined condi-
tions are necessary to determine these coefficients. The six unknown stoichio-
metric coefficients of the acidification of whey were determined after measure-
ment of acid production and COD reduction in the acidification reactor of the
two stage system [52]. In this reactor no methane was produced, and reactions
R2 and R3 were assumed not to proceed under these conditions. The experi-
mentally determined values were: vBu = 0.5; Vpr = 0.12; VAc = 0.15; VCO2 = 0.23;
Vn2 = 0.24; Vn20 = 0.095.
More difficult is the estimation of biomass yield coefficients of individual
reactions. One successful approach uses ATP balancing. The yield of biomass
from ATP (Yx/ATP) is quite constant for many organisms [-4, 481. The average
value given there is YXmTP ~ 10 g tool 1. If the production of ATP from
metabolic reactions is known, the biomass yield can be estimated in a straight-
forward manner. Smith and McCarty [-57] developed these stoichiometric
coefficients from the free energy of the catabolic reactions. Growth was related
to energy availability by calculating the difference between the free energy of the
catabolic reaction and the energy involved in the conversion of the substrates to
pyruvate. It was noted that at least 6 reactions occurred during the experiments,
which involved reduced intermediates that were not included in the model.
Erickson [21] presented an electron balance and discussed it as a means of
estimating biomass yield coefficients in anaerobic processes.
6 Kinetics
~maxi Csi
lai - Cs i + Ks ~ (7)
or with an additional substrate inhibition term (acetic acid)
I-tmax~Csi
gi = Cs i + Ks ~ + Csi2/Kii (8)
The acidic form was taken to be the kinetically determining form, as was done in
previous work [14, 32, 69]. The concentrations of the actual substrates (un-
dissociated carboxylic acids) are pH dependent. At neutral pH-values, only
C O 2 - / H C O 3 - and in sulfur-containing wastewaters HzS/S- have a reasonable
buffer capacity. All other acids mainly exist in their salt form.
In the acetogenic step, acetic acid, hydrogen and carbon dioxide are
produced from propionic and butyric acid. The thermodynamic limits for these
90 E. Heinzle, I. J. D u n n and G. B. Ryhiner
In the pH range of interest (around pH = 7) all strong acids and strong bases are
completely dissociated. Weak acids and bases are only partly dissociated.
AH~A- +H + (10)
The degree of dissociation is determined by the dissociation constant as
CA_H +
K A -- - - (11)
CA
where: C A concentration of the undissociated acid; C A - concentration of the
corresponding base (dissociated acid).
Moderately strong acids occurring in anaerobic systems are: acetic:
KAc ---- 1.73 x 10 - 2 [ m o l m - 3 ] , propionic: Kpr = 1.31 x 10 -2 [mol m - 3 ] ,
butyric: KAc = 1.44x 10 -2 [mol m-3], carbon dioxide: KAc = 3.02x 10 -2
[mol m-3], hydrogen sulfide: KH2s = 1.26 x 10 -4 [mol m -3] and ammonium:
KNH4+ = 5.28 x 10 -7 [mol m-3]. The only moderately strong base is ammonia:
KNH3 = 1.85 X 10 .2 [m01 m-3]. At pH values below 7 ammonia is practically
completely dissociated (cNn3-< CNH4+*0.005). Correspondingly, the buffering
effect of ammonia is negligible. The existing cations in anaerobic degradation
process are therefore cations (K +) from strong bases (e.g. Na +, K +, NH4+,
Ca2 +). In the expected pH range one always has ZCK+ >>CH+,where ZCK+ is the
total cation concentration. Negative ions are mainly from strong acids (e.g. C1-,
SO #-) and from weak acids (Ac-, Pr-, Bu-, HCOj). The concentration of
CO 2- is always much smaller than that of HCO 3 (e.g. at pH = 7;
Cco~-/CHco~ ~ 4 x 1 0 - 4 ) . From this one obtains
KA.
2 K,~Kw -C, ..... + 2 CK+ = 2 KAi~'CH + CA.... ~+ ~CAn- (12)
KB~ + CH+
KA~ Acid dissociation constant (e.g. Knc); KBI - Base dissociation constant (e.g.
-
measured gas concentrations were usually below 2%. The corresponding H2S
concentration in the liquid phase was therefore always < 0.6 [mol m-3].
Neutralization reaction is a very fast reaction reaching equilibrium almost
instantaneously. During pH control strong base is added. The addition of strong
alkali causes an increase in ZCK+ which, by solving Eq. (12), results in a decrease
of CH+.
It is possible to use only the difference between cations and anions
Cz = ~CK+ -- 2CAn (13)
For the numerical solution it is useful to rearrange Eq. (12) and use Eq. (13) in
the form,
Ki KBi
6 = ~ Ki + Cn+ Ctot,~- Z Kw C, ..... - C z (14)
KBi +
CH+
10
pI-I
Titration
Cation :1
Dynamic Cation __.~
Cation uent " - ' - - ~ M a s s Balance
Ion Charge
Balance
(Algebraic
equation)
Acid Feed Dynamic Acid M a s s Balances
(Total m a s s of week acids and
Acid Effluent b a s e s - acetic, CO2, NH3,...)
Production,
Consumption
CA
r=f(CA) [~ Concentration
of free acids
Fig. 3. Schematic of the ion charge balance iteration within a dynamic simulation. C A - free acid
concentration
vdCz
L dt = FL (CZ~ - C Z ) + Ftitr C Z .... (15)
Figure 3 gives an information flow diagram of how the solution of the dynamic
mass balances depends on the reaction kinetics. The kinetics in turn depend on
the undissociated acids, which are given by the total acids and the pH, according
to the equation. The charge balance represents an algebraic loop in the
otherwise sequential dynamic integration.
In those reaction systems where biofilms or flocs are formed (UASB, fluidized
bed, fixed bed) diffusion and direct inter-species transfer of intermediates (e.g.
hydrogen) may be important. The detailed biofilm structure and formation is
Modeling and Control for Anaerobic Wastewater Treatment 95
very complex which can easily be seen in electron microscopic pictures. It can be
assumed that the biofilm structure is also dynamic but with longer time
constants than the hydraulic retention time. Quantitative data on biofilm
structure and inter-species hydrogen transfer are practically not existing [36].
It is also difficult to determine a representative diffusion coefficient. Examples
of biofilm modeling in the literature are: Denac et all [17]; Bryers [8]; Wang
et al. [65].
XiUBA~gp
dni - VLkLa(XiHip - cLi) (16)
dt RT
In addition the time delay of the gaseous components due to the head space was
modeled with a total gas balance. The gas liquid mass transfer coefficient, kLa,
was modeled as a function of fractional gas hold-up, which was a function of the
gas production rate. Mather [33] used a similar concept to model the gas phase
of an anaerobic reactor.
The model of Ryhiner [52] was used to compare simulated results with
experimental data from whey treatment. Steady state and dynamic comparison
was made as given below.
6.2 .8
pH
Gas Production
6.1
._~
-6 E~
6.0
"1- -4 n"
O- 5.9
20
5.8
Y
"O
o
2 IX.
5.7-
5.6 9 . 9 ~ 9 9 9 ~ 9 9 9 ~ 9 . . ~ .
1 2 3 4
Time [hi
Fig. 4. Response of pH and gas production rate on concentration step change (1.4 kg m -3
--* 2.0 kg m 3 whey powder) in single stage reactor (experiment)
Modeling and Control for Anaerobic Wastewater Treatment 97
6.2
pH
6.1, Gas Production
i
6 .c_
E
6.0
"-r
5.9 4 rr
g
5.8
StepC h ! n ~ o
2 13.
5.7 69
(.g
0 1 2 3 4
Time [h]
Fig. 5. Response of pH and gas production rate on feed concentration step change (1.4 kg m -3
2.0 kg m- 3 whey powder) in single stage reactor (simulation)
68 4O
~176
64
36 o'-s
.. 62-
34
O
60 84
32
58
StepChange ~ CH4
CO~
56 9 9 , , , , . 9 9 , 9 9 . 30
O 1 2 3 4 5
Time [h]
68 9 40
66
38
64
-2.'
36
o "6
62
-r-
0 -34 0
60 C)
. 32
58 .l CH4
Step Change ~ C(~
5 6 9 . 9 , . 9 . , 9 9 . , 9 9 , 9 . . 3 0
0 1 2 3 4
Time [h]
Fig. 7. Response of gas composition to a feed concentration step change increase (1.4 kg m -3
--, 2.0 kg m - 3 whey powder) in single stage reactor (simulated)
9e-5
8e-5
E
"6 7e-5
E
-'r-
68-5
t
Step Change
5 e - 5 '
1 2 3 4
Time [h]
Fig.& Response of dissolve hydrogen to a feed concentration step change increase (1.4kg m -3
2.0 kg m- 3 whey powder) in single stage reactor (simulated)
tial, sulfide ion activity, gas production rate, gas composition, hydrogen gas, and
carbon monoxide gas. Bicarbonate alkalinity is determined by acid titration to
defined p H between (4.0 _< pH _< 5.75), as discussed in detail by Powell and
Archer [44] and Di Pinto et al. [43]. Switzenbaum et al. [60] felt that carbon
monoxide and hydrogen gas responded sufficiently fast to overloads to make
them useful for control. A later report by the same authors [27] described their
experimental work on carbon dioxide and hydrogen responses and concluded
that CO gave information on the acetate consumption, while Hz was involved
mainly in the CO2 reduction step. Pauss et al. [40] measured continuously the
Modelingand Controlfor AnaerobicWastewaterTreatment 99
Table 4. (continued)
identical. For the two stage system control the controller gain in the experiment
was set at 40% of the value found suitable by simulation. In spite of this, some
instability was observed, possibly due to the higher buffering in the experiment
(twice the amount of dissolved CO2).
Modeling and Control for Anaerobic Wastewater Treatment 103
,omass . nr, s
Balances Law I
Con,ro,,er~ ~c.,
Fig. 9. Calculation flow diagram for the simulation of on-line control of the anaerobic process
6.0 10
Setpoint 8
5.9
-'T
.E
6 E
-r
Q_ 5.8 0~
4
o
5.7 E
2
pH
FlowRate
5,6 i
0 1
Time [hi
Fig. 10. PI-Control of pH after feed concentration step change (1.4 kg m 3 __,2.8 kg m -3 whey
powder), single stage reactor (experiment). ( K r = 0.7 m 6 mol- 1 h- 1, .q = 0.2 h)
6.0 10
Setpoint
J
'8
5.9,
.=_
.6 E
5.8.
Step C~hang!
I
I
'4 rr
F-
o
5.7. It
pH
Flow Rate
5.6 9 r i 0
1 2 3
Time [h]
Fig. 11. PI-Control of pH after feed concentration step increase (1.4 kg m=3---~ 2.8 kg m -3 whey
powder), single stage reactor, simulation. (Kr = 0.6 m 6 tool- 1 h - 1, Xl = 0.2 h, % = 0.05 h)
50 25
t r, --.-- Acet,0Ac,f0
= = -- -- _ _ -~
<O
15
(g
o 20. rr
o
10 ~-
10,
0 9 9 9 i 9 9 9 i 9 , 9 i , 9 9 i 9 9 9 i 9 9 9 i 9 9 9 5
0 1 2 3 4 5 6 7
Time [h]
Fig. 12. Acetic acid set-point, PID-control of the two stage system, simulated. (K r = 3 • 10 -4,
~ = 0.25, zD = 0.2). Step change concentration increase (1.4 kg m 3 ~ 4.2 kg m -3 whey powder)
was set somewhat high in the experiment and caused oscillations, Similar results
were obtained with propionic acid control.
40 - 35
/
9 Acetic Acid
30 ~ . - Flow Rate 9 30
-- c
20 20
"~ 0~
< 15
0 10'
0'
5
-10 , 9 9 , 9 9 9 ~ 9 9 , ~ 9 9 9 , 9 9 9 ~ 9 9 9 t 9 9 9 0
0 1 2 3 4 5 6 7
Time [hi
Fig. 13. Acetic acid set-point, PID-control of the two stage system, experimental. (K r = 3 x 10 -4,
zl = 0.3, ~D = 0.05). Step change concentration increase (1.4 kg m - 3 ~ 4.2 kg m - 3 whey powder)
le-4 14
H 2
E
6e-5
JL_ _-==== .
Set Point
. . .
10
"-
i
c
a: o ~o
2e-5
4
~ ] Feed Concentration I
0e+0 9 . 9 ~ 9 9 9 ~ 9 9 9 , 9 9 9 ~ 9 9 9 , 9 9 9 , 9 , 9 2
1 2 3 4 5 6
Time [h]
Fig. 14. Simulated response of the PI hydrogen controller to a square wave in feed concentration for
the single stage system. (K r = 2,5 m 6 mol- ~ h- 1 ~ = 0.2 h). Step change concentration increase
(1.4 kg m -3 ~ 2.8 kg m 3 whey powder)
on the methanogenic steps and the neglected propionate reaction was included.
Without control the hydrogen consuming methanogens are washed out of the
reactor. Numerous simulations tested the controller under a range of conditions,
including a comparison with simple P! control. This approach appears to be a
promising way of including some of the model complexities into an adaptive
controller design.
Detailed mechanistic models for adaptive control are generally too complex,
and would require simplification before implementation. However they have the
advantage of containing most information and may be necessary if not all
108 E. Heinzle, I. J. Dunn and G. B. Ryhiner
variables can be measured. When all the important process variables can be
measured, then empirical models of a simple, linear form can be used to
represent the relations between manipulated inputs and measured outputs [493.
Feed
Changes
F ] Reactor ] Measured
~_~Linear I Data f
Process Model ~ Compare
[Updated
IPar~m.
I I Parameter l_
Estimation ]-
F ]Optimizer l- Objective
] ~ Function
Fig. 15. Continual parameter updating allows the use of a simple model. The optimizer calculates
the flow rate required to reach the maximum
Modeling and Control for Anaerobic Wastewater Treatment 109
concentration. In this application the aim was to maximize the rate of anaerobic
degradation without allowing the organic acid concentrations to become ex-
cessively high. Thus the performance index, PX, was formulated as follows:
Methane
Total 1
[Objective~ = production organic acids
- Constant x
L function J rate concentration
PI = Yl - 13 x Y2 (19)
where Yl = rCH4and Y2 = '~CAi' The shape of the performance index was found
by the use of the differential equation simulation model of the process. From this
the influence of feed rate on the performance index at steady state could be
determined. This function, as shown in Fig. 16 is dependent on feed concentra-
tion and feed rate.
Simulation results of the reactor with the estimation (Eqs. 17 and 18)
indicated that the simplest model was sufficient for obtaining good estimates of
Yl and Y2,
yl(t) = b1,1 u(t - 1 - d~) + cl (17)
Y2(t) = b2,1 u(t - I - d2) + c2 (18)
This was therefore used in the experiments. This simple model considers only the
most recent output data directly, and although much of the dynamic informa-
tion was not considered by this simplification, this was found to be of no
disadvantage. Influence of older data was included by variable forgetting factor
for parameter values.
Simple input-output models were used for the adaptive optimization of one
and two stage anaerobic fluidized bed reactors.
15
.1[ " yl : Methane rate [mmol h -1]
10
~= 0
-5 ~. PI = y l - 0.5 y2
9 PI = y l - l . 0 y2
-10 -------13----PI = yl-l.5 y2
o i ~
Feed rate [1h -1]
Fig. 16. The performance index depends on the penalty function constant for the organic acid
concentration, here for a single stage system, simulated (Feed concentration - 42 C-mol m 3)
1 10 E. Heinzle, I. J. D u n n and G. B. Ryhiner
0.007
AC
0.006
NC
7.00e-4 _
0.005
6.00e-4
r
90.004
5.00e-4 _
4.00e4. NC 0.003
3.00e-4,
2.00e-4 : i i
0 20 40 60 80
T i m e [hi
Fig. 17. Simulated optimizer response of the single stage process to a step change (42 to 63
C-mol m - 3) in the feed concentration. AC - adaptive control, N C no control. (~ = 10- 6, 13 = 1.0,
Zo.1 = 10-v, •o,2 = 10 7, ~min = 0.1)
Modelingand Control for Anaerobic WastewaterTreatment 111
0.8
0,7
e-,
1 6
~D
0,6
5 .~=
0.5
E
0.4 l ~D
2
0.3
Feed concentration [ I
0.2
'2 4 ; 8 1'0 1'2
Time [h]
Fig. 18. Controlled response of the single stage process to a step change (42 to 63 C-tool m 3) in the
feed concentration (experiment). (~ = 2 x 10 -6, 13 = 1.0 Z o 1 = 10-6, Zo,2 = 10 -6, ~min : 0.1)
considerable delay. The results obtained using the adaptive controller are
compared to uncontrolled system response. A feed concentration increase by
75% (42 --, 73.5 C-mol m-3) was handled much better by the adaptive control-
ler. Comparison of response to concentration decrease (time = 40 h) is difficult
from Fig. 17 since the starting points were different. In the experiment (Fig. 18)
the controller reacted more quickly in both cases. The flow rate decrease by the
controller corresponded nearly to keeping constant the substrate loading rate of
the reactor. The experimental results for the two stage reactor were similar to the
simulated results but had a significant time delay and correspondingly larger
overshooting.
It is quite remarkable that the detailed mechanistic simulation model was so
useful for designing and testing, the adaptive controller. The controlled process
could be simulated for all controller-types, allowing the controller parameters to
he found. The performance index as defined was satisfactory for control and the
algorithm from the literature [49] was usable.
14 Concluding Remarks
15 References
1. Andrews JF (1978) The development of a dynamic model and control strategies for the
anaerobic digestion process. In: Mathematical models in water pollution control. J Wiley,
Chichester, p 281
2. Archer DB (1983) Enzyme Microb Technol 5:570
3. Archer DB, Hilton MG, Adams P, Wiecko H (1986) Biotechnol Lett 8:570
4. Atkinson B, Mavituna F (1983) Biochemical engineering and biotechnology handbook. Nature
Press, New York
5. Barnes D, Fitzgerald PA (1987) Anaerobic wastewater treatment processes. In: Forster CF,
Wase DAJ (eds) Environmental biotechnology. Ellis Horwood, Chichester, p 57
6. Bastin G, Dochain D (1990) On-line estimation and adaptive control of bioreactors. Elsevier,
Amsterdam
7. Bellgardt KH (1991) personal communication
8. Bryers JD (1985) Biotechnol Bioeng 27:638
9. Colberg PJ (1988) Anaerobic microbial degradation of cellulose, lignin, oligoligninols, and
monoaromatic lignin derivatives. In: Zehnder AJB (ed) Biology of anaerobic microorganisms.
Wiley, New York, p 333
10. Collins AS, Gilliland BE (1974) J Environ Eng Div 100:487
11. Conrad R, Schink B, Phelps TJ (1986) FEMS Microbiology Ecology 38:353
12. Costello DJ, Greenfield PF, Lee PL (1991a) Water Res 25:847
13. Costello DJ, Greenfield PF, Lee PL (1991b) Water Res25:859
14. Dalla Torre A, Stephanopoulos G (1986) Biotechnol Bioeng 28:1106
15. Denac M, Dunn IJ (1988a) Biotechnol Bioeng 32:159
16. Denac M, Miguel A, Dunn IJ (1988b) Biotechnol Bioeng 31:1841
17. Denac M, Lee PL, Newell RB, Greenfield PF (1990) Water Res 24:583
18. Dochain D, Bastin G (1985) Environ Technol Lett 6:584
19. Dochain D, Perrier M, Pauss A (1991) Ind Eng Chem Res 30:129
20. Dolfing J (1988) Acetogenesis. In: Zehnder AJB (ed) Biology of anaerobic microorganisms.
Wiley, New York, p 417
21. Erickson LE (1988) Bioenergetics and yields for anaerobic digestion. In: Erickson LE and Fung
DY-C (eds) Handbook on anaerobic fermentations. Dekker, New York, p 325
22. Erickson LE, Fung DY-C (1988) Handbook on anaerobic fermentations, Dekker, New York
23. Fan T, Erickson LE, Baltes JC, Shah PS (1973) J Water Pollut Control Fed 45:591
24. Golden MP, Ydstie BE (1989) AIChE J 35:1157
25. Gujer W, Zehnder AJB (1983) Water Sci Technol 15:127
26. Harper SR, Pohland FG (1986) Biotechnol Bioeng 28:585
27. Hickey RF, Switzenbaum MS (1991) J Water Pollut Control Fed 63:129
Modeling and Control for Anaerobic Wastewater Treatment l 13
28. Kaspar HF (1977) Untersuchung zur Kopplung von Wasserstoff- und Methanbildung im
Faulschlamm. Diss ETH, Zfirich
29. Kaspar HF, Wuhrmann K (1978) Appl Environ Microbiol 36:1
30. Kissalita WS, Lo KV, Pinder KL (1989) Biotechnol Bioeng 33:623
31. Lloyd D, Whitmore TN (1988) Lett AppI Microbiol 6:179
32. MS~rklH, Mather M, Witty W (1983) Mess- und Regeltechnik bei der anaeroben Abwasserrein-
igung sowie bei Biogasprozessen. In: Rehm H-J, Reed G (eds) Biotechnology Vol. 2,
R. Oldenburg, M/inchen, p 369
33. Mather M (1986) Mathematische Modellierung der Methang~irung. In: Fortschrittsberichte-
VDI, Reihe 14: Landtechnik/Lebensmitteltechnik Vol. 28, VDI Verlag, Dfisseldorf
34. McCarty PL (1981) One hundred years of anaerobic treatment. In: Hughes ef al. (eds) Anaerobic
digestion. Elsevier, Amsterdam, p 1
35. McCarty PL, Smith DP (1986) Environ Sci Technol 20:1200
36. McCarty PL, Mosey FE (1991) Water Sci Technol 24:17
37. McInerney MJ (1988) Anaerobic hydrolysis and fermentation of fats and proteins. In: Zehnder
AJB (ed) Biology of anaerobic microorganisms. Wiley, New York, p 373
38. Moletta R (1989) Environ Technol Lett 10:173
39. Mosey FE (1983) Water Sci Technol 15:209
40. Pauss A,'Beauchemin C, Samson R, Guiot S (1990) Biotechnol Bioeng 35:492
41. Pavlostathis SG, Giraldo-Gomez E (1991) Water Sci Technot 24:1519
42. Petersen JN, Whyatt GA (1990) Biotechnol Bioeng 35:114
43. Di Pinto AC, Limoni N, Passino R, Rozzi A, Tomei MC (1990) Anaerobic process control by
automated bicarbonate monitoring. In: Briggs R (ed) Instrumentation, control and automation
of water and wastewater treatment and transport systems. Pergamon, Oxford, p 5l
44. Powell GE, Archer DB (1989) Biotechnol Bioeng 33:570
45. Da Pra E (1987) Rechnerunterstfitzte Optimierung des Abbaus yon Rfibenschwemmwasser
mittels eines Anaerobfilters Diss. Univ, Zfirich
46. Renard P, Dochain D, Bastin G, Naveau H, Nyns EJ (1988) Biotechnol Bioeng 31:287
47. Ripley LE, Boyle WC, Converse JC (1986) J Water Pollut Control Fed 58:406
48. Roels JA (1983) Energetics and kinetics in biotechnology. Elsevier Biomedical Press, Amsterdam
49. Rolf M J, Lim HC (1984) Chem Eng Commun 29:1236
50. Rozzi A, Eng M (1984) Trans Inst M C 6:153
51. Rozzi A, Di Pinto AC, Brunetti A (1985) Environ Technol Lett 6:594
52. Ryhiner G (1990) Regelung und adaptive On-Line Optimierung von ein- und zweistufigen
anaeroben Fliessbett-Reaktoren. Diss ETH, Zfirich
53. Ryhiner G, Dunn IJ, Heinzle E, Rohani S (1992) J Biotechnol 22:89
54. B61aich JP, Bruschi M, Garcia JL (1990) Microbiology and biochemistry of strict anaerobics
involved in interspecies hydrogen transfer. FEMS Symp. No. 54. Plenum New York
55. Schfirb~scher D, Wandrey Ch (1991) Anaerobic waste water process models. In: Rehm HJ, Reed
G (eds) Biotechnology Vol. 4, VCH, Weinheim, p 441
56. Shi Z, Shimizu K, Watanabe N, Kobayashi T (1989) Biotechnol Bioeng 33:999
57. Smith DP, McCarty PL (1990) J Water Pollut Control Fed 62:175
58. Stephanopoulos G (1988) Modelling mixed culture interactions in anaerobic digestion. In:
Erickson LE and Fung DY-C (eds) Handbook on anaerobic fermentations. Dekker, New York,
p 597
59. Stouthamer AH (1988) Bioenergetics and yields with electron acceptors other than oxygen. In:
Erickson LE and Fung DY-C (eds) Handbook on anaerobic fermentations. Dekker, New York,
p 345
60. Switzenbaum MS, Giraldo-Gomez E, Hickey RF (1990) Enzyme Microb Tecbnol 12:722
61. Thauer, Jungermann, Decker (1977) Bacteriol Rev 41:147
62. Tschui M (1989) Dynamisches Verhalten der mesophilen anaeroben Schlammstabilisierung.
Dissertation. ETH, Zfirich
63. Vogels GD, Keltjens JT, Van der Drift C (1988) Biochemistry of methane production. In:
Zehnder AJB (ed) Biology of anaerobic microorganisms. Wiley, New York, p 707
64. Wandrey C, Aivasidis A, Bastin KH (1986) Ann NY Acad Sci 469:421
65. Wang YT, Suidam MT, Rittman BE (1986) J Environm Eng 112:975
66. Whitmore TN, Lloyd D (1986) Biotechnol Lett 8:203
67. Whitmore TN, Lazzari M, Lloyd D (1985) Biotechnol Lett 7:203
114 E. Heinzle et al.: Modeling and Control for Anaerobic Wastewater Treatment
68. Whitmore TN, Jones G, Lazzari M, Lloyd D (1987) Methanogenesis in mesophilic and
thermophilic anaerobic digesters: Monitoring and control based on dissolved hydrogen. In:
Heinzle E, Reuss M (eds) Mass spectrometry in biotechnological process analysis and control.
Plenum, New York, p 143
69. Wiesmann U (1988) Chem Ing Tech 60:464
70. Zehnder AJB (1988) Biology of anaerobic microorganisms. Wiley, New York
71. Ricker NL, Slater WR, Merigh M (1988) Process monitoring and control strategies for
anaerobic wastewater treatment. DECHEMA Biotechnol Conf Vol. 2, p 413
72. Costello DJ, Lee PL, Greenfield PF (1989) Bioprocess Eng 4:119
Design of Reactors for Plant Cells and Organs
Pauline M. Doran
Department of Biotechnology, University of New South Wales,
P . O . B o x 1, K e n s i n g t o n N S W 2033, A u s t r a l i a
1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 117
2 Suspended Plant Cells . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 118
2.1 Properties of Plant Cultures Relevant to Reactor Design . . . . . . . . . . . . . . . . . . . 118
2.1.1 Rheology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 118
2.1.2 Aggregate F o r m a t i o n . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 120
2.1.3 Shear Sensitivity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 121
2.2 Stirred Versus Airlift Reactors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 123
2.3 Engineering Analysis of Reactors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 124
2.3.1 Mixing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 126
2.3.2 Mass Transfer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 132
2.3.3 Shear . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 135
2.4 Calculation of Time C o n s t a n t s in Plant-Cell Reactors . . . . . . . . . . . . . . . . . . . . . . 137
2.4.1 Oxygen C o n s u m p t i o n . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 139
2.4.2 Mass Transfer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 140
2.4.3 Mixing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 142
2.4.4 Shear . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 142
2.5 C o m p a r i s o n of Reactor Performance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144
2.6 I m p r o v e m e n t of Reactors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 145
2.6.1 Airlift Reactors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 145
2.6.2 Stirred T a n k s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 146
2.7 Alternative Reactor Designs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
3 Immobilised Plant Cells . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
3.1 M a s s Transfer Considerations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 148
3.1.1 External Mass-Transfer Resistance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 148
3.1.2 Internal Mass-Transfer Resistance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 149
3.2 Self-Immobilised Plant Cells . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 151
3.3 Immobilised-Cell Reactors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 151
3.4 Effect of I m m o b i l i s a t i o n on Plant-Cell Metabolism . . . . . . . . . . . . . . . . . . . . . . . . 153
4 Differentiated Plant Tissue . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 154
4.1 H a i r y Roots . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 154
4.2 E m b r y o s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 156
4.3 Shoots, Buds and Plantlets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 157
5 P h o t o a u t o t r o p h i c Plant Cells . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 158
6 Conclusions and F u t u r e O u t l o o k . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 160
7 Nomenclature . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 161
8 References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 163
1 Introduction
Cultured plant cells are capable of producing a range of valuable secondary
metabolites. Research into large-scale plant-tissue culture is aimed mainly at
commercial manufacture of these compounds, although reactor culture of orga-
nised tissue can also be used for plant propagation. The shikonin process developed
by Mitsui Petrochemical Industry Ltd in Japan demonstrated in 1983 that
industrial-scale culture of dispersed plant-cells is technologically feasible. A major
effort is still being made to define the requirements of suspended plant-cell cultures
and their responses to different reactor conditions. This has proven to be a difficult
and challenging area of engineering research. Familiar problems of reactor
hydrodynamics, mixing, and oxygen transfer are coupled with plant-cell shear-
sensitivity and aggregation, formation of non-Newtonian suspensions at relatively
low cell densities, and poor definition of conditions required for secondary
synthesis. Although reactors up to 20000 litres have been used to culture plant
cells, the task of designing an appropriate reactor system for secondary-metabolite
production has not been completed. Yields of most secondary compounds in vitro
are still very low and decline even further after scale-up. In most cases, economics
demands that high cell densities be used to raise volumetric productivity; a reactor
capable of supporting 100 kg m - a dry weight suspended plant-cells without mixing,
mass transfer or shear-damage problems is yet to be found. While immobilised
cells and differentiated plant tissue are rapidly replacing suspended plant-cells as
the focus of research, recent progress with enzyme regulation and genetic
manipulation of plant cells may eventually improve the commercial prospects of
large-scale suspended-cell culture.
Interest in immobilising plant cells followed from the theory that secondary
metabolites are more likely to be produced by slow-growing cells with high cell-cell
contact than by dispersed suspended cells. Early work showed that product yield
and excretion could be greatly improved by immobilisation, with the added benefit
of reduced shear damage to cells. This promoted several studies of immobilisation
methods and bioreactors suitable for large-scale operations. The natural tendency
of plant cells to clump together is being exploited; self-immobilised or aggregated
plant-cells can be grown in packed- and fluidised bed reactors. Surface-attachment
of plant cells in membrane and biofilm reactors is also being investigated. Provision
of adequate oxygen transfer and maintenance of aseptic conditions during reactor
operation are two important factors in immobilised plant-cell reactor design.
Once it was demonstrated that hairy-root cultures produce high levels of
secondary metabolites there was a need for appropriate reactors. Work in this
new area of bioreactor design is progressing rapidly. A major problem with
large-scale culture of roots is poor distribution of biomass. Novel systems are now
being developed to fix the roots in place so the vessel can be completely filled
with tissue. In addition, spray reactors are being used instead of submerged culture;
this mode of culture promotes oxygen transfer and provides low-shear conditions.
Another innovation is use of bioreactors to grow somatic embryos; this technology
could be applied for mass propagation of plants and production of artificial
seed as well as for large-scale phytochemical synthesis. Differentiated shoots and
l 18 P.M. Doran
plantlets have recently been cultivated in reactors. New ideas for providing light
to photoautotrophic plant cells in large-scale vessels are also now appearing.
The aim of this review is to identify problems associated with large-scale culture
of various forms of plant tissue and to summarise the most recent developments
in reactor technology. Emphasis is given to the engineering aspects of plant culture.
More general treatment of applied plant-tissue culture can be found in previous
articles by Fowler [1], Misawa [2], Kurz and Constabel [3], and Rhodes et al. [4].
2.1.1 Rheology
Plant culture rheology remains poorly characterised. The morphology of plant
cells in vitro is highly varied; near-spherical and elongated ceils co-exist in
suspension cultures, and there is usually a wide size-distribution. These properties
mean that high-density plant-cell suspensions become viscous, to the order of
several hundred centipoise. High viscosity leads to longer circulation times in
reactors, and formation of dead zones where the liquid remains stagnant. Reduced
turbulence in viscous suspensions has a pronounced effect on oxygen transfer as
air bubbles are not dispersed as easily as in non-viscous media.
Kato et al. [5] determined that Nicotiana tabacum cultures acted as non-
Newtonian fluids even at low cell concentrations of 0.9 kg m - 3. The culture filtrate
after removal of cells maintained a Newtonian character through the entire culture
period. As the cell concentration grew from 0.9 to 13 kg m-3 dry weight, filtrate
viscosity increased from 0.9 mPa s to only 2.2 mPa s even though apparent viscosity
of the whole cell suspension increased by a factor of 27.5. This small increase in
filtrate viscosity indicates that excretion of extracellular material by plant cells
does not contribute significantly to the overall viscosity; the viscous and non-
Newtonian character of the suspension is due mainly to the solids content.
Several models are available to describe non-Newtonian fluids; those applicable
to microbial cultures are discussed by Roels et al. [6]. Bingham and power-law
models have been used to describe plant-cell suspensions. The equations are:
where ~ is shear stress, % is the yield stress, K is the consistency index, j~ is shear
rate, and n is the flow behaviour index.
By analogy with Newton's law of viscosity, an apparent viscosity gapp can be
defined from Eq. (2) for power-law fluids:
g, pp = K? ~- ~ (3)
so that:
z = Gpp3' 9 (4)
Species n % Reference
(Pa)
2.1.2 A g g r e g a t e F o r m a t i o n
Aggregation is generally a poorly-controlled parameter in plant-cell suspensions;
aggregates and flocs range in size from two cells to clumps of 0.5 cm diameter or
more. Wagner and Vogelmann [9] observed that plant-cell cultures forming flocs
or pellets caused fewer macromixing problems and less shear damage than
suspensions containing mostly single cells. Tendency towards clumping and wall
growth is somewhat species dependent; reducing the Ca 2+ level in the medium
[14, 15] and addition of sorbitol and cell-wall-degrading enzymes such as pectinase
and cellulase [16] have been reported to reduce aggregate formation. In bioreactors,
the size of aggregates depends also on the shear levels present [17]. Control of
clumping may be achieved to some extent by changing reactor geometry
and operating conditions.
The range of aggregate size and specific gravity in several plant cultures is
shown in Table 2. These data can be compared with the density of the suspending
liquid. Murashige and Skoog medium containing 3% sucrose has a specific gravity
of about 1.013; this reduces to about 1.000 at the end of batch culture [10].
Design of Reactors for Plant Cells and Organs 121
Table 2. Size and specific gravity of plant-cell aggregates in suspension culture. (Adapted
from Tanaka [10])
of high apparent viscosity suffered more damage at constant shear rate than
suspensions with low apparent viscosity. Accordingly, shear stress rather than
shear rate appears to be the important variable in assessing shear sensitivity of
cells. When high local velocity-gradients are present, such as at the tip of the
impeller in stirred systems, the extent of cell damage may also depend on the
maximum rather than average shear stress. Maximum shear stresses withstood
by plant cells in laminar flow in a high-shear Couette rheometer have been reported
by Chen et al. [25]. For tobacco cells the critical shear level was 2.5 Pa and for
licorice cells 8 Pa. These figures indicate a significantly lower tolerance to shear
stress than that observed by Meijer [23] with plant cells subjected to stirrer speeds
of 250-1000 rpm.
The actual mechanism of shear damage to plant cells has not been studied.
However, hydrodynamic effects on animal cells have been analysed by several
groups [26-34]; the results may also be applicable to suspended plant-cells. A
complete picture of shear effects in stirred and aerated reactors has not yet been
determined unequivocally. In stirred microcarrier-cultures of animal cells where
there is no entrainment of gas bubbles, interactions between cells and turbulent
eddies are considered most likely to cause shear damage [26, 35]. If the size of the
eddies is approximately the size of the microcarriers, the eddies dissipate their
energy against the surface of the microcarriers and cause damage to the attached
cells. The dimension of the smallest eddies in turbulent flow is given by the
Kolmogorov scale:
n = (~) 1/4
(5)
where 1"1is eddy size, v is fluid kinematic viscosity, and a is local energy dissipation
rate per unit mass of liquid. In large-scale equipment with low-viscosity liquids,
rl is of the order 30-100 gm; smaller eddies are often produced in laboratory-scale
vessels. A relationship between cell damage in microcarrier cultures and T1 has
been demonstrated; decreasing the Kolmogorov scale below 1/2-2/3 the microcar-
rier diameter increases cell damage [26, 28]. Croughan et al. [28] also showed that
increased viscosity, which results in increased Kolmogorov scale, reduces damage.
Since plant cells have dimensions in the vicinity of 30-100 Ixm and can form
clumps of even greater size, it is likely that interaction with turbulent eddies also
causes shear damage in plant suspensions. To date, however, there has been no
experimental investigation of these effects.
In aerated systems, shear damage can occur at much lower impeller speeds than
those producing eddies of the dimensions described above. There is evidence that
in systems with bubble entrainment, significant shear damage occurs at the liquid
surface as a result of bubble disengagement [32]. Data of Handa-Corrigan et al.
[29] relating cell viability to bubble-column height and formation of a stable
protective foam support this theory. Investigations by Kunas and Papoutsakis
[30] have distinguished shear-damage mechanisms in reactors with and without
bubble entrainment. At low agitation speeds between 150 rpm and 600rpm,
damage to hybridoma cells was insignificant if there were no gas-liquid interface
Design of Reactors for Plant Cells and Organs 123
in the vessel. Damage at these stirrer speeds occurred only when bubbles were
allowed to burst at the liquid surface. At 800 rpm, which corresponded to a
Kolmogorov-eddy size comparable with the size of the cells, damage to the cells
increased but was still less significant in the absence of a gas interface.
Studies of shear damage in plant-cell suspensions have concentrated exclusively
on the effects of mechanical agitation. While different species can be expected to
have different shear sensitivities, the extremely wide variation in shear response
from complete cell destruction at 28 rpm [9] to tolerance of 1000 rpm [23] may
be due to intrusion of other factors not considered in these studies. The shear
stress produced by bursting bubbles may be one such factor.
t r
t- - ' l _ _
I I
cl riser
O0 0 O]
o [o -downcomer vdownoomer
oO
~176176
o
i o o o ~
-<
Air Air
Internal loop
External loop
with draft tube
vortexing; a wide variety of impeller sizes and shapes can be used to produce
different flow patterns inside the reactor. In tall vessels, installation of two or more
sets of impellers improves mixing.
In airlift reactors, air sparging is used for both mixing and aeration. The reactor
is divided into two distinct zones. Gas is sparged into the riser section; the increased
gas hold-up and decreased fluid density in this section cause the fluid to rise. Gas
disengages at the top of the reactor leaving heavier liquid which recirculates
through the downcomer. Low and relatively uniform shear fields, low power input
and extended aseptic operation are important advantages of airlift reactors for
plant-cell operations.
As shown in Fig. 1 there are two basic types of airlift reactor:
a) the internal-loop airlift in which riser and downcomer are separated by an
internal baffle or draft tube; and
b) the external- or outer-loop airlift where riser and downcomer are separate
tubes connected by short horizontal sections at top and bottom.
Both types of airlift usually have circular cross-sections. In external-loop reactors
the head-space region at the top of the vessel can be modified to change the
efficiency of gas-liquid separation; this determines whether bubbles are recirculated
or not and can affect gas hold-up, liquid velocity and turbulence. In internal-loop
vessels either the draft tube or annulus may be sparged. Changing the distance
between the lower edge of the baffle and the bottom of the reactor can affect the
pressure drop in this region with consequences for liquid recirculation velocity,
gas hold-up and mixing. Gas-liquid separation depends on the depth of submersion
of the upper edge of the draft tube in the liquid. In internal-loop vessels, since the
downcomer and riser connect directly at the top of the draft tube, gas-liquid
separation is generally not as effective as in external-loop devices and gas bubbles
are often drawn into the downcomer. The capabilities of internal- and external-loop
configurations are compared in Table 3 in terms of various performance indicators
[S0l.
Table 3. Relative performance of internal- and external-loop airlift reactors. (Adapted from
Chisti [50])
Parameter Reactor
Internal-loop ExternaMoop
For plant-cell suspensions it is crucial that the reactor provides good mixing
and suspension of solids, adequate oxygen transfer, and non-damaging levels of
hydrodynamic shear.
2.3.1 Mixing
Good mixing achieves homogeneous distribution of oxygen and other nutrients
in the bulk fluid; dead zones where anoxia or nutrient starvation can occur are
prevented. Mixing also reduces localised build-up of substances such as sugar
which can affect cell metabolism if allowed to accumulate in unmixed regions.
Both heat and mass transfer are greatly influenced by mixing and the intensity of
turbulence in the reactor.
The efficiency of mixing is usually quantified in terms of mixing time. Mixing
time (tmx) is defined as the time required to reach some arbitrary level of uniformity
in the liquid being mixed. In a recirculating-flow mixed reactor, the number of
circulations required for complete mixing depends on the degree of turbulent
diffusion. For a stirred tank with several baffles and small impeller, a satisfactory
relationship between circulation time and mixing time is [51]:
t ~ = 4tc (6)
where tc is the time taken for liquid to travel one complete circulation loop in the
reactor. For internal-loop airlift reactors where gas is sparged into the annulus:
(Ad) ~
tmx = 3.5to \At,/ (7)
(A)OS
tmx = 5.2to \ A r J (8)
10 7
10 6
10 5
10 4
or
r~ 10 3
10 2
IMPELLER TYPES
Fig. 2. Viscosity ranges for different impeller designs. (Adapted from Holland and Chapman
[53])
Metzner and Otto [18] determined that average shear rate (gay) in a stirred tank
is directly proportional to agitation speed:
where the value of k depends on impeller design and type of fluid, and N i is
rotational speed. This relationship has been tested for both Newtonian and
non-Newtonian fluids.
Flow visualisation studies with non-Newtonian fluids have shown that flow
patterns are different from those set up by the same equipment in Newtonian
liquids [54]. For pseudoplastic fluids, apparent viscosity in the impeller region is
rather low because of the higher shear rates; viscosity increases with distance from
the impeller. The recommended geometry for mixing operations is also different
with non-Newtonian fluids. For an ungassed low-viscosity fluid the optimum ratio
of tank diameter to impeller diameter (Dr:D0 is 2.5-3.0: 1. Reduction of this
ratio to 1.5 2.0:1 [55, 56] or use of a close-clearance anchor or helical impeller
gives better mixing with viscous non-Newtonian fluids. At low viscosities the
volume of fluid set in motion by the impeller is of the order of several cubic
impeller-diameters. In viscous fluids this volume decreases markedly. The thickness
of the fluid layer moved by the impeller may become only a fraction of the
128 P. M. Doran
impeller diameter, so that impellers of small diameter are not likely to provide
adequate mixing.
Mixing time in non-aerated Newtonian fluids in baffled cylindrical tanks has
been correlated with system geometry and operating parameters by Norwood and
Metzner [57]. Flat-blade turbine impellers were used. The dimensionless group,
fm:
tmx(NiD2)2/3 gl/6D~/2
fm = (10)
H1/21~3/2
L "L'T
Di2NiQL
Rei - (11)
In Eqs. (10) and (11) tmx is mixing time; Ni is impeller rotational speed; Di is
impeller diameter, g is gravitational acceleration; HL is liquid height; Da- is tank
diameter; QL is fluid density; and gL is fluid viscosity.
Norwood and Metzner [57] proposed that the relationship between fm and Rei
for Newtonian liquids could be used to calculate approximate mixing times for
pseudoplastic fluids at high Reynolds numbers. For non-Newtonian fluids and
viscous cultures of Streptomyces niveus, Wang and Fewkes [58] found that mixing
times measured in stirred tanks were in good qualitative agreement with the
Norwood-Metzner correlation, although the value of fm at Reynolds numbers
Design of Reactors for Plant Cells and Organs 129
above 104 was approx. 20. The discrepancy may be due to aeration effects; aeration
of suspensions increases mixing times for a given Reynolds number [59]. In addition,
Wang and Fewkes calculated the non-Newtonian impeller Reynolds number using
the equation proposed by Calderbank and Moo-Young [60]:
Rel D2 n0,( , )n
- OAK ~ 2 (12)
where K is the consistency index and n is the flow behaviour index for a power-law
fluid. This change should not, however, affect the correlation at high Re i. Ap-
plication of the Norwood-Metzner correlation to non-Newtonian cell suspensions
is discussed further by Charles [7].
Riser 0
o
,o
00
'o
0 ] I ]~
Downcomer e 0
- Riser
o0.0
g.?,'
0r o~
,oOOoOo -- 9
Downeomer
' ~ ~ oo~
"~2 0000200,0
t
GAS b
t
GAS
0 0 0 0 ~~
J~
0 o o
00 0 ~ o O~
'o 0 i
oo o
0 0 0 O0 9
i 0 00 o .
e 0 O0 O0 *
0 00 0 0 4
o~176176
II
II
GAS
t
e GAS d
Fig. 4a-d. Influence of sparger location on gas distribution in airlift reactors. Poor
distribution of gas in a) internal- and b) external-loop reactors. Proper positioning of the
sparger is shown in e) and d). Cross-hatched areas indicate zones which should be filled in
to improve liquid flow and prevent biomass settling. (Adapted from Chisti and Moo-Young
[63])
where u ~ is superficial liquid velocity in the riser, m s-1; uGr is superficial gas
velocity in the riser, m s - 1; Ad is downcomer cross-sectional area, m2; A r is riser
cross-sectional area, m2; and gapp is apparent viscosity, Pa s. Equation (13) was
developed using viscous pseudoplastic suspensions of 0.5-1.5 wt% carboxymethyl
cellulose (CMC) in 0.5 M Na2SOr in external-loop devices with Aa/Ar ratios of
0.111, 0.250 and 0.444. The liquid height in these reactors was 1.88 m. The apparent
viscosity was estimated using the correlation for shear rate in bubble columns
reported by Nishikawa et al. [67]:
~, = 5000u~ (14)
u~ _> 0.04 m s -1. From Eq. (14) the value of •app in Eq. (13) for pseudoplastic
fluids can be taken as:
where ~app is apparent viscosity, Pa s; uar is superficial gas velocity in the riser,
m s-1; n is the flow behaviour index; and K is the fluid consistency index, Pa sn.
When applied to Aspergillus niger suspensions containing up to 20 kg dry weight
per m 3 cells, Eq. (13) over-predicted ULr especially at higher Uar values [68].
Equation (13) was developed for a specified reactor height; however, since liquid
circulation relies on the difference in hydrostatic pressures between riser and
downcomer, liquid velocity also depends on the height of the reactor. Most reported
data on the effects of reactor height are derived from small equipment; these are
difficult to extrapolate because of the predominance of end effects. Rousseau and
Bu'Lock [61] determined for 0.68 < HL(m) _< 1.64 that:
t ~ oc HL~'7 (16)
where H E is liquid height in the reactor, m. This relationship was developed using
air-water in an internal-loop device with Dr/D c between 0.2 and 0.6.
Other equations for liquid velocity have been developed from theoretical analysis
of the hydrodynamics in external- and internal-loop airlift reactors [50]. An energy
balance for non-Newtonian liquids leads to the equation:
+ 2
E KT
( l ~ g , ) z + KB \ A J
]
(1 -- aa)2 = 0
(17)
where uLr is superficial liquid velocity in the riser, m s- 1; APFr is frictional pressure
drop in the riser, Pa; APFa is frictional pressure drop in the downcomer, Pa; EL
is density of the liquid, kg m 3; g is gravitational acceleration, m s-2; HD is
gas-liquid dispersion height, m; ~r is gas hold-up in the riser; ea is gas hold-up in
the downcomer; K T is frictional loss coefficient for the top of the reactor; KB is
frictional loss coefficient for the bottom of the reactor; Ar is riser cross-sectional
area, m/; and Ad is downcomer cross-sectional area, m 2. To use Eq. (17), gas
hold-up in the riser er must be known or estimated from independent correlations,
APFr and APFd are determined from shear stress calculations, and KB and KT are
estimated from correlations involving the area for flow between riser and
downcomer. Solution for uLr requires trial-and-error procedures.
Once uLr is determined, the superficial liquid velocity in the downcomer ULd
can be estimated from continuity:
For internal-loop airlift reactors and for external-loop reactors where the lengths
of horizontal top and bottom sections are small, circulation time can be determined
from ULr, ULe and geometric parameters:
Hr Hd
tc = - - + -- (19)
ULr ULd
where H r is height of the riser and H a is height of the downcomer. Mixing time
can then be calculated from Eq. (7) or (8).
kLaD
kLaL -- (21)
1 -- eT
viscosity increase, the rate of mass transfer can drop to as low as 15% of its initial
value [77]. Prediction and improvement of oxygen mass transfer in reactors relies
on the prediction and enhancement of kLaL.
2.3.2.1 StirredTanks
In aerated stirred reactors one of the functions of the impeller is to break up and
disperse the air bubbles. With pseudoplastic fluids, shear thinning near the impeller
causes gas to channel in toward the impeller rather than into the remainder of
the tank. This non-uniform distribution is difficult to prevent and has adverse
effects on both mass transfer and mixing. Less bubble break-up is observed with
viscous non-Newtonian fluids; kLae values tend to be considerably lower than in
Newtonian systems under the same operating conditions [78].
Correlations for kLaL values in stirred tanks are of two types. Correlations
based on power input are reviewed by Schiigerl [76] who also discusses effects of
gas flow rate, fluid properties and impeller speed on power requirements for
non-Newtonian fluids. Other kLae correlations are presented in terms of dimen-
sionless groups representing system geometry, fluid properties and operating
conditions.
Rates of oxygen desorption in pseudoplastic non-Newtonian fluids were
measured by Yagi and Yoshida [78] in an agitated tank with four baffles and
6-blade turbine impeller. The stirrer speeds tested were 5 s ~ to 10 s z; superficial
gas velocities ranged between 0.002 and 0.008 m s- 1. The correlation obtained for
kLaL in non-Newtonian fluids is:
kLaLD2--O'o6(D2NiQL)l"5(DN2)~176176176
\~app ,/ k,~i--J \ ~ L J
/ N D \o.32
x( i~] [1 + 2.0(ZN0~176 (22)
\uo/
where kLaL is the mass transfer coefficient, s-1; Di is impeller diameter, cm; ~L
is diffusivity in the liquid, cm 2 s - 1; Ni is rotational speed of the impeller, s- 1; ~L
is liquid density, g cm- 3; Papp is apparent viscosity, g cm- 1 s- 1; g is gravitational
acceleration, cm s-2; UG is superficial gas velocity, cm s-~; c~ is surface tension,
g s-Z; and )~ is characteristic material time, s. For non-Newtonian fluids % is a
constant defined as the reciprocal of the shear rate at which the ratio of apparent
viscosity to zero-shear viscosity is 0.67. The correlation of Eq. (22) is limited by
determination of )~; viscosity at zero-shear is difficult to measure accurately.
Apparent viscosity in Eq. (22) is calculated using the Metzner and Otto [18]
relationship for average shear rate as outlined in Eqs. (3) and (9).
An alternative to Eq. (22) is the correlation of Perez and Sandall [79]:
where gg is gas viscosity, g cm- 1 s- 1; the other symbols and units of Eq. (24) are
the same as for Eq. (22). Note that Eq. (23) is a dimensional equation. Apparent
viscosity in Eq. (23) is determined by the method of Calderbank and Moo-Young
[60]:
K (3n + 1~ ~
~app -- (11Ni) 1 - n \ ~ 4 n - n // 9 (24)
Equation (23) was developed for carbon dioxide absorption in a 0.25% Carbopol
solution using a tank equipped with four baffles and a 6-flat-blade turbine impeller
operated at 2810 _< Re i _< 26700; Re i is defined by Eq. (11) with the viscosity
given by Eq. (24). kLaL values calculated using Eq. (23) are for carbon dioxide
and must be corrected for oxygen transfer. The correction based on the two-film
theory of mass transfer is [50]:
where N is diffusivity. The ratio ~o2/Nco2 is equal to approx. 1.2. The applicability
of Eq. (23) is discussed further by Yagi and Yoshida [78].
Popovic and Robinson [84] have proposed a correlation for kta ~ in external-loop
airlift reactors:
where kLao is the mass transfer coefficient based on dispersed volume, s-1; uG~
is superficial gas velocity in the riser, m s- 1; @t is the diffusivity of oxygen in the
liquid, m 2 s- 1; ~i. is liquid density, kg m - 3; Ad is downcomer cross-sectional area,
m2; A~ is riser cross-sectional area, m2; btappis apparent viscosity, Pa s; and cr is
interfacial tension with respect to air, N m - 1. Reactor geometry used to develop
Eq. (27) is the same as that described for Eq. (13); apparent viscosity is calculated
using Eq. (15). Popovic and Robinson [84] tested Eq. (27) with a range of viscous
and non-Newtonian fluids; prediction of kLao was within 13 % of measured values
in each case.
Values of kLaD from Eq. (27) can be converted to kLat using Eq. (21). Gas
hold-up in airlift reactors is not uniform; eT can be evaluated using Eq. (28):
where ~d is gas hold-up in the downcomer and ar is gas hold-up in the riser.
For external-loop vessels ad can be taken as approximately 0.3at; ar is given by:
Equation (29) was developed by Popovic and Robinson [84] for viscous CMC
solutions, but has been shown by Allen and Robinson [68] to apply to AspergilIus
niger suspensions at 0.02 _< uar (m s- 1) < 0.2 and 0.01 < ~app (Pa s) < 0.5.
An alternative equation for kLaL in external-loop airlift reactors has been
proposed by Chisti et al. [83]. This correlation involves solids concentration
explicitly:
2.3.3 Shear
Mechanical damage in plant-cell suspensions depends on the intensity of shear
developed in the reactor. A limited number of correlations for hydrodynamic shear
in stirred and airlift reactors is available.
136 P.M. Doran
Equation (9) is widely accepted as a means for calculating average shear rates in
stirred vessels. Variation of shear has been investigated by Metzner and Taylor
[54] and Oldshue [85]; maximum shear levels occurring at the impeller tip are
significantly higher than those given by Eq. (9).
Wichterle et al. [86] measured maximum shear rates on the front surface of a
Rushton 6-blade turbine in non-Newtonian pseudoplastic fluids. For Re i > 10
the following relationship was determined:
9 / N ? - nD.2,, ",,1/(1+n)
7ma, (1 + 5.3n)1/"~ ' ~ '~L) (31)
Ni
where n is the flow behaviour index for power-law fluids and K is the consistency
index.
At the present time, reliable estimation of average shear levels in airlift reactors
is difficult. Correlations reported for pneumatically agitated columns vary conside-
rably, as summarised in Table 4. Of these relationships, the equation by Nishikawa
et al. [67] is used almost exclusively; it was developed by analysing heat-transfer
coefficients in Newtonian and non-Newtonian liquids in a bubble column at
0.04 _< u~ (m s -1) _< 0.1. The correlation of Henzler [87] was determined by
fitting literature data on mass transfer coefficients obtained with non-Newtonian
liquids. The broad lack of agreement between the shear-rate correlations of Table 4
causes difficulties in estimation of mass transfer coefficients and mixing times;
correlations such as Eqs. (13) and (27) depend on estimation of laavp which, for
non-Newtonian fluids, relies on knowledge of j,.
Constant but relatively low shear-stress is produced in air-driven reactors by
flow of fluid at the walls of the vessel. Tramper et al. [88] estimated the shear stress
associated with frictional wall losses in an airlift reactor using the friction-factor
equation:
"Cw =
1 2
~ ~LULkw
(32)
C 9
trx, - (33)
ro 2
Table 6. Literature values for observed specificoxygen uptake rate qo2 in suspended plant-cell
cultures
Plant %2 Reference
kg 0 2 (kg dry weight)-1 s-x x 10 6
qo2
kg 0 2 (kg fresh weight)- 1 s- 1 x 106
10000
trxn
1000
tmt
tmx \ Jinx
trxn _\
100
(s)
....................... ~_ tmt
10 I q -I-'~ ...... .J
~, trxn
1000
tmt
t rnx
~tmt
trxn
I0 0 "":k
(s)
';:i~......... tmx
".%... ...............................
0 5 10 15 20
b N i (s1 ) )'~
Fig. 5a, b. Mixing, mass transfer and oxygen consumption in bioreactors at a suspended
plant-cell concentration C s = 5 k g m -3 dry weight a) external-loop airlift reactor; b)
stirred-tank reactor
density. The calculated time constants for oxygen uptake are plotted in Figs. 5
and 6. For comparison, time constants for microbial oxygen consumption are of
the order 9-45 s [90].
10000 , , , ,
1000
tmt "",~"............................ ~,trnx
" , , , t rxnJk ................................
tmx
trxn 100
(s)
10 I I P P
l
\
1 1000
10 ~ .............
0 5 10 15 20
b N i (s"1)
Fig. 6a, b. Mixing, mass transfer and oxygen consumption in bioreactors at a suspended
plant-cell concentration Cs = 30kgm -3 dry weight, a) external-loop airlift reactor; b)
stirred-tank reactor
Since data on the viscoelastic properties of plant-cell suspensions are not available
Eq. (23) cannot be easily applied for calculation of kLaL. Instead, with the above
values for K and n, Eqs. (23)-(25) are used to evaluate kLaL for oxygen transfer
in stirred vessels. Results for the mass-transfer time constants are plotted in Figs.
5b and 6b. Depending on the stirrer speed between 0.2 s-1 and 20 s-1, g,vp values
are 88 m P a s to 23 mPa s at the beginning of culture and 850 mPa s to 220 m P a s
at the end.
2.4.3 Mixing
Mixing times for external-loop airlift reactors can be calculated from Eqs. (13),
(16), (18), (19) and (8). The results are plotted in Figs. 5a and 6a.
Mixing times for stirred-tank reactors are obtained from the Norwood-Metzner
correlation of Eq. (10). The value of fm is taken to be 20 as determined by Wang
and Fewkes [58] with Streptomyces niveus cultures. Under turbulent conditions
the mixing time in stirred reactors is largely independent of viscosity; consequently
t~,x values are the same at the end of the plant-cell culture as at the beginning.
The results are plotted in Figs. 5b and 6b.
2.4.4 Shear
10000 , , ,
09
j/j~176176
.1" ~
go 100
~~_-f~ cs=3%
10 I [
0.1 0.2 0.3 0.4 .5
& UGr (ms"1)
10000
c~ = 0.5% .......... :
1000
f - ~ .....3s= " ~
lOO
1 0 i i I r p i r i I I , r i i I r i I i
0 5 10 15
b N i (s"1)
Fig. 7a, b. Liquid shear stress produced in bioreactors containing suspended plant-cell
concentrations C~ = 5 kg m- 3 dry weight (0.5%) and Cs = 30 kg m- 3 (3%). a) airlift reactor;
b) stirred-tank reactor
Fluid shear-stress is maximum in airlift reactors in the top and bottom sections
between riser and downcomer where flow directions change abruptly. The
magnitude of this stress is calculated from Eq. (32). Results for the beginning and
end of plant-cell culture are plotted in Fig. 7a.
In the above analysis, the various functions of bioreactors: mixing, mass transfer
and provision of shear, were considered separately to determine which limits
growth under different culture conditions. Once the limiting process is identified,
measures can be taken to improve it.
and 6.0 x 10- 3 s 1. In the same study an anchor impeller operated at a maximum
speed of 40 rpm produced significantly lower product levels.
Hong et al. [38] installed 25 mm-wide teflon ribbons between the blades of two
turbine impellers in a 7-1itre stirred bioreactor in an attempt to reduce shear and
improve mixing. The modified agitator was tested for culture of strawberry cells;
growth did not occur and the viability of the cells dropped from 80% to less than
25% in 4 days. Better results were obtained by Hooker et al. [19] with a large
flat-blade impeller with and without sail-cloth blades. Higher maximum growth
rates were obtained compared with a smaller flat-blade turbine. Power consump-
tion by large sail-cloth impellers is considerably greater than with Rushton turbines
[511.
alginate beads are added to the medium; Prenosil et al. [119] have postulated that
alginate, being similar to the polycarboxylic acids in cell walls, mediates cell-cell
interactions in the presence of calcium ions to stimulate secondary synthesis.
Enhanced secondary product synthesis has also been reported for plant cells
immobilised in reticulated polyurethane foam [106, 120, 121].
For most efficient operation, the product of interest should be excreted from
immobilised cells into the medium so that the biomass does not have to be
periodically destroyed. There is evidence that immobilisation enhances release of
secondary products from plant cells [101,122]. Alternatively, release of compounds
can be induced by permeabilising agents such as DMSO [123], although limited
success has been achieved using this approach. Once in the medium, products can
be extracted in situ using adsorptive resins such as Amberlite XAD-7 [124] or
immiscible solvents such as vegetable oil [125] or n-hexadecane [126]. Continuous
product extraction encourages further synthesis and minimises degradative meta-
bolism.
Most techniques for immobilising plant cells produce an environment in which
there are substrate, oxygen, hormone and other concentration gradients. This
applies even to methods involving surface attachment; once a monolayer of plant
cells is laid down, other cells readily adhere to it to form a thick biofilm. Designing
large-scale immobilised-cell processes requires knowledge of growth and produc-
tion kinetics. To determine whether mass transfer limits plant-cell activity, analysis
of diffusion processes must be carried out.
Equations for diffusion and reaction are applied to immobilised cells for prediction
of effective reaction rates. Intraparticle diffusion is the main resistance to mass
transfer in entrapped-cell systems, although external boundary-layers surrounding
the particle can also contribute. Despite the relatively slow metabolic rate of plant
cells, effects of diffusion on activity of immobilised plant ceils cannot generally be
neglected.
As well as reducing the overall reaction rate, diffusional limitations mask intrinsic
kinetics. Care must be exercised when comparing the performance of different
immobilised-cell systems; apparently higher rates of growth or product synthesis
in one system may reflect an absence of diffusional effects rather than any intrinsic
metabolic response.
mixing increases the observed reaction velocity [127]. On the other hand, if the
rate of external diffusion is much more rapid than the rate of reaction, the observed
rate should be largely independent of mixing intensity.
Since these parameters vary considerably from system to system, results from one
analysis cannot be readily applied to other immobilised-cell situations.
The molecular species most likely to cause diffusional limitations in immobilised
plant-cells must be determined before the analysis can proceed. Mavituna et al.
[129] analysed diffusion of glucose into reticulated foam containing Capsicum
frutescens cells and determined that mass transfer of this substrate was not an
important problem. Furusaki et al. [130] measured rates of sucrose, oxygen and
codeinone consumption by immobilised Papaver somniferum cells; while effective-
ness factors for these substrates were close to unity, calculations showed that the
concentration of oxygen decreased considerably within the particles. Pras et al.
[131] concluded that limitations in oxygen transport were most likely after
measuring the effective diffusion coefficients and cellular demands for oxygen and
L-DOPA precursors by Mucuna pruriens entrapped in alginate gel.
In the case of oxygen as limiting substrate, for a spherical particle containing
plant cells, if the following assumptions can be made:
a) the cells are identical and their distribution in the particle is uniform;
b) oxygen transport in the particle occurs only by diffusion;
c) oxygen gradients occur only in the radial direction;
d) the effective diffusion coefficient for oxygen is independent of oxygen concentra-
tion;
e) respiration can be described using a Michaelis-Menten kinetic model; and
0 the number of plant cells is constant;
150 P.M. Doran
(1 d (r 2 dC))_ VmaxC
(35)
In Eq. (35), r is radial distance measured from the centre of the particle, m; C is
oxygen concentration within the particle, mol m-3; ~e is the effective diffusion
coefficient for oxygen in the particle containing cells, m 2 s- 1; Vmaxis the maximum
reaction velocity, mo! m - 3 s- 1; Km is the Michaelis-Menten constant, mol m - 3;
Ce is the oxygen concentration at the surface of the particle, mol m - 3 and R is the
radius of the particle, m. Equation (35) can be solved only by numerical integration
to give the oxygen profile in the particle, C as a function of r. A similar though
simpler equation is obtained for a biofilm which can be assumed infinitely long
and in which diffusion of oxygen occurs only in one direction.
Hulst et al. [128] calculated using Eq. (35) the particle size at which the oxygen
concentration in the centre is zero. The analysis was performed for cell aggregates
containing 100% Tagetespatula cells with a dry-weight cell-density of 10 kg m -3,
and intrinsic Michaelis-Menten constants vmax = 4 x 10 .3 m o l m -3 s -1 and
K m = 3 . 5 x l 0 - 2 m o l m -3. ~ e = 1 . 9 x 1 0 - 9 m - Z s -1 and C e = 0 . 2 5 m o l m -3.
Under these conditions, the critical particle diameter for depletion of oxygen at the
centre was 3 mm.
Some restrictions apply to use of the diffusion-reaction analysis represented
by Eq. (35). Equation (35) is a valid description of mass transfer and reaction
under conditions of constant biomass concentration in the immobilised-cell
particle. However, for systems in which growth occurs, the value of Vmaxvaries
according to the quantity of cells present. When cell concentration increases in
alginate beads, reticulated foam and membrane units, oxygen concentration
gradients and diffusional limitations will become more severe with time. In
addition, use of Eq. (35) assumes knowledge of the intrinsic kinetic parameters
for immobilised cells, Vmaxand K m. Experimental determination of these values is
not straightforward [128, 131] and use of suspended-cell parameters may not be
valid.
Although observed rates of oxygen consumption are often lower than in
suspension due to diffusion limitations, immobilised plant-cells still perform
secondary synthesis. The exact response to oxygen starvation depends on the
particular cell system and product of interest. Mavituna et al. [125] report that
capsaicin synthesis by foam-entrapped Capsicumfrutescens is enhanced when the
bulk-liquid dissolved-oxygen concentration is almost zero; if aeration is restored
capsaicin in the medium disappears. Pras et al. [131] showed that production of
Design of Reactors for Plant Cells and Organs 151
tion. The liquid-impelled loop reactor has been tested with free and immobilised
cells of Rubia tinctorum and Morinda citrifolia, with hexadecane as solvent
[145].
Although in vitro culture of excised plant roots was first demonstrated as early
as 1934 [157], slow growth of most roots in liquid medium and the requirement
for exogenous growth factors are major drawbacks for commercial application.
Growth rates can be increased significantly by genetic transformation with
Agrobacteriumrhizogenesbacteria to form 'hairy roots'; these can then be cultured
axenically without external hormones. Hairy-root cultures have been tested for
production of a broad range of secondary products; members of the Solanaceae
family known to synthesise alkaloids in their roots have been studied extensively.
Review articles are available discussing the plant host-range for Agrobacterium
rhizogenes, procedures for infection, methods of hairy-root culture, and the range
of compounds produced in vitro [158-161].
In an effort to identify appropriate reactor configurations, several studies have
been carried out with hairy roots in batch [162, 163] and continuous-flow
air-sparged reactors [164, 165]. Final biomass densities of 10-11 kg m - 3 dry weight
have been obtained [162, 164]. When continuous liquid-flow operation is used,
this figure can be increased to about 40 kg m -3 dry weight [165]. An important
problem with reactor-culture of hairy roots.is poor distribution of biomass in the
vessels. Growth of hairy roots, whether in shake-flasks or in reactors, usually
results in formation of a dense root ball [163, 166]. Nutrient and oxygen limitations
Design of Reactors for Plant Cells and Organs 155
in the centre of the ball reduce cell activity and can lead to necrosis. Poor
distribution of roots also means that the total volume of the reactor is not used
efficiently. To overcome these problems, roots must be more evenly spread and
mass transfer improved. Use of conventional stirred reactors is generally not an
option for freely-suspended hairy roots; root tissue is easily damaged by mechanical
forces and excessive shear in reactors causes callus generation and loss of
productivity [165]. Several modified and novel designs have been proposed recently.
A common feature of these reactors is that the roots are immobilised on supports
within the vessel.
Jung and Tepfer [155] cultured Calystegia sepium hairy roots in a 30-1itre stirred
vessel containing a stainless-steel basket for attachment of the roots. Taya et al.
[162] report culture of horseradish hairy-roots immobilised on polyurethane foam
in a column reactor. The foam stood vertically in the reactor and air was pumped
into the vessel at the bottom. As well as submerged culture, trickling-column
operation and periodic filling and withdrawal of medium were tested. Growth of
carrot hairy-roots and the relationship between growth rate, oxygen mass-transfer
coefficient and reactor design have been investigated using three different bio-
reactors: air-driven, rotating-drum, and modified stirred-tank [167]. Growth in
the air-driven vessel was slow; mass-transfer coefficients were significantly reduced
in the presence of 10 kg m - 3 roots. Roots in the rotating drum were damaged
until a polyurethane-foam sheet was attached to the inner wall of the rotating
drum to provide protection and support. When 10 kg m -3 roots were immobilised
on the foam, mass-transfer coefficients increased significantly compared with those
measured in the absence of cells. In the stirred-tank reactor, a stainless-steel mesh
was installed to separate roots from the turbine impeller; mass transfer in this
vessel was high and unaffected by root growth. In most reactors, oxygen supply
to immobilised hairy-roots depends not only on bulk gas-liquid mass-transfer but
also on diffusion into the tangled root mass.
A stirred reactor in which roots were isolated from the impeller by wire mesh
was also tested by Hilton and Rhodes [168]. In this case a cylindrical mesh cage
was used as support matrix for Datura stramonium hairy roots; this allowed a
more even distribution of biomass up the length of the vessel. With continuous
feeding of medium, packing densities reached 70% (drained weight per volume);
however, mixing under these conditions was not ideal so that localised regions of
low oxygen and nutrient levels are likely to have been formed.
One solution to the problem of oxygen transfer to hairy roots is use of mist
reactors. In these vessels roots are not submerged in nutrient medium but are
sprayed with fine droplets; the liquid drains through the biomass and is collected
for recycle. Hairy roots respond well to direct exposure to air; Wilson et al. [169]
have compared performance of Datura strarnoniurn hairy roots in droplet- and
submerged-culture reactors. A 7-d lag phase observed in submerged culture was
eliminated in the mist reactor and final specific hyoscyamine levels were slightly
higher.
In addition to biomass distribution and mass transfer there are other considera-
tions affecting feasibility of commercial-scale hairy-root culture. These have been
identified by Wilson et al. [169] as:
156 P.M. Doran
4.2 Embryos
Mass culture of somatic embryos can be used for either plant propagation or
production of phytochemicals.
Somatic embryogenesis is a proven means for producing large numbers of plants
in vitro. Suspension cultures in which embryos float freely in the medium are
especially amenable to large-scale mechanical handling. Somatic embryos of some
plants do not respond well in liquid-culture; nevertheless, propagation using
embryo suspensions has been demonstrated for vegetables such as celery and
carrot and for several other species which readily form embryos in a liquid medium
[171-173], The ability to produce and grow embryos on a large scale is necessary
for commercial application of this technology. Theoretically, large bioreactors
would not be necessary since each embryo develops into an entire plant and more
than 100000 embryos can be cultivated in a few litres of medium [174, 175]. In
practice, however, embryo cultures contain a heterogeneous mixture of cells, cell
clusters, embryos in various states of development and deformed embryos. Luckner
and Diettrich [176] report that many Digitalis lanata embryos at the bipolar,
heart-shaped or torpedo stage developed into plantlets with roots and leaves, but
fewer than 1% developed into normal plants after further cultivation. Procedures
for identifying viable embryos and separating them from the remainder of the
culture are required.
Another use of embryos is direct sowing in the field; each embryo is a potential
seedling and can be used as a seed substitute. Encapsulation of somatic embryos
in gel at a saitabie stage of development prevents premature desiccation and
provides growth regulators which control seedling development. The benefits of
artificial seeds include eiimination of seed-borne diseases and development of
improved strains in a relatively short time. Plant Genetics, Inc. in California has
developed techniques for artificial seeds of celery and alfalfa [177]. Somatic embryos
may also prove useful for long-term storage such as in germplasm banks.
Secondary metabolites are synthesised by somatic embryos; it is possible that
high levels are produced because of the effects of cellular organisation. When
torpedo-stage and older embryos produced in celery suspension cultures were
analysed for production of flavour compounds, higher levels of secondary products
Design of Reactors for Plant Cells and Organs 157
Like embryo culture, large-scale plantlet culture can be used for either mass
propagation or secondary-metabolite production. Micropropagation of many
plant species by shoot-culture techniques on agar is well established; however
organs can also be formed from suspended plant-cells. Requirements for shake-flask
and bioreactor cultivation of differentiated plants are being investigated.
Conventional shoot propagation is labour intensive, expensive, and requires
thousands of containers to produce a large number of plants. Use of submerged
158 P.M. Doran
and their characteristics can be found in the reviews of Hfisemann [188], Hfisemann
et al. [189] and Widholm [190]. Potential application of green hairy-roots of Tagetes
and Bidens species for secondary metabolite production is discussed by Flores
et al. [191].
In the laboratory, plant cells can be induced to develop photoautotrophic
metabolism in two-tier culture flasks containing K2CO 3 - K H C O 3 buffer which
produces a CO 2 atmosphere [192]. Formation of chlorophyll and fully functional
chloroplasts in cultured plant ceils has been shown [193] to involve:
a) lowering the sugar content and simultaneously increasing the CO 2 partial-
pressure above ambient level;
b) preventing oxygen and ethylene accumulation;
c) providing high light intensities (6000 8000 lux); and
d) maintaining a balanced ratio of auxin/cytokinin in the medium.
When photoautotrophism is achieved, it is usually in sugar-free medium but in the
presence of air containing 1 2% v/v CO2. The reason why most cultured plant-
cells are only capable of sustained photoautotrophic growth at highly elevated
CO 2 partial pressures is still unknown; however, sustained photoautotrophic
growth under ambient CO2 concentrations has been demonstrated for selected cell
lines of Arachis hypogaea, Daucus carota [194] and Gossypiurn hirsutum [190].
Since some enzymes involved in secondary synthesis are associated with
chloroplasts and chloroplast membranes, it is expected that production rates of
these compounds will be enhanced if chloroplasts are present. Comparative
measurements of secondary-metabolite accumulation in heterotrophic, photo-
mixotrophic and photoautotrophic cell have shown that green cell-cultures are
able to express products partially or totally synthesised by chloroplast-localised
enzymes. Systems in which greening enhances secondary production include
Lupinus polyphyllus [195], Morinda lucida [196], Solanum dulcamara [197] and
Solanurn laciniatum [198]. In many cases, however, transition from heterotrophic
to photoautotrophic conditions is not sufficient to induce secondary metabolism
characteristics of leaves, and production of commercially-interesting secondary
metabolites is either unchanged or decreased in photoautotrophic cells. For
example, in Nicotiana tabacum, nicotine and other secondary metabolites occur
mainly in heterotrophic rather than photoautotrophic cells, although N-methyl-
nicotinic acid formation is totally restricted to photoautotrophic cultures [199, 200].
Hagimori et al. [201] found that digitoxin content in undifferentiated photoauto-
trophic Digitalis purpurea cells was the same as in photomixotrophic cultures; in
shoots, photoautotrophism reduced the alkaloid content to 10% of that in
photomixotrophic cultures. Photoautotrophic suspensions of Catharanthus roseus
synthesise very low levels of vindoline and dimeric alkaloids [202]; production of
these compounds is drastically increased when heterotrophic conditions are
restored [189].
For those compounds produced in quantity in photoautotrophic cells, large-scale
production depends on design of suitable reactors. Important considerations for
maintenance of photoautotrophic cultures include provision of elevated levels of
CO2, reduced levels of Oz, and adequate light intensity.
160 P.M. Doran
7 Nomenclature
Ad downcomer cross-sectional area
aD interfacial area per unit dispersed volume
aL interracial area per unit volume unaerated liquid
Ar riser cross-sectional area
C oxygen concentration within particle
C* oxygen concentration in liquid in equilibrium with air
Ce oxygen concentration at surface of particle
CL oxygen concentration in liquid
C~ solids concentration
D diffusivity
Db bubble diameter
Dc column diameter
Dd downcomer diameter
De effective diffusion coefficient
Di impeller diameter
DE diffusivity in liquid
Dr riser diameter
Dr tank diameter
fm dimensionless mixing group
162 P.M. Doran
g gravitational acceleration
HD gas-liquid dispersion height
Hd height of the downcomer
HL liquid height
H~ height of the riser
K power-law consistency index
k empirical constant in Eq. (9)
KB frictional loss coefficient for bottom of airlift reactor
kL liquid-film mass transfer coefficient
Km Michaelis-Menten constant
Kx frictional loss coefficient for top of airlift reactor
kw resistance coefficient in Eq. (32)
n power-law flow behaviour index
Ni impeller rotational speed
OTR oxygen transfer rate
P power input
APvd frictional pressure drop in the downcomer
APF~ frictional pressure drop in the riser
qo2 observed specific oxygen consumption rate
R radius of particle
r radial distance measured from centre of particle
Re~ impeller Reynolds number
ro 2 volumetric oxygen consumption rate
tc circulation time
trot mass transfer time constant
tmx mixing time constant
trxn oxygen consumption time constant
UB bubble rise velocity
UG superficial gas velocity
UGr superficial gas velocity in the riser
HL liquid velocity
ULd superficial liquid velocity in the downcomer
ULr superficial liquid velocity in the riser
VL liquid volume
Vmax maximum reaction velocity
X cell concentration
exponent in Eq. (26)
shear rate
~{~v average shear rate
~/max maximum shear rate
local energy dissipation rate
8d gas hold-up in the downcomer
~r gas hold-up in the riser
~T total gas hold-up
q eddy size
L characteristic material time
Design of Reactors for Plant Cells and Organs 163
~[app a p p a r e n t viscosity
gg gas viscosity
gL liquid viscosity
v fluid kinematic viscosity
QI~ liquid density
(y surface tension
z shear stress
to yield stress
~w wall shear stress
8 References
1. Fowler MW (1983) In: Mantell SH, Smith H (eds) Plant biotechnology. Cambridge
University Press, Cambridge, p 3
2. Misawa M (1985) Adv Biochem Eng/Biotechnol 31:59
3. Kurz WGW, Constabel F (1985) CRC Crit Rev Biotechnol 2:105
4. Rhodes MJC, Robins RJ, Hamill J and Parr AJ (1986) NZ J Technol 2:59
5, Kato A, Kawazoe S and Soh Y (1978) J Ferment Technol 56:224
6. Roels JA, van den Berg J and Voncken RM (1974) Biotechnol Bioeng 16:181
7. Charles M (1978) Adv Biochem Eng 8:1
8. Metz B, Kossen NWF and van Suijdam JC (1979) Adv Biochem Eng 11:103
9. Wagner F, Vogelmann H (1977) In: Barz W, Reinhard E, Zenk MH (eds) Plant tissue
culture and its bio-technological applications. Springer, Berlin Heidelberg New York,
p 245
10. Tanaka H (1982) Biotechnol Bioeng 24:425
11. Scragg AH, Allan EJ, Bond PA and Smart NJ (1986) In: Morris P, Scragg AH, Stafford
A and Fowler MW (eds) Secondary metabolism in plant cell cultures. Cambridge
University Press, Cambridge, p 178
12. Vogelmann H (1981) In: Moo-Young M, Robinson CW and Vezina C (eds) Advances
in biotechnology, vol 1. Pergamon, NewYork, p 117
13. Laine J and Kuoppam~iki R (1979) Ind Eng Chem Process Des Dev 18:501
14. Takayama S, Misawa M, Ko K and Misato T (1977) Physiol Plant 41:313
15. Sahai OP and Shuler ML (1982) Can J Bot 60:692
16. Street HE (1973) In: Milborrow BV (ed) Biosynthesis and its control in plants. Academic,
London, p 95
17. Tanaka H, Semba H, Jitsufuchi T and Harada H (1988) Biotechnol Lett 10:485
18. Metzner AB and Otto RE (1957) AIChE J 3:3
19. Hooker BS, Lee JM and An G (1990) Biotechnol Bioeng 35:296
20. Scragg AH, Allan EJ and Leckie F (1988) Enzyme Microb Technol 10:361
21. Hooker BS, Lee JM and An G (1989) Enzyme Microb Technol 11:484
22. Allan EJ, Scragg AH and Pugh K (1988) J Plant Physiol 132:176
23. Meijer J (1990) PhD Thesis, Delft University of Technology, The Netherlands
24. Midler M and Finn RK (1966) Biotechnol Bioeng 8:71
25. Chen H, Wang J-J and Liu Y-G (1990) In: Abstracts VIIth Internat Congr Plant Tiss
Cell Cult, Amsterdam, 24-29 June 1990, p 341
26. Cherry RS and Papoutsakis ET (1986) Bioprocess Eng 1:29
27. Tramper J, Williams JB, Joustra D and Vlak JM (1986) Enzyme Microb Technol 8:33
28. Croughan MS, Sayre ES and Wang DIC (1989) Biotechnol Bioeng 33:862
29. Handa-Corrigan A, Emery AN and Spier RE (1989) Enzyme Microb Technol
11:230
30. Kunas KT and Papoutsakis ET (1990) Biotechnol Bioeng 36:476
31. Murhammer DW and Goochee CF (1990) Biotechnol Prog 6:391
164 P.M. Doran
111. Brodelius P (1986) In: Evans DA, Sharp WR and Ammirato PV (eds) Handbook of
plant cell culture, vol 4. Macmillan, New York, p 287
112. Novais JM (1988) In: Pais MSS, Mavituna F and Novais JM (eds) Plant cell
biotechnology. NATO ASI Set H18, Springer, Berlin, Heidelberg, New York, p 353
113. Rosevear A (1988) In: King RD and Cheetham PSJ (eds) Food biotechnology, vol 2.
Elsevier, London, p 83
114. Hulst AC and Tramper J (1989) Enzyme Microb Technol 11:546
115. Brodelius P and Nilsson K (1980) FEBS Lett 122:312
116. Ketel DH, Hulst AC, Gruppen H, Breteler H and Tramper J (1987) Enzyme Microb
Technol 9:303
117. Hulst AC, Tramper J, Brodelius P, Eijkenboom LJC and Luyben KChAM (1985) J
Chem Tech Biotechnol 35B: 198
118. Ayabe S-I, Iida K and Furuya T (1986) Phytochem 25:2803
119. Prenosil JE, Hegglin M, Baumann TW, Frischknecht PM, Kappeler AW, Brodelius P
and Haldimann D (1987) Enzyme Microb Technol 9:450
120. Holden MA, Hall RD, Lindsey K and Yeoman MM (1987) In: Webb C and Mavituna
F (eds) Plant and animal cells: process possibilities. Ellis Horw0od, Chichester, p 45
121. Ishida BK (1988) Plant Cell Rep 7:270
122. Barnabas NJ and David SB (1988) Biotechnol Lett 10:593
123. Brodelius P and Nilsson K (1983) Eur J Appl Microbiol Biotechnol 17:275
124. Robins RJ and Rhodes MJC (1986) Appl Microbiol Biotechnol 24:35
125. Mavituna F, Wilkinson AK, Williams PD and Park JM (1987) In: Moody GW and
Baker PB (eds) Bioreactors and biotransformations. Elsevier, London, p 26
126. Kim DJ and Chang HN (1990) Biotechnol Bioeng 36:460
127. Ford JR, Lambert AH, Cohen W and Chambers RP (1972) Biotechnol Bioeng Syrup
3: 267
128. Hulst AC, Meyer MMT, Breteler H and Tramper J (1989) Appl Microbiot Biotechnol
30:18
129. Mavituna F, Park JM, Wilkinson AK and Williams PD (1987) In: Webb C and
Mavituna F (eds) Plant and animal cells: process possibilities. Ellis Horwood,
Chichester, p 92
130. Furusaki S, Nozawa T, Isohara T and Furuya T (1988) Appl Microbiol Biotechnol
29:437
131. Pras N, Hesselink PGM, Guikema WM and Malingr6 TM (1989) Biotechnol Bioeng
33:1461
132. Kargi F, Ganapathi B and Maricic K (1990) Biotechnol Prog 6:243
133. Fuller KW (1984) Chem Ind 3 Dec, p 825
134. Dainty AL, Goulding KH, Robinson PK, Simpkins I and Trevan MD (1985) Trends
in Biotechnol 3:59
135. Hoekstra SS, Harkes PAA, Verpoorte R and Libbenga KR (1990) Plant Cell Rep 8: 571
136. Tsoulpha P and Doran PM (1991) J Biotechnol 19:99
137. Fuller KW and Bartlett DJ (1982) PCT Patent Application PCT/GB81/00186
138. Panda AK, Mishra S, Bisaria VS and Bhojwani SS (1989) Enzyme Microb Technol
11:386
139. Kobayashi Y, Fukui H and Tabata M (1987) Plant Cell Rep 6:185
140. Kobayashi Y, Fukui H and Tabata M (1988) Plant Cell Rep 7:249
141. Archambault J, Volesky B and Kurz WGW (1989) Biotechnol Bioeng 33:293
142. Archambault J, Volesky B and Kurz WGW (1990) Biotechnol Bioeng 35:702
143. Kurz WGW, Paiva NL and Tyler RT (1990) In: Progress in plant cellular and molecular
biology. Proc VIIth Internat Congr Plant Tiss Cell Cult, Amsterdam, 24-29 June 1990,
Kluwer Academic, Dordrecht, p 682
144. Tramper J, Wolters I and Verlaan P (1987) In: Laane C, Tramper J and Lilly MD
(eds) Biocatalysis in organic media. Elsevier, Amsterdam, p 311
145. Buitelaar RM, Susaeta I and Tramper J (1990) In: Progress in plant cellular and
molecular biology. Proc VIIth Internat. Congr Plant Tiss Cell Cult, Amsterdam, 24-29
June 1990, Kluwer Academic, Dordrecht, p 694
Design of Reactors for Plant Cells and Organs 167
186. Hagimori M, Matsumoto T and Mikami Y (1984) Agric Biol Chem 48:965
187. Hfisemann W, Amino S, Fischer K, Herzbeck H and Callis R (1990) In: Progress in
plant cellular and molecular biology. Proc VIIth Internat Congr Plant Tiss Cell Cult,
Amsterdam, 24-29 June 1990, Kluwer Academic, Dordrecht, p 373
188. Hiisemann W (1985) In: Vasil IK (ed) Cell culture and somatic cell genetics of plants,
vol 2. Academic, Orlando, p 213
189. Hfisemann W, Fischer K, Mittelbach I, Hfibner S, Richter G and Barz W (1989) In:
Kurz WGW (ed) Primary and secondary metabolism of plant cell cultures II. Springer,
Berlin, Heidelberg, New York, p 35
190. Widbolm JM (1989) In: Kurz WGW (ed) Primary and secondary metabolism of plant
cell cultures II. Springer, Berlin, Heidelberg, New York, p 3
191. Flores HE, Hoy MW and Pickard JJ (1987) Trends in biotechnol. 5:64
192. Hfisemann W and Barz W (1977) Physiol Plant 40:77
193. H/isemann W (1988) In: Pais MSS, Mavituna F and Novais JM (eds) Plant cell
biotechnology. NATO ASI Ser H18, Springer, Berlin, Heidelberg, New York, p 179
194. Bender L, Kumar A and Neumann K-H (1985) In: Neumann K-H, Barz W and
Reinhard E (eds) Primary and secondary metabolism of plant cell cultures. Springer,
Berlin, Heidelberg, New York, p 24
195. Wink M and Hartmann T (1980) Planta Med 40:149
196. Igbavboa U, Sieweke H-J, Leistner E, R6wer I, Hfisemann W and Barz W (1985)
Planta 166:537
197. Emke A and Eilert U (1986) Plant Cell Rep 5:31
198. Conner AJ (1987) Phytochem 26:2749
199. Ikemeyer D and Barz W (1989) Plant Cell Rep 8:479
200. Barz W, Ikemeyer D and Beimen A (1989) Planta Med 55:594
201. Hagimori M, Matsumoto T and Mikami Y (1984) Plant Cell Physiol 25:1099
202. Tyler RT, Kurz WGW and Panchuk BD (1986) Plant Cell Rep 3:195
203. Dalton CC (1980) J Exp Bot 31:791
204. Dalton CC and Peel E (1983) Prog Ind Microbiol 17:109
205. Yamada Y, Imaizumi K, Sato F and Yasuda T (1981) Plant Cell Physiol 22:917
206. Peel E (1982) Plant Sci Lett 24:147
207. H/isemann W (1982) Protoplasma 113:214
208. H/isemann W (1983) Plant Cell Rep 2:59
209. Fischer U and Alfermann AW (1989) Planta Med 55:686
210. Ohta Y and Takata T (1990) In: Abstracts VIIth Internat Congr Plant Tiss Cell Cult,
Amsterdam, 24 29 June 1990, p 341
211. Schumpe A and Deckwer W-D (1987) Bioprocess Eng 2:79
212. E1-Temtamy SA, Khalil SA, Nour-E1-Din AA and Gaber A (1984) Appl Microbiol
Biotechnol 19:376
213. Henzler H-J and Kauling J (1985) In: Proc 5th Eur Conf on Mixing, Wiirzburg,
Germany, 10 12 June 1985, p 303
Expert Systems in Bioprocess Control:
Requisite Features
K o n s t a n t i n B. K o n s t a n t i n o v 1., R o b e r t A a r t s 2 a n d T o s h i o m i Y o s h i d a 1
1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 170
2 Basic Structures, Functions and Implementation Schemes . . . . . . . . . . . . . . . . . . . . . . 170
2.1 Structures of Systems for Expert Control . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 170
2.2 Functions of the Knowledge-Based Module of the Control System . . . . . . . . . . . 172
2.3 Implementation Schemes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 173
3 Requisite Features of Expert System Development Tools . . . . . . . . . . . . . . . . . . . . . . 175
3.1 General Features . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 175
3.1.1 Real-Time Capabilities . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 175
3.1.2 Temporal Reasoning . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 176
3.1.3 Integration with External Software Modules . . . . . . . . . . . . . . . . . . . . . . . . . 180
3.1.4 Mechanisms for Knowledge Structuring . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 181
3.1.5 Handling of Various Categories and Levels of Knowledge . . . . . . . . . . . . . 182
3.1.6 Efficient Knowledge Debugging and Integrating . . . . . . . . . . . . . . . . . . . . . . 184
3.1.7 Answering Users' Questions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 184
3.1.8 Advanced User Interface with Graphical Capabilities . . . . . . . . . . . . . . . . . 185
3.1.9 Other Features . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 186
3.2 Specific Features . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 186
3.2.1 Limited "Width" and Enhanced "Depth" . . . . . . . . . . . . . . . . . . . . . . . . . . . 186
3.2.2 Handling of Uncertain, Incomplete and Fuzzy Information . . . . . . . . . . . . 187
3.2.3 Orientation Towards the Processing of Continuous Variables . . . . . . . . . . 187
3.2.4 Capabilities of Intelligent Monitoring . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 187
3.2.5 Availability of C o m m o n "Bioprocess" Knowledge Base . . . . . . . . . . . . . . . 188
3.2.6 Low Prototyping Time . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 188
3.2.7 Low Price . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 189
4 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 189
5 References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 190
The development of intelligent control in biotechnical processes, which only a decade ago
was considered an exciting, but obscure vision, has today become an area of intensive and
realistic research. One of the keys to success in this field is the selection of an adequate
software tool for building the intelligent system. The ideal tool must possess a large set of
features, which reflect both the real-time nature of the control problem, and the peculiarities
of the biotechnical system itself. Here, some of these features are introduced. In addition,
the main concepts and trends in the field of expert control in biotechnical processes are
discussed.
Advancesin BiochemicalEngineering
Biotechnology,Vol. 48
Managing Editor: A. Fiechter
9 Springer-VerlagBerlinHeidelberg 1993
170 K.B. Konstantinov, R. Aarts and T. Yoshida
1 Introduction
Because of the focus on formal analysis and synthesis, many other aspects of
control system design have largely been disregarded in various application areas
[1]. This has been the case with bioprocesses control, even though the limitations
of the traditional control approaches in the field of biological systems have long
been perceived. Nowadays, it is becoming clear that conventional control theory
alone cannot provide the required platform for building high-performance systems
for the control of bioprocesses. Consequently, a trend towards exploitation of new
methods capable of manipulating and utilizing informal knowledge of the biological
plant has emerged. Recent achievements in expert system (ES) technology have
already stimulated research aimed at its application in the field of bioprocesses
control [2].
However, when turning to the ES approach, the bioprocess engineer is likely
to encounter a number of unexpected difficulties, which arise from the still limited
experience in the field of expert control, the fundamental differences from traditional
control methodology, and the development of the software system itself. One of
the tasks which must be solved at an early stage of a project is the selection of
an appropriate expert system development tool (ESDT), also called an "expert
shell". This is crucial to the overall design process; an unsuitable tool may
predestine the research to fruitless effort for years, while selection of the right one
provides a reliable basis for rapid success. Unfortunately, the significance of this
problem is still not sufficiently realized. There is a tendency for people to think
of ESs as one thing, and to consider that any ESDT will suit their particular case
[3]. In fact, current ES technology is leading to the creation of hundreds of products,
which vary tremendously in size, purpose, capability and price. Since a commercial
ESDT specialized for control of bioprocesses has not yet been produced, developers
should make their choice only after careful screening of the tools currently available
on the market. The product selected must provide a large set of features, many
of which are non-trivial and are not supported by the more common ESDTs.
These reflect both the real-time nature of the control problem, and the peculiarities
of the biotechnical system itself.
The principal purpose o f this paper is to introduce a set of ESDT features which
are required in building intelligent systems for the control of bioprocesses. In
addition, the basic structure, functions, and current limitations of ES technology
in the context of problems encountered in the control of bioprocesses, are discussed.
Control system
On-off 1
I ~t~ntrr~llArJ
Ysp _u_> Biotechnical
process
Y
>
Fig. 1. Structure of a system for direct expert control. The knowledge-based (KB) controller
works at the level of the standard controllers, such as PID or "on-off'
whose knowledge-based modules are involved in the control loop (Fig. 1). These
modules operate at the level of standard P I D controllers, and are useful in realizing
more complicated, nonlinear control algorithms. Such modules are known mostly
as fuzzy controllers, because they are often based on fuzzy logic. Recent
achievements in neural networks provide another alternative for the design of
such types of controllers. Generally, the direct expert control scheme, though
useful in some cases, is limited to local and low-level problems. This means that
although its applicability to the control of bioprocesses is possible, it is of limited
significance.
The structure of a system for indirect expert control (also called "supervisory")
is shown in Fig. 2 [4-6]. This is composed of two hierarchical levels providing
clear distribution of the system functions. The standard set of control tasks, such
as measurement, filtration, data acquisition, control, etc., are entrusted to the
lower level, which represents just a conventional control system. Such systems
U IKn~ >
~I supervisory | . Y
| system J =
[__ -r
Conventional 1 Y
control -
system l" O
_J
/
Fig. 2. Structure of a system for indirect
u[ Biotechnical
process
(supervisory) expert control. The knowled-
ge-based module is not involved in the
low-level control loop; instead it supervises
the lower control level by issuing high level
commands
172 K.B. Konstantinov, R. Aarts and T. Yoshida
are included in almost all modern biotechnical equipment, either in the form of
hardware controllers, Or as computer-implemented algorithms. The higher level
represents a knowledge-based superstructure to the conventional control part. In
contrast to the direct scheme, the knowledge-based module, which is indeed a
complete ES, is not explicitly involved in the low-level control (although exceptions
are sometimes possible); it does not generate signals to the control plant, but
helps the lower level to perform its job better. To this end, the higher level issues
supervisory commands which tell the lower level when to do what [7]. As the
techniques for design of low level controllers are relatively well established, the
synthesis of the conventional part of the system is not likely to cause special
difficulties. Thus, the success of the overall system development is dependent
mainly on the design of the knowledge-based module using an appropriate ESDT.
The indirect expert control concept offers new possibilities for development of
high-performance intelligent systems for the control of bioprocesses. This structure
flexibly combines the advantages of the traditional approach with those of ES
technology. It allows enhancement of the control system by the capability of
intelligent decision-making based on informal interpretation of the complex
behavior of the living system. If properly implemented, this structure will be
capable of covering various control problems which usually remain outside the
scope of conventional systems. The subsequent discussion will focus on systems
for indirect expert control.
The features of the ESDT necessary for building the knowledge-based module are
correlated to the functions which this module is expected to perform in real-time.
Although considerable differences are possible according to the application, these
functions can be summarized into four main groups:
- Identification of the state of the cell population. This includes continuous
on-line evaluation of the physiological state of the cell population, informal
interpretation of cell behavior, prediction of future states, detection and diagnosis
of expected (e.g. stage transfers) or unexpected (e.g. deviations from normal
behavior) physiological phenomena, and others. Undoubtedly, this is the most
important, difficult and advanced function of the higher system level [8, 9].
- Identification of the state of the process equipment. The major problem here
is the detection and diagnosis of instrumental failures, such as troubles with sensors
or actuators. It has been pointed out that this function must be explicitly separated
from the previous one [10]. Although this may be not easy in case of biological
plants, it will contribute to the clear modularization of the knowledge base of the ES.
- Supervision of the conventional control part. Identification of the state of the
cell population and the process equipment are passive procedures creating
declarative information on the plant. To become useful, this must be mapped into
high level commands for supervision and synchronization of the work of the
conventional control part. The major purpose is to achieve intelligent handling
and control of the physiological phenomena during the process. Typical supervisory
Expert Systems in Bioprocess Control: Requisite Features 173
According to the method of coupling the knowledge-based part and the conven-
tional control part, two types of implementation scheme are possible:
- Interface scheme. According to this, the knowledge-based part resides on a
dedicated computer, interfaced to the conventional control part (Fig. 3). Both parts
are software and hardware independent of each other. Typically, the higher level
is purchased (as an ESDT), developed and added separately, after the conventional
part has been set in operation.
Today, due to lack of other realistic alternatives, the interfaced scheme is the
standard solution of the bioprocess control problem [4, 14, 15]. This is because
the ESDTs available so far are large programs designed with the interfaced
structure in mind, meant to run on stand-alone computers. Furthermore, modern
bioreactors are equipped with computerized low-level controllers which perform
some of the tasks of the conventional control module, and are easy to link to the ES.
Though popular, an interfaced architecture has several shortcomings: the
software of both system parts is provided by different suppliers, which hampers
integration; communication between computers may cause problems; the informa-
tion storage and the man-machine interface of the system is redundant [10, 16].
- Embedded scheme. In this case, the knowledge-based module is embedded
into the real-time control environment, co-existing with the conventional control
174 K.B. Konstantinov, R. Aarts and T, Yoshida
~Knowledge-based~| "~
supervisory ~/ --~
system~/
Conventional I -~
control ~ y $
system o~
,_1
jf
7
/ / now[edge-base
/ >[ supervisory |4
i I ] modu{e I
f
|~Integrated (hybrid)l| I I . . . .
/I contro[systern 1/ ~---1 t'~
t,.. . . . ~2~ ~l j I
~ 1 INIl ~ I.~l---~- - - ~ -
Biotechnical
process I
Fig. 4. Embedded scheme. The knowledge-based supervisory module coexists with the
conventional control system on the same computer. The knowledge-based module is
embedded into the real-time control environment
Expert Systems in Bioprocess Control: Requisite Features 175
algorithms on the same computer (Fig. 4). The integration between both parts is
much more efficient, and the problems of the interfaced scheme are overcome [17].
Although the embedded scheme is considered to be more advanced, its implementa-
tion is still difficult because the commercial ESDTs available today are not designed
for such integration [18].
i.e. allow the attachment of a specific activation interval to every rule [23].
Furthermore, it is desirable to be able to change dynamically these time intervals
from withing the rule conclusions, for example 1
1 For simplicity, the exemplary rules are written in natural language. In real cases, they
should conform to the ESDT syntax.
Expert Systems in Bioprocess Control: Requisite Features 177
products which do provide some real-time capabilities [28], and it is on these that
the attention of the developed should be entirely focused.
It is worth noting, that to a great extent the real-time capabilities of the ESDT
are an inheritance from those of the operating system; if the operating system
is not real-time, the ESDT will lack real-time features too. Therefore, tools
designed for non-real time operating systems (e.g. DOS), are basically in-
appropriate for control applications. From this viewpoint, the most suitable
ESDTs are those working in real-time UNIX (or UNIX-like) environments [29].
where Time-period defines a period from the process history (often the most recent
one) either explicitly, or in respect to a certain event in the past, Variable is a
particular process variable, and Descriptor represents the expected pattern. The rule
illustrates a simple form of reasoning in respect to the current value of the specific
growth rate SGR and the recent history of the respiration quotient RQ. A slightly
more complicated case is represented by the rule
whose condition considers the recent trend of the SGR. There are many cases
when the information from a biotechnical system can be adequately interpreted
by accounting for the trends of variables. However, this is not always sufficient
to analyze the behavior of bioprocesses in which the temporal shapes of the
variables provide indispensable information about the underlying physiological
phenomena. To achieve better control, these shapes must also be properly
interpreted. This would be quite probable because the decisions of a human
operator are based both on the current values and the recent historical profiles
of the process variables. The rule
178 K.B. Konstantinov, R. Aarts and T. Yoshida
0.10 I
i
I
0.08
0.06
r I
I
0.04 I
n.. I
I
0.02 I
I. At
I
I i
I i
t o- 2 t o- 1 to
Time (h)
Fig. 5. Characteristic time profile of the Ra/g indicating intensive excretion of acetic acid
during cultivation of Escherichia coli
represents such a case. It accounts for the shape of the variable Ra/g (ratio of the
ammonia feed rate to the glucose uptake rate), and will execute its conclusion if
this shape matches the descriptor "'INCREASED CONCAVELY" (Fig. 5). The
shape descriptors may be different, e.g. "PASSED 0 VER MAXIMUM", "STOP-
PED INCREASING", "STOPPED DECREASING", "DECREASING MONO-
TONICALLY", "OSCILLATING", possibly stored in a user-extendible library
of shapes which the ES can recognize [10, 30, 31, 32, 33].
Another useful form of reasoning about the history of variables utilizes the
data over the given time period for calculation of certain statistical markers, and
decision making based on their values [16]:
100
\
9..-. 50
a
o
At
a to- 5 to
Time (rain)
100[ I
I
I
I
where GFR is the glucose feed rate (see Fig. 7). To process this rule, the system
should remember when the continuous glucose feeding began, i.e. to "time-stamp"
~.50- o100I
I
I
I
/DO I GFR
g2s
,,
I
I
O- 01 I
t1- 5 t1
Time (min)
Fig. 7. Typical time profiles of the DO and the GFR upon transfer from batch to continuous
glucose-limited cultivation. If the initial GFR is not sufficient, DO will remain high, which
can be used to correct the GFR
180 K.B. Konstantinov, R. Aarts and T. Yoshida
this event. Such time-stamping of all detected events is another typical feature of
the real-time ES, which always stores facts in the knowledge-bases together with
the time of their appearance.
A slightly different form of temporal reasoning is represented by the introduction
of lower and upper time bounds (permitted interval) for expected events. The
system must then check whether or not the occurrence of the event fulfills these
time constraints, for example
set of rules. They are usually based on some associative criteria, for example
all rules which are related to a particular class of objects or category of problems
[24]. The capabilities of driving the inference engine to narrow its scope of
interest, and to use only part of the knowledge, are sometimes jointly referred to
as "focus-of-attention" feature.
Another consistent way to structure knowledge is to assign priorities to the
rules. The priority represents the importance of the particular rule; a rule will be
checked only after all other rules with higher priority have been processed. This
mechanism is useful because, by assigning different priorities to different groups
of rules, one part of the knowledge base can be forced to work before another [36].
Apart from the logical clarity, there are two important advantages of knowledge
structuring. First, there is a dramatic improvement in the performance speed of
the system due to the reduced number of rules which require processing at any
moment. Second, knowledge structuring simplifies the debugging of the knowledge
base.
If (R~/gis HIGH)
Then (Reduce GFR by A GFR)
which illustrates how to detect and how to prevent the excretion of acetic acid
in a glucose-limited Escherichia coli cultivation [47]. This rule is a typical example
of shallow knowledge. Nevertheless, it can work (and it really does!) quite
satisfactorily. It represents, however, a '"shortcut" through knowledge, which
eliminates all intermediate steps, forming the causal link between the condition
and the conclusion. The underlying physiological phenomena remain hidden, and
the observed fact is mapped directly into the control decision. However, although
very fast and memory-efficient, such rules reduce the capabilities of the ES. They
treat the plant as a black-box, and restrict any possibility of explaining the observed
phenomena. The eliminated intermediate conclusions will not be available in the
knowledge base, and cannot be used by other rules in the system.
To describe the phenomena in more detail, the above rule can be expanded
into a chain of rules. These will show that when the glucose feeding is high the
cell oxidative capacity is exceeded, and cells start to excrete acetic acid, which
results in a raise in the Ra/g value. For more precise quantitative description of
the phenomenon, a mathematical model can be added to the rules. Undoubtedly,
the detailization may continue, involving more and more deep knowledge, down
to the level of enzymes and pathways. The question is whether or not this is useful
or efficient. Generally, excessive overloading of the ES with ballast knowledge is
also harmful. As the control problem itself is usually formulated at a rather high
level of abstraction, extremely detailed description can hardly be useful. Fur-
thermore, the vast amount of knowledge will put a serious load on the computer,
causing an unacceptable slow-down in the real-time environment. On the other
hand, in some cases inclusion of deep knowledge may cause some heuristics to
become redundant, resulting in a net reduction of the size of the knowledge base
[44]. Therefore, the depth of the knowledge must be well balanced with the control
purposes, the desired explanation capabilities, and the time and memory con-
straints.
184 K.B. Konstantinov, R. Aarts and T. Yoshida
This feature is relatively well developed today. The main difficulties ensue from
the necessity to combine simplicity of the man-machine interface with the advanced
performance of the ES. In bioprocess control, it is extremely important to maintain
simplicity because potential users tend to lack specialized computer skills.
Fortunately, modern software techniques for graphical representation, menu-
driven dialogs, icon definitions, simplified data input using a mouse, separation
of the developer interface from the end-user interface, and others, provide a fairly
good balance between simplicity and sophistication.
A major role in the user interface is played by the graphical capabilities of the
ESDT. Since in bioprocesses the main part of the information is represented by
continuous variables, graphical features are of utmost importance [52]. Today,
these are supported by almost all advanced ESDTs, though some companies are
selling the graphics packages as a separate product. Apart from the representation
of process variables, modern systems use graphics for the development and
representation of schemes, tables and diagrams. These help in providing more
realistic visualization o f the situation in the control plant [24], or in representing
the rule base as an easy-to-understand network [36, 37].
186 K.B. Konstantinov, R. Aarts and T. Yoshida
In addition to the general features discussed above, there are some issues related
more specifically to the control of bioprocesses. They reflect the unique charac-
teristics of biological systems, which differ from those of conventional control
plants, and influence the selection of the ESDT.
base will be large. The number of rules will probably range between several dozen
and a few hundred [42, 62, 63]. Thus, it is not likely that ESs for the control of
bioprocesses will grow too much in width.
The reduced size of the knowledge base does not, however, mean simplicity.
To interpret intricate and vague situations, the form of rules and the structure of
the knowledge base will be complex. Consequently, a compact and flexible ESDT,
providing a rich set of functional capabilities (referred to as the "depht" of the
ES) will best suit the problem of interest.
4 Conclusions
Today, there is no single ESDT which combines all the features discussed above,
and since such a perfect tool may not be available soon, the developer should
select the most convenient product from those currently on the marked. The most
appropriate are ESDTs with build-in real-time capabilities, such as NEXPERT
OBJECT, ART, KEE, and particularly those designed for process monitoring and
control, e.g. ESCORT, R'TIME, PICON, MUSE, COMDALE and G2 (the last
being considered to be the most advanced one [16, 28, 64]). Due to their flexibility,
these ESDTs can be customized to a particular application, though such adaption
will not satisfy all specific requirements. Nevertheless, the resulting ES will be
capable of covering a large portion of the problems arising in the control of
bioprocesses.
Certainly, the ideal ESDT in this field would be far more specialized. The lack
of such a product from the market extends the time and effort needed for the
creation of control systems considerably. Since large-scale use of ES technology
in the area of bioprocess control seems inevitable, the development and com-
mercialization of such an ESDT would be an important contribution to the
biotechnology community.
190 K.B. Konstantinov, R. Aarts and T. Yoshida
5 References
35. Haase VH (1990) Preprints l l t h IFAC World Congress, Tallin, USSR, 1:141
36. Payne EC, McArthur RC (1990) Developing expert systems, Wiley, New York
37. Aarts, RJ, Suviranta A, Rauman-Aalto P, Linko P (1989) Food Biotechnology 4:301
38. MontellanoR, Bernier M, CheryA, FarzaM(1990)Preprints llthIFACWorldCongress,
Tallin, USSR, 7:54
39. Johannsen G, Alty JL (1991) Automatica 27:97
40. Kuipers BJ (1986) Artificial Intelligence 29:289
41. Price C, Lee M (1988) Proceed 12th IMACS World Congress, Paris, p 85
42. Liibbert A, Hitzman B, Kracke-Helm H, Sch/igerl K (1988) Proceed IFAC Syrup on
Computer Applications in Fermentation Technology, Cambridge, UK, p 297
43. Voss H (1988) Proceed IFAC Workshop on AI in Real-Time Control, Clyne Castle,
UK, p l
44. Aarts RJ (1989) Proceed International Workshop on Expert Systems in Biotechnology,
Helsinki, Finland, p 108
45. Beaumont SF, Leyval L, Gentil S (1989) Proceed IFAC Conf Advanced Information
Processing in Automatic Control, Nancy, France, p 181
46. Smith JM, Stutely R (1988) SGML: The users guide, Ellis Horwood, Chichester
47. Konstantinov KB, Nishio N, Seki T, Yoshida T (1991) J Ferment Bioeng 71:350
48. Cragun BJ, Stendel HJ (1987) Int J Man-Mach Stud 26:633
49. Walker A (1986) IBM J Res Develop 30 (January), p 2
50. Wick MR, Slagle JR (1989) IEEE Expert 4 (Spring), p 26
51. Chandrasekaran B, Mittal S (1983) Int J Man-Machine Studies (November), p 425
52. Stephanopoulos G (1987) Chemical Engineering Progress 83:44
53. Stephanopoulos G, Stephanopoulos G (1986) Trends in Biotechnology 4:241
54. Linko S, Aarts R (1990) Abstracts of the European Congress in Biotechnology,
Copenhagen, p 326
55. Thibault J, Van Breusegem V, Chwruy A (1990) Biotechnol Bioeng 36:1041
56. Rich E (1990) IEEE Expert 5:5
57. Bielawski L, Lewland R (1991) Intelligent System Design: Integrating Hypermedia and
Expert Systems Technologies, Wiley. New York
58. San Giovanni J, Romans H (1987) Chem Eng Progress 83:52
59. Bristol EH (1987) Proceed IEEE Conf on Decision and Control, Los Angeles, p 1948
60. Moore RL (1985) Control Engineering (April), p 118
61. Satch PA, Paterson AM, Turner MH (1986) Expert Systems 3:22
62. Chen Q, Wang S, Wang J (1988) Proceed IFAC Symp on Computer Applications in
Fermentation Technology, Cambridge, UK, p 253
63. Xu E, Xu G, Zhang S (1989) Proceed IFAC Workshop on A1 in Real-Time Control,
Shenyang, PRC, p 97
64. Eldeib HK (1989) Proc American Control Conference, Pittsburgh, Pennsylvania, p 2637
Author Index Volumes 1-48
Bagnarelli , P., Clementi, M.: Serum-Free Growth of Human Hepatoma Cells. Vol. 34, p. 85
Bailey, J. E.: Host-Vector Interactions in Escherichia coli. Vol. 48, p. 29
Barford, J.-P., see Harbour, C. Vol. 37, p. 1
Barker, A. A., Somers, P. J.: Biotechnology of Immobilized Multienzyme Systems. Vol. 10, p. 27
Barnes, T. J. see Karube, 1. Vol. 46, p. 63
Beardmore, D. H. see Fan, L. T. Vol. 14, p. 101
Bedetti, C., Cantafora, A.: Extraction and Purification of Arachidonic Acid Metabolites from Cell
Cultures. Vol. 35, p. 47
Belfort, G. see Heath, C. Vol. 34, p. 1
Belfort, G. see Heath, C. A. Vol. 47, p. 45
Beker, M. J., Rapoport, A. J.: Conservation of Yeasts by Dehydration. VO1. 35, p. 127
Bell, D. J., Hoare, M., Dunnill, P.: The Formation of Protein Precipitates and their Centrifugal
Recovery. Vol. 26, p. 1
Berlin, J., Sasse, F.: Selection and Screening Techniques for Plant Cell Cultures. Vol. 31, p. 99
Binder, H. see Wiesmann, U. Vol. 24, p. 119
Bjare, M.: Serum-Free Cultivation of Lymphoid Cells. Vol. 34, p. 95
Blanch, H. W., Dunn, I. J.: Modelling and Simulation in Biotechnical Engineering. Vol. 3, p. 127
Blanch. H. W., see Moo-Young, M. Vol. 19, p. 1
Blanch, H. W., see Maiorelle, B. Vol. 20, p. 43
Blaszczyk, R. see Kosaric, N. Vol. 42, p. 27
Blenke, H. see Seipenbusch, R. Vol. 15, p. 1
Blenke, H.: Loop Reactors. Vol. 13, p. 121
Bliem, R. F., Konopitzky, H. W., Katinger, H. W. D.: Industrial Animal Cell Reactor Systems:
Aspects of Selection and Evaluation. Vol. 44, p. 1
194 Author Index Volumes 1 48
Dabora, R. L., Cooney, C. L.: Intracellular Lytic Enzyme Systems and Their Use for Disruption of
Eseherichia colL Vol. 43, p. 1
Daniel, R. M., see Coolbear, T. Vol. 45, p. 57
Daoud, L S. see Atkinson, B. Vol. 4, p. 41
Dos, K. see Ghose, T. K. Vol. 1, p. 55
Davis, P. J. see Smith, R. V. Vol. 14, p. 61
Deckwer, W.-D. see Schumpe, A. Vol. 24, p. 1
Demain, A. L.: Overproduction of Microbial Metabolites and Enzymes due to Alteration of
Regulation. Vol. 1, p. 113
Diamond, A. D., Hsu, J. T.: Aqueous Two- phase Systems for Biomoleculare Separation, Vol. 47, p. 89
Dijkhuizen, L. see de Bes, L. Vol. 41, p. 1
Doelle, H. W., Ewings, K. N., Hollywood, N. W.: Regulation of Glucose Metabolism in Bacterial
Systems. Vol. 23, p. 1
Doelle, H. W. see Johns, M. R. Vol. 44, p. 97
Doran, P. M.: Design of Reactors for Plant Cells and Organs. Vol. 48, p. 115
Dubois, D. see Miller, O. A. Vol. 39, p. 73
Dunn, I. J. see Blanch, H. W. Vol. 3, p. 127
Dunn, I. J. see Heinzle, E. Vol. 48, p. 79
Dunnill, P. see Bell, D. J. Vol. 26, p. 1
Author Index Volumes 1 48 195
Eckenfelder Jr., W. W., Goodman, B. L., Englande, A. J.: Scale-Up of Biological Wastewater
Treatment Reactors. Vol. 46, p. 145
Economidis, I.: An Overview of the Biotechnology research Activities in the European Community.
Vol. 46, p. 225
Einsele, A., Fiechter, A.: Liquid and Solid Hydrocarbons. Vol. 1, p. 169
Electricwala, A. see Griffiths, J. B. Vol. 34, p. 147
Enari, T. M., Markkanen, P.: Production of Cellulolytic Enzymes by Fungi. Vol. 5, p. 1
Enatsu, T., Shinmyo, A.: In Vitro Synthesis of Enzymes. Physiological Aspects of Microbial Enzyme
Production Vol. 9, p. ! 11
Endo, I., Nagamune, T., Tachikava, S., Inaba, H.: A Tubular Bioreactor for High Density
Cultivation of Microorganisms. Vol. 42, p. 1
Endo, I.: A Human Genome Analysis System (HUGA-1) Developed in Japan. Vol. 46, p. 103
Enfors, S-O., Hellbust, H., K6hler, K., Strandberg, L., Veide, A.: Impact of Genetic Engineering on
Downstream Processing of Proteins Produced in E. coli. Vol. 43, p. 31
Engels, J., Uhlmann, E.: Gene Synthesis. Vol 37, p. 73
Englande, A. J. see Eckenfelder Jr., W. W. Vol. 2, p. 145
Eriksson, K. E.: Swedish Developments in Biotechnology Based on Lignocellulose Materials. Vol.
20, p. 193
Eriksson, K.-E. L. see Vallander L. Vol. 42, p. 63
Esser, K.: Some Aspects of Basic Genetic Research on Fungi and Their Practical implications. Vol. 3, p.
69
Esser, K., Lang-Hinrichs, Ch.: Molecular Cloning in Heterologous Systems, Vol. 26, p. 143
Ewings, K. N. see Doelle, H. W. Vol. 33, p. t
Faith, W. T., Neubeck, C. E., Reese, E. T.: Production and Application of Enzymes. Vol. 1, p. 77
Fan, L. S. see Lee, Y. H. Vol. 17, p. 131
Fan, L. T., Lee, Y.-H., Beardmore, D. H.: Major Chemical and Physical Features of Cellulosic
Materials as Substrates for Enzymatic Hydrolysis. Vol. 14, p. 101
Fan, L. T., Lee, Y.-H., Charppuray, M. M.: The Nature of Lignocellulosics and Their Pretreatments
for Enzymatic Hydrolysis. Vol. 23, p. 155
Fan, L. T. see Lee, Y.-H. Vol. 17, p. 101 and p. 131
Faust, U., Sittig, W.: Methanol as Carbon Source for Biomass Production in a Loop Reactor.
Vol. 17, p. 63
Fiechter, A.: Physical and Chemical Parameters of Microbial Growth. Vol. 30, p. 7
Fiechter, A., Gm/inder, F. K.: Metabolic Control of Glucose Degradation in Yeast and Tumor Cells.
Vol. 39, p. 1
Fiechter, A. see Einsele, A. Vol. 1, p. 169
Fiechter, A. see Janshekar, H. Vol. 27, p. 119
Fiechter, A., Sonnleitner, B.: Impacts of Automated Bioprocess Systems on Modern Biological
Research. Vol. 46, p. 143
Finocchiaro, T., Olson, N. F., Richardson, T.: Use of Immobilized Lactase in Milk Systems. Vol. 15,
p. 71
Flaschel, E. see Wandrey, C. Vol. 12, p. 147
Flaschel, E., Wandrey Ch., Kula, M.-R.: Ultrafiltration for the Seperation of Biocatalysts. Vol. 26,
p. 73
Flickinger, M. C. see Gong, Ch.-S. Vol. 20, p. 93
Fowler, H. W. see Atkinson, B. Vol. 3, p. 221
Friehs, K., Reardon, K. F.: Parameters Influencing the Productivity of Recombinant E. coli
Cultivations. Vol. 48, p. 53
Fukui, S., Tanaka, A.: Application of Biocatalysts Immobilized by Prepolymer Methods. Vol. 29, p. 1
Fukui, S., Tanaka, A.: Metabolism of Alkanes by Yeasts. Vol. 19, p. 217
Fukui, S., Tanaka, A.: Production of Useful Compounds from Alkane Media in Japan, Vol. 17, p. 1
Fuller, R. C. see Brandl. H. Vol. 41, p. 77
Furusaki, S. see Chang, H. N. Vol. 44, p. 27
Furusaki, S., Seki, M.: Use and Engineering Aspects of Immobilized Cells in Biotechnology. Vol. 46,
p. 161
196 Author Index Volumes 1-48
Haferberg, D., Hommel, R., Claus, R., Kleber, H.-P.: Extraceltular Microbial Lipids as Bio-
surfactants. Vol. 33, p. 53
HahIbrock, K., Schr6der, J., Vieregge, J.: Enzyme Regulation in Parsley and Soybean Cell Cultures,
Vol. 18, p. 39
I-Iakanson, H. see Mattiasson, B. Vol. 46, p. 81
Haltmeier, Th.: Biomass Utilization in Switzerland. Vol. 20, p. 189
Hampel, W,: Application of Microcomputers in the Study of Microbial Processes. Vol. 13, p. 1
I-Iarbour, C., Barford, J.-P., Low, K>S.: Process Development for Hybridoma Cells. Vol. 37, p. 1
Harder, A., Roels, J. A.: Application of Simple Structures Models in Bioengineering. Vol. 21, p. 55
Hardman, N.: Recent Developments in Enzyme and Microbial Biotechnology - Strategies in
Bioprocess Design. Vol. 40, p. 1
Harrison, D. E. F., Topiwala, H. H.: Transient and Oscillatory States of Continuous Culture. Vol. 3,
p. 167
Heath, C., Belfort, C.: Immobilization of Suspended Mammalian Cells: Analysis of Hollow Fiber
and Microcapsute Bioreactors. Vol. 34, p. 1
Heath, C. A., Belfort, G.: Synthetic Membranes in Biotechnology: Realities and Possibilities. Vol. 47,
p. 45
Hedman, P. see Janson, J.-C. Vol. 25, p. 43
Heinzle, E.: Mass Spectrometry for On-line Monitoring of Biotechnological Processes. Vo. 35, p. 1
Heinzle, E., Dunn, I. J., Ryhiner, G. B.: Modeling and Control for Anaerobic Wastewater Treatment.
Vol. 48, p. 79
Hellebust, H. see Enfors, S.-O. Vol. 43, p. 31
Hermes, H. F. M. see kamphuis, J. Vol. 42, p. 133
Hesolt, H.: Genetics and Genetic Engineering of the Industrial Yeast Yarrowia lipolytica. Vol. 43,
p. 43
Ho, Ch., Smith, M. D., Shanahan, J. F.: Carbon Dioxide Transfer in Biochemical Reactors. Vol. 35,
p. 83
Hoare, M. see Bell, D. J. Vol. 26, p. 1
Hi~ke, H. see Syldatk, Ch. Vol. 41, p. 29
Hofland, A. see Aynsley, M. Vol. 48, p. 1
Hofmann, E. see Kopperschl/iger, G. Vol. 25, p. 101
Hallb, J. see Nyeste, L. Vol. 26, p. 175
Author Index Volumes 1-48 197
Imanaka, T.: Application of Recombinant DNA Technology to the Production of Useful Bio-
materials. Vol. 33, p. 1
lnaba, H. see Endo, I. Vol. 42, p. 1
Inculet, I. I. see Zajic, J. E. Vol. 22, p. 51
Jack, T.-R., Zajic, J. E.: The Immobilization of Whole Cells. Vol. 5, p. 125
Jallageas, J.-C., Arnaud, A., Galzy, P.: Bioconversions of Nitriles and Their Applications. Vol. 14,
p. 1
Jang, C. M., Tsao, G. T.: Packed-Bed Adsorption Theories and Their Applications to Affinity
Chromatography. Vol. 25, p. 1
Jang, C. M., Tsao, G. T.: Affinity Chromatography. Vol. 25, p. 19
Jansen, N. B., Tsao, G. T.: Bioconversion of Pentose to 2,3-Butanediol by Klebsiella pneumonia.
Vol. 27, p. 85
Janshekar, H., Fiechter, A.: Lignin Biosynthesis, Application, and Biodegradation. Vol. 27, p. 119
Janson, J.-C., Hedman, P.: Large-Scale Chromatography of Proteins. Vol. 25, p. 43
Jeffries, Th, W.: Utilization of Xylose by Bacteria. Yeasts, and Fungi, Vol. 27, p. 1
Jiu, J.: Microbial Reactions in Prostaglandin Chemistry, Vol. 17, p. 37
Johns, M. R., Greenfield, P. F., Doelle, H. W.: Byproducts from Zymomonas mobilis. Vol. 44, p. 97
MacLeod, A. J.: The Use of Plasma Protein Fractions as Medium Supplements for Animal Cell
Culture. Vol. 37, p. 41
Magee, R. J., Kosaric, N.: Bioconversion of Hemicellulosics. Vol. 32, p. 61
Maiorella, B., Wilke, Ch. R., Blanch, H. W.: Alcohol Production and Recovery. Vol. 20, p. 43
Mizlek, I.: Present State and Perspectives of Biochemical Engineering: Vol. 3, p. 279
Maleszka, R. See Schneider, H. Vol. 27, p. 57
Author Index Volumes 1-48 199
Nagai, S.: Mass and energy Balances for Microbial Growth Kinetics. Vol. 11, p. 49
Nagamune, T. see Endo, I. Vol. 42, p. 1
Nagatani, M. see Aiba S. Vol. 1, p. 31
Nakajima, H. see Tanaka, A. Vol. 42, p. 97
Nakamura, I. see Kamihara, T. Vo[. 29, p. 35
Neubeck, C. E. see Faith, W. T. Vol. 1, p. 77
Neirinck, L. see Schneider, H. Vol. 27, p. 57
Nguyen, A. L. see Luong, J. H. T. Vol. 47, p. 137
Niederauer, M. Q., Glatz, C. E.: Selective Precipitation. Vol. 47, p. 159
Nielsen, J.: Modelling the Growth of Filamentous Fungi. Vol. 46, p. 187
Nyeste, L., P6cs, M., Sevella, B., Hol16, J.: Production of L-Tryptophan by Microbial Processes,
Vol. 26, p. 175
Nyiri, L. K.: Application of Computers in Biochemical Engineering. Vol. 2, p. 49
Radlett, P. J.: The Use Baby Hamster Kidney (BHK) Suspension Cells for the Production of Foot
and Mouth Disease Vaccines. Vol. 34, p. 129
Radwan, S. S., Mangold, H. K.: Biochemistry of Lipids in Plant Cell Cultures. Vol. 16, p. 109
Ramasubramanyan, K., Venkatasubramanian, K.: Large-Scale Animal Cell Cultures: Design and
Operational Considerations. Vol. 42, p. 13
Ramakrishna, D.: Statistical Models of Cell Populations. Vol. 11, p. 1
Rapoport, A. J. see Beker, M. J. Vol. 35, p. 127
Reardon, K. F. see Friehs, K. Vol. 48, p. 53
Reese, E. T. see Faith, W. T. Vol. 1, p. 77
Reese, E. T., Mandels, M., Weiss, A. H.: Cellulose as a Novel Energy Source. Vol. 2, p. 181
Rehh~ek, Z.: Ergot Alkaloids and Their Biosynthesis. Vol. 14, p. 33
Rehm, H.-J., Reiff, I.: Mechanism and Occurrence of Microbial Oxidation of Long-Chain Alkanes,
Vol. 19, p. 175
Reiff, I. see Rehm, H.-J. Vol. 19, p. 175
Reinhard, E., Alfermann, A. W.: Biotransformation by Plant Cell Cultures. Vol. 16, p. 49
Reiser, J., Glumoff, V., K/ilin, M., Ochsner, U.: Transfer and Expression of Heterologous Genes in
Yeasts Other Than Saccharomyces cerevisiae. Vol. 43, p. 75
Reuveny, S. see Shahar, A. Vol. 34, p. 33
Richardson, T. see Finocchiaro, T. Vol. 15, p. 71
Righelato, R. C. see Pace, G. W. VoL 15, p. 41
Rivera, S. L. see Karim, M. N. Vol. 46, p. 1
Roberts, M. F. see Anderson, L. A. Vol. 31, p. 1
Roels, J. A. see Harder, A. Vol. 21, p. 55
Rogers, P. L.:. Computation in Biochemical Engineering. Vol. 4, p. 125
Rogers, P. L., Lee, K. L., Skotnicki, M. L., Tribe, D. E.: Ethanol Production by Zymomonas Mobilis.
Vol. 23, p. 37
Rolz, C., Humphrey, A.: Microbial Biomass from Renewables: Review of Alternatives. Vol. 21, p. 1
Rosazza, J. P. see Smith, R. V. Vol. 5, p. 69
Rosenberg, M. Z. see Prokop, A. Vol. 39, p. 29
Rosevear, A., Lambe, C. A.: Immobilized Plant Cells. Vol. 31, p. 37
Ryhiner, G. B. see Heinzle, E. Vol. 48, p. 79
Schneider, H., Maleszka, R., Neirinck, L., Veliky, I. A., Chan, Y. K., Wang, P. Y.: Ethanol Production
from D-Xylose and Several Other Carbohydrates by Pachysolen tannophiluys. Vol. 27, p. 57
Shoemaker, H. E. see Kamphuis, J. Vol. 42, p. 133
Schrbder, J. see Hahlbrock, K. Vol. 18, p. 39
Schumpe, A., Quicker, G., Deckwer, W. D.: Gas Solubilities in Microbial Culture Meida. Vol. 24, p. l
SchiigerI, K.: Oxygen Transfer Into Highly Viscous Media. Vol. 19, p. 71
Schfigerl, K.: Characterization and Performance of Single- and Multistage Tower Reactors with
Outer Loop for Cell Mass Production, Vol. 22, p. 93
Schfigerl, K., Oels, U., L/icke, J.: Bubble Column Bioreactors. Vol., p. 1
Schiigerl, K., Lficke, J., Lehmann, J., Wagner, F.: Application of Tower Bioreactors in Cell Mass
Production. Vol. 8, p. 63
Schiigerl, K. see Singh, A. Voi. 45, p. 29
Schwab, H.: Strain Improvement in Industrial Microorganisms by Recombinant DNA
Techniques.Vol. 37, p. 129
Seipenbusch, R., Blenke, H.: The Loop Reactor for Cultivating Yeast on n-Praffin Substrate Vol. 15,
p. 1
Seki, M. see Furusaki, S. Vol. 46, p. 161
Sevella, B. see Nyeste, L. Vol. 26, p. 17
Shahar, A., Reuveny, S,: Nerve and Muscle Cells on Microcarriers Culture. Voi. 34, p. 33
Shanahan, J.-F. see Ho, Ch. S. Vol. 35, p. 83
Shieh, W. K., Keenan, J. D.: Fluidized Bed Biofilm Reactor for Wastewater Treatment. Vol. 33,
p. 131
Shimizu, S. see Yaman~, T. Vol. 30, p. 147
Shinmyo, A. see Enatsu, T. Vol. 9, p. 111
Shioya, S.: Optimization and Control in Fed-Batch Bioreactors. Vol. 46, p. 111
Singh, A., Kumar, P. K. R., Sch/igerl, K.: Bioconversion of Cellulosic Materials to Ethanol by
Filamentous Fungi. Vol. 45, p. 29
Sitti 9, W., see Faust, U. Vol. 17, p. 63
Skotnicki, M. L. see Rogers, P. L. Vol. 23, p. 37
Smith, M. D. see Ho, Ch. S. Vol. 35, p. 83
Smith, R. V., Acosta Jr., D., Rosazza, J. P.: Cellular and Microbial Models in the Investigation of
Mammalian Metabolism of Xenobiotics. Vol. 5, p. 69
Smith, R. V., Davis, P. J.: Induction of Xenobiotic Monooxygenases. Vol. 14, p. 61
Soda, K. see Yonaha, K. Vol. 33, p. 95
Soda, K. see Ohshima, T. Vol. 42, p. 187
Solomon, B.: Starch Hydrolysis by Immobilized Enzymes. Industrial Application. Vol. 10, p. 131
Somers, P. J. see Barker, S. A. Vol. 10, p. 27
Sonnleitner, B. Biotechnology of Thermophilic Bacteria: Growth, Products, and Application.
Vol. 28, p. 69
Sonnleitner, B. see Fiechter, A. Vol. 46, p. 143
Spiegel-Roy, P., Kochba, J.: Embryogenesis in Citrus Tissue Cultures. Vol. 16, p. 27
Spier, R. E.: Rficent Developments in the Large Scale Cultivation of Animal Cells in Monolayers.
Vol. 14, p. 119
Stephanopoulos, G. N. see Grampp, G.E. Vol. 46, p. 35
Stewart, G. G. Kosaric, N. Vol. 20, p. 119
yon Stockar, U., Marison, I. W.: The Use of Calorimetry in Biotechnology. Vol. 40, p. 93
Stohs, S. J.: Metabolism of Steroids in Plant Tissue Cultures. Vol. 16, p. 85
Strandberg, L. see Enfors, S.-O. Vol. 43, p. 31
Sudo, R., Aiby, S.: Role and Function of Protozoa in the Biological Treatment of Poluted Waters.
Vol. 29, p. 117
Suijidam, van, J. C. see Metz, N. W. Vol. 11, p. 103
Sureau, P.: Rabies Vaccine Production in Animal Cell Cultures. Vol. 34, p. 111
Syldatk, Ch., L/iufer, A., Mfiller, R., H6ke, H.: Production of Optically Pure D- and L-m-Amino
Acids by Bioconversion of D,L-5-Monosubstuted Hydantoin Derivatives. Vol. 41, p. 29
Szezesny, T. see Volesky, B. Vol. 27, p. 101