|
|
|
Journal of Bacteriology, April 2002, p . 1888-1894, Vol . 184, No . 7 Effects of Escherichia coli Physiology on Growth of Phage T7 In Vivo and In SilicoLingchong You, Patrick F . Suthers, and John Yin* Department of Chemical Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706-1691 Received 20 September 2001/ Accepted 4 January 2002
The relative simplicity of phage developmental processes, compared with those of microbes or higher organisms, is balanced in part by the complexity of the resources that they need for growth . Phage require at least nucleic acid precursors, protein precursors, and translation machinery from their hosts . Consequently, phage infection processes depend not only on the physicochemical characteristics of their genome-encoded functions but also on the intracellular resources of their hosts, which depend further on the physiological state of their hosts . Studies spanning 60 years have demonstrated this dependence (1, 5-8, 11, 15, 16, 18) . Different stages of phage growth, including the attachment of the phage particle to its host, the penetration of phage DNA into its host, and the synthesis of the phage components, have been found to be sensitive to the physiological state of the host, which has been modulated by its growth medium (5, 15, 16), temperature (11, 18), oxygen tension (18), and pretreatment using chemical agents (1, 5, 15, 16) . These studies have shown that the faster that the host grows at the time of phage infection, the faster that the phage will grow, corresponding to a shorter latent time (or a shorter intracellular eclipse time), faster progeny production rate, and larger burst size . How is phage production influenced by the levels of different essential intracellular resources? The answer is not obvious because of the coupling of resources to host growth . As the growth rate of the host increases, so does its cell size (9), as well as its intracellular levels of genomic DNA, RNA polymerase, ribosomes, nucleoside triphosphates (NTPs), and amino acids (3) . These are all essential resources that, to unknown extents, affect phage growth . In addition, the effects of these factors may be convoluted: the increase in the cell volume as the host cell grows faster will affect the concentrations of various resources, which in turn may affect phage growth . Mutant hosts may help uncouple these effects . For example, the effects of DNA levels can be uncoupled from growth rate by growing a thy mutant under thymine limitation (28) . The one-step growth behavior of phage T4 on this mutant was used to infer the key role of the protein synthesis machinery in phage growth (15) . Such results, while useful, do not exclude possible confounding effects of the mutation on other host resources . Using phage T7 as the model system, we explored the effect of the host physiology on phage growth using both experimental and computational approaches . Phage one-step growth experiments were conducted by infecting Escherichia coli cells growing at different rates, achieved by using a chemostat . These experiments suggested a similar dependence of phage growth rate on the host growth rate as observed for other phage studied: the faster that the host grows, the faster the phage infection . To better quantitatively understand this dependence, we employed a computer simulation for the intracellular growth of phage T7 in E . coli to probe the effect of the host physiology at the molecular level . We have previously developed a simulation to account for the mechanisms and rate of entry of the T7 DNA into the host, the transcription of the T7 genes and the translation of the resulting mRNAs, the assembly of the T7 procapsids, the replication of the T7 genome, the packaging of the new genomes into the procapsids, and the maturation of the T7 progeny (12, 13) . Here we extended this simulation and explored the effects of host resources on phage growth by accounting for the coupling or uncoupling of the host resources to its growth rate .
(ii) Intracellular one-step growth experiments. Established methods were used in the preparation, preservation, and assay of the phage (1, 20, 26) . Bottom agar for plates and soft agar for overlayers were LB broth containing 1.5 and 0.7% Bacto agar (Difco), respectively . The host cells went through at least six to eight doublings after reaching the new steady state for each change in flow rate, as determined by optical density at 600 nm . Sampling of the chemostat was performed under sterile conditions . Ten milliliters of the sample was added to a shaker flask already at 30°C, which was then infected with phage T7 . The multiplicity of infection was 0.05, and all phage growth was carried out at 30°C . At 3 min postinfection, a sample from the infected shaker flask was diluted 500-fold into 30 ml of cell-free LB medium to minimize further binding of phage to bacteria . Samples from the diluted, infected cell culture's flask were treated with chloroform (Sigma) to release intracellular phage and titered on BL21 grown in shaker flasks on LB . The number of infective centers was determined as the difference between the chloroform-treated sample and a sample that was not chloroform treated, both taken immediately after dilution . Phage dilutions for titering were performed in buffer containing 10 mM Tris-HCl (pH 7.5), 1 mM MgCl2, 0.1 M NaCl, 10 mg of gelatin/liter, and 10 mM CaCl2 . Phage growth curves were generated in triplicate at each host growth rate .
(iii) Parameter estimation.
The objective function to extract the eclipse time, rise rate, and burst size from the intracellular one-step growth curves was defined as
where
Simulation . (i) Phage T7 model. The present T7 model, T7v2.5, recasts previous T7 models developed in our group (12, 13) using an object-oriented approach . It extends the earlier models by accounting for the stoichiometric relation of the T7 helicase/primase (gp4A) and the DNA polymerase (gp5) in forming replication complexes, or replisomes, as well as the stoichiometric balance between the number of replication complexes and the maximum number of replication forks that can form on the newly synthesized T7 genomes . In addition, it allocates E . coli RNA polymerases (EcRNAPs) and T7 RNA polymerases (T7RNAPs) to the synthesis of different mRNAs based on the relative strengths of the promoters . Further, it incorporates a module that facilitates the sensitivity analysis of the model with respect to single parameters or two parameters simultaneously . While previous T7 models assumed that host resources were unbounded, the present T7 model incorporates an empirically based host cell model that accounts for experimentally observed correlations between the E . coli growth rate and resources such as the numbers of EcRNAP and ribosomes, the pool sizes of NTPs and amino acids, and the cell volume, shown in Table 1 and reviewed elsewhere (3) . In setting up the host cell model, we treated the host cell as a spatially homogeneous resource reservoir, where the levels of the resources were defined at the initiation of the T7 infection . Moreover, by employing the equations in Table 1, we implicitly assumed the following: (i) cell volume was constant over the course of infection, (ii) there was negligible exchange of metabolites between the infected cell and its extracellular environment during infection, (iii) the model host cell represented an average of a cell population that was growing at the exponential phase immediately prior to phage infection, (iv) the initial NTP pool size was equivalent to the total RNA content of the cell, (v) NTPs were not consumed as the energy source for the reactions, and (vi) the initial amino acid pool size was equivalent to the total protein content of the cell . These assumptions were made based on our still limited understanding of the interplay between phage T7 and its E . coli host; some of them are open to debate . For example, as detailed in Discussion, the value for the NTP pool size was chosen mainly because of the lack of mechanistic understanding on the exact level of NTP in the host cell that is available for phage development . Nevertheless, these assumptions serve as a starting point toward representing the host cell, and they can be refined when more experimental data and mechanisms become available .
(ii) Computer simulations. We performed a sensitivity analysis on all the E . coli physiological parameters that affect the growth of T7 (Table 1) . Here we use the term "sensitivity analysis" in its broad sense, i.e., the determination of changes in the value of model variables in responding to the change in a parameter value . The base case value of the host growth rate was set to 1.5 doublings/h, and the corresponding base case values of the other parameters were calculated using the equations in Table 1 . In each sensitivity analysis, one parameter was varied over a range from 0.1 to 10 times its base case value, while the other parameters were kept constant, and a simulation was conducted for each value of the selected parameter . For scenarios where host physiology was defined by a single growth rate, we coupled host resources to the growth rate using the equations in Table 1 . To provide a basis for comparison with the experiments we conducted, the host growth rate was varied from 0.7 to 1.7 doublings per h . Each simulation was conducted over a period of 60 min with a time step of 0.02 s for numerical integration, except for those where the host cell volume was less than half its base case value . In this case, a time step of 0.001 s was used to avoid computational errors resulting from the integration of stiff equations (e.g., the assembly of procapsid is 4.8th order in terms of the concentration of major capsid protein; significant computational errors may occur when this concentration is high, which may arise for small volumes and large time steps) . The sensitivity of phage growth to both EcRNAP and ribosomes was studied by generating a mesh of 400 (or 20 by 20) nodes ranging from 0.1- to 100-fold of the base case values, and the rise rate was determined at each node . Derivatives with respect to EcRNAP number and ribosome number were then determined at each node, and local changes in these derivatives were used to demark transitions in resource limitation . Specifically, the boundary between the EcRNAP-limited growth and ribosome-limited growth was defined between nodes where the derivative of the rise rate with respect to the EcRNAP numbers changed from positive to either negative or zero . Further, the boundary between growth limitations by protein synthesis and DNA synthesis was defined between nodes where the derivative of the rise rate with respect to both EcRNAP numbers and ribosome numbers changed from positive to either negative or zero .
The number and the elongation rate of ribosomes had the strongest effect on the rise rate (Fig . 3) . They had virtually the same effect until they were large, where the rise rate leveled off at different levels . The increase in ribosome elongation rates increased the rise rate more than did the increase in the ribosome numbers . The dependence of the rise rate on the cell volume was overall moderate and, interestingly, biphasic; the rise rate first increased and then decreased as the cell volume increased from 0.1 to 10 times its base case value (Fig . 3) . The optimal rise rate occurred at the base case cell volume .
To better understand how an increase in the EcRNAP number could slow down T7 growth, we examined the effect of the increase in the EcRNAP number on the production of several phage components . We found that, with the other parameters kept constant at their base case values, an increase in the EcRNAP number reduced the rate of procapsid synthesis (Fig . 4A) . This reduction corresponded with significant increases in the number of ribosomes allocated to the mRNAs for the early T7 genes, such as 0.7 (T7 kinase gene) and 1 (T7RNAP gene), and reductions in the number of ribosomes allocated to the mRNAs for the late T7 genes, such as 9 (scaffold protein gene), 10A (major capsid protein gene), and 19 (DNA maturation protein gene) (Fig . 4B) .
Consistent with our experimental results, the simulations predicted that faster-growing host cells supported faster T7 growth (Fig . 1 and 2B and C) . The simulations give overall better predictions on the rise rate than on the eclipse time . The systematic mismatch between predicted and experimental values for eclipse time (Fig . 2C) may result from the instantaneous nature of several "hard-switches" that we implement in the simulation, where data or mechanisms are not available . These switches include the initiation of phage DNA entry, the initiation of T7RNAP-modulated translocation and transcription, and the initiation of host DNA digestion . For instance, we assume a 90-s delay for the initiation of phage infection to account for the time required for phage adsorption, which can take longer . Increasing this delay by 5.3 min can eliminate the mismatch (Fig . 2C, dashed line) . Nonetheless, this correction does not exclude potential contributions from other factors . Although the intracellular one-step growth curves cannot distinguish between these mechanisms, available experimental data on the intracellular processes of T7 infection at the molecular level may provide some hints . A close examination of earlier data on T7 protein expression (14) suggests that the mismatch is unlikely due to a delay in phage absorption because T7 proteins begin to appear around 2 min after infection . Further, since expression of class II proteins initiates in less than 2 min after the appearance of T7RNAPs (14), the second mechanism cannot completely account for the mismatch either . The last mechanism seems to contribute the major part of this mismatch . In fact, we assumed that the host DNA is digested between 7.5 and 15 min after infection . These time points were based on experiments that measured the degradation of host DNA by T7 infection at 37°C (23) and could be significantly delayed at 30°C, the temperature used in the present study . Similar experiments conducted at 30°C will help to explain the mismatch and improve the simulation by providing more accurate parameters . Previous efforts to identify key host constraints to phage growth have been hindered by the experimental challenge of independently modulating and studying the effects of different host resources on phage growth (15, 16) . We have addressed this challenge by using the simulation to numerically uncouple the host intracellular resources from the host growth rate and then independently examine their effects on phage growth . While such studies may appear to distort relationships among the host resources, they are useful for two reasons . First, they allow us to better understand how uncertainties in parameters may affect the behavior of the simulation . While empirically determined cell compositions and activities are well defined, it is not always clear what may be accessible for phage growth . For example, NTP pool size in the simulation is based on the total host cell RNA content, which includes both stable and unstable RNAs (3) . However, we do not know for certain what fraction of the total NTPs per cell is accessible to T7 processes during infection . On the one hand, since stable RNAs constitute parts of ribosomes that are needed by T7 for protein translation, they are improbable precursors for producing T7 mRNAs . On the other hand, the unstable RNAs alone are probably inadequate to support phage growth . Further, the host cell may have means of providing NTPs for phage growth, for example, by continuing to synthesize NTPs during infection; to our knowledge there has been no evidence that T7 infection shuts down host NTP synthesis . Because of the resultant uncertainty in the accessible NTP pool size, it is useful to explore the effects on phage growth of a broad range of values . Second, studying the sensitivity of the simulated phage growth to extreme changes in single-host parameters may reveal aspects of phage development in cells growing in complex environments or in highly distorted hosts like anucleate minicells (21) . By varying host growth rates in a chemostat, we have considered only a narrow range of conceivable host resource distributions . Hosts growing in complex environments, where spatial gradients and dynamic variations in nutrients, temperatures, and pH can be significant, may well have intracellular resource distributions that drastically differ from our experiments . While we lack detailed information about how resources might be correlated in these cases, determining the sensitivity of phage growth to variations in single or paired host parameters can begin to provide insights into potential constraints under such conditions . The sensitivity analysis with the individual host physiological parameters indicates that T7 growth most strongly depends on the number and the elongation rate of the ribosomes, suggesting that the rate-limiting step of phage growth is the synthesis of phage proteins . However, since the ribosome elongation rate only increases slightly as the cell growth rate increases, its contribution is probably less important than the ribosome number, which increases much more rapidly with the cell growth rate (Table 1) . Thus, our simulation results support the hypothesis that the amount per cell of the host translation machinery is one of the most crucial factors in determining the rate of phage growth (15) . This result is also consistent with the weaker dependence of phage growth on the number and elongation rate of EcRNAPs, the sizes of the NTP and amino acid pools, and the DNA content, although all these factors also play indispensable roles during T7 growth . Because of the complex nature of the viral infection dynamics, increasing the amount of an essential component may actually slow down viral growth, as demonstrated by the counterintuitive observation that an increase in the EcRNAP number can decrease the T7 growth rate (Fig . 3) . One may expect that increasing the EcRNAP number would increase mRNA levels, protein synthesis rates, and finally the T7 growth rate . This expectation holds when the EcRNAP number is sufficiently small or when the ribosome number is sufficiently high (data not shown) . In either scenario, T7 growth is limited by transcription and this limitation can be relieved by an increase in the EcRNAP number . When the EcRNAP number is too large, however, T7 growth will be limited by translation . In this case, increasing the number of EcRNAPs, which are mainly responsible for the transcription of T7 class I genes, will result in a higher level of early mRNAs and divert the ribosomes away from the translation of several late proteins that are needed in large quantities for phage particle formation (Fig . 4B), particularly the scaffold protein (gp9) and the major capsid protein (gp10A) . gp10A is especially important for two reasons: first, compared with other T7 proteins, it is needed in the largest quantity for each T7 particle (25); second, the rate of procapsid assembly is assumed to be proportional to [gp10A]4.8 (12, 13) . Hence, even a slight decrease in gp10A level may significantly decrease the procapsid assembly rate (Fig . 4A) and in turn significantly decrease the rise rate . In summary, when translation is already limited by the ribosome number, an increase in EcRNAP number causes ribosomes to distribute in a manner that is unfavorable for the expression of important late proteins, thus reducing the rise rate (Fig . 3) . Our simulations indicate that T7 growth is more sensitive to the processivity of T7RNAPs than to the processivity of EcRNAPs, although the effect of either factor is much weaker than that of the ribosomes (data not shown) . This is to be expected because the transcription of T7 genes is primarily dependent on T7RNAP, which is more efficient than EcRNAP in internalizing and transcribing T7 DNA . This point further suggests that T7 growth limitation by the host translation machinery may well be because the phage provides its own highly efficient RNA polymerase, thus creating a bottleneck for processing of its mRNAs . The biphasic effect of the cell volume on T7 growth is another demonstration of the complexity of the infection dynamics . A large cell volume may slow down the intracellular interactions by causing a decrease in the concentrations of interacting species . This is particularly true for the assembly of procapsids from major capsid proteins (gp10A) . An increase in the cell volume can result in a decrease in the gp10A concentration, which in turn can cause a major decrease in the procapsid assembly rate . When the cell volume is too small, however, other effects may dominate: the inhibition of EcRNAP by gp0.7 and gp2 and the inhibition of T7RNAP by gp3.5 will be enhanced because of the increased total concentrations for these proteins . Consequently, under such conditions fewer free T7RNAP and EcRNAP molecules are available for the transcription of viral genes and the rate of viral protein expression decreases, which in turn reduces the rate of phage progeny formation . This idea is supported by the observation that the numbers of active EcRNAPs and T7RNAPs decrease with a decreasing cell volume (data not shown) . In summary, extremes in cell volume can, through different mechanisms, have detrimental effects on T7 growth . As the activity of host or phage functions vary, so too may that of bottlenecks to phage production, as shown in the "bottleneck landscape" of Fig . 5 . For the base case parameters, T7 growth is bottlenecked or limited by the rate of protein synthesis, which is determined by the level of ribosomes . If one relieves this bottleneck by increasing the ribosome number, then the phage growth becomes limited by the host transcription rate, determined by the level of EcRNAP . When levels of both ribosomes and EcRNAPs are high, then the bottleneck to phage growth can become the rate of its own DNA synthesis . While such a scenario may be unlikely for wild-type phage growth on a wild-type host, due to the >10-fold increases that would be required for both the EcRNAP and ribosome levels, phage mutants that expressed reduced DNA synthetic capacities could alter the landscape by expanding the DNA synthesis-limiting region . In summary, specific features of the bottleneck landscape will change with genetic or environmental modifications to phage or host functions, but so long as the simulation accounts for the modifications, the corresponding landscapes can be created .
This work was supported by the Office of Naval Research (grant no . N00014-98-1-0226), the National Science Foundation Presidential Early Career Award to J.Y . (grant no . BES-9896067), and Merck Research Laboratories . P.F.S . was supported by a National Science Foundation Graduate Fellowship .
What Is Fermentation?,
What is Food Microbiology?,
What Is Water Purification?,
What Is Molecular Microbiology?,
What Is Botulism?,
i,
Microbe,
r,
Bacteria,
r,
Microbes,
i,
Bacterium,
o,
Bacteriology,
s,
S. cerevisiae,
c,
Escherichia coli,
o,
Escherichia coli,
s,
Escherichia coli,
r,
Erythromycin,
i,
Cell suspensions,
i,
Staphylococcus aureus,
a,
Antibiotics,
i,
Pseudomonas aeruginosa,
n,
Microbiological,
s,
Streptococci,
n,
Bactericidal,
i,
Escherichia coli,
a,
Enterobacteriacea,
s,
Escherichia coli,
n,
Pichia,
c,
Bacteriological,
c,
S. cerevisiae,
s,
Clostridia,
r,
S. cerevisiae,
r,
Escherichia coli
|
© 2005
Transgalactic Ltd (manufacturer of Bioscreen C software) |
Privacy Statement | P.O. Box
1393, 00101 Helsinki, Finland,
Last modified: May 25, 2005
| ||||||