OR WAIT null SECS
Using Bezier curves, an experimental process controller has been developed for biosynthesis applications in which the inactivity of a pH-sensitive enzyme must be decreased. By taking into account various control scenarios of pH and growth rate, as well as the physical and chemical characteristics of the environment, a suitable human-machine interface can be developed.
Biosynthesis may require an active enzyme concentration in a cell-culture growth medium. Reduced biosynthesis capability caused by an inactive enzyme can be an economical concern. Controlling the inactivity of pH-sensitive enzymes can help maintain, or possibly increase, a stable microbial growth environment, thereby preventing a shift in the biosynthesis capabilities of the process harvest.
The biosynthesis of pharmaceutical drug products begins by evaluating charts from the product development laboratory that show a drug's production capability relative to pH, time, substrate concentration, enzyme concentration, and so forth. These informational components were the basis for the "Enzyme Wyeth/RoviSys Controller" (EWrC, RoviSys, Aurora, OH) discussed in this article. The algorithms for the EWrC have been developed and tested on a personal computer. The human–machine interface has not yet been developed.
Figure 1: This reaction mechanism involves six kinetic constants, making it difficult to analyze.
To understand how the information provided in product development laboratory charts can fundamentally improve the approach to process control, it is necessary to briefly review the analytical methods used to evaluate biochemical reactions. Biological cellular growth mediums that make up enzyme-catalyzed reactions are difficult to work with analytically for several reasons (1):
Standard state and relaxation-method analysis. Subsequently, physical biochemists have defined a standard state in which all of the activities and concentrations normally used to analyze reaction kinetics are not required. Hammes developed an analytical method known as the relaxation method to study fast chemical reaction components at equilibrium. This method can be briefly stated as follows:
A chemical system at equilibrium is suddenly perturbed in such a way that it finds itself out of equilibrium under the new conditions. It then relaxes to a state of equilibrium under the new conditions, and the rate of relaxation is governed by the rate constants and mechanism of the process (1).
Relaxation-method analysis and steady-state theory describe the physical model of an enzyme-catalyzed reaction shown in Figure 1, with the following precepts:
Departing from conventional control strategies, the principle objective of the EWrC is to maintain the equilibrium by controlling the reaction perturbation defined by the product-development laboratory charts and supported technically by the relaxation methods' physical model for the enzyme-catalyzed reaction of interest.
Transition-state theory. The transition-state theory, supported by the Hammond postulate for Gibbs free energy of mathematical maximums and minimums for reaction growth diagrams, is used to divide the laboratory charts into various reactant intermediate states (see Figure 2) (2). The transition-state theory for enzyme-catalyzed chemical reactions states the following:
Figure 2: Transition-state theory.
In addition, the following portions of the Hammond postulate are of interest:
The importance of transition-state theory is that it relates the rate of a reaction to the difference in Gibbs free energy between the transition state and the ground state.
Selection of the nonlinear control strategy has been approached theoretically (3). As observed by Antje Kann,
Models are often developed with knowledge of theory and algorithms, but who may not be intimately familiar with the users' situations. Modeling choices made are often guided by practical concerns such as solvability, run-time issues, and ease of implementation in code. However, these choices may impair the ease with which a model can be applied to solve the users' problems (4).
Following this train of thought, a fresh approach in the development of the experimental EWrC was taken. The product-development laboratory charts were the best starting point. Well-documented analysis was available for all the process dynamics, measured variables, responses, and unanswerable process questions for the drug of interest. Classical multivariable, nonlinear control strategies provided well-documented information. In addition, a good curve-generating engine would help convert data from the laboratory charts into useful information.
Bezier curves. Parametric equations allow curves to have the unusual ability to move straight up or even curve back and cross themselves (see Figure 3). Cubic Bezier curves belong to the family of parametric curve equations. Cubic Bezier curves were selected as the first set of parametric curves to use because they do not need a tangent vector. Bezier curves are characterized by defining four coordinate points (origin, destination, values of x, and values of y), one linear variable coordinate point t, and two polynomial equations. The two polynomial equations are used to calculate the values of x and y between the increasing variable t.
Figure 3: Example of a parametric curve crossing itself twice.
As a simple example for the development of a cubic Bezier curve (5), suppose two polynomial equations define the coordinate points on the curve. Both are evaluated for an arbitrary number of values of t. In this example t is incremented according to:
t = t + 0.02, with t ≥ 1
which generates 50 values of t. As increasing values for t are provided to the equations, the point defined by x(t), y(t) moves from the origin to the destination according to:
x(t) = axt3 + bxt2 + cxt + x0
x1 = x0 + cx/3
x2 = x1 + (bx + cx)/3
x3 = x0 + ax + bx + cx
y(t) = ayt3 + byt2 + cyt + y0
y1 = y0 + cy/3
y2 = y1 + (by + cy)/3
y3 = y0 + ay + by + cy
The values of the coefficients a, b, and c can be determined with the following equations:
Σcx = 3(x1 – x0)
bx = 3(x2 – x1) – cx
ax = x3 – x0 – cx – bx
Selecting four control cubic coordinates at random:
origin: xp0 = 3.5 and yp0 = 0
control point 1: xp1 =3.5 and yp1 = 1
control point 2: xp2 = 5 and yp2 = 6
destination: xp3 = 6 and yp3 = 3
These four coordinates evaluate to the coefficients: cx = 0, cy = 3, bx = 4.5, by = 12, ax = –2, ay = –12. The first coordinate pair at t = 0.02, provides x(0.02) = 3.50 and y(0.02) = 0.06. Similarly, at t = 0.04, x(0.04) = 3.50 and y(0.04) = 0.14 (similarly for t = 0.06, etc.). Figure 4 shows the dramatic difference between the results of two selections for control coordinates.
Figure 4: Two examples of cubic Bezier curves. Two different control coordinates were selected for (a) and (b). By using four Cartesian coordinates instead of two to define a point on a curve, calculating correctly dependent variables based upon nonlinear high-gain independent variables is deterministic rather than infinity.
Controller types. Using loosely coupled data from laboratory charts and borrowing concepts from classical controllers with multivariable, nonlinear capabilities generally used for pH regulation produced the following filtered list of controller types:
Summary of information requirements. The EWrC now can be fully developed using the following information about the drug's properties, theories, methods, control strategies, and states:
The following detail requirements will serve as the basis for the controller:
The following measurements will be made and their values compared with the desired values:
Functional specifications for the EWrC (see Figure 5) to be developed are as follows:
Table I: EWrC detail specifications.
Detail specifications for the EWrC are shown in Table I and Figure 5. A block diagram for the EWrC is shown in Figure 6. For clarity, alarms and emergency anomalies are not shown.
Figure 5: Example of the humanâmachine interface for the EWrC (see Table I for descriptions of each screen).
To demonstrate the EWrC, the following are four simulated control scenarios for making a batch with the same set of growth rate curves.
Figure 6: Example of a block diagram for the EWrC.
Scenario 1. The "pH dominate" (pHdom) mode will continuously challenge the current pH measurement and test whether it is possible to increase the pH value ahead of the time period determined by the laboratory data (see Figure 7, scenario 1, screens A and B). Referring to screen B, with the EWrC in pHdom mode, the gain increases at almost a linear slope (see screen B). The desired pH ramp is requested to increase at a much faster rate than normal. The growth rate curve is unaffected in this mode.
Figure 7: Examples of four different control scenarios of the humanâmachine interface for (a) pH dominate, (b) pH normal, (c) rate dominate, and (d) rate normal modes.
EWrCoutput = IF(pHdesired at time x/pHactual at time x)
within Standard Deviation THEN pHcurrent(pHdesired at time (x + 1)/pHestimated at time (x + 1)) ELSE pHcurrent(pHmedian at time x)
Scenario 2. The "pH normal" (pHn) mode does not challenge either the predicted product growth time or the pH of the laboratory data (see Figure 7, scenario 2, screens A and B). Referring to screen B, with the EWrC in pHn mode, the gain tracks the growth-rate curve. The growth-rate curve is unaffected in this mode.
EWrCoutput = IF(pHdesired at time x/pHactual at time x)
within Standard Deviation THEN pHcurrent(pHdesired at time x /pHactual at time x) ELSE pHcurrent(pHmedian at time x)
Scenario 3. The "growth rate dominant" (Rdom) mode will continuously challenge the current pH measurement and tests if it is possible to increase the pH value ahead the of the time period determined by the laboratory data (see Figure 7, scenario 3, screens A and B). Referring to Screen B, with the EWrC in Rdom mode, the estimated growth rate and time interval for pH changes is used in combination to determine the gain output. The growth rate desired increases rapidly and then falls into a constant decline. A general low gain and constant gain value is determined necessary for most of the batch process. The pH curve is unaffected in this mode.
EWrCoutput = IF(Rateestimated at time x/Ratedesired at pH actual at time x)
within Standard Deviation THEN Ratecurrent(Rateestimated at time (x + 1)/Ratedesired at pH estimated at time (x + 1)) ELSE Ratecurrent(Rateestimated at pH median at time x)
Scenario 4. The "growth rate normal" (Raten) mode does not challenge either the predicated product growth time or the pH of the laboratory data (see Figure 7, scenario 4, screens A and B). Referring to Screen B, with the EWrC in Raten mode, the estimated growth rate is used to determine the gain output. A very large gain is determined necessary at first and then falls off rapidly to a constant value until the end of the batch process. The pH curve is unaffected in this mode.
EWrCoutput = IF(Rateestimated at time x/Ratedesired at pH actual at time x)
within Standard Deviation THEN Ratecurrent(Rateestimated at time x/Ratedesired at pHactual at time x) ELSE Ratecurrent(Rateestimated at pH median at time x)
An experimental EWrC has been created and debugged to demonstrate and prove the algorithm on a personal computer (see Figure 8). The human–machine interface has not been created and the controller has not been tested in an actual bioprocess.
Figure 8: Example of a screen shot of the experimental EWrC.
The potenital benefits to product quality, especially for injectable drugs, are the best feature of the EWrC. Product quality associated with pharmaceutical bioprocessing is assessed in terms of the finished product's physical and chemical properties. The EWrC and its control functionality are completely based upon the physical biochemical properties developed in the product development laboratory for the product of interest.
The limitations of the EWrC lie in the statements of physical biochemists listed at the beginning of the article. Reaction intermediates produced in situ accompanying the large volumes used in commercial batches could alter relaxation methods responses markedly. In addition, the forward and back rates of reaction can be affected nonlinearly with perturbations.
The operational and controls benefits of the EWrC functionality depend on the specific physical biochemical properties of the drug of interest. On the basis of the product development laboratory charts, there should be less control transitions because the addition of reagents should follow the relaxation method principles more closely. In addition, the system is less affected by concentration anomalies in reagents; control with either biologics or difficult or latent process measurements is possible; and the system is self-optimizing in pH or rate-dominate control modes.
Production economics benefits can be achieved from the various EWrC modes. For the pH or rate-dominate modes, reagent addition should be at an optimum rate, thus increasing product concentration, decreasing batch perturbations, and decreasing batch-production time intervals. In addition, operating close to equilibrium should optimize utilities and reagent usage. For the pH mode, reagent concentration anomalies should not perturb equilibrium because time-to-pH values is used as a parameter. In addition, operating close to equilibrium should optimize utilities and reagent usage. Finally, the rate mode applies primarily for reactions experiencing enzyme-substrate binding that cause latent pH measurements, and operating close to equilibrium should optimize utilities and reagent usage.
Michael J. Needham is a validation engineer at RoviSys, 1455 Danner Drive, Aurora, OH 44202, firstname.lastname@example.org
1. I. Tinoco, Jr., K. Sauer, and J.C. Wang, Physical Chemistry Principles and Applications in Biological Sciences (Prentice-Hall, Englewood Cliffs, NJ, 2d ed., 1985).
2. Alan Fersht, Structure and Mechanism in Protein Science (W. H. Freeman and Company, New York, NY, 3d ed., 2000).
3. Instrument Society of America, Standards and Practices for Instrumentation (ISA, 6th ed., Research Triangle Park, NC, 1980).
4. Antje Kann, president, Analytics Interactive (San Francisco, CA, www.analyticsinteractive.com).
5. "Bezier curves," http://www.math.ubc.ca/~cass/gfx/bezier.html, accessed, Jan. 16, 2006.