Quantitative relationships for specific growth rates and macromolecular compositions of Mycobacterium tuberculosis, Streptomyces coelicolor A3(2) and Escherichia coli B/r: an integrative theoretical approach

Robert A. Cox

Division of Mycobacterial Research, National Institute for Medical Research, London NW7 1AA, UK

Correspondence
Robert A. Cox
rcox{at}nimr.mrc.ac.uk


   ABSTRACT
TOP
ABSTRACT
INTRODUCTION
THEORETICAL ANALYSES
RESULTS AND DISCUSSION
APPENDIX
REFERENCES
 
Further understanding of the physiological states of Mycobacterium tuberculosis and other mycobacteria was sought through comparisons with the genomic properties and macromolecular compositions of Streptomyces coelicolor A3(2), grown at 30 °C, and Escherichia coli B/r, grown at 37 °C. A frame of reference was established based on quantitative relationships observed between specific growth rates (µ) of cells and their macromolecular compositions. The concept of a schematic cell based on transcription/translation coupling, average genes and average proteins was developed to provide an instantaneous view of macromolecular synthesis carried out by cells growing at their maximum rate. It was inferred that the ultra-fast growth of E. coli results from its ability to increase the average number of rRNA (rrn) operons per cell through polyploidy, thereby increasing its capacity for ribosome synthesis. The maximum growth rate of E. coli was deduced to be limited by the rate of uptake and consumption of nutrients providing energy. Three characteristic properties of S. coelicolor A3(2) growing optimally (µ=0·30 h–1) were identified. First, the rate of DNA replication was found to approach the rate reported for E. coli (µ=1·73 h–1); secondly, all rrn operons were calculated to be fully engaged in precursor-rRNA synthesis; thirdly, compared with E. coli, protein synthesis was found to depend on higher concentrations of ribosomes and lower concentrations of aminoacyl-tRNA and EF-Tu. An equation was derived for E. coli B/r relating µ to the number of rrn operons per genome. Values of µ=0·69 h–1 and µ=1·00 h–1 were obtained respectively for cells with one or two rrn operons per genome. Using the author's equation relating the number of rrn operons per genome to maximum growth rate, it is expected that M. tuberculosis with one rrn operon should be capable of growing much faster than it actually does. Therefore, it is suggested that the high number of insertion sequences in this species attenuates growth rate to still lower values.


Abbreviations: RNAP, RNA polymerase


   INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
THEORETICAL ANALYSES
RESULTS AND DISCUSSION
APPENDIX
REFERENCES
 
Mycobacteria are rod-shaped, Gram-positive, acid-fast, aerobic and non-motile (Wayne & Kubica, 1986). Pathogenic mycobacteria such as Mycobacterium leprae and Mycobacterium tuberculosis have notable properties including slow growth [generation times of 12 days (Shepard, 1960) and 1 day (Wayne, 1994) respectively], an ability to survive within host cells (Armstrong & D'Arcy-Hart, 1971; Ferrari et al., 1999) and an ability to persist in a dormant state (Wayne & Hayes, 1996). Other mycobacteria are pathogenic to man or animals (Wayne & Kubica, 1986; Woods & Washington, 1987). Traditionally, mycobacteria are divided into slow-growers and fast-growers according to the time taken for colonies to appear on a solid medium (Wayne & Kubica, 1986). In general slow-growers were found to have a single rRNA (rrn) operon whereas fast-growers were found to have two rrn operons per genome (Domenech et al., 1994).

Knowledge of mycobacterial physiology has been limited by the slow growth of the major pathogens in the laboratory and other technical problems such as cell aggregation. The availability of the sequence of the genome of M. tuberculosis (Cole et al., 1998) and limited data for the macromolecular composition of its very close relative Mycobacterium bovis bacille Calmette–Guérin (BCG) (Winder & Rooney, 1970) provide the basis for gaining insight into the growth of the tubercle bacillus and other mycobacteria through a comparative approach. In this study a single quantitative framework was sought to relate genomic and macromolecular properties to the rates of protein synthesis of M. bovis BCG and two other representative bacterial species, namely Streptomyces coelicolor A3(2) and Escherichia coli B/r. The three species range in their maximum specific growth rates (µmax) from 0·029 h–1 (M. bovis BCG) to 1·73 h–1 (E. coli).

S. coelicolor A3(2), which is Gram-positive, soil-dwelling and filamentous, is an aerial-mycelium-producing actinomycete (Shahab et al., 1996); the genome, which is linear rather than circular, is almost twice the size of that of M. tuberculosis. E. coli belongs to a group of organisms, Enterobacteriaceae, that are Gram-negative, rod-shaped and capable of ultra-fast growth (Neidhardt, 1996). E. coli was defined as an ultra-fast grower (Cox, 2003) because, when it is growing at its maximum rate the generation time, tD, is less than the time, C, needed to replicate the genome (C>tD), leading to new-born cells possessing more than one genome equivalent. In contrast, M. tuberculosis (a slow-grower) and S. coelicolor (a fast-grower) are thought to replicate their genomes once only during the cell division cycle.

The concept of a virtual or schematic cell was developed in this study to provide an instantaneous view of macromolecular synthesis, particularly protein synthesis, carried out by cells growing at their maximum rate. The quantitative approach not only provides a succinct way of summarizing a wide range of data but also provides models amenable to mathematical manipulation. This integrative approach was used to calculate the specific growth rate of hypothetical E. coli cells having one or two rRNA (rrn) operons per genome as models for mycobacteria and to explore factors limiting the maximum specific growth rate of the three species studied.


   THEORETICAL ANALYSES
TOP
ABSTRACT
INTRODUCTION
THEORETICAL ANALYSES
RESULTS AND DISCUSSION
APPENDIX
REFERENCES
 
Bacterial strains used for the theoretical analyses.
Literature data of Escherichia coli B/r, Streptomyces coelicolor A3(2), Mycobacterium tuberculosis H37Rv and Mycobacterium bovis BCG provided the basis for these theoretical analyses. The genomic properties of these strains are summarized in Table 1.


View this table:
[in this window]
[in a new window]
 
1 Comparisons of features of genomes

 
Definitions, axioms, principles and concepts
The symbols for the variables used in the development of a common theoretical framework are summarized in Table 2.


View this table:
[in this window]
[in a new window]
 
2 Definitions of variables

 
The concept of the average cell.
Schaechter et al. (1958) proposed that for exponentially growing cells an average mass (mx(av)) of constituent, x, per cell may be defined as the amount of x per culture divided by the number of cells per culture. The constituents considered in this study include dry cell mass (mdc(av)), protein (mp(av)), RNA (mRNA(av)) and DNA (mDNA(av)). Average cell volume (v(av)) may be calculated on the basis of the assumptions (Cox, 2003) that 70 % of a cell's mass is water and that a cell's buoyant density is close to 1·09 g cell mass (ml cell volume)–1 [or pg fl–1, or pg µm–3]. It is further assumed that mp(av) constitutes a fixed proportion (40–50 %, according to species) of mdc(av) (Bremer & Dennis, 1996; Cox, 2003). Thus, concentrations of a constituent, for example RNA, may be expressed as mRNA(av)/mp(av), mRNA(av)/mdc(av) or mRNA(av)/v(av).

The concept of the minimal cell.
As the supply of nutrients become scarce E. coli cells become smaller (Bremer & Dennis, 1996) and eventually attain a minimal size and a characteristic macromolecular composition (Jacobsen, 1974, cited by Ingraham et al., 1983). It is proposed that minimal cells are a general feature of bacterial growth. By using the superscript {ddagger}, I indicate the value of the corresponding variable (, etc.) for the minimal cell.

Significance of properties of cells growing at their maximum rate.
The maximum specific growth rate, µmax, is attained when cells are amply supplied with the most favourable nutrients. An asterisk (mp(av)*, etc.) is used to denote a property of a cell growing at its maximum rate. A cell's capacity for growth is defined by the properties of minimal cells and of cells growing at their maximum rate.

The concept of a virtual schematic cell.
The formulation is based on the notions of an average gene, an average mRNA and an average protein (see Table 1) and on a few key equations derived from data for macromolecular compositions of cells growing at different rates. The model cell provides an instantaneous view of macromolecular synthesis, particularly protein synthesis, carried out by an average cell growing at its maximum rate.

Quantitative aspects of cell processes
DNA synthesis.
Except for ultra-fast growth, tD is equal to the sum of periods relating to DNA synthesis [see equation (1)] where B is the period between cell division and the start of DNA replication and D is the period following the completion of DNA synthesis and cell division:

No biochemical significance has been attributed to B, which diminishes as µ increases so that B->0 (Helmstetter, 1996). The D period is thought to be required to allow chromosome segregation and septum formation to take place.

The rate of elongation of DNA per replication fork ({varepsilon}DNA, bp per fork s–1), the size of the genome (lg, bp) and C (h) are, by definition related by equation (2):

The number of genome equivalents per average cell (ng(av)) is given by equation (3) (Helmstetter & Cooper, 1968):

When B->0, (C+D)/tD->1 and D/tD->0·25 then ng->1·6 genome equivalents per average cell, the maximum value for a fast-growing cell.

Transcription/translation coupling.
Coupling between the processes of bacterial transcription and translation has long been accepted (Stent, 1964; Byrne et al., 1964). This coupling was vividly shown by electron microscopy (Miller et al., 1970) because ribosomes were found to be attached to nascent transcripts, demonstrating that translation accompanies transcription. The same technique revealed that, at any instant, the chromosome of an individual E. coli cell is largely transcriptionally inactive and that few, if any, free polyribosomes are found in the cytoplasm. Although the chromosome is largely transcriptionally inactive, a high proportion of its genes are expressed during the lifetime of the cell, as shown by proteome (Tonella et al., 1988) and transcriptome analysis (Bernstein et al., 2002). Transcription/translation coupling is thought to require that the rates of mRNA elongation ({varepsilon}mRNA, nucleotides h–1) and peptide chain elongation ({varepsilon}aa, amino acid residues h–1) are co-ordinated (Bremer & Dennis, 1996); see equation (4), where the factor 3 reflects the number of nucleotides per codon:

Implicit in equation (4) is the notion that the rates of initiation of transcription and translation are co-ordinated. The process of translation is thought to protect mRNA from degradation by RNases (degradosomes) (Grunberg-Manago, 1999; Régnier & Arraiano, 2000; Leroy et al., 2002; Khodursky & Bernstein, 2003).

If the amount of mRNA is represented by nm·nuc(av) nucleotides, and nR(av) is the number of ribosomes and {beta}R the fraction of ribosomes actively synthesizing proteins, then the ratio nm·nuc(av)/{beta}R·nR(av)={sigma} nucleotides per ribosome reflects the organization of an average polyribosome. The minimum value of {sigma} reflects the size of the ‘footprint’ of an initiation form of RNA polymerase (RNAP) bound to a promoter, which is approximately 80 bp (Krummel & Chamberlin, 1989). Accordingly, it is assumed that there is a minimum of 80 nucleotides of mRNA per ribosome. Neglecting secondary structure, this segment corresponds to a length of approximately 27 nm (270 Å), which is comparable to the diameter of a ribosome of 25 nm (Noller & Nomura, 1996). The maximum value of {sigma} corresponds to the length of the region protected by ribosomes plus a region smaller than the number of nucleotides that form a binding site for a degradosome. An arbitrary value of 100 nucleotides (34 nm) was assumed for protected mRNA; that is, unprotected stretches of mRNA of 20 nucleotides or more were considered to be capable of binding degradosomes and, hence, to be very rapidly degraded. Kinetic studies summarized by Bremer & Dennis (1996) led to estimates for {sigma} of 55–80 nucleotides of mRNA per ribosome depending on the growth rate of E. coli (see Table 3).


View this table:
[in this window]
[in a new window]
 
Table 3. Macromolecular properties of ‘average’ cells of E. coli B/r

 
The rate of ribosome synthesis.
It is thought that rRNA synthesis is the rate-limiting step in ribosome synthesis (Bremer & Dennis, 1996). Thus, the number, nrrn(av), of rrn operons per cell is a significant factor in cell proliferation. During ultra-fast growth, both nrrn(av) and ng(av) increase as the specific growth rate increases (see Table 3). The number of rrn operons per average cell is the product of three factors [see equation (5)]; namely nrrn(g), the number of rrn operons per genome, ng(av), the number of genome equivalents per average cell, and f, a constant that reflects the locations of rrn operons within the genome with respect to the origin of replication. For example, f>1·0 when most rrn operons are clustered near to the origin of replication. In this case the average number of rrn operons per cell (nrrn(av)) is underestimated by the average number of genomes per cell.

The specific ribosome synthesis rate, {omega}R(av)(number of ribosomes synthesized per cell h–1), is defined by equation (6a):

leading to equation (6b) because, by definition, {omega}R(av) ribosomes h–1and {omega}rRNA(av) nucleotides h–1/lrrn are equal.

The term iRNAP is the average number of RNAP holoenzymes transcribing an rrn operon, {varepsilon}rRNA(nucleotides) is the elongation rate of nascent precursor-rRNA (pre-rRNA) and lrrn is the length (nucleotides) of a complete transcript of an rrn operon. {varepsilon}rRNA is thought to be approximately equal to 2{varepsilon}mRNA or 6{varepsilon}aa(Bremer & Dennis, 1996), possibly reflecting feedback control of rRNA synthesis. The ratio {omega}R(av)/nrrn(av), which describes the rate of transcription (iRNAP·{varepsilon}rRNA) of a single rrn operon [see equation (6b)], forms the basis of a method for assessing the relationship between µ and nrrn(av)[see Appendix]. The limiting value, ilim, of iRNAP is observed when the rate of initiation is equal to the rate of promoter clearance; then a loading of one RNAP per 80 nucleotides (the footprint of the RNAP initiating complex) would be expected; that is ilim is equal to lrrn/80, or 70–75 RNAPs per rrn operon. The ratio iRNAP/ilim provides a measure of the transcriptional activity of each rrn operon. The minimum number, , required to meet the observed specific ribosome synthesis rate is given by equation (6c); {varepsilon}rRNA is equated with 6 {varepsilon}aa residues h–1.

Equation (7) was derived by relating the specific growth rate, µx, to the average number, xrrn(av), of rrn operons per average cell; {gamma}F is a constant dependent on the growth medium (see Appendix).

Protein synthesis.
The rate-limiting step in protein synthesis is the interaction of a ternary complex of aminoacyl-tRNA, EF-Tu and GTP with the A site of the ribosome. The concentrations of aminoacyl-tRNA and EF-Tu are conveniently expressed as the number of copies of aminoacyl-tRNA per ribosome (naat/R) and the number of copies of EF-Tu per ribosome (nEF/R). It is assumed that both naat/R and nEF/R have characteristic values for a particular species during the growth cycle, irrespective of the specific growth rate (Cox, 2003).

The key equations for describing protein synthesis during exponential growth concern the specific protein synthesis rate, {omega}p(av) (fg protein synthesized per cell h–1), which is defined by equation (8a):

or the equivalent equation (8b), where {omega}aa(av) is the number of amino acids incorporated into protein per cell h–1:

Alternatively, {omega}aa(av) may be related to the number (nR(av)) of ribosomes per cell and the peptide chain elongation rate {varepsilon}aa amino acids incorporated per ribosome h–1 (Bremer & Dennis, 1996); see equation (9), where {beta}R is the fraction of ribosomes actively synthesizing protein:

Equation (10a) may be derived by equating the right-hand sides of equations (8b) and (9) and rearranging to make µ the subject (Bremer & Dennis, 1996):

The term nR(av)/naa(av) is an expression of RNA concentration. Equation (10b) is an alternative form of equation (10a) where {phi}=1x10–4 for E. coli (see Appendix):

Although the terms that define {phi} are species specific, variations in the values of these terms between species are unlikely to be large and so {phi} is expected to remain close to the numerical value of 1x10–4. It was reported for E. coli (Cox, 2003) that specific growth rates were related to RNA concentration by equation (11), where a (the slope) and b (the intercept on the µ axis) are constants.

The ratio b/a is equal to the RNA concentration of minimal cells (); by definition ; . Equation (11) is related to equation (10b) through {phi} and {varepsilon}aa*.

Previously, it was shown (Cox, 2003) for E. coli that {omega}p(av) or {omega}aa(av) is proportional to the third power of the RNA concentration [see, for example, equation (12)] for values of (mRNA(av)/mdc(av))>=b/a.

The slope {theta} is considered to be proportional to the product of kinetic constants, the number of copies of aminoacyl-tRNA per ribosome and the number of copies of EF-Tu per ribosome. The parameters determining {theta} are also thought to be reflected in {varepsilon}aa.

Rate of energy consumption.
Cell proliferation requires the uptake and consumption of energy. Approximately half of the energy is used for the synthesis of macromolecules, of which 90 % or more is used for protein synthesis (Ingraham et al., 1983). The formation of each peptide bond needs the consumption of 4·2 high-energy phosphate bonds. Thus, the rate of energy consumption, {varepsilon}E ATP equivalents h–1, is estimated in terms of ATP equivalents by equation (13):


   RESULTS AND DISCUSSION
TOP
ABSTRACT
INTRODUCTION
THEORETICAL ANALYSES
RESULTS AND DISCUSSION
APPENDIX
REFERENCES
 
The genomic properties of E. coli B/r, S. coelicolor A3(2) and M. bovis BCG (see Table 1) reveal that in each case the average ORF is close to 1000 bp in length, encoding an average protein of approximately 330 amino acid residues. It is noted that the size of the genome of S. coelicolor is almost twice the size of the two other species. The sizes of the rRNA components and the RNA (rrn) operons are very similar.

The macromolecular properties of E. coli B/r grown at 37 °C, in five different media (media A–E), S. coelicolor A3(2) grown at different specific growth rates (µ=0·024 h–1 to µ=0·300 h–1) at 30 °C and M. bovis BCG (µ=0·029 h–1) grown at 37 °C are presented respectively in Table 3 (based on the review of Bremer & Dennis, 1996), Table 4 (based on Shahab et al., 1996) and Table 5 (based on Winder & Rooney, 1970; Cox, 2003).


View this table:
[in this window]
[in a new window]
 
Table 4. Macromoleclar properties of ‘average’ cells of S. coelicolor A3(2) grown at 30 °C

Based on data of Shahab et al. (1996) recalculated using the genome size of S. coelicolor A3(2) reported by Bentley et al. (2002).

 

View this table:
[in this window]
[in a new window]
 
Table 5. Macromolecular properties of ‘average’ cells of M. bovis BCG grown at 37 °C

Data were calculated (Cox, 2003) from the data of Winder & Rooney (1970).

 
Macromolecular compositions and properties of E. coli B/r grown at 37 °C
E. coli grown in medium A was found to replicate its genome once only during the growth cycle (tD>C+D). Cells grown in medium B (tD=1·00 h) have a composition close to that of the transition point (tD=C+D) beyond which ultra-fast growth (C>tD) takes place. The duration of C was found to vary from 1·12 h (µ=0·42 h–1) to 0·70 h (µ=1·73 h–1). The maximum value of {varepsilon}DNA/3600 s h–1 was found to be 933 bp per replication fork s–1.

The average time, trrn, for the transcription of an rrn operon was found to be 1·2 min by direct measurement, although the relation {varepsilon}rRNA=6{varepsilon}aa [see equation (4)] suggests that {varepsilon}rRNA should be dependent on growth rate so that trrn would be expected to decrease from 1·4 min in medium A to 0·8 min in medium E. Irrespective of the specific growth rate, the complement, nrrn(av), of rrn operons was estimated to be in excess of the minimum number, , needed to achieve the observed rate of rrn transcription. In other words, the specific rRNA synthesis rate could have been achieved by using fewer rrn operons more intensively. For discussion, see equations (6b) and (6c) of the Theoretical Analyses section and equations (6A) and (7A) of the Appendix. As shown in Table 3, at high specific growth rates (µ=1·73 h–1) {varepsilon}rRNA was found to be suboptimal ({varepsilon}rRNA<6{varepsilon}aa residues h–1), values of iRNAP/ilim were found to range from 0·06 (µ=0·42 h–1) to 0·90 (µ=1·73 h–1), and [see equation (6c)] was found to range from 0·9 operons (µ=0·42 h–1) to 22·0 operons (µ=1·73 h–1). The fraction , which measures the relative activity of each rrn operon, was found to range from 0·07 (µ=0·42 h–1) to 0·61 (µ=1·73 h–1). Thus, E. coli growing at its maximum rate (µ=1·73 h–1) has a complement of rrn operons in excess of its needs for pre-rRNA synthesis. It was estimated (see Appendix) that the maximum specific growth rates attainable for hypothetical strains with one, two, three and four rrn operons per genome were, respectively, 0·42 h–1, 1·00 h–1, 1·40 h–1 and 1·70 h–1. These conclusions are in accord with the specific growth rates observed on the progressive inactivation of rrn operons present in the E. coli genome (Condon et al., 1993, 1995). It was found that the specific growth rate was maintained by increasing the activities (both {varepsilon}rRNA and iRNAP/ilim) of the remaining operons when the number of functional operons per genome was reduced from seven to four. A reduction from seven to three operons led to a slower growth rate and to smaller cells, in agreement with the estimated value for three rrn operons per genome mentioned above. The unique property of ultra-fast growers is that, beyond the transition point from fast to ultra-fast growth (tD=C*+D*) the number of genome equivalents per cell, g(av), increases to two or more, with a concomitant increase in the number, nrrn(av), of rrn operons per cell [see equation (14)]:

The rate of ribosome synthesis needed to support the maximum specific growth rate of 1·73 h–1 could not be achieved without increasing nrrn(av) from 12·4 to a minimum of 22 operons through polyploidy.

Both the dry cell mass and volume of E. coli were found to increase 4·5-fold as µ increased from 0·42 h–1 to 1·73 h–1 (see Table 3). Empirically, linear plots were obtained by plotting µ2 against mdc(av) or v(av) (see Fig. 1), which are summarized by equations (15a) and (15b):




View larger version (10K):
[in this window]
[in a new window]
 
Fig. 1. Correlations between specific growth rate squared and average dry cell mass. I, {circ}, E. coli (see Table 3); II, {bullet}, S. coelicolor (see Table 4).

 
In equation (15a) the fraction, intercept/slope (q), is equal to the dry cell mass () of minimal cells of E. coli B/r; namely 180 fg, of which 89 fg is protein.

By definition, the slope, q=4·2x10–3, in equation (15a) is equal to , which is a measure of growth­directed metabolic activity of the cell. The relationship between dry cell mass and parameters such as the number of rrn operons per cell (nrrn(av)) and the peptide chain elongation rate ({varepsilon}aa) was made explicit [see equation (16)] by equating the right­hand sides of equations (7) and (15a) and rearranging:

The term {gamma}F relates to a particular growth medium (see Appendix).

The RNA concentration (mRNA(av)/mdc(av)) reflects the ribosome concentration needed to achieve a particular value of µ. The relation between µ and the product of (mRNA(av)/mdc(av)) and {varepsilon}aa is illustrated in Fig. 2 and is described by equation (10b), where {phi}=1x10–4, in accord with the theoretical value. The relations between µ and (mRNA(av)/mdc(av)), and between µ and mRNA(av)/v(av), which are illustrated in Fig. 3 are described by equations (17a) and (17b):




View larger version (14K):
[in this window]
[in a new window]
 
Fig. 2. Correlations between specific growth rate and the product of RNA concentration and the peptide chain elongation rate [see equation (10b)]. {circ}, E. coli (see Table 3); {bullet}, S. coelicolor (see Table 4).

 


View larger version (16K):
[in this window]
[in a new window]
 
Fig. 3. Correlations between specific growth rate and concentration of RNA (mRNA(av)/mdc(av)) or (mRNA(av)/v(av)). I, {circ}, E. coli (see Table 3). II, {bullet}, S. coelicolor (see Table 4). {blacksquare}, M. bovis (see Table 5).

 
As discussed in the Theoretical Analyses section the fractions 0·63/10·3 (0·06 fg RNA/fg dry cell mass) and 0·63/0·03 (21·0 fg RNA/fl) define the RNA concentration of a minimal cell. The correlations between specific protein synthesis rates ({omega}p(av)) and RNA concentrations (mRNA(av)/mdc(av) and mRNA(av)/v(av)), which are presented in Fig. 4, are summarized by equations (18a) and (18b):




View larger version (17K):
[in this window]
[in a new window]
 
Fig. 4. Correlations between specific protein synthesis rate and the third power of RNA concentration [see equation (12)]. (a) I, {circ}, E. coli (see Table 3); {blacksquare}, M. bovis BCG (see Table 5). (b) II, {bullet}, S. coelicolor (see Table 4).

 
Several genomic properties and metabolic activities of E. coli growing optimally in medium E (µ=1·73 h–1) are summarized in Fig. 5(a). The fibril diagram illustrates both the role of coupled transcription/translation and the notion (see Table 1) of an average gene of 950 bp encoding an average protein of 317 amino acid residues. One factor which may limit µmax is a limit to the rate of uptake of the nutrient source of carbon and energy, its active transport across the cytoplasmic membrane and metabolism. E. coli utilizes glucose very efficiently but in order to achieve specific growth rates in excess of 1·04 h–1 it is necessary to supplement a glucose medium by the addition of energy­saving intermediate compounds such as amino acids. The need to supplement a glucose medium in order to achieve µmax=1·73 h–1 suggests that there is a limit to the rate at which E. coli can generate energy, and that this limit defines µmax. The increase in energy generation and consumption needed to support a growth rate greater than 1·73 h–1 (tD=0·40 h) is illustrated by the following example. The properties of a hypothetical cell of µ=2·1 h–1 (tD=0·33 h) were calculated as follows: ng(av)=4·9 genome equivalents [see equation (3)]; nrrn(av)=47·1 operons [see equation (14)]; mdc(av)=1200 fg [see equation (15a)] and hence mp(av)=600 fg; {omega}p(av)=1270 fg h–1 [see equation (8a)], leading to {varepsilon}E/3600 s h–1=1·8x107 ATP equivalents s–1 [see equation (13)]. Thus, a 20 % increase in µ, from 1·73 h–1 to 2·10 h–1, was calculated to require a 60 % increase in the rate of energy consumption from 1·1x107 ATP equivalents s–1 (see Table 3) to 1·8x107 ATP equivalents s–1.



View larger version (35K):
[in this window]
[in a new window]
 
Fig. 5. Virtual schematic cells summarizing the genomic properties and transcription/translation activities of ultra­fast­, fast­ and slow­growing bacteria proliferating at their maximum rates µmax. (a) E. coli B/r (µmax=1·73 h–1); (b) S. coelicolor A3(2) (µmax=0·30 h–1); (c) M. bovis BCG (µmax=0·029 h–1). (i) Summary of genomic properties and the number of genome equivalents per ‘average’ cell. The properties of the genome, namely the number of rrn operons and the number of ‘average’ ORFs (see Table 1), are presented within the square brackets. The number of genome equivalents per ‘average’ cell is indicated outside the brackets on the lower right hand side. Striped bars represent the rrn operon; stippled bars represent the ‘average’ ORF. (ii) Transcription/translation activity of an ‘average’ ORF is represented schematically within the square brackets in the form of a fibril diagram. The size (bp) of an ‘average’ ORF is indicated and also the locations of 5' ends of nascent mRNA transcripts. Nascent polypeptide chains are not shown. The proportion of non­programmed ribosomes per ORF is indicated by the relative number of free ribosomes. The number of ‘average’ ORFs being transcribed/translated at any instant is shown by the number outside the square brackets on the lower right hand side. Stippled bars represent the ‘average’ ORF; vertical black bars along the ORF represent RNAP holoenzyme; black circles represent ribosomes; and the lines joining these circles represent nascent mRNA. (iii) Rate of synthesis of ‘average’ proteins. The rates of synthesis of ‘average’ proteins per fibril is indicated within the square brackets. The number of fibrils synthesizing protein is given outside the square brackets on the lower right­hand side.

 
Macromolecular compositions and properties of S. coelicolor A3(2) grown at 30 °C
The sizes of average genes of E. coli and S. coelicor A3(2) are very similar (approx. 1 kbp). In contrast, the S. coelicolor A3(2) genome approaches twice the size of the E. coli genome. The maximum specific growth rate observed for S. coelicor A3(2) approaches 0·35 h–1 (tD=1·98 h–1); see Shahab et al. (1996). The period, C, of DNA replication is thought to be directly proportional to the size, lg bp, of the genome [see equation (2)]. The duration of C is also dependent on the ambient temperature; for example, DNA synthesis at 30 °C is expected to proceed at approximately 64 % of the rate achieved at 37 °C (Ingraham & Marr, 1996). A minimum value for the rate of DNA replication was calculated for S. coelicolor growing at its maximum rate, as follows. It was supposed that tD*=C*+D*, that is B*->0 h; and, by analogy with the value for E. coli (see Table 3), D*->0·5 h, corresponding to a value for C* of 1·48 h. Hence, the value for {varepsilon}DNA*/3600 s h–1 was calculated to be 812 bp per replication fork s–1, which is 87 % of the maximum rate reported for E. coli growing at 37 °C (see Table 3). The DNA replication rate, adjusted for growth at 37 °C, is 1250 bp per replication fork s–1, which is 135 % of the maximum rate (933 bp per replication fork s–1) observed for E. coli. Thus, the high rate of DNA replication is a notable feature of S. coelicolor A3(2) growing at its maximum rate. The provisional values for C and D allow a provisional value of 1·56 genome equivalents per average cell to be calculated by means of equation (3).

The macromolecular compositions of S. coelicolor A3(2), grown at different specific growth rates are summarized in Table 4. The average numbers of genomes per cell (ng(av)) are provisional and more precise values may lead to an adjustment of dry cell mass (mdc(av)), protein content (mp(av)), etc., but would not affect ratios such as mRNA(av)/mdc(av). Equations (6a) and (6b) were used to analyse features of rRNA synthesis. In cells grown optimally (µmax=0·3 h–1) the rate of rRNA synthesis was found to be slow in comparison with the high rate of DNA synthesis (see Table 4). Briefly, the time taken to synthesize one copy of pre­rRNA was calculated to be 4·8 min, corresponding to {varepsilon}rRNA*/3600 s h–1=19 nucleotides s–1. In contrast with E. coli, the rate of initiation of pre­rRNA synthesis was found to be high (iRNAP/ilim->1), indicating that, when growing optimally, S. coelicolor A3(2) has little spare capacity for pre­rRNA synthesis.

The correlations between µ2 and mdc(av) (or v(av)), which are summarized by equations (19a) and (19b), are presented in Fig. 1.


It was estimated that =320 fg and mdc(av)*=430 fg (µ=0·3 h–1), corresponding to =1·0 fl and v(av)*=1·3 fl. Comparison of the slopes of equations (15a) and (19a) reveals that the growth­directed metabolic activity per fg dry cell mass of S. coelicolor A3(2) is less than one­sixth of the activity of E. coli B/r.

The relation between µ and the product of RNA concentration (mRNA(av)/mdv(av)) and the peptide chain elongation rate ({varepsilon}aa) is illustrated in Fig. 2 and is described by equation (10b), where {phi}=1x10–4, in accord with the theoretical values for E. coli. This result is to be expected if constants such as the sizes of the rRNA species, the fraction of ribosomes actively synthesizing protein, etc., are very similar for both species. The concentration of RNA was found to vary by a little more than twofold over a 12·5­fold increase in µ (see Fig. 2), according to equations (20a) and (20b):


Values of RNA concentrations ( ) for limiting cells were calculated (intercept/slope) to be 0·061 fg RNA per fg dry cell mass and 21 fg per fl cell volume, respectively, similar to values deduced for limiting cells of E. coli [see discussion following equations (17a) and (17b)]. Within experimental error (see Fig. 4), the protein synthesis rate was found to be proportional to the third power of the RNA concentration [see equations (21a) and (21b)]:


Properties of cells of S. coelicolor A3(2) and E. coli B/r synthesizing protein at similar rates (51 and 42 fg h–1, respectively) were compared. It was found that cells of S. coelicolor A3(2) (µ=0·30 h–1) had a higher number (22 000 ribosomes) and a higher concentration (18 200 ribosomes fl–1) of ribosomes than cells of E. coli (µ=0·42 h–1), which had 6800 ribosomes (11 000 ribosomes fl–1). Application of equation (9) reveals that the peptide chain elongation rate {varepsilon}aa/3600 s h–h is lower for S. coelicolor A3(2) than for E. coli B/r; namely 3·17 and 14·0 amino acid residues per ribosome s–1 respectively. These values of {varepsilon}aa are thought to reflect both the different temperatures of growth (30 and 37 °C, respectively) and the rate of interaction between the ternary complex formed between aminoacyl­tRNA, EF­Tu and GTP with the A­site of the ribosome. It is proposed that lower concentrations of aminoacyl­tRNA and of EF­Tu relative to ribosomes (Cox, 2003) would diminish both the rate of ternary complex formation and the rate of peptide bond formation ({varepsilon}aa). Similarly, the 13­fold differences in the numerical values of {theta} [see equation (12a)] found for E. coli B/r [equation (18a)] and S. coelicolor A3(2) [equation (21a)] are attributable to differences in both the temperature of growth and the number of copies per ribosome of aminoacyl­tRNA and EF­Tu. A rough estimate, which was made on the basis of the assumption that the kinetic constants for peptide bond formation are similar for both species, suggests that the numbers of copies of aminoacyl­tRNA per ribosome and EF­Tu per ribosome present in S. coelicolor A3(2) are one­third of the numbers of copies per ribosome present in E. coli.

The transcriptional and translational properties of S. coelicolor A3(2) growing (µ=0·30 h–1) near to its maximum rate (µ{approx}0·35 h–1) are summarized in Fig. 5(b) as a virtual schematic cell.

The analysis of the data for S. coelicolor A3(2) growing optimally leads to provisional values for {varepsilon}DNA, trrn, iRNAP, the number of copies of EF­Tu per ribosome and the number of aminoacyl­tRNA molecules per ribosome. Each of these parameters may be measured by experiment.

Genomic and macromolecular properties of M. bovis BCG
The genomes of M. tuberculosis and E. coli are very similar in their sizes, the numbers of ORFs per genome and the sizes of their ‘average’ proteins (see Table 1). The average macromolecular properties of cells grown at 37 °C (tD=24h, µ=0·029) are summarized in Table 5; data for M. bovis BCG are presented in Figs 3–5. The average cell volume was estimated as 0·96±0·06 fl (or µm3); see Table 5. The concentration of RNA (mRNA(av)/mdc(av)=0·042 fg RNA per fg dry cell mass) was found to be close to the minimum value inferred for viable cells (see Fig. 3).

Equations derived for E. coli were applied to M. bovis BCG. For example the value of the RNA concentration (mRNA(av)/mdc(av)) derived by means of equation (17a) for µ=0·029 h–1 was found to be 0·064 fg RNA per fg dry cell mass compared with the observed value (see Table 5) of 0·042 fg RNA per fg dry cell mass; substitution of {omega}p(av) for M. bovis BCG in equation (18a) yielded a value of mRNA(av)/mdc(av)=0·041 fg RNA per fg dry cell mass. In contrast, calculations based on equations (20a) and (21a) derived for S. coelicolor yielded estimates close to twice the observed value: namely 0·092 and 0·096 fg RNA per fg dry cell mass respectively. In this respect E. coli is a better model for M. bovis BCG than is S. coelicolor. Equation (10b), which was found to apply to both E. coli and S. coelicolor (see Fig. 2), was found to apply also to M. bovis BCG. Substitution of mRNA(av)/mdc(av)=0·042 fg RNA per fg dry cell mass and {varepsilon}aa=2 amino acids incorporated into protein per ribosome s–1 into the equation led to a value of µ=0·030 h–1 compared with the observed value of 0·029 h–1.

The activity of rrn operons is described by {varepsilon}rRNA and iRNAP. A value of {varepsilon}rRNA/3600 s h–1=12 residues s–1 was made on the basis of the assumption that {varepsilon}rRNA=6{varepsilon}aa residues h–1. There is no reported value for {varepsilon}rRNA of M. bovis BCG. However, a very close relative of M. tuberculosis was found (Harshey & Ramakrishnan, 1977) to synthesize pre­rRNA in 7·6 min, corresponding to {varepsilon}rRNA/3600 s h–1=12·2 nucleotides s–1 for pre­rRNA. An average value of iRNAP=12 RNAP complexes per rrn operon was calculated from {omega}R(av)/nrrn(av) for {varepsilon}rRNA/3600 s h–1=12·2 nucleotides s–1. The rate of initiation of transcription of rrn operons of M. bovis BCG is inferred to be close to 12 % of the theoretical maximum (iRNAP/ilim=0·12). The transcriptional and translational properties of M. bovis BCG are summarized in Fig. 5(c). The metabolic activity of the average cell is encapsulated by the parameters C*=10·3 h, trrn*·60 min h–1=7·6 min, tap*·60 min h–1=2·8 min, and {varepsilon}E/3600 s h–1=0·64x105 ATP equivalents s–1. The low rates of macromolecular synthesis require a low rate of nutrient uptake and generation of energy, which may be a valuable asset for growth within macrophages.

The observation that equations (17a) and (18a) obtained for E. coli appear valid for M. bovis BCG suggests that equation (7) relating µ to nrrn(av) may apply not only to E. coli but also to M. bovis BCG and other mycobacteria. Upper limits for the maximum specific growth rates of mycobacteria were equated with the values of µmax calculated for E. coli possessing either one or two rrn operons per genome, namely 0·69 h–1 and 1·00 h–1 respectively. Mycobacteria grow more slowly than these estimated values.

Mycobacterium marinum, which is a very close relative of M. tuberculosis as judged by the relatedness of their 16S rRNA gene sequences, fatty acid profile analysis and DNA–DNA hybridization (Tønjum et al., 1998), has a single rrn operon per genome (C. Helguera­Repetto, R. A. Cox & J. A. Gonzalez­y­Merchand, unpublished work). The maximum specific growth rate of M. marinum was found to be 0·173 h–1 (tD=4 h) at 30 °C (Clark & Shepard, 1963), which compares with 0·69 h–1 calculated for hypothetical cells of E. coli. Mycobacterium chelonae, which also possesses a single rrn operon per genome (Domenech et al., 1994), was found to grow at 30 °C with a µmax of 0·13 h–1 (tD=5·4 h) (M. C. Nuñez & M. J. Garcia, unpublished work). Mycobacterium smegmatis, which is known to have two rrn operons per genome (Bercovier et al., 1986), has a µmax of 0·28 h–1 (tD=2·5 h) at 37 °C (Gonzalez­y­Merchand et al., 1999), compared with tD=1·00 h calculated for hypothetical cells of E. coli.

The tubercle bacillus has a single rrn operon per genome, which was shown to be under growth rate control when transferred to M. smegmatis (Verma et al., 1999). When growing at its maximum rate (see above) M. bovis BCG was found to use about 12 % of the potential activity of its rrn operon. Thus, it is inferred that the slow growth of M. tuberculosis, M. bovis BCG and other members of the complex reflects cell metabolism and not the possession of a single rrn operon per genome. Members of the M. tuberculosis complex have genomes which have gained a large number (56 or more) of insertion sequences through horizontal transfer (Brosch et al., 2000, 2002), which may have attenuated their capacities for growth­directed cell metabolism, leading to slow growth.

Concluding remarks
The framework devised to explore the growth of bacteria was applied to three species with maximum growth rates ranging from µmax=0·029 h–1 to µmax=1·73 h–1 (see Fig. 5). The examination of the properties of E. coli indicates that the requirements for ultra­fast growth include: (i) highly efficient mechanisms for the uptake and metabolism of a nutrient source of carbon and energy; (ii) a capacity for the cell to become polyploid and thereby increase the average number of rrn operons per cell; (iii) a capacity for rapid rRNA synthesis (trrn*·60 min h–1=1·2 min); (iv) a capacity for rapid protein synthesis – for example, an average protein of 330 amino acid residues can be synthesized in 15 s; this rate is achieved through high concentrations of ribosomes, aminoacyl­tRNA and EF­Tu.

S. coelicolor A3(2) growing optimally was found to synthesize macromolecules at contrasting rates. Specifically, DNA synthesis was estimated to proceed at a rate similar to that observed for E. coli (µ=1·73 h–1). In contrast, the predicted rates for rRNA synthesis and for synthesis of an average protein are relatively slow, namely, trrn*·60 min h–1=4·8 min and tap*·60 mins h–1=1·74 min. However, each rrn operon appears to be fully engaged in pre­rRNA synthesis.

Equations derived to describe the growth of E. coli were found to describe the growth of M. bovis BCG, in contrast with equations for S. coelicolor A3(2). Using my equations relating the number of rrn operons per genome to maximum specific growth rate, it is expected that M. tuberculosis with one rrn operon should be capable of growing much faster than it actually does. Therefore, it is suggested that the high number of insertion sequences in this species attenuates growth to much lower values. M. tuberculosis has the ability to persist in the form of a long­term asymptomatic infection (for review see Stewart et al., 2003) and knowledge of the properties of minimal cells of this pathogen is relevant to this phenomenon.

Finally, the framework described in the Theoretical Analyses section is based on the availability of a minimum set of data which includes: (i) the size of genome; (ii) the number of rrn operons per genome; (iii) the DNA : RNA : protein ratios obtained at several growth rates and, preferably, a knowledge of the number of copies of aminoacyl­tRNA and EF­Tu per ribosome. These data provide a range of insights into cell physiology.


   APPENDIX
TOP
ABSTRACT
INTRODUCTION
THEORETICAL ANALYSES
RESULTS AND DISCUSSION
APPENDIX
REFERENCES
 
Correlation between specific growth rate and the average number of rrn operons per cell
When a bacterial culture is growing exponentially both protein content and RNA content increase exponentially. The specific protein synthesis rate, {omega}p(av), of an average cell is defined by equation (1A) (all symbols are defined in Table 2):

Rearrangement leads to equation (2A) [equation (10a) of the Theoretical Analyses section]:

However, nR(av)/naa(av) is proportional to mRNA(av)/mdc(av), leading to equation (3A) [equation (10b) of the Theoretical Analyses section], where {phi}={beta}R·{beta}rRNA·(Mr[aa(av)]/Mr[nuc(av)])/({beta}m·lrRNA); mp(av)={beta}m·mdc(av); naa(av)=mp(av)/(maa(av)/NA) and nR(av)={beta}rRNA·[mRNA(av)/lrRNA·Mr[nuc(av)]/NA)], where Mr[aa(av)] is the average mass of an amino acid residue, NA is Avogadro's constant, {beta}rRNA is the fraction of RNA that is rRNA, lrRNA is the number of nucleotides per ribosome, and Mr[nuc(av)] is the average mass of a nucleotide.

For E. coli the constant {phi} has a numerical value of 1x10–4.

The specific ribosome synthesis rate is defined by equation (4A):

Implicit in this equation is the assumption that all rrn operons are equally active in rRNA synthesis. Rearrangement of equation (4A) leads to equation (5A):

Multiplying equation (2A) by equation (5A) leads to equation (6A):

When particular values, , and , are assigned to iRNAP, {varepsilon}rRNA and naa(av), respectively, equation (6A) reduces to equation (7A), where . µx and xrrn(av) are unknown values of, respectively, µ and nrrn(av):

Equation (7A) was used in the following way. First, appropriate values were assigned to the product iRNAP·{varepsilon}rRNA and to naa(av). A value, µx, was then assigned to µ allowing the specific protein synthesis rate ({omega}=mp(av)·µx) to be calculated. The ratio mRNA(av)/mdc(av) was then obtained by means of equation (12) of the Theoretical Analyses section, allowing {varepsilon}aa to be evaluated by means of equation (3A). Then, the only unknown quantity is xrrn(av). Provisionally, a value of 6{varepsilon}aa residues h–1 was assigned to {varepsilon}rRNA, allowing iRNAP to be calculated.

The validity of equations (6A) and (7A) was established in the following way. The effects of reducing the number of rrn perons per genome on µ were calculated for E. coli B/r grown in medium E. Wild­type E. coli was found to have he following properties (see Table 3): µ=1·73 h–1, nrrn(av)*=35·9 operons, iRNAP*=68 copies/operon, {varepsilon}rRNA/3600 s h–1=85 residues s–1. Please note that the reported value of {varepsilon}rRNA is less than 6{varepsilon}aa, possibly owing to feedback control of rRNA synthesis. An rrn operon functioning at its maximum rate is characterized by iRNAP=75 copies per rrn operon and {varepsilon}rRNA/3600 s h–1=126 residues s–1 [see equation (6c) of the Theoretical Analyses section]. The minimum value of nrrn(av) required for µ=1·73 h–1 was calculated on the basis of the assumption that reducing nrrn(av) would lead to the remaining operons functioning at full capacity; that is iRNAP=75 copies per rrn operon, {varepsilon}rRNA/3600 s h–1=6{varepsilon}aa/3600 s h–1=126 residues s–1, naa(av)=8·5x109 amino acids, {gamma}F=7·25x10–3, µ=1·73. The value of nrrn(av)=19·7 operons, which was obtained by means of equation (7A), is equivalent to four rrn operons per genome [nrrn(av)=(4/7)·35·9=20·5 operons]. A similar calculation revealed that a value of µ=1·5 h–1 (tD=0·47 h) required a minimum value of nrrn(av)=15·5 operons, corresponding to three rrn operons per genome. These estimates are in accord with the results obtained by progressively reducing the number of functional rrn operons present in the E. coli genome (Condon et al., 1993, 1995).


   ACKNOWLEDGEMENTS
 
I thank Simon A. Cox for his help in the preparation of this manuscript, Dr J. A. Gonzalez-y-Merchand, Escuela Nacional de Ciencias Biológicas, IPN, Mexico, for stimulating my interest in Streptomyces, Dra M. C. Nuñez and Dra M. J. Garcia of the Universidad Autonoma de Madrid, and Dr J. A. Gonzalez-y-Merchand and Dra C. Helguera-Repetto of IPN, Mexico, for permission to cite their unpublished work, and my colleague I. D. J. Burdett for his continued interest.


   REFERENCES
TOP
ABSTRACT
INTRODUCTION
THEORETICAL ANALYSES
RESULTS AND DISCUSSION
APPENDIX
REFERENCES
 
Armstrong, J. A. & D'Arcy-Hart, P. (1971). Response of cultured macrophages to Mycobacterium tuberculosis, with observations on fusion of lysosomes with phagosomes. J Exp Med 134, 713–740.

Bentley, S. D., Chater, K. F., Cerdeño-Tárraga, A.-M. & 40 other authors (2002). Complete genome sequence of the model actinomycete Streptomyces coelicolor A3(2). Nature 417, 141–147.[CrossRef][Medline]

Bercovier, H., Kafri, O. & Sela, S. (1986). Mycobacteria possess a surprisingly small number of ribosomal RNA genes in relation to the size of their genome. Biochem Biophys Res Commun 136, 1136–1141.[Medline]

Bernstein, J. A., Khodursky, A. B., Lin, P.-S., Lin-Chao, S. & Cohen, S. N. (2002). Global analysis of mRNA decay and abundance in Escherichia coli at single-gene resolution using two-color fluorescent DNA microarrays. Proc Natl Acad Sci U S A 99, 9697–9702.[Abstract/Free Full Text]

Blattner, F. R., Plunket, G., III, Bloch, C. A. & 14 other authors (1997). The complete genome sequence of Escherichia coli K-12. Science 277, 1453–1462.[Abstract/Free Full Text]

Bremer, H. & Dennis, P. P. (1996). Modulation of chemical composition and other parameters of the cell growth rate. In Escherichia coli and Salmonella: Cellular and Molecular Biology, 2nd edn, pp. 1553–1568. Edited by F. C. Neidhardt and others. Washington DC: American Society for Microbiology.

Brosch, R., Gordon, S. V., Eiglmeier, K., Garnier, T., Tekaia, F., Yeramian, E. & Cole, S. T. (2000). Genomics, biology, and evolution of the Mycobacterium tuberculosis complex. In Molecular Genetics of Mycobacteria, pp. 19–36. Edited by G. F. Hatfull & W. R. Jacobs, Jr. Washington, DC: American Society for Microbiology.

Brosch, R., Gordon, S. V., Marmiesse, M. & 12 other authors (2002). A new evolutionary scenario for the Mycobacterium tuberculosis complex. Proc Natl Acad Sci U S A 99, 3684–3689.[Abstract/Free Full Text]

Byrne, R., Levin, J. G., Bladen, H. A. & Nirenberg, M. W. (1964). The in vitro formation of a DNA–ribosome complex. Proc Natl Acad Sci U S A 52, 140–148.[Medline]

Clark, H. F. & Shepard, C. C. (1963). Effect of environmental temperatures on infection with Mycobacterium marinum (balnei) of mice and a number of poikilothermic species. J Bacteriol 86, 1057–1069.

Cole, S. T., Brosch, R., Parkhill, J. & 39 other authors (1998). Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature 393, 537–544.[CrossRef][Medline]

Condon, C., French, S., Squires, C. & Squires, C. L. (1993). Depletion of functional ribosomal RNA operons in Escherichia coli causes increased expression of the remaining intact copies. EMBO J 12, 4305–4315.[Abstract]

Condon, C., Liveris, D., Squires, C., Schwartz, I. & Squires, C. L. (1995). rRNA operon multiplicity in Escherichia coli and the physiological implications of rrn inactivation. J Bacteriol 177, 4152–4156.[Abstract]

Cox, R. A. (2003). Correlation of the rate of protein synthesis and the third power of the RNA : protein ratio in Escherichia coli and Mycobacterium tuberculosis. Microbiology 149, 729–737.[Abstract/Free Full Text]

Domenech, P., Menendez, M. C. & Garcia, M. J. (1994). Restriction fragment length polymorphism of 16S rRNA genes in the differentiation of fast-growing mycobacterial species. FEMS Microbiol Lett 116, 19–21.[CrossRef][Medline]

Ferrari, G., Langen, H., Naito, M. & Pieters, J. (1999). A coat protein on phagosomes involved in the intracellular survival of mycobacteria. Cell 97, 435–447.[Medline]

Gonzalez-y-Merchand, J. A., Colston, M. J. & Cox, R. A. (1999). Effects of growth conditions on the expression of mycobacterial murA and tyrS genes and contributions of their transcripts to precursor rRNA synthesis. J Bacteriol 181, 4617–4627.[Abstract/Free Full Text]

Grunberg-Manago, M. (1999). Messenger RNA stability and its role in the control of gene expression in bacteria and phages. Annu Rev Genet 33, 193–227.[CrossRef][Medline]

Harshey, R. M. & Ramakrishnan, T. (1977). Rate of ribonucleic acid chain growth in Mycobacterium tuberculosis H37Rv. J Bacteriol 129, 616–622.[Medline]

Helmstetter, C. E. (1996). Timing of synthetic activities in the cell cycle. In Escherichia coli and Salmonella: Cellular and Molecular Biology, 2nd edn, pp. 1627–1639. Edited by F. C. Neidhardt and others. Washington, DC: American Society for Microbiology.

Helmstetter, C. E. & Cooper, S. (1968). DNA synthesis during the division cycle of rapidly growing E. coli B/r. J Mol Biol 31, 507–518.[Medline]

Hiriyanna, K. T. & Ramakrishnan, T. (1986). Deoxyribonucleic acid replication time in Mycobacterium tuberculosis H37Rv. Arch Microbiol 144, 105–109.[Medline]

Ingraham, J. L. & Marr, J. L. (1996). Effect of temperature, pressure, pH, and osmotic stress on growth. In Escherichia coli and Salmonella: Cellular and Molecular Biology, 2nd edn, pp. 1570–1578. Edited by F. C. Neidhardt and others. Washington, DC: American Society for Microbiology.

Ingraham, J. L., Maaløe, O. & Neidhardt, F. C. (1983). Growth of the Bacterial Cell. Sunderland, MA: Sinauer Associates.

Jacobsen, H. (1974). PhD thesis, University of Copenhagen.

Khodursky, A. R. & Bernstein, J. A. (2003). Life after transcription – revisiting the fate of messenger RNA. Trends Genet 19, 113–115.[CrossRef][Medline]

Krummel, B. & Chamberlin, M. J. (1989). RNA chain initiation by Escherichia coli RNA polymerase. Structural transitions of the enzyme in early ternary complexes. Biochemistry 28, 7829–7842.[Medline]

Leroy, A., Vanzo, N. F., Sousa, S., Dreyfus, M. & Carpousis, A. J. (2002). Function in Escherichia coli of the non-catalytic part of RNase E : role in the degradation of ribosome-free mRNA. Mol Microbiol 45, 1231–1243.[CrossRef][Medline]

Miller, O. L., Jr, Hamkalo, B. A. & Thomas, C. A., Jr (1970). Visualization of bacterial genes in action. Science 169, 392–395.[Medline]

Neidhardt, F. C. (1996). The enteric bacterial cell and the age of bacteria. In Escherichia coli and Salmonella: Cellular and Molecular Biology. 2nd edn, pp. 1–3. Edited by F. C. Neidhart and others. Washington, DC: American Society for Microbiology.

Noller, H. E. & Nomura, M. (1996). Ribosomes. In Escherichia coli and Salmonella: Cellular and Molecular Biology, 2nd edn, pp. 167–186. Edited by F. C. Neidhardt and others. Washington, DC: American Society for Microbiology.

Régnier, P. & Arraiano, C. M. (2000). Degradation of mRNA in bacteria: emergence of ubiquitous features. BioEssays 22, 235–244.[CrossRef][Medline]

Schaechter, E., Maaløe, O. & Kjeldgaard, N. O. (1958). Dependence on medium and temperature of cell size and chemical composition during balanced growth of Salmonella typhimurium. J Gen Microbiol 19, 592–606.

Shahab, N., Flett, F., Oliver, S. G. & Butler, P. R. (1996). Growth rate control of protein and nucleic acids in Streptomyces coelicolor A3(2) and Escherichia coli B/r. Microbiology 142, 1927–1935.[Abstract]

Shepard, C. C. (1960). The experimental disease that follows the injection of human leprosy bacilli into the footpads of mice. J Exp Med 112, 445–454.

Stent, G. S. (1964). The operon: on its third anniversary. Science 144, 816–820.[Medline]

Stewart, G. R., Robertson, B. D. & Young, D. B. (2003). Tuberculosis: a problem with persistence. Nat Rev 1, 97–105.[CrossRef]

Tonella, L., Walsh, B. J., Sanchez, J. C. & 14 other authors (1998). '98 Escherichia coli SWISS-2DPAGE database update. Electrophoresis 19, 1960–1971.[Medline]

Tønjum, T., Welty, D. B., Jantzen, E. & Small, P. L. (1998). Differentiation of Mycobacterium ulcerans, M. marinum, and M. haemophilum: mapping of their relationships to M. tuberculosis by fatty acid profile analysis, DNA–DNA hybridization and 16S rRNA gene sequence analysis. J Clin Microbol 35, 918–925.

Verma, A., Sampla, A. K. & Tyagi, J. S. (1999). Mycobacterium tuberculosis rrn promoters: differential usage and growth rate-dependent control. J Bacteriol 181, 4326–4333.[Abstract/Free Full Text]

Wayne, L. G. (1994). Cultivation of Mycobacterium tuberculosis for research purposes. In Tuberculosis: Pathogenesis, Protection and Control, pp. 73–83. Edited by B. Bloom. Washington, DC: American Society for Microbiology.

Wayne, L. G. & Hayes, L. G. (1996). An in vitro model for sequential study of shift down of Mycobacterium tuberculosis through two stages of replicating persistence. Infect Immun 64, 2062–2069.[Abstract]

Wayne, L. G. & Kubica, G. P. (1986). The mycobacteria. In Bergey's Manual of Systematic Bacteriology, vol. 2, pp. 1435–1457. Edited by P. H. A. Sneath, N. S. Mair, M. E. Sharpe & J. G. Holt. London: Williams & Wilkins.

Winder, F. G. & Rooney, S. A. (1970). Effects of nitrogenous components of the medium on the carbohydrate and nucleic acid content of Mycobacterium tuberculosis BCG. J Gen Microbiol 63, 29–39.[Medline]

Woods, G. L. & Washington, J. A. (1987). Mycobacteria other than Mycobacterium tuberculosis: review of microbiologic and clinical aspects. Rev Infect Dis 9, 275–294.[Medline]

Received 12 June 2003; revised 2 December 2003; accepted 20 January 2004.



This Article
Abstract
Full Text (PDF)
Alert me when this article is cited
Alert me if a correction is posted
Citation Map
Services
Email this article to a friend
Similar articles in this journal
Similar articles in PubMed
Alert me to new issues of the journal
Download to citation manager
Google Scholar
Articles by Cox, R. A.
Articles citing this Article
PubMed
PubMed Citation
Articles by Cox, R. A.
Agricola
Articles by Cox, R. A.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS
Copyright © 2004 Society for General Microbiology.