Metabolic investigation of host/pathogen interaction using MS2-infected Escherichia coli
© Jain and Srivastava; licensee BioMed Central Ltd. 2009
Received: 27 March 2009
Accepted: 30 December 2009
Published: 30 December 2009
RNA viruses are responsible for a variety of illnesses among people, including but not limited to the common cold, the flu, HIV, and ebola. Developing new drugs and new strategies for treating diseases caused by these viruses can be an expensive and time-consuming process. Mathematical modeling may be used to elucidate host-pathogen interactions and highlight potential targets for drug development, as well providing the basis for optimizing patient treatment strategies. The purpose of this work was to determine whether a genome-scale modeling approach could be used to understand how metabolism is impacted by the host-pathogen interaction during a viral infection. Escherichia coli/MS2 was used as the host-pathogen model system as MS2 is easy to work with, harmless to humans, but shares many features with eukaryotic viruses. In addition, the genome-scale metabolic model of E. coli is the most comprehensive model at this time.
Employing a metabolic modeling strategy known as "flux balance analysis" coupled with experimental studies, we were able to predict how viral infection would alter bacterial metabolism. Based on our simulations, we predicted that cell growth and biosynthesis of the cell wall would be halted. Furthermore, we predicted a substantial increase in metabolic activity of the pentose phosphate pathway as a means to enhance viral biosynthesis, while a break down in the citric acid cycle was predicted. Also, no changes were predicted in the glycolytic pathway.
Through our approach, we have developed a technique of modeling virus-infected host metabolism and have investigated the metabolic effects of viral infection. These studies may provide insight into how to design better drugs. They also illustrate the potential of extending such metabolic analysis to higher order organisms, including humans.
Viruses are the cause of a variety of diseases ranging from the mildly annoying common cold to the frighteningly lethal ebola. Viruses infect their hosts and hijack the host machinery, using it to produce more progeny viral particles. Viruses, being obligate parasites, require host resources to replicate. Therefore, viral infections lead to alterations in the metabolism of the host, shifting in favor of viral protein production. A systems biology approach for studying these metabolic changes in the host cell could provide new insights to potential drug targets , motivating this study.
Systems biology deals with the studying of organisms that are viewed as one collaborative network of genes, proteins and other metabolites. Recent advancements in high-throughput experimental techniques have brought a flood of information in the form of genomic, transcriptomic, proteomic and metabolomic datasets. Systems biology is answering the growing need for the integration and interpretation of these heterogeneous datasets. One of the methodologies of systems biology, a constraints-based modeling approach , has been successfully demonstrated in identifying potential drug targets . Genome-scale metabolic models of disease-causing organisms have been constructed and evaluated to identify potential drug targets [3–5].
In this study, we have demonstrated the application of the constraints-based, flux balance analysis approach to investigate the effects of host-pathogen interactions on host metabolism. The pathogen-host pair under consideration was a bacterial virus, MS2, and its host, Escherichia coli C-3000. MS2 is a lytic RNA bacteriophage, belonging to the Leviviridae family of viruses and infects F+ Escherichia coli cells .
The Escherichia coli/MS2 system was chosen for a number of reasons. MS2 is harmless to humans and yet possesses many of the same features as its eukaryotic-infecting viral cousins, and as a result may aid in our understanding of RNA viruses in general. It can be cultured quickly, cheaply, and safely, making it easy to work with. Furthermore, the genome-scale metabolic model of E. coli is the most exhaustive one to date. These factors combine to make the E. coli/ MS2 model system ideal for such a study.
The constraints-based modeling approach was used to describe only a part of the infection process. Based on an experimental study tracking the macromolecular synthesis in MS2-infected Escherichia coli , the infection process can be divided into 3 parts - an early transient period where the infection process is initiated, a middle steady state period where the viral protein synthesis has displaced the host protein synthesis and a late transient period where all biosynthesis has diminished and lysis is approaching. This study used the constraints-based modeling approach to investigate the middle steady state period only.
The Escherichia coli genome-scale metabolic model, iAF1260  was used as a basis to represent both - the uninfected cells and the infected cells. C-3000, the Escherichia coli strain used in this study and MG1655, the Escherichia coli strain used to model iAF1260 are both Escherichia coli K-12 derivatives. As a result, iAF1260 was used to represent the uninfected cells. The genome-scale metabolic model of MS2-infected E. coli was based on iAF1260, but was modified to account for the viral infection process. The constraints for these models such as the glucose uptake rate, the oxygen uptake rate and the growth rates of the cells were measured experimentally. Finally, the genome scale metabolomes of the infected and the uninfected cells were compared to find what parts of the metabolic network were activated, upregulated, downregulated or deactivated.
For all in silico reconstructions, linear programming was carried out using the GLPSOL (version 4.38) solver, which is available as a part of the GNU Linear Programming Kit . In addition, to capture the effects of alternate optimal solutions, flux variability analysis  was carried out for all in silico reconstructions.
In silico reconstruction of genome-scale metabolic network of uninfected Escherichia coli
The reconstruction of the uninfected host metabolism was based on the model iAF1260 . In brief, the system was defined as the metabolic network that included extracellular, cytoplasmic and periplasmic metabolites. Exchange reactions were responsible for metabolites entering or leaving the system. Transport reactions accounted for the movement of metabolites across the inner and outer membranes. The remaining reactions accounted for the biochemical interactions among metabolites.
where, v represented a vector of reaction fluxes. An optimization technique, linear programming, was applied to find a specific solution to the underdetermined linear problem. Additional constraints in the form of bounds were placed on each of the fluxes. For all the irreversible reactions, the limits were 0 ≤ v i ≤ ∞ and for the reversible reactions, the limits were -∞ ≤ v i ≤ ∞. For the exchange fluxes, the limits of 0 ≤ v i ≤ ∞ were applied for the metabolites that were not present in the growth medium and the limits of -∞ ≤ v i ≤ ∞ for the metabolites that were provided by the growth medium. In addition, the fluxes of the limiting reactions were constrained to the experimentally measured values.
Based on studies that have systematically evaluated various objective functions [13–15] on experimental systems similar to the one used in this study, maximization of ATP yield was used to model the uninfected cells. Maximization of ATP yield implies the maximization of the accumulation of ATP in the metabolic network. The biological rationale behind this approach is to account for the requirement of ATP in the processes that are not included in the metabolic model, such as mobility, futile cycles  and heat generation. The allowance of accumulation of ATP in the metabolic network would seem to be in contradiction with the basic assumption of steady state of metabolic networks, which implies that there is no accumulation of metabolites. This issue was easily reconciled by including an exchange reaction flux to represent ATP requirements not directly accounted for in the metabolic network as defined. Therefore, the objective function was defined as the maximization of this exchange reaction.
In silico reconstruction of genome-scale metabolic network of infected Escherichia coli
List of viral reactions and metabolites.
(35) ala-L + (29) arg-L + (17) asn-L + (13) asp-L + (3) cys-L + (18) gln-L + (17) glu-L + (28) gly + (5) his-L + (16) ile-L + (36) leu-L + (14) lys-L + (8) met-L + (16) phe-L + (17) pro-L + (36) ser-L + (27) thr-L + (12) trp-L + (16) tyr-L + (30) val-L → mat_v
(14) ala-L + (4) arg-L + (10) asn-L + (4) asp-L + (2) cys-L + (6) gln-L + (5) glu-L + (9) gly + (8) ile-L + (7) leu-L + (6) lys-L + (3) met-L + (4) phe-L + (6) pro-L + (13) ser-L + (9) thr-L + (2) trp-L + (4) tyr-L + (14) val-L → coat_v
(3) ala-L + (8) arg-L + (2) asn-L + asp-L + cys-L + (9) gln-L + (3) glu-L + his-L + (3) ile-L + (12) leu-L + (2) lys-L + met-L + (5) phe-L + (4) pro-L + (6) ser-L + (9) thr-L + (2) tyr-L + (3) val-L → lys_v
(41) ala-L + (42) arg-L + (18) asn-L + (33) asp-L + (7) cys-L + (15) gln-L + (22) glu-L + (44) gly + (10) his-L + (32) ile-L + (50) leu-L + (25) lys-L + (10) met-L + (28) phe-L + (27) pro- L + (48) ser-L + (31) thr-L + (9) trp-L + (21) tyr-L + (32) val- L → rep_v
(834) atp + (933) ctp + (927) gtp + (875) utp → posRNA_v
(875) atp + (927) ctp + (933) gtp + (834) utp → negRNA_v
viral maturation protein
viral coat protein
viral lysis protein
viral replicase protein
positive sense viral RNA genome
negative sense viral RNA genome
The regulated expression of viral genes was incorporated in the infected cell model. This was achieved by fixing the ratios of the fluxes of maturation protein production, replicase protein production and lysis protein production to coat protein production. The coat gene is expressed constitutively and forms the major fraction of the viral proteins that are produced [6, 18, 19]. The other viral proteins are produced in much smaller amounts due to the regulated expression of their corresponding genes . Previous studies have shown that the amounts of each of the maturation and replicase proteins were only 1/10th of that of the coat protein [18, 19]. As for the lysis protein, the amount was found to be only 1/25th of that of the coat protein . These ratios were used in the model to account for the regulated expression of the viral genes. The ratios of fluxes corresponding to viral maturation protein production, viral replicase protein production and viral lysis protein production to viral coat protein production were fixed to 1:10, 1:10 and 1:25, respectively.
Replication of the viral genome was also integrated into the infected cell model by including the positive and negative viral RNA strands as metabolites and the corresponding reactions that produce these metabolites. However, to constrain the fluxes of these reactions, there has been no experimental quantification of the number of positive and negative strands per cell. Hence, the fluxes of the production of these positive and negative strands were constrained based on the composition of MS2 virions. Each MS2 virion consists of 180 copies of the coat protein and one copy of the A protein and a positive strand of the viral genome . Therefore, the ratio of each of the fluxes of positive and negative viral genomic RNA production to the coat protein production was fixed to 1:180.
In addition to the application of constraints of the uninfected model, the exchange flux corresponding to the glucose uptake rate was constrained to the experimentally measured value of 12 mmol/g DW h. The growth rate and the oxygen uptake rate of the host were bound by the limits of 0 ≤ v i ≤ ∞ and -∞ ≤ v i ≤ ∞, respectively. Finally, appropriate objective functions were applied.
The most significant effect of the infection process is the change in the objective of the metabolic network. The objective function used to model the MS2 infected Escherichia coli cells was to maximize the viral reaction fluxes. The underlying reasoning for this approach was that the virus is a parasite and takes over the machinery of its host to produce viral components. Viral protein synthesis dominates over the host protein synthesis . Other objective functions considered were the simultaneous maximization of ATP yield and the viral fluxes, as well as the simultaneous maximization of growth rate and the viral fluxes. These objective functions resulted in no viral protein production and therefore, were not considered further in this study. This result reiterates the fact that host protein synthesis has to be eclipsed for viral protein synthesis to prevail.
The experimental system consisted of the bacteriophage MS2 (ATCC 15597-B1) and its host Escherichia coli C-3000 (ATCC 15597) in minimal M9 media  (supplemented with 2 mM CaCl2 and 10 μg/ml thiamine). The experimental methods used in this study were the same as described previously . In brief, for the uninfected cells, C-3000 was grown in minimal M9 media at 37°C (batch culture) and OD600 was measured using a Biomate 3 Spectrophotometer (Thermo Spectronic, Rochester, NY), glucose concentrations were measured using a YSI 2700 Select Biochemical Analyzer (YSI, Yellow Springs, OH) and dissolved oxygen was measured using a Dissolved Oxygen Meter SM600 (Milwaukee Instruments, Inc, Rocky Mount, NC).
For the infected cells, a C-3000 culture in minimal M9 media at 37°C (batch culture) was infected with MS2 phage at an OD600 of 0.8 and an MOI of 10. OD600 and glucose concentrations were measured. In addition, two kinds of plaque assays were carried out as described in previously published work . The first method was used to quantify the combined sum of infected cells and free phages and started with the mixing of 100 μl of an appropriately diluted sample of the culture with 3 ml E medium (10 g/l of tryptone, 1 g/l of yeast extract and 8 g/l of NaCl; also referred to as "EM") soft agar and 100 μl of the plating culture (an overnight culture of Escherichia coli C-3000). This mixture was poured on pre-warmed EM plates and incubated overnight at 37°C. Plaques were counted the next day. The second method was used to quantify free phages only and was similar to the above technique. The only differences were the use of Escherichia coli C-3000+pBAD102 (ampicillin resistant) as the plating culture and the addition of 100 g/ml ampicillin to EM soft agar and EM plates. The presence of ampicillin did not allow the infected cells to grow as they were ampicillin sensitive. As a result, the plaques that were formed were only due to the free phage. Utilizing both approaches allowed for determination of free phage density and infected cell density.
Determination of glucose and oxygen uptake rates of uninfected Escherichia coli C-3000 cells
where X is the cell density, S is the glucose concentration, μ is the growth rate of the cells and kglu is the glucose uptake rate of the cells.
where X is the cell density, μ is the growth rate of the cells, Soxy is the dissolved oxygen concentration, kLa is the mass transfer coefficient, is the saturated oxygen concentration in the medium and koxy is the oxygen uptake rate of the cells. The first term on the right hand side of (5) relates to the transfer of atmospheric oxygen into the medium.
Determination of glucose uptake rate of MS2 infected Escherichia coli C-3000 cells
where X is uninfected cell density, Z is phage-sensitive uninfected cell density, R is phage-resistant uninfected cell density, Y is phage-infected cell density, P is free phage density and S is substrate (glucose) concentration.
The Berkeley Madonna software package (version 8.3.14)  was used to estimate parameters in the determination of glucose and oxygen uptake rates of uninfected cells and glucose uptake rate of infected cells.
Results and Discussion
The growth rate of the uninfected Escherichia coli C-3000 in minimal M9 media was found to be 0.65 ± 0.02 h-1 (n = 5). The glucose uptake rate of the uninfected cells was 12 ± 0.5 mmol/(g DW h) (n = 3) and the oxygen uptake rate was 27 ± 1 mmol/(g DW h) (n = 2). The glucose uptake rate of the MS2-infected Escherichia coli C-3000 cells in minimal M9 media was the same as that of the uninfected cells, 12 ± 1 mmol/(g DW h) (n = 3). See Additional file 1 for the original data that was used to determine the above mentioned cell properties.
Infected Escherichia coli model
The results of the genome-scale model of the infected cell predicted the flux of the viral coat protein production to be 4.12 × 10-2 mmol/g DW h and the flux variability analysis showed no variation in this value. The experimentally calculated value of the flux, 1.05 × 10-5 mmol/g DW h, was based on the burst size of 10 ± 1 pfu/cell , latent period of 60 ± 2 min  and the composition of the viral capsid containing 180 copies of the coat protein. The discrepancy between the experimentally calculated value and the model prediction may be explained by the fact that the experimentally calculated value only accounts for the plaque forming units (pfu). It does not account for the defective MS2 viral particles that are also produced in host cells [24–27] and are incapable of forming plaques. Although, no literature was found regarding the particle to pfu ratio of MS2, studies have reported particle to pfu ratios of bacteriophage f2, very closely related to MS2, that range from 7 to 39 [28, 29]. In addition, the experimentally measured burst size assumed that all infected cells produced plaque forming units and an average number of progeny virus particles were calculated on a per cell basis. In contrast, studies specific to MS2 bacteriophage and the use of minimal media have reported that the fraction of infected cells producing plaque forming units (infective phages) is as low as 0.01 [24, 30, 31]. This suggests that the average burst size may further be increased by as much as two orders of magnitude. Based on these findings, overall, the burst size could be increased by three orders of magnitude increasing the experimentally calculated value of the coat protein production flux by three orders of magnitude. This may account for the differences between the experimentally measured value and the model prediction of coat protein production flux.
Comparison of the infected and uninfected cell metabolomes
A comparison of the major pathways in the infected cell metabolome was made with those in the uninfected cell metabolome. Flux variability analysis showed no variability in the fluxes of the pathways discussed below - both in the uninfected and the infected models.
The infected cell model predicted that the host cell does not grow. The value of the flux corresponding to the biomass production was zero, which is consistent with the experimental observations of MS2/C-3000 viral dynamics . It was also noted that all the cell wall synthesis and membrane lipid metabolism reactions were deactivated. This result also contributes to the fact that the biomass growth rate was zero and this weakening of the cell wall may further aid the virus in the lysis of the host at the end of the infection cycle.
List of bacterial reactions and metabolites.
Pentose Phosphate Pathway (Figure 2)
g6p + nadp → 6pgl + h + nadph
6pgl + h2o → 6pgc + h
6pgc + nadp → co2 + nadph + ru5p-D
ru5p-D → xu5p-D
ru5p-D → r5p
r5p + xu5p-D → g3p + s7p
e4p + xu5p-D → f6p + g3p
g3p + s7p → e4p + f6p
Glycolytic Pathway (Figure 3)
glc-D + pep → g6p + pyr
glc-D + atp → g6p + adp + h
g6p → f6p
atp + f6p → adp + fdp + h
fdp → dhap + g3p
g3p + nad + pi → 13dpg + h + nadh
13dpg + adp → 3pg + atp
3pg → 2pg
2pg → h2o + pep
TCA Cycle (Figure 4)
accoa + h2o + oaa → cit + coa + h
cit → acon-C + h2o
acon-C + h2o → icit
icit + nadp → akg + co2 + nadph
akg + coa + nad → co2 + nadh + succoa
adp + pi + succoa → atp + coa + succ
q8 + succ → fum + q8h2
fum + h2o → mal-L
mal-L + nad → h + nadh + oaa
Flavin adenine dinucleotide
Flavin adenine dinucleotide - reduced
Nicotinamide adenine dinucleotide
Nicotinamide adenine dinucleotide - reduced
Nicotinamide adenine dinucleotide phosphate
Nicotinamide adenine dinucleotide phosphate - reduced
The changes in the infected cell may be attributed to the up-regulation of biosynthesis of amino acids, as shown in Figure 5. The glucose uptake rates of uninfected and infected cells are similar. A change in the objective function removes the constraints on the usage of ATP that leads to the up-regulation of the hexokinase reaction and the down-regulation of the GLCPTS reaction. However, these changes do not affect the amounts of glucose-6-phosphate, which remains the same in both uninfected and infected cells. The changes start with the division of glucose-6-phosphate. The major portion of glucose-6-phosphate is diverted towards the pentose phosphate pathway, leading to the up-regulation of the pentose phosphate pathway. As a result more NADPH is produced that leads into the up-regulation of biosynthesis of amino acids. Only a minor portion of glucose-6-phosphate feeds the glycolytic pathway leading to the down-regulation of the PGI reaction. Other reactions downstream of the glycolytic pathway are not affected because the up-regulation of the pentose phosphate pathway increases the amounts of the by-products of the pentose phosphate pathway, fructose-6-phosphate and glyceraldehyde 3-phosphate, which are fed back into the glycolytic pathway.
A systems biological approach was taken to investigate the effects of viral infection on the host metabolism. It was seen that the viral infection alters the host metabolism to upregulate the biosynthesis of amino acids. Resources were diverted away from the TCA cycle and from the synthesis of cell wall components in order to upregulate the pentose phosphate pathway. This finally resulted in the production of viral proteins. Interestingly, a similar phenomenon has been noted in cancer cells where glucose is used primarily for anabolic processes via the pentose phosphate pathway [32, 33], suggesting the use of the E. coli/ MS2 system as a potential test bed for the study of drugs targeting metabolism for the treatment of cancer. Clearly, however, any such approach would require far more rigorous study of the commonalities between the viral infection process and cancer metabolism before any real parallels can be drawn.
This study is a proof of principle where techniques of modeling the metabolic system of a virus-infected host were described. Incorporation of regulated viral gene expression, viral genome replication and a suitable objective function to model an infected cell were shown. The greatest advantage of this approach is that the effects can be seen at full genome-scales. Advances in both experimental and computational fields could be used to extend this work to higher order organisms. Efforts to model human cellular metabolism are taking place [34–36] and application of the methods described in the current study may give new insights. Information based on the up-regulation and down-regulation of certain reactions could then be utilized to draw new treatment strategies.
Support for this work was provided by NIH Grant 1R03LM009753-01.
- Ikeda M, Kato N: Modulation of host metabolism as a target of new antivirals. Advanced Drug Delivery Reviews. 2007, 59: 1277-1289. 10.1016/j.addr.2007.03.021View ArticlePubMedGoogle Scholar
- Price ND, Papin JA, Schilling CH, Palsson BO: Genome-scale microbial in silico models: the constraints-based approach. Trends Biotechnol. 2003, 21: 162-169. 10.1016/S0167-7799(03)00030-1View ArticlePubMedGoogle Scholar
- Yeh I, Hanekamp T, Tsoka S, Karp PD, Altman RB: Computational analysis of Plasmodium falciparum metabolism: organizing genomic information to facilitate drug discovery. Genome Res. 2004, 14: 917-924. 10.1101/gr.2050304PubMed CentralView ArticlePubMedGoogle Scholar
- Jamshidi N, Palsson BO: Investigating the metabolic capabilities of Mycobacterium tuberculosis H37Rv using the in silico strain iNJ661 and proposing alternative drug targets. BMC Syst Biol. 2007, 1: 26- 10.1186/1752-0509-1-26PubMed CentralView ArticlePubMedGoogle Scholar
- Schilling CH, Palsson BO: Assessment of the metabolic capabilities of Haemophilus influenzae Rd through a genome-scale pathway analysis. J Theor Biol. 2000, 203: 249-283. 10.1006/jtbi.2000.1088View ArticlePubMedGoogle Scholar
- van Duin J, Tsareva N: The Single-Stranded RNA Bacteriophages. 1988, New York: Plenum Press, 1Google Scholar
- Berzin V, Rosenthal G, Gren EJ: Cellular macromolecule synthesis in Escherichia coli infected with bacteriophage MS2. Eur J Biochem. 1974, 45: 233-242. 10.1111/j.1432-1033.1974.tb03547.xView ArticlePubMedGoogle Scholar
- Feist AM, Henry CS, Reed JL, Krummenacker M, Joyce AR, Karp PD, Broadbelt LJ, Hatzimanikatis V, Palsson BO: A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information. Mol Syst Biol. 2007, 3: 121- 10.1038/msb4100155PubMed CentralView ArticlePubMedGoogle Scholar
- GNU Linear Programming Kit. http://www.gnu.org/software/glpk/
- Mahadevan R, Schilling CH: The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. Metab Eng. 2003, 5: 264-276. 10.1016/j.ymben.2003.09.002View ArticlePubMedGoogle Scholar
- Reed JL, Vo TD, Schilling CH, Palsson BO: An expanded genome-scale model of Escherichia coli K-12 (iJR904 GSM/GPR). Genome Biol. 2003, 4: R54- 10.1186/gb-2003-4-9-r54PubMed CentralView ArticlePubMedGoogle Scholar
- Edwards JS, Covert M, Palsson B: Metabolic modelling of microbes: the flux-balance approach. Environ Microbiol. 2002, 4: 133-140. 10.1046/j.1462-2920.2002.00282.xView ArticlePubMedGoogle Scholar
- Knorr AL, Jain R, Srivastava R: Bayesian-based selection of metabolic objective functions. Bioinformatics. 2007, 23: 351-357. 10.1093/bioinformatics/btl619View ArticlePubMedGoogle Scholar
- Schuetz R, Kuepfer L, Sauer U: Systematic evaluation of objective functions for predicting intracellular fluxes in Escherichia coli. Mol Syst Biol. 2007, 3: 119- 10.1038/msb4100162PubMed CentralView ArticlePubMedGoogle Scholar
- Schuster S, Pfeiffer T, Fell DA: Is maximization of molar yield in metabolic networks favoured by evolution?. J Theor Biol. 2008, 252: 497-504. 10.1016/j.jtbi.2007.12.008View ArticlePubMedGoogle Scholar
- Neijssel OM, Buurman ET, Teixeira de Mattos MJ: The role of futile cycles in the energetics of bacterial growth. Biochim Biophys Acta. 1990, 1018: 252-255. 10.1016/0005-2728(90)90260-BView ArticlePubMedGoogle Scholar
- Birge EA: Bacterial and Bacteriophage Genetics. 1994, New York: Springer-Verlag, 3View ArticleGoogle Scholar
- Nathans D, Oeschger MP, Eggen K, Shimura Y: Bacteriophage-Specific Proteins in E. Coli Infected with an Rna Bacteriophage. Proc Natl Acad Sci USA. 1966, 56: 1844-1851. 10.1073/pnas.56.6.1844PubMed CentralView ArticlePubMedGoogle Scholar
- Nathans D, Oeschger MP, Polmar SK, Eggen K: Regulation of protein synthesis directed by coliphage MS2 RNA. I. Phage protein and RNA synthesis in cells infected with suppressible mutants. J Mol Biol. 1969, 39: 279-292. 10.1016/0022-2836(69)90317-9View ArticlePubMedGoogle Scholar
- Berkhout B, Schmidt BF, van Strien A, van Boom J, van Westrenen J, van Duin J: Lysis gene of bacteriophage MS2 is activated by translation termination at the overlapping coat gene. J Mol Biol. 1987, 195: 517-524. 10.1016/0022-2836(87)90180-XView ArticlePubMedGoogle Scholar
- Sambrook J, Russell DW: Molecular Cloning: A Laboratory Manual. 2001, Cold Spring Harbor, New York: Cold Spring Harbor Laboratory Press, 3Google Scholar
- Jain R, Knorr AL, Bernacki J, Srivastava R: Investigation of bacteriophage MS2 viral dynamics using model discrimination analysis and the implications for phage therapy. Biotechnol Prog. 2006, 22: 1650-1658.View ArticlePubMedGoogle Scholar
- Berkeley Madonna. http://www.berkeleymadonna.com/
- Ricciuti CP, Haywood AM: Production of defective MS2 virions in resting Escherichia coli. Virology. 1974, 58: 164-175. 10.1016/0042-6822(74)90151-2View ArticlePubMedGoogle Scholar
- Renhof RF, Veselova SV, Jansone IV, Gren EJ: Specific RNA fragments in defective particles of bacteriophage MS2. FEBS Lett. 1974, 43: 89-92. 10.1016/0014-5793(74)81112-9View ArticlePubMedGoogle Scholar
- Rohrmann GF, Krueger RG: Physical, biochemical, and immunological properties of coliphage MS-2 particles. J Virol. 1970, 6: 269-279.PubMed CentralPubMedGoogle Scholar
- Shiba T, Saigo K: Packaging of tRNA and 6 S stable RNA into bacteriophage MS2 particles produced in Escherichia coli treated with 5-fluorouracil. Virology. 1982, 119: 209-213. 10.1016/0042-6822(82)90079-4View ArticlePubMedGoogle Scholar
- Lodish HF, Horiuchi K, Zinder ND: Mutants of the bacteriophage f2. V. On the production of noninfectious phage particles. Virology. 1965, 27: 139-155. 10.1016/0042-6822(65)90154-6View ArticlePubMedGoogle Scholar
- Cooper S, Zinder ND: The Growth of an RNA Bacteriophage: The Role of Protein Synthesis. Virology. 1963, 20: 605-612. 10.1016/0042-6822(63)90285-XView ArticlePubMedGoogle Scholar
- Davis JE, Sinsheimer RL: The replication of bacteriophage MS2. 1. Transfer of parental nucleic acid to progeny phage. J Mol Biol. 1963, 6: 203-207. 10.1016/S0022-2836(63)80069-8View ArticlePubMedGoogle Scholar
- Lerner FL: The population biology of male-specific bacteriophages: existence conditions. 1985, University of Massachusetts, ZoologyGoogle Scholar
- Boros LG, Cascante M, Lee WN: Metabolic profiling of cell growth and death in cancer: applications in drug discovery. Drug Discov Today. 2002, 7: 364-372. 10.1016/S1359-6446(02)02179-7View ArticlePubMedGoogle Scholar
- Cascante M, Boros LG, Comin-Anduix B, de Atauri P, Centelles JJ, Lee PW: Metabolic control analysis in drug discovery and disease. Nat Biotechnol. 2002, 20: 243-249. 10.1038/nbt0302-243View ArticlePubMedGoogle Scholar
- Duarte NC, Becker SA, Jamshidi N, Thiele I, Mo ML, Vo TD, Srivas R, Palsson BO: Global reconstruction of the human metabolic network based on genomic and bibliomic data. Proc Natl Acad Sci USA. 2007, 104: 1777-1782. 10.1073/pnas.0610772104PubMed CentralView ArticlePubMedGoogle Scholar
- Mo ML, Jamshidi N, Palsson BO: A genome-scale, constraint-based approach to systems biology of human metabolism. Mol Biosyst. 2007, 3: 598-603. 10.1039/b705597hView ArticlePubMedGoogle Scholar
- Alter O: Genomic signal processing: from matrix algebra to genetic networks. Methods Mol Biol. 2007, 377: 17-60. full_textView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.