JB
Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Rabinovitch, A.
Right arrow Articles by Zaritsky, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Rabinovitch, A.
Right arrow Articles by Zaritsky, A.

 Previous Article  |  Next Article 

Journal of Bacteriology, March 1999, p. 1677-1683, Vol. 181, No. 5
0021-9193/99/$04.00+0
Copyright © 1999, American Society for Microbiology. All rights reserved.

Model for Bacteriophage T4 Development in Escherichia coli

Avinoam Rabinovitch,1 Hilla Hadas,2 Monica Einav,2 Zeev Melamed,3 and Arieh Zaritsky2,*

Departments of Physics1 and Life Sciences,2 and the Computer Center,3 Ben-Gurion University of the Negev, Be'er-Sheva, Israel 84105

Received 3 March 1997/Accepted 10 December 1998


    ABSTRACT
Top
Abstract
Text
References

Mathematical relations for the number of mature T4 bacteriophages, both inside and after lysis of an Escherichia coli cell, as a function of time after infection by a single phage were obtained, with the following five parameters: delay time until the first T4 is completed inside the bacterium (eclipse period, nu ) and its standard deviation (sigma ), the rate at which the number of ripe T4 increases inside the bacterium during the rise period (alpha ), and the time when the bacterium bursts (µ) and its standard deviation (beta ). Burst size [B = alpha (µ - nu )], the number of phages released from an infected bacterium, is thus a dependent parameter. A least-squares program was used to derive the values of the parameters for a variety of experimental results obtained with wild-type T4 in E. coli B/r under different growth conditions and manipulations (H. Hadas, M. Einav, I. Fishov, and A. Zaritsky, Microbiology 143:179-185, 1997). A "destruction parameter" (zeta ) was added to take care of the adverse effect of chloroform on phage survival. The overall agreement between the model and the experiment is quite good. The dependence of the derived parameters on growth conditions can be used to predict phage development under other experimental manipulations.


    TEXT
Top
Abstract
Text
References

Studies on bacteriophage growth and development in the 1940s played a vital role in the history of molecular biology (11, 24, 29). The classical one-step growth experiment (17) defined latent period, rise time, and burst size, and the eclipse period was discovered by disrupting infected bacteria before their spontaneous lysis (14). By the time bacterial physiology was established as a discipline (25, 32), molecular biology had become so attractive that some unsolved questions in phage-host cell interactions have been ignored and never seriously looked at since. The vast amount of knowledge gained during the last 35 years on the biochemistry, genetics, and physiology of bacteria (23, 26, 30) enables a fresh look on these interactions, which may shed light on various cell properties (15, 20).

In a typical one-step growth experiment, a culture of cells is mixed with phage suspension at a low multiplicity of infection to guarantee single infections. Samples are withdrawn periodically and plated on a lawn of sensitive bacteria, and the number of phages is calculated from the number of plaques formed after overnight incubation. This straightforward procedure has recently been used to describe the development of the T4 bacteriophage inside Escherichia coli under varying well-defined physiological states of the host. The dependence of phage growth parameters on cell size, age, and shape, on rates of metabolism and chromosome replication, and on time of lysis was evaluated semiquantitatively (20). In this series of experiments, the parameters obtained were indeed distributed over wide ranges; they were however derived "by eye." To obtain well-defined quantitative values of the parameters, they should be defined rigorously, and the complete time dependence of the process should be calculated and compared with the experiment. Such a quantitative model is developed here, with due note taken of the statistical distributions of the parameters within the populations of both phages and bacteria.

The former model of a one-step growth experiment (3, 12, 17) implicitly assumed that the latent period ends prior to cell burst and that the different numbers of PFU per bacterium (PPBs) obtained by titration as a function of time were therefore due to the different burst times of individual cells only. Since these times are normally distributed, the PPB should have increased as erfc[(q - t)/Delta ], where q is the burst time (latent period) and 2Delta is the width of the normal distribution. Such an assumption would however entail a step function for chloroform titration of bacteria, which is not observed. The present model avoids such a limitation.

The model. All times are measured from the moment of infection of an E. coli cell by a single wild-type bacteriophage T4. If no spread occurs in the times of onset and termination of phage multiplication, the number of PPB is ideal, as shown schematically in Fig. 1a through d. The time nu  (eclipse period; 15 min in the example shown in Fig. 1a) is the average delay between infection and appearance (inside the cell) of the first complete phage. From this time onwards, the PPB is assumed to increase linearly during the rise period with a constant rate alpha  (8 per min; Fig. 1a). The number of PPB (phi ) stops increasing when the bacterium bursts, at the end of the latent period µ (on average, 30 min in the example shown in Fig. 1b), reaching a final value (burst size) of B = alpha (µ - nu ) (120 phages in the example shown in Fig. 1b). The PPB inside an infected cell before its burst (phi 1) is shown in Fig. 1c, and the number which emerges from lysed cells (phi 2) is shown in Fig. 1d. Evidently, phi 1 + phi 2 = phi . The experiments (20) yield both phi  and phi 2 separately; phi 2 is obtained by titrating phages in the suspension (17) while phi  is obtained after lysing the cells artificially (4), such as with chloroform.


View larger version (17K):
[in this window]
[in a new window]
 
FIG. 1.   A schematic ideal situation of the increase in the number of T4 phages per Escherichia coli bacterium (PPB), with parameters as given in the text. (a) PPB if bacteria do not burst; (b) total PPB, in lysed and nonlysed cells; (c) PPB in nonlysed cells; (d) PPB in lysed cells (phages in suspension). Panels e, f, and g display schematically the nonideal cases ("real" situation) of those shown in panels b, c, and d, respectively.

These experiments, however, are performed for a population which is not ideal and for which µ and nu  are distributed around their averages. Thus, in the population, nu  and µ vary from one infected cell to another and from one infecting phage to another; the actual curves for the phi s are thus smoother. These are shown schematically in Fig. 1e through g and are calculated below. In the model, both nu  and µ are assumed to have Gaussian distributions around their averages with coefficients of variation sigma /nu and beta (of 20%, 3 and 6 min, respectively, in the example shown in Fig. 1). The probability that the phages would start to multiply between t' and t' + dt' is thus considered to be
p<SUB>1</SUB>(t′)<UP>d</UP>t′=<FENCE><RAD><RCD>&pgr;</RCD></RAD>&sfgr;</FENCE><SUP><UP>−1</UP></SUP> exp[<UP>−</UP>(t′−&ngr;)<SUP>2</SUP>/&sfgr;<SUP>2</SUP>]<UP>d</UP>t′ (1)
Note that we follow the notation of reference 1 of a Gaussian distribution. Thus, sigma  in our case measures the time increment (positive or negative) with respect to nu , where probability decreases by l/e relative to the peak probability. A different notation is sometimes used (32), in which sigma 1 = sigma 21/2 appears instead of our sigma .

In bacteria where phages started to multiply at time t', the PPB at t > t' will be alpha (t - t'). Therefore, the average PPB at time t is given by
m(t)=&agr;<LIM><OP>∫</OP><LL><UP>−∞</UP></LL><UL>t</UL></LIM> p<SUB>1</SUB>(t′)(t−t′)<UP>d</UP>t′ (2)
provided, of course, that no bacterium bursts up to t. Using well-known mathematical results (1),
m(t)=&agr;&sfgr;ierfc[(&ngr;−t)/&sfgr;]/2 (3)
where ierfc(z) is the integral of the error function erfc(z). Equation 3 can also be written as
m(t)=<FR><NU>&agr;&sfgr;</NU><DE>2</DE></FR> <FENCE><FR><NU>t−&ngr;</NU><DE>&sfgr;</DE></FR> <FENCE>1+erf<FENCE><FR><NU>t−&ngr;</NU><DE>&sfgr;</DE></FR></FENCE></FENCE>+<FR><NU>1</NU><DE><RAD><RCD>&pgr;</RCD></RAD></DE></FR> exp[(t−&ngr;)<SUP>2</SUP>/&sfgr;<SUP>2</SUP>]</FENCE> (4)
Now, to determine the probability that a bacterium bursts between t' and t' + dt' we use
p<SUB>2</SUB>(t′)<UP>d</UP>t′=<FENCE><RAD><RCD>&pgr;</RCD></RAD>&bgr;</FENCE><SUP><UP>−1</UP></SUP> exp[<UP>−</UP>(t′−&mgr;)<SUP>2</SUP>/&bgr;<SUP>2</SUP>]<UP>d</UP>t′ (5)
The average PPB in the lysed bacteria at t is therefore
&phgr;<SUB>2</SUB>(t)=<LIM><OP>∫</OP><LL><UP>−∞</UP></LL><UL>t</UL></LIM> p<SUB>2</SUB>(t′)m(t′)<UP>d</UP>t′ (6)
where the integrand measures PPB for the bacteria that burst between t and t + dt. Using equations 4 and 5, 
&phgr;<SUB>2</SUB>(t)=<FENCE>&agr;&sfgr;/2<RAD><RCD>&pgr;</RCD></RAD>&bgr;</FENCE><LIM><OP>∫</OP><LL><UP>−∞</UP></LL><UL>t</UL></LIM>ierfc[(&ngr;−t′)/&sfgr;]exp[<UP>−</UP>(t′−&mgr;)<SUP>2</SUP>/&bgr;<SUP>2</SUP>]<UP>d</UP>t′ (7)
The average PPB in bacteria which have not lysed by the time t is given by
&phgr;<SUB>1</SUB>(t)=<FENCE>1−<LIM><OP>∫</OP><LL><UP>−∞</UP></LL><UL>t</UL></LIM> p<SUB>2</SUB>(t′)<UP>d</UP>t</FENCE>m(t) (8)
where the relation in the brackets is the probability that the bacterium did not burst by t and m(t) (equation 2) is the PPB if the bacterium indeed did not lyse. Hence,
&phgr;<SUB>1</SUB>(t)=<FR><NU>&agr;&sfgr;</NU><DE>2</DE></FR><FENCE>1−<FR><NU>1</NU><DE>2</DE></FR>erfc<FENCE><FR><NU>&mgr;−t</NU><DE>&bgr;</DE></FR></FENCE></FENCE>ierfc<FENCE><FR><NU>&ngr;−t</NU><DE>&sfgr;</DE></FR></FENCE> (9)
and the total PPB is
&phgr;(t)=&phgr;<SUB>1</SUB>(t)+&phgr;<SUB>2</SUB>(t) (10)
As mentioned before, the measured quantities are phi (t) and phi 2(t).

In order to check these relations, we use the asymptotic forms of the functions appearing here (1). Thus for t >> nu , (nu - t)/sigma , which we denote by -u, is negative, and
ierfc(<UP>−</UP>u) → u(1+erfu)+e<SUP><UP>−</UP>u<SUP>2</SUP></SUP>/<RAD><RCD>&pgr;</RCD></RAD> (11)
For u right-arrow infinity we have erfu right-arrow 1 and e-u2 right-arrow 0; therefore ierfc(-u) right-arrow 2u, and m(t) (equation 2) right-arrow alpha (t - nu ), as it should (Fig. 1). Hence,
&phgr;<SUB>2</SUB>(t)≈<FENCE>&agr;/<RAD><RCD>&pgr;</RCD></RAD>&bgr;</FENCE><LIM><OP>∫</OP><LL><UP>−∞</UP></LL><UL>t</UL></LIM>(t′−&ngr;)exp[<UP>−</UP>(t′−&mgr;)<SUP>2</SUP>/&bgr;<SUP>2</SUP>]<UP>d</UP>t′ (12)
= <FENCE>&agr;/<RAD><RCD>&pgr;</RCD></RAD></FENCE><LIM><OP>∫</OP><LL>−<UP>∞</UP></LL><UL>(t<UP>−&mgr;</UP>)<UP>/&bgr;</UP></UL></LIM> (&bgr;u+&mgr;−&ngr;)exp(<UP>−</UP>u<SUP>2</SUP>)<UP>d</UP>u
where u = (t- µ)/beta . And after some algebra,
&phgr;<SUB>2</SUB>(t)≈&agr;(&mgr;−&ngr;)−(&agr;&bgr;/2){ierfc[(t−&mgr;)/&bgr;]+ (13)
[(t−&ngr;)/&bgr;]erfc[(t−&mgr;)/&bgr;]}
For t >> µ, the second term right-arrow0 and phi 2 right-arrow alpha (µ - nu ), as it should (Fig. 1). Similarly,
1−<FR><NU>1</NU><DE>2</DE></FR> erfc<FENCE><FR><NU>&mgr;−t</NU><DE>&bgr;</DE></FR></FENCE>→ 0 <UP>for  </UP>t >>&mgr;
and thus by equation 9, phi 1 right-arrow 0.

Methods. The model was tested against results (20) obtained with wild-type T4 phage (5) infecting either E. coli B/r (H266) (36) or E. coli K12 (CR34, thr leu thy drm) (35). Bacteria were cultured with vigorous shaking at 37°C in the following media (36): Luria-Bertani broth containing glucose (LBG) (0.4%; doubling time tau  = 23 min); M9 minimal medium supplemented with casein hydrolysate (1%), tryptophan (50 µg/ml), and glucose (GC) (0.4%; tau  = 28 to 30 min) or with 0.4% of either glucose (tau  = 48 min), glycerol (tau  = 70 min), or succinate (tau  = 90 min). For the thyA strain of E. coli K12, 5 µg of thymine/ml was added, either with deoxyguanosine (100 µg/ml) or not (35). Upon achieving a steady state (18, 34) at a concentration of 108 cells ml-1, cells were infected with phage (or treated prior to infection as described previously [20]) at a multiplicity of 0.5 (to guarantee a single infection) in the presence of 2 mM KCN (to synchronize the infective process). Phage development was initiated 4 min after phage addition by dilution (10-4) in the same medium to cease further infections and eliminate the bacteriostatic effect of the cyanide. Samples were withdrawn periodically and plated immediately and through chloroform (after evaporation) after appropriate dilutions, with soft agar using E. coli B/r (H266) as indicator. The number of phages was calculated from the number of plaques formed after overnight incubation at 37°C. The raw data (in PFU per milliliter) were transformed to derive the number of phages per infected cell as a function of time (20). The number of PFU obtained in the chloroform series reflects the concentration of all phages, both free and newly matured inside infected bacteria. PFU obtained without chloroform stem from either ripe phages or infected bacteria (infective centers), whether harboring mature phages or not. The number of unadsorbed phages during the eclipse period is considered the background PFU value and is thus subtracted from all raw data.

To obtain different cell sizes under the same growth rate, parameters which are usually correlated (32), three regimes were employed (20). (i) Synchronous glucose-grown cells, obtained by the "baby-machine" (22), were infected either upon collection ("babies") or after 40 min of growth (almost one mass doubling). The latter, larger synchronous cells indeed supported a slightly faster phage assembly and yielded more phages than the babies. (ii) Thymine limitation of a thyA mutant strain delays cell division due to a slowing down of the chromosome replication rate and thus results in enlarged cells (20, 35). (iii) Low penicillin-G (Pn) concentrations were used to specifically block division without affecting mass growth rate (19). Exposure during about two tau s before infection resulted in ca. fourfold-larger cells (data not shown). Superinfection (SI) (6) was employed to delay cell lysis, thus extending the period during which phages continue to develop (20).

The numerical data processing method used here was developed in the last decade. The data contain values of phi  and phi 2 at t1, t2, ...tn. The independent model parameters were sought in such a way as to minimize the error. An elaborate least-squares method was used here because the problem is complicated due to the functional form (an integral for which no analytic expression was found).

The programs used to evaluate the parameters were Minpack (Argonne National Laboratory, 1980) and Simusolv (Dow Chemical Co., 1990). The algorithm used in Minpack is a modified Newton one (27, 28), where an approximation is built for the Hessian of the Newton method. It can be noted that in many cases Minpack helps to find the global minimum, independent of the starting point. In Simusolv, the movement towards the minimum is accomplished by combining two subgroups, Search and GRG (Generalized Reduced Gradient).

In the processing performed here, the simpler program, Minpack, was used to obtain good first approximations. Note that Simusolv can fit several (here, two [phi and phi 2]) functions simultaneously and even functions given by their differential equations. The latter was used specifically for phi 2.

Results and discussion. Equations 7, 9, and 10 constitute a system whereby the experimental results can be analyzed. The measurement of total PPB (phi ) is performed by applying chloroform. We have found (data not shown) that the addition of chloroform usually causes a reduction of phage ability to form plaques (plating efficiency) by 5 to 20%, resulting in an effective reduction of PPB which does not change with time. Hence, we assumed that the measured PPB is given by zeta  · phi , where zeta  is the chloroform "destruction parameter." Thus, altogether six independent parameters were calculated: µ, nu , alpha , beta , sigma , and zeta .

The model was tested against previously published results (20). Results appear in Table 1, together with the burst size [B = alpha (µ - nu )]. The large variation in derived estimated parameters might have been caused by a small number of data points. To appreciate the problems in analysis, the classical LBG experiment (17) is presented at a higher measured accuracy (minute-by-minute intervals) (Fig. 2). Even under casual observation it is seen that the scatter of points in the chloroform titration results is quite large. Applying the Simusolv procedure as described above yielded the following estimates for the parameters: nu  = 18.2 ± 0.1; sigma  = 2.93 ± 0.1; µ = 24.11 ± 0.1; beta  = 3.8 ± 0.14; zeta  = 0.82 ± 0.02; and B = 250.4 ± 5.2. The r2 values of the fit, shown by the lines in Fig. 2, are however rather poor. For the regular titration, we get 95.2, while for the chloroform titration it is only 90.3. To improve these estimates, we used the moving average method (see, e.g., reference 8), which smooths out the results. The average of the first three points was calculated and taken as the second value, the average of the second to fourth points was used for the third value, etc. The parameters estimated following this averaging process were as follows: nu  = 18.3 ± 0.3; sigma  = 3.37 ± 0.2; µ = 23.6 ± 0.1; beta  = 4.3 ± 0.15; zeta  = 0.84 ± 0.01; and B = 247.1 ± 5.3 (Fig. 3), with an obvious improvement in r2 (regular titration, 97.8; chloroform, 97.4). It is clearly seen that such a smoothing transformation is quite beneficial. In addition, all parameter values other than beta  remained essentially unchanged, indicating that they were quite robust and reliable.

                              
View this table:
[in this window]
[in a new window]
 
TABLE 1.   Parameters, at different growth conditions, calculated by the analysis from experimental dataa


View larger version (14K):
[in this window]
[in a new window]
 
FIG. 2.   Comparison between the model (lines) and results of a classical experiment (17, 20). *, samples titered after chloroform treatment. +, naturally released phage.


View larger version (13K):
[in this window]
[in a new window]
 
FIG. 3.   Comparison between the model (lines) and results of the same experiment as in Fig. 2, with a moving average (3-points) technique. (See text.)

Note that the results (Table 1) were obtained with the same six parameters used for both zeta  · phi  and phi 2(t) in each case. The overall agreement is quite good (see the Statistical Analysis section below). The value of zeta  that is larger than 1.0 cannot be simply explained; more data points in the plateau region may solve this apparent discrepancy.

Figure 4 presents three comparisons between the model and actual experiments previously reported (20) (Table 1). The three examples displayed were selected to cover the whole range of burst sizes (between 9 and 720); thus, a semilog presentation was used. (For clarity, the corresponding experiments with chloroform were not included, and the points at early times were deleted.) Correlations between burst size and cell size were observed in early studies with the T-series bacteriophages (see, e.g., references 13 and 21), but it was too early for them to be accounted for by the physiological parameters of E. coli, which emerged a decade later (25, 32).


View larger version (13K):
[in this window]
[in a new window]
 
FIG. 4.   Comparisons between the model (lines) and results of three experiments previously described (20). open circle  and +, steady-state cultures (17) grown in LBG and succinate-minimal media, respectively. *, cells grown in LBG medium, treated with Pn (75 U/ml) (19) for 60 min and superinfected (6) 6 min after primary infection at a multiplicity of 10.

Figure 5 presents the results of an experiment not previously described, with two samples of a glucose-grown synchronous culture obtained by the "baby machine" (22); one infected upon collection (babies) and the other infected after 40 min of growth (almost one mass doubling). The latter supported a slightly faster phage assembly and yielded more phages (burst size of about 60) than the babies (B = ca. 33). The results are qualitatively consistent with the hypothesis proposed before (20) that the rate of phage synthesis and assembly is proportional to the size of the protein-synthesizing system (9) in the cell upon its infection. This hypothesis will be rigorously tested, and results will be published separately.


View larger version (14K):
[in this window]
[in a new window]
 
FIG. 5.   PPB in synchronous E. coli B/r cells (22) pregrown in glucose-minimal medium, infected (at a multiplicity of 0.5) with bacteriophage T4. + and *, eluted "babies"; open circle  and ×, "old" cells, after 40 min of growth following elution. + and open circle , without chloroform; * and ×, with chloroform.

Statistical analysis. Since the number of points in each experiment is small (of the order of 40 for the chloroform [phi ] and nonchloroform [phi 2] parts combined), the statistical significance of our results is not expected to be excellent. As we presently show, however, some parameters have a higher statistical validity than others. We analyze here two cases (the LBG assay from Fig. 3 and the babies from Fig. 5), while other cases show similar statistical attributes. Consider first the percentage of variation explained, which is a measure equivalent to r2 in our case. As mentioned, for the average LBG, an r2 of ca. 97.5% was obtained, while for Fig. 5 the results were, for regular titration, 96.9% and, for chloroform, 95.1%, with an overall fit of 96.0%. A result of ~96 to 97% is fair and seems higher than would have been expected, since the parameters have to explain both phi  and phi 2.

The standard deviations (Table 2) were calculated from the diagonal elements of the variance-covariance matrix (see below). The results show (i) the much-improved accuracy of the LBG case due to the increased number of measured points and to the averaging procedure and (ii) better accuracy of µ, nu , and zeta  relative to the B and of all these relative to sigma  and beta . For the baby cells, the latter are not really determined. The larger errors in beta  and sigma  (and see the discussion below of the correlation matrix) stems from the fact that these parameters are determined by the curvatures of the PPBs at the beginning and the end of their increase. Since the numbers of experimental points at these regions are small, the inability to get exact values of sigma  and beta  seems obvious. This observation is enhanced by the correlation matrix of the baby cells (Table 2), where it is seen that the covariance between µ and beta  is of the order of -0.9.

                              
View this table:
[in this window]
[in a new window]
 
TABLE 2.   Statistical results of the young "baby" cells (Fig. 5)

To summarize the statistical results, it seems that for the experiments carried out here, the validity of the model is fairly established. The parameters µ, nu , and zeta  have a higher statistical significance than B, and much higher significance than alpha  and beta . The latter is not even robust. Higher accuracies for both were obtained with an increased number of measured points.

Calculation of correlation matrix. To calculate the correlation matrix, the following steps are taken (7, 27, 28, 31). (i) The log likelihood function is calculated by
&PHgr;=<FR><NU>1</NU><DE>2</DE></FR> <LIM><OP>∑</OP><LL>j = 1</LL><UL>r</UL></LIM> <FENCE>n<SUB>j</SUB>[<UP>ln</UP>(2&pgr;)+1]+n<SUB>j</SUB><UP>ln</UP><FENCE><FR><NU>1</NU><DE>n<SUB>j</SUB></DE></FR> <LIM><OP>∑</OP><LL>i = 1</LL><UL>n<SUB>j</SUB></UL></LIM> <FR><NU>(z<SUB>i,j</SUB>−f<SUB>i,j</SUB>)<SUP>2</SUP></NU><DE>f<SUB>i,j</SUB><SUP>&ggr;<SUB>j</SUB></SUP></DE></FR></FENCE>+&ggr;<SUB>j</SUB> <LIM><OP>∑</OP><LL>i = 1</LL><UL>n<SUB>j</SUB></UL></LIM> <UP>ln</UP>f<SUB>i,j</SUB></FENCE>
where zij is the measured value of the j response of the ith data point, fij is the predicted j response of the ith data point, gamma 1 is the heteroscedasticity parameter for the jth response, r is the number of measured response variables and nj is the number of data points of the j response.

(ii) The Hessian matrix Hk1k2 is defined as the matrix of the second partial derivatives of Phi  with respect to each pair of parameters. Here, the Gauss approximation is used:
H<SUB>k<SUB>1</SUB>,k<SUB>2</SUB></SUB>= <LIM><OP>∑</OP><LL>j = 1</LL><UL>r</UL></LIM> <LIM><OP>∑</OP><LL>i = 1</LL><UL>n<SUB>j</SUB></UL></LIM> <FR><NU>∂<SUP>2</SUP>&PHgr;(&thgr;)</NU><DE>∂f<SUB>i,j</SUB><SUP>2</SUP></DE></FR> <FENCE> <LIM><OP>∑</OP><LL>p = 1</LL><UL>r</UL></LIM> <FR><NU>∂f<SUB>i,j</SUB></NU><DE>∂z<SUB>i,p</SUB></DE></FR> <FR><NU>∂z<SUB>i,p</SUB></NU><DE>∂&thgr;<SUB>k<SUB>1</SUB></SUB></DE></FR> </FENCE> <FENCE> <LIM><OP>∑</OP><LL>p = 1</LL><UL>r</UL></LIM> <FR><NU>∂f<SUB>i,j</SUB></NU><DE>∂z<SUB>i,p</SUB></DE></FR> <FR><NU>∂z<SUB>i,p</SUB></NU><DE>∂&thgr;<SUB>k<SUB>2</SUB></SUB></DE></FR> </FENCE>
where theta  is the vector of adjustable parameters.

(iii) The variance-covariance matrix V is estimated from the inverse of the Hessian, V = H-1.

(iv) The correlation matrix is given by the normalized variance-covariance matrix of the parameters estimated by
r<SUB>k<SUB>1</SUB>k<SUB>2</SUB></SUB>=<FR><NU>V<SUB>k<SUB>1</SUB>k<SUB>2</SUB></SUB></NU><DE><RAD><RCD>V<SUB>k<SUB>2</SUB>k<SUB>2</SUB></SUB> V<SUB>k<SUB>1</SUB>k<SUB>1</SUB></SUB></RCD></RAD></DE></FR>
with k1,k2 = 1, 2, ..., m, where m is the number of adjustable parameters.


    ACKNOWLEDGMENTS

We thank Itzhak Fishov for constructive discussions.

This work was partially supported by grant 91-00190/2 from the U.S.-Israel Binational Science Foundation (BSF), Jerusalem (to A.Z.), and by a Ben-Gurion Fellowship administered by the Ministry of Science and the Arts (to H.H.).


    FOOTNOTES

* Corresponding author. Mailing address: Department of Life Sciences, Ben-Gurion University of the Negev, P.O. Box 653, Be'er Sheva 84105, Israel. Phone: 972-7-646.1712. Fax: 972-7-627.8951. E-mail: ariehz{at}bgumail.bgu.ac.il.


    REFERENCES
Top
Abstract
Text
References

1. Abramovitz, M., and I. A. Stegun. 1965. Handbook of mathematical functions. Dover Publishers Inc., New York, N.Y.
2. Adams, M. H. 1959. Bacteriophages. Interscience Publishers Inc., New York, N.Y.
3. Adams, M. H., and F. E. Wassermann. 1951. Frequency distribution of phage release in the one-step growth experiment. Virology 2:96-108.
4. Ang, A., and W. Tang. 1975. Probability concepts in engineering, planning and design. Wiley, New York, N.Y.
5. Benzer, S. 1959. On the topology of the genetic fine structure. Proc. Natl. Acad. Sci. USA 45:1607-1620[Free Full Text].
6. Bode, W. 1967. Lysis inhibition in Escherichia coli infected with bacteriophage T4. J. Virol. 1:948-955[Abstract/Free Full Text].
7. Bord, Y. 1974. Nonlinear estimation. Academic Press, New York, N.Y.
8. Box, G. E. P., and G. M. Jenkins. 1976. Time series analysis: forecasting and control, revised ed. Holden-Day, Inc., San Francisco, Calif.
9. Bremer, H., and P. P. Dennis. 1996. Modulation of chemical composition and other parameters of the cell by growth rate, p. 1553-1569. In F. C. Neidhardt, R. Curtiss III, J. L. Ingraham, E. C. C. Lin, K. B. Low, B. Magasanik, W. S. Reznikoff, M. Riley, M. Schaechter, and H. E. Umbarger (ed.), Escherichia coli and Salmonella: cellular and molecular biology, 2nd ed. American Society for Microbiology, Washington, D.C.
10. Brown, A. 1956. A study of lysis in bacteriophage-infected Escherichia coli. J. Bacteriol. 71:482-490[Free Full Text].
11. Cairns, J., G. S. Stent, and J. D. Watson (ed.). 1966. Phage and the origins of molecular biology. Cold Spring Harbor Laboratory, Cold Spring Harbor, N.Y.
12. Delbrück, M. 1940. Adsorption of bacteriophage under various physiological conditions of the host. J. Gen. Physiol. 23:631-642[Abstract/Free Full Text].
13. Delbrück, M. 1945. The burst size distribution in the growth of bacterial viruses (bacteriophages). J. Bacteriol. 50:131-135[Free Full Text].
14. Doermann, A. H. 1948. The intracellular growth of bacteriophage. Carnegie Inst. Wash. Year Book 47:176-182.
15. Eddy, S. 1992. Introns in the T-even bacteriophages. Ph.D. thesis. University of Colorado, Boulder.
16. Eisenstark, A. 1967. Bacteriophage techniques. Methods Virol. 1:449-524.
17. Ellis, E. L., and M. Delbrück. 1939. The growth of bacteriophage. J. Gen. Physiol. 22:365-384[Abstract/Free Full Text].
18. Fishov, I., A. Zaritsky, and N. B. Grover. 1995. On microbial states of growth. Mol. Microbiol. 15:789-794[Medline].
19. Hadas, H., M. Einav, I. Fishov, and A. Zaritsky. 1995. Division inhibition capacity of penicillin in Escherichia coli is growth-rate dependent. Microbiology 141:1081-1083[Abstract].
20. Hadas, H., M. Einav, I. Fishov, and A. Zaritsky. 1997. Bacteriophage T4 development depends on the physiology of its host Escherichia coli. Microbiology 143:179-185[Abstract].
21. Heden, C. 1951. Studies of the infection of E. coli B with the bacteriophage T2. Acta Pathol. Microbiol. Scand. Suppl. 88:42-62.
22. Helmstetter, C. E. 1969. Methods for studying the microbial division cycle, p. 317-364. In J. R. Norris, and D. W. Ribbons (ed.), Methods in microbiology, vol. 1. Academic Press, New York, N.Y.
23. Ingraham, J. L., O. Maaløe, and F. C. Neidhardt. 1983. Growth of the bacterial cell. Sinauer Associates, Inc., Sunderland, Mass.
24. Karam, J. D. 1994. Molecular biology of bacteriophage T4. American Society for Microbiology, Washington, D.C.
25. Kjeldgaard, N. O., O. Maaløe, and M. Schaechter. 1958. The transition between different physiological states during balanced growth of Salmonella typhimurium. J. Gen. Microbiol. 19:607-616[Medline].
26. Koch, A. L. 1971. The adaptive responses of Escherichia coli to a feast and famine existence. Adv. Microb. Physiol. 6:147-217[Medline].
27. Levenberg, K. 1944. A method for the solution of certain problems in least squares. Q. Appl. Math. 2:164-168.
28. Marquardt, D. W. 1963. An algorithm for least squares estimation of nonlinear parameters. SIAM J. 11:431-441.
29. Mathews, C. K., E. M. Kutter, G. Mosig, and P. B. Berget (ed.). 1983. Bacteriophage T4. American Society for Microbiology, Washington, D.C.
30. Neidhardt, F. C., R. Curtiss III, J. L. Ingraham, E. C. C. Lin, K. B. Low, B. Magasanik, W. S. Reznikoff, M. Riley, M. Schaechter, and H. E. Umbarger (ed.). 1996. Escherichia coli and Salmonella: cellular and molecular biology, 2nd ed. American Society for Microbiology, Washington, D.C.
31. Reilly, P. M., and G. E. Blau. 1974. The use of statistical methods to build mathematical models of chemical reacting systems. Can. J. Chem. Eng. 52:289-299.
32. Schaechter, M., O. Maaløe, and N. O. Kjeldgaard. 1958. Dependency on medium and temperature of cell size and chemical composition during balanced growth of Salmonella typhimurium. J. Gen. Microbiol. 19:592-606[Medline].
33. Stent, G. S. 1963. Molecular biology of bacterial viruses. W. H. Freeman and Co., San Francisco, Calif.
34. Woldringh, C. L., P. G. Huls, and N. O. Vischer. 1993. Volume growth of daughter and parent cells during the cell cycle of Saccharomyces cerevisiae a/alpha as determined by image cytometry. J. Bacteriol. 175:3174-3181[Abstract/Free Full Text].
35. Zaritsky, A., and C. L. Woldringh. 1978. Chromosome replication rate and cell shape in Escherichia coli: lack of coupling. J. Bacteriol. 135:581-587[Abstract/Free Full Text].
36. Zaritsky, A., C. L. Woldringh, and D. Mirelman. 1979. Constant peptidoglycan density in the sacculus of Escherichia coli B/r growing at different rates. FEBS Lett. 98:29-32[Medline].


Journal of Bacteriology, March 1999, p. 1677-1683, Vol. 181, No. 5
0021-9193/99/$04.00+0
Copyright © 1999, American Society for Microbiology. All rights reserved.



This article has been cited by other articles:


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Rabinovitch, A.
Right arrow Articles by Zaritsky, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Rabinovitch, A.
Right arrow Articles by Zaritsky, A.


Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
Appl. Environ. Microbiol. Infect. Immun. Eukaryot. Cell
Mol. Cell. Biol. J. Virol. Microbiol. Mol. Biol. Rev.
ALL ASM JOURNALS